Property gradients

Property gradients#

The first-order nuclear derivatives of properties other than the energy can be related to intensities in vibrational spectroscopies. More specifically, the derivative of the dipole moment and polarizability with respect to nuclear coordinates are needed to determine infrared (IR) intensities and Raman activities, respectively. Their calculation is described in more detail in the following.

IR intensities#

In order to calculate intensities in the IR spectrum, one needs to know how the electric dipole moment μ=(μx,μy,μz) changes along a normal mode. This means that nuclear derivatives of μ have to be calculated, where each component μ can be decomposed into an electronic and a nuclear contribution, μ=μe+μn. The nuclear part of the dipole moment μn is simply given by the classical expression

μn=KZKRK

with the charge ZK and Cartesian coordinate RK of nucleus K. The electronic part μe is calculed quantum-mechanically from the dipole moment integrals in AO basis, μκλ and the one-particle density matrix D,

μe=κλDλκμκλ

The nuclear gradient of the dipole moment can again be calculated either numerically or analytically. Since the dipole moment itself as a first-order property can be considered a first derivative of the energy, its gradient corresponds to a mixed second derivative of the energy (once with respect to an electric field and once with respect to a nuclear coordinate, and thus a second-order property), and its calculation thus resembles the molecular Hessian calculation requiring the solution of a set of coupled-perturbed SCF (CPSCF) equations. While the analytic derivative of the nuclear contribution μn is trivial, the derivative of the electronic part μe with respect to a nuclear coordinate x is given by

dμedx=κλdDλκdxμκλ+κλDλκdμκλdx

where the perturbed density and the derivatives of the dipole integrals are needed.

The IR transition dipole moment is then calculated by taking the dot product of the dipole moment gradient with the Cartesian normal modes lCart, and the IR intensity as the square norm of the corresponding transition moment. The intensities are successively converted from atomic units to the unit of km mol1.

Raman activities#

Intensities in the vibrational Raman spectrum are calculated in an analogous manner, except that the nuclear derivative of the electric-dipole polarizability α(ω) along normal mode Qk is needed. Since the polarizability is already a second-order (or linear-response) property, this is more involved than IR intensities (in the time-independent regime, the polarizability gradient corresponds to a third derivative of the energy). However, it turns that out the analytical nuclear gradient of the polarizability can be calculated in manner analogous to the gradients of the TDHF or TDDFT schemes [RF07], which is similar to what has been discussed for here.

The main differences include the following:

  • all the density matrices (γ and Γ) depend also on two Cartesian components, m,n{x,y,z}, corresponding to the element αmn of the polarizability tensor α, and need to be symmetrized with respect to m and n.

  • the right-hand side of the linear response equation needs to be included in the Lagrangian, corresponding in this case to the contraction of response vector component m with component n of the dipole integrals (also to be symmetrized with respect to m and n).

  • in the dynamic case (ω0), also the term ω(XXYY) needs to be taken into account (for all combinations of vector components m and n).

For randomly oriented molecules and linearly polarized incident light, the Raman differential cross-section is calculated in practice as [Gut19]

dσkdΩ=ωLωS332π2ϵ02c4ωkSk45

where ωL is the frequency of the incident radiation, ωS is the angular frequency of the scattered light, ωk is the angular frequency of vibrational mode Qk, and Sk is the Raman activity of the mode, defined as [Gut19]

Sk=45α¯k2+7γ¯k2

where α¯k2 and γ¯k2 are Raman rotational invariants [Gut19]

α¯k2=19(dαxxdQk+dαyydQk+dαzzdQk)2γ¯k2=3[(dαxydQk+dαxzdQk+dαyzdQk)2]+12[(dαxxdQkdαyydQk)2+12(dαxxdQkdαzzdQk)2+12(dαyydQkdαzzdQk)2]

The rotational invariants can further be used to calculate the parallel (or “polarized”) and perpendicular (or “depolarized”) intensities as Ipol=45α¯k2+4γ¯k2 and Idepol=3γ¯k2, respectively, from which the so-called depolarization ratio ρ can be calculated as

ρ=IdepolIpol