research papers\(\def\hfill{\hskip 5em}\def\hfil{\hskip 3em}\def\eqno#1{\hfil {#1}}\)

Journal logoJOURNAL OF
APPLIED
CRYSTALLOGRAPHY
ISSN: 1600-5767

The linear phase correction of modulation of intensity emerging from zero effort (MIEZE) with magnetic Wollaston prisms

crossmark logo

aNeutron Technologies Division, Oak Ridge National Laboratory, Oak Ridge, TN 37830, USA
*Correspondence e-mail: frankli@ornl.gov

Edited by G. J. McIntyre, Australian Nuclear Science and Technology Organisation, Lucas Heights, Australia (Received 17 May 2021; accepted 10 December 2021)

This study proposes the use of two magnetic Wollaston prisms (MWPs) to correct for the linear Larmor phase aberration of modulation of intensity emerging from zero effort (MIEZE), introduced by the transverse size of the sample. With this approach, the contrast of the intensity modulation can be maximized at any scattering angle of interest such that the same contrast as the direct transmission geometry can be fully recovered. The optimum magnetic fields required for the MWPs depend only on the scattering angle and the frequencies of the radio-frequency flippers, and they are independent of the neutron wavelength and beam divergence, which makes the approach suitable for both pulsed and continuous neutron sources.

1. MIEZE and its Larmor phase aberration

The dynamical properties of materials are often the key to understanding their macroscopic properties. A technique that can reach the time scales with a matching length scale is greatly desired to investigate the dynamics in soft matter materials such as proteins (Rheinstädter et al., 2006[Rheinstädter, M. C., Häußler, W. & Salditt, T. (2006). Phys. Rev. Lett. 97, 048103. ]), polymers (Godfrin et al., 2015[Godfrin, P. D., Hudson, S. D., Hong, K., Porcar, L., Falus, P., Wagner, N. J. & Liu, Y. (2015). Phys. Rev. Lett. 115, 228302. ]) and hard condensed matter materials such as superconductors (Blackburn et al., 2006[Blackburn, E., Hiess, A., Bernhoeft, N., Rheinstädter, M. C., Häußler, W. & Lander, G. H. (2006). Phys. Rev. Lett. 97, 057002.]; Haslbeck et al., 2019[Haslbeck, F., Säubert, S., Seifert, M., Franz, C., Schulz, M., Heinemann, A., Keller, T., Das, P., Thompson, J. D., Bauer, E. D., Pfleiderer, C. & Janoschek, M. (2019). Phys. Rev. B, 99, 014429.]; Keller et al., 2006[Keller, T., Aynajian, P., Habicht, K., Boeri, L., Bose, S. K. & Keimer, B. (2006). Phys. Rev. Lett. 96, 225501.]). The neutron spin echo (NSE) technique satisfies this demand by employing the Larmor labeling of neutron spins in a magnetic field, which decouples the energy resolution of a neutron instrument from the width of the wavelength band and thus provides us with another approach to achieve ultra-high resolution (Mezei, 1972[Mezei, F. (1972). Z. Phys. 255, 146-160.]). To implement NSE, two static magnetic fields are required with one on each side of the sample. The neutron spin will precess continuously inside the two magnetic fields through the sample region, where the neutron spin is flipped and the Larmor phase accumulation will be reversed. Any energy exchange between the neutron and the sample can be encoded as a change in the accumulated neutron Larmor phase. However, certain situations, such as samples that would depolarize the neutron beam or require high magnetic field, would be challenging for NSE. Even though paramagnetic (Murani & Mezei, 1980[Murani, A. P. & Mezei, F. (1980). Neutron Spin Echo, pp. 104-112. Berlin, Heidelberg: Springer.]), ferromagnetic (Mezei et al., 2003[Mezei, F., Pappas, C. & Gutberlet, T. (2003). Neutron Spin Echo Spectroscopy: Basics, Trends and Applications, Lecture Notes in Physics 601, 1st ed. Berlin, Heidelberg: Springer.]) and intensity-modulated NSE (Farago & Mezei, 1986[Farago, B. & Mezei, F. (1986). Phys. B + C, 136, 100-102.]) techniques have been developed, either they cannot be applied due to the depolarizing sample environment or they are too complicated to be routinely operated. All of these difficulties prohibit conventional NSE from accessing the full arena of magnetic phenomena, for example, the spin dynamics in ferromagnetic superconductors (Haslbeck et al., 2019[Haslbeck, F., Säubert, S., Seifert, M., Franz, C., Schulz, M., Heinemann, A., Keller, T., Das, P., Thompson, J. D., Bauer, E. D., Pfleiderer, C. & Janoschek, M. (2019). Phys. Rev. B, 99, 014429.]), the spin-wave fluctuations in ferromagnets (Säubert et al., 2019[Säubert, S., Kindervater, J., Haslbeck, F., Franz, C., Skoulatos, M. & Böni, P. (2019). Phys. Rev. B, 99, 184423.]) and the paramagnetic-to-skyrmion lattice transition in transition metal helimagnets (Kindervater et al., 2019[Kindervater, J., Stasinopoulos, I., Bauer, A., Haslbeck, F. X., Rucker, F., Chacon, A., Mühlbauer, S., Franz, C., Garst, M., Grundler, D. & Pfleiderer, C. (2019). Phys. Rev. X, 9, 041059.]).

An alternative approach to allow for the complicated sample environment that involves high magnetic field and depolarizing sample is to use modulation of intensity emerging from zero effort (MIEZE) (Gähler et al., 1992[Gähler, R., Golub, R. & Keller, T. (1992). Physica B, 180-181, 899-902. ]), which has been demonstrated or routinely operated at several neutron facilities, including Larmor at ISIS (Geerits et al., 2019[Geerits, N., Parnell, S. R., Thijs, M. A., van Well, A. A., Franz, C., Washington, A. L., Raspino, D., Dalgliesh, R. M. & Plomp, J. (2019). Rev. Sci. Instrum. 90, 125101.]), ORNL (Brandl et al., 2012[Brandl, G., Lal, J., Carpenter, J., Crow, L., Robertson, L., Georgii, R., Böni, P. & Bleuel, M. (2012). Nucl. Instrum. Methods Phys. Res. A, 667, 1-4.]; Dadisman et al., 2020[Dadisman, R., Wasilko, D., Kaiser, H., Kuhn, S. J., Buck, Z., Schaeperkoetter, J., Crow, L., Riedel, R., Robertson, L., Jiang, C., Wang, T., Silva, N., Kang, Y., Lee, S.-W., Hong, K. & Li, F. (2020). Rev. Sci. Instrum. 91, 015117. ]; Zhao et al., 2015[Zhao, J., Hamilton, W. A., Lee, S.-W., Robertson, J. L., Crow, L. & Kang, Y. W. (2015). Appl. Phys. Lett. 107, 113508. ]), CMRR (Liu et al., 2020[Liu, B., Wang, Z., Wang, Y., Huang, C. & Sun, G. (2020). J. Instrum. 15, P04007.]), RESEDA of FRM-II (Franz et al., 2019[Franz, C., Säubert, S., Wendl, A., Haslbeck, F. X., Soltwedel, O., Jochum, J. K., Spitz, L., Kindervater, J., Bauer, A., Böni, P. & Pfleiderer, C. (2019). J. Phys. Soc. Jpn, 88, 081002.]) and VIN ROSE at J-PARC (Hino et al., 2013[Hino, M., Oda, T., Kitaguchi, M., Yamada, N. L., Sagehashi, H., Kawabata, Y. & Seto, H. (2013). Phys. Proc. 42, 136-141.]). Depending on the direction of the static magnetic field inside the radio-frequency (RF) flippers with respect to the beam, there are longitudinal and transverse configurations for MIEZE. MIEZE requires two neutron resonance RF spin flippers before the sample and they are operated at different RF frequencies ([\omega_{1,2} = 2\pi f_{1,2}]), as shown in Fig. 1[link]. With such spin flippers, the polarization vector of the neutron beam, which is the average direction of the neutron spin, can be modulated sinusoidally, with the phase determined by the phase angle of the RF signal (Hank et al., 1997[Hank, P., Besenböck, W., Gähler, R. & Köppe, M. (1997). Physica B, 234-236, 1130-1132.]). After the beam has passed through the neutron polarization analyzer, which only transmits the polarization vector along one direction, the polarization modulation can be converted into intensity modulation before the beam interacts with the sample, as shown in Fig. 1[link](a). The intensity modulation takes the form [I(t) = A\cos[2(\omega_{2}-\omega_{1})t]+I_{0}]. Because no spin manipulation or analysis is performed afterwards, the neutron beam is no longer sensitive to the factors that can depolarize it in the sample or sample environment, which allows MIEZE to be used for complicated magnetic sample environments or depolarizing samples. When scattered by the sample quasi-elastically or inelastically, the neutrons' energy may be slightly changed, which can cause a time advance or delay in their arrival at the detector. Therefore, the contrast of the modulation ( A / I0) will be changed, from which the intermediate scattering function of the sample can be deduced (Gähler et al., 1992[Gähler, R., Golub, R. & Keller, T. (1992). Physica B, 180-181, 899-902. ]; Golub et al., 1994[Golub, R., Gähler, R. & Keller, T. (1994). Am. J. Phys. 62, 779-788. ]).

[Figure 1]
Figure 1
(a) The conventional transverse MIEZE setup. (b) The MIEZE setup combined with magnetic Wollaston prisms (MWPs) to correct for the linear Larmor phase from the sample. (c) A schematic of the combined intensity modulation in the time and space domains. (d) The intensity modulation in direct transmission and scattering directions. The lines denote the wavefront, and [2S^{\prime}] is the footprint of the sample projected towards the scattering direction, which denotes the effective source size of the intensity modulation. The frequencies of the RF flippers are [\omega_{1,2} = 2\pi f_{1,2}] and the thickness of the RF flipper is δ. The magnetic fields inside the MWPs are B1,2, respectively. L1,2 and D1,2 are the distances from the RF flippers and MWPs to the sample, respectively. The intensity modulation of the neutron beam is developed from the RF flippers in the time domain and a snapshot of the modulation is shown at the sample and detector position for a parallel beam, for the case of without (a) and with (b) the MWPs. s and h are the positions where the neutrons are scattered and captured at the sample and detector, respectively, which have a full width of 2S and 2H, respectively. w is the transverse offset of the neutron's trajectory at the entrance of the slit located at x = L with a full size of 2W. For demonstration purposes, the dashed trajectories are the ideal case where both the incoming and outgoing neutron beams are parallel to the center trajectory. The solid trajectory shows the actual case in the experiment, which is used in the following calculations. The red dots along the neutron trajectory are the intersection points with various components of the setup, and the corresponding Larmor phase and time at these points are denoted as [\Phi_{i}] and ti, respectively.

However, one of the problems of MIEZE is the variation of neutron time of flight (TOF) from the position of scattering at the sample to the position of capture at the detector. This can directly introduce a Larmor phase variation and hence phase aberration on the detector, especially for large scattering angles. To resolve this problem, the contribution of both the sample and detector geometry to the neutron TOF needs to be investigated. For the detector side, there have been extensive discussions about using post data reduction to reconstruct the contrast of the intensity modulations (Schober et al., 2019[Schober, A., Wendl, A., Haslbeck, F. X., Jochum, J. K., Spitz, L. & Franz, C. (2019). J. Phys. Commun. 3, 103001. ]; Oda et al., 2020[Oda, T., Hino, M., Endo, H., Seto, H. & Kawabata, Y. (2020). Phys. Rev. Appl. 14, 054032.]) by shifting the modulation on each pixel of the detector with a known phase. However, the application of such a reduction method is still limited, especially for a sample with large transverse size, as it cannot differentiate the locations where the neutrons are scattered at the sample position. For the sample side, as stated by Brandl et al. (2011[Brandl, G., Georgii, R., Häußler, W., Mühlbauer, S. & Böni, P. (2011). Nucl. Instrum. Methods Phys. Res. A, 654, 394-398.]) and Weber et al. (2013[Weber, T., Brandl, G. R., Georgii, R., Häußler, W., Weichselbaumer, S. & Böni, P. (2013). Nucl. Instrum. Methods Phys. Res. A, 713, 71-75.]), rotating the sample with respect to the beam can also change the resolution function of the MIEZE because it changes the neutron TOF from the sample to the detector. However, rotating a flat sample is only practical within limits and it will also significantly restrict the choices of sample. It has also been suggested by Martin (2018[Martin, N. (2018). Nucl. Instrum. Methods Phys. Res. A, 882, 11-16.]) that the detector can be simply configured to be parallel to the sample surface such that the neutron TOF from the sample to the detector is equal for all the neutron trajectories. However, this approach provides a complete correction only in the ideal case when the scattered neutron trajectories are parallel to the center line shown in Fig. 1[link](a), and there is a one-to-one correspondence between the scattering location at the sample and the pixel of the capture on the detector. In reality, the neutrons scattered by a specific location of the sample can be captured by any pixel of the detector and vice versa. The effectiveness of having the detector parallel to the sample surface has also been investigated with McStas simulations (Dadisman et al., 2022[Dadisman, R., Ehlers, G. & Li, F. (2022). Rev. Sci. Instrum. 93, 013301. ]); this shows that, to minimize the Larmor phase aberration, the detector needs to be perpendicular to the scattered neutrons instead of the incoming beam. Another approach is to physically rotate the RF flippers (Dadisman et al., 2022[Dadisman, R., Ehlers, G. & Li, F. (2022). Rev. Sci. Instrum. 93, 013301. ]) such that the wavefront of the intensity modulation can be shaped spatially to match the scattering angle without changing the neutron TOF from the sample to the detector, but its performance is sensitive to the beam divergence because the space focusing condition is coupled with the time focusing condition and they cannot be independently tuned to be focused towards different planes.

In this report, I propose to use two magnetic Wollaston prisms (MWPs) in addition to the conventional MIEZE setup such that we can precisely manipulate the Larmor phase of the neutron beam spatially to match the scattering angle of interest. Therefore, the resolution function of MIEZE can be maximized at any scattering angle of interest without any physical rotation of the sample or the RF flippers. In the following demonstration, for simplicity, both the incoming and scattered neutrons are parallel to the beam center line with no wavelength dispersion. The neutron intensity modulation will be taken as waves, which originate from the RF flippers, are scattered by the sample and propagate towards the detector, similar to ripples in water. As shown in Fig. 1[link](a), the wavefront of the incident modulation is perpendicular to the beam direction at the sample position, as well as the detector position in the direct transmission geometry. However for the detector at an angle ([2\theta]), the arrival time of a neutron will vary depending on the position of the scattering event at the sample position (s). The consequent wavefront of the scattered neutrons will be tilted and misaligned with the detector plane. As discussed previously, for a realistic beam, the pixelated intensity modulations could only be partially corrected by post data reduction. Without knowing where the neutrons are scattered from the sample, the TOF variation and hence phase variation caused by the transverse size of the sample cannot be corrected. Also, the phase variation will change with scattering angle, which means that the resolution of MIEZE, especially at large scattering angles, will be greatly limited. Without resolving this problem, it would be challenging to further increase the energy resolution, especially at large scattering angles, unless one could compromise the size of the sample and thus neutron intensity.

MIEZE, with two RF flippers operated at different frequencies [\omega_{1,2}], generates a Larmor phase gradient in the time domain along the x direction, as shown in Fig. 1[link](c). With this phase gradient in the time domain, the neutron spin can pick up a Larmor phase as it moves forward. By the time it is captured at the detector position, the Larmor phase it accumulates is [\Phi(t) = 2\Delta\omega t = 2(\omega_{2}-\omega_{1})t], where t is the clock time of the detector. When a neutron travels in a magnetic field, its spin will pick up a Larmor phase that is proportional to the neutron wavelength and the magnetic field integral along the neutron trajectory. For an MWP constructed using two oppositely directed, triangular cross section magnetic fields separated by a planar interface, a gradient in the magnetic field integral can be generated when neutrons are transmitting through (Li et al., 2014[Li, F., Parnell, S. R., Hamilton, W. A., Maranville, B. B., Wang, T., Semerad, R., Baxter, D. V., Cremer, J. T. & Pynn, R. (2014). Rev. Sci. Instrum. 85, 053303.]). Therefore, a linear Larmor phase gradient and hence modulation along the transverse direction (y) will be generated, as shown in Fig. 1[link](c). If both the RF flippers and the MWPs were employed, the neutron spin would experience a Larmor phase gradient in both the time domain (along the x direction) and the space domain (along the y direction) simultaneously, which is not the case in conventional MIEZE. Therefore, with the addition of the MWPs, the gradient direction of the total Larmor phase can be changed from the x direction to any direction in the xy plane by choosing the right combination of the phase gradients in the two directions (x and y), as shown in Figs. 1[link](b), 1[link](c) and 1[link](d). Since the direction of the phase gradient determines the wavefront of the intensity modulation, it can be adjusted such that the wavefront is always perpendicular to the average neutron trajectory and thus parallel to the detector plane. The phase aberrations can thus be minimized.

2. Larmor phase of MIEZE with magnetic Wollaston prisms

Before showing the total Larmor phase a neutron picks up through the setup in Fig. 1[link](b), I will briefly discuss the Larmor phase contributed by an individual RF flipper and MWP. As given by Golub & Gähler (1987[Golub, R. & Gähler, R. (1987). Phys. Lett. A, 123, 43-48.]), the Larmor phase generated by a neutron RF flipper with an angular frequency of ω is given by [\Phi_{f} = 2\Phi_{\rm rf}(t_{i})+\omega t_{\pi}-\phi_{i}]. [\Phi_{\rm rf} = \omega t_{i}] is the phase of the RF flipper upon the entrance of the neutron at ti. [t_{\pi}] is the time the neutron stays inside the flipper, and it is given by [t_{\pi} = {\delta} / (v\cos\psi)] with ψ being the divergence angle, δ the thickness of the flipper and v the neutron speed. [\phi_{i}] is the phase the neutron spin carries before entering the RF flipper. For a MWP, as given by Li et al. (2014[Li, F., Parnell, S. R., Hamilton, W. A., Maranville, B. B., Wang, T., Semerad, R., Baxter, D. V., Cremer, J. T. & Pynn, R. (2014). Rev. Sci. Instrum. 85, 053303.], 2016[Li, F., Parnell, S. R., Bai, H., Yang, W., Hamilton, W. A., Maranville, B. B., Ashkar, R., Baxter, D. V., Cremer, J. T. & Pynn, R. (2016). J. Appl. Cryst. 49, 55-63.]) and Li & Pynn (2014[Li, F. & Pynn, R. (2014). J. Appl. Cryst. 47, 1849-1854.]), the Larmor phase generated ([\Phi_{m}]) depends on the magnetic field integral along the neutron's trajectory [\Phi_{m}(B,y,\psi) =] [ ({2\gamma} / {v})By[\cot\beta+(\cot\beta)^{2}\psi]], where γ is the neutron gyromagnetic ratio, B is the magnitude of the magnetic field of the MWP, β is the inclination angle of the hypotenuse of the MWP with respect to the beam and y is the transverse coordinate of the intersection point of the beam with the midplane of the MWP, as denoted by the red dots in Fig. 1[link](b). With [\Phi_{m}], the MWP can be taken as a device with no thickness along the beam. With each neutron's transverse coordinates, divergence and velocity known, the Larmor phase picked up inside the MWP can be calculated. In the following calculations, the inclination angle of the MWP is assumed to be 45°. For a neutron trajectory incident upon the sample with an offset of s and divergence of ψ as indicated by the solid black line in Fig. 1[link](b), the intersection points with each device are denoted by the red dots. The coordinates of the intersection in time and space and the corresponding Larmor phase have been calculated and are shown in Table 1[link]; using these values the Larmor phase [\Phi_{i}] can be calculated at each point. By propagating the Larmor phase through all the components in Fig. 1[link], the phase [\Phi(t)] at the detector can be obtained, which is given in equation (1)[link]:

[\eqalignno{\Phi(t) &= \underbrace {2(\omega_{2}-\omega_{1})t}_{1}\,+\,{{1} \over {v}}\Bigg(\underbrace{2(\omega_{1} L_{1}-\omega_{2}L_{2})(1+\psi^{2} )^{1/2}}_{2}&\cr &\quad+(\omega_{1}-\omega_{2})\bigg\{-\delta\sec\psi &\cr &\quad+2[ \underbrace{L_{s}^{2}}_{2}+\underbrace{h^{2}+s^{2}-2s(h\cos 2\theta+L_{s} \sin 2\theta)}_{3} ]^{1/2}\bigg\}&\cr &\quad +2\gamma\bigg[\underbrace{(B_{1}-B_{2})s}_{4}+\underbrace{ (B_{1}D_{1}-B_{2}D_{2})\psi}_{5}\bigg](1+\psi)\Bigg)&\cr &\simeq 2(\omega_{2}-\omega_{1})t+{{2} \over { v}}\Bigg\{\underline{\omega_{1}(L_{1}+L_{s})-\omega_{2}(L_{2}+L_{s})}&\cr &\quad + \underline{\gamma\psi(B_{1}D_{1}-B_{2}D_{2})}&\cr &\quad+s\bigg[\underline{\gamma(B_{1}-B_{2})-(\omega_{1}-\omega_{2}) \sin 2\theta}\bigg]\Bigg\}.&(1)}]

Table 1
The propagation of the Larmor phase in the proposed MIEZE setup

[\Phi_{i}] and ti are, respectively, the Larmor phase of neutron spin and time at the intersection points of the neutron trajectory and the spin manipulation components. Fig. 1[link] provides detailed definitions of the parameters.

Position Time ( ti) Neutron phase ([\Phi_{i}])
Entrance of RF1 t0 [\Phi_{0} = 0]
Exit of RF1 t1 [\displaystyle\Phi_{1} = 2\omega_{1}\left (t_{0}+{{\delta} \over {2v\cos{\psi}}}\right)-\Phi_{0}]
Center of MWP1 t2 [\displaystyle\Phi_{2} = \Phi_{m}(B_{1},y_{1},\psi)+\Phi_{1}]
Center of MWP2 t3 [\Phi_{3} = \Phi_{m}(B_{2},y_{2},\psi)+\Phi_{2}]
Entrance of RF2 t4 [\Phi_{4} = \Phi_{3}]
Exit of RF2 t5 [\displaystyle\Phi_{5} = 2\omega_{2}\left(t_{4}+{{\delta} \over {2v\cos{\psi}}}\right)-\Phi_{4}]
Detector t [\Phi = \Phi_{5}]
[t_{0} = t-({1} / {v})[L_{1}(1+\psi^{2} )^{1/2}+L^{\prime}_{s}]]. [y_{1}\simeq -D_{1}\psi-s]. [y_{2}\simeq -D_{2}\psi-s]. [t_{4} = t-({1} / {v})] × [[L_{2}(1+\psi^{2} )^{1/2}+L^{\prime}_{s}]]. [L^{\prime}_{s}] = [[h^{2}+L_{s}^{2}+s^{2}-2s(h\cos{2\theta}+L_{s}\sin 2\theta )]^{1/2}].

In this equation, Bi denotes the magnitude of the magnetic field inside the two MWPs. Li, Di and Ls are the distances from the centers of the RF flippers, MWPs and detector to the sample, respectively. w, s and h are the transverse displacement of the neutron trajectory at the slit, sample and detector, respectively. In the calculation, the divergence angle can be explicitly defined as [\psi\simeq({w-s}) / {L}], where L is the distance from the slit to the sample position. The terms labeled with different numbers in equation (1)[link] represent the major terms that contribute to the Larmor phase, whereas the terms that are not labeled are higher-order terms and thus are negligible. The term labeled as 1 denotes the intensity modulation produced by the RF flippers. The Larmor phase aberration, due to the dispersion of neutron speed (v), the transverse offset of scattering (s) and the distribution of scattering angles ([2\theta]), is explicitly given as 2 and 3. The Larmor phase generated by the MWP is provided by 4 and 5, which are time independent. Term 1 is the main MIEZE signal of interest. Terms 2 and 3 are the aberration terms contributed from the sample size, scattering angle, neutron speed etc., which need to be minimized by the correction. In a MIEZE setup, the sample size (∼2 cm) and the area of interest on the detector (∼10 × 10 cm) are considerably smaller than the length of the scattering arm ( Ls > 2 m), and also the divergence of the incident beam is small. Therefore, equation (1)[link] can be further expanded in s, h and ψ. With the lowest orders summarized in equation (1)[link], for the speed-dependent terms, I will minimize their contributions by setting each of the underlined combinations to zero, yielding the following:

[{\rm Time\ focusing\colon }\ \omega_{1}(L_{1}+L_{s}) = \omega_{2}(L_{2}+L_ {s}),\eqno(2)]

[{\rm Space\ focusing\colon }\ B_{1}D_{1} = B_{2}D_{2},\eqno(3)]

[{\rm Steering\ condition\colon }\ B_{2}-B_{1} = {{1} \over {\gamma}}(\omega_{2 }-\omega_{1})\sin 2\theta.\eqno(4)]

Equation (2)[link] is termed the time focusing condition or the MIEZE condition. In a conventional MIEZE setup without MWPs, it can remove the speed dependence of the Larmor phase for the directly transmitted neutron beam. Equation (3)[link] can ensure the Larmor phase the neutron spin accumulates through the MWPs is independent of the beam divergence. With such a condition, the wavefront of neutron spin can be tilted with respect to the incident beam direction by the same amount regardless of the divergence angle (ψ). For a given scattering angle (2θ) and RF frequency ([\omega_{1,2}]), to ensure the wavefront can be tilted by the right amount, equation (4)[link] needs to be satisfied. Similar to the phased array radar, equation (4)[link] can ensure the scattered beam is modulated with high contrast towards the detector at a specific angle; I call it the steering condition. Equations (3)[link] and (4)[link] combined yield the optimum field required for the MWPs to correct for the linear phase aberrations caused by the transverse dimension of the sample in a regular MIEZE. In combination with equation (2)[link], the Larmor phase of neutron spin at the detector position would be independent of the neutron speed and beam divergence to the first order, which is the prerequisite condition to maximize the contrast of the MIEZE signal.

To understand the principle of this method more intuitively, note that the introduction of the two MWPs can generate a magnetic field integral gradient along the transverse direction (y) as shown in Fig. 1[link](c), which is perpendicular to the Larmor phase gradient direction generated by the RF flippers (x). With the two independent phase gradients, the combined gradient direction of the total Larmor phase can be steered to any direction of interest. Since the phase gradient direction is always perpendicular to the wavefront of the intensity modulation, it can be tuned such that it is perpendicular to the surface of the detector. In this case, the wavefront will be parallel to the detector and the phase aberrations can be minimized. Additionally, when the scattering angle is small ([2\theta\simeq 0]), no field is required for the MWPs ( B1,2 = 0). As also shown by equation (4)[link], the field required for the MWPs is independent of neutron wavelength, which makes the setup capable of accepting a broad wavelength band at both pulsed and continuous neutron sources.

The employment of MWPs to steer the wavefront of the intensity-modulated neutron beam is very similar to a phased array radar, which can create a beam of radio waves that can be electronically steered to point in different directions without moving the antennas. With such analogy, the approach of rotating the RF flippers is similar to the conventional radar, which rotates physically and steadily to sweep the airspace with a narrow beam (Dadisman et al., 2022[Dadisman, R., Ehlers, G. & Li, F. (2022). Rev. Sci. Instrum. 93, 013301. ]). As illustrated in Fig. 1[link](d), as the scattering angle increases, the footprint of the sample or the effective sample size ([2S^{\prime} = 2S\cos{2\theta}]) will decrease. Hence as the scattering angle is approaching [2\theta] = 90°, the intensity modulations get more `coherent', which means fewer aberrations will be contributed from the transverse size of the sample. Ultimately, with [2\theta ] = 90°, the sample becomes a coherent point source of intensity modulations. These effects will be shown and discussed in the following sections.

3. The phase correction calculations of MIEZE

With the Larmor phase known in equation (1)[link], the time modulations measured by the detector can be calculated by integrating over the beam size (w) and sample size (s) and thus divergence ([w = L\psi+s]). Ultimately, the detector will measure a three-dimensional data set where one dimension denotes the time and the other two denote the pixel coordinates of the detector, and each pixel sees intensity modulations in the time domain. With the data reduction procedure outlined by Schober et al. (2019[Schober, A., Wendl, A., Haslbeck, F. X., Jochum, J. K., Spitz, L. & Franz, C. (2019). J. Phys. Commun. 3, 103001. ]) and Oda et al. (2020[Oda, T., Hino, M., Endo, H., Seto, H. & Kawabata, Y. (2020). Phys. Rev. Appl. 14, 054032.]), the modulations can be shifted by a known phase over all the pixels before they are integrated spatially across the detector along h, yielding a one-dimensional intensity oscillation in the time domain, as given in equation (5)[link]. The contrast of the modulation can therefore be given by the amplitude of the one-dimensional oscillation. A higher contrast will yield a higher sensitivity to a small change in neutron energy. Therefore, the dependence of the contrast as a function of the scattering angle yields the resolution function of the MIEZE setup. To calculate the resolution function, the detector thickness is assumed to be zero:

[I(2\theta,t) = \textstyle\int\limits_{-H}^{H}\int\limits_{-S}^{S}\int\limits_{-W}^{W}\cos\left[\Phi(2 \theta,\omega_{1},B_{1},w,s,h,t)\right]\,{\rm d}w\,{\rm d}s\,{\rm d}h.\eqno(5)]

3.1. The phase aberrations of a regular MIEZE

First, to demonstrate the phase aberration problem of a conventional MIEZE, calculations have been performed by setting the field in the MWPs to zero in equation (1)[link]. In the calculation, the parameters are chosen such that L = 5 m, L1 = 4 m, L2 = 2 m, Ls = 2 m, 2H = 100 mm and δ = 30 mm, which are very close to the parameters of the existing RESEDA instrument at FRM-II (Franz et al., 2019[Franz, C., Säubert, S., Wendl, A., Haslbeck, F. X., Soltwedel, O., Jochum, J. K., Spitz, L., Kindervater, J., Bauer, A., Böni, P. & Pfleiderer, C. (2019). J. Phys. Soc. Jpn, 88, 081002.]) and the VIN ROSE instrument at J-PARC (Hino et al., 2013[Hino, M., Oda, T., Kitaguchi, M., Yamada, N. L., Sagehashi, H., Kawabata, Y. & Seto, H. (2013). Phys. Proc. 42, 136-141.]). For demonstration purposes, a parallel beam is used for this calculation. The contrast of the intensity modulation has been plotted in Fig. 2[link] as a function of scattering angle ([2\theta]) for various frequencies ([\Delta f]) and sample sizes (2S). The peak contrast of the resolution function is always centered at zero scattering angle and the contrast of the modulation decreases rapidly as scattering angle increases, especially for large sample size (2S) and high RF frequency differences ([\Delta f]). For example, for [\Delta f] = 100 kHz and 2S = 10 mm, it is impossible to measure the scattering at a scattering angle of 20° with a vanishing contrast.

[Figure 2]
Figure 2
The contrast of the modulation as a function of the scattering angle [2\theta] when changing the frequency difference between the two RF flippers ([\Delta f]) and the size of the sample (2S). The two RF flippers are maintained in time focusing.

To understand the results intuitively, I have calculated the modulations observed on the detector at a scattering angle of [2\theta] = 20°. Note that the calculations performed in the following sections take a realistic beam with the divergence angle defined as [\psi = ({w-s}) / {L}]. As shown in Fig. 3[link], without the correction provided by the MWPs, the contrast of the modulations in the time domain is low across all the pixels of the detector. Setting the magnetic fields inside the MWPs according to equation (4)[link] allows the contrast of the modulations to be fully recovered on each pixel. Along the transverse direction of the detector, the phase of the modulation is monotonically shifted towards the edges of the detector due to the monotonic increase in the neutron TOF. Because of the `coherent' scattered beam prepared by the MWPs, the phase variation measured by the detector can be corrected by post data reduction as outlined by Schober et al. (2019[Schober, A., Wendl, A., Haslbeck, F. X., Jochum, J. K., Spitz, L. & Franz, C. (2019). J. Phys. Commun. 3, 103001. ]) and Oda et al. (2020[Oda, T., Hino, M., Endo, H., Seto, H. & Kawabata, Y. (2020). Phys. Rev. Appl. 14, 054032.]), which is applied in all the following calculations.

[Figure 3]
Figure 3
The intensity modulations observed on the detector along the transverse direction (a) with no MWPs, (b) with MWPs at a scattering angle of [2\theta] = 20°. In this calculation, λ = 5.5 Å, Δf = 100 kHz, 2S = 20 mm, L = 5 m and 2W = 40 mm.

3.2. The scattering angle and frequency dependence of the phase correction

The resolution function of the MIEZE setup with the MWPs has been calculated and plotted in Fig. 4[link] by scanning the scattering angle and the difference of the magnetic field between the two MWPs, where the time focusing and space focusing are maintained for the RF flippers and MWPs, respectively. As one can see, for any given scattering angle, the contrast of the modulation can be maximized by adjusting the difference of the magnetic fields inside the two MWPs ([\Delta B]). Namely, the resolution function of the MIEZE setup can be optimized by the MWPs such that the highest sensitivity can be fully recovered at any scattering angle of interest. This would allow the MIEZE instrument to reach any scattering angle without compromising the contrast of the intensity modulation or the size of the sample. It also shows that the steering condition obtained for the MWPs from the first-order approximation [equation (4)[link]] agrees well with the optimum field calculated from equation (1)[link]. As the frequency of the modulation increases, the window of solid angle with high-contrast modulations gets narrower. Therefore, for a given detector occupying a finite solid angle, it is critical to apply the post data reduction procedure, especially for high-frequency oscillations. While the contrast of the modulation is high at low scattering angles, as shown in Fig. 4[link](c), it will further increase as the scattering angle increases. As shown in Fig. 1[link](c), this could be explained by the diminishing effective sample size and thus increasing coherence as the scattering angle approaches [2\theta ] = 90°, which has been discussed in the previous sections.

[Figure 4]
Figure 4
The contrast of the modulation as a function of the scattering angle (2θ) and the difference of the magnetic fields between the MWPs ([\Delta B = B_{2}-B_{1}]) for three different frequencies: (a) [\Delta f] = 25 kHz, (b) [\Delta f] = 100 kHz and (c) [\Delta f] = 400 kHz. The dashed black lines denote the optimum setting of the two MWPs directly obtained from the steering condition given in equation (4)[link], which show great agreement. In this calculation, λ = 5.5 Å, 2S = 20 mm, L = 5 m, 2W = 40 mm and 2H = 100 mm, which corresponds to a divergence of 0.34° in full width at half-maximum (FWHM) for the incoming beam.

3.3. Wavelength dependence of the phase correction

It has been shown by equation (4)[link] that the optimum fields of the MWPs are wavelength independent. The applicability of this method at different neutron wavelengths is now calculated, where a continuous spectrum of neutrons is scanned simultaneously. The calculations are performed with two RF frequencies ([\Delta f] = 200, 400 kHz). The scattering angle is fixed to be [2\theta] = 20°. For each setting, the fields of the MWPs are optimized and the neutron wavelength is swept to obtain the resolution function. For comparison, the resolution function of conventional MIEZE in direct transmission geometry has also been calculated and included in Fig. 5[link]. With the correction, the contrast of the modulations at the 20° scattering angle can be fully recovered as in the direct transmission geometry. As the modulation frequency and wavelength increase, the total Larmor phase of neutron spin, as well as its phase aberration, gets higher. Therefore, the peak contrast will slightly drop towards longer wavelength and higher frequency. The ability to maintain a high contrast on a wide-bandwidth neutron spectrum means the correction method proposed here is suitable for accepting a broad wavelength band at both pulsed and continuous neutron sources.

[Figure 5]
Figure 5
The variation of the peak value of the resolution function as a function of wavelength for various frequencies. For each case, the field in the MWP has been optimized according to equation (4)[link]. The parameters for the calculation are 2W = 40 mm, 2H = 10 mm and 2S = 20 mm. For comparison, the calculation of MIEZE in transmission geometry has also been included.

3.4. Divergence dependence of the phase correction

According to equation (3)[link], the magnetic field required for the MWPs is independent of the beam divergence to the first order. In this part, I will calculate and investigate the divergence dependence of the correction efficiency. As shown in Fig. 6[link], when the beam divergence increases, the contrast of the modulations will decrease slightly below 2° but dramatically above 2°. To determine the origin of such an effect, I have also calculated the contrast of the conventional MIEZE setup without the MWPs in the transmission direction using a sample with the same effective size [2S′ = 2S cos (20°)], the results of which are denoted by the dashed horizontal lines in Fig. 6[link]. With the correction provided by the MWPs, the peak contrast of the MIEZE setup at a scattering angle of 20° is the same as for the direct geometry MIEZE for a sample with the same effective size. Basically, a higher beam divergence will introduce a larger variation of a neutron's TOF prior to its interaction with the sample even in the direct transmission geometry, where no phase correction is required. In a scattering geometry, the introduction of the MWPs would correct the phase aberration such that the same contrast as in the transmission geometry can be fully recovered with no additional divergence dependence.

[Figure 6]
Figure 6
The variation of the resolution function as a function of field difference between the two MWPs ([\Delta B]) for various beam divergences. The parameters for the calculation are [2\theta] = 20°, Δf = 100 kHz and 2S = 20 mm. The dashed lines indicate the polarization of the conventional MIEZE setup in direct transmission with no MWPs, where a sample with the same effective size has been used as in the scattering beam [2S′ = 2S cos (20°)]. The beam divergence angle is given in FWHM.

4. Discussion and conclusions

Using two magnetic Wollaston prisms, I have presented an approach to correct for the first-order Larmor phase aberration of MIEZE caused by the transverse size of the sample. With this approach, the resolution function of MIEZE can be modified such that the contrast of the intensity modulation can be maximized at any scattering angle of interest. In contrast to the approach of rotating the RF flippers (Dadisman et al., 2022[Dadisman, R., Ehlers, G. & Li, F. (2022). Rev. Sci. Instrum. 93, 013301. ]), the approach with MWPs can shape the wavefront at the sample position precisely such that the correction efficiency is independent of the beam divergence and neutron wavelength. Such a correction approach with MWPs will be demonstrated in a future experiment. Though this approach allows MIEZE to work at a wide scattering angle, MIEZE is still limited to use at a finite solid angle. As discussed by Weber et al. (2013[Weber, T., Brandl, G. R., Georgii, R., Häußler, W., Weichselbaumer, S. & Böni, P. (2013). Nucl. Instrum. Methods Phys. Res. A, 713, 71-75.]), Brandl et al. (2011[Brandl, G., Georgii, R., Häußler, W., Mühlbauer, S. & Böni, P. (2011). Nucl. Instrum. Methods Phys. Res. A, 654, 394-398.]) and Martin (2018[Martin, N. (2018). Nucl. Instrum. Methods Phys. Res. A, 882, 11-16.]), the geometry of the sample would affect the resolution function of MIEZE and thus the contrast of the modulations. Unfortunately, such limitation still exists for the method proposed here unless one can shape the wavefront to any arbitrary shape to match the geometry of the sample. As the thickness of the sample increases, higher phase aberrations would be expected due to the larger variation of the TOF from the sample to the detector, regardless of the phase correction mechanism provided by the MWPs. To achieve the best efficiency of correction and thus contrast, a flat and thin sample is always preferred for this method.

While the calculation is based on a conventional MIEZE instrument with no additional neutron optics in between the RF flippers, for some MIEZE instruments, special optics are embedded in between the two RF flippers to either increase the neutron flux or minimize the phase aberrations; the performance of the correction methodology proposed here would need to be reevaluated according to the specific design of the instrument. Currently, only transverse magnetic Wollaston prisms have been developed (Li et al., 2014[Li, F., Parnell, S. R., Hamilton, W. A., Maranville, B. B., Wang, T., Semerad, R., Baxter, D. V., Cremer, J. T. & Pynn, R. (2014). Rev. Sci. Instrum. 85, 053303.]; Pynn et al., 2009[Pynn, R., Fitzsimmons, M., Lee, W., Stonaha, P., Shah, V. R., Washington, A., Kirby, B., Majkrzak, C. & Maranville, B. (2009). Physica B, 404, 2582-2584.]), which means it would be easier to correct for the phase aberration of MIEZE in the transverse configuration, such as at the VIN ROSE instrument at J-PARC. Although one could adiabatically change a longitudinal magnetic field to the transverse direction such that the transverse MWP can also be used for longitudinal MIEZE like at the RESEDA instrument, the magnetic fields would need to be carefully designed and simulated to ensure the uniformity of magnetic field integral across the beam. Therefore, for longitudinal MIEZE, the approach of rotating the RF flippers might be preferable. With the magnetic fields in the MWP well confined, the MWPs can be easily coupled to other spin manipulation components, such as the RF flippers or guide fields. If both the MWPs and RF flippers are superconducting (Dadisman et al., 2020[Dadisman, R., Wasilko, D., Kaiser, H., Kuhn, S. J., Buck, Z., Schaeperkoetter, J., Crow, L., Riedel, R., Robertson, L., Jiang, C., Wang, T., Silva, N., Kang, Y., Lee, S.-W., Hong, K. & Li, F. (2020). Rev. Sci. Instrum. 91, 015117. ]), they can share the same vacuum chamber and cryogenics, which means the whole setup would be very compact. With the RF flippers and MWPs combined, in a quantum picture, the two neutron spin states can be separated in both time and space domains. Such a setup could also be potentially used for a space–time interferometer to measure the mesoscopic space and time correlation function G(r, t) (Van Hove, 1954[Van Hove, L. (1954). Phys. Rev. 95, 249-262.]).

Acknowledgements

I would like to express my gratitude and appreciation to Steven Parnell (Delft University of Technology), Roger Pynn, Steve Kuhn (Indiana University Bloomington) and Georg Ehlers (Oak Ridge National Laboratory) for proofreading the manuscript.

Funding information

This work was sponsored by the Laboratory Directed Research and Development Program of Oak Ridge National Laboratory, managed by UT-Battelle, LLC, for the US Department of Energy. This material is based upon work supported by the US Department of Energy, Office of Science, Office of Basic Energy Sciences under contract No. DE-AC05-00OR22725.

References

First citationBlackburn, E., Hiess, A., Bernhoeft, N., Rheinstädter, M. C., Häußler, W. & Lander, G. H. (2006). Phys. Rev. Lett. 97, 057002.  CrossRef PubMed Google Scholar
First citationBrandl, G., Georgii, R., Häußler, W., Mühlbauer, S. & Böni, P. (2011). Nucl. Instrum. Methods Phys. Res. A, 654, 394–398.  CrossRef CAS Google Scholar
First citationBrandl, G., Lal, J., Carpenter, J., Crow, L., Robertson, L., Georgii, R., Böni, P. & Bleuel, M. (2012). Nucl. Instrum. Methods Phys. Res. A, 667, 1–4.  Web of Science CrossRef CAS Google Scholar
First citationDadisman, R., Ehlers, G. & Li, F. (2022). Rev. Sci. Instrum. 93, 013301.   Google Scholar
First citationDadisman, R., Wasilko, D., Kaiser, H., Kuhn, S. J., Buck, Z., Schaeperkoetter, J., Crow, L., Riedel, R., Robertson, L., Jiang, C., Wang, T., Silva, N., Kang, Y., Lee, S.-W., Hong, K. & Li, F. (2020). Rev. Sci. Instrum. 91, 015117.   Google Scholar
First citationFarago, B. & Mezei, F. (1986). Phys. B + C, 136, 100–102.  CrossRef CAS Google Scholar
First citationFranz, C., Säubert, S., Wendl, A., Haslbeck, F. X., Soltwedel, O., Jochum, J. K., Spitz, L., Kindervater, J., Bauer, A., Böni, P. & Pfleiderer, C. (2019). J. Phys. Soc. Jpn, 88, 081002.  Web of Science CrossRef Google Scholar
First citationGähler, R., Golub, R. & Keller, T. (1992). Physica B, 180–181, 899–902.   Google Scholar
First citationGeerits, N., Parnell, S. R., Thijs, M. A., van Well, A. A., Franz, C., Washington, A. L., Raspino, D., Dalgliesh, R. M. & Plomp, J. (2019). Rev. Sci. Instrum. 90, 125101.  Web of Science CrossRef PubMed Google Scholar
First citationGodfrin, P. D., Hudson, S. D., Hong, K., Porcar, L., Falus, P., Wagner, N. J. & Liu, Y. (2015). Phys. Rev. Lett. 115, 228302.   Google Scholar
First citationGolub, R. & Gähler, R. (1987). Phys. Lett. A, 123, 43–48.  CrossRef CAS Web of Science Google Scholar
First citationGolub, R., Gähler, R. & Keller, T. (1994). Am. J. Phys. 62, 779–788.   CrossRef Google Scholar
First citationHank, P., Besenböck, W., Gähler, R. & Köppe, M. (1997). Physica B, 234–236, 1130–1132.  CrossRef CAS Web of Science Google Scholar
First citationHaslbeck, F., Säubert, S., Seifert, M., Franz, C., Schulz, M., Heinemann, A., Keller, T., Das, P., Thompson, J. D., Bauer, E. D., Pfleiderer, C. & Janoschek, M. (2019). Phys. Rev. B, 99, 014429.  CrossRef Google Scholar
First citationHino, M., Oda, T., Kitaguchi, M., Yamada, N. L., Sagehashi, H., Kawabata, Y. & Seto, H. (2013). Phys. Proc. 42, 136–141.  CrossRef CAS Google Scholar
First citationKeller, T., Aynajian, P., Habicht, K., Boeri, L., Bose, S. K. & Keimer, B. (2006). Phys. Rev. Lett. 96, 225501.  CrossRef PubMed Google Scholar
First citationKindervater, J., Stasinopoulos, I., Bauer, A., Haslbeck, F. X., Rucker, F., Chacon, A., Mühlbauer, S., Franz, C., Garst, M., Grundler, D. & Pfleiderer, C. (2019). Phys. Rev. X, 9, 041059.  Google Scholar
First citationLi, F., Parnell, S. R., Bai, H., Yang, W., Hamilton, W. A., Maranville, B. B., Ashkar, R., Baxter, D. V., Cremer, J. T. & Pynn, R. (2016). J. Appl. Cryst. 49, 55–63.  Web of Science CrossRef CAS IUCr Journals Google Scholar
First citationLi, F., Parnell, S. R., Hamilton, W. A., Maranville, B. B., Wang, T., Semerad, R., Baxter, D. V., Cremer, J. T. & Pynn, R. (2014). Rev. Sci. Instrum. 85, 053303.  Web of Science CrossRef PubMed Google Scholar
First citationLi, F. & Pynn, R. (2014). J. Appl. Cryst. 47, 1849–1854.  Web of Science CrossRef CAS IUCr Journals Google Scholar
First citationLiu, B., Wang, Z., Wang, Y., Huang, C. & Sun, G. (2020). J. Instrum. 15, P04007.  Google Scholar
First citationMartin, N. (2018). Nucl. Instrum. Methods Phys. Res. A, 882, 11–16.  Web of Science CrossRef CAS Google Scholar
First citationMezei, F. (1972). Z. Phys. 255, 146–160.  CrossRef CAS Web of Science Google Scholar
First citationMezei, F., Pappas, C. & Gutberlet, T. (2003). Neutron Spin Echo Spectroscopy: Basics, Trends and Applications, Lecture Notes in Physics 601, 1st ed. Berlin, Heidelberg: Springer.  Google Scholar
First citationMurani, A. P. & Mezei, F. (1980). Neutron Spin Echo, pp. 104–112. Berlin, Heidelberg: Springer.  Google Scholar
First citationOda, T., Hino, M., Endo, H., Seto, H. & Kawabata, Y. (2020). Phys. Rev. Appl. 14, 054032.  Web of Science CrossRef Google Scholar
First citationPynn, R., Fitzsimmons, M., Lee, W., Stonaha, P., Shah, V. R., Washington, A., Kirby, B., Majkrzak, C. & Maranville, B. (2009). Physica B, 404, 2582–2584.  CrossRef CAS Google Scholar
First citationRheinstädter, M. C., Häußler, W. & Salditt, T. (2006). Phys. Rev. Lett. 97, 048103.   Google Scholar
First citationSäubert, S., Kindervater, J., Haslbeck, F., Franz, C., Skoulatos, M. & Böni, P. (2019). Phys. Rev. B, 99, 184423.  Google Scholar
First citationSchober, A., Wendl, A., Haslbeck, F. X., Jochum, J. K., Spitz, L. & Franz, C. (2019). J. Phys. Commun. 3, 103001.   Google Scholar
First citationVan Hove, L. (1954). Phys. Rev. 95, 249–262.  CrossRef CAS Web of Science Google Scholar
First citationWeber, T., Brandl, G. R., Georgii, R., Häußler, W., Weichselbaumer, S. & Böni, P. (2013). Nucl. Instrum. Methods Phys. Res. A, 713, 71–75.  CrossRef CAS Google Scholar
First citationZhao, J., Hamilton, W. A., Lee, S.-W., Robertson, J. L., Crow, L. & Kang, Y. W. (2015). Appl. Phys. Lett. 107, 113508.   Google Scholar

This article is published by the International Union of Crystallography. Prior permission is not required to reproduce short quotations, tables and figures from this article, provided the original authors and source are cited. For more information, click here.

Journal logoJOURNAL OF
APPLIED
CRYSTALLOGRAPHY
ISSN: 1600-5767
Follow J. Appl. Cryst.
Sign up for e-alerts
Follow J. Appl. Cryst. on Twitter
Follow us on facebook
Sign up for RSS feeds