An Update to The Demagnetizing Factor Dataset Calculated for The General Ellipsoid by Osborn

: The exact formulae for calculating the demagnetizing factors of a general ellipsoid along the three main axes a ≥ b ≥ c have been long known. According to these formulae, the demagnetizing factors depend only on the axial ratios b / a and c / a . Although the calculation of the demagnetizing factors is a straightforward task, the calculation itself is not a simple one. Therefore, tabular and graphical representations of these demagnetizing factor data have also been presented which can then be used for approximating the demagnetizing factors of a rectangular ferromagnetic slab with the same axial ratios. It turned out in our recent study, however, that, in some ranges of axial ratios (e.g., for very small c / a values), the available tables and graphs do not provide sufficient resolution for obtaining the demagnetizing factors with reasonable accuracy. It was decided to calculate these missing values, and they are presented here in both tabular and graphical form by giving instructions for how to obtain conveniently further interpolated data. In addition, the previous and current demagnetizing factor data have been replotted and fitted to a polynomial function with high accuracy. The functional form of these fitting polynomials is presented in a table for the whole range of the axial ratios b / a and c / a . By graphically displaying these functions, one can obtain, in a relatively simple manner, the demagnetizing factors of a general ellipsoid with known axial ratios without the need to directly calculate through the exact formulae. This may be helpful in obtaining a quick estimate for the demagnetizing factors of any rectangular ferromagnetic slab of interest.


Introduction
The determination of the anisotropic magnetoresistance (AMR) [1,2] requires the measurement of the field dependence of the resistivity in the configurations when the magnetic field H is parallel to the measuring current (this is called longitudinal magnetoresistance, LMR) and when H is perpendicular to the current (this is the transverse magnetoresistance, TMR).The measurement of the MR(H) curves can be conveniently carried out on thin stripshaped samples with a current flowing along the long axis of the strip, and the magnetic field H is oriented in the strip plane either parallel or perpendicular to the current flow direction [3].
It turned out, in some recent MR studies of various ferromagnetic materials [4][5][6][7][8], that the shape and width of the low-field LMR(H) and TMR(H) curves were similar for some of the samples whereas they were distinctly different for other samples.The strip-shaped foil samples for the MR(H) measurements were typically 5 to 10 mm long, 1 to 2 mm wide, and 10 to 50 µm thick.In order to reveal the origin of the observed differences in the LMR(H) and TMR(H) curves, magnetization (M) measurements have been subsequently performed [9] on similar strip-shaped samples.These M(H) studies [9] revealed that the differences between the longitudinal and transverse configurations can be well explained by demagnetizing field effects.Two examples are shown in Figure 1 where the magnetization is displayed as a function of the external magnetic field for a cold-rolled Ni 50 Co 50 strip and an electrodeposited nanocrystalline (nc) Ni strip, with the magnetic field oriented in the strip plane parallel (M L ) and transverse (M T ) to the longest edge.These samples have the largest and smallest thickness/length ratio, respectively, showing the huge effect of sample geometry on the transverse magnetization curves.The detailed evaluation of the M(H) study to be published elsewhere [9] necessitated a thorough consideration of the demagnetizing factors for such strip-shaped samples, and this was a major motivation for the present work.where the magnetization is displayed as a function of the external magnetic field for a cold-rolled Ni 50 Co 50 strip and an electrodeposited nanocrystalline (nc) Ni strip, with the magnetic field oriented in the strip plane parallel (M L ) and transverse (M T ) to the longest edge.These samples have the largest and smallest thickness/length ratio, respectively, showing the huge effect of sample geometry on the transverse magnetization curves.The detailed evaluation of the M(H) study to be published elsewhere [9] necessitated a thorough consideration of the demagnetizing factors for such strip-shaped samples, and this was a major motivation for the present work.From the viewpoint of the demagnetizing effects, the strip-shaped thin foil samples can be considered as a rectangular ferromagnetic slab which, on the other hand, can be approximated by a general ellipsoid.The exact formulae for calculating the demagnetizing factors of a general ellipsoid along the three main axes a ≥ b ≥ c have been long known [10,11].According to these formulae, the demagnetizing factors depend only on the axial ratios b/a and c/a.Osborn [10] presented a tabular and graphical representation of the demagnetizing factor data for the general ellipsoid which was then complemented with some further data by Cronemeyer [12].
It should be noted that, following the work of Osborn [10], Stoner [11], and Cronemeyer [12], numerous further studies (see, e.g., Refs.[13][14][15][16][17][18][19][20][21]) have been carried out for calculating the demagnetizing factors of both the general ellipsoid and its specific forms and also various rectangular or circular ferromagnetic slabs, such as rods and disks, by accounting in numerous cases also for the fact that, in not perfectly ellipsoidal objects, the magnetization distribution within the specimen is not homogeneous, in contrast to the general ellipsoid.However, by looking at these reports [13][14][15][16][17][18][19][20][21], it turned out that, whereas they are really useful for the specific cases considered, e.g., for infinite cylinders or rods, the results cannot be straightforwardly applied to the strip-shaped samples with the axial ratios of our interest.
Therefore, we have found it more useful to follow the scheme to rely on the formulae of Osborn [10] and Stoner [11] to directly calculate the required demagnetizing factors of interest for our study.Namely, it turned out that the approximation of a rectangular From the viewpoint of the demagnetizing effects, the strip-shaped thin foil samples can be considered as a rectangular ferromagnetic slab which, on the other hand, can be approximated by a general ellipsoid.The exact formulae for calculating the demagnetizing factors of a general ellipsoid along the three main axes a ≥ b ≥ c have been long known [10,11].According to these formulae, the demagnetizing factors depend only on the axial ratios b/a and c/a.Osborn [10] presented a tabular and graphical representation of the demagnetizing factor data for the general ellipsoid which was then complemented with some further data by Cronemeyer [12].
It should be noted that, following the work of Osborn [10], Stoner [11], and Cronemeyer [12], numerous further studies (see, e.g., Refs.[13][14][15][16][17][18][19][20][21]) have been carried out for calculating the demagnetizing factors of both the general ellipsoid and its specific forms and also various rectangular or circular ferromagnetic slabs, such as rods and disks, by accounting in numerous cases also for the fact that, in not perfectly ellipsoidal objects, the magnetization distribution within the specimen is not homogeneous, in contrast to the general ellipsoid.However, by looking at these reports [13][14][15][16][17][18][19][20][21], it turned out that, whereas they are really useful for the specific cases considered, e.g., for infinite cylinders or rods, the results cannot be straightforwardly applied to the strip-shaped samples with the axial ratios of our interest.
Therefore, we have found it more useful to follow the scheme to rely on the formulae of Osborn [10] and Stoner [11] to directly calculate the required demagnetizing factors of interest for our study.Namely, it turned out that the approximation of a rectangular ferromagnetic slab with a general ellipsoid having the same axial ratios is a fairly good approach since direct experimental determinations of the demagnetizing factors of rectan-gular slab samples [9,14] have demonstrated a not-too-large deviation.As noted above, such a deviation arises due to the fact that, even in a homogeneous external magnetic field, the magnetization orientation is homogeneous only in a body having the shape of a general ellipsoid.Therefore, in a rectangular slab, the magnetization orientation is inhomogeneous to some extent around the edges [14] as a result of which the effective demagnetizing factor will be somewhat different [9,14] than that calculated for a general ellipsoid having the same axial ratios.The recent magnetic measurements [9] revealed that the experimental demagnetizing factor N b when the magnetic field is oriented along the slab edge b is about 8% larger only than the N b value calculated for the corresponding general ellipsoid.
It is, therefore, useful to have available detailed data for the demagnetizing factors of the general ellipsoid for the whole range of axial ratios.It turned out, however, in our recent study [9], that, in some ranges of axial ratios (e.g., for very small c/a values), the available tables and graphs [10,12] do not provide sufficient resolution for obtaining the demagnetizing factors with reasonable accuracy.The two samples shown in Figure 1 correspond to the smallest and largest c/a values (c/a = 0.00074 for nc-Ni and c/a = 0.01533 for Ni 50 Co 50 ) in our recent work [9].Thus, the c/a values for our samples fall just in the range where the demagnetizing factors cannot be obtained with reasonable accuracy from previously existing sources.
It was decided, therefore, to calculate these missing values, and they are presented here in both tabular and graphical form by giving instructions for how to conveniently obtain further interpolated data.In addition, the previous and current demagnetizing factor data [10,12] have been replotted and fitted to a polynomial function with high accuracy.The functional form of these fitting polynomials is presented in a table for the whole range of the axial ratios b/a and c/a.By graphically displaying these functions, one can obtain, in a relatively simple manner, the demagnetizing factors of a general ellipsoid with known axial ratios without the need to directly calculate through the exact formulae.This may be helpful in obtaining a quick estimate for the demagnetizing factors of any rectangular ferromagnetic slab of interest.

Demagnetizing Factors of the General Ellipsoid
The magnetic induction inside a ferromagnetic specimen is given by the general expression [22] B = H − NM + 4πM (1) where H is the external magnetic field and N is the demagnetizing factor (we use here the CGS system).Strictly speaking, this expression is valid only for a homogeneously magnetized specimen in the form of a general ellipsoid.In this case, the demagnetizing field H d is also uniform within the specimen and is proportional to the magnetization: H d = −NM; and this appears in Expression (1).The demagnetizing factor is a scalar quantity along the three main axes a ≥ b ≥ c of the ellipsoid, and for the three demagnetizing factors, the relation holds.Osborn [10] and Stoner [11] have provided exact formulae for the calculation of the demagnetizing factors of the general ellipsoid.
As outlined in the Introduction, a general ellipsoid will be used for approximating the demagnetizing factors of a rectangular slab shown in Figure 2. In this approach, we will approximate the slab with a corresponding inscribed general ellipsoid having the same full axis lengths as the edges of the rectangular slab.(Since the demagnetizing factors depend only on the axial ratios b/a and c/a, any general ellipsoid with the same axial ratios is an equally appropriate approximation, the term "inscribed" does not represent any restriction, and the same is valid whether the parameters a, b and c are defined as full axes or semiaxes.)ratios is an equally appropriate approximation, the term "inscribed" does not represent any restriction, and the same is valid whether the parameters a, b and c are defined as full axes or semiaxes.)

New Calculated Data for the Demagnetizing Factors of the General Ellipsoid
As noted in the Introduction, Osborn [10] and Cronemeyer [12] presented tabular values and graphical representations of the dependence of the demagnetizing factors of the general ellipsoid on the axial ratios c/a and b/a for a wide range of the axial ratios.Unfortunately, no useful data were included in these works [10,12] for axial ratios c/a < 0.1.For completing these missing data, by using formulae (2.1)-(2-3) from Osborn [10], we have calculated a set of N/4π values in this c/a range for b/a values of integer multiples of 0.1, and these values are collected in Table 1.
Table 1.Demagnetizing factors of the general ellipsoid as calculated in the present work from formulae (2.1)-(2.3)from Osborn [10] for selected b/a and c/a values.The N/4π values for c/a = 0.1 were taken from Refs.[10,12] in the case of b/a = 0.1 to 0.9.For b/a = 1.0, we obtained the N/4π values from the formulae for oblate spheroids [10,11,22].

New Calculated Data for the Demagnetizing Factors of the General Ellipsoid
As noted in the Introduction, Osborn [10] and Cronemeyer [12] presented tabular values and graphical representations of the dependence of the demagnetizing factors of the general ellipsoid on the axial ratios c/a and b/a for a wide range of the axial ratios.Unfortunately, no useful data were included in these works [10,12] for axial ratios c/a < 0.1.For completing these missing data, by using formulae (2.1)-(2-3) from Osborn [10], we have calculated a set of N/4π values in this c/a range for b/a values of integer multiples of 0.1, and these values are collected in Table 1.
Table 1.Demagnetizing factors of the general ellipsoid as calculated in the present work from formulae (2.1)-(2.3)from Osborn [10] for selected b/a and c/a values.The N/4π values for c/a = 0.1 were taken from Refs.[10,12] in the case of b/a = 0.1 to 0.9.For b/a = 1.0, we obtained the N/4π values from the formulae for oblate spheroids [10,11,22]

A Useful Graphical Representation of the Demagnetizing Factor Data for the General Ellipsoid
In order to provide help for obtaining N values for more c/a and b/a values than available in the published tables and graphs [10,12] and in the present Table 1, it was found useful to replot the available data.Therefore, we have prepared graphs of N a /4π and N b /4π as a function of c/a for all values of b/a for which data are available.The datasets N/4π vs. c/a obtained in this manner were fitted for each available value of b/a to a polynomial under the constraint that the polynomial should be zero (N = 0) at c/a = 0.
Examples of these graphs are shown in Figure 3 for the demagnetizing factor N a /4π and in Figure 4 for N b /4π.These graphs demonstrate that the available N a and N b data for b/a = 0.1, 0.2, and 0.3 can be fitted to a fourth-order polynomial with a very good fit quality.The same order of polynomial was sufficient to obtain similarly good fits for N a and N b with b/a = 0.4 to 1.0.The fit quality R2 was at least 0.9999 or even higher for all the fits.The parameters of the polynomial fitting functions are collected in Table 2.
By using these fitting functions, one can now display the demagnetizing factors as a function of c/a for several fixed values of b/a which fall close to the b/a value of the rectangular slab sample of interest.Such plots can now be easily created even for the range of extremely small c/a values which cannot be resolved properly in the graphs presented by Osborn [10].The good quality of fits ensures that we can reliably estimate the N data in the whole range of c/a values.
Since the b/a values of actual samples are always somewhat different from the fixed b/a values selected as integral multiples of 0.1 only, one should make a non-linear interpolation for actual b/a values in between the displayed polynomial curves in order to obtain a good N value for a given c/a value.We can see, for example, from Table 1, that for c/a = 0.   [10] and Cronemeyer [12] or given in Table 1.The solid lines are fourth-order polynomial fits to the displayed N a values under the constraint that N a should be zero at c/a = 0.
Since the b/a values of actual samples are always somewhat different from the fixed b/a values selected as integral multiples of 0.1 only, one should make a non-linear interpolation for actual b/a values in between the displayed polynomial curves in order to obtain a good N value for a given c/a value.We can see, for example, from Table 1, that for c/a = 0.   [10] and Cronemeyer [12] or given in Table 1.The solid lines are fourth-order polynomial fits to the displayed N a values under the constraint that N a should be zero at c/a = 0.

Summary
In the present paper, demagnetizing factor data are reported for the general ellipsoid, which were missing in previous works [10,12].In addition, all available demagnetizing factor data of the general ellipsoid were displayed as a function of the axial ratio c/a and fitted by a polynomial for each b/a value.The parameters of the fitting polynomials are provided here from which a graphical representation of the demagnetizing factors can be easily created.Some hints were also provided for obtaining interpolated values

Summary
In the present paper, demagnetizing factor data are reported for the general ellipsoid, which were missing in previous [10,12].In addition, all available demagnetizing factor data of the general ellipsoid were displayed as a function of the axial ratio c/a and fitted by a polynomial for each b/a value.The parameters of the fitting polynomials are provided here from which a graphical representation of the demagnetizing factors can be easily created.Some hints were also provided for obtaining interpolated values not displayed.An advantage of the suggested procedure is that one can obtain a good value for the demagnetizing factor values also in axial ratio ranges not available in previous reports [10,12].It is hoped that the present paper will be useful for quickly obtaining a reliable estimate of the demagnetizing factors for a rectangular slab of interest.

Figure 1 .
Figure 1.Magnetization M as a function of external magnetic field H for a cold-rolled Ni 50 Co 50 strip and an electrodeposited nanocrystalline (nc) Ni strip, with in-plane magnetic field oriented parallel (M L , black and red curves for Ni 50 Co 50 and nc-Ni, respectively) and transverse (M T , green and blue curves for Ni 50 Co 50 and nc-Ni, respectively) to the longest edge.

Figure 1 .
Figure 1.Magnetization M as a function of external magnetic field H for a cold-rolled Ni 50 Co 50 strip and an electrodeposited nanocrystalline (nc) Ni strip, with in-plane magnetic field oriented parallel (M L , black and red curves for Ni 50 Co 50 and nc-Ni, respectively) and transverse (M T , green and blue curves for Ni 50 Co 50 and nc-Ni, respectively) to the longest edge.

Figure 2 .
Figure 2. The dimensional parameters of a rectangular slab.The demagnetizing factors N a , N b, and N c correspond to the different orientations of the magnetic field (actually, the magnetization) as indicated.

Figure 2 .
Figure 2. The dimensional parameters of a rectangular slab.The demagnetizing factors N a , N b , and N c correspond to the different orientations of the magnetic field (actually, the magnetization) as indicated.
1, we find that N b /4π = 0.4899 for b/a = 0.1, N b /4π = 0.3183 for b/a = 0.2, and N b /4π = 0.2322 for b/a = 0.3.One easily finds that the difference in the N b /4π values between b/a = 0.1 and b/a = 0.2 is almost precisely a factor of two larger than that between b/a = 0.2 and b/a = 0.3 (0.1716 and 0.0861, respectively).This implies that the N b /4π value for b/a = 0.15 can be safely obtained by taking a subdivision at a ratio of 2:1 between the N values for b/a = 0.1 and b/a = 0.2, and then we end up with N b /4π = 0.3755 for b/a = 0.15 at c/a = 0.1.With a similar non-linear interpolation, one can straightforwardly obtain N b data for any further intermediate b/a values by using Figure 4.The same procedure can also be used for interpolating N a data with the help of graphs like Figure 3.

Figure 3 .
Figure 3. Plots of the demagnetizing factor N a /4π of the general ellipsoid as a function of c/a for three selected values of b/a.The symbols are the values tabulated by Osborn[10] and Cronemeyer[12] or given in Table1.The solid lines are fourth-order polynomial fits to the displayed N a values 1, we find that N b /4π = 0.4899 for b/a = 0.1, N b /4π = 0.3183 for b/a = 0.2, and N b /4π = 0.2322 for b/a = 0.3.One easily finds that the difference in the N b /4π values between b/a = 0.1 and b/a = 0.2 is almost precisely a factor of two larger than that between b/a = 0.2 and b/a = 0.3 (0.1716 and 0.0861, respectively).This implies that the N b /4π value for b/a = 0.15 can be safely obtained by taking a subdivision at a ratio of 2:1 between the N values for b/a = 0.1 and b/a = 0.2, and then we end up with N b /4π = 0.3755 for b/a = 0.15 at c/a = 0.1.With a similar non-linear interpolation, one can straightforwardly obtain N b data for any further intermediate b/a values by using Figure 4.The same procedure can also be used for interpolating N a data with the help of graphs like Figure 3.

Figure 3 .
Figure 3. Plots of the demagnetizing factor N a /4π of the general ellipsoid as a function of c/a for three selected values of b/a.The symbols are the values tabulated by Osborn[10] and Cronemeyer[12] or given in Table1.The solid lines are fourth-order polynomial fits to the displayed N a values under the constraint that N a should be zero at c/a = 0.

Magnetism 2024, 4 , 7 Figure 4 .
Figure 4. Plots of the demagnetizing factor N b /4π of the general ellipsoid as a function of c/a for three values of b/a.For an explanation of the symbols and lines, see Figure 3.Table 2. Fit functions y = f(x) = k 4 x 4 + k 3 x 3 + k 2 x 2 + k 1 xobtained by fitting a fourth-order polynomial to the N/4π vs. c/a data from Refs.[10,12] as well as from Table1and graphs shown in Figures3 and 4for various fixed values of b/a.The N c /4π values can then be obtained from Equation (2).DemagnetizingFitted Polynomials for N /4π and N /4π vs. c/a Data

Figure 4 .
Figure 4. Plots of the demagnetizing factor N b /4π of the general ellipsoid as a function of c/a for three values of b/a.For an explanation of the symbols and lines, see Figure 3.

N
a and in Figure 6 for N b .With the help of the gridlines, from these plots, one can make a fairly good reading of all four coefficients for the actual b/a value of the sample of interest.Then, by inserting these values into the fitting equation y = f(x) = k 4 x 4 + k 3 x 3 + k 2 x 2 + k 1 x together with the actual x = c/a value of the sample, one can directly obtain the N a and N b values for any b/a and c/a values with sufficient accuracy.

Figure 5 .
Figure 5. Plots of the four coefficients of the fitting equation y = f(x) = k 4 x 4 + k 3 x 3 + k 2 x 2 + k 1 x for the demagnetizing factor N a /4π of the general ellipsoid as a function of b/a.The symbols are the fitted coefficient values as given in Table2, and the lines are approximate trend lines through the data points.Note that the extremely high values of some of the coefficients for low b/a values were omitted for better visibility of the overall trends.

Figure 5 . 9 Figure 6 .
Figure 5. Plots of the four coefficients of the fitting equation y = f (x) = k 4 x 4 + k 3 x 3 + k 2 x 2 + k 1 x for the demagnetizing factor N a /4π of the general ellipsoid as a function of b/a.The symbols are the fitted coefficient values as given in Table 2, and the lines are approximate trend lines through the data points.Note that the extremely high values of some of the coefficients for low b/a values were omitted for better visibility of the overall trends.Magnetism 2024, 4, FOR PEER REVIEW 9

Figure 6 .
Figure 6.The same as Figure 5 but for the demagnetizing factor N b /4π.