|
|
Sample directory: adf/e_ESR_TiF3/
You calculate Electron Spin Resonance properties with the keywords ESR and QTENS. ESR is a block-type key and is used to compute the G-tensor or the Nuclear Magnetic Dipole Hyperfine interaction. QTENS is a simple key and invokes the computation of the Nuclear Electric Quadrupole Hyperfine interaction.
Proper usage of the key ESR requires that you do either
(a) A Spin-Orbit calculation, spin-restricted, with exactly one unpaired
electron, or
(b) No Spin-Orbit terms and spin-unrestricted.
In case (a) you obtain the G-tensor. In case (b) you get the Magnetic Dipole Hyperfine interaction.
Note: in case (a) the program also prints a Magnetic Dipole Hyperfine interaction data, but these have then been computed without the terms from the spin-density at the nucleus.
After the preliminary calculations (DIRAC, to get the relativistic TAPE12 file with relativistic potentials, and the Create runs), we first calculate the Dipole Hyperfine interaction: a spin-unrestricted calculation without Spin-Orbit coupling.
title TiF3 relativistic unrestricted
noprinf sfo, frag, functions
DEFINE
RTIF = 1.780
RY = RTIF*SQRT(3)/2
END
esr
end
qtens
atoms
Ti 0 0 0
F RTIF 0 0
F -RTIF/2 RY 0
F -RTIF/2 -RY 0
end
fragments
Ti t21.ti
F t21.f
end
xc
GGA Becke Perdew
end
charge 0 1
unrestricted
relativistic scalar zora
Corepotentials t12.rel &
Ti 1
F 2
end
end input
eor
Then, for the same molecule, we compute the G-tensor in a Spin-Orbit run (spin-restricted).
The here-computed and printed Dipole Hyperfine interaction misses the terms from the spin-density at the nucleus: compare with the outcomes from the first calculation.
In each of the calculations, the QTENS key invokes the computation of the Electric Quadrupole Hyperfine interaction.
Note that an all-electron calculation is carried out. This is relevant for the computation of the A-tensor, the nuclear magnetic dipole hyperfine interaction, where an accurate value of the spin-polarization density at the nucleus is important. For the G-tensor (and also for the Q-tensor) this plays a minor role, but for reasons of consistency both calculations use the same basis set and (absence of) frozen core.
title TiF3 relativistic spinorbit
noprinf sfo, frag, functions
DEFINE
RTIF = 1.780
RY = RTIF*SQRT(3)/2
END
esr
end
qtens
atoms
Ti 0 0 0
F RTIF 0 0
F -RTIF/2 RY 0
F -RTIF/2 -RY 0
end
fragments
Ti t21.ti
F t21.f
end
xc
GGA Becke Perdew
end
relativistic spinorbit zora
Corepotentials t12.rel &
Ti 1
F 2
end
end input
eor