π-Stacking Isomerism in Polycyclic Aromatic Hydrocarbons: The 2-Naphthalenethiol Dimer

Competition between O-H and S-H Intermolecular Interactions in Conformationally Complex Systems: The 2-Phenylethanethiol and 2-Phenylethanol Dimers

Noncovalent interactions involving sulfur centers play a relevant role in biological and chemical environments. Yet, detailed molecular descriptions are scarce and limited to very simple model systems. Here we explore the formation of the elusive S-H···S hydrogen bond and the competition between S-H···O and O-H···S interactions in pure and mixed dimers of the conformationally flexible molecules 2-phenylethanethiol (PET) and 2-phenylethanol (PEAL), using the isolated and size-controlled environment of a jet expansion. The structure of both PET-PET and PET-PEAL dimers was unraveled through a comprehensive methodology that combined rotationally resolved microwave spectroscopy, mass-resolved isomer-specific infrared laser spectroscopy, and quantum chemical calculations. This synergic experimental-computational approach offered unique insights into the potential energy surface, conformational equilibria, molecular structure, and intermolecular interactions of the dimers. The results show a preferential order for establishing hydrogen bonds following the sequence S-H···S < S-H···O ≲ O-H···S < O-H···O, despite the hydrogen bond only accounting for a fraction of the total interaction energy.

Rotational Spectroscopy and Conformational Flexibility of 2-Phenylethanethiol: The Dominant S-H⋅⋅⋅π Intramolecular Hydrogen Bond

We present a rotational-computational investigation of the aromatic mercaptan 2-phenylethanethiol, addressing its potential energy surface, conformational equilibrium, internal dynamics and intramolecular interactions. The experiment used broadband chirped-pulse Fourier transform microwave spectroscopy in a supersonic jet expansion, recording the rotational spectrum in the 2-8 GHz frequency region. Two different conformers were detected in the spectrum. The most intense transitions correspond to a skew (gauche-gauche) conformation, identified as the global minimum. The spectra of ten different isotopologues were assigned for this species, leading to accurate effective and substitution structures. The weaker spectrum presents small tunnelling doublings caused by the torsional motion of the thiol group, which are only compatible with an antiperiplanar skeleton and a gauche thiol. The larger stability of the global minimum is attributed to an intramolecular S-H⋅⋅⋅π weak hydrogen bond. A comparison of the intramolecular interactions in the title molecule and 2-phenylethanol, similarly stabilized by a O-H⋅⋅⋅π hydrogen bond, shows the different strength of these interactions. Density functional (B3LYP-D3, B2PLYP-D3) and ab initio (MP2) calculations were conducted for the molecule.

Molecular mechanisms underlying nanowire formation in pristine phthalocyanine

Understanding the molecular processes of nanowire self-assembly is crucial for designing and controlling nanoscale structures that could lead to breakthroughs in functional materials. In this work, we focus on pristine phthalocyanines as a representative example of mesogenic supramolecular assemblies and have analyzed the formation of nanowires using classical molecular dynamics simulations. In the simulations, the molecules spontaneously form multi-columnar structures resembling supramolecular polymers that subsequently grow into more ordered aggregates. These self-assemblies are concentration dependent, leading to the formation of multi-columnar, dynamic aggregates at higher concentrations and nanowires at lower concentrations. The multi-columnar assemblies on a whole are more disordered than the nanowires, but have locally ordered domains of parallel facing molecules that can fluctuate while maintaining their overall shape. The nanowire formation at lower concentrations involves the initial interaction and clustering of randomly oriented phthalocyanine molecules, followed by the merging of small clusters into elongated segments and the eventual formation of a stable nanowire. We observe three main conformers in these self-assemblies, the parallel, T-shaped and edge-to-edge stacking of the phthalocyanine dimers. We calculate the underlying free energy landscape and show that the parallel conformers form the most stable configuration which is followed by the T-shaped and edge-to-edge dimer configurations. The findings provide insights into the mechanisms and pathways of nanowire formation and a step towards the understanding of self-assembly processes in supramolecular mesogens.

The Energetic Origins of Pi-Pi Contacts in Proteins

Accurate potential energy models of proteins must describe the many different types of noncovalent interactions that contribute to a protein's stability and structure. Pi-pi contacts are ubiquitous structural motifs in all proteins, occurring between aromatic and nonaromatic residues and play a nontrivial role in protein folding and in the formation of biomolecular condensates. Guided by a geometric criterion for isolating pi-pi contacts from classical molecular dynamics simulations of proteins, we use quantum mechanical energy decomposition analysis to determine the molecular interactions that stabilize different pi-pi contact motifs. We find that neutral pi-pi interactions in proteins are dominated by Pauli repulsion and London dispersion rather than repulsive quadrupole electrostatics, which is central to the textbook Hunter-Sanders model. This results in a notable lack of variability in the interaction profiles of neutral pi-pi contacts even with extreme changes in the dielectric medium, explaining the prevalence of pi-stacked arrangements in and between proteins. We also find interactions involving pi-containing anions and cations to be extremely malleable, interacting like neutral pi-pi contacts in polar media and like typical ion-pi interactions in nonpolar environments. Like-charged pairs such as arginine-arginine contacts are particularly sensitive to the polarity of their immediate surroundings and exhibit canonical pi-pi stacking behavior only if the interaction is mediated by environmental effects, such as aqueous solvation.

Deciphering the non-covalent interactions in the furan⋯hexane complex using rotational spectroscopy and theoretical analyses

The rotational spectrum of a binary complex formed between furan and n-hexane was investigated using a chirped pulse Fourier transform microwave spectrometer in the range of 2-6 GHz. While furan has only one conformer, n-hexane exists in multiple conformations. The conformational landscape of the binary complex was systematically explored by using a semiempirical conformational search tool, namely CREST. The CREST conformational candidates were subjected to further geometry optimization and harmonic frequency calculations at the B3LYP-D3BJ/def2-TZVP level of theory, resulting in 34 minima within an energy window of 5 kJ mol-1. The three most stable furan⋯hexane minima all contain the most stable n-hexane conformer subunit and are separated by relatively low conformational conversion barriers. Additional calculations were carried out to support the conclusive identification of the global minimum structure responsible for the set of assigned rotational transitions. These include calculations at the B3LYP-D3BJ level with the aug-cc-pVTZ and 6-311++G(d,p) basis sets and the MP2/def2-TZVP level, as well as the single point energy calculations at the CCSD(T)-F12/cc-pVDZ level. Further non-covalent interaction and principal interacting orbital analyses show that the synergy of the πfuran → σ*hexane and σhexane → π*furan interactions plays an important role in stabilizing the observed furan-hexane conformer.