Full-dimensional Quantum Dynamics Research Articles

Mode specificity in the HCl + OH → Cl + H2O reaction: Polanyi's rules vs sudden vector projection model.

The dynamics and mode specificity of the HCl + OH → Cl + H2O reaction are investigated using a full-dimensional quantum dynamics method on an accurate global potential energy surface. It is shown that the vibrational excitation of the HCl reactant greatly enhances the reactivity while the OH vibrational excitation has little effect. The surprising HCl vibrational enhancement of this early barrier reaction contradicts a naive extension of Polanyi's rules, but can be explained by the sudden vector projection model, which attributes the promotional effect of the HCl vibration to its strong coupling with the reaction coordinate at the transition state. In addition, it is found that the fundamental and overtone excitations of the HCl reactant change the reaction mechanism from a direct barrier crossing process to a capture-like process.

Loss of Memory in H + CH4 → H2 + CH3 State-to-State Reactive Scattering.

State-to-state reaction probabilities for the H + CH4→ H2 + CH3 reaction are calculated by accurate full-dimensional quantum dynamics calculations employing the multilayer multiconfigurational time-dependent Hartree approach and the quantum transition-state concept. Reactions starting from different vibrational and rotational states of the methane reactant are investigated for vanishing total angular momentum. The vibrational state distributions of the products are found to be essentially independent of the initial rovibrational state of the reactants. The reaction products only show vibrational excitation in the methyl umbrella mode. No excitation in H2 vibration or another CH3 vibration is observed. Analyzing the results, the observed loss of vibrational memory can be explained by a transition-state-based view of the reaction process.

Mode specificity in the HF + OH → F + H2O reaction.

Full-dimensional quantum dynamics and quasi-classical trajectory calculations are reported for the title reaction on a recently constructed ab initio based global potential energy surface. Strong mode specificity was found, consistent with the prediction of the sudden vector projection model. Specifically, the HF vibration strongly promotes the reaction while the OH vibration has little effect. Rotational excitations of both reactants slightly enhance the reaction.

Communication: Ro-vibrational control of chemical reactivity in H+CH₄→ H₂+CH₃: full-dimensional quantum dynamics calculations and a sudden model.

The mode-selective chemistry of the title reaction is studied by full-dimensional quantum dynamics simulation on an accurate ab initio potential energy surface for vanishing total angular momentum. Using a rigorous transition state based approach and multi-configurational time-dependent Hartree wave packet propagation, initial state-selected reaction probabilities for many ro-vibrational states of methane are calculated. The theoretical results are compared with experimental trends seen in reactions of methane. An intuitive interpretation of the ro-vibrational control of the chemical reactivity provided by a sudden model based on the quantum transition state concept is discussed.

A full-dimensional quantum dynamical study of H2+H2 collisions: Coupled-states versus close-coupling formulation

Collision-induced energy transfer involving H2 molecules plays an important role in many areas of physics. Kinetic models often require a complete set of state-to-state rate coefficients for H2+H2 collisions in order to interpret results from spectroscopic observations or to make quantitative predictions. Recent progress in full-dimensional quantum dynamics using the numerically exact close-coupling (CC) formulation has provided good agreement with existing experimental data for low-lying states of H2 and increased the number of state-to-state cross sections that may be reliably determined over a broad range of energies. Nevertheless, there exist many possible initial states (e.g., states with high rotational excitation) that still remain elusive from a computational standpoint even at relatively low collision energies. In these cases, the coupled-states (CS) approximation offers an alternative full-dimensional formulation. We assess the accuracy of the CS approximation for H2+H2 collisions by comparison with benchmark results obtained using the CC formulation. The results are used to provide insight into the orientation effects of the various internal energy transfer mechanisms. A statistical CS approximation is also investigated and cross sections are reported for transitions which would otherwise be impractical to compute.

A full dimensional investigation of infrared spectroscopy of the RbCs dimer using the multi-configuration time-dependent Hartree method

The geometry and infrared absorption spectrum of (RbCs)2 have been studied by full dimensional quantum dynamics simulations. For this purpose, the potential energy and dipole moment surfaces are generated by means of a cluster expansion with all two and three mode correlations, and fitted to analytical expressions with negligible deviations. Accordingly, the ground state (RbCs)2 has a diamond geometry with D(2h) symmetry. The infrared spectrum with frequencies up to 120 cm(-1), exhibits rich details of the fundamentals, overtones, and combination bands; the highest fundamental frequency of (RbCs)2 is only 40.26 cm(-1). The present study unravels important details of the interactions between the widely investigated ultracold RbCs molecules.

Competition between abstraction and exchange channels in H + HCN reaction: Full-dimensional quantum dynamics

Dynamics of the title reaction is investigated on an ab initio based potential energy surface using a full-dimensional quantum wave packet method within the centrifugal sudden approximation. It is shown that the reaction between H and HCN leads to both the hydrogen exchange and hydrogen abstraction channels. The exchange channel has a lower threshold and larger cross section than the abstraction channel. It also has more oscillations due apparently to quantum resonances. Both channels are affected by long-lived resonances supported by potential wells. Comparison with experimental cross sections indicates underestimation of the abstraction barrier height.

Resonances in the Entrance Channel of the Elementary Chemical Reaction of Fluorine and Methane

Extending the fully quantum-state-resolved description of elementary chemical reactions beyond three or four atom systems is a crucial issue in fundamental chemical research. Reactions of methane with F, Cl, H or O are key examples that have been studied prominently. In particular, reactive resonances and nonintuitive mode-selective chemistry have been reported in experimental studies for the F+CH4 →HF+CH3 reaction. By investigating this reaction using transition-state spectroscopy, this joint theoretical and experimental study provides a clear picture of resonances in the F+CH4 system. This picture is deduced from high-resolution slow electron velocity-map imaging (SEVI) spectra and accurate full-dimensional (12D) quantum dynamics simulations in the picosecond regime.

Excitation of helium Rydberg states and doubly excited resonances in strong extreme ultraviolet fields: Full-dimensional quantum dynamics using exponentially tempered Gaussian basis sets

Recently optimized exponentially tempered Gaussian basis sets [P. R. Kapralova-Zdanska and J. Smydke, J. Chem. Phys. 138, 024105 (2013)] are employed in quantitative simulations of helium absorption cross-sections and two-photon excitation yields of doubly excited resonances. Linearly polarized half-infinite and Gaussian laser pulses at wavelengths 38-58 nm and large intensities up to 100 TW/cm(2) are considered. The emphasis is laid on convergence of the results with respect to the quality of the Gaussian basis sets (typically limited by a number of partial waves, density, and spatial extent of the basis functions) as well as to the quality of the basis set of field-free states (typically limited by the maximum rotational quantum number and maximum excitation of the lower electron). Particular attention is paid to stability of the results with respect to varying complex scaling parameter. Moreover, the study of the dynamics is preceded by a thorough check of helium energies and oscillator strengths as they are obtained with the exponentially tempered Gaussian basis sets, being also compared with yet unpublished emission wavelengths measured in electric discharge experiments.

Coupled potential energy surface for the F(2P) + CH4 → HF + CH3 entrance channel and quantum dynamics of the CH4 · F− photodetachment

An approach to construct vibronically and spin-orbit coupled diabatic potential energy surfaces (PESs) which describe all three relevant electronic states in the entrance channels of the X(P) + CH4 →HX + CH3 reactions (with X=F((2)P), Cl((2)P), or O((3)P)) is introduced. The diabatization relies on the permutational symmetry present in the methane molecule and results in diabatic states which transform as the three p orbitals of the X atom. Spin-orbit coupling is easily and accurately included using the atomic spin-orbit coupling matrix of the isolated X atom. The method is applied to the F + CH4 system obtaining an accurate PES for the entrance channel based on ab initio multi-reference configuration interaction (MRCI) calculations. Comparing the resulting PESs with spin-orbit MRCI calculations, excellent agreement is found for the excited electronic states at all relevant geometries. The photodetachment spectrum of CH4·F(-) is investigated via full-dimensional (12D) quantum dynamics calculations on the coupled PESs using the multi-layer multi-configurational time-dependent Hartree approach. Extending previous work [J. Palma and U. Manthe, J. Chem. Phys. 137, 044306 (2012)], which was restricted to the dynamics on a single adiabatic PES, the contributions of the electronically excited states to the photodetachment spectrum are calculated and compared to experiment. Considering different experimental setups, good agreement between experiment and theory is found. Addressing questions raised in the previous work, the present dynamical calculations show that the main contribution to the second peak in the photodetachment spectrum results from electron detachment into the electronically excited states of the CH4F complex.

Mode specificity in the H + H2O → H2 + OH reaction: A full-dimensional quantum dynamics study

The initial state-selected time-dependent wave packet approach to an atom-triatom reaction is employed to study the H + H2O → H2 + OH reaction without the centrifugal sudden approximation. The total reaction probabilities and integral cross sections, which are the exact coupled-channel results, are calculated for the H2O reactant initially in the ground and several vibrationally excited states, including bending excited states, first and second stretching excited states, and simultaneous excitations of both bending and stretching modes. The reactivity enhancements from different initial states of the H2O reagent are presented and discussed in detail. The thermal rate constant for the title reaction and the contributions to this coefficient from individual vibrational states of H2O are also obtained and compared with the previous theoretical and experimental data.

Vibration-vibration and vibration-translation energy transfer in H2-H2 collisions: A critical test of experiment with full-dimensional quantum dynamics

Quantum scattering calculations of vibration-vibration (VV) and vibration-translation (VT) energy transfer for non-reactive H2-H2 collisions on a full-dimensional potential energy surface are reported for energies ranging from the ultracold to the thermal regime. The efficiency of VV and VT transfer is known to strongly correlate with the energy gap between the initial and final states. In H2(v = 1, j = 0) + H2(v = 0, j = 1) collisions, the inelastic cross section at low energies is dominated by a VV process leading to H2(v = 0, j = 0) + H2(v = 1, j = 1) products. At energies above the opening of the v = 1, j = 2 rotational channel, pure rotational excitation of the para-H2 molecule leading to the formation of H2(v = 1, j = 2) + H2(v = 0, j = 1) dominates the inelastic cross section. For vibrationally excited H2 in the v = 2 vibrational level colliding with H2(v = 0), the efficiency of both VV and VT process is examined. It is found that the VV process leading to the formation of 2H2(v = 1) molecules dominates over the VT process leading to H2(v = 1) + H2(v = 0) products, consistent with available experimental data, but in contrast to earlier semiclassical results. Overall, VV processes are found to be more efficient than VT processes, for both distinguishable and indistinguishable H2-H2 collisions confirming room temperature measurements for v = 1 and v = 2.

First principles determination of the NH2/ND2($\skew4\tilde A, \skew3\tilde X$Ã,X̃) branching ratios for photodissociation of NH3/ND3 via full-dimensional quantum dynamics based on a new quasi-diabatic representation of coupled ab initio potential energy surfaces

The A-band photodissociation of ammonia is an archetypical and long studied example of non-adiabatic dynamics in a polyatomic system. Due to a well-known conical intersection seam, electronically excited NH(3) can produce either the ground (X̃(2)B(1)) state or the excited (Ã(2)A(1)) state of the NH(2) fragment. In this work, the non-adiabatic dynamics is investigated using a six-dimensional wave packet method and an improved version of a newly developed diabatic Hamiltonian based on high quality ab initio data. The Ã(2)A(1)/X̃(2)B(1) branching ratios are in excellent agreement with experimental estimates, thus validating the non-adiabatically coupled Hamiltonian.

Quantum Dynamics of the HO + CO → H + CO2 Reaction on an Accurate Potential Energy Surface.

Full-dimensional quantum dynamics of the HO + CO → H + CO2 reaction is investigated on a recent global potential energy surface based on a large number of ab initio points. The J = 0 reaction probability is small and essentially a monotonically increasing function with energy, superimposed by overlapping resonances. The reactivity is considerably enhanced by OH vibrational excitation while relatively insensitive to CO vibrational excitation. The rate constant estimated by the J-shifting approximation indicates a much better agreement with experiment than that obtained on a previous potential energy surface.

Theoretical Investigation of the He4Br2 Conformers

Full dimensional quantum dynamics calculations of the three lowest isomers of the He(4)Br(2) van der Waals molecule in its ground electronic state are reported. The calculations are performed using the multiconfiguration time-dependent Hartree (MCTDH) method and a realistic potential form that includes the sum of three body ab initio coupled-cluster single double triple [CCSD(T)] He-Br(2) interactions plus the He-He and Br-Br interactions. This potential exhibits several multiple minima, with the three lowest ones lying very close in energy, just within 2 cm(-1). Such small differences are also found in the calculated binding energies of the three most stable conformers, indicating the floppiness of the system and, thus, the need of accurate potential forms and quantum full dynamics methods to treat this kind of complexes. The 12 dimensional results reported in this work present benchmark data and, thus, can serve to evaluate approximate methods aiming to describe higher order rare gas-dihalogen (N > 4) complexes. A comparison with previous studies using different potential forms and approaches to the energetics for the He(4)Br(2) cluster is also presented.

Full-dimensional quantum dynamics study of exchange processes for the D + H2O and D + HOD reactions

The exchange processes of D + H(2)O and D + HOD reactions are studied using initial state-selected time-dependent wave packet approach in full dimension. The total reaction probabilities for different partial waves, together with the integral cross sections, are obtained both by the centrifugal sudden (CS) approximation and exact coupled-channel (CC) calculations, for the H(2)O(HOD) reactant initially in the ground rovibrational state. In the CC calculations, small resonance peaks in the reaction probabilities and quick diminishing of the resonance peaks with the increase of total angular momenta J do not lead to clear step-like features just above the threshold in the cross sections for the title reactions, which are different in other isotopically substituted reactions where the hydrogen atom was included as the reactant instead of the deuterium atom [B. Fu, Y. Zhou, and D. H. Zhang, Chem. Sci. 3, 270 (2012); B. Fu and D. H. Zhang, J. Phys. Chem. A 116, 820 (2012)]. It is interesting that the shape resonance-induced features resulting from the reaction tunneling are significantly diminished accordingly in the reactions of the deuterium atom and H(2)O or HOD, owing to the weaker tunneling capability of the reagent deuterium atom in the title reactions than the reagent hydrogen atom in other reactions. In the CS calculations, the resonance peaks persist in many partial waves but cannot survive the partial-wave summations. The cross sections for the D(') + H(2)O → D(')OH + H and D(') + HOD → D(')OD + H reactions are substantially larger than those for the D(') + HOD → HOD(') + D reaction, indicating that the D(')/H exchange reactions are much more favored than the D(')/D exchange.

Rovibrational energy transfer in collisions of H2 with D2: a full-dimensional wave packet propagation study

We present full-dimensional quantum dynamics calculations for the process of rovibrational energy transfer in collisions between H2 and D2 molecules, specifically H2(ν = 1) + D2(ν = 0) → H2(ν = 0) + D2(ν = 1). Rotationally resolved state-to-state cross-sections are obtained for collision energies up to 0.5 eV. From these we calculate rotationally averaged thermal rate coefficients in the temperature range from 100 to 500 K, and compare them with available experimental data. For some transitions, we found it numerically advantageous to compute cross-sections of the reverse collision process and then use microscopic reversibility to obtain the originally sought cross-sections. We employ the Multi-Configuration Time-Dependent Hartree (MCTDH) method for propagating wave packets, and calculate the cross-sections from transition probabilities obtained by the correlation function formalism introduced by Tannor and Weeks. Computations are performed with a potential energy surface that is based on the six-dimensional surface from Boothroyd et al. but reduced in anisotropy, as suggested by Pogrebnya and Clary. The expression of the kinetic energy operator in terms of internal curvilinear coordinates allows us to treat the kinematics of the system exactly, without any decoupling approximations.

Quantum-dynamics study of the ${\rm H}_5^+$H5+ cluster: Full dimensional benchmark results on its vibrational states

A full-dimensional quantum dynamics study is carried out for the highly fluxional H(5)(+) cation on a recent reference potential energy surface by using the multi configuration time-dependent Hartree method. With five equivalent light atoms and shallow barriers between various low-lying stationary points on the surface, the spectroscopic characterization of H(5)(+) represents a huge challenge for accurate quantum dynamics simulations. The present calculation is the first such a study on this cation, which together with its isotope analogies are of primary importance in the interstellar chemistry. The vibrational ground state properties and several vibrationally excited states corresponding to low vibrational frequency motions, not yet directly observable by the experiment, are presented and analyzed.

Shape resonance in the H + D2O → D + HOD reaction: a full-dimensional quantum dynamics study

We present exact coupled-channel (CC) results of full-dimensional quantum mechanical calculations for the H + D2O → D + HOD exchange reaction which provide the first evidence for a shape resonance in this reaction, employing the initial state-selected time-dependent wave packet method with both OD bonds in the D2O reactant treated as reactive bonds. Due to the shallow C3v minimum along the reaction path, the signature of shape resonance is manifested as the existence of a distinct step-like feature in the CC integral cross sections of the exchange reaction just above the threshold. Our study also shows that the current CC integral cross sections are in very good agreement with the experimental results, and a bit larger than those calculated in the previous quantum dynamics study, where the one reactive bond and centrifugal sudden (CS) approximations were invoked, revealing that the one reactive bond approximation overestimated, while the CS approximation underestimated the integral cross sections.

Photodissociation of methyl iodide embedded in a host-guest complex: A full dimensional (189D) quantum dynamics study of CH3I@resorc[4]arene

Accurate full dimensional quantum dynamics calculations studying the photodissociation of CH(3)I@resorc[4]arene on an ab initio based potential energy surface (PES) model are reported. The converged 189D quantum dynamics calculations are facilitated by the multilayer multi-configurational time-dependent Hartree (ML-MCTDH) approach combined with the correlation discrete variable representation (CDVR) for the evaluation of potential energy matrix elements. The potential employed combines an established ab initio PES describing the photodissociation of methyl iodide in the A band with a harmonic description of the resorc[4]arene host and a bilinear modeling of the host-guest interaction. All potential parameters required in the description of the vibrations of the host molecule and the host-guest interaction are derived from ab initio calculations on the host-guest complex. Absorption spectra at 0 K and 300 K are calculated and the electronic population dynamics during the bond breaking process occurring in the first 20-30 fs after the photoexcitation is investigated. Weak but significant effects resulting from the host-guest interaction on this time scale are found and interpreted. The present study demonstrates that accurate fully quantum mechanical dynamics calculations can be preformed for systems consisting of more than 50 atoms using the ML-MCTDH/CDVR approach. Utilizing an efficient statistical approach for the construction of the ensemble of initial wavepackets, these calculations are not restricted to zero temperature but can also study the dynamics at 300 K.