Absolute polarization measurement by elastic scattering

Leave a comment

The magnitude of proton polarization can be measured by NMR technique with a reference. Because the NMR gives the free-induction decay signal, which is a voltage or current. For Boltzmann polarization using strong magnetic field and low temperature, the polarization can be calculated. However, when a reference point is not available, the absolute magnitude of proton polarization can be measured using proton-proton elastic scattering. The principle is the nuclear spin-orbital coupling. That creates left-right asymmetry on the scattering cross section.

Because of spin-orbital interaction:

V_{ls}(r) = f(r) \vec{l} \cdot \vec{s} ,

where f(r) is the distance function, \vec{l} is the relative angular momentum, \vec{s} is the spin of the incident proton. In the following picture, the spin of the incident proton can be either out of the plane (\uparrow ) or into the plan (\downarrow). When the proton coming above, the angular momentum is into the plane (\downarrow ). The 4 possible sign of the spin-orbital interaction is shown. We can see, when the spin is up, the spin-orbital force repulses the proton above and attracts the proton below. That creates an asymmetry in the scattering cross section.



The cross section is distorted and characterized using analysing power A_y. Analyzing power is proportional to the difference between left-right cross-section. By symmetry (parity, time-reversal) consideration, A_y = 1 + P sin(2\theta) (why?), in center of mass frame. In past post, the transformation between difference Lorentz frame. The angle in the A_y has to be expressed in lab angle. The cross section and A_y can be obtained from http://gwdac.phys.gwu.edu/ .

In scattering experiment, the number of proton (yield) is counted in left and right detectors. The yield should be difference when either proton is polarized. The yield is

Y(\theta, \phi) = L \epsilon \sigma_0 (1 + cos(\phi)A_y(\theta) P) ,

where L is the luminosity, \epsilon is the detector efficiency, \sigma_0 is the integrated cross-section of un-polarized beam and target of the detector, P is the polarization of either the target or beam. When both target and the beam are polarized, the cross section is

\sigma = \sigma_0 (1 + (P + P_T)A_y + P P_T C_yy),

where C_yy is spin-spin correlation due to spin-spin interaction of nuclear force.

Using the left-right yield difference, the absolute polarization of the target or the beam can be found using,

\displaystyle A_y P = \frac{Y_L - Y_R}{Y_L + Y_R} ,

where Y_L = Y(\phi =0) and Y_R = Y(\phi=\pi) .




Magnetic Dipole Moment & Gyromagnetic Ratio

Leave a comment

I always confuses on the definition, and wiki did not have any summary. so,

The Original definition is the Hamiltonian of a magnetic dipole under external magnetic field \vec{B},

H = -\vec{\mu}\cdot \vec{B},

where \vec{\mu} is magnetic dipole moment (MDM). It is

\vec{\mu} = g \frac{q}{2 m} \vec{J} = g \frac{\mu}{\hbar} \vec{J} = \gamma \vec{J}.

Here, the g is the g-factor, \mu is magneton, and \vec{J} is the total spin, which has a intrinsic factor m\hbar / 2 inside. \gamma is gyromegnetic ratio.

We can see, the g-factor depends on the motion or geometry of the MDM. For a point particle, the g-factor is exactly equal to 2. For a charged particle orbiting, the g-factor is 1.

Put everything into the Hamiltonian,

H = -\gamma \vec{J}\cdot \vec{B} = -\gamma J_z B = -\gamma \hbar \frac{m}{2} B [J],

Because energy is also equal E = \hbar f , thus, we can see the \gamma has unit of frequency over Tesla.


Take electron as an example, the MDM is Bohr magneton \mu_{e} = e\hbar/(2m_e). The MDM is,

\vec{\mu_e} = g_e \frac{e}{2 m_e} \vec{S} = g_e \frac{\mu_e}{\hbar}\vec{S} = \gamma_e \vec{S}.

The magnitude of MDM is,

|\vec{\mu_e}|= g_e \frac{e}{2 m_e} \frac{\hbar}{2} = \gamma_e \frac{\hbar}{2} [JT^{-1}],

The gyromagnetic ratio is,

\gamma_e = g_e \frac{\mu_e}{\hbar} [rad s^{-1} T^{-1}].

Since using rad s^{-1} is not convenient for experiment. The gyromagnetic ratio usually divided by 2\pi,

\gamma_e = g_e \frac{\mu_e}{2\pi\hbar} [Hz T^{-1}].


To evaluate the magnitude of  MDM of  single particle state, which has orbital angular momentum and spin, the total spin \vec{J} = \vec{L} + \vec{S}. However, the g-factor for \vec{L} is difference from that for \vec{S}. Thus, the MDM is not parallel to total spin. We have to use Landé Formula,

\left< JM|\vec{V}|JM'\right> = \frac{1}{J(J+1)} \left< JM|(\vec{J}\cdot\vec{V})|JM\right> \left<JM|\vec{J}|JM'\right>

or see wiki, sorry for my laziness.

The result is


For J = L \pm 1/2,

g = J(g_L \pm \frac{g_S-g_L}{2L+1})


Leave a comment

Q-factor is Not Q-value. in Q-factor, the Q is for Quality.

Q-factor is a dimensionless factor for showing the degree of reasonace of a coil or an Oscillator.

Higer value means a better coil. The definition is the resonance frequency over the Full Width Half Maximum (FWHM).

Q = f/\Delta f

High value means:

  1. high sensitivity
  2. noise reduction due to narrow band of absorption.
  3. over damped with long decay time.
  4. higher energy stored
  5. lower energy loss

For complex electric circuit, the Q factor is:

Q = | im Z/ re Z|

It seem that there is a conflict between impedance matching.

Finding a 90 degree pulse of the NMR system

Leave a comment

Tuning of NMR System check the pdf… so tired to rewrite again….




T1 and T2 measurement

Leave a comment

Measuring  T1 in NMR, we apply follow  pulse sequence:

\pi_x \longrightarrow \tau \longrightarrow (\pi/2)_x

according to previous post on density matrix in operator form, we can evaluate the polarization. suppose the magnetization is pointing with the external B-field.

\rho_R = \rho = k\sigma_z

\downarrow \pi_x

\rho_R = - k \sigma_z

\downarrow \tau

\rho_R = -k \left( 1- 2 e^{-\tau/T_1} \right) \sigma_z

\downarrow (\pi/2)_x

\rho_R = k \left( 1- 2 e^{-\tau/T_1} \right) \sigma_y

in lab frame:

\rho = k \left( 1-2 e^{-\tau/T_1} \right) (\sigma_y cos(\omega_0 t) - \sigma_x sin(\omega_0 t) )

if the NMR coil is placed alone with x-axis. the magnetization is proportional to :

\left<\sigma_x\right> = 2 k\left( 1- 2 e^{-\tau/T_1} \right) sin(\omega_0 t)

the amplitude of the magnetization is only a function of \tau . bu measuring the amplitude with different \tau , we can determine the T_1 .

To measure the T2, we use follow pulse sequence:

(\pi/2)_x \longrightarrow \tau \longrightarrow (\pi)_y \longrightarrow \tau

again we use the same initial state, but this time, we are going watch it as lab frame.

\rho = k\sigma_z

\downarrow (\pi/2)_x

\rho = - k \sigma_y

\downarrow \tau

\rho = -k (\sigma_y cos( \omega_0 \tau ) - \sigma_x sin(\omega_0 \tau) ) e^{- \tau/T_2}

\downarrow \pi_y

\rho = -k (\sigma_y cos( \omega_0 \tau ) + \sigma_x sin(\omega_0 \tau) ) e^{- \tau/T_2}

\downarrow \tau

\rho = - k \sigma_ye^{- 2 \tau/T_2}

at  the last step, the free-induction decay is not decay but a revert process of decay. we see that the state back to its beginning state! thus, this method also called the spin echo. this can be see in pictorial  way. after the (\pi/2)_x pulse, the spin go to – y axis. now, due to the incoherence of each spin, some spin are faster and some are slower. after a time \tau , the \pi_y pulse flips all the spin by 180 degree. now, the “slower” spin become ahead of the “faster” spin. After the same time interval, the “faster” spin will catch up the “slower” spin and all the spin becomes coherence at that moment again! thus, the amplitude of the signal will become large again. However, since the spins are not at same Larmor frequency, some will flip more than 180 degrees while some flip less, thus, the decay are still there.


the FID ( free induction decay) rate is not exactly equal to T2, since there are many other way to make the spin polarization decay. the FID decay rate T_2^* should be :

\frac{1}{T_2 ^*} = \frac{1}{T_2} + \frac{1}{T_{others}}

and the spin echo method can eliminate the others.

[Pol. p target] Matching Impedance of NMR coil


i had played with the impedance of coax cable. connect the signal generator at 50Ω output and a 4meter long coax to a CRO with 50Ω input with a tee. although the web said, the coax impedance is 52, there is no observed different. the peak-to peak signal does not depend on the frequency.

then connect a 12.8MHz coax to the tee and open at the other end. in theory, it should be zero input-impedance and he CRO reading should be

V_{in}= 1 - cos( 2\pi L/\lambda)

where the CRO reading is the input voltage. when the frequency adjusted to 12.8MHz, which is quata-wavelength, the input voltage is equal to 1.

when connect the end of the coax with a 50Ω resistor, the input voltage does not depends on frequency, as expected.

i haven’t try to short the end of the coax.

after that. i going to matching the impedance of the tuner and coil. we use a short coax to connect the input to the tuner before, i replaced it with a 12.8MHz coax. the reason is, the little bit mismatching of the impedance can be saved by the length of the cable, such that:

Z_{in} = Z_0^2/Z_L \approx Z_0

and i opened the tuner to see the circuit inside:

The above is the circuit diagram. there is a fixed capacitor with 2.2pF in parallel. i cannot identify the type of the 2 variable capacitor.  the coax cable can be put in port 1 or port 2 and the coil put in the other port. different configuration has different behavior.

i found that, the input is in port 2 and the coil is in port 1, and the box doesn’t ground but just wrapped with metal sticker. i took it out and grounded to NMR system.

For the input at port 2, the total impedance of the coil and tuner is:

Z_L= - i /(\omega (C_0+C_p + C_s/(1-\omega^2 L C_s)))

The first things to notice is the impedance depends on the frequency. which mean, the impedance matching can only on particular frequency. when the Cs adjusted to matching the frequency. the impedance solely depends on Cp. when the driving voltage gone, the LC circuit will oscillate at natural frequency:

\omega_n = 1/\sqrt{ L ( C_s + C_0 C_p / (C_0+C_p) ) }

Thus, i tired to measure to inductance of the coil by a parallel resistance. but i cannot find any suitable wire to convert BNC cable to wires. after a long time finding, i gave up and wait unit work with my partner.

the input in port 2 is not a common config, so, i changed it to port 1. and the impedance is :

Z_L=i\omega (C_s +L/(1-\omega^2L(C_0+ C_p)))

The impedance also depends on frequency. and the natural frequency when the driving voltage gone is:

\omega_n = 1/\sqrt{L(C_0+C_p)}

there is one way to tune the Cp to match the input frequency and make load impedance solely depends on Cs. by using a pulse signal. and measure the natural frequency of the LC loop, such that the natural frequency is same as input frequency. However, the Low Pass Filter only let frequency less then 1MHz pass and out frequency is 12.8MHz. can i use other pulse source? may be, if i have a mixer.

so, i matching the impedance by very naive way. i use the method on testing the impedance matching. i use a continuous signal source and fixed the frequency at 12.8MHz, then connect it with a coax cable to 50Ω CRO input by a tee, then connect th tee with a 12.8MHz coax cable. the other end of the cable connect to a 50Ω resistor. this setup should be matched impedance. so, i record the input voltage on the CRO. and replace the resistor with the tuner, which port 2 connected to the coil. then adjust the capacitors (both) so that the CRO voltage is same as 50Ω resistor.

by solving the load impedance formula of port 1 configuration, there are mulitple solution for Cp and Cs to give 50Ω. and i think, as long as the impedance is 50Ω at 12.8MHz, any configuration can do the job.

later, i try to find the water NMR signal. although i cannot find any. but the noise level reduced to ±5mV. more or less equal to the background.

i played with the NMR program. the record data is counted by point, so the CRO horizontal setting should set to 5000 points over the screen.

and just before i leave, i don’t know what wrong, the program doesn’t read the CRO signal…

Dynamic Nuclear Polarization

Leave a comment

DOI : 10.1103/RevModPhys.34.173

The Dynamic Nuclear Polarization (DNP) means we has a pumping source to change the population of nuclear spin, then create a polarization. in contrast, Static Nuclear Polarization (SNP) means thermal equilibrium of nuclear state population.

the introduction of the paper gives 7 applications on polarized nuclear spin.i only list some below:

  1. the angular distribution on radiations can serve as a test on the theory of nuclear interaction
  2. Polarized target can be used in scatter experiment
  3. obtain detail information on static and dynamic interaction between nuclear spin and its environment.
  4. increase the sensitivity of NMR

this paper focus on a general system and represents them by graphs ( called chart in the paper ). the graphs are based on electron spin ½ and nuclear spin also ½.

on section II, it give out the Spin Hamiltonian and use it for the discussion on the population distribution. by that, the author used the rate equations to related the population in each state. Then, he defined the Enhancement of polarization, which is the ratio between the population with saturating radiation to the thermal thermal distribution.

on section III, it mention about the first 2 successful dynamic nuclear polarization experiments around 1953-4. one group polarized the 6Li nucleus in metallic lithium. the other group polarized the 1H in solid DPPH.

The paper gives conditions for DNP, which is coupling between nuclear spin and an unpaired electron spin. the paramagnetic environment can be archived by

  1. the conduction electron in metals or metal ammonia solution
  2. the donor or acceptor electrons in semi-conductor
  3. paramagnetic ions in diamagnetic solid
  4. paramagnetic ions in solution
  5. free radical
  6. color centers

the detection of DNP can be via:

  1. NMR
  2. shift of EPR frequency
  3. the β asymmetry or γ anisotropy from an oriented radioisotope

[Pol. p target] Looking for a signal

Leave a comment

the NMR signal was calibrated by water, since it has contains proton. the pule frequency is 12.78 MHz. this frequency is also equal to the Larmor frequency under a static B field.

2 \pi f = \mu_p B /2

the time of the pule is 5 μs. this should be the π/2 pulse. thus

\mu_p B_R \tau_{\pi/2} /2 = \pi / 2

we did a parameter search for any possible NMR signal. the parameters are the microwave power and the delay time.

(the following is a rough understanding, many places are unclear. )

the dynamic nuclear polarization requires a magnetic field and an microwave. the microwave is function as a power source or photo source. when nuclear spin flips, it requires a photo to provide the necessary energy. since our microwave frequency is fixed by the microwave cavity, so, in order to matching the correct frequency of absorption, we use a Sweeping magnetic field, when the magnetic field is at right magnitude, than the energy gap between spin up and spin down will equal to the photo energy and then absorbed.

after the electrons are all polarized by optical pumping, which means the electron polarization is saturated. the applied sweeping B-field and microwave will able to transfer the polarization to the proton.

thus, when the laser stopped by a chopper, the sweeping field and microwave will be turned on. the time delay between the laser and the sweeping field is called delay time.

the upper line is the chopper signal ( channel 4) . the next one is the sweeping field ( channel 2 ). the 3rd one is the reflected microwave, we can see a very small change on it, the base line are off screen, since the signal is very small compare to the overall signal ( channel 1). the last one is the RF Amp trigger signal, which control the microwave on-off ( channel 3 ).

[Pol. p target] NMR system and Control Unit

Leave a comment

The NMR system was re-programmed. Not by me.

the trigger signal is from the laser beam. there are 2 ways, 1 is from the chopper diode, another is by the power meter. the chopper signal is typical square wave. But the power meter has phase delay and relaxation time.

the field sweep is a amplifier connected to the magnetic coil. the amplifier is APEX-PA05 model. we use it as a Voltage-Controlled-Current-Supply. since the magnetic field is depends on Current and the coil resistance is so small that applying a voltage on it will very sensitive.

the RF Amp control is usual, we connect a voltage output in to the “Trig IN”. but the RF Amp just like turned off. when the control changed back to local, every things work well. so, the problem is on the Remote setting of the RF Amp. we checked the manual and it said that to complete control from remote, we may have to connect it with a computer. Oh….

Analyzing power for proton elastic scattering from the neutron-rich 6He nucleus

Leave a comment

DOI : 10.1103/PhysRevC.82.021602

This paper is based on the solid proton polarized target, and make improvement on the experiment. The proton polarization was monitored by NMR method and the absolute polarization is known. Therefore, the Analyzing power ( or the spin asymmetry) can be determined. The experimental result is different form the prediction of the t-matrix folding model.


The 1st and 2nd paragraphs explain why the spin-asymmetry is important in nuclear research. a simple reason is, by manipulating the spin, we can reveal the spin-dependence of the nuclear-structure. For example, the orbital-spin coupling, which is the major factor in the nuclear potential.

the 3rd, 4th and 5th paragraphs talk about the solid proton polarization target and the advantage of operating under weak magnetic field. (already discussed in Here )

The 6th paragraph report the result from the 2003 experiment. and say that the unsatisfiable result on the spin-asymmetry.

The 7th to 11th paragraphs talk about the experiment procedures and conditions.

the 12th paragraph explains the angular distribution in center of mass frame of the differential cross section. the data of 6He is similar to the 6Li for angle less then 50 degrees, which is the forward angle in the lab frame. at the backward angle or angle more then 50 degree in C.M. frame, there are some different and it may be due to the halo structure of 6He.

the next paragraph turns the focus on the analyzing power, which can be determined at this time. the data is very different from 6Li data. the t-matrix folding method predicted that the analyzing power should be positive but the experiment result is different. this indicated that there is other mechanism to explain the nuclear structure of 6He.

The t-matrix stands for transition matrix.

the other mechanism may be the g-matrix folding model…..i am sorry, i don’t understand the following….

Older Entries