Search

This study presents a parallel algorithm for high-dimensional quantum dynamics simulations in poly atomic reactions, integrating distributed- and shared-memory models. The distributions of the wave function and potential energy matrix across message passing interface processes are based on bundled radial and angular dimensions, with implementations featuring either two- or one-sided communication schemes. Using realistic parameters for the H + NH3 reaction, performance assessment reveals linear scalability, exceeding 90% efficiency with up to 600 processors. In addition, owing to the universal and concise structure, the algorithm demonstrates remarkable extensibility to diverse reaction systems, as demonstrated by successes with six-atom and four-atom reactions. This work establishes a robust foundation for high-dimensional dynamics studies, showcasing the algorithm's efficiency, scalability, and adaptability. The algorithm's potential as a valuable tool for unraveling quantum dynamics complexities is underscored, paving the way for future advancements in the field. [ABSTRACT FROM AUTHOR]

We present a comprehensive investigation into the dissociative chemisorption of HOD on a rigid Ni(100) surface using an approximate full-dimensional (9D) quantum dynamics approach, which was based on the time-dependent wave-packet calculations on a full-dimensional potential energy surface obtained through neural network fitting to density functional theory energy points. The approximate-9D probabilities were computed by averaging the seven-dimensional (7D) site-specific dissociation probabilities across six impact sites with appropriate relative weights. Our results uncover a distinctive bond-selective effect, demonstrating that the vibrational excitation of a specific bond substantially enhances the cleavage of that excited bond. The product branching ratios are substantially influenced by which bond undergoes excitation, exhibiting a clear preference for the product formed through the cleavage of the excited bond over the alternative product. [ABSTRACT FROM AUTHOR]

Polanyi's rules predict that a late-barrier reaction yields vibrationally cold products; however, experimental studies showed that the H2 product from the late-barrier H + H2O(|04⟩−) and H + HOD(vOH = 4) reactions is vibrationally hot. Here, we report a potential-averaged five-dimensional state-to-state quantum dynamics study for the H + HOD(vOH = 0–4) → H2 + OD reactions on a highly accurate potential energy surface with the total angular momentum J = 0. It is found that with the HOD vibration excitation increasing from vOH = 1 to 4, the product H2 becomes increasingly vibrationally excited and manifests a typical characteristic of an early barrier reaction for vOH = 3 to 4. Analysis of the scattering wave functions revealed that vibrational excitation in the breaking OH bond moves the location of dynamical saddle point from product side to reactant side, transforming the reaction into an early barrier reaction. Interestingly, pronounced oscillatory structures in the total and product vibrational-state-resolved reaction probabilities were observed for the H + HOD(vOH = 3, 4) reactions, in particular at low collision energies, which originate from the Feshbach resonance states trapped in the bending/torsion excited vibrational adiabatic potential wells in the entrance region due to van der Waals interactions. [ABSTRACT FROM AUTHOR]

The 6D time-dependent wave packet calculations were performed to explore H2 dissociation on Ag, Au, and two AgAu alloy surfaces, using four newly fitted potential energy surfaces based on the neural network fitting to density functional theory energy points. The ligand effect resulting from the Ag–Au interaction causes a reduction in the barrier height for H2+Ag/Au(111) compared to H2+Ag(111). However, the scenario is reversed for H2+Au/Ag(111) and H2+Au(111). The 6D dissociation probabilities of H2 on Ag/Au(111) surfaces are significantly higher than those on the pure Ag(111) surface, but the corresponding results for H2 on Au/Ag(111) surfaces are substantially lower than those on the pure Au(111) surface. The reactivity of H2 on Au(111) is larger than that on Ag(111), despite Ag(111) having a slightly lower static barrier height. This can be attributed to the exceptionally small dissociation probabilities at the hcp and fcc regions, which are at least 100 times smaller compared to those at the bridge or top site for H2+Ag(111). Due to the late barrier being more pronounced, the vibrational excitation of H2 on Ag(111) is more effective in promoting the reaction than on Au(111). Moreover, a high degree of alignment dependence is detected for the four reactions, where the H2 dissociation has the highest probability at the helicopter alignment, as opposed to the cartwheel alignment. [ABSTRACT FROM AUTHOR]

Tunneling splittings observed in molecular rovibrational spectra are significant evidence for tunneling motion of hydrogen nuclei in water clusters. Accurate calculations of the splitting sizes from first principles require a combination of high-quality inter-atomic interactions and rigorous methods to treat the nuclei with quantum mechanics. Many theoretical efforts have been made in recent decades. This Perspective focuses on two path-integral based tunneling splitting methods whose computational cost scales well with the system size, namely, the ring-polymer instanton method and the path-integral molecular dynamics (PIMD) method. From a simple derivation, we show that the former is a semiclassical approximation to the latter, despite that the two methods are derived very differently. Currently, the PIMD method is considered to be an ideal route to rigorously compute the ground-state tunneling splitting, while the instanton method sacrifices some accuracy for a significantly smaller computational cost. An application scenario of such a quantitatively rigorous calculation is to test and calibrate the potential energy surfaces of molecular systems by spectroscopic accuracy. Recent progress in water clusters is reviewed, and the current challenges are discussed. [ABSTRACT FROM AUTHOR]

The mode-specific dynamics for the dissociative chemisorption of H2O on rigid Ni(100) is investigated by approximate nine-dimensional (9D) quantum dynamics calculations. The vibrational state-specific 9D dissociation probabilities are obtained by site-averaging the site-specific seven-dimensional results based on an accurate full-dimensional potential energy surface newly developed by neural network fitting to density functional theory energy points with the revised version of the Perdew, Burke, and Ernzerhof functional. The mode specificity of H2O/Ni(100) is very different from that of H2O/Ni(111) or H2O/Cu(111) whose reactivity enhancement by vibrational excitations is quite efficient. For H2O/Ni(100), it is found that the excitation in the symmetric stretching mode is more efficacious than increasing the translational energy in promoting the reaction, while the excitations in the asymmetric stretching mode and bending mode are less efficacious than the translational energy at low collision energies. These interesting observations can be attributed to the near central-barrier reaction for H2O/Ni(100), as well as large discrepancies between the site-specific mode specificities at different impact sites. The mode-specific dynamics obtained in this study is different from that obtained with the reaction path Hamiltonian approach, indicating the importance of full-dimensional quantum dynamics for gas–surface reactions. [ABSTRACT FROM AUTHOR]

Potential-averaged fivedimensional (PA5D) time-dependent wave packet calculations have been carried out to obtain the initial state-specific reaction probabilities and rate constants for a series of reactant rotationally excited states of the title reaction on one of our recently developed potential energy surface (F12a PES). It is found that individual excitation of the reactant HC1 to j1=2 or OH to j2=1 yields the largest total reaction probability except for low collision energies. In the low energy region, the resonance enhancement of reaction probability decreases rapidly with the increase of OH rotational excitation, resulting in a rapid decline in the rate constant below 250 K. The rate constant of the simultaneous HC1 and OH rotational excitation is not a simple superposition of the results of the separate excitation. Finally, we obtain the thermal rate constants of the reaction by Boltzmann averaging the rate constants for initial rotational states up to (5, 5), which are in reasonable agreement with the experimental results. [ABSTRACT FROM AUTHOR]

The reactivity and selectivity of bimetallic surfaces are of fundamental importance in industrial applications. Here, we report the first six-dimensional (6D) quantum dynamics study for the role of surface strain and ligand effects on the reactivity of HCl on a strained pseudomorphic monolayer of Au deposited onto a Ag(111) substrate, with the aid of accurate machine learning-based potential energy surfaces. The substitute of Au into Ag changes the location of the transition state; however, the static barrier height remains roughly the same as pure Au(111). The 6D quantum dynamics calculations reveal that the surface strain due to lattice expansion slightly enhances the reactivity. The ligand effect due to electronic structure interactions between Au and Ag substantially suppresses the reactivity of HCl in the ground vibrational state but promotes the reactivity via vibrational excitation at high kinetic energies. This finding can be attributed to more close interaction with Ag atoms at the transition state close to the fcc site, as well as the tight transition-state region, making the vibrational excitation highly efficient in enhancing the reactivity. Our study quantitatively unravels the dynamical origin of reactivity control by two metals, which will ultimately provide valuable insight into the selectivity of the catalyst. [ABSTRACT FROM AUTHOR]

In this work, the permutation invariant polynomial neural network (PIP-NN) approach is employed to construct a quasi-diabatic Hamiltonian for system with non-Abelian symmetries. It provides a flexible and compact NN-based diabatic ansatz from the related approach of Williams, Eisfeld, and co-workers. The example of H 3 + is studied, which is an (E + A) × (e + a) Jahn–Teller and Pseudo-Jahn–Teller system. The PIP-NN diabatic ansatz is based on the symmetric polynomial expansion of Viel and Eisfeld, the coefficients of which are expressed with neural network functions that take permutation-invariant polynomials as input. This PIP-NN-based diabatic ansatz not only preserves the correct symmetry but also provides functional flexibility to accurately reproduce ab initio electronic structure data, thus resulting in excellent fits. The adiabatic energies, energy gradients, and derivative couplings are well reproduced. A good description of the local topology of the conical intersection seam is also achieved. Therefore, this diabatic ansatz completes the PIP-NN based representation of DPEM with correct symmetries and will enable us to diabatize even more complicated systems with complex symmetries. [ABSTRACT FROM AUTHOR]

A full-dimensional ab initio potential energy surface (PES) and dipole moment surface (DMS) for the ammonia dimer (NH3)2 are reported. The database of the PES consists of 27 736 ab initio energy points and all of these points were calculated at the UCCSD(T)-F12a/AVTZ level. The PES was fitted by using the fundamental invariant neural network (FI-NN) method that satisfies the permutational symmetry of identical atoms, and the root mean square fitting error for the PES is very small as low as 0.562 meV. The geometries for the (NH3)2 DMS are the same as those used for the PES and are calculated at the XYG3/AVTZ level. This PES can describe a variety of internal floppy motions, including all kinds of vibrational modes no matter intermolecular or intramolecular. The CCSD(T)-PES can dissociate correctly to two NH3 monomers, with De = 1135.55 cm−1 (13.58 kJ/mol) which agrees accurately with the 13.5 ± 0.3 kJ/mol predicted by previous work. [ABSTRACT FROM AUTHOR]

A novel flying cascaded acceleration mechanism is proposed to generate energetic proton beams with multi-peak energy spectra using a circularly polarized (CP) Laguerre–Gaussian (LG) laser pulse in three-dimensional particle-in-cell simulations. Simulations show that the protons are initially accelerated and compressed into the beam center via the radiation pressure of the CP LG (σz = −1) laser pulse. Then, they are tailored by flying dipolar electric fields in this LG laser, resulting in a multi-peak energy spectrum. Each shaped proton peak exhibits a narrow energy spread of ∼5% and high flux of ∼2 × 108 protons/MeV at giga-electron volts energy. Such a flying cascaded acceleration mechanism extends the energy spectra of proton beams from monoenergetic to multi-peak structure, thereby potentially enhancing the generation efficiency of monoenergetic proton beams for various applications, such as proton-induced spallation reactions, proton radiography, and proton therapy. [ABSTRACT FROM AUTHOR]

Collimation control of energetic plasma beams is crucial in the laser–plasma field. In this paper, we report on a self-collimated acceleration scheme for a plasma beam using an ultra-intense Laguerre–Gaussian (LG) laser irradiating a solid target. Three-dimensional (3D) particle-in-cell simulations show that a plasma beam with a high current density is stably formed by the radiation pressure of the hollow LG laser. The initial interaction of LG laser with solid target can be approximately researched by a deformable mirror model. Under the effect of the ponderomotive force of the LG laser, the plasma converges in the center axis to form a narrow beam. An elongated strong-magnetic tunnel (B ∼ 2 kT) is self-generated around the plasma beam, capable of trapping some electrons in a region with a radius of less than 500 nm (r < 500 nm). Compared with the case driven by the conventional Gaussian laser, the beam radius size is dramatically reduced from the microscale to hundreds on the nanoscale. The beam density is increased by at least ten times. Such an interesting scheme can provide a feasible and efficient way to achieve and enhance the collimation of energetic particle beams, which may benefit the general applications of fast ignition in inertial fusion, radiotherapy, realization of high-energy density states, and so on. [ABSTRACT FROM AUTHOR]

A full-dimensional quantum dynamical study for the bimolecular reactions of hydrogen molecules with amino radicals for different isotopologues is reported. The nonreactive amino radical is described by two Radau vectors that are very close to the valence bond coordinates. Potential-optimized discrete variable representation basis is used for the vibrational coordinates of the amino radical. Starting from the reaction H2 + NH2, we study the isotope effects for the two reagents separately, i.e., H2 + NH2/ND2/NHD and H2/D2/HD + NH2. The effects of different vibrational mode excitations of the reagents on the reactivities are studied. Physical explanations about the isotope effects are also provided thoroughly including the influence of vibrational energy differences between the different isotopologues and the impact of the tunneling effect. [ABSTRACT FROM AUTHOR]

The collisions transferring large portions of energy are often called supercollisions. In the H + C2H2 reactive system, the rovibrationally cold C2H2 molecule can be activated with substantial internal excitations by its collision with a translationally hot H atom. It is interesting to investigate the mechanisms of collisional energy transfer in other important reactions of H with hydrocarbons. Here, an accurate, global, full-dimensional potential energy surface (PES) of H + C2H4 was constructed by the fundamental invariant neural network fitting based on roughly 100 000 UCCSD(T)-F12a/aug-cc-pVTZ data points. Extensive quasi-classical trajectory calculations were carried out on the full-dimensional PES to investigate the energy transfer process in collisions of the translationally hot H atoms with C2H4 in a wide range of collision energies. The computed function of the energy-transfer probability is not a simple exponential decay function but exhibits large magnitudes in the region of a large amount of energy transfer, indicating the signature of supercollisions. The supercollisions among non-complex-forming nonreactive (prompt) trajectories are frustrated complex-forming processes in which the incoming H atom penetrates into C2H4 with a small C–H distance but promptly and directly leaves C2H4. The complex-forming supercollisions, in which either the attacking H atom leaves (complex-forming nonreactive collisions) or one of the original H atoms of C2H4 leaves (complex-forming reactive trajectories), dominate large energy transfer from the translational energy to internal excitation of molecule. The current work sheds valuable light on the energy transfer of this important reaction in the combustion and may motivate related experimental investigations. [ABSTRACT FROM AUTHOR]

We present a new approach to incorporate the geometric phase in the time-dependent wave packet calculations based on the analytic diabatic potential energy matrices for two-state systems connecting via a conical intersection. The approach only requires information on the location of the conical intersection and the adiabatic potential energy surface of the ground electronic state and merely takes the same computational cost as a diabatic calculation. Demonstrations of the benchmark H + H2/HD reactions show that the new approach can accurately include the geometric phase in dynamics calculation and can be easily extended to the cold regime where the GP effects become more pronounced. Due to its simplicity and numerical efficiency, the new approach has the potential to extend the dynamics study of the geometric effects to a wide range of reaction systems. [ABSTRACT FROM AUTHOR]

Symmetry adaptation is crucial in representing a permutationally invariant potential energy surface (PES). Due to the rapid increase in computational time with respect to the molecular size, as well as the reliance on the algebra software, the previous neural network (NN) fitting with inputs of fundamental invariants (FIs) has practical limits. Here, we report an improved and efficient generation scheme of FIs based on the computational invariant theory and parallel program, which can be readily used as the input vector of NNs in fitting high-dimensional PESs with permutation symmetry. The newly developed method significantly reduces the evaluation time of FIs, thereby extending the FI-NN method for constructing highly accurate PESs to larger systems beyond five atoms. Because of the minimum size of invariants used in the inputs of the NN, the NN structure can be very flexible for FI-NN, which leads to small fitting errors. The resulting FI-NN PES is much faster on evaluating than the corresponding permutationally invariant polynomial-NN PES. [ABSTRACT FROM AUTHOR]

A full-dimensional time-dependent wave packet study using mixed polyspherical Jacobi and Radau coordinates for the title reaction has been reported. The non-reactive moiety CH3 has been described using three Radau vectors, whereas two Jacobi vectors have been used for the bond breaking/formation process. A potential-optimized discrete variable representation basis has been employed to describe the vibrational coordinates of the reagent CH4. About one hundred billion basis functions have been necessary to achieve converged results. The reaction probabilities for some initial vibrational states are given. A comparison between the present approach and other methods, including reduced and full-dimensional ones, is also presented. [ABSTRACT FROM AUTHOR]

The dependence of reactivity on different facets of a surface is an interesting subject in dynamics at gas-surface interfaces. Here, we constructed new six-dimensional (6D) potential energy surfaces (PESs) for the dissociative chemisorption of HCl on rigid Ag(100) and Ag(110) surfaces, using the neural network method based on extensive density functional theory (DFT) calculations with the Perdew-Burke-Ernzerhof (PBE) functional, and compared the two PESs with the previously fitted PES of HCl/Ag(111). Time-dependent wave packet calculations show that the new PESs are very well converged with respect to the fitting procedure as well as to the number of DFT data points. The 6D dissociation probabilities for HCl initially in the ground rovibrational state decrease gradually for HCl/Ag(110), HCl/Ag(100), and HCl/Ag(111), consistent with the increasing barrier heights for the three reactions. The validity of the site-averaging approximation for HCl/Ag(110) does not hold well as compared with HCl/Ag(100) and HCl/Ag(111), in particular, at low kinetic energies, due to the strong steering effect this reaction exhibits if it is modeled with the semilocal PBE functional, which results in a low reaction barrier and a deep physisorption well. [ABSTRACT FROM AUTHOR]

We study the photodissociation dynamics of the hydrogen sulfide cations (H2S+) using the time-sliced velocity map ion imaging (VMI) technique and high-accuracy calculations. High-resolution ion images of the S+(4S) products were measured at four photolysis wavelengths of 325.158, 325.200, 325.243, 325.307 nm, which correspond to the excitation to the A2A1(0,13,0) K=1 state of H2S+. Rotational state-resolved total kinetic energy releases and angular distributions have been derived as a function of the photolysis wavelengths. Notably, photolysis wavelength dependent product rotational state and anisotropy parameter distributions have been clearly observed. Full-dimensional potential energy surface characterization suggests that nonadiabatic coupling between A2A1 and B2B2 states at C2v configurations, as well as relaxation of the symmetry to Cs in the conical intersection region between the two states, plays a key role in the photodissociation process. [ABSTRACT FROM AUTHOR]

The collision between hydrogen and ammonia is a benchmark system to study chemical elementary reactions with five atoms. In this work, we present a description of the system based on mixed Jacobi and Radau coordinates combined with the time-dependent wave packet method to study the H + NH3 reaction. The Radau coordinates are used to describe the reactive moiety NH2. A salient feature of this approach is that the present coordinates have a great advantage that a very small number of basis set functions can be used to describe the NH2 group. Potential-optimized discrete variable representation basis is applied for the vibrational coordinates of the reagent NH3. The reaction probabilities for several initial vibrational states are presented in this paper. The role of the different vibrational excitations on the reactivity is thoroughly described. [ABSTRACT FROM AUTHOR]

We present a theoretical study on plasmonic properties of an Au based nanoshell structure in which the localized surface plasmon resonances (LSPRs) and the inter-band electronic transitions (IBTs) are presented. The optical properties of Au nanoshell are depicted by a corrected dielectric function of Drude-Lorentz mode, and Mie theory is applied to calculate the absorption efficiencies. It shows that there are two localized surface plasmon absorption peaks, which are induced from the outside LSPRs and the inside LSPRs, respectively. The positions of these two peaks can shift by changing core radius and shell thickness. By varying the core radius and thickness of the shells, the LSPRs are modulated with the IBTs, leading to the strong interactions between the LSPRs and the IBTs. It is found that the LSPRs and IBTs can couple with each other and the new resonance absorption peaks can be observed. Our results show that the LSPRs are enhanced by the IBTs. This enhancement can be applied to improve biosensor signals and advanced plasmonic applications. [ABSTRACT FROM AUTHOR]

We carried out six-dimensional quantum dynamics calculations for the dissociative chemisorption of HCl on a rigid Ag(111) surface, employing three potential energy surfaces (PESs) which were recently constructed using the neural network approach based on extensive density functional theory calculations with Perdew-Burke-Ernzerhof, Perdew-Wang91, and revised Perdew-Burke-Ernzerhof functionals, respectively. The vibrational excitation of HCl enhances the reactivity substantially, and the dissociation is most favored for HCl molecules colliding with rotation in a plane parallel to the Ag(111) surface (helicopter alignment). The influence of rotational excitation on the dissociation probability is much more complicated, with different trends at high and at low kinetic energies. The usage of three different PESs does not change the effects of vibrational excitation, rotational excitation, and rotational-alignment qualitatively, but it does change the magnitude of dissociation probabilities quantitatively due to the different barrier heights. [ABSTRACT FROM AUTHOR]

The rate coefficients of the H + H2O2 → H2 + HO2 reaction are calculated using the ring polymer molecular dynamics (RPMD), quasi-classical trajectory (QCT), and canonical variational transition state theory (CVT) with small curvature tunneling (SCT) correction, in conjunction with the recently constructed fundamental invariant-neural network (FI-NN) potential energy surface (PES) [X. Lu et al., Phys. Chem. Chem. Phys. 20, 23095 (2018)]. In RPMD calculations, 32, 16, and 8 beads are used for computing the rate coefficients at 200 K ≤ T ≤ 400 K, 500 K ≤ T ≤ 700 K, and 700 K < T ≤ 1000 K, respectively. Given that the previous experimental rate coefficients vary widely, in particular, at low temperatures, the present RPMD rate coefficients agree well with most of the experimental results. In addition, comparing with some experimental values, the present QCT and CVT/SCT calculations on the FI-NN PES also predict accurate results at some temperatures. These results strongly support the accuracy of the present dynamics calculations as well as the full-dimensional FI-NN PES. [ABSTRACT FROM AUTHOR]

A micrometer-scale "plasma lens" self-constructed by the prepulse and main pulse of the Laguerre–Gaussian (LG) laser is realized to enhance the collimation and acceleration of proton beams in a target normal sheath field acceleration mechanism. Hydrodynamic FLASH and particle-in-cell simulations are carried out and find that a collimated proton source with beam divergence ∼2.7° is generated by the LG laser, which is smaller than the case driven by the traditional Gaussian laser. It demonstrates that the curved sheath field on the "plasma lens" plays an important role in the beam collimation. Such an approach considerably relaxes the constraints of complex design for the target fabrication and auxiliary laser pulse, opening new doors for high-repetition-rate collimated proton accelerations for innovative applications in upcoming high-repetition-rate petawatt laser systems. [ABSTRACT FROM AUTHOR]

Transition state wave packet calculations have been carried out to compute cumulative reaction probabilities for the H2 + OH reaction on the NN1 potential energy surface, as given in the work of Chen et al., from which well converged thermal rate constants for the reaction up to a temperature of 1000 K were obtained. It was found that both the centrifugal sudden approximation and the "J-K-shifting" approximation slightly underestimated the thermal rate constants, while the ring-polymer molecular dynamics overestimated the rates in the low temperature region. After considering the correction of the barrier height by spin-orbit coupling effect and the more accurate level of theory, the calculated rate constants were in good consistency with experimental measurements in the entire temperature region for this benchmark reaction. [ABSTRACT FROM AUTHOR]

We first constructed six-dimensional potential energy surfaces (PESs) for the dissociative chemisorption of HCl on rigid Ag(111). Three global PESs were calculated using the neural network method based on extensive density functional theory (DFT) calculations with Perdew-Burke-Ernzerhof, Perdew-Wang91, and revised Perdew-Burke-Ernzerhof functionals, respectively. The resulting three PESs are all accurately fitted and smooth, based on the small fitting errors and good agreements between the fitted PESs and the direct DFT calculations. Time-dependent wave packet calculations show that the PESs are very well converged with respect to the fitting process, as well as the number of DFT data points. Comparisons of three PESs and time-dependent quantum dynamics results on the three PESs were given and discussed. The validity of site-averaging approximation holds for the title reaction, where the 25 site-averaged four-dimensional dissociation probability can accurately reproduce the six-dimensional dissociation probability on the current three PESs. We expect that the current theoretical investigations can stimulate the experimental work for this gas-surface reaction. [ABSTRACT FROM AUTHOR]

We developed a systematic theoretical study of nonlinear optical properties of semiconductors. The eight-band k⋅p model and the energy-balance equation are employed to calculate the transmission and optical absorption coefficients in the presence of both the linear one-photon absorption and the nonlinear two-photon absorption (TPA) processes. A substantial reduction of the optical transmittance far below the band-gap can be observed under relatively high-intensity radiation fields due to the nonlinear TPA. The TPA-induced optical transmittance decreases with increasing intensity of the radiation fields. Our theoretical results are in line with those observed experimentally. The theoretical approach can be applied to understand the nonlinear optical properties of semiconductors under high-field conditions. [ABSTRACT FROM AUTHOR]

A new set of diabatic potential energy surfaces (PESs) for LiFH is constructed with artificial neural networks (NNs). The adiabatic PESs of the ground state and the first excited state are directly fitted with NNs. Meanwhile, the adiabatic-to-diabatic transformation (ADT) angles (mixing angles) are obtained by simultaneously fitting energy difference and interstate coupling gradients. No prior assumptions of the functional form of ADT angles are used before fitting, and the ab initio data including energy difference and interstate coupling gradients are well reproduced. Converged dynamical results show remarkable differences between adiabatic and diabatic PESs, which suggests the significance of non-adiabatic processes. [ABSTRACT FROM AUTHOR]

The dissociative chemisorption of N2 is the rate-limiting step for ammonia synthesis in industry. Here, we investigated the role of initially vibrational excitation and rotational excitation of N2 for its reactivity on the Fe(111) surface, based on a recently developed six-dimensional potential energy surface. Six-dimensional quantum dynamics study was carried out to investigate the effect of vibrational excitation for incidence energy below 1.6 eV, due to significant quantum effects for this reaction. The effects of vibrational and rotational excitations at high incidence energies were revealed by quasiclassical trajectory calculations. We found that raising the translational energy can enhance the dissociation probability to some extent, however, the vibrational excitation or rotational excitation can promote dissociation more efficiently than the same amount of translational energy. This study provides valuable insight into the mode-specific dynamics of this heavy diatom-surface reaction. [ABSTRACT FROM AUTHOR]

We performed extensive quasiclassical trajectory calculations for the H+C2D2→HD+C2D/D2+C2H reaction based on a recently developed, global and accurate potential energy surface by the fundamental-invariant neural network method. The direct abstraction pathway plays a minor role in the overall reactivity, which can be negligible as compared with the roaming pathways. The acetylene-facilitated roaming pathway dominates the reactivity, with very small contributions from the vinylidene-facilitated roaming. Although the roaming pathways proceed via the short-lived or long-lived complex forming process, the computed branching ratio of product HD to D2 is not far away from 2:1, implying roaming dynamics for this reaction is mainly contributed from the long-lived complex-forming process. The resulting angular distributions for the two product channels are also quite different. These computational results give valuable insights into the significance and isotope effects of roaming dynamics in the biomolecular reactions. [ABSTRACT FROM AUTHOR]

The dissociative chemisorption of HCl on the Au(111) surface has recently been an interesting and important subject, regarding the discrepancy between the theoretical dissociation probabilities and the experimental sticking probabilities. We here constructed an accurate full-dimensional (sixdimensional (6D)) potential energy surface (PES) based on the density functional theory (DFT) with the revised Perdew-Burke-Ernzerhof (RPBE) functional, and performed 6D quantum mechanical (QM) calculations for HCl dissociating on a rigid Au(111) surface. The effects of vibrational excitations, rotational orientations, and site-averaging approximation on the present RPBE PES are investigated. Due to the much higher barrier height obtained on the RPBE PES than on the PW91 PES, the agreement between the present theoretical and experimental results is greatly improved. In particular, at the very low kinetic energy, the QM-RPBE dissociation probability agrees well with the experimental data. However, the computed QM-RPBE reaction probabilities are still markedly different from the experimental values at most of the energy regions. In addition, the QM-RPBE results achieve good agreement with the recent ab initio molecular dynamics calculations based on the RPBE functional at high kinetic energies. [ABSTRACT FROM AUTHOR]

The time-dependent wave packet method has been employed to calculate state-to-state differential cross sections for the title reaction in full dimensions. It is found that the majority of H2O is produced in the first stretching excited states, with a large fraction of available energy for the reaction ending up as product internal motion. The differential cross sections for collision energy up to 0.4 eV are all peaked in the backward direction, but the width of the angular distribution increases considerably as the increase of collision energy. The isotope effect was also examined by comparing the scattering angular distribution for the title reaction with those for the HD + OH and D2 + OH reactions obtained in our previous work. [ABSTRACT FROM AUTHOR]

The time-dependent wave packet method has been employed to calculate the state-to-state integral cross sections and differential cross sections (DCSs) for three initial states of the title reaction on the recently constructed neural network potential energy surface. It is found that the product HBr(υ′ = 2, 3, 4) states have the dominated population in the entire energy region considered here, indicating an inverted HBr vibrational state distribution. More than half of the available energy ends up as product internal motion, and most of which goes into the vibrational motion. Our calculations show that initial rotational excitation of Br2 has little effect on the product ro-vibrational state distributions and DCSs of the reaction. While the initial vibrational excitation has some influences. The initial vibrational excitation to υ0 = 5 obviously enhance the product vibrational excitation in the low energy region. The DCSs for collision energy up to 0.5 eV at the ground and rotationally excited state are peaked in the backward direction, but the width of the angular distribution increases considerably with the increase of collision energy. For the vibrationally excited state, the DCSs are rather complicated with some strong forward scattering peaks for highly vibrationally excited products. [ABSTRACT FROM AUTHOR]

A more flexible neural network (NN) method using the fundamental invariants (FIs) as the input vector is proposed in the construction of potential energy surfaces for molecular systems involving identical atoms. Mathematically, FIs finitely generate the permutation invariant polynomial (PIP) ring. In combination with NN, fundamental invariant neural network (FI-NN) can approximate any function to arbitrary accuracy. Because FI-NN minimizes the size of input permutation invariant polynomials, it can efficiently reduce the evaluation time of potential energy, in particular for polyatomic systems. In this work, we provide the FIs for all possible molecular systems up to five atoms. Potential energy surfaces for OH3 and CH4 were constructed with FI-NN, with the accuracy confirmed by full-dimensional quantum dynamic scattering and bound state calculations. [ABSTRACT FROM AUTHOR]

Full quantum mechanical calculations of vibrational energies of methane and fluoromethane are carried out using a polyspherical description combining Radau and Jacobi coordinates. The Hamiltonian is built in a potential-optimized discrete variable representation, and vibrational energies are solved using an iterative eigensolver. This new approach can be applied to a large variety of molecules. In particular, we show that it is able to accurately and efficiently compute eigenstates for four different molecules : CH4, CHD3, CH2D2, and CH3F. Very good agreement is obtained with the results reported previously in the literature with different approaches and with experimental data. [ABSTRACT FROM AUTHOR]

This work performs a time-dependent wavepacket study of the H2 + C2H â?' H + C2H2 reaction on a new ab initio potential energy surface (PES). The PES is constructed using neural network method based on 68 478 geometries with energies calculated at UCCSD(T)-F12a/aug-cc-pVTZ level and covers H2 + C2H â?" H + C2H2, H + C2H2 â?' HCCH2, and HCCH2 radial isomerization reaction regions. The reaction dynamics of H2 + C2H â?' H + C2H2 are investigated using full-dimensional quantum dynamics method. The initial-state selected reaction probabilities are calculated for reactants in eight vibrational states. The calculated results showed that the H2 vibrational excitation predominantly enhances the reactivity while the excitation of bending mode of C2H slightly inhibits the reaction. The excitations of two stretching modes of C2H molecule have negligible effect on the reactivity. The integral cross section is calculated with J-shift approximation and the mode selectivity in this reaction is discussed. The rate constants over 200-2000 K are calculated and agree well with the experimental measured values. [ABSTRACT FROM AUTHOR]

To fast and accurately compute rate coefficients of the H/D + CH4→H2/HD + CH3 reactions, we propose a segmented strategy for fitting suitable potential energy surface (PES), on which ring-polymer molecular dynamics (RPMD) simulations are performed. On the basis of recently developed permutation invariant polynomial neural-network approach [J. Li et al., J. Chem. Phys. 142, 204302 (2015)], PESs in local configuration spaces are constructed. In this strategy, global PES is divided into three parts, including asymptotic, intermediate, and interaction parts, along the reaction coordinate. Since less fitting parameters are involved in the local PESs, the computational efficiency for operating the PES routine is largely enhanced by a factor of ~20, comparing with that for global PES. On interaction part, the RPMD computational time for the transmission coefficient can be further efficiently reduced by cutting off the redundant part of the child trajectories. For H + CH4, good agreements among the present RPMD rates and those from previous simulations as well as experimental results are found. For D + CH4, on the other hand, qualitative agreement between present RPMD and experimental results is predicted. [ABSTRACT FROM AUTHOR]

Understanding the role of reactant ro-vibrational degrees of freedom (DOFs) in reaction dynamics of polyatomic molecular dissociation on metal surfaces is of great importance to explore the complex chemical reaction mechanism. Here, we present an expensive quantum dynamics study of the dissociative chemisorption of CH4 on a rigid Ni(111) surface by developing an accurate nine-dimensional quantum dynamical model including the DOF of azimuth. Based on a highly accurate fifteen-dimensional potential energy surface built from first principles, our simulations elucidate that the dissociation probability of CH4 has the strong dependence on azimuth and surface impact site. Some improvements are suggested to obtain the accurate dissociation probability from quantum dynamics simulations. [ABSTRACT FROM AUTHOR]

In the present work, we develop a highly accurate, fifteen-dimensional potential energy surface (PES) of CH4 interacting on a rigid flat Ni(111) surface with the methodology of neural network (NN) fit to a database consisted of about 194 208 ab initio density functional theory (DFT) energy points. Some careful tests of the accuracy of the fitting PES are given through the descriptions of the fitting quality, vibrational spectrum of CH4 in vacuum, transition state (TS) geometries as well as the activation barriers. Using a 25-60-60-1 NN structure, we obtain one of the best PESs with the least root mean square errors: 10.11 meV for the entrance region and 17.00 meV for the interaction and product regions. Our PES can reproduce the DFT results very well in particular for the important TS structures. Furthermore, we present the sticking probability S0 of ground state CH4 at the experimental surface temperature using some sudden approximations by Jackson's group. An in-depth explanation is given for the underestimated sticking probability. [ABSTRACT FROM AUTHOR]

A collimated electron sheet is generated by using an intense circularly polarized (CP) Laguerre–Gaussian (LG) laser to irradiate the plasma in three-dimensional particle-in-cell simulations. It is found that the CP L G p l (p = 0 , l = 1 , σ z = 1) laser provides a more stable vortex force to ensure the concentrated transport of electrons as compared to a linearly polarized LG laser. Herein, the effects of the laser intensity, target density, and preplasma on collimated sheet formation are discussed in detail. Such collimated electron sheets may be potentially used in various applications such as the generation of vortex radiation sources and electron microscopy. [ABSTRACT FROM AUTHOR]

The ring polymer molecular dynamics (RPMD) calculations are performed to calculate rate constants for the title reaction on the recently constructed potential energy surface based on permutation invariant polynomial (PIP) neural-network (NN) fitting [J. Li et al., J. Chem. Phys. 142, 204302 (2015)]. By inspecting convergence, 16 beads are used in computing free-energy barriers at 300 K ≤ T ≤ 1000 K, while different numbers of beads are used for transmission coefficients. The present RPMD rates are in excellent agreement with quantum rates computed on the same potential energy surface, as well as with the experimental measurements, demonstrating further that the RPMDis capable of producing accurate rates for polyatomic chemical reactions even at rather low temperatures. [ABSTRACT FROM AUTHOR]

We report a permutationally invariant global potential energy surface (PES) for the H + CH4 system based on ~63 000 data points calculated at a high ab initio level (UCCSD(T)-F12a/AVTZ) using the recently proposed permutation invariant polynomial-neural network method. The small fitting error (5.1 meV) indicates a faithful representation of the ab initio points over a large configuration space. The rate coefficients calculated on the PES using tunneling corrected transition-state theory and quasi-classical trajectory are found to agree well with the available experimental and previous quantum dynamical results. The calculated total reaction probabilities (Jtot = 0) including the abstraction and exchange channels using the new potential by a reduced dimensional quantum dynamic method are essentially the same as those on the Xu-Chen-Zhang PES. [ABSTRACT FROM AUTHOR]

The potential energy surfaces for the H + H2 exchange reaction are determined, using the standard Becke-3_parameter-Lee-Yang-Parr (B3LYP) hybrid exchange-correlation functional and the recently developed XYG3 type (xDH) doubly hybrid functionals. Quantum dynamical results calculated on these density functional theory surfaces revealed that the accuracy for the surfaces from the xDH functionals is quite satisfactory as compared to the benchmark Boothroyd-Keogh-Martin-Peterson-2 surface, whereas B3LYP is wholly inadequate in describing this simplest reaction despite the fact that it has been widely used to study reactions in complex systems. It is shown that further improvements can be achieved if the xDH functional parameters are fine-tuned to be reaction specific. [ABSTRACT FROM AUTHOR]

We employ the initial state-selected time-dependent wave packet approach to an atom-triatom reaction to study the H + HOD → OH + HD/OD + H2 reaction without the centrifugal sudden approximation, based on an accurate potential energy surface which was recently developed by neural network fitting to high level ab initio energy points. The total reaction probabilities and integral cross sections, which are the exact coupled-channel results, are calculated for the HOD reactant initially in the ground and several vibrationally excited states, including the bending excited state, OD stretching excited states, OH stretching excited states, and combined excitations of them. The reactivity enhancements from different initial states of HOD are presented, which feature strong bond-selective effects of the reaction dynamics. The current results for the product branching ratios, reactivity enhancements, and relative cross sections are largely improved over the previous calculations, in quantitatively good agreement with experiment. The thermal rate constant for the title reaction and the contributions from individual vibrational states of HOD are also obtained. [ABSTRACT FROM AUTHOR]

A three dimensional potential energy surface for the F + H2 → HF + H reaction has been computed by the spin unrestricted coupled cluster method with singles, doubles, triples, and perturbative quadruples [UCCSDT(2)Q] using the augmented correlation-consistent polarised valence quadruple zeta basis set for the fluorine atom and the correlation-consistent polarised valence quadruple zeta basis set for the hydrogen atom. All the calculations are based on the restricted open-shell Hartree-Fock orbitals, together with the frozen core approximations, and the UCCSD(T)/complete basis set (CBS) correction term was included. The global potential energy surface was calculated by fitting the sampled ab initio points without any scaling factor for the correlation energy part using a neutral network function method. Extensive dynamics calculations have been carried out on the potential energy surface. The reaction rate constants, integral cross sections, product rotational states distribution, and forward and backward scattering as a function of collision energy of the F + HD → HF + D, F + HD → DF + H, and F + H2 reaction, were calculated by the time-independent quantum dynamics scattering theory using the new surface. The satisfactory agreement with the reported experimental observations previously demonstrates the accuracy of the new potential energy surface. [ABSTRACT FROM AUTHOR]

Eight-dimensional (8D) transition-state wave packet simulations have been performed on two latest potential energy surfaces (PES), the Zhou-Fu-Wang-Collins-Zhang (ZFWCZ) PES [Y. Zhou, B. Fu, C. Wang, M. A. Collins, and D. H. Zhang, J. Chem. Phys. 134, 064323 (2011)] and the Xu-Chen- Zhang (XCZ)-neural networks (NN) PES [X. Xu, J. Chen, and D. H. Zhang, Chin. J. Chem. Phys. 27, 373 (2014)]. Reaction rate constants for both the H+CH4 reaction and the H2+CH3 reaction are calculated. Simulations of the H+CH4 reaction based on the XCZ-NN PES show that the ZFWCZ PES predicts rate constants with reasonable high accuracy for low temperatures while leads to slightly lower results for high temperatures, in line with the distribution of interpolation error associated with the ZFWCZ PES. The 8D H+CH4 rate constants derived on the ZFWCZ PES compare well with full-dimensional 12D results based on the equivalent m-ZFWCZ PES, with a maximum relative difference of no more than 20%. Additionally, very good agreement is shown by comparing the 8D XCZ-NN rate constants with the 12D results obtained on the ZFWCZ-WM PES, after considering the difference in static barrier height between these two PESs. The reaction rate constants calculated for the H2+CH3 reaction are found to be in good consistency with experimental observations. [ABSTRACT FROM AUTHOR]

A new finding of the site-averaging approximation was recently reported on the dissociative chemisorption of the HCl/DCl+Au(111) surface reaction [T. Liu, B. Fu, and D. H. Zhang, J. Chem. Phys. 139, 184705 (2013); J. Chem. Phys. 140, 144701 (2014)]. Here, in order to investigate the dependence of new site-averaging approximation on the initial vibrational state of H2 as well as the PES for the dissociative chemisorption of H2 on Cu(111) surface at normal incidence, we carried out six-dimensional quantum dynamics calculations using the initial state-selected time-dependent wave packet approach, with H2 initially in its ground vibrational state and the first vibrational excited state. The corresponding four-dimensional site-specific dissociation probabilities are also calculated with H2 fixed at bridge, center, and top sites. These calculations are all performed based on two different potential energy surfaces (PESs). It is found that the site-averaging dissociation probability over 15 fixed sites obtained from four-dimensional quantum dynamics calculations can accurately reproduce the six-dimensional dissociation probability for H2 (v = 0) and (v = 1) on the two PESs. [ABSTRACT FROM AUTHOR]

Seven-dimensional time-dependent wave packet calculations have been carried out for the title reaction to obtain reaction probabilities and cross sections for CHD3 in JO= 1, 2 rotationally excited initial states with kO= 0 - JO (the projection of CHD3 rotational angular momentum on its C3 axis). Under the centrifugal sudden (CS) approximation, the initial states with the projection of the total angular momentum on the body fixed axis (KO) equal to kO are found to be much more reactive, indicating strong dependence of reactivity on the orientation of the reagent CHD3 with respect to the relative velocity between the reagents H and CHD3. However, at the coupled-channel (CC) level this dependence becomes much weak although in general the K0 specified cross sections for the KO= kO initial states remain primary to the overall cross sections, implying the Coriolis coupling is important to the dynamics of the reaction. The calculated CS and CC integral cross sections obtained after KO averaging for the JO= 1, 2 initial states with all different k0 are essentially identical to the corresponding CS and CC results for the JO= 0 initial state, meaning that the initial rotational excitation of CHD3 up to JO= 2, regardless of its initial kO, does not have any effect on the total cross sections for the title reaction, and the errors introduced by the CS approximation on integral cross sections for the rotationally excited JO= 1, 2 initial states are the same as those for the JO= 0 initial state. [ABSTRACT FROM AUTHOR]

We present a potential energy surface fitting scheme based on multiplicative artificial neural networks. It has the sum of products form required for efficient computation of the dynamics of multidimensional quantum systems with the multi configuration time dependent Hartree method. Moreover, it results in analytic potential energy matrix elements when combined with quantum dynamics methods using Gaussian basis functions, eliminating the need for a local harmonic approximation. Scaling behavior with respect to the complexity of the potential as well as the requested accuracy is discussed. [ABSTRACT FROM AUTHOR]

The origin of the anomalous frequency dispersion in accumulation capacitance of metal-insulatorsemiconductor devices on InGaAs and InP substrates is investigated using modeling, electrical characterization, and chemical characterization. A comparison of the border trap model and the disorder induced gap state model for frequency dispersion is performed. The fitting of both models to experimental data indicate that the defects responsible for the measured dispersion are within approximately 0.8 nm of the surface of the crystalline semiconductor. The correlation between the spectroscopically detected bonding states at the dielectric/III-V interface, the interfacial defect density determined using capacitance-voltage, and modeled capacitance-voltage response strongly suggests that these defects are associated with the disruption of the III-V atomic bonding and not border traps associated with bonding defects within the high-k dielectric. [ABSTRACT FROM AUTHOR]