SciELO - Scientific Electronic Library Online

 
vol.68 número3Hydrophobization of paper intended for packagingBreak-up effect of the weakly bound 6Li ions scattered by 208Pb target índice de autoresíndice de materiabúsqueda de artículos
Home Pagelista alfabética de revistas  

Servicios Personalizados

Revista

Articulo

Indicadores

Links relacionados

  • No hay artículos similaresSimilares en SciELO

Compartir


Revista mexicana de física

versión impresa ISSN 0035-001X

Rev. mex. fis. vol.68 no.3 México may./jun. 2022  Epub 14-Abr-2023

https://doi.org/10.31349/revmexfis.68.031101 

Research

Medical Physics

Proton wave function in a water molecule: breakdown of degeneracy caused by interactions with the magnetic field of a magnetic resonance imaging device

C. H. Zepeda Fernándeza  b  * 

J. L. Aguilar Cuevasb 

E. Moreno Barbosab 

aCátedra CONACyT, 03940, CDMX Mexico

bFacultad de Ciencias Físico Matemáticas, Benemérita Universidad Autónoma de Puebla, Av. San Claudio y 18 Sur, Ciudad Universitaria 72570, Puebla, Mexico.


Abstract

The concept of a Magnetic Resonance Imaging (MRI) device is based on the emission of radio waves produced by the protons of the hydrogen atoms in water molecules when placed in a constant magnetic field after they interact with a pulsed radio frequency (RF) current. When the RF field is turned on, the protons are brought to a spin excited state. When the RF field is turned off, the MRI sensors are able to detect the energy released as the protons realign their spins with the magnetic field. In this work we provide a simple model to describe the basic physical mechanism responsible for the operation of MRI devices. We model the water molecule in terms of a central force problem, where the protons move around the (unstructured) doubly negatively charged oxygen atom. First, we employ an analytical treatment to obtain the system’s wave function as well as its energy levels, which we show are degenerate. Next, the energy levels from the water molecule are studied in the presence of a uniform external magnetic field. As a result, they get shifted and the degeneration is lifted. We provide numerical results for a magnetic field strength commonly used in MRI devices.

Keywords: Magnetic resonance imaging; water molecule; energy levels; quantum mechanics

1 Introduction

Magnetic Resonance Imaging (MRI) is a powerful technique to obtain three dimensional images of the structure and composition of human body parts [1,2,3]. The technique is not invasive, in the sense that it does not require the use of radiation. It is constantly evolving as the procedure to obtain the images improves [4,5,6,7]. An alternative technique to obtain three dimensional images is the Computed Tomography (CT). This technique is however invasive as it requires using x-rays [8,9,10] that can cause damage to body tissues [11,12,13]. In addition of not being invasive, MRI is preferred over CT because it can produce clearer images, particularly for the brain, the nerves, the muscles, the spinal cord, among other body parts.

The MRI main components are the magnet and the radio frequency coils (RFC) [14]. MRI works by employing a strong magnetic field, usually in the range of 0.5-7 T, namely, six orders of magnitude larger than the magnetic field strength on the surface of the Earth. For this reason, patients must inform to MRI operators of any metallic implant they may have, as these can be affected by the strong magnetic field. A typical MRI procedure consists of laying down the patient within the device that contains the magnet. The magnetic field, interacts with the protons in the hydrogen atoms that make up the water molecules, abundantly present in the human body. The interaction is through the proton spin that can be oriented either parallel or anti parallel with respect to the magnetic field direction. A RF field produced by the RFC is used to reverse the orientation of the spin. When this RF is turned off, the proton spin returns to their original orientation producing radio waves that are detected by the MRI device to produce the images [14]. The human body is of course made out of different molecules and atoms and not only water. Moreover, the water molecule contains also protons and neutrons within the oxygen nucleus as well as electrons, all of which are subject to interactions with the magnetic field. It is therefore useful to describe why the most relevant interaction of the MRI magnetic field is the one with the protons in the hydrogen that makes up the water molecule. In this work we use a simple model to describe the basis of the interaction between the MRI magnetic field and the protons in the water molecule. The work is organized as follows: In Sec. 2 we first explain the reason to consider the spin of the protons in the water molecule to be the source of the main interaction for the MRI technique. In Sec.3 we model the water molecule in terms of a central force problem and find the solution of the Schrödinger equation describing the protons wave function in the absence of a magnetic field. In Sec. 4 we introduce a constant magnetic field and find the breakdown of degeneration. In particular, we compute the shift of the energy levels with respect to the case in the absence of the magnetic field. We also compute the energy eigenvalues for typical field strengths employed in MIR devices. The interaction of the proton’s spin with the magnetic field is not taken into account for these calculations, then, only the interaction of the proton’s angular momentum with the magnetic field is considered, leaving the interaction with the spin for a future work. We finally discuss our results and conclude in Sec.5.

2 Basic principles of the MRI technique

The human body is mainly composed of water. 80% of human tissues consist of water molecules [15] which in turn are made of two hydrogen and one oxygen atoms. Let us first consider the protons in the hydrogen atoms. This is the principal assumption used for the physics in the MRI [16,17,18]. Under normal circumstances, the orientation of the protons spins is arbitrary. When a human body is placed in a constant magnetic field, some proton spins in the water molecules become aligned whereas some others become anti-aligned with respect to the direction of the magnetic field. The latter correspond to excited energy states. If a RF is turned on, the effect is to reverse the spin alignment with respect to the magnetic field: parallel becomes anti-parallel and vice versa. When the RF is turned off and the proton spins return to the ground state, they release the excess energy in the form of electromagnetic radiation.

Recall that the proton’s Bohr magneton is given by

μNp=e2mp=3.14×10-8eV/T, (1)

and that the proton’s Larmor frequency -a quantity which depends of the magnetic field strenght B- is given by

νp=2μpBh, (2)

where μp=2.79 μNp is the proton’s magnetic moment, that depends on the sum of the magnetic moments of the quarks that make up the proton [19]. The above means that, νp=42.595B MHz/T. Using a typical value for the field strength in a MRI device of B = 7 T, one obtains that the protons emit radiation with a frequency νp=297.21 MHz. In general, for typical values of the field strength used in these devices, this radiation is in the radio wave part of the spectrum.

The water molecule is also composed of electrons. Their corresponding Larmor frequency is given by νe=27.93B GHz/T. Using B = 7 T, one obtains that νe=195.53 GHz. In general, for the typical field strengths used in MRI devices, νe is in the GHz range of the spectrum, namely in microwaves (MW). The other nucleons (including those from the oxygen atom) could also be subject to magnetic resonance [20]. However, they are mostly spin-paired (their total spin is zero). This means in particular that they are not subject to a significant interaction with the magnetic field. Since the RFC is tuned to detect signals in the radio part of the spectrum, the MW radiation is signaled out. As a consequence, the magnetic field in MRI devices only detects transitions produced by the spin flip of protons in the water molecules.

After having spelled out these general remarks, we now proceed to study the proton wave function and energy levels when these are taken as being part of a water molecule. For simplicity we refer to these as the protons.

3 Model for the water molecule

Finding analytical solutions of Schrödinger’s equation for a system of many interacting particles (MIP) is nearly impossible. In most cases, this can be achieved only after resorting to some simplifying assumptions. The case of a water molecule is an example of a MIP system that has attracted considerable attention since long ago to find solutions using either suitable approximations [21, 22, 23, 24] or numerical methods [25,26].

As a first approximation, the water molecule is commonly described as having a di-polar structure, whereby the two protons of the hydrogen atoms represent the positive side, whereas the oxygen atom together with the two other electrons represent the negative side. Due to electrostatic repulsion, these three components form a triangle [27]. It is known that the distance d p between the protons of the two hydrogen atoms is, roughly speaking, constant, namely d p ∼ 151.05 pm. On the same footing, the distance d o between the protons and the oxygen nucleus is about d o ∼ 95.60 pm. One can then consider that the protons move as if they were a rigid body, with a fixed distance between them. The electrons of the two hydrogen atoms have a covalent bond with the oxygen atom. Therefore, the oxygen atom can be regarded as being a negative, doubly charged system that provides a central force for the motion of the two protons while its internal structure can be ignored, for the purpose of this study. The problem is then reduced to considering the motion of the two-proton system around a charged center with a charge Q = -2e.

3.1 Hamiltonian

The origin of the reference system is taken as the position of the charge that produces the central force (oxygen nucleus). As a first approximation the protons are considered as distinguishable. Since the distances between the protons and the oxygen nucleus are constant, the potential energy is also constant, yielding a system’s Hamiltonian given by

H^0=H^1+H^2=L^122mdo2+V^1(do,dp)+L^222mdo2+V^2(do,dp) (3)

where 1 and 2 label the protons. The potential is given by

V^(do,dp)=V^1(do,dp)+V^2(do,dp)=14πϵ0-2e2do+e22dp+14πϵ0-2e2do+e22dp=14πϵ0-2e2do-2e2do+e2dp=14πϵ0e2dp-4e2do. (4)

Substituting the values for dp and do , we obtain V = -23.4 eV.

3.2 Protons wave function

Putting together all the above considerations, the Schödinger’s equation can be separated to be written as two equations, one for each proton, such that

L^122mdo2+V^1ψ1(θ1,ϕ1)=ϵ1ψ1(θ1,ϕ1), (5)

L^222mdo2+V^2ψ2(θ2,ϕ2)=ϵ2ψ2(θ2,ϕ2), (6)

where L^i2 is the orbital angular momentum operator for the protons, i = 1,2, given by

L^i2=-21sinθiθisinθiθi+1sin2θi2ϕi2 (7)

The total system’s energy is given by E=ϵ1+ϵ2. Since V < 0, then ϵi-V>0 for both i = 1,2. As a consequence, the solutions of Eqs. (5) and (6) represent the eigenfunctions of the orbital angular momentum, which are the spherical harmonics, namely,

ψ1(θ1,ϕ1)=Yl1m1(θ1,ϕ1)=(-1)m12l1+14π(l1-m1)!(l1+m1)!Pl1m1(cosθ1),ψ2(θ2,ϕ2)=Yl2m2(θ2,ϕ2)=(-1)m22l2+14π(l2-m1)!(l2+m2)!Pl2m2(cosθ2) (8)

where Plimi(cosθi) for i = 1,2, are the associated Legendre polynomials. Therefore, the total system’s eigenfunction is

Ψl1,m1;l2,m2(θ1,ϕ1;θ2ϕ2)=Yl1m1(θ1,ϕ1)Yl2m2(θ2,ϕ2), (9)

where mi=-li,-li+1,...,li-1,li; for i = 1,2. As an example, Fig. 1 shows four probability distribution obtained from the corresponding proton wave functions in the water molecule. The energy of the system is given by

El1,l2=22mdo2l1(l1+1)+22mdo2l2(l2+1)+V. (10)

Figure 1 Examples of probability distributions obtained from the corresponding wave functions of protons in the water molecule. The dimensions are arbitrary.  

We notice that for given values of l 1 and l 2 the degeneracy is (2l 1 + 1) ( (2l 2 + 1).

3.3 Wave function for identical protons

To finish the analysis of the wave function, as mentioned, the protons were treated as distinguishable. In order to formulate the model considering the protons as indistinguishable, it is necessary to antisymmetrize the spatial part of the wave function, namely,

Ψl1,m1;l2,m2(θ1,ϕ1;θ2ϕ2)=12Yl1m1(θ1,ϕ1)Yl2m2(θ2,ϕ2)±Yl2m2(θ1,ϕ1)Yl1m1(θ2,ϕ2). (11)

Since the proton spin is s = 1/2, they can be either in a singlet or a triplet state. Since the total wave function must be antisymmetric, when the spatial wave function is symmetric, the spin part must be antisymmetric (singlet) whereas if the spatial wave function is antisymmetric, the spin part is symmetric (triplet).

4 Water molecule in a magnetic field

For this analysis, the interaction of the proton’s spin with the magnetic field is not considered, then, only the interaction of the angular momentum with magnetic field is considered.

To introduce the effects of a magnetic field B, we make the minimal substitution pp-qcA, where A is the magnetic potential and B=×A. Therefore, the Hamiltonian for each proton becomes

H^i=12mp^i-qcA^2;    i=1,2, (12)

resulting in the total system’s Hamiltonian

H^=12mp^1-qcA^2+12mp^2-qcA^2+V=H^0-q2mcA^p^1+p^1A^-q2mcA^p^2+p^2A^+q2mc2A^2 (13)

where H^0 is the Hamiltonian with B = 0. Working in the Coulomb gauge A^=0, one obtains that piA^-A^p^i=0. Finally, using that A^=(1/2)B×r^ and L^=r^×p^, the Hamiltonian is given by

H^=H^0-q2mcBL^1-q2mcBL^2+q2mc2A2. (14)

Since the magnetic field is directed along the z^-axis, Eq. (14) reduces to

H^=H^0-q2mcBL^z1-q2mcBL^z2+q24mc2x12+y12+q24mc2x22+y22 (15)

It is easy to note that H^, H^0, L^z1 and L^z2 commute with each other. Therefore, the common set of eigenfunctions is given by |l1m1l2m2>=Ψl1,m1;l2,m2. Using this basis, one can compute the energy eigenvalues

El1,m1;l2,m2=<l1m1l2m2|H^|l1m1l2m2>=El1,l2-m1Bμp-m2Bμp+q24mc2<l1m1l2m2|(x12+y12)|l1m1l2m2>+q24mc2<l1m1l2m2|(x22+y22)|l1m1l2m2>=22md2l1(l1+1)+l2(l2+1)+V-m1Bμp-m2Bμp+q24mc2<l1m1l2m2|(x12+y12)|l1m1l2m2>+q24mc2<l1m1l2m2|(x22+y22)|l1m1l2m2>, (16)

where, μp=e/2mc. Defining ωp=μpB the energy can be written as

El1,m1;l2,m2=22md2l1(l1+1)+l2(l2+1)+V-m1ωp-m2ωp=El1,l2-m1ωp-m2ωp+q24mc2<l1m1l2m2|(x12+y12)|l1m1l2m2>+q24mc2<l1m1l2m2|(x22+y22)|l1m1l2m2> . (17)

The two last terms of Eq. (17) can be written using spherical coordinates, xi2+yi2=do2-do2cos2θi, for i=1,2; as:

<l1m1l2m2|(xi2+yi2)|l1m1l2m2>=do2(1-<l1m1l2m2|cos2θi|l1m1l2m2>)=do2(1-<limi|cos2θi|limi>). (18)

The cos2θi can be written in terms of the spherical harmonics as 16π/5  Y20(θi,ϕi)+1, then, we obtained an integral of three spherical harmonics (Appendix A), whose result is in terms of the Wigner 3j-symbol as:

<limi|cos2θi|limi>=1cm(2li+1)(2li+1)(5)4π112000112mimi0, (19)

where

112000=215

Then, the value of Eq. (18) is,

<l1m1l2m2|(xi2+yi2)|l1m1l2m2>=do21-23215lili2mimi0+1, (20)

For i = 1,2. Finally, Eq. (17) becomes,

El1,m1;l2,m2=22md2l1(l1+1)+l2(l2+1)+V-m1ωp-m2ωp=El1,l2-m1ωp-m2ωp+q24mc2do21-23215l1l12m1m10+1+q24mc2do21-23215l2l22m2m20+1 (21)

Notice that the degeneracy in the quantum numbers m1 and m2, obtained for the B0, is now lifted. As an example, Fig. 2 shows the first two energy levels comparing the cases when B = 0 and B0.

Figure 2 Left: For the B = 0 case, the energy levels are degenerate in the quantum numbers l1, l2, m1, m2. Right: For the B0 case, the degeneracy with respect to the quantum numbers m1, m2 is lifted. The system is however still degenerate with respect to l1, l2. The notation (l1, l2; m1, m2) represents the quantum numbers associated to the wave function Ψl1,m1;l2,m2.  

In a MRI study, the magnetic field strength can be anywhere in the range of 0.5 T and 7 T. Using Eq. (21), it is simple to compute the proton energy levels for any magnetic field strength. As an example, we take B = 7 T. Table I shows the values of the energy levels of the protons in the water molecule for the ground state and the first and second excited states. The values of the last two terms are not written for the best appreciation of the Table, because they are very small -on the order of 10-15. As a comment, the values of the Wigner 3j-symbol are:

112000=215;  112-1-10=112110=0. (22)

Table I Energy levels for the ground, first and second state of the protons in the water molecule.  

State (l1, l2; m1, m2) Energy (eV))
Ground (0,0;0,0) -23.4
First (0,1;0,-1)=(1,0;-1,0) -23.39545379
(0,1;0,0)=(0,1;0,0) -23.397727
(1,0:1,0)=(0,1;0,1) -23.39545421
Second (1,1;-1,-1) -23.39090758
(1,1;0,-1)=(1,1;-1,0) -23.39090779
(1,1;0,0)=(1,1;-1,1)=(1,1;1,-1) -23.390908
(1,1;0,1)=(1,1;1,0) -23.39090821
(1,1;1,1) -23.39090842

We notice that the degeneracy in the quantum numbers m1 and m2 is removed, as it is illustrated also in Fig. 2.

5 Discussion and conclusions

The working principle behind MRI devices is the interaction of the constant magnetic field with the spin of the protons (hydrogen nuclei) in the water molecule, we worked on this idea to know and understand their energy levels. By using a simple model to describe the water molecule, we formulated and analytically solved the Schödinger’s equation. In this model, since the protons are described as having a constant separation between them and with respect to the oxygen nucleus, they move as a rigid body around the oxygen nucleus, which is regarded as a double negatively charged source of a central force. Notice that the result of this simple description of the interaction between the magnetic field and the protons in the water molecule is similar to the Normal Zeeman effect (S=0), where a lifting of the degeneration is also induced by the magnetic field. The similarity comes from the model of the water molecule that we used, which is the analog to the hydrogen atom system with two particles interacting in a central potential. The wave function for the case of B = 0 is a product of the spherical harmonics, and thus it is described by the four quantum numbers l1, l2, m1 and m2. We have obtained the explicit energy degeneracy for all quantum numbers. This degeneracy is lifted when the water molecule is immersed in a constant magnetic field. However, the degeneracy persists with respect to the quantum numbers l1 and l2. Comparing with the case of the Anomalous Zeeman effect, it is expected that this last degeneracy is lifted when the spin of both protons is also accounted for. The interaction of the magnetic field with the total angular momentum ( J=L+S) is considering for a future work.

With this analysis, we show the energy levels for the protons interacting with a magnetic field hoping these results can be used to improve the techniques to obtain the image from a MRI, having more precision.

Acknowledgments

The authors thank A. Ayala for a thorough reading and suggestions to improve this work.

References

[1] Care at Mayo Clinic. https://www.mayoclinic.org/tests-procedures/mri/about/pac-20384768. [ Links ]

[2] Web MD. https://www.webmd.com/a-to-z-guides/what-is-an-mri. [ Links ]

[3] Medical News Today https://www.medicalnewstoday.com/articles/146309. [ Links ]

[4] J. A. Fessler. Optimization Methods for Magnetic Resonance Image Reconstruction: Key Models and Optimization Algorithms, IEEE Signal Processing Magazine, 37 (2020) 33, https://doi.org/10.1109/MSP.2019.2943645. [ Links ]

[5] F. Knoll et al., Deep Learning Methods for Parallel Magnetic Resonance Image Reconstruction. arXiv:1904.01112v1 [eess.SP]. [ Links ]

[6] X. Qu, et al., Magnetic resonance image reconstruction from undersampled measurements using a patch-based nonlocal operator. Medical Image Analysis 18 (2013) 843, https://doi.org/10.1016/j.media.2013.09.007. [ Links ]

[7] B. Wen, et al., Transform Learning for Magnetic Resonance Image Reconstruction: From Model-Based Learning to Building Neural Networks, IEEE Signal Processing Magazine 37 (2020) 41, https://doi.org/10.1109/MSP.2019.2951469. [ Links ]

[8] Mayo Clinic, https://www.mayoclinic.org/tests-procedures/ct-scan/about/pac-20393675Links ]

[9] National Institute of Biomedical Imaging and Bioengineering, https://www.nibib.nih.gov/science-education/science-topics/computed-tomography-ct. [ Links ]

[10] RaiologyInfo.org For patients, https://www.radiologyinfo.org/en/ctscan. [ Links ]

[11] Medical News Today, https://www.medicalnewstoday.com/articles/219970. [ Links ]

[12] Harvard Health Publishing Harvard Medical School, https://www.health.harvard.edu/cancer/radiation-risk-from-medical-imaging. [ Links ]

[13] U.S. Food and drugs Administration. https://www.fda.gov/radiation-emitting-products/medical-imaging/medical-x-ray-imaging. [ Links ]

[14] S.D. Serai, et al., Components of a magnetic resonance imaging system and their relationship to safety and image quality. Pediatr. Radiol. 51 (2021) 716, https://doi.org/10.1007/s00247-020-04894-9. [ Links ]

[15] J. Zhang, G. Shim, S. M. de Toledo, and E. I. Azzam, The Translationally Controlled Tumor Protein and the Cellular Response to Ionizing Radiation-Induced DNA Damage, in TCTP/tpt1 - Remodeling Signaling from Stem Cell to Disease. Results and Problems in Cell Differentiation, edited by A. Telerman and R. Amson (Springer, Cham, 2017), https://doi.org/10.1007/978-3-319-67591-6_12. [ Links ]

[16] R. A. Pooley, Fundamental Physics of MR Imaging, The AAPM/RSNA physics tutorial for residents 25 (2005) 1087, https://doi.org/10.1148/rg.254055027. [ Links ]

[17] A. Pai, R. Shetty, Y. S. Chowdhury, Magnetic Resonance Imaging Physics. (StatPearls2021) https://www.ncbi.nlm.nih.gov/books/NBK564320/. [ Links ]

[18] V. S. Vassiliou et al., Magnetic resonance imaging: Physics basics for the cardiologist. JRSM Cardiovascular Dis. 7 (2018) 1, https://doi.org/10.1177/2048004018772237. [ Links ]

[19] J.H. Sanders and K. C. Turberfield, The Magnetic Moment of the Proton. I. The Value in Nuclear Magnetons. Proc. R. Soc. Lond. 272 (1963) 79, https://doi.org/10.1098/rspa.1963.0043. [ Links ]

[20] Britannica, The Editors of Encyclopaedia. Magnetic resonance. Encyclopedia Britannica, 25 Jan. 2011, https://www.britannica.com/science/magnetic-resonance . Accessed 14 September 2021. [ Links ]

[21] R. McWeeny and K. A. Ohno, A quantum-mechanical study of eater molecule, Proc. R. Soc. Lond. A 255 (1960) 367, https://doi.org/10.1098/rspa.1960.0072. [ Links ]

[22] J. Chem, Multiconfiguration Wavefunctions for the Water Molecule, J. Chem. Phys. 55 (1971) 1720, https://doi.org/10.1063/1.1676302. [ Links ]

[23] A. Sprague Coolidge, A Quantum Mechanics Treatment of the Water Molecule, Phys. Rev. 42 (1932) 189. [ Links ]

[24] S. Aung, Thesis of Doctor of Philosophy. Part I Approximate Hartree-Fock Wavefunctions One-electron Properties and Electronic Structure of the W ater Molecule. Part II Perturbation-Variational Calculation of the Nuclear Spin-Spin Isotropic Coupling Constant in HD California Institute of Technology Pasadena California (1969). [ Links ]

[25] B. Clark et al., Computing the energy of a water molecule using MultiDeterminants: A simple, efficient algorithm. J. Chem. Phys. 135 (2011) 244105, https://doi.org/10.1063/1.3665391. [ Links ]

[26] M. Johansson, Introduction to Electronic Structure Theory. CSC/PRACE Spring School in Coputational Chemistry 2018. http://iki.fi/mpjohans. [ Links ]

[27] S. S. Zumdahl water Encyclopedia Britannica, (2021). https://www.britannica.com/science/water. [ Links ]

Appendix A

The integral of three spherical harmonics is given in terms of the Wigner 3j-symbol, as:

02π0πYl1m1θ,ϕYl2m2θ,ϕYl3m3θ,ϕsinθdθdϕ=2l1+12l2+12l3+14π×l1l2l3000l1l2l3m1m2m3

For this work the values of l1 = l2 (then m1 = m2), l3 = 2 and m3 = 0 are used.

Received: September 27, 2021; Accepted: October 26, 2021

Creative Commons License This is an open-access article distributed under the terms of the Creative Commons Attribution License