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.

MARVEL analysis of high-resolution rovibrational spectra of 16O12C18O.

Empirical rovibrational energy levels are presented for the third most abundant, asymmetric carbon dioxide isotopologue, 16O12C18O, based on a compiled dataset of experimental rovibrational transitions collected from the literature. The 52 literature sources utilized provide 19,438 measured lines with unique assignments in the wavenumber range of 2-12,676 cm-1. The MARVEL (Measured Active Rotational-Vibrational Energy Levels) protocol, which is built upon the theory of spectroscopic networks, validates the great majority of these transitions and outputs 8786 empirical rovibrational energy levels with an uncertainty estimation based on the experimental uncertainties of the transitions. Issues found in the literature data, such as misassignment of quantum numbers, typographical errors, and misidentifications, are fixed before including them in the final MARVEL dataset and analysis. Comparison of the empirical energy-level data of this study with those in the line lists CDSD-2019 and Ames-2021 shows good overall agreement, significantly better for CDSD-2019; some issues raised by these comparisons are discussed.

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.

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.

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 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.

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.

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.

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.

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.

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.

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.

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.

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.

Comment on "Wigner numbers" [J. Chem. Phys. 151, 244122 (2019)].

An alternative combinatorial expression is presented for the Wigner numbers Wm,n J that leads to a generating function and to several new identities. An extended class of Wigner numbers, Wigner numbers of the ℓth order, Wm,n J,ℓ, is also introduced.

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.

Rovibrational quantum dynamics of the vinyl radical and its deuterated isotopologues.

Rotational-vibrational states up to 3200 cm-1, beyond the highest-lying stretching fundamental, are computed variationally for the vinyl radical (VR), H2Cß[double bond, length as m-dash]CαH, and the following deuterated isotopologues of VR: CH2[double bond, length as m-dash]CD, CHD[double bond, length as m-dash]CH, and CD2[double bond, length as m-dash]CD. The height of the CαH tunneling rocking barrier of VR, partially responsible for the complex nuclear dynamics of VR and its isotopologues, is determined to be 1641 ± 25 cm-1 by the focal-point analysis approach. The definitive nuclear-motion computations performed utilize two previously published potential energy hypersurfaces and reveal interesting energy-level and tunneling patterns characterizing the internal motions of the four isotopologues. A full assignment, including symmetry labels, of the vibrational states computed for CH2[double bond, length as m-dash]CH is provided, whenever feasible, based on the analysis of wave functions and the related one- and two-mode reduced density matrices. The computed vibrational states of CH2[double bond, length as m-dash]CD and CD2[double bond, length as m-dash]CD are characterized up to slightly above the top of the barrier. Interestingly, it is the interplay of the ν6 (formally CH2 rock) and ν7 (formally CH rock) modes that determines the tunneling dynamics; thus, the description of tunneling in VR needs, as a minimum, the consideration of two in-plane bending motions at the two ends of the molecule. When feasible, the computed results are compared to their experimental counterparts as well as to previous computational results. Corrections to the placement of the ν4 and ν6 fundamentals of VR are proposed. Tunneling switching, a unique phenomenon characterizing tunneling in slightly asymmetric effective double-well potentials, is observed and discussed for CHD[double bond, length as m-dash]CH. Despite the extensive tunneling dynamics, the rotational energy-level structure of VR exhibits rigid-rotor-type behavior.

Accurate empirical rovibrational energies and transitions of H216O.

Several significant improvements are proposed to the computational molecular spectroscopy protocol MARVEL (Measured Active Rotational-Vibrational Energy Levels) facilitating the inversion of a large set of measured rovibrational transitions to energy levels. The most important algorithmic changes include the use of groups of transitions, blocked by their estimated experimental (source segment) uncertainties, an inversion and weighted least-squares refinement procedure based on sequential addition of blocks of decreasing accuracy, the introduction of spectroscopic cycles into the refinement process, automated recalibration, synchronization of the combination difference relations to reduce residual uncertainties in the resulting dataset of empirical (MARVEL) energy levels, and improved classification of the lines and energy levels based on their accuracy and dependability. The resulting protocol, through handling a large number of measurements of similar accuracy, retains, or even improves upon, the best reported uncertainties of the spectroscopic transitions employed. To show its advantages, the extended MARVEL protocol is applied for the analysis of the complete set of highly accurate H216O transition measurements. As a result, almost 300 highly accurate energy levels of H216O are reported in the energy range of 0-6000 cm-1. Out of the 15 vibrational bands involved in accurately measured rovibrational transitions, the following three have definitely highly accurate empirical rovibrational energies of 8-10 digits of accuracy: (v1v2v3) = (0 0 0), (0 1 0), and (0 2 0), where v1, v2, and v3 stand for the symmetric stretch, bend, and antisymmetric stretch vibrational quantum numbers. The dataset of experimental rovibrational transitions and empirical rovibrational energy levels assembled during this study, both with improved uncertainties, is considerably larger and more accurate than the best previous datasets.

Molecular dimers of methane clathrates: ab initio potential energy surfaces and variational vibrational states.

Motivated by the energetic and environmental relevance of methane clathrates, highly accurate ab initio potential energy surfaces (PESs) have been developed for the three possible dimers of the methane and water molecules: (H2O)2, CH4·H2O, and (CH4)2. While only a single monomer geometry was used for each monomer in the ab initio calculations, the PES parameterization makes it possible to produce distinct surfaces for all isotopologues within the rigid-monomer approximation. The PESs were fitted to computations at the frozen-core coupled-cluster level with single, double, and non-iterative triple excitations, employing basis sets of augmented triple- and quadruple-zeta quality plus bond functions, followed by extrapolations to the complete basis set limit. The long-range parts of the PESs are computed using the asymptotic version of symmetry-adapted perturbation theory based on a density-functional description of the monomers. All PESs are polarizable, i.e., in cluster or condensed-phase applications they approximate many-body effects by the induced dipole polarization model. The PESs were developed in a fully automated procedure applying the autoPES method, which is used for the first time to generate near-spectroscopic quality surfaces. The stationary points (SPs) on the PESs have been determined and compared with literature data. For CH4·H2O, previously unknown SPs have been identified and the first detailed study of the (CH4)2 potential energy landscape has been carried out. The PESs were used in variational quantum nuclear motion computations. For the water dimer, the resulting vibrational transitions are in excellent agreement with available high-resolution spectroscopic data. For (CH4)2, the intermonomer vibrational states are reported for the first time.