⚛️ quantum physics

Chiral Torsional Double Wells: Tunneling and Polarizability

This study demonstrates that while quantum tunneling rates between chiral conformers in hydrogen peroxide, disulfide, and diselenide decrease significantly across the series, the intrinsic transition polarizability remains robust or even increases because the two properties depend on different regions of the torsional potential landscape.

Original authors: Shah Ishmam Mohtashim, Manas Sajjan, Sabre Kais

Published 2026-11-02

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Shah Ishmam Mohtashim, Manas Sajjan, Sabre Kais

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

In the microscopic world of molecules, shape is destiny. Some molecules are chiral, meaning they exist in two forms that are mirror images of each other, much like a left hand and a right hand. These two versions cannot be superimposed; they are distinct entities. For many years, scientists have been fascinated by how these mirror-image forms can switch places. In certain simple molecules, the atoms are not locked in a rigid structure but can twist and turn, allowing the molecule to flip from its left-handed form to its right-handed form and back again. This flipping happens through a process called quantum tunneling, a phenomenon where a particle passes through a barrier it seemingly cannot climb over, simply because the laws of quantum mechanics allow it to exist on the other side. The speed of this flip determines how long the molecule stays in one shape before becoming its mirror image. Understanding this speed is crucial for fields ranging from chemistry to the study of life's origins, as it tells us how stable these shapes are and how they might interact with light.

A team of researchers set out to investigate a specific question about this flipping process in three related molecules: hydrogen peroxide, hydrogen disulfide, and hydrogen diselenide. These molecules share a similar structure, with two hydrogen atoms attached to a chain of two other atoms (oxygen, sulfur, or selenium). The researchers wanted to know if slowing down the tunneling process—making the flip take longer—would also weaken the molecule's ability to interact with light in a specific way. In physics, this interaction is often measured by something called polarizability, which describes how easily a molecule's electron cloud can be distorted by an electric field. The prevailing intuition might suggest that if the molecule is sluggish and takes a long time to flip, its internal response to light should also be weak. However, the researchers found that this intuition is incorrect. They discovered that even as the tunneling time slowed down dramatically, the molecule's ability to respond to light remained strong, and in some cases, even grew stronger.

To uncover this relationship, the scientists built a detailed computer model of these three molecules. They started by calculating the exact energy landscape the molecules travel through as they twist. This landscape looks like a valley with two dips, separated by a hill. One dip represents the left-handed shape, and the other represents the right-handed shape. The molecule sits in one dip, but thanks to quantum mechanics, it can tunnel through the hill to the other dip. The researchers used advanced computational methods to map out this terrain with high precision, accounting for the masses of the atoms and the forces holding them together. They then calculated the specific quantum states the molecule can occupy within this valley and determined how fast it would tunnel from one side to the other. Crucially, they also calculated the molecule's polarizability, specifically looking at how it responds to light as it twists, to see if this property changed in step with the tunneling speed.

The results revealed a striking disconnect between the speed of the flip and the strength of the light interaction. As the researchers moved from hydrogen peroxide to hydrogen disulfide and then to hydrogen diselenide, the tunneling time slowed down by a factor of more than six million. The molecule went from flipping in a few trillionths of a second to taking several microseconds, which is a vast difference in the molecular world. Despite this massive slowdown, the strength of the polarizability interaction did not fade away. Instead, it increased by nearly sixty percent. The researchers found that the tunneling speed depends on how easily the molecule can penetrate the barrier between the two shapes, a process that becomes much harder as the atoms get heavier. In contrast, the polarizability is determined by what happens deep inside the valleys where the molecule spends most of its time. Because the shape of these valleys and the behavior of the electrons within them remained favorable, the molecule retained its strong ability to interact with light, regardless of how slowly it was flipping.

To confirm that this was truly a matter of mass and not a change in the molecule's electronic structure, the team performed a controlled test using hydrogen peroxide. They replaced the light hydrogen atoms with heavier deuterium atoms, a heavier version of hydrogen, while keeping the rest of the molecule exactly the same. This change slowed the tunneling time by nearly seven times. Yet, the polarizability interaction changed by less than two percent. This experiment proved that the slowing of the flip and the strength of the light response are governed by different parts of the molecule's physics. The tunneling time is sensitive to the barrier crossing, while the polarizability is a property of the stable shapes themselves. The study also looked at higher energy states where the molecules were more excited, finding that even when the tunneling time varied by orders of magnitude, the polarizability remained remarkably stable.

These findings suggest that the ability of chiral molecules to interact with light is not tied to how fast they flip between shapes. This is a significant insight for scientists who hope to use light to control or distinguish between these mirror-image forms. If the interaction strength faded as the molecules slowed down, it would be much harder to manipulate them with lasers or other optical tools. The fact that the interaction remains strong even for slow-flipping molecules means that researchers have a robust physical mechanism to work with. The study confirms that while the clock of molecular flipping can be slowed down by changing the atoms, the molecule's intrinsic response to light remains a constant, powerful feature. This separation of timescales and strengths provides a clearer picture of how chiral molecules behave, offering a solid foundation for future experiments aimed at controlling molecular handedness with light.

Technical Summary: Chiral Torsional Double Wells: Tunneling and Polarizability

Problem Statement
Chiral molecules such as hydrogen peroxide ( $H_2O_2$ ), hydrogen disulfide ( $H_2S_2$ ), and hydrogen diselenide ( $H_2Se_2$ ) possess enantiomeric conformers that interconvert via internal rotation (torsion) around the X–X bond. This process occurs through quantum tunneling between mirror-related potential wells. The timescale for this coherent interconversion is determined by the energy splitting of the lowest even and odd torsional eigenstates. While it is established that tunneling rates vary dramatically across this chemical series (spanning nearly seven orders of magnitude), it remains unclear whether a reduction in tunneling speed (slower interconversion) correlates with a weakening of the intrinsic transition polarizability. Understanding this relationship is critical for assessing the feasibility of controlling chiral torsional motion via light-matter interactions, specifically Raman processes, which rely on polarizability rather than dipole moments.

Methodology
The authors employed a rigorous computational framework to decouple nuclear dynamics from electronic response:

Potential Energy Surfaces: Relaxed one-dimensional torsional paths were constructed using the CCSD(T)//MP2 method. Geometries were optimized at the MP2 level under $C_2$ symmetry for each dihedral angle ( $\phi$ ), while single-point energies were calculated using CCSD(T). The path spans from cis ( $\phi=0$ ) to trans ( $\phi=\pm\pi$ ).
Kinetic Model: To isolate torsional motion, overall molecular translation and rotation were removed. An effective moment of inertia, $I_{eff}(\phi)$ , was calculated by projecting rigid rotations from the internal displacement vectors. A coordinate transformation to a signed mass-weighted arc length ( $\chi$ ) was applied to yield a Hamiltonian with a constant kinetic coefficient ( $I_{arc}$ ).
Quantum Dynamics: The torsional Hamiltonian was solved using a periodic Fourier discrete-variable representation (DVR) with $N=384$ grid points. This allowed for the precise calculation of eigenstates, energy splittings ( $\Delta_n$ ), and tunneling times ( $t_{L\to R}$ ).
Electronic Properties: Static polarizabilities were computed using the PBE0 functional at the MP2 geometries. The electric dipole and polarizability tensors were transformed into a body-fixed molecular frame aligned with the torsional path.
Symmetry Analysis: The study explicitly analyzed selection rules based on $C_2$ rotation and reflection symmetry ( $\sigma_{xz}$ ) to determine which molecular properties can couple opposite-parity torsional states.
Comparative Cases: The study examined three chemical substitutions ( $O \to S \to Se$ ) and one isotopic substitution ( $H \to D$ in $H_2O_2$ ) to distinguish between changes in the electronic surface and changes in nuclear mass.

Key Contributions and Results

Symmetry Selection Rules: The authors demonstrated that in the defined molecular frame, the electric dipole component ( $\mu_x$ ) is even under well exchange, while the off-diagonal polarizability component ( $\alpha_{yz}$ ) is odd. Consequently, the transition matrix element for the dipole between the lowest even and odd states vanishes ( $\langle + | \hat{\mu}_x | - \rangle = 0$ ). However, the matrix element for the polarizability, $A_{yz} = \langle + | \hat{\alpha}_{yz} | - \rangle$ , is symmetry-allowed and non-zero.
Divergence of Tunneling and Coupling: Across the chemical series ( $H_2O_2 \to H_2S_2 \to H_2Se_2$ ), the ground-state tunneling splitting decreases by a factor of $6.3 \times 10^6$ , increasing the coherent transfer time from 1.5 ps to 9.57 $\mu$ s. Contrary to the hypothesis that slower tunneling implies weaker coupling, the magnitude of the transition polarizability matrix element $|A_{yz}^{01}|$ increases by 59% (from 1.57 to 2.49 a.u.).
Mechanism of Decoupling: The study reveals that tunneling rates are governed by wavefunction penetration through the potential barriers, whereas the polarizability matrix element is formed primarily within the potential wells. Because the electronic response (polarizability) is concentrated in the wells, it remains robust even as the barrier penetration (and thus tunneling rate) diminishes significantly.
Excited States and Isotopes:
- Excited Doublets: In $H_2S_2$ and $H_2Se_2$ , higher torsional doublets exhibit tunneling times that are orders of magnitude faster than the ground state, yet their squared coupling ( $|A_{yz}|^2$ ) retains over 90% of the ground-state value.
- Deuteration: Replacing hydrogen with deuterium in $H_2O_2$ slows the tunneling time by a factor of 6.75 but alters the transition polarizability matrix element by only 1.01%. This confirms that the decoupling of tunneling speed and polarizability strength is driven by nuclear mass effects on the kinetic term, not changes in the electronic surface.
Benchmarking: Calculated torsional energies for $H_2O_2$ and $D_2O_2$ agree with experimental assignments within a few percent, validating the relaxed path model despite its one-dimensional approximation.

Significance
The paper establishes that the slowing of coherent interconversion between chiral conformers does not necessitate a weakening of the intrinsic transition polarizability. This finding challenges the assumption that slower tunneling renders chiral molecules less responsive to optical control. Specifically, it suggests that the odd polarizability component $\alpha_{yz}$ provides a viable, symmetry-allowed pathway for Raman-based manipulation of chiral torsional motion, even in heavy analogues where tunneling is extremely slow. The work provides a physical basis for exploring chiroptical control mechanisms that rely on polarizability rather than dipole transitions, distinguishing the spatial origins of tunneling (barrier penetration) and polarizability coupling (well localization).

More like this