Neutron imaging for magnetization inside an operating inductor


HomeHome / Blog / Neutron imaging for magnetization inside an operating inductor

Aug 30, 2023

Neutron imaging for magnetization inside an operating inductor

Scientific Reports volume 13,

Scientific Reports volume 13, Article number: 9184 (2023) Cite this article

1 Altmetric

Metrics details

Magnetic components are key parts of energy conversion systems, such as electric generators, motors, power electric devices, and magnetic refrigerators. Toroidal inductors with magnetic ring cores can be found inside such electric devices that are used daily. For such inductors, magnetization vector M is believed to circulate with/without distribution inside magnetic cores as electric power was used in the late nineteenth century. Nevertheless, notably, the distribution of M has never been directly verified. Herein, we measured a map of polarized neutron transmission spectra for a ferrite ring core assembled on a familiar inductor device. The results showed that M circulates inside the ring core with a ferrimagnetic spin order when power is supplied to the coil. In other words, this method enables the multiscale operando imaging of magnetic states, allowing us to evaluate the novel architectures of high-performance energy conversion systems using magnetic components with complex magnetic states.

Magnetic components are key parts of energy conversion systems, such as electric generators, motors, power electric devices, and magnetic refrigerators. Therefore, they have been a mainstay of modern society since the late nineteenth century1. For example, toroidal inductors with ferrite cores and copper coils (Fig. 1a) are found in various electric devices in everyday life. According to Ampère's circuital law, a circumferential magnetic field H is generated inside an inductor coil when an electric current is supplied to the inductor1. Ampère's circuital law predicts that the amplitude of H in the inner circumference is 1.5 times greater than that in the outer one because of the difference in perimeters (Fig. 1b). If magnetization M is induced in a direction parallel to H and its magnitude is proportional to H, M also goes around circumferentially inside the ferrite core with 1.5 times different magnitudes between the inner and outer sides. Although such simple assumptions are not always valid for actual magnetic components used in energy conversion systems, the distribution of M inside it has never been directly verified. Because of magnetic saturation, nonlinear magnetic responses are frequently expected in a homogeneous large H. Moreover, demagnetizing fields generated at the corners or magnetic anisotropy tilt the direction of M form H in general magnetic components. In other words, actual magnetic components do not satisfy the well-known relationship with the following magnetic flux B = μ0(H + M) = μ0(1 + χ)H, where μ0 is the vacuum permeability and χ is the susceptibility. Thus, the distribution of M is not the same as the distribution of H or B. However, over the centuries, magnetic components have been designed using only information from magnetization curves averaged as a whole (Fig. 1c) because the distribution of M cannot be measured unless disassembling the system.

Toroidal inductor with manganese a zinc ferrite ring core. (a) Photograph prior to the measurement. According to Ampère's circuital law, the amplitude H at the inner and outer peripheries, Hin and Hout, is calculated in (b). (c) Average magnetization is shown as a function of average magnetic fields at 296 K. The square symbols in (c) denote the points for measuring neutron transmission spectra.

Currently, magneto–optical Kerr microscopy or spin-polarized scanning electron microscopy is used to elucidate the fine distribution of M on the bare surface of bulky magnetic components2, whereas the H distribution outside the component can be accurately measured using fluxgate or Hall magnetic sensors. On the other hand, the means to non-destructively observe the distributions inside bulky components assembled in energy conversion systems have not yet been established. For example, H inside the magnetic core was approximated to the measured value using a search coil placed inside holes drilled into the core3. Consequently, we usually infer the internal distributions of H and M by comparing the surface or outside information with electromagnetic simulations4,5. This indirect evaluation worked for current designs using existing simple magnetic components. However, future sustainable societies will require a higher efficiency of energy conversion enabled by highly sophisticated magnetic components, where M or its original spin orientations are designed to be nonparallel and nonproportional to H at multiscale. For example, microscopically, some spins are antiparallel to H orientations (so-called ferrimagnetic) in a permanent magnet of (Nd1−xDyx)2Fe14B6, soft magnet of (Mn1−xZnx)Fe2O47, magnetocaloric material of ErCo28, and spintronic material of GdFeCo9, where the antiparallel spins play an important role in their magnetic performance. The tilted spins in the soft-magnetic phase contribute to increasing the energy product in exchange-spring composite magnets at the mesoscale10. Macroscopically, functionally graded magnetic materials11 and multimaterial components12 are used, where the magnetic properties are designed to vary from one place to another inside a single component unit. The magnetic components must be complex on multiscales. It is difficult to evaluate the internal distributions of H and M (otherwise B and M) in such advanced materials using the information obtained from the surface or outside. The lack of a useful evaluation method is a barrier to the advancement of energy conversion systems.

Neutrons are highly penetrative compared with electrons and X-rays. Additionally, the spins of neutrons exhibit precession motion around B, and some of the neutrons are scattered by electron spins that are microscopic M. Therefore, imaging techniques for the precession of the initially polarized neutrons have been developed recently to visualize the distribution of B13,14,15. Meanwhile, neutron scattering has played a critical role in the fundamental study of complex spin arrangements for a long time, where diffractometry detects scattered neutrons with high accuracy. In this context, the M distribution can be easily visualized using neutron diffractometry. However, neutrons scattered at different positions intersect each other when the scattering angles are not the same (See supplementary information). Alternatively, narrowed incident beams must be scanned for such mapping using diffractometry. Consequently, this method is time consuming. However, when such scattering occurs, the intensity of transmission neutrons decreases. This type of decrease observed in the neutron transmission spectra is known as Bragg edges16,17. Recently, we have demonstrated that spin arrangements can be evaluated by analyzing Bragg edges18, indicating that the simultaneous mapping of internal spin arrangements in a wide area is possible using a pair of parallel large-diameter rectilinearity beams and a two-dimensional (2D) detector because the straight trajectories of transmitted neutrons are not crossed over. Herein, we used neutron transmission spectroscopy to evaluate the potential value for imaging of M and microscopic spin arrangements inside magnetic components. We verified the distribution of M and microscopic spin arrangements inside the inductor (Fig. 1) as the simplest magnetic component suitable for the first experiment.

Figure 2 shows the neutron transmissions Tr of the left part of the test inductor with manganese zinc (Mn–Zn) ferrite ring core as a function of the time of flight (TOF) of neutrons. In this figure, the standby state represents the 0 A current (I) to the coil and the operating state represents the 2 A current, corresponding to the remanent and almost saturated states respectively (Fig. 1c). The polarization vector of the incident neutron P was kept upward P(↑) or rotated downward P(↓) by turning the neutron spin flipper off or on19 (Fig. 2a). It is observed that as TOF increases, Tr gradually decreases with edge-like fine structures. Tr for neutrons with P(↑) and P(↓) differs considerably in the operating state, whereas they are the same in the standby state. In general, magnetic scattering at j-th atom is constructive for nuclear scattering at the same atom when P is antiparallel to m⊥j and destructive for nuclear scattering when P is parallel to m⊥j, where m⊥j = (m⊥xj, m⊥yj, m⊥zj) is the vector projection of the magnetic moment mj = (mxj, myj, mzj) of j-th atom on the plane perpendicular to the scattering vector, q. As stated above, the variation in the intensity of transmission has an inverse relationship with that in the total magnitude of scatterings with variously oriented q at the atoms. Therefore, the result that Tr is smaller for P(↓) than that for P(↑) can be attributed to the state that the sum of mj, i.e., M is primarily directed upward through the transmission path (details are analysed later). Conversely, when Tr for P(↓) exceeds that for P(↑), we can consider the downward direction of M. The contour maps of the difference in Tr for P(↓) to that for P(↑) in each pixel of the detector, where Tr is averaged from 10 to 30 ms TOF is shown in Fig. 3. In the operating state, Tr in the left part of the test inductor decreases, as P changes from P(↑) to P(↓), whereas Tr in the right part of the test inductor increases with the change (Fig. 3b), indicating that when H is generated, M is directed upward in the left part and downward in the right part. However, M in the top and bottom parts of the test inductor seems not direct neither upward nor downward. These results are consistent with the simple model, where M circulates in a clockwise direction inside the ferrite ring core during the operating state. The insets in Fig. 3 show the difference in Tr for P(↓) to that for P(↑) along a horizontal line at the vertical centre on the right (the grey striped region). In the operating state, the difference is almost constant in the ranges where the ring core exists. In brief, the magnitude M is constant inside the ferrite core. For the ring core, the perimeter of the outer circle is 1.5 times longer than that of the inner one. Consequently, the amplitude H at the inner periphery is 1.5 times larger than that at the outer one because of the difference in the wire number density (Fig. 1b). If M is proportional to H, the magnitude of M induced on the inner side of the core must be 1.5 times greater than that on the outer side. Alternatively, when the core is magnetically saturated, the magnitude of M is almost constant over the entire core. The feature observed here seems consistent with the latter case. This is reasonable because the magnetization curve is almost saturated when H is applied (Fig. 1c). The current analysis shows that M circulates peripherally inside the ring core with almost equivalent magnitudes between the inner and outer sides in contrast to the naïve prediction based on the relationship M = μ0χH. Let us examine edge-like fine structures of the obtained spectra to elucidate the microscopic magnetic states next.

Transmissions Tr of the test inductor with manganese zinc ferrite ring core for neutrons with upward/downward polarizations as functions of the time of flight (TOF) of neutrons. (a) shows the schematic chart of the transmission experiment, and (b) show the results obtained in the standby state with a 0 A current I and the operating state with a 2 A current I, where Tr is averaged in the coloured part in the inset.

Contour maps of the difference in Tr for P(↓) to that for P(↑) in each pixel of the detector in (a) the standby state and (b) the operation state. The insets show the difference along a horizontal line at the vertical centre on the right. (the grey striped region).

The neutron transmission is represented by an exponential decay \(Tr = e^{ - A\left( \lambda \right)}\), where A(λ) is the absorbance. For the multicomponent test inductor, A(λ) is expressed as the sum of every attenuation effect17:

where i represents the ferrite core (Fe), copper wire (Cu), and resin coating (Re). The elastic Bragg scattering cross section \(\sigma_{{{\text{Bragg}},{\text{i}}}}^{{{\text{ela}}}}\), elastic diffuse one \(\sigma_{{{\text{diffuse}},{\text{ i}}}}^{{{\text{ela}}}}\), elastic incoherent one \(\sigma_{{{\text{coh}},{\text{ i}}}}^{{{\text{ela}}}}\), inelastic one \(\sigma_{{\text{i}}}^{{{\text{inela}}}}\), and absorption cross section \(\sigma_{{\text{i}}}^{{{\text{abs}}}}\), contribute for each component, ni and ti represent the total number of unit cells in the unit volume and the effective thickness of the i-th component, respectively. Among these terms, only the Bragg scattering for the periodic structures of the nucleus at the position of (xj, yj, zj) and magnetic moments mj form fine structures in the spectra because the Laue condition is satisfied at a specific λ. In polycrystalline materials, such Bragg scattering occurs at various λ that are less than twice the interplanetary spacing, dhkl, of the {hkl} planes. Consequently, we can observe an edge-like structure for absorbance, known as the Bragg edge, at 2dhkl17. For thermal neutrons, the other terms have monotonous dependencies on λ; hence, they are discussed in the supplementary section, and we shall concentrate on \(\sigma_{{{\text{Bragg}},{\text{i}}}}^{{{\text{ela}}}}\) described as follows:

where v0 is the unit cell volume2,17. The complementary error function, March–Dollase orientation distribution function, and Sabine's primary extinction function were used for the resolution function Rhkl, the preferred orientation function Phkl, and primary extinction function Ehkl with crystallite size Rc, respectively. The crystal and magnetic structure factors for nonspin flip and spin-flip scatterings of neutrons, FN(hkl), FMnsf(h′k′l′), and FMsf(h″k″l″) are expressed as follows:

where rm is the factor of magnitude (5.39 fm); Biso is the temperature factor; and oj, bj, and fj are the site occupancy, nucleus scattering length, and magnetic form factor of the j-th atom, respectively2. The negative and positive signs in Eq. (3b) correspond to the cases for neutrons with upward polarized states P(↑) and with downward state P(↓), respectively.

Let us begin the investigation of the edge-like fine structures with analysing the transmission spectra in the standby state using the 0 A current I (Fig. 4). Because the spectrum is invariant when P is rotated from upward P(↑) to downward P(↓) by the neutron spin flipper, macroscopic M does not direct either upward or downward directions. Additionally, we cannot find extra Bragg edges compared with the cubic-spinel structure of ferrite and the FCC structure of copper. Therefore, the magnetic structure is collinear ferrimagnetic, as previously observed for Mn1−xZnxFe2O4 ferrites (x < 0.5): {h′k′l′} = {h″k″l″} = {hkl}20. For the ferrite core, we assume that all Zn ions occupy the tetrahedral interstitial site (A site) of the cubic-spinel structure21 and that each crystallite is randomly oriented: Phkl = 1. Furthermore, Biso for the ferrite and copper is set at 0.4 Å221. The observed absorbance can be well reproduced with the curve calculated using Eqs. (1–3) (Fig. 4), and the parameters are shown in Table 1, where the agreement factor R is 0.3% for the edge width of 0.02 Å (see details in the supplementary section) The distribution ratio of Mn ions of the A site to octahedral (B) site is 0.267–0.124, which can be compared with that in a previous report22. The magnetic moment mB averaged over the ions of Mn0.062Fe0.938 at the B site is 1.6 μB, and it is antiparallel to the magnetic moment mA of − 1.4 μB averaged over the ions of Zn0.498Mn0.267Fe0.235 at the A site. The total magnetic moment for a unit chemical formula is roughly 2 μB, which is almost the same as the value estimated from the bulk magnetization at room temperature (see details in the supplementary section). We can confirm that the microscopic magnetic states in the standby state agree with the previously reported results for bare Mn–Zn ferrites that have not yet been assembled in an inductor device20,21,22,23. In the subsequent section, we will consider how the ferrite core becomes microscopically magnetized when an exciting current is supplied to the copper wires.

Wavelength dependence of the absorbance A(λ) in the standby state using the 0 A current: the black circles show the observed results, and the red curve shows the fitted line. {hkl} shows the diffraction planes of Mn–Zn ferrite and Cu.

As discussed above, there are variations in the transmission at the right and left parts of the test inductor when P was rotated from upward P(↑) to downward P(↓) using the neutron spin flipper in the operating state with the exciting current (Fig. 2). Looking the spectra in more details, variations of some Bragg edge heights are considerable, while those of the other are insignificant. Considering Eqs. 2 and 3, such differences between the edges come from variety of \(F_{{\text{M}}}^{{{\text{nsf}}}} \left( {h^{\prime } k^{\prime } l^{\prime } } \right)\) with positive and negative signs in the term of \(\left( {F_{{\text{N}}} \left( {hkl} \right) + F_{{\text{M}}}^{{{\text{nsp}}}} \left( {h^{\prime } k^{\prime } l^{\prime } } \right)} \right)^{2}\). In other words, we can, in principle, estimate the magnitude and sign of \(m_{ \bot j}^{{{z}}}\) of j-th atom in each crystalline site from the analyses on fine variations of the spectra between P(↑) and P(↓). For simplicity, we shall concentrate on the variations at Bragg edges because q at the edges λ = 2dhkl becomes antiparallel to the incident direction (the backscattering). In the case, m⊥j is just written as (mxj, 0, mzj) for neutrons from the Y axis as q is given by (0, − qy, 0). Consequently, the difference in A(λ) between neutrons with P(↑) and P(↓) is expressed as follows:

at λ = 2dhkl.

In addition to these scatterings, the precession motion around B occurs, as stated earlier. Therefore, we must also consider this effect. In general soft-magnets, M is expected to be parallel to B on average: 〈M〉//H. On the condition, B is expressed as μ0(H + 〈M〉) + μ0ΔM, where ΔM is the local fluctuation of M. If ΔM is negligible, the polarization P starts to process around μ0(H + 〈M〉) after the neutrons enter the ferrite core. During the precession motion, the angle between P and H + 〈M〉 is constant. Consequently, the angle between microscopic mj and P is also preserved because of the nature of M//mj in collinear ferrimagnets. In other words, the above discussion on scattering cross sections is applicable even if the precession motion occurs.

Finally, we shall consider the case that effects of locally fluctuated ΔM are not negligible. As known as the depolarization effect, numerous tiny random rotations around inhomogeneous μ0ΔM cause the reduction of the length of P. In previous studies, the reduction rate per unit transmission path is described as \(\alpha = cR_{f} (\mu_{0} \Delta M)^{2} \lambda^{2}\), where c = 2.15 × 1029 m−4 T−2 is a constant, relating the Larmor precession and Rf is the typical size of the local magnetic inhomogeneities, respectively24. In this analysis, the depolarization effect was considered because it was observed in a Mn–Zn ferrite even when the magnetization is almost saturated25. Because this type of depolarization is dominant, Eq. (4) has been revised.

where the details are discussed in the supplementary section. In the evaluation of the magnitude and sign of \({m}_{\text{z}j}\), Rc estimated in the above analyses is used as Rf because the prior study on polycrystalline Mn–Zn ferrites showed that Rf is almost identical to the grain sizes26. In addition, we set ΔM2 at Ms2 − Mz2, where the saturation magnetization Ms is estimated using the law of approximation to magnetic saturation and Mmean at I = 2 A is used as Mz (see the supplementary section). Because FN(hkl), Phkl, and Ehkl are unalterable by the current, we can use them as estimated in the standby state. We estimate mzj at each pixel by fully analysing the variations in the Bragg edge heights in the {311}, {511}, {440}, and {531} planes of the spinel ferrite, where the spin arrangement is assumed collinear as observed in the standby state because it was reported that the collinear structure holds in Mn1−xZnxFe2O4 ferrites (x < 0.5) at the high magnetic field of 4 MA/m20.

The estimated mzj of the ions at A and B sites, mzA and mzB, are mapped in Fig. 5. The magnitude of mzA is ~ 1 μB and its sign is positive in the right part of the ferrite ring core, indicating that the direction of \({\varvec{m}}_{{\text{A}}}\) of A-site ions in the right part is in the upper hemisphere, whereas \({\varvec{m}}_{{\text{A}}}\) of A-site ions in the left part of the ring core is directed in the lower hemisphere because of the estimated negative sign of \(m_{{\text{zA}}}\) in the part. These results are consistent with the interpretation that the direction of \({\varvec{m}}_{{\text{A}}}\) rotates anticlockwise inside the ferrite ring core during the operating state, against the clockwise circulation of macroscopic M (Fig. 3). Conversely, the sign of mzB is negative in the right parts and positive in the left part. These show that \({\varvec{m}}_{{\text{B}}}\) of B-site ions directs downward and upward in the right and left parts, respectively. In other words, the direction of \({\varvec{m}}_{{\text{B}}}\) appear to rotate clockwise, as seen for macroscopic M. The observed antiparallel relationship between \({\varvec{m}}_{{\text{A}}}\) and \({\varvec{m}}_{{\text{B}}}\) is reasonable because it has been known that \({\varvec{m}}_{{\text{B}}}\) is antiferromagnetically coupled with \({\varvec{m}}_{{\text{A}}}\) on the microscopic scale22. Notably, the result obtained here seemed qualitatively invariant with respect to the detailed conditions of the assumption, although many assumptions were employed as stated above.

Distribution of microscopic magnetic state. (a) shows the Z-component of magnetic moment, \({\varvec{m}}_{{\text{A}}}\) estimated for ions in A site and (b) exhibits Z-component of \({\varvec{m}}_{{\text{B}}}\) in B site, where they are averaged in the region of 5 × 5 pixels.

We described the magnetic states inside an inductor device using a Mn–Zn ferrite ring core because no study has directly validated the internal distribution of magnetization vectors M even in such simple inductors. In the experiment, polarized neutron transmission spectra were measured at each pixel of a 2D time-resolved detector. The results showed that M circulates in the clockwise direction inside the ring core when power is supplied to the coil. Detailed analyses on the edge-like fine structures in the spectra can also clarify the microscopic magnetic state and its distribution. For example, magnetic moments in the A-site orient in the direction opposite to that of macroscopic M. The results obtained here are consistent with the simple prediction for the inductor device using a Mn–Zn ferrite ring core, as expected for over a century. This is the first direct verification of the prediction for the magnetic core and an indication that polarized neutron transmission spectroscopic imaging has potential to image the distribution of the macroscopic magnetic state in highly sophisticated magnetic devices with complex structures. We also hope that material/site selective magnetic imaging derived from the edge-like fine structures would be helpful in evaluating the distributions of microscopic magnetic states not only in a homogeneous ferrimagnet employed for the demonstration here but also in functionally graded materials or in materials unevenly cooled/ pressed during a bulky manufacturing process.

However, there remain many issues in the current stage. In respect of the spatial resolution, it is important to selectively use the neutron transmission spectroscopic imaging and other imaging techniques such as X-ray magnetic imaging according to the size of measuring object. In the generally used energy range, the penetration depth of X-ray on magnetic materials is less than 1 mm. Therefore, we expect the neutron magnetic imaging to measure the object thicker than 1 mm. If a minimum spatial resolution of 1% is essential, it corresponds to 0.01 mm for the parts measuring a few millimetres unit per side. The actual spatial resolution obtained herein, 1 mm or more (see Fig. 3), is apparently insufficient for the purpose. Thus, further progress on detectors with high spatial resolutions is necessary for the experiments on smaller objects. Meanwhile, the maximum size depends on the penetration depth of neutron on the measuring object. Considering the actual transmission was 20–30% for the present inductor with a total thickness of 15 mm (Fig. 2), the upper limit of total thickness is a few tens of millimetres for the magnetic devices (not containing neutron absorbing elements such as B, Cd, and Gd). Therefore, this technique is suitable for evaluating small- to medium-sized power conversion devices. With respect to time resolution, it is presently possible to measure only steady states and their periodic responses that can be stroboscopically imaged, as evidenced by the fact that 36 ks was needed for the present measurement. The observations of time-dependent variations such as aging effects require improving the counting rate of the detector because we could not effectively use the powerful pulsed neutron source in J-PARC to avoid miscounting in the detector (see "Materials and Methods" section). With respect to evaluating microscopic states, further advances in analysis methods for the depolarization effects would be highly required to estimate the magnitude of the magnetic moments with high accuracy. Finally, we must mention that it is inevitable that irradiated objects are more or less radio-activated. Especially, it is difficult to bring back the magnetic devices containing elements with high neutron activations such as Co and Eu even after a year. As discussed here, neutron transmission spectroscopic imaging has many issues to be resolved and it is still inconvenient to apply to the evaluation of a wide variety of novel magnetic components with complex magnetic states at the multiscale; however, we can say that this new method is worth developing further because the evaluation of such complex magnetic states designed by novel architectures will be essential for improving the performance of next generation energy conversion systems.

Manganese zinc (Mn–Zn) ferrite core was commercially supplied by TDK Co. and used as the main component of the test inductor, where the outer and inner diameters of the ring core, thickness, and density are 44.5, 30.0, 13.0 mm, 5000 kg/m3, respectively. Inductively coupled plasma emission analysis shows that the compositional ratios of Mn, Zn, and Fe are 0.13:0.166:0.704, respectively. The M–H loop was measured using a B–H analyzer (SY8219, IWATSU Electric Co). The coil was made by winding polyester resin-coated copper wire with a 0.5-mm diameter on the ferrite core with 480 turns (Fig. 1).

The neutron transmission spectra of the inductor were measured as functions of the TOF at beamline 22 (BL22) RADEN in J-PARC27. A magnetic mirror was used to polarize a 50 × 50 mm neutron beam. Consequently, the polarization P = (I+ − I−)/(I+ + I−) became approximately Poff ~ 0.9, where I+ and I− are the intensities of neutrons in spin-up and spin-down states19. Subsequently, the polarities of neutron spins were reversed by a gradient- RF type neutron spin flipper. Because the efficiency of the spin flipper was very close to one, the polarization changed to Pon ~  − 0.9. These polarized states were maintained using guide magnets until the neutrons entered the inductor. A time-resolved neutron gas electron multiplier n-GEM 2D-detector with an effective detection area of 100 × 100 mm and pixels of size 0.8 × 0.8 mm was used to count the transmitted neutrons at a distance of 18.5 m from the source. The effective peak count rate of the detector was 180 kcps; hence, the incident neutron flux density was reduced to the magnitude of 2 × 102 n/s/mm2 in average to avoid miscounting in the detector. Consequently, the measuring time of each image was 36 ks. The velocity estimated from TOF was used to calculate λ. The experiments were performed in the standby state where no current I (0 A) was supplied to the coil and in the operating state with a current I of 2 A, where a magnetic field ranging from 6.87 to 10.2 kA/m was generated inside the coil (Fig. 1b).

All data needed to evaluate the conclusions in the paper are present in the paper and/or the Supplementary Materials.

DuBois, H. Magnetic Circuit in Theory and Practice (G. LONGMANS, and Co. London, 1896).

Google Scholar

Zhu, Y. Modern Techniques for Characterizing Magnetic Materials (Springer, 2005).

Book Google Scholar

Nakata, T., Fujiwara, K., Takahashi, N., Nakano, M. & Okamoto, N. An improved numerical analysis of flux distributions in anisotropic materials. IEEE Trans. Magn. 30, 3395–3398 (1994).

Article ADS Google Scholar

Sykulski, J. K. Computational Magnetics (Springer, 1995).

Book Google Scholar

Li, L. et al. Distributed magnetic flux density on the cross-section of a transformer core. Electronics 8, 297 (2019).

Article Google Scholar

Coey, J. M. D. Magnetism and Magnetic Materials (Cambridge University Press, 2010).

Google Scholar

Van Groenou, A. B., Bongers, P. F. & Stuyts, A. L. Magnetism, microstrueture and crystal chemistry of spinel ferrites. Mater. Sci. Eng. 3, 317–392 (1968).

Article Google Scholar

Tang, X. et al. Magnetic refrigeration material operating at a full temperature range required for hydrogen liquefaction. Nat. Commun. 13, 1817 (2022).

Article ADS CAS PubMed PubMed Central Google Scholar

Kim, K.-J. et al. Fast domain wall motion in the vicinity of the angular momentum compensation temperature of ferrimagnets. Nat. Mater. 16, 1187–1192 (2017).

Article ADS CAS PubMed Google Scholar

Zeng, H., Li, J., Liu, J. P., Wang, Z. L. & Sun, S. Exchange-coupled nanocomposite magnets by nanoparticle self-assembly. Nature 420, 395–398 (2002).

Article ADS CAS PubMed Google Scholar

Jkasai, S., Namikawa, M. & Hiratani, T. Recent progress of high silicon electrical steel in JFE steel. JFE Tech. Rep. 21, 14 (2016).

Google Scholar

Dall’Olio, S. et al. Novel design of a high efficiency multi-bed active magnetic regenerator heat pump. Int. J. Refrig. 132, 243–254 (2021).

Article CAS Google Scholar

Kardjilov, N. et al. Three-dimensional imaging of magnetic fields with polarized neutrons. Nat. Phys. 4, 399–403 (2008).

Article CAS Google Scholar

Hiroi, K. et al. Magnetic field imaging of a model electric motor using polarized pulsed neutrons at J-PARC/MLF. J. Phys. Conf. Ser. 862, 012008 (2017).

Article Google Scholar

Hilger, A. et al. Tensorial neutron tomography of three-dimensional magnetic vector fields in bulk materials. Nat. Commun. 9, 4023 (2018).

Article ADS CAS PubMed PubMed Central Google Scholar

Fermi, E., Sturm, W. J. & Sachs, R. G. The transmission of slow neutrons through microcrystalline materials Phys. Rev. 71, 589–594 (1947).

CAS Google Scholar

Sato, H., Kamiyama, T. & Kiyanagi, Y. Rietveld-type analysis code for pulsed neutron Bragg-edge transmission imaging and quantitative evaluation of texture and microstructure of a welded α-iron plate. Mater. Trans. 52, 1294–1302 (2011).

Article CAS Google Scholar

Mamiya, H. et al. Magnetic Bragg dip and Bragg edge in neutron transmission spectra of typical spin superstructures. Sci. Rep. 7, 15516 (2017).

Article ADS PubMed PubMed Central Google Scholar

Shinohara, T. et al. Polarization analysis for magnetic field imaging at RADEN in J-PARC/MLF. J. Phys. Conf. Ser. 862, 012025 (2017).

Article Google Scholar

Morrish, A. H. & Clark, P. E. High-field Mössbauer study of manganese–zinc ferrites. Phys. Rev. B 11, 278–286 (1975).

Article ADS CAS Google Scholar

Mamiya, H., Terada, N., Kitazawa, H., Hoshikawa, A. & Ishigaki, T. Structural and magnetic properties of dilute spinel ferrites: neutron diffractometry and magnetometry investigations. J. Magn. 16, 134–139 (2011).

Article Google Scholar

Krupicka, P. "Novak ‘oxide spinels" in Ferromagnetic materials, E. P., Wohlfarth, Ed. (North-Holland, Amsterdam, 1982), vol 3.

König, U., Gros, Y. & Chol, G. Study of the sublattices magnetisations Mn0.6Zn0.4Fe2O4 by neutron diffractometry and mössbauer spectrometry. Phys. Stat. Sol. 33, 811 (1969).

Article ADS Google Scholar

Grigoriev, S. V., Maleyev, S. V., Okorokov, A. I., Eckerlebe, H. & van Dijk, N. H. Critical scattering of polarized neutrons in the invar Fe65Ni35 alloy. Phys. Rev. B 69, 134417 (2004).

Article ADS Google Scholar

Rekveldt, MTh. The magnetization process in hot-pressed ferrites studied by neutron depolarization. J. Phys. Colloq. 38, C1–C23 (1977).

Article Google Scholar

van der Zaag, P. J., Johnson, M. T., Noordermeer, A., Por, P. T. & Rekveldt, MTh. Relation between grain size and domain size in MnZn ferrite studied by neutron depolarisation. J. Magn. Magn. Mater. 99, L1 (1991).

Article Google Scholar

Shinohara, T. et al. The energy-resolved neutron imaging system. RADEN. Rev. Sci. Instrum. 91, 043302 (2020).

Article ADS CAS PubMed Google Scholar

Download references

The experiments at BL22 of J-PARC were performed under user programs (2017A0042 2018A0062.) This work was partially supported by MEXT-Program for Creation of Innovative Core Technology for Power Electronics (Grant Number JPJ009777), JSPS KAKENHI (Grant No.19H04400) and JST-Mirai Program, Japan (Grant No. JPMJMI18A3).

Yojiro Oba

Present address: Toyohashi University of Technology, Toyohashi, 441-8580, Japan

National Institute for Materials Science, Tsukuba, 305-0047, Japan

Hiroaki Mamiya, Noriki Terada & Tadakatsu Ohkubo

Japan Atomic Energy Agency, Tokai, 319-1195, Japan

Yojiro Oba, Kosuke Hiroi & Takenao Shinohara

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

H.M., Y.O., K.H., and T.S. conducted the neutron imaging experiments. H.M. and N.T. significantly contributed to data analysis and interpretation. T.O. supervised the conduct of this study. All authors critically reviewed and revised the manuscript draft and approved the final version for submission.

Correspondence to Hiroaki Mamiya.

The authors declare no competing interests.

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit

Reprints and Permissions

Mamiya, H., Oba, Y., Terada, N. et al. Neutron imaging for magnetization inside an operating inductor. Sci Rep 13, 9184 (2023).

Download citation

Received: 03 February 2023

Accepted: 02 June 2023

Published: 06 June 2023


Anyone you share the following link with will be able to read this content:

Sorry, a shareable link is not currently available for this article.

Provided by the Springer Nature SharedIt content-sharing initiative

By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.