Älykkyystesti: osaatko selittää magneettikuvauksen toimintamekanismin matemaattisesti?
Tuskin.
MRI uses the properties of hydrogen nuclei when they are exposed to a high magnetic field and a radio-frequency field, so it is important to analyze the physical characteristics of the hydrogen atom 11.
The hydrogen atom is the simplest atom since, there is a proton in its core and an electron orbiting it experimenting an attractive Coulomb potential.
V(r)=−e24πϵ0r. (1)
To perform the analysis, the Schrödinger equation indepen- dent of time is used:
HΨ(r)=(p22μ+V(r))Ψ(r)=EΨ(r), (2)
where µ is the reduced mass of the system and p, the momentum operator
μ=mp⋅memp+me, (3)
p= −iℏ▽. (4)
Because there is a central potential V (r), natural coordinates are spherical coordinates, and Eq. (2) becomes:
−ℏ22μr2senθ[senθ∂∂r(r2∂∂r)+∂∂θ(senθ∂∂θ)+1senθ∂2∂ϕ2[Ψ(r)−e24πϵ0rΨ(r)=EΨ(r), (5)
re-arranging the equation above:
r2∂2Ψ∂r2+2r∂Ψ∂r+2μℏ2(re24πϵ0+Er2)Ψ+∂2Ψ∂θ2+cosθsenθ∂Ψ∂θ+1sin2θ∂2Ψ∂ϕ2=0. (6)
To solve Eq. (6) the separation of variables is:
Ψ(r,θ,ϕ)=R(r)Y(θ,ϕ), (7)
replacing the solution proposed in Eq. (6) and dividing by R(r)Y (θ,φ) we obtain:
r2R(r)∂2R(r)∂r2+2rR(r)∂R(r)∂r+2μℏ2(re24πϵ0+Er2)=−1Y(θ,ϕ)∂2Y(θ,ϕ)∂θ2−cosθY(θ,ϕ)senθ∂Y(θ,ϕ)∂θ−1Y(θ,ϕ)sen2θ∂2Y(θ,ϕ)∂ϕ2.
Kommentit (26)
Performing the above changes, Eq. (16) is written as:
(1−x2)d2ydx2−2xdxdy+[l(l+1)−m21−x2]y=0.(1-x2)d2ydx2-2xdxdy+l(l+1)-m21-x2y=0. (17)
Equation (17) is known as Associated Legendre Differential Equation, and also reveals why we chose the constant above to be l(l+1). That is, its solutions are given by the Associated Legendre Polynomials:
Pml(x)=(−1)m2ll!(1−x2)m/2dl+mdxl+m(x2−1)l.Plm(x)=(-1)m2ll!(1-x2)m/2dl+mdxl+m(x2-1)l. (18)
The solution Y (θ,φ) = Θ(θ) · Φ(φ), comprising the angular part of the Eq. (8) is of the form:
Yl,m(θ,ϕ)=NPml(cosθ)eimϕ;Yl,m(θ,ϕ)=NPlm(cosθ)eimϕ; (19)
where N is a normalization constant defined as:
N=(2l+1)!(l−1)!4π(l+1)!−−−−−−−−√,N=(2l+1)!(l-1)!4π(l+1)!, (20)
and thus it ensures that the functions Y l,m (θ,φ) are orthonormal. Equation (19) is known as Spherical Harmonics. Returning to the radial part of Schrödinger Eq. (10), rearranging it to be:
d2R(r)dr2+2rdRdr+[2μℏ2(e24πϵ0r+E)−l(l+1)r2[R(r)=0,d2R(r)dr2+2rdRdr+[2μℏ2(e24πϵ0r+E)-l(l+1)r2[R(r)=0, (21)
whole solution is given by the Associated Laguerre Polynomials, which meet the following normalization condition:
∫0∞e−ρρ2l[L2l+1n+l(ρ)[2ρ2dρ=2n[(n+l)!]3(n−l−1)!;∫0∞e-ρρ2l[Ln+l2l+1(ρ)[2ρ2dρ=2n(n+l)!3n-l-1!; (22)
so the solution of Eq. (8) for the radial part is expressed as:
R(r)=(n−l−1)!2n[(n+l)![3−−−−−−−√×(2na0)3/2e−ρ/2ρlL2l+1n−l−1(ρ),R(r)=(n-l-1)!2n[(n+l)![3×(2na0)3/2e-ρ/2ρlLn-l-12l+1(ρ), (23)
Since we have the solutions of each of the equations depending only on one variable, we proceed to build the exact solution of the Schrödinger equation independent of time for the hydrogen atom,
Ψ(r,θ,ϕ)=(n−l−1)!2n[(n+l)![3−−−−−−−√(2na0)3/2×e−ρ/2ρlL2l+1n−l−1(ρ)Yl,m(θ,ϕ);Ψ(r,θ,ϕ)=(n-l-1)!2n[(n+l)![3(2na0)3/2×e-ρ/2ρlLn-l-12l+1(ρ)Yl,m(θ,ϕ); (24)
explicitly writing in terms of the spherical harmonics:
Ψ(r,θ,ϕ)=(2l+1)!(l−1)!4π(l+1)!−−−−−−−−√(n−l−1)!2n[(n+l)![3−−−−−−−√(2na0)3/2×e−ρ/2ρlL2l+1n−l−1(ρ)Pml(cosθ)eimϕ,Ψ(r,θ,ϕ)=(2l+1)!(l-1)!4π(l+1)!(n-l-1)!2n[(n+l)![3(2na0)3/2×e-ρ/2ρlLn-l-12l+1(ρ)Plm(cosθ)eimϕ, (25)
where the principal, azimuthal, and magnetic (n,l,m) quantum numbers take the following values:
nlm=1,2,3,…=0,1,2,…,n−1=−l,−l+1,…,0,…,l−1,ln=1,2,3,…l=0,1,2,…,n-1m=-l,-l+1,…,0,…,l-1,l
The wave function Ψ(r,θ,φ) by itself has no physical meaning, but the expression |Ψ|2 calculated for a place and a given time is proportional to the probability of finding an electron in that place, at that instant.
3. Nuclear magnetic resonance
Nuclear magnetic resonance is a physical phenomenon associated with the intrinsic angular momentum of the spin and the magnetic properties of atomic nuclei. When a nucleus is placed in a magnetic field an interaction occurs between the magnetic moment of the nucleus and the field, resulting in an energy splitting. By the absorption and emission of photons with the right frequency, transitions between these energy states may occur.
Local magnetic changes produced by the absorption and emission of photons are detected by an antenna that sends the signal to a computer for decoding and image generating.
3.1. Nuclear Spin
The particles that make up the atomic nuclei (protons and neutrons) have the intrinsic quantum mechanical property of spin. In nuclear physics the total angular momentum that the nucleus has is called nuclear spin 18, though the term should not be confused with spin of each nucleon or total spin as the sum of all nucleons, because this is just one of the two contributions to the nuclear spin, the other is the angular momentum of the nucleons. So each nucleon will have a net angular momentum such that:
ji=li+si,ji=li+si, (26)
then the nuclear spin will be:
J=∑iji=∑i(li+si),J=∑iji=∑i(li+si),
(27)
or just:
J=L+S.J=L+S. (28)
Each of these vectors have a similar quantum number. The magnitude of orbital angular momentum satisfies:
L2=l(l+1)ℏ2 con l=0,1,2,…L2=l(l+1)ℏ2 con l=0,1,2,… (29)
The direction of the angular momentum vector L is associated with the quantum number m. For a given quantum orbital momentum number l, there are 2l +1 integral magnetic quantum numbers m ranging from −l to l. Along the z axis the component of angular momentum is given as:
Lz=mlℏ.Lz=mlℏ. (30)
For the vector of spin angular momentum S, its quantum number s can take integer or half-integer values:
S2=s(s+1)ℏ2 con s=0,1/2,1,3/2,…S2=s(s+1)ℏ2 con s=0,1/2,1,3/2,… (31)
and for the z component, we find that sz = msℏ, there are 2s + 1 values of ms. To exhibit the property of MRI, the nucleus must have a non-zero value of s. In medical applications, the proton (1H) is the nucleus of most interest, due to its high natural abundance, however other nuclei have been studied.
The proton is a fermion whose value of s is 1/2, therefore there will be two possible values for Sz, i.e. ±ℏ/2. The eigenfunctions describing the nucleus can be written as |+1/2〉 and |−1/2〉, and since in quantum mechanics every physical observable has an associated operator, we can write an eigenvalue equation to describe the observation of the spin state as:
Sz|ms⟩=ms|ms⟩,Sz|ms⟩=ms|ms⟩, (32)
where Sz is the operator describing the measurement of angular momentum along the z axis.
3.2. Voxel magnetization
To measure the energy of the spin system it is necessary to construct a Hamiltonian operator. The shape of the Hamiltonian can be derived from classical electromagnetism, for the energy of a magnetic moment placed in a magnetic field 15. The nuclei have a magnetic moment μ⃗ μ⃗ which is proportional to the spin angular momentum 4:
μ⃗ =γS,μ⃗=γS, (33)
where γ is a constant of proportionality called gyromagnetic ratio, which for the proton it has a value of 2.675×108 rad/s·T. When this magnetic moment is placed in a magnetic field B its ground state is degenerated and the energy of a particular level will be:
E=−μ⃗ ⋅B;E=-μ⃗⋅B; (34)
combining Eqs. (6), (7) the Hamiltonian is obtained:
H=−ℏγS⋅B,H=-ℏγS⋅B, (35)
and since the field B is oriented in the z direction, the Hamiltonian becomes:
H=−ℏγSzBz;H=-ℏγSzBz; (36)
now using Schrödinger’s equation, the energy eigenstates are found
H|ms⟩=−ℏγBzSz|ms⟩=−ℏγBzms|ms⟩,H|ms⟩=-ℏγBzSz|ms⟩=-ℏγBzms|ms⟩, (37)
and they will be:
E= −ℏγBzmsE= -ℏγBzms (38)
ΔE= E(ms=−1/2)−E(ms=1/2)=ℏγBz.ΔE= Ems=-1/2-Ems=1/2=ℏγBz. (39)
For the proton ms = ±1/2, a transition between the two states represents a change in energy. This is called Zeeman splitting17 and is shown in Fig. 1, where it shows that as the intensity of the magnetic field increases, the energy difference between the states is bigger. The two possible states are more commonly known as “spin−up” and “spin−down”, where the latter has a higher energy state than the “spin−up” state. Transitions between the two states can be induced by absorption or emission of a photon such that:
ΔE=ℏγBz=ℏω.ΔE=ℏγBz=ℏω. (40)
Thus obtaining the Larmor equation which underpins the whole phenomenon of Nuclear Magnetic Resonance.
ω0=γB0.ω0=γB0. (41)
The characteristic frequency ω0 is called Larmor frequency and B0 = Bz. In a real system there is not only an isolated nucleus, but many nuclei which may occupy a particular spin state. This means that our development must be extended to consider an ensemble of spins. To do this we define a Ψ eigenstate, which is a linear combination of possible spin states for a single nucleus as:
∣∣Ψ⟩=∑msams|ms⟩.|Ψ⟩=∑msams|ms⟩. (42)
When performing a measurement on the system, the expected value of the operation on this superposition of states is:
⟨Ψ∣∣Sz∣∣Ψ⟩=∑ms|ams|2ms,⟨Ψ|Sz|Ψ⟩=∑ms|ams|2ms, (43)
where the value |ams|2 represents the probability of finding a single nucleus in the ms state. For the case of a proton, which is the particle of interest, with two spin states we have:
∣∣Ψ⟩=a+12∣∣+12⟩+a−12∣∣−12⟩.|Ψ⟩=a+12|+12⟩+a-12|-12⟩. (44)
The ratio of the population in the two energy states follows the Boltzmann statistics 4,
∣∣a−1/2∣∣2∣∣a+1/2∣∣2=exp(−ΔEkBT)=exp(−ℏγB0kBT),|a-1/2|2|a+1/2|2=exp-ΔEkBT=exp-ℏγB0kBT, (45)
and since kBT ≫ ℏγB0, expanding the first order exponential we have
∣∣a−1/2∣∣2∣∣a+1/2∣∣2≈1−ℏγB0kBT,|a-1/2|2|a+1/2|2≈1-ℏγB0kBT, (46)
and the difference in the number of spins between the spin-up state and the spin-down state will be:
∣∣a+1/2∣∣2−∣∣a−1/2∣∣2≈∣∣a+1/2∣∣2⋅ℏγB0kBT≈ℏγB02kBT.|a+1/2|2-|a-1/2|2≈|a+1/2|2⋅ℏγB0kBT≈ℏγB02kBT. (47)
The magnitude of the magnetization of a voxel containing a nucleus density ρ is obtained by multiplying the density by the difference between the number of spins in both states and the magnetic moment of the nuclei µz = ±(γℏ/2)
M=ργ2ℏ2B04kBT.M=ργ2ℏ2B04kBT. (48)
The magnetization of the volume element also has a direction which is equal to the orientation of the main magnetic field for that reason it should be considered as a vector quantity, as it shown in Fig. 2.
3.3. Excitation of the Magnetization by radio-frequency pulse
A volume element of the sample is defined as a voxel. Having placed a transmitting antenna in the direction of maximum emission to the voxel in the vertical plane and changing the transmission frequency, when we emit the frequency of precession (Larmor frequency), the nuclei are able to absorb energy, i.e. enter resonance.
When the nuclei of the voxel are entering into resonance, the magnetization vector M travels performing a spiral rotational movement relative to the direction of the field B 0. Each nucleus entering into resonance at a specific frequency specified by the Larmor Law, depending on B 0 perceived and biochemical environment in which it is located. Therefore, the radio-frequency emission contains an approximate bandwidth of 100 kHz.
Separation of M with respect to its equilibrium position is determined by the angle of inclination or FLIP angle. Its value depends on the power and time of radio-frequency emission, among other factors.
The process of the radio-frequency emission is of the order of microseconds called radio-frequency pulse and quantified by the value α. A 90° pulse moves the magnetization vector about the xy plane, and a 180° pulse inverts the magnetization respect to its equilibrium position.
After a 90° pulse, the longitudinal component of the magnetization vector is zero, in this position the number of nuclei in the lower energy state equals the number of nuclei in the higher energy state, it is called a state of saturation.
3.4. MRI selectivity
If we have multiple voxels placed under different magnetic fields, we can selectively excite either by simply changing the transmission frequency of the antenna.
A magnetic gradient enables atomic nuclei to perceive a distinct magnetic field. In addition to changes in biochemical environment, one can selectively resonate all nuclei within positions that are being excited by the frequency bandwidth used in the pulse emitter. Therefore all voxels contained in a plane perpendicular to the direction of the gradient plane and whose thickness depends, once defined the value of the gradient, on the bandwidth used on the pulse emitter, will be excited.
Magnetic resonance images are obtained by sending pulses to different spaced time intervals with suitable values, which is a pulse sequence.
3.5. Evolution of magnetization under temporal variation of the magnetic field
Taking a volume element (Voxel) that has a magnetization vector M, the set of equations describing the temporal evolution of said vector are called Bloch equations1 and are defined as follows:
∂Mx∂t= γ(M×B)x−MxT2∂Mx∂t= γM×Bx-MxT2 (49)
∂My∂t= γ(M×B)y−MyT2∂My∂t= γM×By-MyT2 (50)
∂My∂t= γ(M×B)z−Mz−M0T1∂My∂t= γM×Bz-Mz-M0T1 (51)
where γ is the gyromagnetic ratio of the proton and M0 = M is the magnitude of the magnetization in equilibrium, T1 and T2 are time constants that will be discussed later. The above equations can be solved with the appropriate conditions; for example, immediately after when the radio-frequency pulse is turned off:
Mx(t)= [Mx(0)cos(ω0t)+My(0)sen(ω0t)]⋅e−t/T2,Mxt= Mx(0)cos(ω0t)+My(0)sen(ω0t)⋅e-t/T2, (52)
My(t)= [My(0)cos(ω0t)−Mx(0)sen(ω0t)]⋅e−t/T2,Myt= My(0)cos(ω0t)-Mx(0)sen(ω0t)⋅e-t/T2, (53)
Mz(t)= Mz(0)e−t/T1+M0[1−e−t/T1].Mzt= Mz(0)e-t/T1+M01-e-t/T1. (54)
3.6. Nuclear relaxation
Once the radio emission is completed, the magnetization vector, M, returns to its initial position by releasing energy through a process called relaxation. Relaxation occurs because nuclei emit the excess energy absorbed when entering resonance. Relaxation ends when the proportion of nuclei voxel between higher energy states and lower energy states match the Boltzmann equilibrium.
The return to the equilibrium position of M produces magnetic field modifications that may be collected by a receiving antenna, since the magnetic field variations induce an electrical signal called Free Induction Decay (FID) 16, which is a damped sinusoidal signal. The frequency of the sinusoid is the precession frequency imposed by the value of the magnetic field during relaxation.
Given the definition for the electromotive force (fem) and the magnetic flux through a coil:
fem=−dΦdt Φ=∫SB⋅ds.fem=-dΦdt Φ=∫SB⋅ds. (55)
On the other hand, for the magnetization vector M of a sample, there is an associated magnetic field from the current
density
JM(r,t)=▽⃗ ×M(r,t),JM(r,t)=▽⃗×M(r,t), (56)
the vector potential derived from a source is defined as:
A(r)=μ04π∫J(r)|r−r'|dr,'A(r)=μ04π∫J(r)|r-r'|dr,' (57)
with this, the magnetic field can be calculated as B=▽⃗ ×AB=▽⃗×A, and we can write the flux through the coil as:
Φ=∫sB⋅ds=∫(▽⃗ ×A)⋅ds=∮A⋅dl.Φ=∫sB⋅ds=∫▽⃗×A⋅ds=∮A⋅dl. (58)
Substituting the vector potential in the term for the current density and later in the line integral, we have:
ΦM=∮dl⋅[μ04π∫▽⃗ '×M(r')|r−r'|dr'],ΦM=∮dl⋅μ04π∫▽⃗'×M(r')|r-r'|dr', (59)
simplifying the above equation and performing integration by parts:
ΦM=μ04π∫dr'M(r')⋅[▽⃗ ×(∮dl|r−r'|)].ΦM=μ04π∫dr'M(r')⋅▽⃗×∮dl|r-r'|. (60)
The expression in the parentheses of Eq. (60) can be compared with the vector potential A, obtaining an expression for circuits:
A(r')=μ04π∮Idl|r−r'|,A(r')=μ04π∮Idl|r-r'|, (61)
thus the curl of the line integral is indeed the magnetic field per unit current that would be produced by the coil at point r'
Br(r')=B(r')I=▽⃗ '×(μ04π∮dl|r−r'|).Br(r')=B(r')I=▽⃗'×μ04π∮dl|r-r'|. (62)
The magnetic flux through the coil can be written as:
ΦM=∫muestraBr(r')⋅M(r',t)dr'.ΦM=∫muestraBr(r')⋅M(r',t)dr'. (63)
Finally, the electromotive force is expressed as:
fem=−ddt∫muestraBr(r')⋅M(r',t)dr'.fem=-ddt∫muestraBr(r')⋅M(r',t)dr'. (64)
Two voxels placed under different magnetic fields in the moment of relaxation, will have different relaxation frequencies and, therefore, their signals may be differentiated through a frequency analysis such as Fourier analysis. By studying signals of relaxation, one can obtain information on the density (D) of hydrogen nuclei in the voxel and information related with the medium by the parameters T1, T2, and T2 *.
Figure 3 The recovery of Mz to M0 is controlled by an exponential. T1 is the necessary time to recover the 63% of the equilibrium value of the magnetization 4.
The values of T1 increases with the value of the magnetic field.
3.8. Spin-spin relaxation (T2 Relaxation)
We can get information related to the biochemical structure of the medium, by studying the variations over time of the component for the magnetization vector in the vertical plane (x,y) during the transverse relaxation.
When Mx,y is zero it implies that the magnetization vector is aligned on the z-axis with the main magnetic field. Immediately after excitation, part of the spins precess synchronously, these spins have a 0° phase and are said to be in phase, this state is called coherent phase.
The coherent phase is gradually lost as the spins advance and others are delayed on their way precession. In other words the transverse relaxation is the decay of transverse magnetization due to loss of coherence in the spins.
The transverse relaxation differs from longitudinal ralaxation in that the spins do not dissipate energy to the surroundings but, rather, exchange energy with each other.
The evolution of the transverse magnetization over time, until it vanishes, corresponds to a sinusoid with a dampened relaxation frequency generating an exponential decay (see Fig. 4). This exponential decay of the surroundings is regulated by a parameter T2* 2. When we take into consideration all factors that influence the asynchronism of the nuclei, or T2 if neither are influenced by the external magnetic field’s inhomogeneities or local magnetic variations acting permanently on the nuclei. Therefore T2 indicates asynchronism of nuclei of the voxel during relaxation due to spin-spin random influences which depend on the composition and structure of the tissue. T2 is the time that must elapse before the transverse magnetization loses 63% of its value. The time constant T2 is called the Transverse Relaxation Time.
4. CONCLUSIONS
MRI is a technique for visualizing tissues that takes use of the physical phenomenon of nuclear magnetic resonance, which is the union of quantum mechanics with classical electrodynamics, that uses the quantum-mechanical properties of the hydrogen atom to produce high resolution images that help with medical diagnosis. Nowadays, there are many MRI machines used for both medicine and research, where, for the second case it is essential to have knowledge of the facets at a quantum level of the physical process involved. Likewise there is a wide variety of research topics such as digital processing of images acquired during a study, magnetic resonance spectroscopy, designing new radio frequency coils, pulse frequency design, etc., thus the broad field for the development and monitoring of research.
Vierailija kirjoitti:
En mutta fysikaalisesti osaan. Tai ainakin osasin tentissä joskus.
Tuskin osaat oikeasti.
Just meidän kuusivuotiaan kanssa mietittiin samaa ja kyllä hän sen ihan oikein oli esittänyt.
Vierailija kirjoitti:
Just meidän kuusivuotiaan kanssa mietittiin samaa ja kyllä hän sen ihan oikein oli esittänyt.
Noobelin palkintoa pukkaa viimestään seitsemänvuotiaana!
Kun somejulkkis työntää pakaraimplanttinsa kuvauslaitteeseen, niin posket alkavat resonoida ja kohta ulos putkahtaa kuva. Sitä on MRI.
Ei mittaa puhtaasti älyä, vaan ensisijaisesti perehtyneisyyttä fysiikkaan. Myönnän kyllä, että aivan tomppeli ei noita yhtälöitä ottaisi lukemallakaan haltuun. Itsekään en suoralta kädeltä pääse kiinni vaikka Schrödingeriä ja Boltzmannia on tullut pyöriteltyä yliopistossa.
Vierailija kirjoitti:
Just meidän kuusivuotiaan kanssa mietittiin samaa ja kyllä hän sen ihan oikein oli esittänyt.
Ei kuusivuotias hyödy tuollaisesta taidosta mitenkään jos ei samalla osaa arvostaa sitä ettei hänen äidillään ole silikoneja.
Vierailija kirjoitti:
Just meidän kuusivuotiaan kanssa mietittiin samaa ja kyllä hän sen ihan oikein oli esittänyt.
Sellaista se on noiden erityislasten kanssa. Kehitys laahaa selkeästi jäljessä. Normaali lapsihan osaa kertoa tuon viimeistään kun pääsee synnäriltä kotiin.
Tuo matematiikka on hankalahkoa, mutta ei se toimintaperiaate mitenkään monimutkainen ole. Ydinten spin suuntautuu ulkoisen magneettikentän mukaan, ja kun kenttää tuodaan voimakas poikkeutusimpulssi, atomien spin käy "resonoimaan", eli palautuu ikään kuin vaappuen siihen alkutilaan. Tämän perusteella voidaan sitten nähdä protonien määrän erojen perusteella eri kudokset.
We can separate the Eq. (8) into two with a separation constant. For reasons that will eventually become clear, we will write this separation constant as l(l + 1), i.e.
1Y(θ,ϕ)∂2Y(θ,ϕ)∂θ2+cosθY(θ,θ)senθ∂Y(θ,ϕ)∂θ +1Y(θ,ϕ)sen2θ∂2Y(θ,ϕ)∂ϕ2=−l(l+1),1Y(θ,ϕ)∂2Y(θ,ϕ)∂θ2+cosθY(θ,θ)senθ∂Y(θ,ϕ)∂θ +1Y(θ,ϕ)sen2θ∂2Y(θ,ϕ)∂ϕ2=-l(l+1), (9)
r2R(r)∂2R(r)∂r2+2rR(r)∂R(r)∂r+2μℏ2(re24πϵ0+Er2)=l(l+1).r2R(r)∂2R(r)∂r2+2rR(r)∂R(r)∂r+2μℏ2re24πϵ0+Er2=l(l+1). (10)
For the solution of Eq. (9) the separation of variables method is used again, this time making Y (θ,φ) = Θ(θ) · Φ(φ), to get:
1Θ(θ)d2Θdθ2+cosθΘ(θ)senθdΘdθ+1Φ(ϕ)sen2θd2Φdϕ2=−l(l+1);1Θ(θ)d2Θdθ2+cosθΘ(θ)senθdΘdθ+1Φ(ϕ)sen2θd2Φdϕ2=-l(l+1); (11)
multiplying the above equation by sin2 θ: sen
sen2θΘ(θ)d2Θdθ2+senθcosθΘ(θ)dΘdθ+1Φ(ϕ)d2Φdϕ2=−l(l+1)sen2θ,sen2θΘ(θ)d2Θdθ2+senθcosθΘ(θ)dΘdθ+1Φ(ϕ)d2Φdϕ2=-l(l+1)sen2θ, (12)
separating the azimuth part,
1Φ(ϕ)d2Φdϕ2=−m.1Φ(ϕ)d2Φdϕ2=-m. (13)
The solution to Eq. (13) is given as:
Φ(ϕ)=Ae−imϕ+Beimϕ,Φ(ϕ)=Ae-imϕ+Beimϕ, (14)
where A and B are constants. Equation (12) is as follows:
sen2θΘ(θ)d2Θdθ2+senθcosθΘ(θ)dΘdθ−m2=−l(l+1)sen2θ,sen2θΘ(θ)d2Θdθ2+senθcosθΘ(θ)dΘdθ-m2=-l(l+1)sen2θ, (15)
multiplying the equation by Θ(θ) we get:
sen2θd2Θdθ2+senθcosθdΘdθ+[l(l+1)sen2θ−m2[Θ(θ)=0.sen2θd2Θdθ2+senθcosθdΘdθ+[l(l+1)sen2θ-m2[Θ(θ)=0. (16)
Making the following change in the Eq. (16) cos θ → x and Θ → y, to get:
dΘdθd2Θdθ2sen2θ=−senθdΘdx;=sen2θd2Θdx2−cosθdθdx;=1−x2.