Processing math: 100%

4293

An analytic expression for the ultimate intrinsic SNR in a uniform sphere
Hong-Hsi Lee1, Daniel K. Sodickson1, and Riccardo Lattanzi1

1Center for Biomedical Imaging, New York University, New York, NY, United States

Synopsis

Ultimate intrinsic SNR (UISNR) is the theoretically highest SNR for given geometry and electrical properties, independent of the coil design. Here, we introduce an analytic exact expression to calculate the UISNR at the sphere center, enabling to directly analyze the dependence on main magnetic field, sample geometry and electric properties. The analytic expression can approximate the UISNR near the center with < 5% error. This work can enable people without access to the full simulation code to calculate UISNR and use it, for example, as an absolute reference to assess the performance of head coils with spherical phantoms.

Purpose

Ultimate intrinsic signal-to-noise ratio (UISNR) is the theoretically largest SNR for given sample geometry and electrical properties, independent of coil design [1-4]. UISNR can be calculated via a current mode expansion and dyadic Green's functions (DGF) [5,6]. Only few modes contribute to the UISNR in the central region of a uniform sphere, compared to thousands of modes at other voxel positions [4,7]. Based on this previous observation, the aim of this work was to derive an analytic expression for the UISNR ζr0 at the center of a sphere.

Theory and Methods

All equations are defined in Fig.1. The receive sensitivity Bx(r)iBy(r) (Eq.(1)) inside a dielectric sphere of radius a<b, where b is the radius of the surface where the current modes are defined, is determined by the vector wave functions Ml,m and Nl,m (Eq.(2)).

At the sphere center (r0), we found that Ml,m(r0)=0 (l,m) and Nl,m(r0)=0 (l,m), except N1,1, N1,0, and N1,1 (Eq.(3)). Substituting Eq.(3) into Eq.(1), we see that only N1,1 contributes to the receive sensitivity at the center, meaning that the corresponding UISNR (ζr0) is completely determined by one divergence-free current mode. This explains why the corresponding shape of the ideal current patterns [5] does not change with field strength (Fig.2), and indicates that loops are the optimal coil elements to maximize central SNR.

The analytic expression for ζr0 is given in Eq.(4), which explicitly shows the dependence of ζr0 on radius and wave number (electrical properties) of the dielectric sphere.

To validate the analytic solution, we investigated the dependence of ζr0 on main magnetic field strength B0 and compared the results with corresponding DGF simulation results, with expansion order lmax=55. Average brain tissue was used for the dielectric properties of the sphere, ϵr and σ [4]. We also tested the case with conductivity σ equals to 105(Ω1m1) to approximate the lossless condition as in Ref. [4].

We used the analytic expression to study the dependence of ζr0 on the sphere radius (a) and the radius of the surface where the current distribution is defined (b).

Finally, we investigated how the UISNR (ζr) varies with the distance r from the center, either along the x-y plane or z-axis (B0 direction), and determined the range for which the UISNR can be approximated by the value at the center (ζr0).

Results and Discussion

Fig.3 shows the double-logarithmic plots of ζr0vs.B0 for B0=0.5-12 T, replicating the plots in fig.6 of Ref. [4]. The ζr0 based on the analytic solution (solid lines) was consistent with the results in Ref. [4], which used a multipole expansion to calculate UISNR, and also with values calculated with full DGF simulations (data points). Fig.3 shows that ζr0 is approximately linear with respect to B0 at low B0 and increase nonlinearly ~B0n with an exponent n>1 at high B0, compatible to published results [4].

The sensitivity of the central UISNR with respect to changes in sphere radius (aζr0) and current radius (bζr0) is shown in Fig.4. The analytic results (solid lines) agree with those from DGF simulations (data points). The UISNR at the center grows monotonically with sphere radius until reaching a plateau (Fig.4a), except at 11 T, where its value oscillates for small sphere radii, likely due to wavelength effects. Intuitively, since there are no conductive losses in air, ζr0 should not depend on the radius at which the current distribution is defined (b). This is, in fact, confirmed by the analytic expression (Eq.(4)) and by Fig. 4b, where bζr0=0.

Fig.5 shows that using the UISNR at the center to approximate the UISNR for a voxel at an intermediate position (0<r/a<1) on the x-y plane (Fig. 5a), or z-axis (Fig. 5b) yields an error <5% if r/a<10-20%.

While an approximate solution, with very limited validity, had been proposed for the UISNR [8], our analytic solution of ζr0 is exact and consistent to full simulations. Calculation speed is 800 times faster than for full simulations.

Conclusion

We introduced an analytic expression to calculate the UISNR at the center of a dielectric sphere, which is rapid, accurate, and enables to directly analyze the dependence on B0, sample geometry and electric properties. The analytic expression can approximate the UISNR near the center with an acceptable error: for a typical spherical phantom with radius = 8.5cm, the UISNR would be valid within a central volume of ~(2cm)3. This work can enable people without access to the full simulation code to calculate UISNR and use it, for example, as an absolute reference to assess the performance of head coils with spherical phantoms [9].

Acknowledgements

This work was performed under the rubric of the Center for Advanced Imaging Innovation and Research (CAI2R, www.cai2r.net), a NIBIB Biomedical Technology Resource Center (NIH P41 EB017183).

References

[1] Ocali, O. & Atalar, E. (1998) MRM 39,462–473.

[2] Schnell, W., et al. (2000) IEEE Trans. Antennas. Propag. 48(3),418-428

[3] Ohliger, M.A., et al. (2003) MRM 50,1018-1030

[4] Wiesinger, F., et al. (2004) MRM, 52:376-390.

[5] Lattanzi, R. & Sodickson, D.K. (2012) MRM, 68:286-304.

[6] Vaidya, M.V., et al. (2014) Concepts Magn. Reson. B Magn. Reson. Eng., 44(3):53-65.

[7] Lattanzi, R., et al. (2012) Proc. ISMRM, 20, 2817.

[8] Kopanoglu, E., et al. (2011) MRM, 66:846-858.

[9] Lattanzi, R., et al. (2010) NMR Biomed. 23(2),142-151.

Figures

Fig.1. Analytic solution for the UISNR at the center of a uniform dielectric sphere (ζr0). Ml,m and Nl,m are vector wave functions; k0 and kin are complex wave numbers outside and inside the sphere, respectively, ω is the Larmor frequency, μ0 and ϵ0 are magnetic permeability and electric permittivity in free space, ϵr and σ are dielectric constant and conductivity of the dielectric sphere, jl is a spherical Bessel function of order l, and Yml is the spherical harmonic function of order (l,m).

Fig.2. A time snapshot of the ideal current patterns associated with the UISNR at the center of a sphere, for B0 = 1.5 T, 3 T, 7 T, and 11 T. Ideal current are the same in all cases, except for a phase shift [5], which is clearly visible in the 2D “unwrapped” view of the spherical surface (bottom row).

Fig.3. UISNR at the sphere center (ζr0) as a function of main magnetic field strength B0, calculated by analytic solution (solid lines) and DGF simulations (data points). Dielectric sphere radius a was set to [5,15,25]cm. Surface current radius b was set to 27.5cm. The dashed line is a reference line with a slope n=1. Dielectric properties of the sphere were set to (a) average brain tissue, or (b) approximately lossless conditions (σ=105Ω1m1)[4]. UISNR (ζr0)wasnormalizedby(FOV)^3$$$ to validate plots against those in fig.6 of previous work that calculated UISNR with multipole expansions [4]. The two results are identical.

Fig.4. (a) aζr0 with a fixed b = 27.5 cm and (b) bζr0 with a fixed a = 5 cm based on analytic solutions (solid lines) and DGF simulations (data points). The dielectric properties were set to values in the average brain tissue [4]. Corresponding analytic equations are shown at the bottom.

Fig.5. The error ζrζr0ζr0 associated with using the analytic solution for the central voxel to estimate UISNR at positions away from the center. The error is shown as a function of the normalized distance r/a from the sphere center. The dielectric properties were set to values in the average brain tissue. The dielectric sphere radius a was set to 8.5 cm, and the surface current radius b was set to 9.5 cm. The normalized distance r/a varies either along the (a) x-y plane or (b) z-axis. The solid line indicates an error of 5%.

Proc. Intl. Soc. Mag. Reson. Med. 25 (2017)
4293