Physics:Avoided crossing

Short description: Quantum physics and chemistry phenomenon

In quantum physics and quantum chemistry, an avoided crossing (sometimes called intended crossing, non-crossing or anticrossing) is the phenomenon where two eigenvalues of a Hermitian matrix representing a quantum observable and depending on N continuous real parameters cannot become equal in value ("cross") except on a manifold of N-3 dimensions. The phenomenon is also known as the von Neumann–Wigner theorem. In the case of a diatomic molecule (with one parameter, namely the bond length), this means that the eigenvalues cannot cross at all. In the case of a triatomic molecule, this means that the eigenvalues can coincide only at a single point (see conical intersection).

This is particularly important in quantum chemistry. In the Born–Oppenheimer approximation, the electronic molecular Hamiltonian is diagonalized on a set of distinct molecular geometries (the obtained eigenvalues are the values of the adiabatic potential energy surfaces). The geometries for which the potential energy surfaces are avoiding to cross are the locus where the Born–Oppenheimer approximation fails.

Avoided crossing also occur in the resonance frequencies of undamped mechanical systems, where the stiffness and mass matrices are real symmetric. There the resonance frequencies are the square root of the generalized eigenvalues.

In two-state systems

Emergence

Study of a two-level system is of vital importance in quantum mechanics because it embodies simplification of many of physically realizable systems. The effect of perturbation on a two-state system Hamiltonian is manifested through avoided crossings in the plot of individual energy vs energy difference curve of the eigenstates. The two-state Hamiltonian can be written as

[math]\displaystyle{ H= \begin{pmatrix}E_{1}&0\\0&E_{2}\end{pmatrix} \,\! }[/math]

The eigenvalues of which are [math]\displaystyle{ \textstyle E_{1} }[/math] and [math]\displaystyle{ \textstyle E_{2} }[/math] and the eigenvectors, [math]\displaystyle{ \textstyle \begin{pmatrix}1\\0\end{pmatrix} }[/math] and [math]\displaystyle{ \textstyle \begin{pmatrix}0\\1\end{pmatrix} }[/math]. These two eigenvectors designate the two states of the system. If the system is prepared in either of the states it would remain in that state. If [math]\displaystyle{ \textstyle E_{1} }[/math] happens to be equal to [math]\displaystyle{ E_{2} }[/math] there will be a twofold degeneracy in the Hamiltonian. In that case any superposition of the degenerate eigenstates is evidently another eigenstate of the Hamiltonian. Hence the system prepared in any state will remain in that forever.

Avoided crossing in two-state system. The energy level crossing is avoided with increasing the parameter [math]\displaystyle{ \textstyle w (= |W| ) }[/math]. In the absence of external perturbation the levels would have crossed if the original energy states were degenerate, i.e. [math]\displaystyle{ \textstyle \Delta E = 0 }[/math]

However, when subjected to an external perturbation, the matrix elements of the Hamiltonian change. For the sake of simplicity we consider a perturbation with only off diagonal elements. Since the overall Hamiltonian must be Hermitian we may simply write the new Hamiltonian

[math]\displaystyle{ H' = H + P= \begin{pmatrix}E_{1}&0\\0&E_{2}\end{pmatrix} + \begin{pmatrix}0&W\\W^{*}&0\end{pmatrix} = \begin{pmatrix}E_{1}&W\\W^{*}&E_{2}\end{pmatrix} \,\! }[/math]

Where P is the perturbation with zero diagonal terms. The fact that P is Hermitian fixes its off-diagonal components. The modified eigenstates can be found by diagonalising the modified Hamiltonian. It turns out that the new eigenvalues are,

[math]\displaystyle{ E_{+}=\frac{1}{2}(E_{1}+E_{2})+\frac{1}{2}\sqrt{(E_{1}-E_{2})^{2}+4|W|^{2}} }[/math]

[math]\displaystyle{ E_{-}=\frac{1}{2}(E_{1}+E_{2})-\frac{1}{2}\sqrt{(E_{1}-E_{2})^{2}+4|W|^{2}} }[/math]

If a graph is plotted varying [math]\displaystyle{ \textstyle (E_{1}-E_{2}) }[/math] along the horizontal axis and [math]\displaystyle{ \textstyle E_{+} }[/math] or [math]\displaystyle{ \textstyle E_{-} }[/math] along the vertical, we find two branches of a hyperbola (as shown in the figure). The curve asymptotically approaches the original unperturbed energy levels. Analyzing the curves it becomes evident that even if the original states were degenerate (i.e. [math]\displaystyle{ \textstyle E_{1}=E_{2} }[/math] ) the new energy states are no longer equal. However, if [math]\displaystyle{ \textstyle W }[/math] is set to zero we may find at [math]\displaystyle{ \textstyle (E_{1}-E_{2})=0 }[/math], [math]\displaystyle{ \textstyle E_{+}=E_{-} }[/math] and the levels cross. Thus with the effect of the perturbation these level crossings are avoided.

Quantum resonance

The immediate impact of avoided level crossing in a degenerate two state system is the emergence of a lowered energy eigenstate. The effective lowering of energy always correspond to increasing stability.(see: Energy minimization) Bond resonance in organic molecules exemplifies the occurrence of such avoided crossings. To describe these cases we may note that the non-diagonal elements in an erstwhile diagonalised Hamiltonian not only modify the energy eigenvalues but also superpose the old eigenstates into the new ones. These effects are more prominent if the original Hamiltonian had degeneracy. This superposition of eigenstates to attain more stability is precisely the phenomena of chemical bond resonance.

Our earlier treatment started by denoting the eigenvectors [math]\displaystyle{ \textstyle \begin{pmatrix}1\\0\end{pmatrix} }[/math] and [math]\displaystyle{ \textstyle \begin{pmatrix}0\\1\end{pmatrix} }[/math] as the matrix representation of eigenstates [math]\displaystyle{ \textstyle |\psi_{1} \rangle }[/math] and [math]\displaystyle{ \textstyle |\psi_{2} \rangle }[/math] of a two-state system. Using bra–ket notation the matrix elements of [math]\displaystyle{ H' }[/math] are actually the terms

[math]\displaystyle{ H'_{ij}=\langle \psi_{i}|H'|\psi_{j} \rangle }[/math] with [math]\displaystyle{ i,j \in \left\{ {1,2}\right\} }[/math]

where [math]\displaystyle{ H'_{11}=H'_{22}=E }[/math] due to the degeneracy of the unperturbed Hamiltonian and the off-diagonal perturbations are [math]\displaystyle{ H'_{12}=W }[/math] and [math]\displaystyle{ H'_{21}=W^{*} }[/math].

The new eigenstates [math]\displaystyle{ \textstyle |\psi_{+} \rangle }[/math] and [math]\displaystyle{ \textstyle |\psi_{-} \rangle }[/math] can be found by solving the eigenvalue equations [math]\displaystyle{ H'|\psi_{+}\rangle=E_{+}|\psi_{+}\rangle }[/math] and [math]\displaystyle{ H'|\psi_{-}\rangle=E_{-}|\psi_{-}\rangle }[/math]. From simple calculations it can be shown that

[math]\displaystyle{ |\psi_{+}\rangle = \frac{1}{\sqrt{2}}\begin{pmatrix}e^{i\phi}\\1\end{pmatrix}= \frac{1}{\sqrt{2}} (e^{i\phi}| \psi_{1}\rangle +|\psi_{2}\rangle) }[/math] and

[math]\displaystyle{ |\psi_{-}\rangle = \frac{1}{\sqrt{2}}\begin{pmatrix}-e^{i\phi}\\1\end{pmatrix}= \frac{1}{\sqrt{2}} (-e^{i\phi}| \psi_{1}\rangle +|\psi_{2}\rangle) }[/math] where [math]\displaystyle{ e^{i\phi}=W/|W| }[/math]

It is evident that both of the new eigenstates are superposition of the original degenerate eigenstates and one of the eigenvalues (here [math]\displaystyle{ E_{-} }[/math]) is less than the original unperturbed eigenenergy. So the corresponding stable system will naturally mix up the former unperturbed eigenstates to minimize its energy. In the example of benzene the experimental evidences of probable bond structures give rise to two different eigenstates, [math]\displaystyle{ \textstyle |\psi_{1} \rangle }[/math] and [math]\displaystyle{ \textstyle |\psi_{2} \rangle }[/math]. The symmetry of these two structures mandates that [math]\displaystyle{ \langle \psi_{1}|H|\psi_{1}\rangle=\langle \psi_{2}|H|\psi_{2}\rangle=E }[/math].

However it turns out that the two-state Hamiltonian [math]\displaystyle{ H }[/math] of benzene is not diagonal. The off-diagonal elements result into lowering of energy and the benzene molecule stabilizes in a structure which is a superposition of these symmetric ones with energy [math]\displaystyle{ E_{-}\lt E }[/math]. For any general two-state system avoided level crossing repels the eigenstates [math]\displaystyle{ |\psi_{+}\rangle }[/math] and [math]\displaystyle{ |\psi_{-}\rangle }[/math] such that it requires more energy for the system to achieve the higher energy configuration.

Resonances in avoided crossing

In molecules, the nonadiabatic couplings between two adiabatic potentials build the avoided crossing (AC) region. The rovibronic resonances in the AC region of two-coupled potentials are very special, since they are not in the bound state region of the adiabatic potentials, and they usually do not play important roles on the scatterings and are less discussed. Yu Kun Yang et al studied this problem in the New J. Phys. 22 (2020). Exemplified in particle scattering, resonances in the AC region are comprehensively investigated. The effects of resonances in the AC region on the scattering cross sections strongly depend on the nonadiabatic couplings of the system, it can be very significant as sharp peaks, or inconspicuous buried in the background. More importantly, it shows a simple quantity proposed by Zhu and Nakamura to classify the coupling strength of nonadiabatic interactions, can be well applied to quantitatively estimate the importance of resonances in the AC region.

General avoided crossing theorem

The above illustration of avoided crossing however is a very specific case. From a generalised view point the phenomenon of avoided crossing is actually controlled by the parameters behind the perturbation. For the most general perturbation [math]\displaystyle{ \textstyle P=\begin{pmatrix}W_{1}&W\\W&W_{2}\end{pmatrix} }[/math] affecting a two-dimensional subspace of the Hamiltonian [math]\displaystyle{ H }[/math], we may write the effective Hamiltonian matrix in that subspace as

[math]\displaystyle{ \begin{pmatrix}E_{1}&0\\0&E_{2}\end{pmatrix} + \begin{pmatrix}W_{1}&W\\W&W_{2}\end{pmatrix} =\begin{pmatrix}E_{1}+W_{1}&W\\W&E_{2}+W_{2}.\end{pmatrix}. }[/math]

Here the elements of the state vectors were chosen to be real so that all the matrix elements become real. Now the eigenvalues of the system for this subspace are given by

[math]\displaystyle{ E_{\pm}=\frac{1}{2}(E_{1}+E_{2}+W_{1}+W_{2}) \pm \frac{1}{2}\sqrt{(E_{1}-E_{2}+W_{1}-W_{2})^{2}+4W^{2}} }[/math]

The terms under the square root are squared real numbers. So for these two levels to cross we simultaneously require

[math]\displaystyle{ (E_{1}-E_{2}+W_{1}-W_{2})=0 }[/math]

[math]\displaystyle{ W=0. }[/math]

Now if the perturbation [math]\displaystyle{ P }[/math] has [math]\displaystyle{ k }[/math] parameters [math]\displaystyle{ { \alpha_{1},\alpha_{2},\alpha_{3}.....\alpha_{k} } }[/math] we may in general vary these numbers to satisfy these two equations.

[math]\displaystyle{ (E_{1}-E_{2}+W_{1}-W_{2})=F_{1}(\alpha_{1},\alpha_{2},\alpha_{3}.....\alpha_{k})=0 }[/math]

[math]\displaystyle{ W=F_{2}(\alpha_{1},\alpha_{2},\alpha_{3}.....\alpha_{k})=0. }[/math]

If we choose the values of [math]\displaystyle{ \alpha_{1} }[/math] to [math]\displaystyle{ \alpha_{k-1} }[/math] then both of the equations above have one single free parameter. In general it is not possible to find one [math]\displaystyle{ \alpha_{k} }[/math] such that both of the equations are satisfied. However, if we allow another parameter to be free, both of these two equations will now be controlled by the same two parameters

[math]\displaystyle{ F_{1}(\alpha_{k-1},\alpha_{k})|_{\alpha_{1},\alpha_{2},...,\alpha_{k-2} \, fixed}=0 }[/math]

[math]\displaystyle{ F_{2}(\alpha_{k-1},\alpha_{k})|_{\alpha_{1},\alpha_{2},...,\alpha_{k-2} \, fixed}=0. }[/math]

And generally there will be two such values of them for which the equations will be simultaneously satisfied. So with [math]\displaystyle{ k }[/math] distinct parameters [math]\displaystyle{ k-2 }[/math] parameters can always be chosen arbitrarily and still we can find two such [math]\displaystyle{ \alpha_{k} }[/math]s such that there would be crossing of energy eigenvalues. In other words, the values of [math]\displaystyle{ E_{+} }[/math] and [math]\displaystyle{ E_{-} }[/math] would be the same for [math]\displaystyle{ k-2 }[/math] freely varying co-ordinates (while the rest of the two co-ordinates are fixed from the condition equations). Geometrically the eigenvalue equations describe a surface in [math]\displaystyle{ k }[/math] dimensional space.

[math]\displaystyle{ E_{\pm}=E_{\pm}(\alpha_{1},\alpha_{2},\alpha_{3}.....\alpha_{k}). }[/math]

Since their intersection is parametrized by [math]\displaystyle{ k-2 }[/math] coordinates, we may formally state that for [math]\displaystyle{ k }[/math] continuous real parameters controlling the perturbed Hamiltonian, the levels (or surfaces) can only cross at a manifold of dimension [math]\displaystyle{ k-2 }[/math]. However the symmetry of the Hamiltonian has a role to play in the dimensionality. If the original Hamiltonian has asymmetric states, [math]\displaystyle{ \langle \psi_{1}|P|\psi_{2}\rangle \neq \langle \psi_{2}|P|\psi_{1}\rangle }[/math], the off-diagonal terms vanish automatically to ensure hermiticity. This allows us to get rid of the equation [math]\displaystyle{ W=0 }[/math]. Now from similar arguments as posed above, it is straightforward that for an asymmetrical Hamiltonian, the intersection of energy surfaces takes place in a manifold of dimension [math]\displaystyle{ k-1 }[/math].

In polyatomic molecules

In an N-atomic polyatomic molecule there are 3N-6 vibrational coordinates (3N-5 for a linear molecule) that enter into the electronic Hamiltonian as parameters. For a diatomic molecule there is only one such coordinate, the bond length r. Thus, due to the avoided crossing theorem, in a diatomic molecule we cannot have level crossings between electronic states of the same symmetry. However, for a polyatomic molecule there is more than one geometry parameter in the electronic Hamiltonian and level crossings between electronic states of the same symmetry are not avoided.

See also

• Geometric phase
• Christopher Longuet-Higgins
• Conical intersection
• Vibronic coupling
• Adiabatic theorem
• Bond hardening
• Bond softening
• Landau–Zener formula
• Level repulsion