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.

Gas-Phase Characterization of Adipic Acid, 6-Hydroxycaproic Acid, and Their Thermal Decomposition Products by Rotational Spectroscopy.

We report the spectroscopic investigation of two bifunctional aliphatic carboxylic acids, namely, adipic acid and 6-hydroxycaproic acid, in the gas phase by combining high-resolution rotational spectroscopy and supersonic expansions. Their pure rotational spectra were successfully identified and characterized. However, due to the low thermal stability of these two chemicals, the measured rotational spectra were significantly congested with transitions corresponding to their decomposition products upon heating. We observed cyclopentanone and adipic anhydride in the spectrum of adipic acid and ε-caprolactone and its monohydrate in the spectrum of 6-hydroxycaproic acid. On the basis of the distinct fingerprints of both carboxylic acids and a series of their decomposition products, the spectra were analyzed in a time-segmented manner. This provides valuable insights into the thermal decomposition mechanisms of these two samples over time, which highlights the robustness of microwave spectroscopy as a potent tool for analyzing complex chemical mixtures in a species-, isomer-, and conformer-selective way.

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.

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.

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

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.

Interaction Types in C6H5(CH2)nOH-CO2 (n = 0-4) Determined by the Length of the Side Alkyl Chain.

C6H5(CH2)nOH-CO2 complexes have been investigated using rotational spectroscopy (n = 0-2) complemented by quantum chemical calculations (n = 0-4), which implies that the side alkyl chain length can determine the types of intermolecular interactions. Unlike the in-plane C···O tetrel bond in phenol-CO2, the π*CO2···πaromatic interaction has been shown to link CO2 to phenylmethanol and 2-phenylethanol, which is, to the best of our knowledge, the first time it has been demonstrated by rotational spectroscopy. Further elongations of the side alkyl chain gradually increase the energies of intramolecular hydrogen bonds in 3-phenylpropanol and 4-phenylbutanol so that CO2 cannot break it. CO2 will be pushed farther from the monomers and link with the -OH group through a dominating C···O tetrel bond. Our observations would allow, with the choice of the proper length of the side alkyl chain, new strategies for engineering C···πaromatic-centered noncovalent bonding schemes for the capture, utilization, and storage of CO2.

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.

Internal rotation and chlorine nuclear quadrupole coupling in 2-chloro-4-fluorotoluene explored by microwave spectroscopy and quantum chemistry.

2-Chloro-4-fluorotoluene was investigated using a combination of molecular jet Fourier transform microwave spectroscopy in the frequency range from 5 to 21 GHz and quantum chemistry. The molecule experiences an internal rotation of the methyl group, which causes fine splittings of all rotational transitions into doublets with separation on the order of a few tens of kHz. In addition, hyperfine effects originating from the chlorine nuclear quadrupole moment coupling its nuclear spin to the end-over-end rotation of the molecule are observed. The torsional barrier was derived using both the rho and the combined-axis-method, giving a value of 462.5(41) cm-1. Accurate rotational constants and quadrupole coupling constants were determined for the 35Cl and 37Cl isotopologues and compared with Bailey's semi-experimental quantum chemical predictions. The gas phase molecular structure was deduced from the experimental rotational constants supplemented with those calculated by quantum chemistry at various levels of theory. The values of the methyl torsional barrier and chlorine nuclear quadrupole coupling constants were compared with the theoretical predictions and with those of other chlorotoluene derivatives.

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.

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.

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.

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.

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.