Tunneling molecular dynamics in the light of the corpuscular-wave dualism theory

This paper presents the experimental demonstration of the corpuscular-wave dualism theory. The correlation between the de Broglie wavelength related to the thermal motion and the potential barrier width and height is reported. The stochastic jumps of light atoms (hydrogen, deuterium) between two equilibrium sites A and B (identical geometry) occur via different pathways; one pathway is over the barrier (classical dynamics), and the other one is through the barrier (tunneling). On the over-the-barrier pathway, there are no obstacles for the de Broglie waves, and this pathway exists from high to low temperatures up to 0 K because the thermal energy is subjected to the Maxwell distribution and a certain number of particles owns enough energy for the hopping over the barrier. On the tunneling pathway, the particles pass through the barrier, or they are reflected from the barrier. Only particles with the energy lower than barrier heights are able to perform a tunneling hopping. The de Broglie waves related to these energies are longer than the barrier width. The Schr?dinger equation is applied to calculate the rate constant of tunneling dynamics. The Maxwell distribution of the thermal energy has been taken into account to calculate the tunneling rate constant. The equations for the total spectral density of complex motion derived earlier by us together with the expression for the tunneling rate constant, derived in the present paper, are used in analysis of the temperature dependence of deuteron spin-lattice relaxation of the ammonium ion in the deuterated analogue of ammonium hexachloroplumbate ((ND4)2PbCl6). It has been established that the equation CpTtun = EH (thermal energy equals activation energy), where Cp is the molar heat capacity (temperature-dependent, known from literature), determines directly the low temperature Ttun at which the de Broglie wavelength, lambdadeBroglie, related to the thermal energy, CpT, is equal to the potential barrier width, L. Above Ttun, the lambdadeBroglie wavelength related to the CpT energy is shorter than the potential barrier width and not able to overcome the barrier. The activation energy EH equals 7.5 kJ/mol, and therefore, the Ttun temperature for deuterons in ((ND4)2PbCl6 is 55.7 K. The agreement between the potential barrier width following from the simple geometrical calculations (L = 0.722 A) and de Broglie wavelength at Ttun (L = 0.752 A) is good. The temperature plots of the deuteron correlation times for (ND4)2PbCl6 reveal comparable values of the correlation times of the tunneling, (tau(T)), and over-the-barrier jumps (tau(H)) near 34.8 K. Matsuo, on the basis of the molar heat capacity study, found the first-order phase transition at this temperature.     

The ground state tunneling splitting and the zero point energy of malonaldehyde: a quantum Monte Carlo determination

Quantum dynamics calculations of the ground state tunneling splitting and of the zero point energy of malonaldehyde on the full dimensional potential energy surface proposed by Yagi et al. [J. Chem. Phys. 1154, 10647 (2001)] are reported. The exact diffusion Monte Carlo and the projection operator imaginary time spectral evolution methods are used to compute accurate benchmark results for this 21-dimensional ab initio potential energy surface. A tunneling splitting of 25.7+/-0.3 cm-1 is obtained, and the vibrational ground state energy is found to be 15 122+/-4 cm-1. Isotopic substitution of the tunneling hydrogen modifies the tunneling splitting down to 3.21+/-0.09 cm-1 and the vibrational ground state energy to 14 385+/-2 cm-1. The computed tunneling splittings are slightly higher than the experimental values as expected from the potential energy surface which slightly underestimates the barrier height, and they are slightly lower than the results from the instanton theory obtained using the same potential energy surface.     

Tunneling time in driven double-well system using entangled molecular dynamics method

We calculate quantum tunneling time in the double-well system without perturbation and with symmetry-breaking driven force using the entangled trajectory molecular dynamics method in the present article. Without perturbation, quantum tunneling time decreased with the increase of the energy, and the different contributions of the barrier traversal time and the intrinsic decay time have been shown. The tunneling time dependence on the amplitude and frequency of the symmetry-breaking driven force are present. In the case of weak driven force, tunneling time has a minimum value in the resonant frequency. For strong driven force, chaos brings a huge change to quantum dynamics, tunneling time significantly becomes short. Finally, we directly show the enhancement of quantum tunneling process by chaotic behavior of entangled trajectories and indicate that it was caused by quantum effect.     

Applicability of the WKB method in asymmetric double wells with degenerate and nondegenerate minima

To test the applicability of the WKB semiclassical approximation to correctly describe the tunneling frequencies and energy levels in asymmetric double wells, we have considered the cases where the asymmetry comes from the nondegeneracy of the two minima and when it comes from the asymmetry in the shape of the barrier between both minima. To do this, we compare the tunneling frequencies and energy levels obtained through the WKB method for symmetric and asymmetric cases with exact results obtained by a basis set procedure. Our results show that the semiclassical WKB approximation has to be used with caution for asymmetric double wells. © 1994 by John Wiley & Sons, Inc.     

An ab initio close-coupling calculation of the lower vibrational energies of the HCl dimer

We have previously determined an analytical ab initio six-dimensional potential energy surface for the HCl dimer, and in the present paper we use this potential, with the HCl bond lengths held fixed, in a full (four-dimensional) close-coupling calculation to determine the energies of the lowest 24 vibrational states. These vibrational states involve the intermolecular stretch ν₄, the trans-bend tunneling vibration ν₅, and the torsion ν₆. The highest of the 24 levels is the (ν₄ν₅ν₆)=(111) state, for which we calculate an energy of 200 cm⁻¹ above the (000) state. As well as determining tunneling energies up to 5ν₅=183 cm⁻¹, we determine ν₄=49 cm⁻¹, 2ν₄=93 cm⁻¹, 3ν₄=134 cm⁻¹, 4ν₄=172 cm⁻¹, ν₆=137 cm⁻¹ and ν₄+ν₆=178 cm⁻¹, together with tunneling energies in all these states. Making allowance for the HCl stretching zero-point energy we determine the dissociation energy D₀ as 390 cm⁻¹ on this analytical surface. We determine that below 300 cm⁻¹ there are 72 vibrational (J=K=0) states, and below dissociation there are 162 vibrational (J=K=0) states, for this potential surface.     

Energy level and band alignment for GaAs-alkylthiol monolayer-Hg junctions from electrical transport and photoemission experiments

A series of p- and n-GaAs-S-C(n)H(2n+1) || Hg junctions are prepared, and the electronic transport through them is measured. From current-voltage measurements, we find that, for n-GaAs, transport occurs by both thermionic emission and tunneling, with the former dominating at low forward bias and the latter dominating at higher forward bias. For p-GaAs, tunneling dominates at all bias voltages. By combining the analysis of the transport data with results from direct and inverse photoemission spectroscopy, we deduce an energy band diagram of the system, including the tunnel barrier and, with this barrier and within the Simmons tunneling model, extract an effective mass value of 1.5-1.6m(e) for the electronic carriers that cross the junctions. We find that transport is well-described by lowest unoccupied and highest occupied states at 1.3-1.4 eV above and 2.0-2.2 eV below the Fermi level. At the same time, the photoemission data indicate that there are continua of states from the conduction band minimum and the valence band maximum, the density of which varies with energy. On the basis of our results, it appears likely that, for both types of junctions, electrons are the main carrier type, although holes may contribute significantly to the transport in the p-GaAs system.     

化学动力学中隧道效应的研究（Ⅱ）──甲醛单分子反应渠道间的相互作用研究

对甲醛单分子热反应模拟给出了分子渠道Ｅｃｋａｒｔ垒，计算了该渠道在不同温度下的隧道校正因子，研究了分子渠道与游离基渠道间的相互作用，计算了存在反应渠道间相互作用时各渠道的Ｋ＿（ｕｎｉ）（ｉ）、Ｅａ（ｉ），分子渠道考虑隧道效应，而游离基渠道无隧道效应，计算结果与实验结果一致。     

The importance of tunneling in the first hydrogenation step in ammonia synthesis over a Ru(0001) surface

The hydrogenation of nitrogen (N(ads)+H(ads)-->NH(ads)) on metal surfaces is an important step in ammonia catalysis. We investigate the reaction dynamics of this hydrogenation step by time independent scattering theory and variational transition state theory (VTST) including tunneling corrections. The potential energy surface is derived by hybrid density functional theory on a model cluster composed of 12 ruthenium atoms resembling a Ru(0001) surface. The scattering calculations are performed on a reduced dimensionality potential energy hypersurface, where two dimensions are treated explicitly and all others are included implicitly by the zero-point correction. The VTST calculations include quantum effects along the reaction coordinate by applying the small curvature tunneling scheme. Even at room temperature (where ruthenium already shows catalytic activity) we find rate enhancement by tunneling by a factor of approximately 70. Inspection of the reaction probabilities shows that the major contribution to reactivity comes from the vibrational ground state of the reactants into vibrationally excited product states. The reaction rates are higher than determined in previous studies, and are compatible with experimental overall rates for ammonia synthesis.     

Including anharmonicity in the calculation of rate constants. 1. The HCN/HNC isomerization reaction

A method for calculating anharmonic vibrational energy levels in asymmetric top and linear systems that is based on second-order perturbation theory in curvilinear coordinates is extended to the bound generalized normal modes at nonstationary points along a reaction path. Explicit formulas for the anharmonicity coefficients, x(ij), and the constant term, E0, are presented, and the necessary modifications for resonance cases are considered. The method is combined with variational transition state theory with semiclassical multidimensional tunneling approximations to calculate thermal rate constants for the HCN/HNC isomerization reaction. Although the results for this system are not very sensitive to the choice of coordinates, we find that the inclusion of anharmonicity leads to a substantial improvement in the vibrational energy levels. We also present detailed comparisons of rate constants computed with and without anharmonicity, with various approximations for incorporating tunneling along the reaction path, and with a more practical approach to calculating the vibrational partition functions needed for larger systems.     

Very fast tunneling in the early stage of reaction dynamics

The difference between quantum and classical survival probabilities for molecular dissociation dynamics in the time domain, which arises mainly from quantum mechanical tunneling, has interesting characteristics that are not noticed through the counterpart in energy domain. It is shown that the early stage undergoes a fast tunneling, while the later stage is characterized with a long-lasting slow tunneling. The mechanism of this behavior is analyzed in terms of a quasi-semiclassical theory featuring the geometrical distribution of the so-called tunneling points. In particular, the role of dynamical tunneling is discussed as a phenomenon that typifies the time dependence of tunneling dynamics. It is predicted that these tunneling characteristics will be reflected in the isotope effect and should be experimentally observable.     

Tunneling and classical paths for proton transfer in an enzyme reaction dominated by tunneling: oxidation of tryptamine by aromatic amine dehydrogenase

Proton tunneling dominates the oxidative deamination of tryptamine catalyzed by the enzyme aromatic amine dehydrogenase. For reaction with the fast substrate tryptamine, a H/D kinetic isotope effect (KIE) of 55 +/- 6 has been reported-one of the largest observed in an enzyme reaction. We present here a computational analysis of this proton-transfer reaction, applying combined quantum mechanics/molecular mechanics (QM/MM) methods (PM3-SRP//PM3/CHARMM22). In particular, we extend our previous computational study (Masgrau et al. Science 2006, 312, 237) by using improved energy corrections, high-level QM/MM methods, and an ensemble of paths to estimate the tunneling contributions. We have carried out QM/MM molecular dynamics simulations and variational transition state theory calculations with small-curvature tunneling corrections. The results provide detailed insight into the processes involved in the reaction. Transfer to the O2 oxygen of the catalytic base, Asp128beta, is found to be the favored reaction both thermodynamically and kinetically, even though O1 is closer in the reactant complex. Comparison of quantum and classical models of proton transfer allows estimation of the contribution of hydrogen tunneling in lowering the barrier to reaction in the enzyme. A reduction of the activation free energy due to tunneling of 3.1 kcal mol-1 is found, which represents a rate enhancement due to tunneling by 2 orders of magnitude. The calculated KIE of 30 is significantly elevated over the semiclassical limit, in agreement with the experimental observations; a semiclassical value of 6 is obtained when tunneling is omitted. A polarization of the C-H bond to be broken is observed due to the close proximity of the catalytic aspartate and the (formally) positively charged imine nitrogen. A comparison is also made with the related quinoprotein methylamine dehydrogenase (MADH)-the much lower KIE of 11 that we obtain for the MADH/methylamine system is found to arise from a more endothermic potential energy surface for the MADH reaction.     

Interface and molecular electronic structure vs tunneling characteristics of CH3- and CF3-terminated thiol monolayers on Au(111)

By means of density functional theory calculations, we investigate work functions, energy level alignments, charge transfers, and tunneling characteristics of CH3- and CF3-terminated alkane- and diphenylthiol monolayers on Au(111). While the alignments of the energy levels and the charge transfers at the metal-molecule interface are found to be determined by the value of the clean Au surface work function relative to the HOMO ionization potential (IP) at the thiolate end of the monolayer, the change of work function for the modified Au(111) surface is dominated by the properties of the thiolate monolayer, including the character, saturated or conjugated, of the molecule and the chemical nature and orientation of the terminal group. The tunneling currents through the adsorbed molecular monolayers are calculated using the Tersoff-Hamann approach. The computed difference between the I-V characteristics for the CH3- and CF3-terminated alkanethiol monolayers agree well with available experimental data. The energy barrier at the metal-molecule interface, the molecular electronic structure, and the IP of the terminal group are the key parameters which determine the tunneling properties.     

Correlated ab initio study of the ground electronic state of the O2-HF complex

In this paper, we present the first correlated ab initio investigations on the ground electronic state of the O(2)-HF complex. Calculations were performed using the CCSD(T) method with the aug-cc-pVDZ and aug-cc-pVTZ basis sets. The results show that there are two equivalent minimum energy hydrogen-bonded structures of planar bent geometry, where the minima correspond to exchange of the oxygen atoms. For each minimum the length of the O-H hydrogen bond is 2.16 A. The best calculated value of D(e) of the equivalent minima is 271 cm(-1). The T-shaped geometry of the complex, with oxygen perpendicular to the axis connecting the center of masses of O(2) and the HF molecule, represents a barrier to tunneling between the equivalent minima. The best estimated value of that barrier height is 217 cm(-1). The linear O-O-HF geometry of the complex represents a saddle point. The calculated geometrical parameters of the minimum energy structure of the complex are in reasonable agreement with the previously reported spectroscopic results. However, results of the current calculations suggest that a full understanding of the fine structures of the observed infrared spectrum of the complex requires the development of an effective Hamiltonian that takes the effects of tunneling into account.     

Site-directed electronic tunneling in a dissipative molecular environment

The ability to control electronic tunneling in complex molecular networks of multiple donor/acceptor sites is studied theoretically. Our past analysis, demonstrating the phenomenon of site-directed transport, was limited to the coherent tunneling regime. In this work we consider electronic coupling to a dissipative molecular environment including the effect of decoherence. The nuclear modes are classified into two categories. The first kind corresponds to the internal molecular modes, which are coupled to the electronic propagation along the molecular bridges. The second kind corresponds to the external solvent modes, which are coupled to the electronic transport between different segments of the molecular network. The electronic dynamics is simulated within the effective single electron picture in the framework of the tight binding approximation. The nuclear degrees of freedom are represented as harmonic modes and the electronic-nuclear coupling is treated within the time-dependent Redfield approximation. Our results demonstrate that site-directed tunneling prevails in the presence of dissipation, provided that the decoherence time is longer than the time period for tunneling oscillations (e.g., at low temperatures). Moreover, it is demonstrated that the strength of electronic coupling to the external nuclear modes (the solvent reorganization energy) controls the coherent intramolecular tunneling dynamics at short times and may be utilized for the experimental control of site-directed tunneling in a complex network.     

Stereomutation tunneling switching dynamics and parity violation in chlorineperoxide Cl-O-O-Cl

In a search for efficient spectroscopic avenues toward experiments on molecular parity violation, we investigate the stereomutation tunneling processes in the axially chiral chlorine isotopomers of Cl2O2 by the quasi-adiabatic channel reaction path Hamiltonian (RPH) approach and the corresponding parity violating potentials by means of quantum chemical calculations including our recently developed Multiconfiguration linear response (MC-LR) approach to electroweak quantum chemistry. The calculated ground-state torsional tunneling splittings for all isotopomers of Cl2O2 are much smaller than the parity violating energy differences Delta(pv)E between the enantiomers of these molecules and therefore parity violation is predicted to dominate the quantum dynamics of stereomutation at low energies. We also compare these with torsional ground-state tunneling splittings and parity violating energy differences of the whole series of axially chiral HXYH(+) isotopomers (with X, Y= Cl(+), O, S, Se, Te). A comparison with our previous results for the homologous molecule Cl2S2 shows that for Cl2O2 a spectroscopic high-resolution analysis should be easier and the energy region of large tunneling splittings should be more easily accessible by IR excitation. We thus propose a scheme using "tunneling switching" with vibrational excitation in order to carry out the measurement of time-dependent parity violation in superposition states of initially well-defined parity. We discuss the advantages and drawbacks of such an experiment that can be carried out entirely in the IR spectral range (for Cl2O2 or related molecules).     

Large tunneling effect on the hydrogen transfer in bis(μ-oxo)dicopper enzyme: a theoretical study

Type-III copper-containing enzymes have dicopper centers in their active sites and exhibit a novel capacity for activating aliphatic C-H bonds in various substrates by taking molecular oxygen. Dicopper enzyme models developed by Tolman and co-workers reveal exceptionally large kinetic isotope effects (KIEs) for the hydrogen transfer process, indicating a significant tunneling effect. In this work, we demonstrate that variational transition state theory allows accurate prediction of the KIEs and Arrhenius parameters for such model systems. This includes multidimensional tunneling based on state-of-the-art quantum-mechanical calculations of the minimum-energy path (MEP). The computational model of bis(μ-oxo)dicopper enzyme consists of 70 atoms, resulting in a 204-dimensional potential energy surface. The calculated values of E(a)(H) - E(a)(D), A(H)/A(D), and the KIE at 233 K are -1.86 kcal/mol, 0.51, and 28.1, respectively, for the isopropyl ligand system. These values agree very well with experimental values within the limits of experimental error. For the representative tunneling path (RTP) at 233 K, the pre- and post-tunneling configurations are 3.3 kcal/mol below the adiabatic energy maximum, where the hydrogen travels 0.54 ? by tunneling. We found that tunneling is very efficient for hydrogen transfer and that the RTP is very different from the MEP. It is mainly heavy atoms that move as the reaction proceeds from the reactant complex to the pretunneling configuration, and the hydrogen atom suddenly hops at that point.     

Effect of Coulomb interaction between the electrons on two-electron redox-mediated tunneling

Two-electron non-adiabatic redox-mediated tunneling through a symmetric electrochemical contact with a bridge molecule having one electron energy level participating in tunneling is considered under ambient conditions. It is shown that the current/overpotential dependence in this system can disclose two distinct or overlapping clear-cut maxima depending on the value of the effective Coulomb repulsion energy. This new effect is due to the opening of the channel for tunneling of second electron with the variation of the electrode potential. The system manifests also a rectification effect in the current/bias voltage curve which depends on the value of the effective Coulomb repulsion energy.     

Trajectory on-the-fly molecular dynamics approach to tunneling splitting in the electronic excited state: A case of tropolone

The semiclassical tunneling method is applied to evaluate the tunneling splitting of tropolone due to the intramolecular proton transfer in the electronic excited state, first time, in a framework of the trajectory on-the-fly molecular dynamics (TOF-MD) approach. To prevent unphysical zero-point vibrational energy transfer among the normal modes of vibration, quantum zero-point vibrational energies are assigned only to the vibrational modes related to intramolecular proton transfer, whereas the remaining modes are treated as bath modes. Practical ways to determine the tunnel-initiating points and tunneling path are introduced. It is shown that the tunneling splitting decreases as the bath-mode energy increases. The experimental tunneling splitting value is well reproduced by the present TOF-MD approach based on the Wentzel-Kramers-Brillouin (WKB) approximation.     

Intramolecular proton transfer in malonaldehyde: accurate multilayer multi-configurational time-dependent Hartree calculations

Full-dimensional (multilayer) multi-configurational time-dependent Hartree calculations studying the intramolecular proton transfer in malonaldehyde based on a recent potential energy surface (PES) [Wang et al., J. Chem. Phys. 128, 224314 (2008)] are presented. The most accurate calculations yield a ground state tunneling splitting of 23.8 cm(-1) and a zero point energy of 14,678 cm(-1). Extensive convergence tests indicate an error margin of the quantum dynamics calculations for the tunneling splitting of about 0.2 cm(-1). These results are to be compared with the experimental value of the tunneling splitting of 21.58 cm(-1) and results of Monte Carlo calculations of Wang et al. on the same PES which yielded a zero point energy of 14,677.9 cm(-1) with statistical errors of 2-3 cm(-1) and a tunneling splitting of 21.6 cm(-1). The present data includes contributions resulting from the vibrational angular momenta to the tunneling splitting and the zero point energy of 0.2 cm(-1) and 2.4 cm(-1), respectively, which have been computed using a perturbative approach.