Accurate analytic potentials for Li2(X Σ1g+) and Li2(A Σ1u+) from 2 to 90 Å, and the radiative lifetime of Li(2p)

Accurate Determination of Blackbody Radiation Shifts in a Strontium Molecular Lattice Clock

Molecular lattice clocks enable the search for new physics, such as fifth forces or temporal variations of fundamental constants, in a manner complementary to atomic clocks. Blackbody radiation (BBR) is a major contributor to the systematic error budget of conventional atomic clocks and is notoriously difficult to characterize and control. Here, we combine infrared Stark-shift spectroscopy in a molecular lattice clock and modern quantum chemistry methods to characterize the polarizabilities of the Sr_{2} molecule from dc to infrared. Using this description, we determine the static and dynamic blackbody radiation shifts for all possible vibrational clock transitions to the 10^{-16} level. This constitutes an important step toward millihertz-level molecular spectroscopy in Sr_{2} and provides a framework for evaluating BBR shifts in other homonuclear molecules.

Galván I, Alavi A, Aleotti F, Aquilante F, Autschbach J, Avagliano D, Baiardi A, Bao JJ, Battaglia S, Birnoschi L, Blanco-González A, Bokarev SI, Broer R, Cacciari R, Calio PB, Carlson RK, Carvalho Couto R, Cerdán L, Chibotaru LF, Chilton NF, Church JR, Conti I, Coriani S, Cuéllar-Zuquin J, Daoud RE, Dattani N, Decleva P, de Graaf C, Delcey M, De Vico L, Dobrautz W, Dong SS, Feng R, Ferré N, Filatov(Gulak) M, Gagliardi L, Garavelli M, González L, Guan Y, Guo M, Hennefarth MR, Hermes MR, Hoyer CE, Huix-Rotllant M, Jaiswal VK, Kaiser A, Kaliakin DS, Khamesian M, King DS, Kochetov V, Krośnicki M, Kumaar AA, Larsson ED, Lehtola S, Lepetit MB, Lischka H, López Ríos P, Lundberg M, Ma D, Mai S, Marquetand P, Merritt ICD, Montorsi F, Mörchen M, Nenov A, Nguyen VHA, Nishimoto Y, Oakley MS, Olivucci M, Oppel M, Padula D, Pandharkar R, Phung QM, Plasser F, Raggi G, Rebolini E, Reiher M, Rivalta I, Roca-Sanjuán D, Romig T, Safari AA, Sánchez-Mansilla A, Sand AM, Schapiro I, Scott TR, Segarra-Martí J, Segatta F, Sergentu DC, Sharma P, Shepard R, Shu Y, Staab JK, Straatsma TP, Sørensen LK, Tenorio BNC, Truhlar DG, Ungur L, Vacher M, Veryazov V, Voß TA, Weser O, Wu D, Yang X, Yarkony D, Zhou C, Zobel JP, Lindh R. The OpenMolcas Web: A Community-Driven Approach to Advancing Computational Chemistry

The developments of the open-source OpenMolcas chemistry software environment since spring 2020 are described, with a focus on novel functionalities accessible in the stable branch of the package or via interfaces with other packages. These developments span a wide range of topics in computational chemistry and are presented in thematic sections: electronic structure theory, electronic spectroscopy simulations, analytic gradients and molecular structure optimizations, ab initio molecular dynamics, and other new features. This report offers an overview of the chemical phenomena and processes OpenMolcas can address, while showing that OpenMolcas is an attractive platform for state-of-the-art atomistic computer simulations.

MLRNet: Combining the Physics-Motivated Potential Models with Neural Networks for Intermolecular Potential Energy Surface Construction

A physics-based machine learning model called MLRNet has been developed to construct the high-accuracy two-body intermolecular potential energy surface (IPES). The outputs of the neural network are integrated into the physically realistic Morse/long-range (MLR) function, which ensures that the MLRNet has meaningful extrapolation at both short and long ranges and solves the asymptotic problem in common neural network potential (NNP) models. The neural network representation of the MLR parameters is more flexible and more efficient than the polynomial expansion in the conventional mdMLR model, especially for systems containing nonrigid monomer(s). The present work illustrates the basic framework of the current MLRNet model, including (i) how to combine the physically meaningful MLR function with different possible NN structures, (ii) the preservation of permutation symmetry, and (iii) the predetermination of the long-range function u_LR. We choose two realistic systems to demonstrate the performance of MLRNet: the three-dimensional IPES of CO₂-He including the CO₂ antisymmetric vibration Q₃ and the six-dimensional IPES of the H₂O-Ar system. In both cases, the fitting errors of the MLRNet are several times smaller than those of the conventional mdMLR model. Both short-range and long-range extrapolation tests were performed to illustrate the extrapolation ability of the MLRNet and its damping function version. Moreover, for the 6-D H₂O-Ar system, the MLRNet only needs 1596 trainable parameters, which is almost equal to the number needed for the 5-D mdMLR model (1509) and half that needed for the PIP-NN model (3501) within similar accuracy, which illustrates the model efficiency in high-dimensional IPES fitting.

The ab initio potential energy curves of atom pairs and transport properties of high-temperature vapors of Cu and Si and their mixtures with He, Ar, and Xe gases

The potential energy curves (PECs) for the homonuclear He-He, Ar-Ar, Cu-Cu, and Si-Si dimers, as well as heteronuclear Cu-He, Cu-Ar, Cu-Xe, Si-He, Si-Ar, and Si-Xe dimers, are obtained in quantum Monte Carlo (QMC) calculations. It is shown that the QMC method provides the PECs with an accuracy comparable with that of the state-of-the-art coupled cluster singles and doubles with perturbative triples corrections [CCSD(T)] calculations. The QMC data are approximated by the Morse long range (MLR) and (12-6) Lennard-Jones (LJ) potentials. The MLR and LJ potentials are used to calculate the deflection angles in binary collisions of corresponding atom pairs and transport coefficients of Cu and Si vapors and their mixtures with He, Ar, and Xe gases in the range of temperature from 100 K to 10 000 K. It is shown that the use of the LJ potentials introduces significant errors in the transport coefficients of high-temperature vapors and gas mixtures. The mixtures with heavy noble gases demonstrate anomalous behavior when the viscosity and thermal conductivity can be larger than that of the corresponding pure substances. In the mixtures with helium, the thermal diffusion factor is found to be unusually large. The calculated viscosity and diffusivity are used to determine parameters of the variable hard sphere and variable soft sphere molecular models as well as parameters of the power-law approximations for the transport coefficients. The results obtained in the present work include all information required for kinetic or continuum simulations of dilute Cu and Si vapors and their mixtures with He, Ar, and Xe gases.

More about properties of Morse oscillator

Using Birge-Sponer extrapolation we have analyzed the approximation of the potential of a real diatomic molecule by the Morse model, which implies a constant value of anharmonicity ωx. The real values of ωx*(v) for each vibrational level are estimated from transition frequencies between neighboring levels. The dependence of ωx* on the vibrational quantum number v up to dissociation is calculated from the literature data for the ground electronic state of H₂, O₂, Be₂, Li₂, ArXe, Xe₂, Kr₂ and the excited state of Li₂. Characteristic features of deviations of the anharmonicity parameter x* - x from the Morse model are described.

Reaction microscope for investigating ionization dynamics of weakly bound alkali dimers

We report on the implementation of a far-off-resonant, optical dipole force trap in a reaction microscope combined with a magneto-optical trap. Kinematically complete multi-photon ionization experiments were performed on optically trapped ⁶Li atoms and photo-associated ⁶Li₂ molecules in their highest vibrational state. The apparatus allows us to distinguish different ionization mechanisms related to the presence of the IR field of the optical dipole trap that can occur during ionization of ⁶Li and ⁶Li₂ in strong fields. In a series of proof-of-principle experiments, we detect weakly bound dimers via three-photon ionization with femtosecond pulses (τ = 30 fs) at a central wavelength of 780 nm and measure directly the momenta of the photoelectrons in coincidence with recoil ions.

Accurate semiempirical potential energy curves for the a3Σ+-state of NaCs, KCs, and RbCs

A five parameter semiempirical Tang-Toennies type model is used to describe the potential curves of the a³Σ⁺-state of the heteronuclear polar molecules NaCs, KCs, and RbCs. These molecules are of current interest in experiments at ultra-cold conditions to explore the effects of the strong dipole-dipole forces on the collective many-body quantum behavior. New quantum phenomena are also anticipated in systems consisting of atomic species with different fermion/boson statistics. The model parameters are obtained by simultaneously fitting all five of the parameters to the extensive LIF-Fourier transform spectroscopy published by Tiemann and collaborators [e.g., Docenko et al. J. Phys. B: At., Mol. Opt. Phys. 39, S929-S943 (2006)], who also report best fit potential curves. Although the new potentials are in good agreement with the earlier potentials, they have the advantage that they are continuous over the entire range of internuclear distances and have the correct long-range behavior. The scattering lengths for all isotope combinations show good agreement with dedicated experiments where available. The new potentials are also in excellent agreement with combining rules based on the potentials of the homonuclear systems.

Transition Strength Measurements to Guide Magic Wavelength Selection in Optically Trapped Molecules

Optical trapping of molecules with long coherence times is crucial for many protocols in quantum information and metrology. However, the factors that limit the lifetimes of the trapped molecules remain elusive and require improved understanding of the underlying molecular structure. Here we show that measurements of vibronic line strengths in weakly and deeply bound ^{88}Sr_{2} molecules, combined with ab initio calculations, allow for unambiguous identification of vibrational quantum numbers. This, in turn, enables the construction of refined excited potential energy curves, informing the selection of magic wavelengths that facilitate long vibrational coherence. We demonstrate Rabi oscillations between far-separated vibrational states that persist for nearly 100 ms.

Explicitly correlated ab initio potential energy surface and predicted rovibrational spectra for H2O-N2 and D2O-N2 complexes

An ab initio intermolecular potential energy surface (PES) for the van der Waals complex of H₂O-N₂ that explicitly incorporates the intramolecular Q₂ bending normal mode of the H₂O monomer is presented. The electronic structure computations have been carried out at the explicitly correlated coupled cluster theory [CCSD(T)-F12] with an augmented correlation-consistent triple zeta basis set and an additional bond function. Analytic five-dimensional intermolecular PESs for ν₂(H₂O) = 0 and 1 are obtained by fitting to the multi-dimensional Morse/long-range potential function form. These fits to 40 890 points have the root-mean-square (rms) discrepancy of 0.88 cm^-1 for interaction energies less than 2000.0 cm^-1. The resulting vibrationally averaged PESs provide good representations of the experimental microwave and infrared data: for microwave transitions of H₂O-N₂, the rms discrepancy is only 0.0003 cm^-1, and for infrared transitions of the A₁ symmetry of the H₂O(ν₂ = 1 ← 0)-N₂, the rms discrepancy is 0.001 cm^-1. The calculated infrared band origin shifts associated with the ν₂ bending vibration of water are 2.210 cm^-1 and 1.323 cm^-1 for H₂O-N₂ and D₂O-N₂, respectively, in good agreement with the experimental values of 2.254 cm^-1 and 1.266 cm^-1. The benchmark tests and comparisons of the predicted spectral properties are carried out between CCSD(T)-F12a and CCSD(T)-F12b approaches.

Analytic intermolecular potential energy surface and first-principles prediction of the rotational profiles for a symmetric top ion-atom complex: A case study of H3O+-Ar

We presented the first three-dimensional (3D) ab initio intermolecular potential energy surface (PES) for the H₃O⁺-Ar complex. The electronic structure computations were carried out at the explicitly correlated coupled cluster theory-F12 with an augmented correlation-consistent triple zeta basis set. Analytic 3D PES was obtained by least-squares fitting the multi-dimensional Morse/Long-Range (mdMLR) potential model to interaction energies, where the mdMLR function form was applied to the nonlinear ion-atom case for the first time. The 3D PES fitting to 1708 points has root-mean-square deviations of 0.19 cm^-1 with only 108 parameters for interaction energies less than 500 cm^-1. With the 3D PES of the H₃O⁺-Ar complex, we employed the combined radial discrete variable representation/angular finite basis representation method and Lanczos algorithm to calculate rovibrational energy levels. The rotational profiles of the O-H anti-stretching vibrational bands of v₃ ⁺(S)←0⁺ and v₃ ^-(A)←0^- for the H₃O⁺-Ar complex were predicted and were in good agreement with the experimental results.

Full quantum calculation of the rovibrational states and intensities for a symmetric top-linear molecule dimer: Hamiltonian, basis set, and matrix elements

The rovibrational energy levels and intensities of the CH₃F-H₂ dimer have been obtained using our recent global intermolecular potential energy surface [X.-L. Zhang et al., J. Chem. Phys. 148, 124302 (2018)]. The Hamiltonian, basis set, and matrix elements are derived and given for a symmetric top-linear molecule complex. This approach to the generation of energy levels and wavefunctions can readily be utilized for studying the rovibrational spectra of other van der Waals complexes composed of a symmetric top molecule and a linear molecule, and may readily be extended to other complexes of nonlinear molecules and linear molecules. To confirm our method, the rovibrational levels of the H₂O-H₂ dimer have been computed and shown to be in good agreement with experiment and with previous theoretical results. The rovibrational Schrödinger equation has been solved using a Lanczos algorithm together with an uncoupled product basis set. As expected, dimers containing ortho-H₂ are more strongly bound than dimers containing para-H₂. Energies and wavefunctions of the discrete rovibrational levels of CH₃F-paraH₂ complexes obtained from the direct vibrationally averaged 5-dimensional potentials are in good agreement with the results of the reduced 3-dimensional adiabatic-hindered-rotor (AHR) approximation. Accurate calculations of the transition line strengths for the orthoCH₃F-paraH₂ complex are also carried out, and are consistent with results obtained using the AHR approximation. The microwave spectrum associated with the orthoCH₃F-orthoH₂ dimer has been predicted for the first time.

Force Field for Water over Pt(111): Development, Assessment, and Comparison

Metal/water interfaces are key in many natural and industrial processes, such as corrosion, atmospheric, or environmental chemistry. Even today, the only practical approach to simulate large interfaces between a metal and water is to perform force-field simulations. In this work, we propose a novel force field, GAL17, to describe the interaction of water and a Pt(111) surface. GAL17 builds on three terms: (i) a standard Lennard-Jones potential for the bonding interaction between the surface and water, (ii) a Gaussian term to improve the surface corrugation, and (iii) two terms describing the angular dependence of the interaction energy. The 12 parameters of this force field are fitted against a set of 210 adsorption geometries of water on Pt(111). The performance of GAL17 is compared to several other approaches that have not been validated against extensive first-principles computations yet. Their respective accuracy is evaluated on an extended set of 802 adsorption geometries of H₂O on Pt(111), 52 geometries derived from icelike layers, and an MD simulation of an interface between a c(4 × 6) Pt(111) surface and a water layer of 14 Å thickness. The newly developed GAL17 force field provides a significant improvement over previously existing force fields for Pt(111)/H₂O interactions. Its well-balanced performance suggests that it is an ideal candidate to generate relevant geometries for the metal/water interface, paving the way to a representative sampling of the equilibrium distribution at the interface and to predict solvation free energies at the solid/liquid interface.

Analytic Morse/long-range potential energy surfaces and "adiabatic-hindered-rotor" treatment for a symmetric top-linear molecule dimer: A case study of CH3F-H2

A first effective six-dimensional ab initio potential energy surface (PES) for CH₃F-H₂ which explicitly includes the intramolecular Q₃ stretching normal mode of the CH₃F monomer is presented. The electronic structure computations have been carried out at the explicitly correlated coupled cluster level of theory [CCSD(T)-F12a] with an augmented correlation-consistent triple zeta basis set. Five-dimensional analytical intermolecular PESs for ν₃(CH₃F) = 0 and 1 are then obtained by fitting the vibrationally averaged potentials to the Morse/Long-Range (MLR) potential function form. The MLR function form is applied to the nonlinear molecule-linear molecule case for the first time. These fits to 25 015 points have root-mean-square deviations of 0.74 cm^-1 and 0.082 cm^-1 for interaction energies less than 0.0 cm^-1. Using the adiabatic hindered-rotor approximation, three-dimensional PESs for CH₃F-paraH₂ are generated from the 5D PESs over all possible orientations of the hydrogen monomer. The infrared and microwave spectra for CH₃F-paraH₂ dimer are predicted for the first time. These analytic PESs can be used for modeling the dynamical behavior in CH₃F-(H₂)_N clusters, including the possible appearance of microscopic superfluidity.

Characteristic footprints of an exceptional point in the dynamics of Li dimer under a laser field

Non-hermitian quantum mechanics is a formalism that excels in describing time-dependent states such as resonances. As one, it opens up a window to explore new and undiscovered phenomena. Under this formalism coalescence of two eigenstates and a deficient spectrum are a possible situation. These situations are unique and can occur solely in specific conditions known as Exceptional Points (EPs). An EP holds unique characteristics. One of which is a switch-like behavior: upon adiabatically changing the conditions in a closed loop around the EP, the population of one resonance can be transferred completely to another resonance. The phenomenon was not experimentally observed in an atomic or molecular system so far, although experiments involving nonlinear PT symmetry optics and microwave cavities have already indicated its existence. In this work, we demonstrate and confirm that the switch-like behavior exists in the spectrum of a lithium dimer taking into account both the rotations and the vibrations of the system. Moreover, a footprint of the EP is also shown to exist in the photo-association process of the lithium dimer. In this process, the EP's resonances serve as the mean to associate two free lithium atoms into a dimer. Based on this, we suggest a corresponding experiment to demonstrate for the first time the EP phenomenon in a molecular system.

A Density Functional with Spherical Atom Dispersion Terms

A new hybrid density functional, APF, is introduced, which avoids the spurious long-range attractive or repulsive interactions that are found in most density functional theory (DFT) models. It therefore provides a sound baseline for the addition of an empirical dispersion correction term, which is developed from a spherical atom model (SAM). The APF-D empirical dispersion model contains nine adjustable parameters that were selected based on a very small training set (15 noble gas dimers and 4 small hydrocarbon dimers), along with two computed atomic properties (ionization potential and effective atomic polarizability) for each element. APF-D accurately describes a large portion of the potential energy surfaces of complexes of noble gas atoms with various diatomic molecules involving a wide range of elements and of dimers of small hydrocarbons, and it reproduces the relative conformational energies of organic molecules. The accuracy for these weak interactions is comparable to that of CCSD(T)/aug-cc-pVTZ calculations. The accuracy in predicting the geometry of hydrogen bond complexes is competitive with other models involving DFT and empirical dispersion.

A new four-dimensional ab initio potential energy surface for N2O-He and vibrational band origin shifts for the N2O-He(N) clusters with N = 1-40

A new four-dimensional ab initio potential energy surface for N(2)O-He is constructed at the CCSD(T) level with an aug-cc-pVQZ basis set together with bond functions. The vibrational coordinates Q(1) and Q(3) of N(2)O are explicitly included, due to the strong coupling between the symmetric and asymmetric stretches of N(2)O. A global potential energy surface is obtained by fitting the original potential points to a four-dimensional Morse∕long range (MLR) analytical form. In the fitting, the ab initio noise in the long range region of the potential is smoothed over by theoretically fixed long range parameters. Two-dimensional intermolecular potentials for both the ground and the excited υ(3) states of N(2)O are then constructed by vibrationally averaging the four-dimensional potential. Based on the two-dimensional potentials, we use the path integral Monte Carlo algorithm to calculate the vibrational band origin shifts for the N(2)O-He(N) clusters using a first order perturbation theory estimate. The calculated shifts agree reasonably well with the experimental values and reproduce the evolution tendency from dimer to large clusters.