# Perturbed angular correlation

The perturbed γ-γ angular correlation, PAC for short or PAC-Spectroscopy, is a method of nuclear solid-state physics with which magnetic and electric fields in crystal structures can be measured. In doing so, electrical field gradients and the Larmor frequency in magnetic fields as well as dynamic effects are determined. With this very sensitive method, which requires only about 10–1000 billion atoms of a radioactive isotope per measurement, material properties in the local structure, phase transitions, magnetism and diffusion can be investigated. The PAC method is related to nuclear magnetic resonance and the Mössbauer effect, but shows no signal attenuation at very high temperatures. Today only the time-differential perturbed angular correlation (TDPAC) is used.

## History and development

PAC goes back to a theoretical work by Donald R. Hamilton from 1940. The first successful experiment was carried out by Brady and Deutsch in 1947. Essentially spin and parity of nuclear spins were investigated in these first PAC experiments. However, it was recognized early on that electric and magnetic fields interact with the nuclear moment, providing the basis for a new form of material investigation: nuclear solid-state spectroscopy.
Step by step the theory was developed. After Abragam and Pound published their work on the theory of PAC in 1953 including extra nuclear fields, many studies with PAC were carried out afterwards. In the 1960s and 1970s, interest in PAC experiments sharply increased, focusing mainly on magnetic and electric fields in crystals into which the probe nuclei were introduced. In the mid-1960s, ion implantation was discovered, providing new opportunities for sample preparation. The rapid electronic development of the 1970s brought significant improvements in signal processing. From the 1980s to the present, PAC has emerged as an important method for the study and characterization of materials, e.g. for the study of semiconductor materials, intermetallic compounds, surfaces and interfaces, and a number of applications have also appeared in biochemistry.
While until about 2008 PAC instruments used conventional high-frequency electronics of the 1970s, in 2008 Christian Herden and Jens Röder et al. developed the first fully digitized PAC instrument that enables extensive data analysis and parallel use of multiple probes. Replicas and further developments followed.

## Measuring principle

PAC uses radioactive probes, which have an intermediate state with decay times of 2 ns to approx. 10 μs, see example $^{111}$In in the picture on the right. After electron capture (EC), indium transmutates to cadmium. Immediately thereafter, the $^{111}$cadmium nucleus is predominantly in the excited 7/2+ nuclear spin and only to a very small extent in the 11/2- nuclear spin, the latter should not be considered further. The 7/2+ excited state transitions to the 5/2+ intermediate state by emitting a 171 keV γ-quantum. The intermediate state has a lifetime of 84.5 ns and is the sensitive state for the PAC. This state in turn decays into the 1/2+ ground state by emitting a γ-quantum with 245 keV. PAC now detects both γ-quanta and evaluates the first as a start signal, the second as a stop signal.
Now one measures the time between start and stop for each event. This is called coincidence when a start and stop pair has been found. Since the intermediate state decays according to the laws of radioactive decay, one obtains an exponential curve with the lifetime of this intermediate state after plotting the frequency over time. Due to the non-spherically symmetric radiation of the second γ-quantum, the so-called anisotropy, which is an intrinsic property of the nucleus in this transition, it comes with the surrounding electrical and/or magnetic fields to a periodic disorder (hyperfine interaction). The illustration of the individual spectra on the right shows the effect of this disturbance as a wave pattern on the exponential decay of two detectors, one pair at 90° and one at 180° to each other. The waveforms to both detector pairs are shifted from each other. Very simply, one can imagine a fixed observer looking at a lighthouse whose light intensity periodically becomes lighter and darker. Correspondingly, a detector arrangement, usually four detectors in a planar 90 ° arrangement or six detectors in an octahedral arrangement, "sees" the rotation of the core on the order of magnitude of MHz to GHz.
According to the number n of detectors, the number of individual spectra (z) results after z=n²-n, for n=4 therefore 12 and for n=6 thus 30. In order to obtain a PAC spectrum, the 90° and 180° single spectra are calculated in such a way that the exponential functions cancel each other out and, in addition, the different detector properties shorten themselves. The pure perturbation function remains, as shown in the example of a complex PAC spectrum. Its Fourier transform gives the transition frequencies as peaks.
$R(t)$, the count rate ratio, is obtained from the single spectra by using:
    $R(t)=2{\frac {W(180^{\circ },t)-W(90^{\circ },t)}{W(180^{\circ },t)+2W(90^{\circ },t)}}$
Depending on the spin of the intermediate state, a different number of transition frequencies show up. For 5/2 spin, 3 transition frequencies can be observed with the ratio ω$_1$+ω$_2$=ω$_3$. As a rule, a different combination of 3 frequencies can be observed for each associated site in the unit cell.
PAC is a statistical method: Each radioactive probe atom sits in its own environment. In crystals, due to the high regularity of the arrangement of the atoms or ions, the environments are identical or very similar, so that probes on identical lattice sites experience the same hyperfine field or magnetic field, which then becomes measurable in a PAC spectrum. On the other hand, for probes in very different environments, such as in amorphous materials, a broad frequency distribution or no is usually observed and the PAC spectrum appears flat, without frequency response. With single crystals, depending on the orientation of the crystal to the detectors, certain transition frequencies can be reduced or extinct, as can be seen in the example of the PAC spectrum of zinc oxide (ZnO).

## Instrumental setup

In the typical PAC spectrometer, a setup of four 90° and 180° planar arrayed detectors or six octahedral arrayed detectors are placed around the radioactive source sample. The detectors used are scintillation crystals of BaF$_2$ or NaI. For modern instruments today mainly LaBr$_3$:Ce or CeBr$_3$ are used. Photomultipliers convert the weak flashes of light into electrical signals generated in the scintillator by gamma radiation. In classical instruments these signals are amplified and processed in logical AND/OR circuits in combination with time windows the different detector combinations (for 4 detectors: 12, 13, 14, 21, 23, 24, 31, 32, 34, 41, 42, 43) assigned and counted. Modern digital spectrometers use digitizer cards that directly use the signal and convert it into energy and time values and store them on hard drives. These are then searched by software for coincidences. Whereas in classical instruments, "windows" limiting the respective γ-energies must be set before processing, this is not necessary for the digital PAC during the recording of the measurement. The analysis only takes place in the second step. In the case of probes with complex cascades, this makes it makes it possible to perform a data optimization or to evaluate several cascades in parallel, as well as measuríng different probes simultaneously. The resulting data volumes can be between 60 and 300 GB per measurement.

## Sample materials

As materials for the investigation (samples) are in principle all materials that can be solid and liquid. Depending on the question and the purpose of the investigation, certain framework conditions arise. For the observation of clear perturbation frequencies it is necessary, due to the statistical method, that a certain proportion of the probe atoms are in a similar environment and e.g. experiences the same electric field gradient. Furthermore, during the time window between the start and stop, or approximately 5 half-lives of the intermediate state, the direction of the electric field gradient must not change. In liquids, therefore, no interference frequency can be measured as a result of the frequent collisions, unless the probe is complexed in large molecules, such as in proteins. The samples with proteins or peptides are usually frozen to improve the measurement.
The most studied materials with PAC are solids such as semiconductors, metals, insulators, and various types of functional materials. For the investigations, these are usually crystalline. Amorphous materials do not have highly ordered structures. However, they have close proximity, which can be seen in PAC spectroscopy as a broad distribution of frequencies. Nano-materials have a crystalline core and a shell that has a rather amorphous structure. This is called core-shell model. The smaller the nanoparticle becomes, the larger the volume fraction of this amorphous portion becomes. In PAC measurements, this is shown by the decrease of the crystalline frequency component in a reduction of the amplitude (attenuation).

## Sample preparation

The amount of suitable PAC isotopes required for a measurement is between about 10 to 1000 billion atoms (10$^{10}$-10$^{12}$). The right amount depends on the particular properties of the isotope. 10 billion atoms are a very small amount of substance. For comparison, one mol contains about 6.22x10$^{23}$ particles. 10$^{12}$ atoms in one cubic centimeter of beryllium give a concentration of about 8 nmol/L (nanomol=10$^{−9}$ mol). The radioactive samples each have an activity of 0.1-5 MBq, which is in the order of the exemption limit for the respective isotope.
How the PAC isotopes are brought into the sample to be examined is up to the experimenter and the technical possibilities. The following methods are usual:

### Implantation

During implantation, a radioactive ion beam is generated, which is directed onto the sample material. Due to the kinetic energy of the ions (1-500 keV) these fly into the crystal lattice and are slowed down by impacts. They either come to a stop at interstitial sites or push a lattice-atom out of its place and replace it. This leads to a disruption of the crystal structure. These disorders can be investigated with PAC. By tempering these disturbances can be healed. If, on the other hand, radiation defects in the crystal and their healing are to be examined, unperseived samples are measured, which are then annealed step by step.
The implantation is usually the method of choice, because it can be used to produce very well-defined samples.

### Evaporation

In a vacuum, the PAC probe can be evaporated onto the sample. The radioactive probe is applied to a hot plate or filament, where it is brought to the evaporation temperature and condensed on the opposite sample material. With this method, e.g. surfaces are examined. Furthermore, by vapor deposition of other materials, interfaces can be produced. They can be studied during tempering with PAC and their changes can be observed. Similarly, the PAC probe can be transferred to sputtering using a plasma.

### Diffusion

In the diffusion method, the radioactive probe is usually diluted in a solvent applied to the sample, dried and it is diffused into the material by tempering it. The solution with the radioactive probe should be as pure as possible, since all other substances can diffuse into the sample and affect thereby the measurement results. The sample should be sufficiently diluted in the sample. Therefore, the diffusion process should be planned so that a uniform distribution or sufficient penetration depth is achieved.

### Added during synthesis

PAC probes may also be added during the synthesis of sample materials to achieve the most uniform distribution in the sample. This method is particularly well suited if, for example, the PAC probe diffuses only poorly in the material and a higher concentration in grain boundaries is to be expected. Since only very small samples are necessary with PAC (about 5 mm), micro-reactors can be used. Ideally, the probe is added to the liquid phase of the sol-gel process or one of the later precursor phases.

### Neutron activation

In neutron activation, the probe is prepared directly from the sample material by converting very small part of one of the elements of the sample material into the desired PAC probe or its parent isotope by neutron capture. As with implantation, radiation damage must be healed. This method is limited to sample materials containing elements from which neutron capture PAC probes can be made. Furthermore, samples can be intentionally contaminated with those elements that are to be activated. For example, hafnium is excellently suited for activation because of its large capture cross section for neutrons.

### Nuclear reaction

Rarely used are direct nuclear reactions in which nuclei are converted into PAC probes by bombardment by high-energy elementary particles or protons. This causes major radiation damage, which must be healed. This method is used with PAD, which belongs to the PAC methods.

## Laboratories

The currently largest PAC laboratory in the world is located at ISOLDE in CERN with about 10 PAC instruments, that receives its major funding form BMBF. Radioactive ion beams are produced at the ISOLDE by bombarding protons from the booster onto target materials (uranium carbide, liquid tin, etc.) and evaporating the spallation products at high temperatures (up to 2000 °C), then ionizing them and then accelerating them. With the subsequent mass separation usually very pure isotope beams can be produced, which can be implanted in PAC samples. Of particular interest to the PAC are short-lived isomeric probes such as: $^{111m}$Cd, $^{199m}$Hg, $^{204m}$Pb, and various rare earth probes.

## Theory

The first $\gamma$-quantum ($\gamma _{1},k_{1}$) will be emitted isotropically. Detecting this quantum in a detector selects a subset with an orientation of the many possible directions that has a given. The second $\gamma$-quantum ($\gamma _{2},k_{2}$) has an anisotropic emission and shows the effect of the angle correlation. The goal is to measure the relative probability $W(\Theta ){\textrm {d}}(\Omega )$ with the detection of $\gamma _{2}$ at the fixed angle $\Theta$ in relation to $\gamma _{1}$. The probability is given with the angle correlation (perturbation theory):
    $W(\Theta )=\sum _{k}^{k_{max}}A_{kk}P_{k}cos(\Theta )$
For a $\gamma$-$\gamma$-cascade, $k$ is due to the preservation of parity:
    $0<k<{\textrm {min}}(2I_{S},I_{i}+I'_{i})$
Where $I_{S}$ is the spin of the intermediate state and $I_{i}$ with $i=1;2$ the multipolarity of the two transitions. For pure multipole transitions, is $I_{i}=I'_{i}$.
$A_{kk}$ is the anisotropy coefficient that depends on the angular momentum of the intermediate state and the multipolarities of the transition.
The radioactive nucleus is built into the sample material and emits two $\gamma$-quanta upon decay. During the lifetime of the intermediate state, i.e. the time between $\gamma _{1}$ and $\gamma _{2}$, the core experiences a disturbance due to the hyperfine interaction through its electrical and magnetic environment. This disturbance changes the angular correlation to:
    $W(\Theta )=\sum _{k}^{k_{max}}A_{kk}G_{kk}$
$G_{kk}$ is the perturbation factor. Due to the electrical and magnetic interaction, the angular momentum of the intermediate state $I_{i}$ experiences a torque about its axis of symmetry. Quantum-mechanically, this means that the interaction leads to transitions between the M states. The second $\gamma$-quantum ($\gamma _{2}$) is then sent from the intermediate level. This population change is the reason for the attenuation of the correlation.
The interaction occurs between the magnetic core dipole moment ${\vec {\nu }}$ and the intermediate state $I_{S}$ or/and an external magnetic field ${\vec {B}}$. The interaction also takes place between nuclear quadrupole moment and the off-core electric field gradient $V_{zz}$.

### Magnetic dipole interaction

For the magnetic dipole interaction, the frequency of the precession of the nuclear spin around the axis of the magnetic field ${\vec {B}}$ is given by:
    $\omega _{L}={\frac {g\cdot u_{N}\cdot B}{\hbar }}$
    $\Delta E=\hbar \cdot \omega _{L}=-g\cdot u_{N}\cdot B$
$g$ is the Landé g-factor und $u_{N}$ is the nuclear magneton.
With $N=M-M'$ follows:
    $E_{magn}(M)-E_{magn}(M')=-(M-M')g\mu _{N}B_{z}=N\hbar \omega _{L}$
From the general theory we get:
    $G_{k_{1}k_{2}}^{NN}={\sqrt {(2k_{1}+1)(2k_{2}+1)}}\cdot e^{-iN\omega _{L}t}\times \sum _{M}{\begin{pmatrix}I&I&k_{1}\\M'&-M&N\\\end{pmatrix}}{\begin{pmatrix}I&I&k_{2}\\M'&-M&N\\\end{pmatrix}}$
For the magnetic interaction follows:
    $G_{k_{1}k_{2}}^{NN}=e^{\left({-iN\omega _{L}t}\right)}$

### Static electric quadrupole interaction

The energy of the hyperfine electrical interaction between the charge distribution of the core and the extranuclear static electric field can be extended to multipoles. The monopole term only causes an energy shift and the dipole term disappears, so that the first relevant expansion term is the quadrupole term:
    $E_{Q}=\sum _{ij}Q_{ij}V_{ij}$ ij=1;2;3
This can be written as a product of the quadrupole moment $Q_{ij}$ and the electric field gradient $V_{ij}$. Both [tensor]s are of second order. Higher orders have too small effect to be measured with PAC.
The electric field gradient is the second derivative of the electric potential $\Phi ({\vec {r}})$ at the core:
    $V_{ij}={\frac {\partial ^{2}\Phi ({\vec {r}})}{\partial x_{i}\partial x_{j}}}={\begin{pmatrix}V_{xx}&0&0\\0&V_{yy}&0\\0&0&V_{zz}\\\end{pmatrix}}$
$V_{ij}$ becomes diagonalized, that:
    $|V_{zz}|\geq |V_{yy}|\geq |V_{xx}|$
The matrix is free of traces in the main axis system (Laplace equation)
    $V_{xx}+V_{yy}+V_{zz}=0$
Typically, the electric field gradient is defined with the largest proportion $V_{zz}$ and $\eta$:
    $\eta ={\frac {V_{yy}-V_{xx}}{V_{zz}}}$, $0\leq \eta \leq 1$
In cubic crystals, the axis parameters of the unit cell x, y, z are of the same length. Therefore:
    $V_{zz}=V_{yy}=V_{xx}$ and $\eta =0$
In axisymmetric systems is $\eta =0$.
For axially symmetric electric field gradients, the energy of the substates has the values:
    $E_{Q}={\frac {eQV_{zz}}{4I(2I-1)}}\cdot (3m^{2}-I(I+1))$
The energy difference between two substates, $M$ and $M'$, is given by:
    $\Delta E_{Q}=E_{m}-E_{m'}={\frac {eQV_{zz}}{4I(2I-1)}}\cdot 3|M^{2}-M'^{2}|$
The quadrupole frequency $\omega _{Q}$ is introduced. The formulas in the colored frames are important for the evaluation:
       $\omega _{Q}={\frac {eQV_{zz}}{4I(2I-1)\hbar }}={\frac {2\pi eQV_{zz}}{4I(2I-1)h}}={\frac {2\pi \nu _{Q}}{4I(2I-1)}}$
       $\nu _{Q}={\frac {eQ}{h}}V_{zz}={\frac {4I(2I-1)\omega _{Q}}{2\pi }}$
The publications mostly list $\nu _{Q}$. $e$ as elementary charge and $h$ as Planck constant are well known or well defined. The nuclear quadrupole moment $Q$ is often determined only very inaccurately (often only with 2-3 digits). Because $\nu _{Q}$ can be determined much more accurately than $Q$, it is not useful to specify only $V_{zz}$ because of the error propagation. In addition, $\nu _{Q}$ is independent of spin! This means that measurements of two different isotopes of the same element can be compared, such as $^{199m}$Hg(5/2−), $^{197m}$Hg(5/2−) and $^{201m}$Hg(9/2−). Further, $\nu _{Q}$ can be used as finger print method.
For the energy difference then follows:
    $\Delta E_{Q}=\hbar \omega _{Q}\cdot 3|m^{2}-m'^{2}|$
If $\eta =0$, then:
    $\omega ^{n}=n\cdot \omega _{Q}^{0}$
with:
    $\omega _{Q}^{0}={\textrm {min}}\left({\frac {\Delta E_{Q}}{\hbar }}\right)$
For integer spins applies:
    $\omega _{Q}^{0}=3\cdot \omega _{Q}$ und $n=|M^{2}-M'^{2}|$
For half integer spins applies:
    $\omega _{Q}^{0}=6\cdot \omega _{Q}$ und $n={\frac {1}{2}}|M^{2}-M'^{2}|$
The perturbation factor is given by:
    $G_{k_{1}k_{2}}^{NN}=\sum _{n}s_{nN}^{k_{1}k_{2}}\cos {(n\omega _{Q}^{0}t)}$
With the factor for the probabilities of the observed frequencies:
    $s_{nN}^{k_{1}k_{2}}={\sqrt {(2k_{1}+1)(2k_{2}+1)}}\cdot \sum _{M,M'}{\begin{pmatrix}I&I&k_{1}\\M'&-M&N\\\end{pmatrix}}{\begin{pmatrix}I&I&k_{2}\\M'&-M&N\\\end{pmatrix}}$
As far as the magnetic dipole interaction is concerned, the electrical quadrupole interaction also induces a precision of the angular correlation in time and this modulates the quadrupole interaction frequency. This frequency is an overlap of the different transition frequencies $\omega _{n}$. The relative amplitudes of the various components depend on the orientation of the electric field gradient relative to the detectors (symmetry axis) and the asymmetry parameter $\eta$. For a probe with different probe nuclei, one needs a parameter that allows a direct comparison: Therefore, the quadrupole coupling constant $\nu _{Q}$ independent of the nuclear spin ${\vec {I}}$ is introduced.

### Combined interactions

If there is a magnetic and electrical interaction at the same time on the radioactive nucleus as described above, combined interactions result. This leads to the splitting of the respectively observed frequencies. The analysis may not be trivial due to the higher number of frequencies that must be allocated. These then depend in each case on the direction of the electric and magnetic field to each other in the crystal. PAC is one of the few ways in which these directions can be determined.

### Dynamic interactions

If the hyperfine field fluctuates during the lifetime $\tau _{n}$ of the intermediate level due to jumps of the probe into another lattice position or from jumps of a near atom into another lattice position, the correlation is lost. For the simple case with an undistorted lattice of cubic symmetry, for a jump rate of $\omega _{s}<0.2\cdot \nu _{Q}$ for equivalent places $N_{s}$, an exponential damping of the static $G_{22}(t)$-terms is observed:
    $G_{22}^{dyn}(t)=e^{-\lambda _{d}t}G_{22}(t)$ $\lambda _{d}=(N_{s}-1)\omega _{s}$
Here $\lambda _{d}$ is a constant to be determined, which should not be confused with the decay constant $\lambda ={\frac {1}{\tau }}$. For large values of $\omega _{s}$, only pure exponential decay can be observed:
    $G_{22}^{dyn}(t)=e^{-\lambda _{d}t}$
The boundary case after Abragam-Pound is $\lambda _{d}$, if $\omega _{s}>3\cdot \nu _{Q}$, then:
    $\lambda _{d}\approx {\frac {2,5\nu _{Q}^{2}}{N_{s}\omega _{s}}}$

### After effects

Cores that transmute beforehand of the $\gamma$-$\gamma$-cascade usually cause a charge change in ionic crystals (In$^{3+}$) to Cd$^{2+}$). As a result, the lattice must respond to these changes. Defects or neighboring ions can also migrate. Likewise, the high-energy transition process may cause the Auger effect, that can bring the core into higher ionization states. The normalization of the state of charge then depends on the conductivity of the material. In metals, the process takes place very quickly. This takes considerably longer in semiconductors and insulators. In all these processes, the hyperfine field changes. If this change falls within the $\gamma$-$\gamma$-cascade, it may be observed as an after effect.
The number of nuclei in state (a) in the image on the right is depopulated both by the decay after state (b) and after state (c):
    $\mathrm {d} N_{a}=-N_{a}\left(\Gamma _{r}+{\frac {1}{\tau _{7/2}}}\right)\mathrm {d} t$
mit: $\tau _{7/2}={\frac {120{\textrm {ps}}}{\ln {2}}}$
From this one obtains the exponential case:
    $N_{a}(t)=N_{a_{0}}\cdot e^{\left({-(\Gamma _{r}+{\frac {1}{\tau _{7/2}}})t}\right)}$
For the total number of nuclei in the static state (c) follows:
    $N_{c}(t)=\Gamma _{r}\int \limits _{0}^{t}N_{a}(t)\mathrm {d} t=N_{0}{\frac {\Gamma _{r}\tau _{7/2}}{\Gamma _{r}\tau _{7/2}+1}}\left(1-e^{-(\Gamma _{r}+{\frac {1}{\tau _{7/2}}})t}\right)$
The initial occupation probabilities $\rho$ are for static and dynamic environments:
    $\rho _{stat}={\frac {\Gamma _{r}\tau _{7/2}}{\Gamma _{r}\tau _{7/2}+1}}$
    $\rho _{dyn}={\frac {1}{\Gamma _{r}\tau _{7/2}+1}}$

### General theory

In the general theory for a transition $M_{i}\rightarrow M_{f}$ is given:
    $W(M_{i}\rightarrow M_{f})=\left|\sum _{M}\langle M_{f}|{\mathcal {H}}_{2}|M\rangle \langle M|{\mathcal {H}}_{1}|M_{i}\rangle \right|^{2}$
    $W({\vec {k}}_{1},{\vec {k}}_{2})=\sum _{M_{i},M_{f},\sigma _{1},\sigma _{2}}\left|\sum _{M}\langle M_{f}|{\mathcal {H}}_{2}|M\rangle \langle M|{\mathcal {H}}_{1}|M_{i}\rangle \right|^{2}$
    $W({\vec {k}}_{1},{\vec {k}}_{2})=W(\Theta )=\sum _{k_{gerade}}^{k_{max}}A_{k}(1)A_{k}(2)P_{k}(\cos {\Theta })$
    $0\leq k\leq$ Minimum von $(2I,l_{1}+l_{1}',l_{2}+l_{2}')$
    $W(\Theta ,t)=\sum _{k=2,4}A_{kk}P_{k}(\cos {\Theta })$
    $|M_{a}\rangle \rightarrow \Lambda (t)|M_{a}=\sum _{M_{b}}|M_{b}\rangle \langle M_{b}|\Lambda (t)|M_{a}\rangle$
    $W({\vec {k}}_{1},{\vec {k}}_{2},t)=\sum _{M_{i},M_{f},\sigma _{1},\sigma _{2}}\left|\sum _{M_{a}}\langle M_{f}|{\mathcal {H}}_{2}\Lambda (t)|M_{a}\rangle \langle M_{a}|{\mathcal {H}}_{1}|M_{i}\rangle \right|^{2}=\langle \rho ({\vec {k}}_{2})\rangle _{t}$
    $W({\vec {k}}_{1},{\vec {k}}_{2},t)=\sum _{k_{1},k_{2},N_{1},N_{2}}A_{k_{1}}(1)A_{k_{2}}(2){\frac {1}{\sqrt {(2k_{1}+1)(2k_{2}+1)}}}\times Y_{k_{1}}^{N_{1}}(\Theta _{1},\Phi _{1})\cdot Y_{k_{2}}^{N_{2}}(\Theta _{2},\Phi _{2})G_{k_{1}k_{2}}^{N_{1}N_{2}}(t)$
with:
    $G_{k_{1}k_{2}}^{N_{1}N_{2}}=\sum _{M_{a},M_{b}}(-1)^{2I+M_{a}+M_{b}}{\sqrt {(2k_{1}+1)(2k_{2}+)}}\times \langle M_{b}|\Lambda (t)|M_{a}\rangle \langle M_{b}'|\Lambda (t)|M_{a}'\rangle ^{*}\times {\begin{pmatrix}I&I&k_{1}\\M_{a}'&-M_{a}&N_{1}\end{pmatrix}}{\begin{pmatrix}I&I&k_{2}\\M_{b}'&-M_{b}&N_{2}\end{pmatrix}}$