---
title: Qx–Qy Vibronic Mixing
url: https://www.emergentmind.com/topics/qx-qy-vibronic-mixing
type: topic
---

# Qx–Qy Vibronic Mixing

Searching arXiv for the specified papers and closely related work on Qx–Qy vibronic mixing.
Qx–Qy vibronic mixing denotes the entanglement of two orthogonal degrees of freedom through nuclear motion, but the meaning of “Qx” and “Qy” depends on the literature. In Jahn–Teller and pseudo–Jahn–Teller treatments of solids, $Q_x$ and $Q_y$ are usually the two components of an $E_g$ vibrational doublet, and their direction in the distortion plane fixes the admixture of an electronic doublet. In porphyrins, chlorins, and related aggregates, by contrast, Qx and Qy usually denote two orthogonally polarized electronic transitions within the Q band, and vibronic mixing refers to nuclear-coordinate-dependent coupling between those electronic manifolds. Across both usages, the common structure is a two-component electronic sector, a symmetry-selected vibrational sector, and an adiabatic potential-energy surface whose topology controls intensity borrowing, pseudorotation, order–disorder behavior, or energy-flow pathways [1510.06671] [2409.08095] [2105.04725] [2507.11007].

## 1. Terminological scope and definitions

The term is not tied to a single physical system. It designates a pattern of coupling in which two orthogonal components become mixed by vibronic interaction, but those components may be vibrational coordinates, electronic basis states, or both.

| Context | Meaning of $Q_x,Q_y$ | Character of mixing |
|---|---|---|
| CuInP$_2$S$_6$ | $Q_x \equiv Q_1$, $Q_y \equiv Q_2$ of a $\Gamma_5(E_g)$ mode | An $E_g$ distortion mixes a $\Gamma_5$ electronic doublet |
| Ba$_2$CaReO$_6$ | $Q_x \equiv q_{x^2-y^2}$, $Q_y \equiv q_{z^2}$ | An $E_g$ Jahn–Teller mode mixes pseudoorbital states in the SOC-adapted manifold |
| LHCII and porphyrin aggregates | Qx and Qy are orthogonally polarized electronic transitions | Herzberg–Teller activity and vibronic resonance mix Qx and Qy manifolds |
| SrOPh / SrOPh-d$_5$ | No true Qx–Qy degeneracy | Two distinct non-degenerate bends generate a second-order analogue rather than degenerate-mode Qx–Qy mixing |

In CuInP$_2$S$_6$, the near-gap copper $d$-derived states form an elementary energy band whose $\Gamma$-point doublets transform as $\Gamma_5(E_g)$ and $\Gamma_6(E_u)$ of $D_{3d}$, and the active vibrational mode is likewise $\Gamma_5(E_g)$. In Ba$_2$CaReO$_6$, the Jahn–Teller-active $E_g$ coordinates are written as $q_{x^2-y^2}$ and $q_{z^2}$, which can be mapped directly to $Q_x$ and $Q_y$. In LHCII and porphyrin nanotubes, Qx and Qy are the two orthogonally polarized transitions of the Q band, and the mixing is assigned to Herzberg–Teller activity rather than to a purely Franck–Condon mechanism. In SrOPh and SrOPh-d$_5$, the analogy is explicitly limited: the in-plane and out-of-plane bends are distinct $b_2$ and $b_1$ modes, so the observed state mixing is not a true degenerate Qx–Qy problem [1510.06671] [2409.08095] [2105.04725] [2510.22388].

## 2. Symmetry structure and model Hamiltonians

In the trigonal protostructure of CuInP$_2$S$_6$, the $\Gamma_5$ electronic doublet is represented in an orthonormal basis $\{|e_1\rangle,|e_2\rangle\}$, with Pauli matrices $\tau_x$ and $\tau_z$. The symmetry content of the symmetric square is
$$
\Gamma_5 \otimes_s \Gamma_5 = \Gamma_1 \oplus \Gamma_5,
$$
so the scalar $\Gamma_1$ contributes the harmonic term and the $\Gamma_5$ component contributes the linear and quadratic anisotropy terms. The resulting vibronic matrix is
$$
H_{\mathrm{vib}}(Q_x,Q_y)=\frac{1}{2}\omega^2(Q_x^2+Q_y^2)I+V(Q_x\tau_x+Q_y\tau_z)+W[(Q_x^2-Q_y^2)\tau_z+(Q_xQ_y)\tau_x].
$$
This is a crystal realization of the linear $E\otimes e$ problem supplemented by symmetry-allowed warping terms, and it is explicitly tied to the elementary energy band induced from the copper Wyckoff position $d(2/3,1/3,1/4)$ [1510.06671].

In Ba$_2$CaReO$_6$, the onsite model is written as
$$
H = H_{\mathrm{lat}} + H_{\mathrm{SOC}} + V_{\mathrm{JT}} + V_{\mathrm{vib}},
$$
with
$$
H_{\mathrm{lat}}=\frac{\hbar\omega}{2}\left(p_{x^2-y^2}^2+q_{x^2-y^2}^2+p_{z^2}^2+q_{z^2}^2\right),\qquad
H_{\mathrm{SOC}}=\lambda\,\mathbf l\cdot \mathbf s,
$$
and a linear $E_g$ Jahn–Teller interaction. In the standard $E\otimes e$ form, the coupling is written as $H_{\mathrm{JT}}=\tilde g(Q_x\tau_x+Q_y\tau_z)$. The authors use $\lambda=0.32$ eV and $\hbar\omega=66$ meV, and emphasize that SOC quenches the Jahn–Teller coupling in the $j_{\mathrm{eff}}=3/2$ states by half [2409.08095].

For molecular excited-state manifolds, a broader linear vibronic coupling representation writes
$$
W(\mathbf Q)=\sum_i\left(E_i^0+\sum_\alpha \kappa_{i\alpha}Q_\alpha\right)|i\rangle\langle i|
+\sum_{i\neq j}\sum_\alpha \lambda_{ij,\alpha}Q_\alpha |i\rangle\langle j|.
$$
Here $\kappa_{i\alpha}$ are intrastate couplings and $\lambda_{ij,\alpha}$ are interstate vibronic couplings. The mode-selection rule is
$$
\Gamma_n\otimes \Gamma_Q\otimes \Gamma_m \supset \Gamma_A,
$$
which separates gap-tuning totally symmetric modes from non-totally symmetric modes that provide off-diagonal interstate mixing [1803.11360].

## 3. Adiabatic potentials, mixing angles, and geometric phase

Diagonalization of the CuInP$_2$S$_6$ vibronic matrix gives
$$
E_\pm(Q_x,Q_y)=\frac{1}{2}\omega^2(Q_x^2+Q_y^2)\pm
\sqrt{[VQ_y+W(Q_x^2-Q_y^2)]^2+[VQ_x+W(Q_xQ_y)]^2},
$$
with electronic mixing angle
$$
\tan[2\chi(Q_x,Q_y)] =
\frac{VQ_x+W(Q_xQ_y)}{VQ_y+W(Q_x^2-Q_y^2)}.
$$
In the linear limit, the lower sheet has the circular trough characteristic of $E\otimes e$ coupling. The quadratic term $W$ warps that trough and, in the $D_{3d}$ protostructure, the lower adiabatic surface acquires six equivalent minima with $60^\circ$ angular spacing. In the monoclinic paraphase, by contrast, the active vibronic instability is reduced to a one-dimensional double-well,
$$
\epsilon(Q)=\frac{1}{2}Q^2\pm(\Delta^2+V^2Q^2)^{1/2},
$$
and in the ferriphase further symmetry lowering selects one well, which is the basis for the discussion of spontaneous polarization [1510.06671].

In Ba$_2$CaReO$_6$, the adiabatic $E\otimes e$ form is
$$
E_\pm(Q)=\frac{K}{2}Q^2\pm g_e Q,\qquad Q=\sqrt{Q_x^2+Q_y^2},
$$
and the local pseudoorbital orientation is fixed by
$$
\tan\theta=\frac{Q_y}{Q_x}.
$$
The ideal linear problem has a circular trough at $Q_0=g_e/K$, but the work distinguishes three regimes: a strongly warped classical limit without SOC, a shallow almost circular trough with strong SOC, and an anisotropy-lifted trough under static elastic fields. The Berry phase around the conical intersection imposes antiperiodic boundary conditions and gives the vibronic sequence $\Gamma_8,\Gamma_7,\Gamma_6,\dots$ [2409.08095].

The same geometry reappears in the explicit $E\otimes e$ treatment of cavity-coupled Jahn–Teller molecules, where
$$
Q_x=\rho\cos\theta,\qquad Q_y=\rho\sin\theta,
$$
and the lower adiabatic surface is
$$
E_-(\rho)=\frac{1}{2}M\omega^2\rho^2-g\rho.
$$
This produces the conical intersection at $\rho=0$, the Mexican-hat minimum at $\rho_0=g/(M\omega^2)$, and half-integer vibronic angular momentum associated with pseudorotation. This suggests that Qx–Qy vibronic mixing is governed not only by coupling strength but also by the topology of the adiabatic surface [2511.07880].

## 4. Crystalline and oxide realizations

In CuInP$_2$S$_6$, Qx–Qy vibronic mixing is embedded in an order–disorder phase-transition scenario. The near-gap elementary energy band is induced from the copper site-symmetry group, and the copper $d$ states at the valence-band top provide the $\Gamma_5$ and $\Gamma_6$ doublets relevant for the Jahn–Teller analysis. The trigonal protostructure supports six equivalent minima on the lower adiabatic surface, the monoclinic paraphase reduces this to a two-minimum double-well associated with order–disorder between Cu1u and Cu1d sites, and the ferriphase selects one minimum and breaks inversion symmetry. The discussion connects this sequence directly to the possibility that spontaneous polarization arises normal to the layers in CuInP$_2$S$_6$ [1510.06671].

In Ba$_2$CaReO$_6$, Iwahara et al. identify spectroscopic signatures that cannot be reproduced by a purely classical lattice description. At the Re $L_3$ edge, the main inelastic feature appears at $\Delta E\approx0.52$ eV and is accompanied by a high-energy tail and vibronic satellites at $0.59$, $0.67$, and $0.73$ eV. At the O $K$ edge, a distinct low-energy peak occurs at $\Delta E\approx0.10$ eV with a broad tail up to $\sim0.4$ eV. These features persist across the structural transition at $T_s\approx130$ K. The analysis attributes them to orbital–lattice entangled states generated by Qx–Qy vibronic mixing within the $j_{\mathrm{eff}}=3/2$ manifold. The same work reports dynamic vibronic stabilization of approximately $30$ meV, static Jahn–Teller stabilization of approximately $14$ meV, and cooperative elastic fields of order several meV to $\sim10$ meV, concluding that the dynamic effect remains dominant even below the onset of multipolar order [2409.08095].

## 5. Molecular and aggregate realizations

In photosynthetic light harvesting, Qx–Qy vibronic mixing refers to coupling between orthogonally polarized electronic transitions rather than to a doubly degenerate normal mode. In LHCII, the relevant mechanism is assigned to Herzberg–Teller activity. The note on LHCII identifies a higher-energy side-band in linear absorption and 2DEV spectra that appears only with Herzberg–Teller activity in the model, reproduces higher-lying excitonic states composed mainly of chlorophyll $b$ in the $15200$–$15600$ cm$^{-1}$ range, and is tracked through a ground-state bleach at $1690$ cm$^{-1}$. Follow-up work cited there attributes green-light absorption across nearly a $3000$ cm$^{-1}$ span of the electronic spectrum to mixed vibronic Qy–Qx states. The note therefore revises earlier interpretations by assigning the dominant role to Herzberg–Teller rather than realistic Franck–Condon-only coupling [2105.04725].

In porphyrin nanotube aggregates, polarization-controlled 2DES resolves mixed Qx–Qy pathways at room temperature. The cross-peak-selective polarization sequence $(0^\circ,0^\circ,+60^\circ,-60^\circ)$ removes isotropic contributions and suppresses purely Qx–Qx or Qy–Qy pathways. Cross-peaks between the main Q band and vibronic shoulders are already present at $T=30$ fs, the lower cross-peak is about $10\times$ stronger than the upper one at that time, and the ratio $\mathrm{DP}_L/\mathrm{CP}_L$ changes from $3.2$ in all-parallel 2D to $0.8$ in the cross-peak-selective experiment. Low-frequency beats at $236$ and $309$ cm$^{-1}$ are removed by the polarization filter, whereas anisotropic beats at $432$ and $880$ cm$^{-1}$ survive and are assigned to out-of-plane pyrrole-ring deformations that modulate Qx–Qy mixing. Static energetic disorder with $\sigma\approx200$ cm$^{-1}$ spreads Qx character across the Q band and increases simulated $0$–$1$ mixing from approximately $5\%$ to $19\%$ for the $880$ cm$^{-1}$ mode and from approximately $20\%$ to $55\%$ for the $440$ cm$^{-1}$ mode [2507.11007].

A more restrictive conclusion emerges from a heterodimer study of vibronic excitons in photosynthetic energy transfer. There, a low-frequency intramolecular vibration with $\omega_v=180$ cm$^{-1}$ and $S=0.025$ generates robust long-lived beating at cryogenic temperature, but environmentally induced fluctuations eradicate the mixing at physiological temperature. The same work argues that such electronic–vibrational quantum mixtures do not necessarily play a significant role in electronic energy-transfer dynamics even when they enhance long-lived beating in 2D spectra. This suggests that the spectroscopic visibility of Qx–Qy vibronic mixing and its dynamical functionality need not coincide [1505.05281].

## 6. Computational parameterization, diagnostics, and limits of the label

One route to quantitative parameterization is the overlap-based extraction of interstate vibronic coupling constants. For a displaced geometry along mode $Q_i$, the protocol computes
$$
S_i^{(nm)}=\langle \Psi_n(0)|\Psi_m(\delta Q_i)\rangle,
$$
constructs a Löwdin-orthogonalized transformation $U$, forms
$$
W(\delta Q_i)=U^{\mathrm T}\,\mathrm{diag}[E_n(\delta Q_i)]\,U,
$$
and estimates the linear coupling through
$$
\lambda_i^{(nm)}\approx \frac{W_{nm}(\delta Q_i)}{\delta Q_i}.
$$
Because it uses overlaps between auxiliary excited-state wavefunctions rather than only excited-state energies, this protocol was proposed as an alternative to Hessian-based extraction in dense excited-state manifolds [1803.11360].

For aggregates, the dimensionality of the vibronic problem can be reduced by exact construction of aggregate normal modes within each equal-frequency set. In the trimer example, the orthonormal coordinates
$$
Q_k^+=\frac{q_{Ak}+q_{Bk}+q_{Ck}}{\sqrt3},\qquad
Q_k^-=\frac{q_{Ak}-2q_{Bk}+q_{Ck}}{\sqrt6},\qquad
Q_k^{AC}=\frac{q_{Ak}-q_{Ck}}{\sqrt2}
$$
separate a spectator mode from promoter modes that tune energy gaps. The mode $Q_k^-$ tunes the $AB$ and $BC$ gaps simultaneously, whereas $Q_k^+$ does not tune any pairwise gap. This reduction was used to show how vibronic resonance can convert an electronically uncoupled “trap” into a “conduit” for transport, and it provides a mode-resolved language for identifying which coordinates actually promote Qx–Qy mixing in aggregates [2011.03717].

The spectroscopic diagnosis of vibronic mixing in 2D coherent spectroscopy is likewise symmetry- and resonance-sensitive. In the Holstein-type dimer analysis, off-resonant cases preserve a clean distinction between electronic beats and vibrational beats, while near-resonant cases with $\omega_0\approx2J$ generate avoided crossings, multiple oscillation frequencies, and congested Fourier maps. The effective two-level description uses
$$
\tan(2\theta)=\frac{2V_{\mathrm{eff}}}{\Delta_{\mathrm{eff}}},
$$
and complete exciton–vibronic mixing occurs when $|2V_{\mathrm{eff}}|\gtrsim|\Delta_{\mathrm{eff}}|$. The same work formulates vibrational lifetime borrowing as the mechanism by which mixed coherences can outlive purely electronic ones [1310.1343].

A recurrent terminological limit is that not every orthogonal two-mode interaction is a genuine Qx–Qy problem. In SrOPh and SrOPh-d$_5$, the observed mixing is between $B,v0$ and the combination level $A,v21\,v33$, with $v33$ a $b_2$ in-plane bend and $v21$ a $b_1$ out-of-plane bend. The effective $A$–$B$ coupling is deduced to be about $0.5$ cm$^{-1}$, the direct linear $a_2$ coupling is vanishingly small, and the dominant pathway is second order through the higher $C(A_1)$ state. The work explicitly states that this is not Qx–Qy degenerate-mode mixing. This clarification is important because it separates true two-component vibronic angular-momentum physics from combination-mode symmetry pathways that only resemble it at a coarse level [2510.22388].

Source: https://www.emergentmind.com/topics/qx-qy-vibronic-mixing