Weiner's theory for exactly solvable Schrödinger equation with symmetric double well potential

This paper develops a modified Weiner's theory based on an exactly solvable Schrödinger equation with a symmetric trigonometric double-well potential to accurately calculate proton transfer rates, demonstrating its application to the ammonia dimer cation by capturing the transition from thermal activation to quantum tunneling and vibrationally enhanced tunneling.

The Gist

The Big Picture: A Proton's Mountain Climb

Imagine a proton (a tiny, positively charged particle) trapped in a valley. But this isn't just any valley; it's a double-well potential. Think of it as a "W" shape. The proton is sitting in the left dip of the W, but it wants to get to the right dip.

To get there, it has to climb over the middle hump (the barrier). In the world of chemistry, this is called Proton Transfer, and it's the engine behind many reactions, including how enzymes in our bodies work.

The author, A.E. Sitnitsky, is trying to solve a very old problem: How fast does the proton get across?

The Old Way: The "Rough Sketch" (Weiner's Theory)

For decades, scientists used a method called Weiner's Theory (WT) to guess the speed.

• The Problem: The old method treated the mountain (the potential barrier) like a jagged, blocky structure made of flat pieces of wood glued together.
• The Flaw: To make the math work with these blocky mountains, the scientists had to make some very "rough" guesses (approximations). It was like trying to calculate the speed of a car driving over a bumpy road by pretending the road is perfectly flat, then just adding a little "bumpiness" factor at the end.
• The Result: It gave decent answers, but it was messy, relied on shaky assumptions, and sometimes broke down when the mountain wasn't perfectly symmetrical.

The New Way: The "Smooth Slide" (Modified Weiner's Theory)

The author proposes a Modified Weiner's Theory (mWT).

• The Upgrade: Instead of a blocky, jagged mountain, the author uses a Trigonometric Double-Well Potential (TDWP). Imagine the mountain is now a perfectly smooth, curved slide made of silk.
• The Magic: Because this "smooth slide" is mathematically "exactly solvable" (meaning we have a perfect map of it), the author doesn't need to make those rough guesses anymore. They can calculate the proton's journey with high precision.
• The Tool: The author uses a powerful computer software (Mathematica) that knows the "language" of these smooth curves (called Spheroidal Functions) to do the heavy lifting.

The Two Ways to Cross the Mountain

The paper explains that the proton can cross the barrier in two ways, depending on the temperature:

1. Thermal Activation (The Hiker): At high temperatures, the proton is energetic. It gets enough "heat energy" to climb over the top of the mountain. This follows the standard rules of physics (Arrhenius behavior).
2. Quantum Tunneling (The Ghost): At low temperatures, the proton is cold and lazy. It doesn't have the energy to climb over. However, because it's a quantum particle, it can act like a ghost and tunnel through the mountain instead of going over it.

The "Crossover": The paper shows exactly where the proton switches from being a "Hiker" (climbing over) to a "Ghost" (tunneling through). For the specific molecule they studied (the Ammonia Dimer Cation), this switch happens around 60 Kelvin (very cold, about -213°C).

The "Super-Charge" Effect (Vibrationally Enhanced Tunneling)

This is the most exciting part of the paper. The author looks at a phenomenon called Vibrationally Enhanced Tunneling (VET).

• The Analogy: Imagine the proton is trying to tunnel through a thick wall. Usually, it's hard. But, imagine the wall is attached to a giant speaker. If you play a specific musical note (a specific vibration frequency) that matches the wall's natural rhythm, the wall starts to shake violently.
• The Result: When the wall shakes at just the right frequency, the proton doesn't just tunnel; it gets a massive "boost."
• The Math: The author shows that at this "resonant" frequency, the speed of the reaction can jump by 26 orders of magnitude. That's like going from walking at a snail's pace to moving faster than the speed of light (in a relative sense).
• Why it matters: This explains how enzymes work so incredibly fast. Enzymes might be "vibrating" the chemical bonds in just the right way to help protons tunnel through barriers instantly, speeding up life-saving reactions by factors of a trillion or more.

The Real-World Test: The Ammonia Dimer

To prove their new theory works, the author tested it on a specific molecule: the Proton-Bound Ammonia Dimer Cation ($N_2H_7^+$).

• They took real data from experiments (infrared spectroscopy) and computer simulations (quantum chemistry) to build their "smooth mountain."
• They ran the numbers.
• The Verdict: The new theory (mWT) successfully predicted the transition from "climbing over" to "tunneling through" and matched the behavior of the molecule perfectly. It also showed that the old theory (WT) would have struggled with this specific molecule because the mountain wasn't perfectly symmetrical.

Summary

In simple terms, this paper says:

"We found a better way to calculate how fast protons move between molecules. Instead of using a jagged, approximate map (the old way), we used a perfectly smooth, mathematically precise map (the new way). This allows us to see exactly when protons switch from climbing over barriers to tunneling through them, and it explains how vibrations can supercharge these reactions, helping us understand how enzymes work so efficiently."

It's a move from guessing the speed of a particle to calculating it with the precision of a master cartographer.

Technical Summary

1. Problem Statement

The paper addresses the theoretical challenges in modeling proton transfer (PT) reactions in hydrogen bonds (HB), particularly within the context of enzymatic hydrogen transfer (EHT) and vibrationally enhanced tunneling (VET).

• Limitations of Existing Models: Traditional reaction rate theories (e.g., Transition State Theory, Kramers' theory) often struggle with quantum effects like tunneling. The original Weiner's Theory (WT) is a phenomenological approach for barrier crossing but relies on severe approximations:
  • It uses piecewise quadratic double-well potentials (DWP), which introduce artificial "sewing" errors at potential boundaries.
  • It requires quasi-classical (WKB) approximations for wave functions and transmission coefficients, leading to rough estimates (e.g., assuming $\psi_n \approx \psi_{n+1}$).
  • It assumes non-normalizable right- and left-moving states, creating conceptual difficulties in defining boundary conditions and flux.
  • It struggles with asymmetric or split doublets where one level of a doublet is below the barrier and the other is above.
• The Need: There is a need for a mathematically rigorous, exactly solvable framework that avoids these approximations to accurately describe the transition from thermal activation (Arrhenius behavior) to quantum tunneling, especially for low-barrier hydrogen bonds.

2. Methodology: Modified Weiner's Theory (mWT)

The author proposes a Modified Weiner's Theory (mWT) based on an exactly solvable Schrödinger equation (SE) using a Trigonometric Double-Well Potential (TDWP).

• The Potential: The study utilizes a symmetric TDWP defined as:
  $$U(x) = \left(\frac{m^2 - 1}{4}\right)\tan^2 x - p^2 \sin^2 x$$
  where $m$ is an integer and $p$ is a real number. This potential is smooth (not piecewise) and defined on a finite interval $[-\pi/2, \pi/2]$.
• Exact Solutions: Unlike the original WT, the SE for TDWP has exact analytic solutions in terms of normalized angular prolate spheroidal functions (SF), which are implemented in standard software (Mathematica).
  • Wave Functions: $\psi_q(x) = \sqrt{\cos x} \, \bar{S}^{m}_{(q+m)}(p; \sin x)$
  • Energy Levels: Derived from the eigenvalue spectrum of the spheroidal functions.
• Derivation of Flux and Transmission:
  • Instead of using WKB approximations, the author decomposes the known exact wave function into right-moving ($\psi_{q;r}$) and left-moving ($\psi_{q;l}$) components using a phase function $\chi_q(x)$.
  • A new relationship is derived for the probability flux $J_{q;r}$ that avoids divergent integrals by introducing a regularization constant $\mu_q^2$.
  • The transmission coefficient ($|T_q|^2$) is calculated by matching scattering states at the barrier turning points, utilizing the exact wave functions rather than approximations.
• Rate Constant Calculation: The PT rate constant $k(\beta)$ is obtained via Boltzmann averaging of the product of the flux and transmission coefficient for all energy levels (both below and above the barrier).

3. Key Contributions

1. Elimination of Severe Approximations: By using the exact analytic solution of the TDWP, mWT removes the need for piecewise potentials and the rough WKB approximations inherent in the original WT.
2. Handling Split Doublets: The theory successfully handles cases where a doublet is "split" by the barrier top (one level below, one above), a scenario where the original WT fails or requires ad-hoc fixes.
3. Analytic Formula for Rate Constants: The author derives a closed-form analytic expression (Eq. 34) for the rate constant that can be directly implemented in computational software (Mathematica), facilitating parameter scanning.
4. Vibrationally Enhanced Tunneling (VET) Framework: The formalism is shown to be compatible with VET. By coupling the proton coordinate to an external oscillator (simulating enzyme dynamics or polariton fields), the model predicts resonant increases in flux.

4. Results and Applications

The theory was applied to two specific systems:

• Proton-Bound Ammonia Dimer Cation ($N_2H_7^+$):

  • Parameters: Extracted from IR spectroscopy and quantum chemical calculations ($R_{NN} = 3.15$ Å).
  • Findings: The model successfully reproduces the transition from Arrhenius-like exponential temperature dependence (thermal activation) to a temperature-independent plateau (quantum tunneling).
  • Crossover Temperature ($T_c$): The calculated crossover temperature is $T_c \approx 60$ K, which aligns well with the Goldanskii criterion estimate ($\approx 50$ K).
  • Comparison: mWT results for $N_2H_7^+$ show good agreement with WT where WT is applicable but provide a more rigorous foundation for split doublets.

• Zundel Ion ($H_5O_2^+$) and VET:

  • The model was extended to a 2D scenario (reduced to effective 1D) to study VET.
  • Resonance Effect: At a specific resonant frequency of an external oscillator, the probability flux $J_{q;r}$ increases dramatically.
  • Magnitude: The model predicts a flux increase of 26 orders of magnitude at resonance compared to the non-resonant plateau. This magnitude is sufficient to explain the massive rate accelerations observed in enzymatic reactions (e.g., factors of $10^{12}$ to $10^{26}$).

5. Significance and Conclusion

• Theoretical Robustness: mWT provides a self-consistent, mathematically rigorous alternative to the original WT, grounded in exact solutions rather than phenomenological approximations.
• Chemical Relevance: It offers a viable tool for interpreting experimental data in enzymatic hydrogen transfer and polariton chemistry, specifically explaining how vibrational modes can catalyze proton transfer via resonance.
• Limitations: While mWT alleviates many technical flaws of WT, the author notes that the fundamental conceptual issue of non-normalizable scattering states in the decomposition of the wave function remains, though its consequences are mitigated.
• Future Work: The approach is applicable to various hydrogen-bonded systems (intramolecular and intermolecular) provided that reliable IR or quantum chemical data exist to parameterize the potential.

In summary, Sitnitsky presents a refined theoretical framework that bridges the gap between exact quantum mechanics and reaction rate theory, offering a powerful tool to investigate the quantum nature of proton transfer in complex biological and chemical environments.