---
title: Herzberg-Teller Approximation in Vibronic Spectroscopy
url: https://www.emergentmind.com/topics/herzberg-teller-approximation
type: topic
---

# Herzberg-Teller Approximation in Vibronic Spectroscopy

The Herzberg–Teller approximation is a first-order vibronic expansion of the electronic transition dipole moment with respect to nuclear coordinates. In place of the Condon approximation, which assumes a coordinate-independent dipole, it retains the leading linear dependence of the dipole on normal modes and thereby provides a systematic description of non-Condon effects, intensity borrowing, and vibronic activation of electronically weak or symmetry-forbidden transitions. In the cited treatments, it is central to vibrationally resolved absorption and emission spectroscopy, coherent multidimensional spectroscopy, and open-system descriptions of excitonic dynamics [2001.08414] [1812.04400] [2601.09346].

## 1. Formal definition and relation to the Condon approximation

In the Condon approximation, the electronic transition dipole moment is taken to be independent of nuclear geometry, for example
\[
\mu(Q)\approx \mu_0
\]
or, in the notation of one of the semiclassical formulations,
\[
\vec{\mu}(q)\approx \vec{\mu}(q_0).
\]
Under this assumption, vibronic intensities are controlled solely by Franck–Condon overlaps between vibrational states on the initial and final electronic surfaces [2001.08414] [1812.04400].

The Herzberg–Teller approximation relaxes that constraint by expanding the transition dipole in a Taylor series around a reference geometry and retaining the linear term. In normal-mode coordinates,
\[
\mu(Q)\approx \mu_0 + \sum_k \left(\frac{\partial \mu}{\partial Q_k}\right)_0 Q_k + \text{higher orders},
\]
and, equivalently, in generalized coordinates,
\[
\vec{\mu}(q)\approx \vec{\mu}(q_0)+\mathrm{grad}_q\vec{\mu}|_{q_0}^{\top}\cdot (q-q_0).
\]
The first derivative term is the Herzberg–Teller, or first-order non-Condon, contribution. Higher-order terms may also be written, but the cited works consistently emphasize the linear truncation as the operative HT approximation [2001.08414] [1812.04400] [1609.03685].

The physical content of the linear term is that vibronic motion modulates the electronic transition moment. This enables intensity borrowing from bright states, activates vibronic sidebands, and can make electronically forbidden transitions vibronically allowed when the zeroth-order dipole vanishes but its gradient does not [1812.04400] [1609.03685].

## 2. Vibronic transition moments, intensity borrowing, and selection rules

For a transition between vibrational states on two electronic surfaces, the vibronic transition moment in the first-order HT approximation is
\[
\mathcal{M}_{v'v''}
= \langle \chi_{v'}(e')| \mu(Q) |\chi_{v''}(e'')\rangle
\approx
\mu_0 \langle \chi_{v'}|\chi_{v''}\rangle
+
\sum_k \left(\frac{\partial \mu}{\partial Q_k}\right)_0
\langle \chi_{v'}|Q_k|\chi_{v''}\rangle.
\]
When the amplitude is squared, the intensity separates into a Franck–Condon term, a pure Herzberg–Teller term, and an FC–HT cross term,
\[
|\mathcal{M}|^2
\approx
|\mu_0\langle\chi_{v'}|\chi_{v''}\rangle|^2
+
\left|\sum_k \mu_k' \langle\chi_{v'}|Q_k|\chi_{v''}\rangle\right|^2
+
2\,\mathrm{Re}\!\left[
\mu_0\langle\chi_{v'}|\chi_{v''}\rangle
\left(\sum_k \mu_k' \langle\chi_{v'}|Q_k|\chi_{v''}\rangle\right)^{\!*}
\right].
\]
This decomposition makes explicit that HT terms both add new intensity channels and interfere with the FC contribution [2001.08414].

Because the coordinate operator connects different vibrational states, HT terms relax FC selection rules. In the harmonic-mode picture they activate transitions involving odd vibronic quanta and, more generally, transitions with nonzero matrix elements such as \(\langle \chi_e|Q_i|\chi_g\rangle\). The cited works therefore describe HT coupling as a route by which electronically weak or symmetry-forbidden bands gain measurable intensity [2205.06413] [2601.09346].

One exact identity used to clarify intensity borrowing relates the gradient of the transition dipole to nonadiabatic couplings and dipole couplings:
\[
\frac{\partial}{\partial q_i}\vec{\bm{\mu}}
=
[\vec{\bm{\mu}},\mathbf{F}_i]
+
\left(e\sum_{j=1}^N Z_j \frac{\partial \vec{R}_j}{\partial q_i}\right)\mathbf{1}.
\]
For off-diagonal transition elements near the ground-state equilibrium geometry, this reduces to the statement that a nonzero dipole gradient can arise from vibronic coupling to intermediate bright states. In that sense, the HT gradient is a direct formal expression of intensity borrowing [1812.04400] [1901.05208].

The symmetry analysis of pyrrole provides an explicit selection-rule example. For the \(1^1A_2(\pi\sigma^*) \leftarrow \tilde{X}^1A_1(\pi\pi)\) transition in \(C_{2v}\), the direct dipole moment vanishes by symmetry, but first-order HT terms activate the transition when the dipole symmetry and vibrational symmetry combine to contain \(A_2\). In the notation given there, \(B_1\) vibrational distortions activate \(\mu_y\), \(B_2\) distortions activate \(\mu_x\), and \(A_2\) ring modes activate \(\mu_z\) [1712.02339].

## 3. Time-dependent spectral formulations

A major virtue of the HT approximation is that it enters naturally in time-domain correlation-function formulations of spectra. For zero-temperature absorption, one cited formulation writes
\[
\sigma_{\mathrm{abs}}(\epsilon,\omega)
=
\frac{4\pi \omega}{\hbar c}
\mathrm{Re}\int_0^\infty dt\, C(\epsilon,t)\, e^{i(\omega+\omega_{1,g})t},
\]
with
\[
C(\epsilon,t)=\langle \phi(0)|\phi(t)\rangle,\qquad
|\phi(0)\rangle=\hat{\mu}|1,g\rangle,\qquad
|\phi(t)\rangle=e^{-i\hat{H}_2 t/\hbar}|\phi(0)\rangle.
\]
For emission from \(|2,g\rangle\), an analogous expression propagates the initial dipole-prepared state on the ground-state surface [2001.08414].

In the Condon limit, \(|\phi(0)\rangle\propto \mu_0 \psi(0)\), so the autocorrelation reflects only nuclear dynamics on the final surface. With HT coupling,
\[
|\phi(0)\rangle \propto [\mu_0+\mu'^{T}\hat{Q}]\,\psi(0),
\]
and the correlation function acquires pure HT and FC–HT cross contributions involving coordinate operators at \(t=0\) and during propagation. This is the basic time-domain origin of non-Condon structure in vibronic spectra [2001.08414].

The same viewpoint is used in an explicitly tensorial orientational average beyond the Condon approximation. For coordinate-dependent dipoles, the isotropic average is written as
\[
\overline{C_{\mu\mu}(t)}=\frac{1}{3}\mathrm{Tr}[C_{\vec{\mu}\vec{\mu}}(t)],
\qquad
\overline{\sigma(\omega)}=\frac{1}{3}\mathrm{Tr}[\sigma_{\vec{\mu}\vec{\mu}}(\omega)].
\]
Under Condon conditions this reduces to the textbook \((1/3)\) rule based on \(|\vec{\mu}|^2\), but the tensor formulation remains valid in the HT case [1812.04400].

At finite temperature, the correlation function becomes
\[
C(t)=\mathrm{Tr}\!\left(\hat{\mu}^\dagger e^{-i\hat{H}_2 t/\hbar}\hat{\mu}\hat{\rho}e^{i\hat{H}_1 t/\hbar}\right),
\]
and can be mapped, via thermo-field dynamics, to an autocorrelation in a doubled space,
\[
C(t)=\langle \bar{\phi}_0|\bar{\phi}_t\rangle,
\qquad
i\hbar \dot{\bar{\phi}}_t(\bar{q})=\bar{H}(\bar{q})\bar{\phi}_t(\bar{q}),
\]
with \(\bar{H}=H_2(q)-H_1(q')\). In the cited implementation, the HT dipole acts only on the physical coordinates, while the thermal purification is carried by the doubled Gaussian initial state [2005.09126].

In coherent multidimensional spectroscopy, the same non-Condon mechanism appears as multiplicative oscillatory factors in the time-domain response. For the displaced harmonic oscillator model with linear HT coupling,
\[
\mu = \mu_0 \otimes \mathcal{I} + \mu_1 \otimes (a+a^\dagger),
\]
the response functions factor into an electronic term, a Franck–Condon vibrational factor, and finite sums of oscillatory HT corrections. Upon Fourier transformation, these oscillations generate replicas of FC multidimensional spectral features shifted by integer multiples of the vibrational frequency [2601.09346].

## 4. Semiclassical, wavepacket, and ab initio realizations

Several cited works implement HT effects through wavepacket propagation. In the thawed Gaussian approximation, the nuclear wavepacket is propagated in a local harmonic approximation,
\[
\psi_t(q)=N_0 \exp\!\left\{-(q-q_t)^{\top}A_t(q-q_t)+\frac{i}{\hbar}[p_t^{\top}(q-q_t)+\gamma_t]\right\},
\]
with center, width, and phase evolving according to the local gradient and Hessian of the potential [1812.04400] [1901.05208].

The extended thawed Gaussian approximation propagates a Gaussian multiplied by a linear polynomial,
\[
\phi_t(q)=[a_t+b_t^{\top}(q-q_t)]\psi_t(q),
\]
or, in the notation used for non-Condon spectra of azulene,
\[
\phi(q,t)=\big[\mu(q_0)+\mu'(q_0)^{T}Q_0Q_t^{-1}(q-q_t)\big]\psi(q,t).
\]
Within the local harmonic potential this functional form is preserved, so a single Gaussian dynamics yields both FC and HT contributions analytically and efficiently [2001.08414] [1812.04400].

To reduce cost further, one implementation replaces the time-dependent Hessian in the width equations by a fixed adiabatic reference Hessian evaluated once at a final-state optimized geometry. This “single-Hessian” construction keeps the Gaussian center on the true force field while using constant curvature for widths. In the cited azulene study, the HT propagation formula remains valid under this approximation, and the method treats anharmonicity and mode mixing efficiently through a single on-the-fly ab initio trajectory [2001.08414].

A related variant, the three thawed Gaussians approximation, represents the initial HT contribution not as a linear prefactor on one Gaussian but as an antisymmetric superposition of two displaced Gaussians, with a third Gaussian carrying the Condon term. The motivation stated there is that the HT packet has a node at the center and can be viewed as split into lobes in anharmonic potentials. Each Gaussian is propagated independently, and the spectra are assembled from their coherent superposition [1901.05208].

For finite temperature, the doubled-space ETGA carries over with nearly no additional computational cost because the ground-state branch is typically stationary and the formal doubling does not require new on-the-fly electronic-structure calculations. The cited benzene application uses this framework to separate finite-temperature, anharmonicity, and HT effects in a symmetry-forbidden spectrum [2005.09126].

The wavepacket perspective also underlies the multiconfigurational time-dependent Hartree treatment of pyrrole. There, coordinate-dependent transition dipole functions constructed from a linear HT expansion prepare the initial wavepacket, and the subsequent propagation on a locally diabatic Hamiltonian yields frequency-resolved absorption spectra and dissociation time scales for an optically dark \(\pi\sigma^*\) state [1712.02339].

## 5. Open-system and excitonic formulations

In excitonic open-system models, HT coupling is often represented by promoting the transition dipole to an operator that depends on bath coordinates. One cited DEOM formulation writes
\[
\hat{V}=\sum_a \mu_a(X_a)\hat{D}_a,
\qquad
\mu_a(X_a)=\mu_a+\mu_a' \hat{F}_a/\lambda,
\]
where \(\hat{F}_a\) is the collective bath operator of an underdamped Brownian oscillator mode. In the Franck–Condon limit \(\mu_a'=0\); with \(\mu_a'\neq 0\), the field couples directly to the vibrational coordinate during electronic excitation [1609.03685].

For linear absorption of a two-level monomer in that framework, the cited analytical result is
\[
I(\omega)=\big(\mu+\mu' \Delta\omega/\lambda\big)^2 I_0(\omega).
\]
The cross term and quadratic term explicitly redistribute intensity: increasing \(\mu'\) suppresses the zero-phonon peak and enhances vibronic sidebands, especially on the blue side. The same study emphasizes that non-Condon coupling intensifies dynamical electronic–vibrational energy transfer, enhances total system-and-bath quantum coherence, and strengthens long-lived vibronic oscillations in 2D spectroscopy [1609.03685].

A closely related open-system construction in CODDE and DEOM uses a total dipole
\[
\hat{\mu}_{\mathrm{tot}}
=
\hat{\mu}\Big(1+\sum_a \nu_a \hat{F}_a\Big),
\]
which is the linear HT form expressed in system–bath language. In that formalism, non-Condon response functions are obtained through one-dissipaton auxiliary variables and correlated driving kernels, and the linear absorption correlation function contains both the reduced density operator and first-tier hybrid bath objects [2205.06413].

HEOM treatments extend the same idea to both coordinate-dependent transition dipoles and coordinate-dependent interstate couplings. The cited dimer Hamiltonian includes first- and second-order HT terms in the resonance coupling,
\[
\hat{H}^{(1)}=\sum_{m}\sum_{n\neq m} J_{mn}^{(1)}(\hat{q}_m+\hat{q}_n)\hat{B}_m^\dagger \hat{B}_n,
\]
\[
\hat{H}^{(2)}=\sum_{m}\sum_{n\neq m}\frac{1}{2}J_{mn}^{(2)}(\hat{q}_m+\hat{q}_n)^2\hat{B}_m^\dagger \hat{B}_n,
\]
as well as a first-order non-Condon dipole
\[
\hat{\boldsymbol{\mu}}^{(1)}
=
\sum_m \boldsymbol{\mu}_m^{(1)} \hat{q}_m(\hat{B}_m^\dagger+\hat{B}_m).
\]
In HEOM space these terms introduce additional couplings between adjacent auxiliary density operators, and for sufficiently strong excitonic coupling the resulting linear absorption spectra can be captured very well by effective Huang–Rhys factors [1807.07475].

## 6. Representative molecular applications and spectroscopic signatures

The cited applications show that HT effects are highly system dependent. For the phenyl radical, on-the-fly semiclassical spectra with and without HT are very similar and the HT contribution only slightly broadens peaks; the conclusion stated there is that anharmonicity outweighs HT corrections. By contrast, for benzene’s symmetry-forbidden \(\tilde{\mathrm{A}}^1\mathrm{B}_{2u}\leftarrow \tilde{\mathrm{X}}^1\mathrm{A}_{1g}\) transition, the Condon spectrum is zero because \(\vec{\mu}(q_0)=0\), whereas the HT spectrum agrees well with experiment and reproduces intensities of the main progression within about \(5\%\) of experimental values [1812.04400].

Azulene provides a different non-Condon regime. In the cited study, \(S_1\leftarrow S_0\) absorption is FC-dominated and HT does not improve the band appreciably, while \(S_2\leftarrow S_0\) absorption and \(S_2\rightarrow S_0\) emission require HT to reproduce the experimental vibronic structure and intensity. The same work states that the breaking of absorption–emission mirror symmetry for \(S_2\) arises predominantly from HT terms interacting with Duschinsky mode mixing, and that turning off either HT or mode mixing largely restores symmetry [2001.08414].

Pyrrole illustrates HT-enabled access to dissociative dark states. The \(1^1A_2(\pi\sigma^*) \leftarrow \tilde{X}^1A_1(\pi\pi)\) transition is electric-dipole forbidden at the \(C_{2v}\) equilibrium geometry, but becomes weakly allowed through vibronic coupling. In the 15-dimensional treatment cited there, the total absorption maximum is at \(4.70\ \mathrm{eV}\), the envelope width is \(0.61\ \mathrm{eV}\), and the polarization-resolved intensities satisfy \(\sigma_y > \sigma_z > \sigma_x\). The spectrum is diffuse, structured by ring-mode progressions together with obligatory H-bend excitation, and the associated photodissociation is predominantly direct with short survival times [1712.02339].

In strongly coupled vibronic dimers treated by HEOM, first-order HT coupling enhances the vibrational progression in one exciton band and diminishes it in the other, while non-Condon dipoles redistribute intensities without changing vibronic spacings. The same study reports that effective Huang–Rhys factors provide a compact description of these HT-induced spectral changes when excitonic coupling is sufficiently strong [1807.07475].

Across these applications, a recurring diagnostic pattern is that HT terms become decisive when a transition is electronically forbidden or weak at the reference geometry, when absorption–emission mirror symmetry is broken, or when vibronic sidebands and intensity patterns cannot be reproduced within FC-only models [1812.04400] [2001.08414].

## 7. Validity regime, approximations, and practical use

All cited implementations emphasize that the standard HT approximation is a truncation at first order in the dipole expansion. It is accurate when the transition dipole varies moderately over the support of the nuclear wavepacket. Strongly nonlinear \(\mu(Q)\) requires higher orders, and several works explicitly note that second-order vibronic coupling or quadratic HT terms are omitted from the first-order treatments [2001.08414] [1812.04400] [2005.09126].

The wavepacket and semiclassical realizations impose further assumptions. The ETGA and related thawed-Gaussian methods assume a locally harmonic potential around the trajectory center, and single-trajectory versions cannot describe strong wavepacket splitting or strongly anharmonic long-time dynamics. The single-Hessian variant is most reliable in near-harmonic regions, although it still captures part of the anharmonicity through propagation of the packet center on the full ab initio surface [2001.08414] [1812.04400].

Open-system approaches impose different restrictions. The DEOM and CODDE formulations use harmonic baths and Gaussian statistics; the HEOM implementation discussed for HT and non-Condon effects relies on exponential decompositions of bath correlation functions and requires increasing hierarchy depth for convergence, especially for undamped oscillators and second-order HT terms [1609.03685] [2205.06413] [1807.07475].

Within these limits, the practical guidance given in the cited studies is consistent. HT terms should be included when electronically weak or symmetry-forbidden bands are of interest, when mirror-image symmetry between absorption and emission is broken, or when nearby bright states and vibronic couplings suggest strong intensity borrowing. Mode mixing should be treated accurately, \(\mu_0\) and its gradient should be evaluated at the appropriate reference geometry, and in single-surface spectral simulations the Born–Oppenheimer assumption should be validated independently when nonadiabatic effects may be significant [2001.08414] [1812.04400].

Source: https://www.emergentmind.com/topics/herzberg-teller-approximation