New Ab Initio Potential Energy Surfaces for NH3 Constructed from Explicitly Correlated Coupled-Cluster Methods

Rovibrational Line Lists of Triplet and Singlet Methylene

The methylene molecule (CH2) is a short-lived radical with lacking data on its spectral line intensities. Although the lifetime of CH2 is extremely short under Earth's conditions, it exists in a free form in interstellar media. CH2 is an important intermediate species in chemical reactions associated with the formation and destruction of complex hydrocarbons. We present the first rovibrational line lists of CH2 in its ground triplet and first excited singlet electronic state. To this end, our previously developed accurate ab initio potential energy surface (PES) was used for the ground electronic triplet state [Egorov et al. J. Comp. Chem. 2024. V. 45. (2). P. 83] while a new PES for the singlet state was constructed in this work using the single-reference coupled cluster approach [CCSD(T)] combined with the extrapolation to the complete basis set (CBS) limit based on the correlation-consistent orbital basis sets with the core-valence electron correlation effects [aug-cc-pCVXZ, X = T, Q, 5, and 6]. In addition, the contributions to the correlation energy from highly excited Slater determinants [CC(n), n = 3-5] were included as well as the scalar relativistic effects and DBOC. The most accurate description of the infrared band origins of singlet CH2 was thus achieved for the energy range where the impact of the nonadiabatic coupling due to the Renner-Teller effect can be neglected. To obtain the probabilities of the rovibrational transitions, new ab initio DMSs were constructed both for the triplet and singlet CH2 using the CCSD(T)/aug-cc-pCVQZ approach. Finally, the absorption spectra of triplet and singlet methylene were predicted from the variationally computed line lists.

Gravitationally-induced wave function collapse time for molecules

The Diósi-Penrose model states that the wave function collapse ending a quantum superposition occurs due to the instability of coexisting gravitational potentials created by distinct geometric conformations of the system in different states. The Heisenberg time-energy principle can be invoked to estimate the collapse time for the energy associated with this instability, the gravitational self-energy. This paper develops atomistic models to calculate the Diósi-Penrose collapse time. It applies them to a range of systems, from small molecules to large biological structures and macroscopic systems. An experiment is suggested to test the Diósi-Penrose hypothesis, and we critically examine the model, highlighting challenges from an atomistic perspective, such as gravitational self-energy saturation and limited extensivity.

Direct flipping dynamics and quantized enrichment of chirality at single-molecule resolution

Chirality is an important aspect of nature, and numerous macroscopic methods have been developed to understand and control chirality. For the chiral tertiary amines, their flexible flipping process makes it possible to achieve high chiral controllability without bond formation and breaking. Here, we present a type of stable chiral single-molecule devices formed by tertiary amines, using graphene-molecule-graphene single-molecule junctions. These single-molecule devices allow real-time, in situ, and long-time measurements of the flipping process of an individual chiral nitrogen center with high temporal resolution. Temperature- and bias voltage-dependent experiments, along with theoretical investigations, revealed diverse chiral intermediates, indicating the regulation of the flipping dynamics by energy-related factors. Angle-dependent measurements further demonstrated efficient enrichment of chiral states using linearly polarized light by a symmetry-related factor. This approach offers a reliable means for understanding the chirality's origin, elucidating microscopic chirality regulation mechanisms, and aiding in the design of effective drugs.

Accurate ab initio potential energy surface, rovibrational energy levels and resonance interactions of triplet ( X ~ 3 B1 ) methylene

In this work, we report rovibrational energy levels for four isotopologues of methylene (CH2 , CHD, CD2 , and 13 CH2 ) in their ground triplet electronic state (X ~ 3 B1 ) from variational calculation up to ~10,000 cm-1 and using a new accurate ab initio potential energy surface (PES). Triplet methylene exhibits a large-amplitude bending vibration and can reach a quasilinear configuration due to its low barrier (~2000 cm-1 ). To construct the ab initio PES, the Dunning's augmented correlation-consistent core-valence orbital basis sets were employed up to the sextuple-ζ quality [aug-cc-pCVXZ, X = T, Q, 5, and 6] combined with the single- and double-excitation unrestricted coupled cluster approach with a perturbative treatment of triple excitations [RHF-UCCSD(T)]. We have shown that the accuracy of the ab initio energies is further improved by including the corrections due to the scalar relativistic effects, DBOC and high-order electronic correlations. For the first time, all the available experimental rovibrational transitions were reproduced with errors less than 0.12 cm-1 , without any empirical corrections. Unlike more "traditional" nonlinear triatomic molecules, we have shown that even the energies of the ground vibrational state (000) with rather small rotational quantum numbers are strongly affected by the very pronounced rovibrational resonance interactions. Accordingly, the polyad structure of the vibrational levels of CH2 and CD2 was analyzed and discussed. The comparison between the energy levels obtained from the effective Watson A-reduced Hamiltonian, from the generating-function approach and from a variational calculation was given.

A numerical-tensorial "hybrid" nuclear motion Hamiltonian and dipole moment operator for spectra calculation of polyatomic nonrigid molecules

The analysis and modeling of high-resolution spectra of nonrigid molecules require a specific Hamiltonian and group-theoretical formulation that differs significantly from that of more familiar rigid systems. Within the framework of Hougen-Bunker-Johns (HBJ) theory, this paper is devoted to the construction of a nonrigid Hamiltonian based on a suitable combination of numerical calculations for the nonrigid part in conjunction with the irreducible tensor operator method for the rigid part. For the first time, a variational calculation from ab initio potential energy surfaces is performed using the HBJ kinetic energy operator built from vibrational, large-amplitude motion, and rotational tensor operators expressed in terms of curvilinear and normal coordinates. Group theory for nonrigid molecules plays a central role in the characterization of the overall tunneling splittings and is discussed in the present approach. The construction of the dipole moment operator is also examined. Validation tests consisting of a careful convergence study of the energy levels as well as a comparison of results obtained from independent computer codes are given for the nonrigid molecules CH2, CH3, NH3, and H2O2. This work paves the way for the modeling of high-resolution spectra of larger nonrigid systems.

New Theoretical Infrared Line List for the Methyl Radical with Accurate Vibrational Band Origins from High-Level Ab Initio Calculations

In the present work, high-level ab initio calculations were carried out for the ground electronic state of the methyl radical (CH₃). Dunning's augmented correlation-consistent orbital basis sets were employed up to the quintuple-ζ valence quality with the core-valence electron correlation [aug-cc-pCV5Z] combined with the single- and double-excitation unrestricted coupled-cluster approach with a perturbative treatment of triple excitations [RHF-UCCSD(T)]. The explicitly correlated version of the coupled-cluster approach [RHF-UCCSD(T)-F12x{x = a, b}] was additionally applied with the core-valence cc-pCVQZ-F12 basis set in order to study convergence with respect to the basis set size. The contributions beyond the coupled-cluster level of the theory like Douglas-Kroll-Hess scalar relativistic Hamiltonian, diabatic Born-Oppenheimer corrections, and high-order electronic correlations have been included into the ab initio potential energy surfaces (PESs). It is shown that the theoretical band origins of CH₃ converge gradually to the experimental values when applying the ab initio PESs using the aug-cc-pCVXZ [X = T, Q, and 5] basis sets. For the first time, all available experimental band origins of the gaseous CH₃ are reproduced within an accuracy of 0.2 cm^-1 using a newly developed PES extrapolated to the complete basis set limit [CBS(TQ5Z)]. The reached accuracy is one order of magnitude better than that of the best available calculations. A new theoretical infrared line list was generated for astrophysical applications using an ab initio dipole moment surface computed at the RHF-UCCSD(T)/aug-cc-pCVQZ level of the theory. The manifestation of a large-amplitude motion in CH₃ is also discussed.

Novel methodology for systematically constructing global effective models from ab initio-based surfaces: A new insight into high-resolution molecular spectra analysis

In this paper, a novel methodology is presented for the construction of ab initio effective rotation-vibration spectroscopic models from potential energy and dipole moment surfaces. Non-empirical effective Hamiltonians are obtained via the block-diagonalization of selected variationally computed eigenvector matrices. For the first time, the derivation of an effective dipole moment is carried out in a systematic way. This general approach can be implemented quite easily in most of the variational computer codes and turns out to be a clear alternative to the rather involved Van Vleck perturbation method. Symmetry is exploited at all stages to translate first-principles calculations into a set of spectroscopic parameters to be further refined on experiment. We demonstrate on H₂CO, PH₃, CH₄, C₂H₄, and SF₆ that the proposed effective model can provide crucial information to spectroscopists within a very short time compared to empirical spectroscopic models. This approach brings a new insight into high-resolution spectrum analysis of polyatomic molecules and will be also of great help in the modeling of hot atmospheres where completeness is important.