Potential-Energy Surfaces for Ring-Puckering Motions of Flexible Cyclic Molecules through Cremer-Pople Coordinates: Computation, Analysis, and Fitting

Deciphering the conformational preference and ionization dynamics of tetrahydrofuran: Insights from advanced spectroscopic techniques

Tetrahydrofuran (THF) has garnered significant attention due to its pivotal role in biological and chemical processes. The diverse array of conformations exhibited by THF profoundly impacts its reactivity and interactions with other molecules. Understanding these conformational preferences is crucial for comprehending its molecular behavior. In this study, we utilize infrared (IR) resonant vacuum ultraviolet photoionization/mass-analyzed threshold ionization (VUV-PI/MATI) mass spectroscopies to capture distinctive vibrational spectra of individual conformers, namely, "twisted" and "bent," within THF. Our conformer-specific vibrational spectra provide valuable insights into the relative populations of these two conformers. The analysis reveals that the twisted (C2) conformer is more stable than the bent (CS) conformer by 17 ± 15 cm-1. By precisely tuning the VUV photon energy to coincide with vibrational excitation via IR absorption, we selectively ionize specific conformers, yielding two-photon IR + VUV-PI/MATI spectra corresponding to the twisted and bent conformers. This investigation conclusively affirms that both the twisted and bent conformers coexist in the neutral state, while only the twisted conformer exists in the cationic state. These findings not only bridge gaps in existing knowledge but also provide profound insights into the behavior of this pivotal molecule in the realms of biology and medicine.

Quantum chemistry meets high-resolution spectroscopy for characterizing the molecular bricks of life in the gas-phase

Computation of accurate geometrical structures and spectroscopic properties of large flexible molecules in the gas-phase is tackled at an affordable cost using a general exploration/exploitation strategy. The most distinctive feature of the approach is the careful selection of different quantum chemical models for energies, geometries and vibrational frequencies with the aim of maximizing the accuracy of the overall description while retaining a reasonable cost for all the steps. In particular, a composite wave-function method is used for energies, whereas a double-hybrid functional (with the addition of core-valence correlation) is employed for geometries and harmonic frequencies and a cheaper hybrid functional for anharmonic contributions. A thorough benchmark based on a wide range of prototypical molecular bricks of life shows that the proposed strategy is close to the accuracy of state-of-the-art composite wave-function methods, and is applicable to much larger systems. A freely available web-utility post-processes the geometries optimized by standard electronic structure codes paving the way toward the accurate yet not prohibitively expensive study of medium- to large-sized molecules by experimentally-oriented researchers.

Photodissociation dynamics of tetrahydrofuran at 193 nm

Tetrahydrofuran (THF), a cyclic ether with the chemical formula C4H8O, can be considered the simplest analog of the deoxyribose backbone component of deoxyribonucleic acid (DNA). As such, it provides a useful model for probing the photochemistry of such biomolecular motifs. We present a velocity-map imaging study into the ultraviolet dissociation of THF at a wavelength of 193 nm. Excitation to the S1 state occurs via a 3s ← n transition involving a lone-pair electron on the oxygen atom, and has been shown by other authors to result in rapid ring opening via cleavage of one of the C-O bonds to form a ring-opened C4H8O diradical, followed by C-C bond cleavage over a longer timescale to form either OCH2 + C3H6 products (Channel 1a), HOCH2 + C2H5 products (Channel 1b), or OCH2CH2 + C2H4 products (Channel 2). The C2H4O products formed via Channel 2 are unstable on the timescale of our experiment and dissociate further to form CH3 and CHO. We also observe a number of minor products resulting from H or H2 loss from the primary photofragments. The speed distributions observed for all photofragments are broad, indicating excitation of a range of rotational and vibrational states of the products. The angular distributions of the photofragments show an interesting speed dependence: the slowest products have almost isotropic angular distributions, but the magnitude of the recoil anisotropy increases monotonically with photofragment speed. The fastest products exhibit highly anisotropic angular distributions, with the recoil anisotropy parameter β approaching its limiting value of -1 (-0.75 for Channel 1 and -0.5 for Channel 2). This behaviour is attributed to the range of timescales over which the diradical intermediate dissociates into the observed photofragments. Rapid dissociation leads to fast photofragments which retain the correlation between the transition dipole moment for the S1 ← S0 excitation (which lies perpendicular to the ring) and the photofragment velocities (which lie predominantly in the plane of the ring). Slow dissociation results in a high degree of energy redistribution into internal modes, slower photofragments, and loss of correlation between the photofragment velocities and the transition dipole. The higher barrier associated with dissociation via Channel 2 suggests somewhat longer lifetimes for the diradical intermediate and is consistent with a corresponding reduction in the maximum observed value for β.

DFT Meets Wave-Function Methods for Accurate Structures and Rotational Constants of Histidine, Tryptophan, and Proline

A new computational strategy has been applied to the conformational and spectroscopic properties in the gas phase of amino acids with very distinctive features, ranging from different tautomeric forms (histidine) to ring puckering (proline), and heteroaromatic structures with non-equivalent rings (tryptophan). The integration of modern double-hybrid functionals and wave-function composite methods has allowed us to obtain accurate results for a large panel of conformers with reasonable computer times. The remarkable agreement between computations and microwave experiments allows an unbiased interpretation of the latter in terms of stereoelectronic effects.

Potential Energy Surfaces Sampled in Cremer-Pople Coordinates and Represented by Common Force Field Functionals for Small Cyclic Molecules

The complex conformations of the cyclic moieties impact the physical and chemical properties of molecules. In this work, we chose 22 molecules of four-, five-, and six-membered rings and performed a thorough conformational sampling using Cremer-Pople coordinates. With consideration of symmetries, we obtained a total of 1504 conformational structures for four-membered, 5576 for five-membered, and 13509 for six-membered rings. All well-known and many less well-known conformers for each molecule were identified. We represented the potential energy surfaces (PESs) by fitting the data to common analytical force field (FF) functional forms. We found that the general features of PESs can be described by the essential FF functional forms; however, the accuracy of representation can be improved remarkably by including the torsion-bond and torsion-angle coupling terms. The best fit yields R-squared (R²) values close to 1.0 and mean absolute errors in energy less than 0.3 kcal/mol.

Low-temperature heat capacity and pseudorotation in 2 methyltetrahydrofuran

The heat capacity and phase transitions of 2-methyltetrahydrofuran in the temperature range from 7 to 350 K were measured by adiabatic calorimetry. The smoothed molar thermodynamic functions in the condensed...

Optical Activity in Saturated Cyclic Amines: Untangling the Roles of Nitrogen-Inversion and Ring-Puckering Dynamics

The dispersive optical activity of two saturated cyclic amines, (R)-2-methylpyrrolidine (R-2MPY) and (S)-2-methylpiperidine (S-2MPI), has been interrogated under isolated and solvated conditions to elucidate the roles of large-amplitude motion associated with nitrogen-center inversion and ring-puckering dynamics. Experimental optical rotatory dispersion profiles were almost mirror images of one another and displayed parallel solvent dependencies. Quantum-chemical analyses built on density-functional and coupled-cluster methods revealed four low-lying conformers for each molecule, which are distinguished by axial/equatorial orientations of their amino hydrogens and methyl substituents. Chiroptical signatures predicted for these species were combined through an independent-conformer ansatz to simulate the ensemble-averaged response, with a polarizable continuum model (PCM) being used to treat implicit solute-solvent interactions. The intrinsic behavior observed for isolated (gaseous) R-2MPY and S-2MPI was reproduced best by merging coupled-cluster (CCSD) estimates of rotatory powers with thermal population fractions deduced from complete basis set (CBS-APNO) free-energy calculations. Although prior claims of sizable chiroptical contributions arising from helically twisted (chiral) heterocyclic frameworks could be discounted, less satisfactory agreement between experiment and theory was realized for solution phases. Response properties sustained modest isomer-dependent changes in the presence of PCM solvation, but the corresponding energy metrics showed systematic trends, whereby structures having larger electric-dipole moments were stabilized preferentially in media of high polarity. Despite the fact that R-2MPY conformations were predicted to undergo a progressive reordering of their relative energies across the six solvents of interest, S-2MPI was found to exhibit more pronounced solvent-induced perturbations at long wavelengths (viz., in regions far removed from electronic resonances). Experimental results are discussed in terms of the distinct ring-puckering mechanisms for R-2MPY and S-2MPI, which are expected to be dominated by hindered pseudorotation among envelope/twist motifs and semi-inversion between chairlike antipodes, respectively.

Interplay of Rotational and Pseudorotational Motions in Flexible Cyclic Molecules

Solutions to the time-independent nuclear Schrödinger equation associated with the pseudorotational motion of three flexible cyclic molecules are presented and discussed. Structural relaxations related to the pseudorotational motion are described as functions of a pseudorotation angle ϕ which is formulated according to the definition of ring-puckering coordinates originally proposed by Cremer and Pople ( J. Am. Chem. Soc. 1975, 97 (6), 1354-1358). In order to take into account the interplay between pseudorotational and rotational motions, the rovibrational Hamiltonian matrices are formulated for the rotational quantum numbers J = 0 and J = 1. The rovibrational Hamiltonian matrices are constructed and diagonalized using a Python program developed by the authors. Suitable algorithms for (i) the construction of one-dimensional cuts of potential energy surfaces along the pseudorotation angle ϕ and (ii) the assignment of the vibrorotational wave functions (which are needed for the automatic calculation of rotational transition energies J = 0 → J = 1) are described and discussed.

Impact of Deuteration and Temperature on Furan Ring Dynamics

Despite significant progress in conformational analysis of cyclic molecules, the number of computational studies is still limited while most of that available in the literature data have been obtained long time ago with outdated methods. In present research, we have studied temperature driven conformational changes of the furan ring at three different temperatures. Additionally, the effect of deuteration on the ring dynamics is discussed; in addition, the aromaticity indices following the Bird and HOMA schemes are computed along all trajectories. Our ab initio molecular dynamic simulations revealed that deuteration has changed the furan ring dynamics and the obvious consequences; in addition, the shape and size of molecule are expected to be different.

Understanding Ring Puckering in Small Molecules and Cyclic Peptides

The geometry of a molecule plays a significant role in determining its physical and chemical properties. Despite its importance, there are relatively few studies on ring puckering and conformations, often focused on small cycloalkanes, 5- and 6-membered carbohydrate rings, and specific macrocycle families. We lack a general understanding of the puckering preferences of medium-sized rings and macrocycles. To address this, we provide an extensive conformational analysis of a diverse set of rings. We used Cremer-Pople puckering coordinates to study the trends of the ring conformation across a set of 140 000 diverse small molecules, including small rings, macrocycles, and cyclic peptides. By standardizing using key atoms, we show that the ring conformations can be classified into relatively few conformational clusters, based on their canonical forms. The number of such canonical clusters increases slowly with ring size. Ring puckering motions, especially pseudo-rotations, are generally restricted and differ between clusters. More importantly, we propose models to map puckering preferences to torsion space, which allows us to understand the inter-related changes in torsion angles during pseudo-rotation and other puckering motions. Beyond ring puckers, our models also explain the change in substituent orientation upon puckering. We also present a novel knowledge-based sampling method using the puckering preferences and coupled substituent motion to generate ring conformations efficiently. In summary, this work provides an improved understanding of general ring puckering preferences, which will in turn accelerate the identification of low-energy ring conformations for applications from polymeric materials to drug binding.

Systematic description of molecular deformations with Cremer-Pople puckering and deformation coordinates utilizing analytic derivatives: Applied to cycloheptane, cyclooctane, and cyclo[18]carbon

The conformational properties of ring compounds such as cycloalkanes determine to a large extent their stability and reactivity. Therefore, the investigation of conformational processes such as ring inversion and/or ring pseudorotation has attracted a lot of attention over the past decades. An in-depth conformational analysis of ring compounds requires mapping the relevant parts of the conformational energy surface at stationary and also at non-stationary points. However, the latter is not feasible by a description of the ring with Cartesian or internal coordinates. We provide in this work, a solution to this problem by introducing a new coordinate system based on the Cremer-Pople puckering and deformation coordinates. Furthermore, analytic first- and second-order derivatives of puckering and deformation coordinates, i.e., B-matrices and D-tensors, were developed simplifying geometry optimization and frequency calculations. The new coordinate system is applied to map the potential energy surfaces and reaction paths of cycloheptane (C₇H₁₄), cyclooctane (C₈H₁₆), and cyclo[18]carbon (C₁₈) at the quantum chemical level and to determine for the first time all stationary points of these ring compounds in a systematic way.

Describing Polytopal Rearrangements of Fluxional Molecules with Curvilinear Coordinates Derived from Normal Vibrational Modes: A Conceptual Extension of Cremer-Pople Puckering Coordinates

In this work a new curvilinear coordinate system is presented for the comprehensive description of polytopal rearrangements of N-coordinate compounds (N = 4-7) and systems containing an N-coordinate subunit. It is based on normal vibrational modes and a natural extension of the Cremer-Pople puckering coordinates ( J. Am. Chem. Soc. 1975, 97, 1354) together with the Zou-Izotov-Cremer deformation coordinates ( J. Phys. Chem. A 2011, 115, 8731) for ring structures to N-coordinate systems. We demonstrate that the new curvilinear coordinates are ideal reaction coordinates describing fluxional rearrangement pathways by revisiting the Berry pseudorotation and the lever mechanism in sulfur tetrafluoride, the Berry pseudorotation and two Muetterties' mechanisms in pentavalent compounds, the chimeric pseudorotation in iodine pentafluoride, Bailar and Ray-Dutt twists in hexacoordinate tris-chelates as well as the Bartell mechanism in iodine heptafluoride. The results of our study reveal that this dedicated curvilinear coordinate system can be applied to most coordination compounds opening new ways for the systematic modeling of fluxional processes.

Chemical promenades: Exploring potential-energy surfaces with immersive virtual reality

The virtual-reality framework AVATAR (Advanced Virtual Approach to Topological Analysis of Reactivity) for the immersive exploration of potential-energy landscapes is presented. AVATAR is based on modern consumer-grade virtual-reality technology and builds on two key concepts: (a) the reduction of the dimensionality of the potential-energy surface to two process-tailored, physically meaningful generalized coordinates, and (b) the analogy between the evolution of a chemical process and a pathway through valleys (potential wells) and mountain passes (saddle points) of the associated potential energy landscape. Examples including the discovery of competitive reaction paths in simple A + BC collisional systems and the interconversion between conformers in ring-puckering motions of flexible rings highlight the innovation potential that augmented and virtual reality convey for teaching, training, and supporting research in chemistry.

Toward Fully Unsupervised Anharmonic Computations Complementing Experiment for Robust and Reliable Assignment and Interpretation of IR and VCD Spectra from Mid-IR to NIR: The Case of 2,3-Butanediol and trans-1,2-Cyclohexanediol

The infrared (IR) and vibrational circular dichroism (VCD) spectra of 2,3-butanediol and trans-1,2-cyclohexanediol from 900 to 7500 cm^–1 (including mid-IR, fundamental CH and OH stretchings, and near-infrared regions) have been investigated by a combined experimental and computational strategy. The computational approach is rooted in density functional theory (DFT) computations of harmonic and leading anharmonic mechanical, electrical, and magnetic contributions, followed by a generalized second-order perturbative (GVPT2) evaluation of frequencies and intensities for all the above regions without introducing any ad hoc scaling factor. After proper characterization of large-amplitude motions, all resonances plaguing frequencies and intensities are taken into proper account. Comparison of experimental and simulated spectra allows unbiased assignment and interpretation of the most interesting features. The reliability of the GVPT2 approach for OH stretching fundamentals and overtones is confirmed by the remarkable agreement with a local mode model purposely tailored for the latter two regions. Together with the specific interest of the studied molecules, our results confirm that an unbiased assignment and interpretation of vibrational spectra for flexible medium-size molecules can be achieved by means of a nearly unsupervised reliable, robust, and user-friendly DFT/GVPT2 model.

Unbiased Determination of Absolute Configurations by vis-à-vis Comparison of Experimental and Simulated Spectra: The Challenging Case of Diplopyrone

A new experimental-computational strategy for the determination of the absolute configuration (AC) of complex chiral molecules is proposed by combining diverse experimental spectroscopies with quantum-mechanical simulations well beyond the current computational practice. Key features are the conformer search and relative stability evaluation performed by a new stochastic two-level tool followed by a vis-à-vis comparison of experimental and computed spectra without any ad hoc adjustment. The entire computational procedure is embedded in the user-friendly VMS software, and its reliability is granted by the inclusion of mechanic/electric/magnetic anharmonicity as well as ro-vibrational and vibronic couplings by means of generalized perturbation theory in conjunction with double-hybrid functionals combined with empirical dispersion contributions and suitable basis sets. To test and validate the new approach, the puzzling case of diplopyrone, a fungal phytotoxic metabolite, has been chosen: the close match between new experimental and simulated infrared absorption and vibrational circular dichroism spectra has led to the unbiased evaluation of its AC.