Vibronic coupling (also called nonadiabatic coupling or derivative coupling) in a molecule involves the interaction between electronic and nuclear vibrational motion. The term "vibronic" originates from the combination of the terms "vibrational" and "electronic", denoting the idea that in a molecule, vibrational and electronic interactions are interrelated and influence each other. The magnitude of vibronic coupling reflects the degree of such interrelation. In theoretical chemistry, the vibronic coupling is neglected within the Born–Oppenheimer approximation. Vibronic couplings are crucial to the understanding of nonadiabatic processes, especially near points of conical intersections. The direct calculation of vibronic couplings used to be uncommon due to difficulties associated with its evaluation, but has recently gained popularity due to increased interest in the quantitative prediction of internal conversion rates, as well as the development of cheap but rigorous ways to analytically calculate the vibronic couplings, especially at the TDDFT level.
Definition Vibronic coupling describes the mixing of different electronic states as a result of small vibrations.
f k ′ k ≡ ⟨ χ k ′ ( r ; R ) | ∇ ^ R χ k ( r ; R ) ⟩ ( r ) {\displaystyle \mathbf {f} _{k'k}\equiv \langle \,\chi _{k'}(\mathbf {r} ;\mathbf {R} )\,|\,{\hat {\nabla }}_{\mathbf {R} }\chi _{k}(\mathbf {r} ;\mathbf {R} )\rangle _{(\mathbf {r} )}}
Evaluation The evaluation of vibronic coupling often involves complex mathematical treatment.
Numerical gradients The form of vibronic coupling is essentially the derivative of the wave function. Each component of the vibronic coupling vector can be calculated with numerical differentiation methods using wave functions at displaced geometries. This is the procedure used in MOLPRO. First order accuracy can be achieved with forward difference formula:
( f k ′ k ) l ≈ 1 d [ γ k ′ k ( R | R + d e l ) − γ k ′ k ( R | R ) ] {\displaystyle (\mathbf {f} _{k'k})_{l}\approx {\frac {1}{d}}\left[\gamma ^{k'k}(\mathbf {R} |\mathbf {R} +d\mathbf {e} _{l})-\gamma ^{k'k}(\mathbf {R} |\mathbf {R} )\right]}
Second order accuracy can be achieved with central difference formula:
( f k ′ k ) l ≈ 1 2 d [ γ k ′ k ( R | R + d e l ) − γ k ′ k ( R | R − d e l ) ] {\displaystyle (\mathbf {f} _{k'k})_{l}\approx {\frac {1}{2d}}\left[\gamma ^{k'k}(\mathbf {R} |\mathbf {R} +d\mathbf {e} _{l})-\gamma ^{k'k}(\mathbf {R} |\mathbf {R} -d\mathbf {e} _{l})\right]}
Here, e l {\displaystyle \mathbf {e} _{l}} is a unit vector along direction l {\displaystyle l} . γ k ′ k {\displaystyle \gamma ^{k'k}} is the transition density between the two electronic states.
… excerpt ends here. Continue reading the full article.
