Hydrogen-iodine scattering. II. Rovibronic analysis and collisional dynamics.

Our recently published [Weike et al., J. Chem. Phys. 159, 244119 (2023)] spin-orbit coupled diabatic potential energy model for HI is used in a thorough analysis of bound and quasi-bound states as well as elastic and inelastic processes in H + I collisions. The potential energy model, designed explicitly for studying scattering, accurately describes the various couplings in the system, which lead to complex dynamics. Ro-vibronic bound and quasi-bound states related to the adiabatic electronic ground state and an excited electronic state are analyzed. Calculations using the full 104 × 104 diabatic matrix model or a single adiabatic state are compared in order to investigate approximations in the latter. Elastic and inelastic scattering cross sections as well as thermal rates between the ground and first excited fine structure levels of iodine are computed for collision energies up to 12 500 cm-1. Resonances related to the quasi-bound states are analyzed in terms of their energy, width, lifetime, and decay probabilities. The effect of different resonances on the thermal rates is discussed. Resonances between 30 000 and 40 000 cm-1 are also studied for selected values of the total angular momentum, in particular their decay probabilities into different final states of iodine and hence their potential effect on branching ratios in photodissociation of HI.

Hydrogen-iodine scattering. I. Development of an accurate spin-orbit coupled diabatic potential energy model.

The scattering of H by I is a prototypical model system for light-heavy scattering in which relativistic coupling effects must be taken into account. Scattering calculations depend strongly on the accuracy of the potential energy surface (PES) model. The methodology to obtain such an accurate PES model suitable for scattering calculations is presented, which includes spin-orbit (SO) coupling within the Effective Relativistic Coupling by Asymptotic Representation (ERCAR) approach. In this approach, the SO coupling is determined only for the atomic states of the heavy atom, and the geometry dependence of the SO effect is accounted for by a diabatization with respect to asymptotic states. The accuracy of the full model, composed of a Coulomb part and the SO model, is achieved in the following ways. For the SO model, the extended ERCAR approach is applied, which accounts for both intra-state and inter-state SO coupling, and an extended number of diabatic states are included. The corresponding coupling constants for the SO operator are obtained from experiments, which are more accurate than computed values. In the Coulomb Hamiltonian model, special attention is paid to the long range behavior and accurate c6 dispersion coefficients. The flexibility and accuracy of this Coulomb model are achieved by combining partial models for three different regions. These are merged via artificial neural networks, which also refine the model further. In this way, an extremely accurate PES model for hydrogen iodide is obtained, suitable for accurate scattering calculations.

Simulation of the photodetachment spectra of the nitrate anion (NO3-) in the B̃ 2E' energy range and non-adiabatic electronic population dynamics of NO3.

The photodetachment spectrum of the nitrate anion (NO3-) in the energy range of the NO3 second excited state is simulated from first principles using quantum wave packet dynamics. The prediction at 10 K and 435 K relies on the use of an accurate full-dimensional fully coupled five state diabatic potential model utilizing an artificial neural network. The ability of this model to reproduce experimental spectra was demonstrated recently for the lower energy range [A. Viel, D. M. G. Williams and W. Eisfeld, J. Chem. Phys. 2021, 154, 084302]. Analysis of the spectra indicates a weaker Jahn-Teller coupling compared to the first excited state. The detailed non-adiabatic dynamics is studied by computing the population dynamics. An ultra-fast non-statistical radiationless decay is found only among the Jahn-Teller components, which is followed by a slow statistical non-radiative decay among the different state manifolds. The latter is reproduced perfectly by a simple first order kinetics model. The dynamics in the second excited state is not affected by the presence of a conical intersection with the first excited state manifold.

Accurate quantum dynamics simulation of the photodetachment spectrum of the nitrate anion (NO3 -) based on an artificial neural network diabatic potential model.

The photodetachment spectrum of the nitrate anion (NO3 -) is simulated from first principles using wavepacket quantum dynamics propagation and a newly developed accurate full-dimensional fully coupled five state diabatic potential model. This model utilizes the recently proposed complete nuclear permutation inversion invariant artificial neural network diabatization technique [D. M. G. Williams and W. Eisfeld, J. Phys. Chem. A 124, 7608 (2020)]. The quantum dynamics simulations are designed such that temperature effects and the impact of near threshold detachment are taken into account. Thus, the two available experiments at high temperature and at cryogenic temperature using the slow electron velocity-map imaging technique can be reproduced in very good agreement. These results clearly show the relevance of hot bands and vibronic coupling between the X̃ 2A2 ' ground state and the B̃ 2E' excited state of the neutral radical. This together with the recent experiment at low temperature gives further support for the proper assignment of the ν3 fundamental, which has been debated for many years. An assignment of a not yet discussed hot band line is also proposed.

Diabatic neural network potentials for accurate vibronic quantum dynamics-The test case of planar NO3.

A recently developed scheme to produce high-dimensional coupled diabatic potential energy surfaces (PESs) based on artificial neural networks (ANNs) [D. M. G. Williams and W. Eisfeld, J. Chem. Phys. 149, 204106 (2019)] is tested for its viability for quantum dynamics applications. The method, capable of reproducing high-quality ab initio data with excellent accuracy, utilizes simple coupling matrices to produce a basic low-order diabatic potential matrix as an underlying backbone for the model. This crude model is then refined by making its expansion coefficients geometry-dependent by the output neurons of the ANN. This structure, strongly guided by a straightforward physical picture behind nonadiabatic coupling, combines structural simplicity with high accuracy, reproducing ab initio data without introducing unphysical artifacts to the surface, even for systems with complicated electronic structure. The properties of diabatic potentials obtained by this method are tested thoroughly in the present study. Vibrational/vibronic eigenstates are computed on the X̃ and à states of NO3, a notoriously difficult Jahn-Teller system featuring strong nonadiabatic couplings and complex spectra. The method is investigated in terms of how consistently it produces dynamics results for PESs of similar (fitting) quality and how the results depend on the ANN size and ANN topography. A central aspect of this work is to understand the convergence properties of the new method in order to evaluate its predictive power. A previously developed, high-quality model utilizing a purely (high-order) polynomial ansatz is used as a reference to showcase improvements of the overall quality which can be obtained by the new method.

Quantum dynamics and geometric phase in E⊗e Jahn-Teller systems with general Cnv symmetry.

E ⊗ e Jahn-Teller (JT) systems are considered the prototype of symmetry-induced conical intersections and of the corresponding geometric phase effect (GPE). For decades, this has been analyzed for the most common case originating from C3v symmetry and these results usually were generalized. In the present work, a thorough analysis of the JT effect, vibronic coupling Hamiltonians, GPE, and the effect on spectroscopic properties is carried out for general Cnv symmetric systems (and explicitly for n = 3-8). It turns out that the C3v case is much less general than often assumed. The GPE due to the vibronic Hamiltonian depends on the leading coupling term of a diabatic representation of the problem, which is a result of the explicit n, α, and ß values of a CnvEα ⊗ eß system. Furthermore, the general existence of n/m (m∈N depending on n, α, and ß) equivalent minima on the lower adiabatic sheet of the potential energy surface (PES) leads to tunneling multiplets of n/m states (state components). These sets can be understood as local vibrations of the atoms around their equilibrium positions within each of the local PES wells symmetrized over all equivalent wells. The local vibrations can be classified as tangential or radial vibrations, and the quanta in the tangential mode together with the GPE determine the level ordering within each of the vibronic multiplets. Our theoretical predictions derived analytically are tested and supported by numerical model simulations for all possible Eα ⊗ eß cases for Cnv symmetric systems with n = 3-8. The present interpretation allows for a full understanding of the complex JT spectra of real systems, at least for low excitation energies. This also opens a spectroscopic way to show the existence or absence of GPEs.

Large amplitude motion within acetylene-rare gas complexes hosted in helium droplets.

Near-infrared spectroscopy of the C2H2-Ar, Kr complexes was performed in the spectral region overlapping the ν3/ν2 + ν4 + ν5 Fermi-type resonance of C2H2. The experiment was conducted along the HElium NanoDroplet Isolation (HENDI) technique in order to study the coupling dynamics between a floppy molecular system (C2H2-Ar and C2H2-Kr) and a mesoscopic quantum liquid (the droplet). Calculations were performed using a spectral element based close-coupling program and state-of-the-art 2-dimensional potential energy surfaces to determine the bound states of the C2H2-Ar and C2H2-Kr complexes and simulate the observed spectra. This furnished a quantitative basis to unravel how the superfluid and non-superfluid components of the droplet affect the rotation and the deformation dynamics of the hosted complex.

Application of the spectral element method to the solution of the multichannel Schrödinger equation.

We apply the spectral element method to the determination of scattering and bound states of the multichannel Schrödinger equation. In our approach, the reaction coordinate is discretized on a grid of points whereas the internal coordinates are described by either purely diabatic or locally diabatic (diabatic-by-sector) bases. Bound levels and scattering matrix elements are determined with spectral accuracy using relatively small number of points. The scattering problem is cast as a linear system solved using state-of-the-art sparse matrix non-iterative packages. Boundary conditions can be imposed so as to compute a single column of the matrix solution. A comparison with log-derivative propagators customarily used in molecular physics is performed. The same discretization scheme can also be applied to bound levels that are computed using direct scalable sparse-matrix solvers.

Vibronic eigenstates and the geometric phase effect in the 2E″ state of NO3.

The 2E″ state of NO3, a prototype for the Jahn-Teller effect, has been an enigma and a challenge for a long time for both experiment and theory. We present a detailed theoretical study of the vibronic quantum dynamics in this electronic state, uncovering the effects of tunnelling, geometric phase, and symmetry. To this end, 45 vibronic levels of NO3 in the 2E″ state are determined accurately and analyzed thoroughly. The computation is based on a high quality diabatic potential representation of the two-sheeted surface of the 2E″ state developed by us [W. Eisfeld et al., J. Chem. Phys. 140, 224109 (2014)] and on the multi-configuration time dependent Hartree approach. The vibrational eigenstates of the NO3- anion are determined and analyzed as well to gain a deeper understanding of the symmetry properties of such D3h symmetric systems. To this end, 61 eigenstates of the NO3- anion ground state are computed using the single sheeted potential surface of the 1A1 state published in the same reference quoted above. The assignments of both the vibrational and vibronic levels are discussed. A simple model is proposed to rationalize the computed NO3 spectrum strongly influenced by the Jahn-Teller couplings, the associated geometric phase effect, and the tunnelling. Comparison with the available spectroscopic data is also presented.

Full-dimensional diabatic potential energy surfaces including dissociation: the ²E″ state of NO3.

A scheme to produce accurate full-dimensional coupled diabatic potential energy surfaces including dissociative regions and suitable for dynamical calculations is proposed. The scheme is successfully applied to model the two-sheeted surface of the (2)E″ state of the NO3 radical. An accurate potential energy surface for the NO3⁻ anion ground state is developed as well. Both surfaces are based on high-level ab initio calculations. The model consists of a diabatic potential matrix, which is expanded to higher order in terms of symmetry polynomials of symmetry coordinates. The choice of coordinates is key for the accuracy of the obtained potential energy surfaces and is discussed in detail. A second central aspect is the generation of reference data to fit the expansion coefficients of the model for which a stochastic approach is proposed. A third ingredient is a new and simple scheme to handle problematic regions of the potential energy surfaces, resulting from the massive undersampling by the reference data unavoidable for high-dimensional problems. The final analytical diabatic surfaces are used to compute the lowest vibrational levels of NO3⁻ and the photo-electron detachment spectrum of NO3⁻ leading to the neutral radical in the (2)E″ state by full dimensional multi-surface wave-packet propagation for NO3 performed using the Multi-Configuration Time Dependent Hartree method. The achieved agreement of the simulations with available experimental data demonstrates the power of the proposed scheme and the high quality of the obtained potential energy surfaces.

Low-temperature rate coefficients for vibrational relaxation of (3)Σ(u)(+)Rb2 molecules by (3)He and (4)He atoms.

We present quantum-scattering calculations of (4)He and (3)He colliding with (87)Rb2. For both helium isotopes, the elastic and inelastic rate coefficients are strongly influenced by the J = 1 partial wave. For the lighter isotope, a strong resonance feature of the J = 1 partial wave is responsible for an extremely efficient vibrational relaxation process. We also perform bound-state calculations of the Rb2He complex for even Rb permutation symmetry and nonzero total angular momentum. The global Rb2He (3)Σu(+) potential-energy surface used supports four bound states for (4)He and a single one for (3)He. We propose an analysis of the (87)Rb2(4)He spectrum separating the contributions of Rb2 rotation and helium motion.

Full dimension Rb2He ground triplet potential energy surface and quantum scattering calculations.

We have developed a three-dimensional potential energy surface for the lowest triplet state of the Rb(2)He complex. A global analytic fit is provided as in the supplementary material [see supplementary material at http://dx.doi.org/10.1063/1.4709433 for the corresponding Fortran code]. This surface is used to perform quantum scattering calculations of (4)He and (3)He colliding with (87)Rb(2) in the partial wave J = 0 at low and ultralow energies. For the heavier helium isotope, the computed vibrational relaxation probabilities show a broad and strong shape resonance for a collisional energy of 0.15 K and a narrow Feshbach resonance at about 17 K for all initial Rb(2) vibrational states studied. The broad resonance corresponds to an efficient relaxation mechanism that does not occur when (3)He is the colliding partner. The Feshbach resonance observed at higher collisional energy is robust with respect to the isotopic substitution. However, its effect on the vibrational relaxation mechanism is faint for both isotopes.

Excited Li and Na in He(n): influence of the dimer potential energy curves.

The X(2)Σ ground and the A(2)Π and B(2)Σ first two excited states of Li-He and Na-He are determined using high level complete active space self-consistent field-multireference configuration interaction ab initio method. The obtained potentials differ from the ones proposed by Pascale [Phys. Rev. A 28, 632 (1983)], more strongly for the ground than for the excited states. Quantum diffusion Monte Carlo studies of small Li(∗)He(n) and Na(∗)He(n) with n ≤ 5 are performed using a diatomics-in-molecule approach to model the non-pair additive interaction potential. The sensitivity of our results to the A(2)Π and B(2)Σ potentials used is assessed by an analysis of the structure and of the energetics of the clusters. For these small clusters, the physical conclusions are essentially independent of the diatomic curves employed.

Theoretical study of Rb2 in He(N): potential energy surface and Monte Carlo simulations.

An analytical potential energy surface for a rigid Rb2 in the ³Σ(u)⁺ state interacting with one helium atom based on accurate ab initio computations is proposed. This 2-dimensional potential is used, together with the pair approximation approach, to investigate Rb2 attached to small helium clusters He(N) with N = 1-6, 12, and 20 by means of quantum Monte Carlo studies. The limit of large clusters is approximated by a flat helium surface. The relative orientation of the dialkali axis and the helium surface is found to be parallel. Dynamical investigations of the pendular and of the in-plane rotation of the rigid Rb2 molecule on the surface are presented.

Electronically excited rubidium atom in helium clusters and films. II. Second excited state and absorption spectrum.

Following our work on the study of helium droplets and film doped with one electronically excited rubidium atom Rb(∗) ((2)P) [M. Leino, A. Viel, and R. E. Zillich, J. Chem. Phys. 129, 184308 (2008)], we focus in this paper on the second excited state. We present theoretical studies of such droplets and films using quantum Monte Carlo approaches. Diffusion and path integral Monte Carlo algorithms combined with a diatomics-in-molecule scheme to model the nonpair additive potential energy surface are used to investigate the energetics and the structure of Rb(∗)He(n) clusters. Helium films as a model for the limit of large clusters are also considered. As in our work on the first electronic excited state, our present calculations find stable Rb(∗)He(n) clusters. The structures obtained are however different with a He-Rb(∗)-He exciplex core to which more helium atoms are weakly attached, preferentially on one end of the core exciplex. The electronic absorption spectrum is also presented for increasing cluster sizes as well as for the film.

Multiconfigurational time-dependent Hartree calculations for tunneling splittings of vibrational states: Theoretical considerations and application to malonaldehyde.

Full-dimensional multiconfigurational time-dependent Hartree calculations on the tunneling splitting of the vibrational ground state and the low lying excited states of malonaldehyde are presented. Methodological developments utilizing the symmetry of double well systems for the efficient calculation of tunneling splittings are described and discussed. Important aspects of the theory underlying the previously communicated results for the ground state tunneling splitting [M. D. Coutinho-Neto et al., J. Chem. Phys. 121, 9207 (2004)] are detailed and further developments facilitating the calculation of tunneling splittings for vibrationally excited states are introduced. Utilizing these developments, the 14 lowest vibrational states of malonaldehyde, i.e., seven tunneling splittings, have been computed. The tunneling splittings are found to vary significantly depending on the particular vibrational excitation. This results in a complex pattern of vibrational levels. Studying the dependence of the tunneling splittings on the vibrational excitation, good agreement with available experimental results is found and intuitive interpretations of the results can be given.

Electronically excited rubidium atom in a helium cluster or film.

We present theoretical studies of helium droplets and films doped with one electronically excited rubidium atom Rb( *) ((2)P). Diffusion and path integral Monte Carlo approaches are used to investigate the energetics and the structure of clusters containing up to 14 helium atoms. The surface of large clusters is approximated by a helium film. The nonpair additive potential energy surface is modeled using a diatomic in molecule scheme. Calculations show that the stable structure of Rb( *)He(n) consists of a seven helium atom ring centered at the rubidium, surrounded by a tirelike second solvation shell. A very different structure is obtained when performing a "vertical Monte Carlo transition." In this approach, a path integral Monte Carlo equilibration starts from the stable configuration of a rubidium atom in the electronic ground state adsorbed to the helium surface after switching to the electronically excited surface. In this case, Rb( *)He(n) relaxes to a weakly bound metastable state in which Rb( *) sits in a shallow dimple. The interpretation of the results is consistent with the recent experimental observations [G. Aubock et al., Phys. Rev. Lett. 101, 035301 (2008)].

Linewidths of C2H2 perturbed by H2: experiments and calculations from an ab initio potential.

In this work we present a theoretical and experimental study of the acetylene-hydrogen system. A potential surface considering rigid monomers has been obtained by ab initio quantum chemistry methods. This 4-dimensional potential is further employed to compute, using the close-coupling approach and the coupled-states approximation, pressure broadening coefficients of C2H2 isotropic Raman Q lines over a temperature range of 77 to 2000 K. Experimental data for the acetylene nu2 Raman lines broadened by molecular hydrogen are obtained using stimulated Raman spectroscopy. The comparison of theoretical values with experimental data at 143 K is promising. Approximations to increase the computational efficiency are proposed.

Blueshift and intramolecular tunneling of NH3 umbrella mode in 4He n clusters.

We present diffusion Monte Carlo calculations of the ground and first excited vibrational states of NH(3) (4)He(n) for n< or =40. We use the potential energy surface developed by one of us [M. P. Hodges and R. J. Wheatley, J. Chem. Phys. 114, 8836 (2001)], which includes the umbrella mode coordinate of NH(3). Using quantum Monte Carlo calculations of excited states, we show that this potential is able to reproduce qualitatively the experimentally observed effects of the helium environment, namely, a blueshift of the umbrella mode frequency and a reduction of the tunneling splittings in ground and first excited vibrational states of the molecule. These basic features are found to result regardless of whether dynamical approximations or exact calculations are employed.

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

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