Wave Packet Method Research Articles

State-resolved quantum dynamics and reaction mechanisms of the Ne + HD+(v = 1) reaction

In the work, accurate converged theoretical simulations of the Ne + HD + ( v = 1 ) reaction are presented with the time-dependent wave packet (TDWP) method using the accurate potential energy surface by Lv (J. Chem. Phys. 2010, 132, 014303). In particular, the effect of intramolecular substitution and the translational collision energy on the reactive dynamics is investigated in detail to discern the reaction mechanisms. Total and state-to-state integral cross sections (ICSs) and differential cross sections (DCSs) of the two channels are derived and compared over a wide range of collision energies up to 1.2 eV. The analysis of the ICSs in terms of the total angular momentum contributions shows that the contribution of the J partial wave varies in the two reactive channels at different collision energies. The product state-resolved DCSs reveal that two opposite mechanisms co-exist in the two reactive channels with different ratios at low and high collision energies. Overall, NeH + + D channel appears to form the products via collision complex formation, which survives long enough leading to a forward–backward symmetric total DCS at low collision energies. In comparison, the NeD + + H proceeds through a direct stripping reaction, with predominantly forward scattering in the DCSs.

Study on Quantum Dynamics of the Na + Na2 Exchanged Reaction and Lifetime Prediction of Na3 Complex Based on the Neural Network Potential Energy Surface.

A high-precision global potential energy surface (PES) is constructed for the Na3 system based on high-level ab initio calculations and the fundamental invariant neural network (FI-NN) method. The root-mean-square error (RMSE) of the PES is 2.88 cm-1. The state-resolved quantum dynamics of the ground-state Na + Na2 (v = 0, j = 0) → Na2 (v', j') + Na reaction is studied using the time-dependent wave packet (TDWP) method on the new PES. Analysis of the relevant integral cross sections revealed a complicated energy-transfer mechanism during collisions. Similarly, the characteristics of the differential cross sections indicate that the complex-forming mechanism plays a dominant role in the reaction, providing conditions for a comprehensive exploration of the lifetimes of the complexes. Based on the Rice-Ramsperger-Kassel-Marcus (RRKM) theory, the calculated lifetime of the Na3 complex is approximately 3.9 ns.

Ab Initio Neural Network Potential Energy Surface and Quantum Dynamics Calculations on Na(2S) + H2 → NaH + H Reaction.

The collisions between Na atoms and H2 molecules are of great significance in the field of chemical reaction dynamics, but the corresponding dynamics results of ground-state reactions have not been reported experimentally or theoretically. Herein, a global and high-precision potential energy surface (PES) of NaH2 (12A') is constructed by the neural network model based on 21,873 high-level ab initio points. On the newly constructed PES, the quantum dynamics calculations on the Na(2S) + H2(v0 = 0, j0 = 0) → NaH + H reaction are carried out using the time-dependent wave packet method to study the microscopic reaction mechanism at the state-to-state level. The calculated results show that the low-vibrational products are mainly formed by the dissociation of the triatomic complex; whereas, the direct reaction process dominates the generation of the products with high-vibrational states. The reaction generally follows the direct H-abstraction process, and there is also the short-lived complex-forming mechanism that occurs when the collision energy exceeds the reaction threshold slightly. The PES could be used to further study the stereodynamics effects of isotope substitution and rovibrational excitations on the title reaction, and the presented dynamics data would provide an important reference on the corresponding experimental research at a higher level.

Non-adiabatic dynamics studies of the C+(2P1/2, 3/2) + H2 reaction: Based on global diabatic potential energy surfaces of CH2.

Global diabatic potential energy surfaces (PESs) of CH2+ are constructed using the neural network method with a specific function based on 18 213 abinitio points. The multi-reference configuration interaction method with the aug-cc-pVQZ basis set is adopted to perform the abinitio calculations. The topographical properties of the diabatic PESs are examined in detail. In general, the diabatic PESs provide an accurate quasi-diabatic representation. To validate the diabatic PESs, the dynamics studies of the C+(2P1/2, 3/2) + H2 (v0 = 0, j0 = 0) → H + CH+(X1Σ+) reaction are performed using the time-dependent wave packet method. The reaction probabilities, integral cross sections, differential cross sections, and rate constants are calculated and compared with the experimental and theoretical results. Non-adiabatic dynamics results are in good agreement with experimental data. In addition, the non-adiabatic effect in the C+(2P1/2, 3/2) + H2 reaction is significant due to the non-adiabatic results being obviously larger than adiabatic values. The reasonable non-adiabatic dynamics results indicate that present diabatic PESs can be recommended for any type of dynamics study.

A Globally Accurate Neural Network Potential Energy Surface and Quantum Dynamics Studies on Be+(2S) + H2/D2 → BeH+/BeD+ + H/D Reactions.

Chemical reactions between Be+ ions and H2 molecules have significance in the fields of ultracold chemistry and astrophysics, but the corresponding dynamics studies on the ground-state reaction have not been reported because of the lack of a global potential energy surface (PES). Herein, a globally accurate ground-state BeH2+ PES is constructed using the neural network model based on 18,657 ab initio points calculated by the multi-reference configuration interaction method with the aug-cc-PVQZ basis set. On the newly constructed PES, the state-to-state quantum dynamics calculations of the Be+(2S) + H2(v0 = 0; j0 = 0) and Be+(2S) + D2(v0 = 0; j0 = 0) reactions are performed using the time-dependent wave packet method. The calculated results suggest that the two reactions are dominated by the complex-forming mechanism and the direct abstraction process at relatively low and high collision energies, respectively, and the isotope substitution has little effect on the reaction dynamics characteristics. The new PES can be used to further study the reaction dynamics of the BeH2+ system, such as the effects of rovibrational excitations and alignment of reactant molecules, and the present dynamics data could provide an important reference for further experimental studies at a finer level.

Atomistic wave packet investigation of phonon scattering at rough surfaces

Investigating the phonon-surface scattering mechanism is essential for the evaluation of thermal transport in nanostructures. The theoretical description to quantify this mechanism remains to be developed at the atomic scale. This work presents a phonon wave packet method to study the phonon-surface scattering behavior at rough surfaces. We obtain the specularity distribution dependent on phonon polarization, wavelength, and surface roughness. The reflection of the diffuse wave packet is primarily attributed to the surface shadowing effect at a higher incident angle, surface disorder, and surface-induced localized modes. Taking the wavevector-dependent specularity data as input, the thermal conductivity of silicon nanowires is calculated based on the Boltzmann Transport Equation. Our specularity model provides an accurate evaluation for predicting thermal conductivity. This work offers an atomic-level analysis for phonon-surface interaction, which is helpful for the understanding of thermal transport in nanostructures.

Product state-resolved reactive scattering calculations using stair shaped grids in hyperspherical coordinates for the quantum wave packet method

Product state-resolved reactive scattering calculations using stair shaped grids in hyperspherical coordinates for the quantum wave packet method

A Neural Network Potential Energy Surface and Quantum Dynamics Study of Ca(1S) + H2(v 0 = 0, j 0 = 0) → CaH + H Reaction.

The reactive collision between Ca and H2 molecules has attracted great interest experimentally due to the key role of the product CaH molecule in the field of astrophysics and cold molecules. However, quantum dynamics calculations for this system have not been reported due to the lack of a global potential energy surface (PES). Herein, a globally accurate PES of the ground-state CaH2 is developed by combining 11365 high-level ab initio points and permutation invariant polynomial neural network method. Based on the newly constructed PES, the state-to-state quantum dynamics calculations for the Ca(1S) + H2 (v 0 = 0, j 0 = 0) → CaH + H reaction are carried out using the time-dependent wave packet method. The dynamic results reveal that the reaction follows the complex-forming mechanism near the reactive threshold, whereas both the indirect insertion mechanism and direct abstraction mechanism have effects at higher collision energies. The newly constructed PES can be used to further study the influence of isotope substitution, rovibrational excitation, and spatial orientation of reactant molecules on reaction dynamics.

Quantum Dynamics Studies of the Li + Na2 (V = 0, j = 0) → Na + NaLi Reaction on a New Neural Network Potential Energy Surface.

The ultracold reaction offers a unique opportunity to elucidate the intricate microscopic mechanism of chemical reactions, and the Na2Li system serves as a pivotal reaction system in the investigation of ultracold reactions. In this work, a high-precision potential energy surface (PES) of the Na2Li system is constructed based on high-level ab initio energy points and the neural network (NN) method, and a proper asymptotic functional form is adopted for the long-range interaction, which is suitable for the study of cold or ultracold collisions. Based on the new NN PES, the dynamics of the Li + Na2 (v = 0, j = 0) → Na + NaLi reaction are studied in the collision energy range of 10-7 to 80 cm-1. In the high collision energy range of 8 to 80 cm-1, the dynamics of the reaction is studied using the time-dependent wave packet method and the statistical quantum mechanical (SQM) method. Comparing the results of the two methods, it is found that the SQM method provides a rough description of the product ro-vibrational state distribution but overestimates the integral cross-section values. With the decrease of collision energy, the reaction differential cross section gradually changes from forward-backward symmetric scattering to predominant forward scattering. In the low collision energy range from 10-7 to 8 cm-1, the SQM method is used to study the reaction dynamics, and the rate constant in the Wigner threshold region is estimated to be 2.87 × 10-10 cm3/s.

Quantum stereodynamics studies for the Li + HF (v = 0, j = 1) → H + LiF reaction

In reactive collisions, the colliding direction has an important effect on the reaction dynamics, and the study of reaction stereodynamics which controls the initial orientation or alignment of reactants is of great significance for controlling chemical reactions. In this work, the influence of different initial reactant orientations on the Li + HF (v0 = 0, j0 = 1) → LiF + H reaction dynamics are studied by using the quantum wave packet method. We found that the initial orientation of HF molecule has a great influence on the reactivity and scattering direction of product, but has little influence on the product rotational-vibrational state distribution.

Quantum dynamics of the P(2D) + H2 (X1[formula omitted]) → PH (X3Σ-) + H(2S) reaction

The title reaction is thought to be important for the astronomical PO formation, which may contribute to the generation of bioavailable H3PO3. Using quasi-classical trajectory (QCT) method, the state-specified (v = 0, j = 0) rate constant of this reaction was calculated. The comparison with previous temperature-averaged QCT rate constant shows that the effect of ro-vibrational excitation is relatively small below 2000 K. Then, the time-dependent wave packet (TDWP) method was applied to investigate the quantum reaction mechanisms, obtaining the state-specified (v = 0, j = 0) reaction probability, integral cross sections and rate constants. The number of partial waves and the integral interval of collision energy considered were rigorously tested to ensure the convergence of the TDWP rate constant below 2000 K. The TDWP rate constant agrees better to the latest experimental results at 291 to 685 K than those from the two QCT calculations. Also, the TDWP calculations can produce relatively reliable rate constants at 700 to 2000 K, where the Arrhenius formula fitted by experimental results ignores the temperature dependence of the activation energy. The calculated rate constants are helpful for modelling the phosphorus chemistry and determining the abundances of the P-bearing species in interstellar media (ISM).

Quantum study of the CH3+ photodissociation in full-dimensional neural network potential energy surfaces.

C H 3 + , a cornerstone intermediate in interstellar chemistry, has recently been detected for the first time by using the James Webb Space Telescope. The photodissociation of this ion is studied here. Accurate explicitly correlated multi-reference configuration interaction abinitio calculations are done, and full-dimensional potential energy surfaces are developed for the three lower electronic states, with a fundamental invariant neural network method. The photodissociation cross section is calculated using a full-dimensional quantum wave packet method in heliocentric Radau coordinates. The wave packet is represented in angular and radial grids, allowing us to reduce the number of points physically accessible, requiring to push up the spurious states appearing when evaluating the angular kinetic terms, through projection technique. The photodissociation spectra, when employed in astrochemical models to simulate the conditions of the Orion bar, result in a lesser destruction of CH3+ compared to that obtained when utilizing the recommended values in the kinetic database for astrochemistry.

Influence of spatial coherence on phonon transmission across aperiodically arranged interfaces

In this study, we employ the atomistic wave-packet method to directly simulate coherent phonon transport and scattering dynamics in an aperiodic superlattice structure with aperiodically arranged interfaces. Our investigation reveals that interference dynamics, contingent upon the relative magnitudes of phonon wavelength, spatial coherence length, and average interface spacing, profoundly influence transmission spectra and mode-conversion behavior. Particularly, longer-wavelength phonons possessing larger coherence lengths exhibit more pronounced destructive interference and mode-conversion, leading to generally lower transmission. The insights gained into coherent phonon wave physics from this study can significantly augment the analysis and design of phononic devices.

A comparative study on adiabatic and nonadiabatic dynamics of the H(2S) + NaH(X1Σ+) reaction

To reveal the influence of nonadiabatic couplings on the H(2S) + NaH(X1Σ+) reaction, a comparative study on adiabatic and nonadiabatic dynamics is conducted by using the time-dependent wave packet method. The dynamics results show that nonadiabatic couplings reduce the reactivity of H2(X1Σg+) forming channel and vibrational state density of the product, while the exchange reaction channel is almost unaffected. Although the nonadiabatic couplings have a great influence on the product state distribution, it has little influence on the scattering direction of the product. For the H2(X1Σg+) forming channels, products tend to forward scattering. For the exchange reaction channel, the NaH(X1Σ+) product exhibits a maximum scattering peak in the backward direction. The study reveals that the effect of nonadiabatic couplings on the ro-vibration distribution of the H2(X1Σg+) forming channel in the H(2S) + NaH(X1Σ+) reaction is significantly different from the effect in the H(2S) + LiH(X1Σ+) reaction.

Quantum Dynamics of O(3P)+HBr→OH+Br Reaction: Integral Cross Sections and Rate Constants and their Initial State Dependence.

The integral cross sections (ICSs) and rate constants (RCs) for the title reaction are calculated by means of accurate quantum wave packet method employing the recent ab initio potential energy surface developed by Peterson and co-workers. Hundreds of partial wave contributions (up to J=150) are calculated explicitly taking Coriolis coupling into account. The ICSs are found to increase monotonically with energy and their energy dependences are smooth except in the energy range below the potential energy barrier, where several resonance peaks are observed. The temperature dependences of the RCs obey in general the Arrhenius law. These behaviors manifest the dominance of abstract mechanism. Initial rotational state (ji) excitations are found to greatly enhance the reactivity starting from ji=4, while exhibiting a complicated dependence for ji<4. From Boltzmann average of the initial state specified RCs (for ji=0-15) we obtained thermal RCs which are in reasonable agreement with available experimental results.

The influence of stereodynamical control on the nonadiabatic effects in the D + HD (v = 1, j = 2) → D2 + H reaction.

Stereodynamics is a field that studies the influence of the alignment or orientation of colliding partners on the results of collisions. At present, the intersection of nonadiabatic effects and stereodynamics remains to be explored. In this study, we theoretically demonstrate significant stereodynamical effects in the D + HD (v = 1, j = 2) → D2 + H reaction within the collision energy range of 0.01-2.99eV by using the time-dependent wave packet method. It is found that the stereodynamical control not only facilitates the reaction but also allows precise control of the products over a range of different scattering angles. The analysis at the state-to-state level reveals that the nonadiabatic effects are stronger in the parallel configuration than in the perpendicular configuration. By topological approach to separate the two reaction pathways at the conical intersection, the scattering amplitude of the roaming pathway in the parallel configuration is larger than that of the perpendicular configuration, which leads to more dramatic nonadiabatic features in the collision with parallel configuration.

Reactant-Product Decoupling Technique Using the Intermediate Coordinate Method.

Although the reactant-product decoupling (RPD) technique was proposed over two decades ago, it remains an efficient approach for calculating product state-resolved information on some simple direct reactions using the quantum wave packet method. In the past, usually the RPD technique employed the collocation method to transform the wave function between reactant and product arrangements, which requires quite large computational efforts. In this work, the intermediate coordinate (IC) method is employed to realize the RPD technique. Numerical examples demonstrate that this new IC RPD (IRPD) technique has superior computational efficiency compared with the original method employing the collocation method. Especially, the new IRPD technique significantly saves disk space and computer memory. To illustrate the features of our new method, the total reaction probabilities of the H + H2, H + Br2, and F + H2 reactions with J = 0 and the differential cross sections of the H + H2 and F + H2 reactions at a series of collision energy are calculated and presented. With this efficient and effective new RPD technique, the Li + HF reaction, which involves sharp resonances with long-range wave functions in the van der Waals wells in both the reactant and product arrangements, is also calculated with several J at the product state-resolved level to reveal the ability of the RPD technique for describing resonance wave functions. With these numerical examples, it is found that, for the reaction with resonances, the RPD approach should be applied carefully. Otherwise, it is very possible that the resonances could disappear with the application of the RPD technique.

A new global potential energy surface of CaH2 system constructed with neural network method and dynamics studies of the Ca(4s21S) + H2 reaction

A new global potential energy surface (PES) for the ground state of CaH2 system was constructed by using neural network method. This paper provides a comprehensive discussion on the topographical characteristics of the PES and the spectral properties of stationary points. Furthermore, a time-dependent wave packet method was used to perform the dynamics studies on the Ca(4 s2 1S) + H2 reaction based on the newly constructed PES. Some meaningful dynamics results are reported, including reaction probabilities, integral cross sections and differential cross sections. In addition, the reaction mechanism of the Ca(4 s2 1S) + H2 reaction was discussed in detail.

Product State-Resolved Reactive Scattering Studies of the H + Cl2 (v0 = 1-3, j0 = 0) → HCl + Cl Reaction by the Time-Dependent Wave Packet Method.

The typical hydrogen atom plus halogen molecule reaction H + Cl2 → HCl + Cl has implications across many fields. In this paper, product state-resolved quantum dynamics calculations for the vibrationally excited reaction H + Cl2 (v0 = 1-3, j0 = 0) → HCl + Cl were conducted using the time-dependent wave packet method on a newly developed accurate potential energy surface. Numerical results indicate that the initial vibrational excitation of Cl2 does enhance the reactivity for this early barrier reaction, although less than the enhancement of the translational energy. The calculated product vibrational state-resolved integral cross sections and rate constants reveal that the product vibrational state distribution and the initial vibrational state of Cl2 are highly correlated. The thermal rate constant in the temperature range from 100 to 1000 K was given and is found to be in reasonable agreement with the experimental measurements.

A globally accurate neural network potential energy surface and quantum dynamics study of Mg+(2S) + H2 → MgH+ + H reaction

A globally accurate ground-state MgH2+ potential energy surface (PES) is constructed by neural network method based on 12,056 high-level ab initio points. On the PES, the dynamics calculations of the Mg+(2S) + H2/HD (v0 = 0, j0 = 0) reaction are performed using the time-dependent wave packet method. The branching ratio (MgD+/MgH+) in the Mg+ + HD reaction is much larger than one in most studied collision energies. The dynamics results imply that the title reaction follows the complex-forming mechanism near the threshold, and the direct abstraction process gradually plays the dominant role as the collision energy increases.