Dispersion Control over Molecule Cohesion: Exploiting and Dissecting the Tipping Power of Aromatic Rings.

ConspectusWe have learned over the past years how London dispersion forces can be effectively used to influence or even qualitatively tip the structure of aggregates and the conformation of single molecules. This happens despite the fact that single dispersion contacts are much weaker than competing polar forces. It is a classical case of strength by numbers, with the importance of London dispersion forces scaling with the system size. Knowledge about the tipping points, however difficult to attain, is necessary for a rational design of intermolecular forces. One requires a careful assessment of the competing interactions, either by sensitive spectroscopic techniques for the study of the isolated molecules and aggregates or by theoretical approaches. Of particular interest are the systems close to the tipping point, when dispersion interactions barely outweigh or approach the strength of the other interactions. Such subtle cases are important milestones for a scale-up to realistic multi-interaction situations encountered in the fields of life and materials science. In searching for examples that provide ideal competing interactions in complexes and small clusters, aromatic systems can offer a diverse set of molecules with a variation of dispersion and electrostatic forces that control the dominant and peripheral interactions. Our combined spectroscopic and theoretical investigations provide valuable insights into the balance of intermolecular forces because they typically allow us to switch the aromatic substituent on and off. High-resolution rotational spectroscopy serves as a benchmark for molecular structures, as correct calculations should be based on correct geometries. When discussing the competition with other noncovalent interactions, obvious competitors are directional hydrogen bonds. As a second counterweight to aryl interactions, we will discuss aurophilic/metallophilic interactions, which also have a strong stabilization with a small number of atoms involved. Vibrational spectroscopy is most sensitive to interactions of light atoms, and the competition of OH hydrogen bonds with dispersion forces in a molecular aggregate can be judged well by the OH stretching frequency. Experiments in the gas phase are ideal for gauging the accuracy of quantum chemical predictions free of solvent forces. A tight collaboration utilizing these three methods allows experiment vs experiment vs theory benchmarking of the overall influence of dispersion in molecular structures and energetics.

Weak hydrogen bonding to halogens and chirality communication in propanols: Raman and microwave spectroscopy benchmark theory.

Constitutional and conformational isomers of bromopropanol are vibrationally and rotationally characterised with parallels drawn to the structural chlorine analogues. A previous microwave spectroscopic study of the chloropropanols is re-examined and all systems are explored by Raman jet spectroscopy. For bromine, the entire nuclear quadrupole coupling tensors are accurately determined and compared to their chlorine counterparts. Tensor asymmetry parameters are determined and linked with the hydrogen bond strength as indicated by the downshift of the OH-stretching frequency. The spectroscopic constants derived from the observed transitions are used as benchmarks for a large variety of electronic structure methods followed by harmonic and anharmonic rovibrational treatments. The CCSD(T) electronic structure calculations provide the best performance, in particular once anharmonic and relativistic corrections are applied or implied. Standard DFT approaches vary substantially with respect to their systematic error cancellation across the investigated species, and cost-effective compromises for the different observables are proposed.

Quantum Tunneling Facilitates Water Motion across the Surface of Phenanthrene.

Quantum tunneling is a fundamental phenomenon that plays a pivotal role in the motion and interaction of atoms and molecules. In particular, its influence in the interaction between water molecules and carbon surfaces can have significant implications for a multitude of fields ranging from atmospheric chemistry to separation technologies. Here, we unveil at the molecular level the complex motion dynamics of a single water molecule on the planar surface of the polycyclic aromatic hydrocarbon phenanthrene, which was used as a small-scale carbon surface-like model. In this system, the water molecule interacts with the substrate through weak O-H···π hydrogen bonds, in which phenanthrene acts as the hydrogen-bond acceptor via the high electron density of its aromatic cloud. The rotational spectrum, which was recorded using chirped-pulse Fourier transform microwave spectroscopy, exhibits characteristic line splittings as dynamical features. The nature of the internal dynamics was elucidated in great detail with the investigation of the isotope-substitution effect on the line splittings in the rotational spectra of the H218O, D2O, and HDO isotopologues of the phenanthrene-H2O complex. The spectral analysis revealed a complex internal dynamic showing a concerted tunneling motion of water involving its internal rotation and its translation between the two equivalent peripheral rings of phenanthrene. This high-resolution spectroscopy study presents the observation of a tunneling motion exhibited by the water monomer when interacting with a planar carbon surface with an unprecedented level of detail. This can serve as a small-scale analogue for water motions on large aromatic surfaces, i.e., large polycyclic aromatic hydrocarbons and graphene.

Benchmarking the quadrupolar coupling tensor for chlorine to probe weak-bonding interactions.

Rotational spectroscopy relies on quantum chemical calculations to interpret observed spectra. Among the most challenging molecules to assign are those with additional angular momenta coupling to the rotation, contributing to the complexity of the spectrum. This benchmark study of computational methods commonly used by rotational spectroscopists targets the nuclear quadrupole coupling constants of chlorine containing molecules and the geometry of its complexes and clusters. For each method, the quality of both structural and electronic parameter predictions is compared with the experimental values. Ab initio methods are found to perform best overall in predicting both the geometry of the complexes and the coupling constants of chlorine with moderate computational cost. This cost can be reduced by combining these methods with density functional theory structure optimization, which still yields adequate predictions. This work constitutes a first step in expanding Bailey's quadrupole coupling data set to encompass molecular clusters. [W. C. Bailey, Calculation of Nuclear Quadrupole Coupling Constants in Gaseous State Molecule, 2019, https://nqcc.wcbailey.net/].

Hydrogen-Atom Tunneling in a Homochiral Environment.

The role-exchanging concerted torsional motion of two hydrogen atoms in the homochiral dimer of trans-1,2-cyclohexanediol was characterized through a combination of broadband rotational spectroscopy and theoretical modeling. The results reveal that the concerted tunneling motion of the hydrogen atoms leads to the inversion of the sign of the dipole moment components along the a and b principal axes, due to the interchange motion that cooperatively breaks and reforms one intermolecular hydrogen bond. This motion is also coupled with two acceptor switching motions. The energy difference between the two ground vibrational states arising from this tunneling motion was determined to be 29.003(2) MHz. The corresponding wavefunctions suggest that the two hydrogen atoms are evenly delocalized on two equivalent potential wells, which differs from the heterochiral case where the hydrogen atoms are confined in separate wells, as the permutation-inversion symmetry breaks down. This intriguing contrast in hydrogen-atom behavior between homochiral and heterochiral environments could further illuminate our understanding of the role of chirality in intermolecular interactions and dynamics.

The many forms of alpha-methoxy phenylacetic acid in the gas phase: flexibility, internal dynamics, and their intramolecular interactions.

We present a rotational spectroscopy study of alpha-methoxy phenylacetic acid in the gas phase. This acid is a derivative of mandelic acid and is used in various organic reactions. The conformational landscape of alpha-methoxy phenylacetic acid was explored to gain insight into its intramolecular dynamics. A rich rotational spectrum was obtained using chirped-pulse Fourier transform microwave spectroscopy in the 2-8 GHz range. Five conformers out of six calculated low-energy forms were identified in the spectrum, and the assignment of the 13C singly substituted isotopologues for the lowest-energy conformer led to its accurate structure determination. Splitting patterns were analyzed and attributed to the internal rotation of a methyl top. The analysis of the non-covalent interactions within the molecule highlights the subtle balance in the stabilization of the different conformers. We thus provide high-level structural and intramolecular dynamics information that is also used to benchmark the performance of quantum-chemical calculations.

Rotational analysis of naphthol-aromatic ring complexes stabilized by electrostatic and dispersion interactions.

For complexes involving aromatic species, substitution effects can influence the preferred geometry. Using broadband rotational spectroscopy, we report the structures of three naphthol-aromatic ring complexes with different heteroatoms (furan and thiophene) and alkyl groups (2,5-dimethylfuran). The aim was to analyze the influence of the presence of heteroatoms or alkyl groups on the structure of the complex and the kind of intermolecular forces that control it. Face or edge arrangements can take place in these complexes via π-π or O-H⋯O/O-H⋯π interactions, respectively. All the experimentally observed complexes present O-H⋯O/O-H⋯π interactions with the hydroxyl group, with different structures and intermolecular interactions depending on the heteroatom present in the five-membered aromatic rings, yielding different symmetries in the experimental structure. Structures are experimentally identified through the use of planar moments of inertia. Further results from SAPT calculations show that dispersion and electrostatic interactions contribute similarly to the stabilization of all the studied complexes. These new spectroscopic results shed light on the influence of dispersion and hydrogen bonding in molecular aggregation of systems with substituted aromatic residues.

Competitive tetrel bond and hydrogen bond in benzaldehyde-CO2: characterization via rotational spectroscopy.

The rotational spectrum of the 1 : 1 benzaldehyde-CO2 complex has been investigated using pulsed-jet Fourier transform microwave spectroscopy complemented with quantum chemical calculations. Two isomers, both characterized by one C⋯O tetrel bond (n → π* interaction) and one C-H⋯O hydrogen bond (n → σ* interaction), have been observed in the pulsed jet. Competition between the tetrel bond and the hydrogen bond has been disclosed by natural bond orbital analysis: isomer I is characterized by one dominating OCCO2⋯O tetrel bond (12.6 kJ mol-1) and a secondary (C-H)formyl⋯O hydrogen bond (2.2 kJ mol-1); by contrast, in isomer II the (C-H)phenyl⋯O hydrogen bond (7.6 kJ mol-1) becomes the dominant bond, while the OCCO2⋯O tetrel bond (5.8 kJ mol-1) becomes much weaker with respect to that of isomer I. Using intensity measurements the relative population ratio of the two isomers was estimated to be NI/NII ≈ 2/1.

Van der Waals interactions of the disulfide bond revealed: A microwave spectroscopic study of the diethyl disulfide-argon complex.

The van der Waals complex formed between diethyl disulfide (DEDS) and an argon atom was investigated by pulsed-jet Fourier transform microwave spectroscopy in conjunction with quantum chemical computations. One set of transition lines belonging to the configuration of the global potential energy minimum was measured and assigned. The rotational constants A, B, and C were accurately determined to be 1262.5758(1) MHz, 845.402 12(9) MHz, and 574.006 38(8) MHz, respectively. The distance between the argon atom and the center of mass of the DEDS subunit is 4.075(16) Å. Quantum theory of atoms in molecules and non-covalent interaction analyses reveal that the interactions take place between the argon atom and four sites of the DEDS subunit. Furthermore, the usage of the energy decomposition analysis approach provides further understanding of the characteristics of the van der Waals interactions. Additionally, ab initio calculations and symmetry-adapted perturbation theory analysis of the binary complexes of DEDS with He, Ne, Kr, and Xe atoms were carried out to get further insight into the characteristics of the van der Waal interactions of the disulfide bond.

The Characteristics of Disulfide-Centered Hydrogen Bonds.

The disulfide-centered hydrogen bonds in the three different model systems of diethyl disulfide⋅⋅⋅H2 O/H2 CO/HCONH2 clusters were characterized by high-resolution Fourier transform microwave spectroscopy and quantum chemical computations. The global minimum energy structures for each cluster are experimentally observed and are characterized by one of the three different S-S⋅⋅⋅H-C/N/O disulfide-centered hydrogen bonds and two O⋅⋅⋅H-C hydrogen bonds. Non-covalent interaction and natural bond orbital analyses further confirm the experimental observations. The symmetry-adapted perturbation theory (SAPT) analysis reveals that electrostatic is dominant in diethyl disulfide⋅⋅⋅H2 O/HCONH2 clusters being consistent with normal hydrogen bonds, whilst dispersion takes over in diethyl disulfide⋅⋅⋅H2 CO cluster. Our study gives accurate structural parameters for the disulfide bond involved non-covalent clusters providing important benchmarking data for the theoretical evaluation of more complex systems.

XeOCS: relatively straightforward?

We report a benchmark-quality equilibrium-like structure of the XeOCS complex, obtained from microwave spectroscopy. The experiments are supported by a wide array of highly accurate calculations, expanding the analysis to the complexes of He, Ne, Ar, Kr, Xe, and Hg with OCS. We investigate the trends in the structures and binding energies of the complexes. The assumption that the structure of the monomers does not change significantly upon forming a weakly bound complex is also tested. An attempt at reproducing the r structure of the XeOCS complex with correlated wavefunction theory is made, highlighting the importance of relativistic effects, large basis sets, and inclusion of diffuse functions in extrapolation recipes.

The puzzling hyper-fine structure and an accurate equilibrium geometry of succinic anhydride.

An accurate semiexperimental equilibrium structure of succinic anhydride has been determined from a combination of experiment and theory. The cm-wave and mm-wave rotational spectra of succinic anhydride, 3,4-dihydrofuran-2,5-dione, were recorded in a pulsed supersonic jet using Fourier-transform microwave spectroscopy and in a free-jet using mm-wave absorption spectroscopy. Many lines in the cm-wave spectrum show fine structure and after eliminating all other possibilities the origin of this fine structure is determined to be from spin-spin interaction. Accurate rotational and quartic centrifugal distortion constants are determined. Assignments of 13C and 18O singly substituted isotopologues in natural abundance were used to obtain a substitution geometry for the heavy atoms of succinic anhydride. Theoretical approaches permitted the calculation of a Born-Oppenheimer ab initio structure and the determination of a semiexperimental equilibrium structure in which computed rovibrational corrections were utilized to convert vibrational ground state rotational constants into equilibrium constants. The agreement between the semiexperimental structure and the Born-Oppenheimer ab initio structure is excellent. Succinic anhydride has been shown to have a planar heavy atom equilibrium structure with the effects of a large amplitude vibration apparent in the resultant rotational constants.

New findings from old data: A semi-experimental value for the eQq of the nitrogen atom.

Nuclear quadrupole coupling arises from the interaction of the nuclear quadrupole moment with the electric field gradient. Thus, it is associated with electron occupancy and the electronic structure of molecules. We demonstrate a simple method for planar molecules based on a direct correlation between the out-of-plane quadrupole coupling constant and the electron occupancy in the p orbital perpendicular to the molecular plane. This method is applied to 98 molecular systems containing a 14N quadrupolar nucleus using data from more than 40 years of rotational spectroscopy and comparing the performance of three levels of theory from quantum-chemical computations. From this extensive dataset, we have analyzed chemical properties of molecules, such as the hybridization of the atom, and we could quantify the extent of polarization and resonance processes as well as physical characteristics of the quadrupolar nucleus, such as eQq. This is a constant, which represents the interaction in the hypothetical case of having a single electron in an electronic orbital at the isolated nucleus, and its value has been under debate for a long time. Here, the eQq value has been determined for the 14N nucleus, and the methodology to calculate it for other nuclei is provided.

The first microsolvation step for furans: New experiments and benchmarking strategies.

The site-specific first microsolvation step of furan and some of its derivatives with methanol is explored to benchmark the ability of quantum-chemical methods to describe the structure, energetics, and vibrational spectrum at low temperature. Infrared and microwave spectra in supersonic jet expansions are used to quantify the docking preference and some relevant quantum states of the model complexes. Microwave spectroscopy strictly rules out in-plane docking of methanol as opposed to the top coordination of the aromatic ring. Contrasting comparison strategies, which emphasize either the experimental or the theoretical input, are explored. Within the harmonic approximation, only a few composite computational approaches are able to achieve a satisfactory performance. Deuteration experiments suggest that the harmonic treatment itself is largely justified for the zero-point energy, likely and by design due to the systematic cancellation of important anharmonic contributions between the docking variants. Therefore, discrepancies between experiment and theory for the isomer abundance are tentatively assigned to electronic structure deficiencies, but uncertainties remain on the nuclear dynamics side. Attempts to include anharmonic contributions indicate that for systems of this size, a uniform treatment of anharmonicity with systematically improved performance is not yet in sight.

Tetrel bonds and conformational equilibria in the formamide-CO2 complex: a rotational study.

The rotational spectra of the complex formamide-CO2 have been measured and assigned by pulsed jet Fourier transform microwave spectroscopy. Two isomers of the complex have been detected where a CO tetrel bond dominates the interactions, and either N-HO or C-HO forms a secondary linkage. Bader's quantum theory of atoms in molecules and Johnson's non-covalent interaction analyses were applied to unveil the intermolecular binding sites and energetic properties in the complex. Relative intensity measurements on a set of µa-type transitions allowed estimating the relative population of the observed two isomers as NI/NII ≈ 18/1.

Ab initio conformational analysis of 1,2,3,4-tetrahydroquinoline and the high-resolution rotational spectrum of its lowest energy conformer.

The saturated part of the 1,2,3,4-tetrahydroquinoline (THQ) molecule allows for the possibility of multiple conformers' existence. High-resolution microwave spectroscopy, supported by high-level quantum chemistry calculations, was used to determine the precise molecular structures of the conformers of THQ. Via the MP2 calculations, we were able to discriminate four stable conformations, i.e. two pairs of energetically equivalent enantiomorphic conformers. The results of the calculations also indicate that energetically non-equivalent conformers are separated by a low energy barrier (104 cm-1) that allows for conformational cooling to occur. The high resolution rotational spectrum with resolved hyperfine structure in the frequency range of 7-20 GHz was obtained using both the In-phase/quadrature-phase-Modulation Passage-Acquired-Coherence Technique (IMPACT) and the coaxially oriented beam resonator arrangement (COBRA) to perform Fourier transform microwave (FTMW) spectroscopy. The precise values of the rotational constants, 14N nuclear hyperfine coupling parameters and centrifugal distortion parameters were determined from the measured transition frequencies. Based on our experimental results, only the most stable enantiomeric pair of THQ contributes to the rotational spectrum under the conditions of our experiment as the less stable conformers seem to efficiently relax to the lower energy conformers. Thus the experimentally evaluated molecular constants unambiguously define the lowest energy conformer of 1,2,3,4-tetrahydroquinoline.

Unveiling the Sulfur-Sulfur Bridge: Accurate Structural and Energetic Characterization of a Homochalcogen Intermolecular Bond.

By combining rotational spectroscopy in supersonic expansion with the capability of state-of-the-art quantum-chemical computations in accurately determining structural and energetic properties, the genuine nature of a sulfur-sulfur chalcogen bond between dimethyl sulfide and sulfur dioxide has been unveiled in a gas-jet environment free from collision, solvent and matrix perturbations. A SAPT analysis pointed out that electrostatic S⋅⋅⋅S interactions play the dominant role in determining the stability of the complex, largely overcoming dispersion and C-H⋅⋅⋅O hydrogen-bond contributions. Indeed, in agreement with the analysis of the quadrupole-coupling constants and of the methyl internal rotation barrier, the NBO and NOCV/CD approaches show a marked charge transfer between the sulfur atoms. Based on the assignment of the rotational spectra for 7 isotopologues, an accurate semi-experimental equilibrium structure for the heavy-atom backbone of the molecular complex has been determined, which is characterized by a S⋅⋅⋅S distance (2.947(3) Å) well below the sum of van der Waals radii.

The covalent interaction between dihydrogen and gold: A rotational spectroscopic study of H2-AuCl.

The pure rotational transitions of H2-AuCl have been measured using a pulsed-jet cavity Fourier transform microwave spectrometer equipped with a laser ablation source. The structure was found to be T-shaped, with the H-H bond interacting with the gold atom. Both 35Cl and 37Cl isotopologues have been measured for both ortho and para states of H2. Rotational constants, quartic centrifugal distortion constants, and nuclear quadrupole coupling constants for gold and chlorine have been determined. The use of the nuclear spin-nuclear spin interaction terms Daa, Dbb, and Dcc for H2 were required to fit the ortho state of hydrogen, as well as a nuclear-spin rotation constant Caa. The values of the nuclear quadrupole coupling constant of gold are χaa=-817.9929(35) MHz, χbb=504.0(27) MHz, and χcc=314.0(27). This is large compared to the eQq of AuCl, 9.63 312(13) MHz, which indicates a strong, covalent interaction between gold and dihydrogen.

A Study of 2-Iodobutane by Rotational Spectroscopy.

Rotational transitions belonging to 2-iodobutane (sec-butyl-iodide, CH3CHICH2CH3) were measured over the frequency range 5.5-16.5 GHz via jet-pulsed Fourier transform microwave spectroscopy. The complete nuclear quadrupole coupling tensor of iodine, χ, was obtained for the gauche (g)-, anti (a)-, and gauche' (g')-conformers as well as the four (13)C isotopologues of the gauche species. Rotational constants, centrifugal distortion constants, quadrupole coupling constants, and nuclear spin-rotation constants were determined for each species. Changes in χ of the iodine nucleus, resulting from conformational and isotopic differences, are discussed. Isotopic substitution of g-2-iodobutane allowed for an rs structure to be determined for the carbon backbone. Additionally, isotopic substitution in conjunction with an ab initio structure allowed for a fit of various r0 structural parameters belonging to g-2-iodobutane.

Microwave Spectra, Structure, and Ring-Puckering Vibration of Octafluorocyclopentene.

The rotational spectra of octafluorocyclopentene (C5F8) has been measured for the first time using pulsed jet Fourier transform microwave spectroscopy in a frequency range of 6 to 16 GHz. As in the molecule cyclopentene, the carbon ring is nonplanar, and inversion through the plane results in an inversion pair of ground state vibrational energy levels with an inversion splitting of 18.4 MHz. This large amplitude motion leads to the vibration-rotation coupling of energy levels. The symmetric double minimum ring-puckering potential function was calculated, resulting in a barrier of 222 cm-1. The rotational constants A0 = 962.9590(1) MHz, B0 = 885.1643(4) MHz, C0 = 616.9523(4) MHz, A1 = 962.9590(1) MHz, B1 = 885.1643(4) MHz, C1 = 616.9528(4) MHz, and two centrifugal distortion constants for each state were determined for the parent species and all 13C isotopologues. A mixed coordinate molecular structure was determined from a least-squares fit of the ground state rotational constants of the parent and each 13C isotopologue combined with the equilibrium bond lengths and angles from quantum chemical calculations.