Coupling polyatomic molecules to lossy nanocavities: Lindblad vs Schrödinger description.

The use of cavities to impact molecular structure and dynamics has become popular. As cavities, in particular plasmonic nanocavities, are lossy and the lifetime of their modes can be very short, their lossy nature must be incorporated into the calculations. The Lindblad master equation is commonly considered an appropriate tool to describe this lossy nature. This approach requires the dynamics of the density operator and is thus substantially more costly than approaches employing the Schrödinger equation for the quantum wave function when several or many nuclear degrees of freedom are involved. In this work, we compare numerically the Lindblad and Schrödinger descriptions discussed in the literature for a molecular example where the cavity is pumped by a laser. The laser and cavity properties are varied over a range of parameters. It is found that the Schrödinger description adequately describes the dynamics of the polaritons and emission signal as long as the laser intensity is moderate and the pump time is not much longer than the lifetime of the cavity mode. Otherwise, it is demonstrated that the Schrödinger description gradually fails. We also show that the failure of the Schrödinger description can often be remedied by renormalizing the wave function at every step of time propagation. The results are discussed and analyzed.

All paths lead to hubs in the spectroscopic networks of water isotopologues H216O and H218O.

Network theory has fundamentally transformed our comprehension of complex systems, catalyzing significant advances across various domains of science and technology. In spectroscopic networks, hubs are the quantum states involved in the largest number of transitions. Here, utilizing network paths probed via precision metrology, absolute energies have been deduced, with at least 10-digit accuracy, for almost 200 hubs in the experimental spectroscopic networks of H216O and H218O. These hubs, lying on the ground vibrational states of both species and the bending fundamental of H216O, are involved in tens of thousands of observed transitions. Relying on the same hubs and other states, benchmark-quality line lists have been assembled, which supersede and improve, by three orders of magnitude, the accuracy of the massive amount of data reported in hundreds of papers dealing with Doppler-limited spectroscopy. Due to the omnipresence of water, these ultraprecise line lists could be applied to calibrate high-resolution spectra and serve ongoing and upcoming space missions.

MARVEL analysis of high-resolution rovibrational spectra of 13 C 16 O 2 .

A set of empirical rovibrational energy levels, obtained through the MARVEL (measured active rotational-vibrational energy levels) procedure, is presented for the 13 C 16 O 2 isotopologue of carbon dioxide. This procedure begins with the collection and analysis of experimental rovibrational transitions from the literature, allowing for a comprehensive review of the literature on the high-resolution spectroscopy of 13 C 16 O 2 , which is also presented. A total of 60 sources out of more than 750 checked provided 14,101 uniquely measured and assigned rovibrational transitions in the wavenumber range of 579-13,735 cm - 1 . This is followed by a weighted least-squares refinement yielding the energy levels of the states involved in the measured transitions. Altogether 6318 empirical rovibrational energies have been determined for 13 C 16 O 2 . Finally, estimates have been given for the uncertainties of the empirical energies, based on the experimental uncertainties of the transitions. The detailed analysis of the lines and the spectroscopic network built from them, as well as the uncertainty estimates, all serve to pinpoint possible errors in the experimental data, such as typos, misassignment of quantum numbers, and misidentifications. Errors found in the literature data were corrected before including them in the final MARVEL dataset and analysis.

Verification labels for rovibronic quantum-state energy uncertainties.

Transition wavenumbers contained in line-by-line rovibronic databases can be compromised by errors of various nature. When left undetected, these errors may result in incorrect quantum-state energies, potentially compromising a large number of derived spectroscopic data. Spectroscopic networks treat the complete set of line-by-line spectroscopic data as a large graph, and through a least-squares refinement the measured line positions are converted into empirical quantum-state energies. Spectroscopic networks also offer a highly useful framework to develop mathematical tools helping to identify possible errors and conflicts within the dataset. For example, wavenumber errors can be detected by checking for violations of the law of energy conservation. This paper describes a new graph-theory tool, which results in so-called verification labels for the quantum states. Verification labels help to express the vulnerability of a calculated empirical energy value and its uncertainty against possible wavenumber errors, providing complementary information to simple statistical uncertainties.

Temperature-Dependent Line-Broadening Effects in CO2 Caused by Ar.

This computational study of line-broadening effects is based on an accurate, analytical representation of the intermonomer potential energy surface (PES) of the CO2 ⋅ Ar van der Waals (vdW) complex. The PES is employed to compute collisional broadening coefficients for rovibrational lines of CO2 perturbed by Ar. The semiclassical computations are performed using the modified Robert-Bonamy approach, including real and imaginary terms, and the exact trajectory model. The lines investigated are in the 10001←00011, 01101←00001, 00011←00001, and 00031←00001 vibrational bands and the computations are repeated at multiple temperatures. The computed results are in good agreement with the available experimental values, validating both the intermonomer PES developed and the methodology used. For lines in the 01101←00001 band of CO2 , temperature-dependent Ar-broadening coefficients are reported for the first time. The parameters presented should prove useful, among other applications, for the accurate experimental determination of CO2 and Ar abundances in planetary atmospheres.

Quantum Chemical Investigation of the Cold Water Dimer Spectrum in the First OH-Stretching Overtone Region Provides a New Interpretation.

Intramolecular vibrational transition wavenumbers and intensities were calculated in the fundamental HOH-bending, fundamental OH-stretching, first OH-stretching-HOH-bending combination, and first OH-stretching overtone (ΔvOH = 2) regions of the water dimer's spectrum. Furthermore, the rotational-vibrational spectrum was calculated in the ΔvOH = 2 region at 10 K, corresponding to the temperature of the existing jet-expansion experiments. The calculated spectrum was obtained by combining results from a full-dimensional (12D) vibrational and a reduced-dimensional vibrational-rotational-tunneling model. The ΔvOH = 2 spectral region is rich in features due to contributions from multiple vibrational-rotational-tunneling sub-bands. Origins of the experimental vibrational bands depend on the assignment of the observed sub-bands. Based on our calculations, we assign the observed sub-bands, and our reassignment leads to new values for the vibrational band origins of the free donor and antisymmetric acceptor OH-stretching first overtones of ∼7227 and ∼7238 cm-1, respectively. The observed bands with origins at 7192.34 and ∼7366 cm-1 are assigned to the symmetric acceptor OH-stretching first overtone and the OH-stretching combination of the donor, respectively.

Unusual Dynamics and Vibrational Fingerprints of van der Waals Dimers Formed by Linear Molecules and Rare-Gas Atoms.

Detailed structural, dynamical, and vibrational analyses have been performed for systems composed of linear triatomic molecules solvated by a single rare-gas atom, He, Ne, or Ar. Among the chromophores of these van der Waals (vdW) dimers, there are four neutral molecules (CO2, CS2, N2O, and OCS) and six molecular cations (HHe2+, HNe2+, HAr2+, HHeNe+, HHeAr+, and HNeAr+), both of apolar and polar nature. Following the exploration of bonding preferences, high-level four-dimensional (4D) potential energy surfaces (PESs) have been developed for 24 vdW dimers, keeping the two intramonomer bond lengths fixed. For these 24 complexes, over 1500 bound vibrational states have been obtained via quasi-variational nuclear-motion computations, employing exact kinetic-energy operators together with the accurate 4D PESs and their 2D/3D cuts. The reduced-dimensional (2D to 4D) dimer models have been compared with full-dimensional (6D) ones in the cases of the neutral CO2·Ar and charged HHe2+·He dimers, corroborating the high accuracy of the 2D to 4D vibrational energies. The reduced-dimensional models suggest that (a) while the equilibrium structures are T-shaped and planar, the effective ground-state structures are nonplanar, (b) certain bound states belong to collinear molecular structures, even when they are not minima, (c) the vdW vibrations are heavily mixed and many states have amplitudes corresponding to both the T-shaped and collinear structures, (d) there are a few dimers, for which even some of the vdW fundamentals lie above the first dissociation limit, and (e) the vdW vibrations are almost fully decoupled from the intramonomer bending motion.

On the 12C2H2 near-infrared spectrum: absolute transition frequencies and an improved spectroscopic network at the kHz accuracy level.

Lamb dips of twenty lines in the P, Q, and R branches of the ν1 + ν3 + ν41 vibrational band of 12C2H2, in the spectral window of 7125-7230 cm-1, have been measured using an upgraded comb-calibrated frequency-stabilized cavity ring-down spectrometer, designed for extensive sub-Doppler measurements. Due to the large number of carefully executed Lamb-dip experiments, and to the extrapolation of absolute frequencies to zero pressure in each case, the combined average uncertainty of the measured line-center positions is 15 kHz (5 × 10-7 cm-1) with a 2-σ confidence level. Selection of the twenty lines was based on the theory of spectroscopic networks (SN), ensuring that a large number of transitions, measured previously by precision-spectroscopy investigations, could be connected to the para and ortho principal components of the SN of 12C2H2. The assembled SN contains 331 highly precise transitions, 119 and 121 of which are in the ortho and para principal components, respectively, while the rest remain in floating components. The para- and ortho-12C2H2 energy-level lists, determined during the present study, contain 82 and 80 entries, respectively, with an accuracy similar to that of the lines. Based on the newly assembled lists of para- and ortho-12C2H2 empirical energy levels, a line list, called TenkHz, has been generated. The TenkHz line list contains 282 entries in the spectral range of 5898.97-7258.87 cm-1; thus far, only 149 of them have been measured directly via precision spectroscopy. The TenkHz line list includes 35 intense lines that are missing in the HITRAN2020 database.

Quantum Nuclear Delocalization and its Rovibrational Fingerprints.

Quantum mechanics dictates that nuclei must undergo some delocalization. In this work, emergence of quantum nuclear delocalization and its rovibrational fingerprints are discussed for the case of the van der Waals complex HHe 3 + ${{\rm{HHe}}_3^ + }$ . The equilibrium structure of HHe 3 + ${{\rm{HHe}}_3^ + }$ is planar and T-shaped, one He atom solvating the quasi-linear He-H+ -He core. The dynamical structure of HHe 3 + ${{\rm{HHe}}_3^ + }$ , in all of its bound states, is fundamentally different. As revealed by spatial distribution functions and nuclear densities, during the vibrations of the molecule the solvating He is not restricted to be in the plane defined by the instantaneously bent HHe 2 + ${{\rm{HHe}}_2^ + }$ chomophore, but freely orbits the central proton, forming a three-dimensional torus around the HHe 2 + ${{\rm{HHe}}_2^ + }$ chromophore. This quantum delocalization is observed for all vibrational states, the type of vibrational excitation being reflected in the topology of the nodal surfaces in the nuclear densities, showing, for example, that intramolecular bending involves excitation along the circumference of the torus.

Quantum-Chemical and Quantum-Graph Models of the Dynamical Structure of CH5.

Experimental and computational results about the structure, dynamics, and rovibrational spectra of protonated methane have challenged a considerable number of traditional chemical concepts. Hereby theoretical and computational results are provided about the dynamical structure of CH5+. It is shown that the ground vibrational state investigated thus far by computations, forbidden by nuclear-spin statistics, has a structure similar to the first allowed vibrational state and, in fact, the structures of all vibrational states significantly below 200 cm-1 are highly similar. Spatial delocalization of the nuclei, determined by nuclear densities computed from accurate variational vibrational wave functions, turns out to be limited when viewed in the body-fixed frame, confirming that the effective structure of CH5+ is well described as a CH3+ tripod with a H2 unit on top of it. The interesting and unusual qualitative aspects of the sophisticated state-dependent variational results receive full explanation via simple quantum-graph models.

Parity-pair-mixing effects in nonlinear spectroscopy of HDO.

A non-linear spectroscopic study of the HDO molecule is performed in the wavelength range of 1.36-1.42 µm using noise-immune cavity-enhanced optical-heterodyne molecular spectroscopy (NICE-OHMS). More than 100 rovibrational Lamb dips are recorded, with an experimental precision of 2-20 kHz, related to the first overtone of the O-H stretch fundamental of HD16O and HD18O. Significant perturbations, including distortions, shifts, and splittings, have been observed for a number of Lamb dips. These spectral perturbations are traced back to an AC-Stark effect, arising due to the strong laser field applied in all saturation-spectroscopy experiments. The AC-Stark effect mixes parity pairs, that is pairs of rovibrational states whose assignment differs solely in the Kc quantum number, where Kc is part of the standard J K a,K c asymmetric-top rotational label. Parity-pair mixing seems to be especially large for parity pairs with Ka ≥ 3, whereby their energy splittings become as small as a few MHz, resulting in multi-component asymmetric Lamb-dip profiles of gradually increasing complexity. These complex profiles often include crossover resonances. This effect is well known in saturation spectroscopy, but has not been reported in combination with parity-pair mixing. Parity-pair mixing is not seen in H2 16O and H2 18O, because their parity pairs correspond to ortho and para nuclear-spin isomers, whose interaction is prohibited. Despite the frequency shifts observed for HD16O and HD18O, the absolute accuracy of the detected transitions still exceeds that achievable by Doppler-limited techniques.

Analysis of measured high-resolution doublet rovibronic spectra and related line lists of 12CH and 16OH.

Detailed understanding of the energy-level structure of the quantum states as well as of the rovibronic spectra of the ethylidyne (CH) and the hydroxyl (OH) radicals is mandatory for a multitude of modelling efforts within multiple chemical, combustion, astrophysical, and atmospheric environments. Accurate empirical rovibronic energy levels, with associated uncertainties, are reported for the low-lying doublet electronic states of 12CH and 16OH, using the Measured Active Rotational-Vibrational Energy Levels (MARVEL) algorithm. For 12CH, a total of 1521 empirical energy levels are determined in the primary spectroscopic network (SN) of the radical, corresponding to the following seven electronic states: X 2Π, A 2Δ, B 2Σ-, C2 Σ+, D 2Π, E 2Σ+, and F 2Σ+. The energy levels are derived from 6348 experimentally measured and validated transitions, collected from 29 sources. For 16OH, the lowest four doublet electronic states, X 2Π, A 2Σ+, B 2Σ+, and C 2Σ+, are considered, and a careful analysis and validation of 15 938 rovibronic transitions, collected from 45 sources, results in 1624 empirical rovibronic energy levels. The large set of spectroscopic data presented should facilitate the refinement of line lists for the 12CH and 16OH radicals. For both molecules hyperfine-resolved experimental transitions have also been considered, forming SNs independent from the primary SNs.

Structure, energetics, and spectroscopy of the chromophores of HHe+n, H2He+n, and He+n clusters and their deuterated isotopologues.

The linear molecular ions H2He+, HHe+2, and He+3 are the central units (chromophores) of certain He-solvated complexes of the H2He+n, HHe+n, and He+n families, respectively. These are complexes which do exist, according to mass-spectrometry studies, up to very high n values. Apparently, for some of the H2He+n and He+n complexes, the linear symmetric tetratomic H2He+2 and the diatomic He+2 cations, respectively, may also be the central units. In this study, definitive structures, relative energies, zero-point vibrational energies, and (an)harmonic vibrational fundamentals, and, in some cases, overtones and combination bands, are established mostly for the triatomic chromophores. The study is also extended to the deuterated isotopologues D2He+, DHe+2, and D2He+2. To facilitate and improve the electronic-structure computations performed, new atom-centered, fixed-exponent, Gaussian-type basis sets called MAX, with X = T(3), Q(4), P(5), and H(6), are designed for the H and He atoms. The focal-point-analysis (FPA) technique is employed to determine definitive relative energies with tight uncertainties for reactions involving the molecular ions. The FPA results determined include the 0 K proton and deuteron affinities of the 4He atom, 14 875(9) cm-1 [177.95(11) kJ mol-1] and 15 229(8) cm-1 [182.18(10) kJ mol-1], respectively, the dissociation energies of the He+2 → He+ + He, HHe+2 → HHe+ + He, and He+3 → He+2 + He reactions, 19 099(13) cm-1 [228.48(16) kJ mol-1], 3948(7) cm-1 [47.23(8) kJ mol-1], and 1401(12) cm-1 [16.76(14) kJ mol-1], respectively, the dissociation energy of the DHe+2 → DHe+ + He reaction, 4033(6) cm-1 [48.25(7) kJ mol-1], the isomerization energy between the two linear isomers of the [H, He, He]+ system, 3828(40) cm-1 [45.79(48) kJ mol-1], and the dissociation energies of the H2He+ → H+2 + He and the H2He+2 → H2He+ + He reactions, 1789(4) cm-1 [21.40(5) kJ mol-1] and 435(6) cm-1 [5.20(7) kJ mol-1], respectively. The FPA estimates of the first dissociation energy of D2He+ and D2He+2 are 1986(4) cm-1 [23.76(5) kJ mol-1] and 474(5) cm-1 [5.67(6) kJ mol-1], respectively. Determining the vibrational fundamentals of the triatomic chromophores with second-order vibrational perturbation theory (VPT2) and vibrational configuration interaction (VCI) techniques, both built around the Eckart-Watson Hamiltonian, proved unusually challenging. For the species studied, VPT2 has difficulties yielding dependable results, in some cases even for the fundamentals of the H-containing molecular cations, while carefully executed VCI computations yield considerably improved spectroscopic results. In a few cases unusually large anharmonic corrections to the fundamentals, on the order of 15% of the harmonic value, have been observed.

Reduced-dimensional vibrational models of the water dimer.

A model based on the finite-basis representation of a vibrational Hamiltonian expressed in internal coordinates is developed. The model relies on a many-mode, low-order expansion of both the kinetic energy operator and the potential energy surface (PES). Polyad truncations and energy ceilings are used to control the size of the vibrational basis to facilitate accurate computations of the OH stretch and HOH bend intramolecular transitions of the water dimer (H2 16O)2. Advantages and potential pitfalls of the applied approximations are highlighted. The importance of choices related to the treatment of the kinetic energy operator in reduced-dimensional calculations and the accuracy of different water dimer PESs are discussed. A range of different reduced-dimensional computations are performed to investigate the wavenumber shifts in the intramolecular transitions caused by the coupling between the intra- and intermolecular modes. With the use of symmetry, full 12-dimensional vibrational energy levels of the water dimer are calculated, predicting accurately the experimentally observed intramolecular fundamentals. It is found that one can also predict accurate intramolecular transition wavenumbers for the water dimer by combining a set of computationally inexpensive reduced-dimensional calculations, thereby guiding future effective-Hamiltonian treatments.

Normal-Mode Vibrational Analysis of Weakly Bound Oligomers at Constrained Stationary Points of Arbitrary Order.

Following the full realization of the importance of noncovalent interactions, finding and characterizing stationary points (SP), of various order, for weakly bound oligomers have become important tasks for computational chemists. An efficient algorithm and an associated computer code, called oligoCGO, are described, facilitating constrained geometry optimization of oligomers of arbitrary structure and complexity and normal-mode vibrational analysis at nonstationary geometries. To minimize the adverse effects of nonzero forces on harmonic vibrational analyses at constrained stationary points (cSP), two residual gradient correction (RGC) schemes are proposed. RGC1, for which a rigorous justification is given, is based on ignoring the remaining forces in internal-coordinate space. RGC2 modifies the geometry of the cSP in a single Newton step and recalculates the Cartesian Hessian at this updated geometry. As demonstrated by 10 examples related to the water-water, water-methane, and methane-methane dimers as well as the methane trimer, without RGC the harmonic analysis of cSPs may result in even qualitatively incorrect results when compared to reference values obtained at the nearby unconstrained SPs (uSP). Both RGC protocols work exceedingly well, and the corrected harmonic wavenumbers of the cSPs are very close to their uSP counterparts.

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.

The rovibrational Aharonov-Bohm effect.

Another manifestation of the Aharonov-Bohm effect is introduced to chemistry, in fact to nuclear quantum dynamics and high-resolution molecular spectroscopy. As demonstrated, the overall rotation of a symmetric-top molecule influences the dynamics of an internal vibrational motion in a way that is analogous to the presence of a solenoid carrying magnetic flux. To a good approximation, the low-energy rovibrational energy-level structure of the quasistructural molecular ion H+5 can be understood entirely in terms of this effect.

Selecting lines for spectroscopic (re)measurements to improve the accuracy of absolute energies of rovibronic quantum states.

Improving the accuracy of absolute energies associated with rovibronic quantum states of molecules requires accurate high-resolution spectroscopy measurements. Such experiments yield transition wavenumbers from which the energies can be deduced via inversion procedures. To address the problem that not all transitions contribute equally to the goal of improving the accuracy of the energies, the method of Connecting Spectroscopic Components (CSC) is introduced. Using spectroscopic networks and tools of graph theory, CSC helps to find the most useful target transitions and target wavenumber regions for (re)measurement. The sets of transitions suggested by CSC should be investigated by experimental research groups in order to select those target lines which they can actually measure based on the apparatus available to them. The worked-out examples, utilizing extensive experimental spectroscopic data on the molecules H[Formula: see text]O, [Formula: see text]S[Formula: see text]O[Formula: see text], H[Formula: see text]C[Formula: see text]O, and [Formula: see text]NH[Formula: see text], clearly prove the overall usefulness of the CSC method and provide suggestions how CSC can be used for various tasks and under different practical circumstances.

Exactly solvable 1D model explains the low-energy vibrational level structure of protonated methane.

A new one-dimensional model is proposed for the low-energy vibrational quantum dynamics of CH5+ based on the motion of an effective particle confined to a 60-vertex graph Γ60 with a single edge length parameter. Within this model, the quantum states of CH5+ are obtained in analytic form and are related to combinatorial properties of Γ60. The bipartite structure of Γ60 gives a simple explanation for curious symmetries observed in numerically exact variational calculations on CH5+.

Understanding the structure of complex multidimensional wave functions. A case study of excited vibrational states of ammonia.

Generalization of an earlier reduced-density-matrix-based vibrational assignment algorithm is given, applicable for systems exhibiting both large-amplitude motions, including tunneling, and degenerate vibrational modes. The algorithm developed is used to study the structure of the excited vibrational wave functions of the ammonia molecule, 14NH3. Characterization of the complex dynamics of systems with several degenerate vibrations requires reconsidering the traditional degenerate-mode description given by vibrational angular momentum quantum numbers and switching to a symmetry-based approach that directly predicts state degeneracy and uncovers relations between degenerate modes. Out of the 600 distinct vibrational eigenstates of ammonia obtained by a full-dimensional variational computation, the developed methodology allows for the assignment of about 500 with meaningful labels. This study confirms that vibrationally excited states truly have modal character recognizable up to very high energies even for the non-trivial case of ammonia, a molecule which exhibits a tunneling motion and has two two-dimensional normal modes. The modal characteristics of the excited states and the interplay of the vibrational modes can be easily visualized by the reduced-density matrices, giving an insight into the complex modal behavior directed by symmetry.