An efficient and flexible approach for computing rovibrational polaritons from first principles.

A theoretical framework is presented for the computation of the rovibrational polaritonic states of a molecule in a lossless infrared (IR) microcavity. In the proposed approach, the quantum treatment of the rotational and vibrational motions of the molecule can be formulated using arbitrary approximations. The cavity-induced changes in electronic structure are treated perturbatively, which allows using the existing polished tools of standard quantum chemistry for determining electronic molecular properties. As a case study, the rovibrational polaritons and related thermodynamic properties of H2O in an IR microcavity are computed for varying cavity parameters, applying various approximations to describe the molecular degrees of freedom. The self-dipole interaction is significant for nearly all light-matter coupling strengths investigated, and the molecular polarizability proved important for the correct qualitative behavior of the energy level shifts induced by the cavity. On the other hand, the magnitude of polarization remains small, justifying the perturbative approach for the cavity-induced changes in electronic structure. Comparing results obtained using a high-accuracy variational molecular model with those obtained utilizing the rigid rotor and harmonic oscillator approximations revealed that as long as the rovibrational model is appropriate for describing the field-free molecule, the computed rovibropolaritonic properties can be expected to be accurate as well. Strong light-matter coupling between the radiation mode of an IR cavity and the rovibrational states of H2O leads to minor changes in the thermodynamic properties of the system, and these changes seem to be dominated by non-resonant interactions between the quantum light and matter.

A quantum-chemical perspective on the laser-induced alignment and orientation dynamics of the CH3 X (X = F, Cl, Br, I) molecules.

Motivated by recent experiments, the laser-induced alignment-and-orientation (A&O) dynamics of the prolate symmetric top CH3 X (X = F, Cl, Br, I) molecules is investigated, with particular emphasis on the effect of halogen substitution on the rotational constants, dipole moments, and polarizabilities of these species, as these quantities determine the A&O dynamics. Insight into possible control schemes for preferred A&O dynamics of halogenated molecules and best practices for A&O simulations are provided, as well. It is shown that for accurate A&O -dynamics simulations it is necessary to employ large basis sets and high levels of electron correlation when computing the rotational constants, dipole moments, and polarizabilities. The benchmark-quality values of these molecular parameters, corresponding to the equilibrium, as well as the vibrationally averaged structures are obtained with the help of the focal-point analysis (FPA) technique and explicit electronic-structure computations utilizing the gold-standard CCSD(T) approach, basis sets up to quintuple-zeta quality, core-correlation contributions and, in particular, relativistic effects for CH3 Br and CH3 I. It is shown that the different A&O behavior of the CH3 X molecules in the optical regime is mostly caused by the differences in their polarizability anisotropy, in other terms, the size of the halogen atom. In contrast, the A&O dynamics of the CH3 X series induced by an intense few-cycle THz pulse is mostly governed by changes in the rotational constants, due to the similar dipole moments of the CH3 X molecules. The A&O dynamics is most sensitive to the B rotational constant: even the difference between its equilibrium and vibrationally-averaged values results in noticeably different A&O dynamics. The contribution of rotational states having different symmetry, weighted by nuclear-spin statistics, to the A&O dynamics is also studied.

Nonadiabatic phenomena in molecular vibrational polaritons.

Nonadiabatic phenomena are investigated in the rovibrational motion of molecules confined in an infrared cavity. Conical intersections (CIs) between vibrational polaritons, similar to CIs between electronic polaritonic surfaces, are found. The spectral, topological, and dynamic properties of the vibrational polaritons show clear fingerprints of nonadiabatic couplings between molecular vibration, rotation, and the cavity photonic mode. Furthermore, it is found that for the investigated system, composed of two rovibrating HCl molecules and the cavity mode, breaking the molecular permutational symmetry, by changing 35Cl to 37Cl in one of the HCl molecules, the polaritonic surfaces, nonadiabatic couplings, and related spectral, topological, and dynamic properties can deviate substantially. This implies that the natural occurrence of different molecular isotopologues needs to be considered when modeling realistic polaritonic systems.

Robust field-dressed spectra of diatomics in an optical lattice.

The absorption spectra of the cold Na2 molecule dressed by a linearly polarized standing laser wave is investigated with a theoretical model incorporating translational, electronic, vibrational as well as rotational degrees of freedom. In such a situation a light-induced conical intersection (LICI) can be formed (J. Phys. B: At. Mol. Opt. Phys., 2008, 41, 221001). To measure the spectra a weak field is used whose propagation direction is perpendicular to the direction of the dressing field but has identical polarization direction. Although LICIs are present in our model, the simulations demonstrate a very robust absorption spectrum, which is insensitive to the intensity and the wavelength of the dressing field and which does not reflect clear signatures of light-induced nonadiabatic phenomena related to the strong mixing between the electronic, vibrational, rotational and translational motions. However, by widening artificially the very narrow translational energy level gaps, the fingerprint of the LICI appears to some extent in the spectrum.

Rotational-vibrational resonance states.

Resonance states are characterized by an energy that is above the lowest dissociation threshold of the potential energy hypersurface of the system and thus resonances have finite lifetimes. All molecules possess a large number of long- and short-lived resonance (quasibound) states. A considerable number of rotational-vibrational resonance states are accessible not only via quantum-chemical computations but also by spectroscopic and scattering experiments. In a number of chemical applications, most prominently in spectroscopy and reaction dynamics, consideration of rotational-vibrational resonance states is becoming more and more common. There are different first-principles techniques to compute and rationalize rotational-vibrational resonance states: one can perform scattering calculations or one can arrive at rovibrational resonances using variational or variational-like techniques based on methods developed for determining bound eigenstates. The latter approaches can be based either on the Hermitian (L2, square integrable) or non-Hermitian (non-L2) formalisms of quantum mechanics. This Perspective reviews the basic concepts related to and the relevance of shape and Feshbach-type rotational-vibrational resonance states, discusses theoretical methods and computational tools allowing their efficient determination, and shows numerical examples from the authors' previous studies on the identification and characterization of rotational-vibrational resonances of polyatomic molecular systems.

Spectroscopic signatures of HHe2+ and HHe3.

Using two different action spectroscopic techniques, a high-resolution quantum cascade laser operating around 1300 cm-1 and a cryogenic ion trap machine, the proton shuttle motion of the cations HHe2+ and HHe3+ has been probed at a nominal temperature of 4 K. For HHe3+, the loosely bound character of this complex allowed predissociation spectroscopy to be used, and the observed broad features point to a lifetime of a few ps in the vibrationally excited state. For He-H+-He, a fundamental linear molecule consisting of only three nuclei and four electrons, the method of laser-induced inhibition of complex growth (LIICG) enabled the measurement of three accurate rovibrational transitions, pinning down its molecular parameters for the first time.

Complex rovibrational dynamics of the Ar·NO+ complex.

Rotational-vibrational states of the Ar·NO+ cationic complex are computed, below, above, and well above the complex's first dissociation energy, using variational nuclear motion and close-coupling scattering computations. The HSLH potential energy surface used in this study (J. Chem. Phys., 2011, 135, 044312) is characterized by a first dissociation energy of D0 = 887.0 cm-1 and supports 200 bound vibrational states. The bound-state vibrational energies and the corresponding wave functions allow the interpretation of the scarcely available experimental results about the intermonomer vibrational motion of the complex. A very large number of long-lived quasibound combination states of the three vibrational modes, exhibiting a very similar energy-level structure as that of the bound states, are found embedded in the continuum. Additional short-lived resonance states are also identified and their properties are analyzed.

A general variational approach for computing rovibrational resonances of polyatomic molecules. Application to the weakly bound H2He+ and H2⋠CO systems.

The quasi-variational quantum chemical protocol and code GENIUSH [E. Mátyus et al., J. Chem. Phys. 130, 134112 (2009) and C. Fábri et al., J. Chem. Phys. 134, 074105 (2011)] has been augmented with the complex absorbing potential (CAP) technique, yielding a method for the determination of rovibrational resonance states. Due to the effective implementation of the CAP technique within GENIUSH, the GENIUSH-CAP code is a powerful tool for the study of important dynamical features of arbitrary-sized molecular systems with arbitrary composition above their first dissociation limit. The GENIUSH-CAP code has been tested and validated on the H2He+ cation: the computed resonance energies and lifetimes are compared to those obtained with a previously developed triatomic rovibrational resonance-computing code, D2FOPI-CCS [T. Szidarovszky and A. G. Császár Mol. Phys. 111, 2131 (2013)], utilizing the complex coordinate scaling method. A unique feature of the GENIUSH-CAP protocol is that it allows the simple implementation of reduced-dimensional dynamical models. To prove this, resonance energies and lifetimes of the H2⋅CO van der Waals complex have been computed utilizing a four-dimensional model (freezing the two monomer stretches), and a related potential energy surface, of the complex.

Toward accurate thermochemistry of the (24)MgH, (25)MgH, and (26)MgH molecules at elevated temperatures: corrections due to unbound states.

The total partition functions QT and their first two moments Q(')T and Q(″)T, together with the isobaric heat capacities CpT, are computed a priori for three major MgH isotopologues on the temperature range of T = 100-3000 K using the recent highly accurate potential energy curve, spin-rotation, and non-adiabatic correction functions of Henderson et al. [J. Phys. Chem. A 117, 13373 (2013)]. Nuclear motion computations are carried out on the ground electronic state to determine the (ro)vibrational energy levels and the scattering phase shifts. The effect of resonance states is found to be significant above about 1000 K and it increases with temperature. Even very short-lived states, due to their relatively large number, have significant contributions to QT at elevated temperatures. The contribution of scattering states is around one fourth of that of resonance states but opposite in sign. Uncertainty estimates are given for the possible error sources, suggesting that all computed thermochemical properties have an accuracy better than 0.005% up to 1200 K. Between 1200 and 2500 K, the uncertainties can rise to around 0.1%, while between 2500 K and 3000 K, a further increase to 0.5% might be observed for Q(″)T and CpT, principally due to the neglect of excited electronic states. The accurate thermochemical data determined are presented in the supplementary material for the three isotopologues of (24)MgH, (25)MgH, and (26)MgH at 1 K increments. These data, which differ significantly from older standard data, should prove useful for astronomical models incorporating thermodynamic properties of these species.

Grid-based empirical improvement of molecular potential energy surfaces.

A grid-based method designed to refine adiabatic potential energy surfaces (PES) of molecules via minimizing a suitable objective function is described. The objective function contains deviations from the reference (experimental) (ro)vibrational energy levels and is based on PES correction values determined at the grid points within a discrete-variable-representation nuclear-motion algorithm and first-order perturbation theory (PT). The proposed PES refinement technique is tested on the ground electronic state of the MgH molecule. The large number of numerical test results obtained suggest the following: (1) first-order PT is able to yield accurate correction values at the grid points representing the PES, and for practical cases there seems to be no need to go to higher orders of PT; (2) with the number of grid points greatly exceeding the number of experimental energy levels included in the refinement procedure, terms additional to the "obs-calc" term, including numerical first and second derivatives of the correction surface, are necessary in the objective function to arrive at a physically meaningful, "smooth" correction surface; (3) for a given J rotational quantum number, the corrected PES is able to reproduce experimental (ro)vibrational energies to within tenths of cm(-1) if they are included in the refinement or interpolated between states that are involved in the optimization, whereas extrapolated states tend to have somewhat larger remaining discrepancies; (4) the PES refined only for the J = 0 states introduces a minor systematic error for J > 0 states, with discrepancies growing with J; (5) when the number of experimental energies included in the refinement greatly exceeds the number of grid points upon which the PES is optimized, the systematic error of treating states with different J rotational quantum numbers can be reduced and an impressive average accuracy can be achieved for all rovibrational states; and (6) in the case of quasibound (also known as resonance) rovibrational states, energies can be computed to accuracies similar to those of the bound states and excellent lifetimes (widths) can also be determined. Changes in thermochemical functions upon inclusion of quasibound states during direct summation is discussed.

Modelling non-adiabatic effects in H3⁺: solution of the rovibrational Schrödinger equation with motion-dependent masses and mass surfaces.

Introducing different rotational and vibrational masses in the nuclear-motion Hamiltonian is a simple phenomenological way to model rovibrational non-adiabaticity. It is shown on the example of the molecular ion H3(+), for which a global adiabatic potential energy surface accurate to better than 0.1 cm(-1) exists [M. Pavanello, L. Adamowicz, A. Alijah, N. F. Zobov, I. I. Mizus, O. L. Polyansky, J. Tennyson, T. Szidarovszky, A. G. Császár, M. Berg et al., Phys. Rev. Lett. 108, 023002 (2012)], that the motion-dependent mass concept yields much more accurate rovibrational energy levels but, unusually, the results are dependent upon the choice of the embedding of the molecule-fixed frame. Correct degeneracies and an improved agreement with experimental data are obtained if an Eckart embedding corresponding to a reference structure of D(3h) point-group symmetry is employed. The vibrational mass of the proton in H3(+) is optimized by minimizing the root-mean-square (rms) deviation between the computed and recent high-accuracy experimental transitions. The best vibrational mass obtained is larger than the nuclear mass of the proton by approximately one third of an electron mass, m(opt,p)((v))=m(nuc,p)+0.31224m(e). This optimized vibrational mass, along with a nuclear rotational mass, reduces the rms deviation of the experimental and computed rovibrational transitions by an order of magnitude. Finally, it is shown that an extension of the algorithm allowing the use of motion-dependent masses can deal with coordinate-dependent mass surfaces in the rovibrational Hamiltonian, as well.

MARVEL analysis of the rotational-vibrational states of the molecular ions H2D+ and D2H+.

Critically evaluated rotational-vibrational line positions and energy levels, with associated critically reviewed labels and uncertainties, are reported for two deuterated isotopologues of the H3(+) molecular ion: H2D(+) and D2H(+). The procedure MARVEL, standing for Measured Active Rotational-Vibrational Energy Levels, is used to determine the validated levels and lines and their self-consistent uncertainties based on the experimentally available information. The spectral ranges covered for the isotopologues H2D(+) and D2H(+) are 5.2-7105.5 and 23.0-6581.1 cm(-1), respectively. The MARVEL energy levels of the ortho and para forms of the ions are checked against ones determined from accurate variational nuclear motion computations employing the best available adiabatic ab initio potential energy surfaces of these isotopologues. The number of critically evaluated, validated and recommended experimental (levels, lines) are (109, 185) and (104, 136) for H2D(+) and D2H(+), respectively. The lists of assigned MARVEL lines and levels and variational levels obtained for H2D(+) and D2H(+) as part of this study are deposited in the ESI to this paper.

Precision measurements and computations of transition energies in rotationally cold triatomic hydrogen ions up to the midvisible spectral range.

First-principles computations and experimental measurements of transition energies are carried out for vibrational overtone lines of the triatomic hydrogen ion H(3)(+) corresponding to floppy vibrations high above the barrier to linearity. Action spectroscopy is improved to detect extremely weak visible-light spectral lines on cold trapped H(3)(+) ions. A highly accurate potential surface is obtained from variational calculations using explicitly correlated Gaussian wave function expansions. After nonadiabatic corrections, the floppy H(3)(+) vibrational spectrum is reproduced at the 0.1 cm(-1) level up to 16600 cm(-1).

The fourth age of quantum chemistry: molecules in motion.

Developments during the last two decades in nuclear motion theory made it possible to obtain variational solutions to the time-independent, nuclear-motion Schrödinger equation of polyatomic systems as "exact" as the potential energy surface (PES) is. Nuclear motion theory thus reached a level whereby this branch of quantum chemistry started to catch up with the well developed and widely applied other branch, electronic structure theory. It seems to be fair to declare that we are now in the fourth age of quantum chemistry, where the first three ages are principally defined by developments in electronic structure techniques (G. Richards, Nature, 1979, 278, 507). In the fourth age we are able to incorporate into our quantum chemical treatment the motion of nuclei in an exact fashion and, for example, go beyond equilibrium molecular properties and compute accurate, temperature-dependent, effective properties, thus closing the gap between measurements and electronic structure computations. In this Perspective three fundamental algorithms for the variational solution of the time-independent nuclear-motion Schrödinger equation employing exact kinetic energy operators are presented: one based on tailor-made Hamiltonians, one on the Eckart-Watson Hamiltonian, and one on a general internal-coordinate Hamiltonian. It is argued that the most useful and most widely applicable procedure is the third one, based on a Hamiltonian containing a kinetic energy operator written in terms of internal coordinates and an arbitrary embedding of the body-fixed frame of the molecule. This Hamiltonian makes it feasible to treat the nuclear motions of arbitrary quantum systems, irrespective of whether they exhibit a single well-defined minimum or not, and of arbitrary reduced-dimensional models. As a result, molecular spectroscopy, an important field for the application of nuclear motion theory, has almost black-box-type tools at its disposal. Variational nuclear motion computations, based on an exact kinetic energy operator and an arbitrary PES, can now be performed for about 9 active vibrational degrees of freedom relatively straightforwardly. Simulations of high-resolution spectra allow the understanding of complete rotational-vibrational spectra up to and beyond the first dissociation limits. Variational results obtained for H(2)O, H, NH(3), CH(4), and H(2)CCO are used to demonstrate the power of the variational techniques for the description of vibrational and rotational excitations. Some qualitative features of the results are also discussed.

The role of axis embedding on rigid rotor decomposition analysis of variational rovibrational wave functions.

Approximate rotational characterization of variational rovibrational wave functions via the rigid rotor decomposition (RRD) protocol is developed for Hamiltonians based on arbitrary sets of internal coordinates and axis embeddings. An efficient and general procedure is given that allows employing the Eckart embedding with arbitrary polyatomic Hamiltonians through a fully numerical approach. RRD tables formed by projecting rotational-vibrational wave functions into products of rigid-rotor basis functions and previously determined vibrational eigenstates yield rigid-rotor labels for rovibrational eigenstates by selecting the largest overlap. Embedding-dependent RRD analyses are performed, up to high energies and rotational excitations, for the H(2) (16)O isotopologue of the water molecule. Irrespective of the embedding chosen, the RRD procedure proves effective in providing unambiguous rotational assignments at low energies and J values. Rotational labeling of rovibrational states of H(2) (16)O proves to be increasingly difficult beyond about 10,000 cm(-1), close to the barrier to linearity of the water molecule. For medium energies and excitations the Eckart embedding yields the largest RRD coefficients, thus providing the largest number of unambiguous rotational labels.

Calibration-quality adiabatic potential energy surfaces for H3(+) and its isotopologues.

Calibration-quality ab initio adiabatic potential energy surfaces (PES) have been determined for all isotopologues of the molecular ion H(3)(+). The underlying Born-Oppenheimer electronic structure computations used optimized explicitly correlated shifted Gaussian functions. The surfaces include diagonal Born-Oppenheimer corrections computed from the accurate electronic wave functions. A fit to the 41,655 ab initio points is presented which gives a standard deviation better than 0.1 cm(-1) when restricted to the points up to 6000 cm(-1) above the first dissociation asymptote. Nuclear motion calculations utilizing this PES, called GLH3P, and an exact kinetic energy operator given in orthogonal internal coordinates are presented. The ro-vibrational transition frequencies for H(3)(+), H(2)D(+), and HD(2)(+) are compared with high resolution measurements. The most sophisticated and complete procedure employed to compute ro-vibrational energy levels, which makes explicit allowance for the inclusion of non-adiabatic effects, reproduces all the known ro-vibrational levels of the H(3)(+) isotopologues considered to better than 0.2 cm(-1). This represents a significant (order-of-magnitude) improvement compared to previous studies of transitions in the visible. Careful treatment of linear geometries is important for high frequency transitions and leads to new assignments for some of the previously observed lines. Prospects for further investigations of non-adiabatic effects in the H(3)(+) isotopologues are discussed. In short, the paper presents (a) an extremely accurate global potential energy surface of H(3)(+) resulting from high accuracy ab initio computations and global fit, (b) very accurate nuclear motion calculations of all available experimental line data up to 16,000 cm(-1), and (c) results suggest that we can predict accurately the lines of H(3)(+) towards dissociation and thus facilitate their experimental observation.

Achieving high molecular alignment and orientation for CH[Formula: see text]F through manipulation of rotational states with varying optical and THz laser pulse parameters.

Increasing interest in the fields of high-harmonics generation, laser-induced chemical reactions, and molecular imaging of gaseous targets demands high molecular "alignment" and "orientation" (A&O). In this work, we examine the critical role of different pulse parameters on the field-free A&O dynamics of the CH[Formula: see text]F molecule, and identify experimentally feasible optical and THz range laser parameters that ensure maximal A&O for such molecules. Herein, apart from rotational temperature, we investigate effects of varying pulse parameters such as, pulse duration, intensity, frequency, and carrier envelop phase (CEP). By analyzing the interplay between laser pulse parameters and the resulting rotational population distribution, the origin of specific A&O dynamics was addressed. We could identify two qualitatively different A&O behaviors and revealed their connection with the pulse parameters and the population of excited rotational states. We report here the highest alignment of [Formula: see text] and orientation of [Formula: see text] for CH[Formula: see text]F molecule at 2 K using a single pulse. Our study should be useful to understand different aspects of laser-induced unidirectional rotation in heteronuclear molecules, and in understanding routes to tune/enhance A&O in laboratory conditions for advanced applications.

Gas-phase and Ar-matrix SQM scaling factors for various DFT functionals with basis sets including polarization and diffuse functions.

Scaling factors for Pulay's scaled quantum mechanical (SQM) scheme have been determined for four different widely used DFT functionals (PBE, B3LYP, B3PW91, and M06-2X) and for two basis sets (6-31++G** and aug-cc-pVTZ) by fitting computed results to 347 fundamental experimental vibrational frequencies of 33 molecules. Measurements in the gas phase and in solid argon matrices were used independently in the fitting procedure in order to provide a simple method of estimating matrix shifts. The accuracy of the new scaling factors is demonstrated on test molecules including hydrogen-bonded systems and molecules containing chlorine and sulfur atoms.

On the efficiency of treating singularities in triatomic variational vibrational computations. The vibrational states of H(+)3 up to dissociation.

Several techniques of varying efficiency are investigated, which treat all singularities present in the triatomic vibrational kinetic energy operator given in orthogonal internal coordinates of the two distances-one angle type. The strategies are based on the use of a direct-product basis built from one-dimensional discrete variable representation (DVR) bases corresponding to the two distances and orthogonal Legendre polynomials, or the corresponding Legendre-DVR basis, corresponding to the angle. The use of Legendre functions ensures the efficient treatment of the angular singularity. Matrix elements of the singular radial operators are calculated employing DVRs using the quadrature approximation as well as special DVRs satisfying the boundary conditions and thus allowing for the use of exact DVR expressions. Potential optimized (PO) radial DVRs, based on one-dimensional Hamiltonians with potentials obtained by fixing or relaxing the two non-active coordinates, are also studied. The numerical calculations employed Hermite-DVR, spherical-oscillator-DVR, and Bessel-DVR bases as the primitive radial functions. A new analytical formula is given for the determination of the matrix elements of the singular radial operator using the Bessel-DVR basis. The usually claimed failure of the quadrature approximation in certain singular integrals is revisited in one and three dimensions. It is shown that as long as no potential optimization is carried out the quadrature approximation works almost as well as the exact DVR expressions. If wave functions with finite amplitude at the boundary are to be computed, the basis sets need to meet the required boundary conditions. The present numerical results also confirm that PO-DVRs should be constructed employing relaxed potentials and PO-DVRs can be useful for optimizing quadrature points for calculations applying large coordinate intervals and describing large-amplitude motions. The utility and efficiency of the different algorithms is demonstrated by the computation of converged near-dissociation vibrational energy levels for the H molecular ion.

Assigning quantum labels to variationally computed rotational-vibrational eigenstates of polyatomic molecules.

A procedure is investigated for assigning physically transparent, approximate vibrational and rotational quantum labels to variationally computed eigenstates. Pure vibrational wave functions are analyzed by means of normal-mode decomposition (NMD) tables constructed from overlap integrals with respect to separable harmonic oscillator basis functions. Complementary rotational labels J(K(a)K(c)) are determined from rigid-rotor decomposition (RRD) tables formed by projecting rotational-vibrational wave functions (J not equal 0) onto products of symmetrized rigid-rotor basis functions and previously computed (J=0) vibrational eigenstates. Variational results for H(2)O, HNCO, trans-HCOD, NCCO, and H(2)CCO are presented to demonstrate the NMD and RRD schemes. The NMD analysis highlights several resonances at low energies that cause strong mixing and cloud the assignment of fundamental vibrations, even in such simple molecules. As the vibrational energy increases, the NMD scheme documents and quantifies the breakdown of the normal-mode model. The RRD procedure proves effective in providing unambiguous rotational assignments for the chosen test molecules up to moderate J values.