Mutual influence of parallel, CH/O, OH/π and lone pair/π interactions in water/benzene/water system

First-Principles Characterization of Noncovalent Interactions in the Enantioselective Oxetane Ring Opening via SPHENOL-Based Chiral Phosphoric Acid Catalyst

The asymmetric ring-opening reaction of 3,3-substituted oxetanes catalyzed by chiral phosphoric acid (CPA) derived from a newly developed SPHENOL (2,2',3,3'-tetrahydro-1,1'-spirobi[phenalene]-9,9'-diol) framework was investigated by performing density functional calculations. The plausible transition states were traced by considering the oxetane activation mode. The energy profiles obtained for various substrates provide a rational understanding of the reaction at the atomic level. The truncation models help to attribute the reaction initiation and origin of enantioselectivity to three types of noncovalent interactions between the catalyst and the substrate. The Quantum theory of atoms in molecules (QTAIM), Noncovalent Interactions (NCI) Plots, and Wiberg Bond Index (WBI) provide conclusive evidence of the origin of stereoselectivity for the intramolecular desymmetrization of the five substrates considered here. The current study establishes that the SPHENOL-derived CPA catalyst forges enantioselective desymmetrization of 3,3-substituted oxetanes with higher enantiomeric excess in comparison to previously known SPINOL- and BINOL-derived CPA catalysts.

Atmospheric microplastics and nanoplastics as vectors of primary air pollutants - A theoretical study on the polyethylene terephthalate (PET) case

Polyethylene terephthalate (PET) microplastics and nanoplastics are ubiquitously present in the atmosphere as atmospheric and airborne forms (PET-aMPs). Using first-principles calculations, we analyze the uptake of primary air pollutants onto PET-aMPs, focusing on their stabilities, adsorption mechanisms, and thermochemistry. The results show that PET-aMPs are selective for the spontaneous adsorption of CO, CO₂, NO, N₂O, NO₂, NH_3, and SO₂, reaching stable adsorption energies of 6-20 kcal/mol per molecule, with comparable uptake ability than carbon-based materials, metals/metalloids, and metal oxide surfaces. Then, PET-aMPs become a vector for coexisting air pollutants in the atmosphere, which adsorb by inner or outer adsorption depending on the molecular polarity (dipole moment) and atomic constitution (electronegativity) of gaseous molecules. Also, atmospheric H₂O and O₂ are not competitive molecules, and ozone could enhance adsorption due to surface oxidation and structure breakdown. The interplay of electrostatic (46-61%) and dispersion forces (21-58%) drives the adsorption mechanism, where low-polar pollutants display almost a balanced electrostatic vs. dispersion contribution, while high polar molecules display a higher electrostatic stabilization. The outer adsorption is reached by strong dispersion, hydrogen bonding, and dipole-dipole-induced pairs, while lone-pair-π interactions appear in the inner adsorption regime. These results expand the understanding of the hazards and risks of atmospheric and airborne microplastics/nanoplastics, their impacts, co-transport ability, and interaction with the environment.

Comparative oxidative ability of mononuclear and dinuclear high-valent iron-oxo species towards the activation of methane: does the axial/bridge atom modulate the reactivity?

Over the years, mononuclear Fe^IVO species have been extensively studied, but the presence of dinuclear Fe^IVO species in soluble methane monooxygenase (sMMO) has inspired the development of biomimic models that could activate inert substrates such as methane. There are some successful attempts; particularly the [(Por)(m-CBA) Fe^IV(μ-N)Fe^IV(O)(Por˙⁺)]^- species has been reported to activate methane and yield decent catalytic turnover numbers and therefore regarded as the closest to the sMMO enzyme functional model, as no mononuclear Fe^IVO analogues could achieve this feat. In this work, we have studied a series of mono and dinuclear models using DFT and ab initio DLPNO-CCSD(T) calculations to probe the importance of nuclearity in enhancing the reactivity. We have probed the catalytic activities of four complexes: [(HO)Fe^IV(O)(Por)]^- (1), [(HO)Fe^IV(O)(Por˙⁺)] (2), μ-oxo dinuclear iron species [(Por)(m-CBA)Fe^IV(μ-O)Fe^IV(O) (Por˙⁺)]^- (3) and N-bridged dinuclear iron species [(Por)(m-CBA)Fe^IV(μ-N)Fe^IV(O)(Por˙⁺)]^- (4) towards the activation of methane. Additionally, calculations were performed on the mononuclear models [(X)Fe^IV(O)(Por˙⁺)]ⁿ {X = N 4a (n = -2), NH 4b (n = -1) and NH₂4c (n = 0)} to understand the role of nuclearity in the reactivity. DFT calculations performed on species 1-4 suggest an interesting variation among them, with species 1-3 possessing an intermediate spin (S = 1) as a ground state and species 4 possessing a high-spin (S = 2) as a ground state. Furthermore, the two Fe^IV centres in species 3 and 4 are antiferromagnetically coupled, yielding a singlet state with a distinct difference in their electronic structure. On the other hand, species 2 exhibits a ferromagnetic coupling between the Fe^IV and the Por˙⁺ moiety. Our calculations suggest that the higher barriers for the C-H bond activation of methane and the rebound step for species 1 and 3 are very high in energy, rendering them unreactive towards methane, while species 2 and 4 have lower barriers, suggesting their reactivity towards methane. Studies on the system reveal that model 4a has multiple FeN bonds facilitating greater reactivity, whereas the other two models have longer Fe-N bonds and less radical character with steeper barriers. Strong electronic cooperativity is found to be facilitated by the bridging nitride atom, and this cooperativity is suppressed by substituents such as oxygen, rendering them inactive. Thus, our study unravels that apart from enhancing the nuclearity, bridging atoms that facilitate strong cooperation between the metals are required to activate very inert substrates such as methane, and our results are broadly in agreement with earlier experimental findings.

Association Complexes of Calix[6]arenes with Amino Acids Explained by Energy-Partitioning Methods

Intermolecular complexes with calixarenes are intriguing because of multiple possibilities of noncovalent binding for both polar and nonpolar molecules, including docking in the calixarene cavity. In this contribution calix[6]arenes interacting with amino acids are studied with an additional aim to show that tools such as symmetry-adapted perturbation theory (SAPT), functional-group SAPT (F-SAPT), and systematic molecular fragmentation (SMF) methods may provide explanations for different numbers of noncovalent bonds and of their varying strength for various calixarene conformers and guest molecules. The partitioning of the interaction energy provides an easy way to identify hydrogen bonds, including those with unconventional hydrogen acceptors, as well as other noncovalent bonds, and to find repulsive destabilizing interactions between functional groups. Various other features can be explained by energy partitioning, such as the red shift of an IR stretching frequency for some hydroxy groups, which arises from their attraction to the phenyl ring of calixarene. Pairs of hydrogen bonds and other noncovalent bonds of similar magnitude found by F-SAPT explain an increase in the stability of both inclusion and outer complexes.

Switching a HO···π Interaction to a Nonconventional OH···π Hydrogen Bond: A Completed Crystallographic Puzzle

In this article, we present crystallographic and spectroscopic evidence of a tunable system wherein a HO···π interaction switches incrementally to a nonconventional OH···π hydrogen bonding (HB) interaction. More specifically, we report the synthesis of substituted forms of model system 1 to study the effects of aryl ring electronic density on the qualitative characteristics of OH···π hydrogen bonds therein. The OH stretch in experimental infrared data, in agreement with density-functional theory (DFT) calculations, shows continuous red-shifts as the adjacent ring becomes more electron rich. For example, the OH stretch of an amino-substituted analogue is red-shifted by roughly 50 cm^-1 compared to the same stretch in the CF₃ analogue, indicating a significantly stronger HB interaction in the former. Moreover, DFT calculations (ωB97XD/6-311+G**) predict that increasing electronic density on the adjacent top ring reduces the aryl π-OH σ* energy gap with a concomitant enhancement of the OH n-π* energy gap. Consequently, a dominant π-σ* interaction in the amino substituted analogue locks the system in the in-form while a favorable n-π* interaction reverses the orientation of the oxygen-bound hydrogen in its protonated form. Additionally, the ¹H NMR data of various analogues reveal that stronger OH···π interactions in systems with electron-rich aromatic rings slow exchange of the alcoholic proton, thereby revealing coupling with the geminal proton. Finally, X-ray crystallographic analyses of a spectrum of analogues clearly visualize the three distinct stages of "switch"-starting with exclusive HO···π, to partitioned HO···π/OH···π, and finally to achieving exclusive OH···π forms. Furthermore, the crystal structure of the amino analogue reveals an interesting feature in which an extended HB network, involving two conventional (NH···O) and two nonconventional (OH···π) HBs, dimerizes and anchors the molecule in the unit cell.

Probing the ionization potentials of the formaldehyde dimer

In this work, we present a computational investigation on the ionization potentials (IPs) of the formaldehyde dimer, (H₂CO)₂. Twelve lowest lying IPs (corresponding to the entire valence orbitals) for both C_2h and C_s symmetry conformers have been computed at the coupled cluster level of theory using large correlation consistent basis sets with extrapolation to the complete basis set limit and consideration of core electron correlation effects. Specifically, the equation-of-motion ionization potential coupled-cluster with single and double (EOMIP-CCSD) excitations method with the aug-cc-pVXZ and aug-cc-pCVXZ (X = T, Q, and 5) basis sets combined with the Feller-Peterson-Dixon approach was employed, as well as CCSD with perturbative triples [CCSD(T)] with the aug-cc-pVTZ basis sets. In general, excellent agreement was observed from the comparison between the results obtained through the use of these approaches. In addition, the IPs for the formaldehyde monomer were also obtained using such methodologies and the results compared with existing experimental data; excellent agreement was also observed in this case. To the best of our knowledge, this work represents the first of its kind to determine the IPs for all these systems using a high level theory approach and is presented to motivate experimental investigations, e.g., studies involving photoionization, particularly for the formaldehyde dimer. The equilibrium binding energy of the C_2h dimer is calculated in this work at the CCSD(T)/aug-cc-pVTZ level of theory to be -4.71 kcal/mol. At this same level of theory, the equilibrium isomerization energy between C_2h and C_s conformers is 0.76 kcal/mol (C_s conformer being more stable).

Strong CH/O interactions between polycyclic aromatic hydrocarbons and water: Influence of aromatic system size

Energies of CH/O interactions between water molecule and polycyclic aromatic hydrocarbons with a different number of aromatic rings were calculated using ab initio calculations at MP2/cc-PVTZ level. Results show that an additional aromatic ring in structure of polycyclic aromatic hydrocarbons significantly strengthens CH/O interactions. Calculated interaction energies in optimized structures of the most stable tetracene/water complex is -2.27 kcal/mol, anthracene/water is -2.13 kcal/mol and naphthalene/water is -1.97 kcal/mol. These interactions are stronger than CH/O contacts in benzene/water complex (-1.44 kcal/mol) while CH/O contacts in tetracene/water complex are even stronger than CH/O contacts in pyridine/water complexes (-2.21 kcal/mol). Electrostatic potential maps for different polycyclic aromatic hydrocarbons were calculated and used to explain trends in the energies of interactions.

Non-covalent interaction of benzene with methanol and diethyl ether solid surfaces

We have investigated the interactions involved at the interface of binary, layered ices (benzene on methanol and on diethyl ether) by means of laboratory experiments and ab initio calculations on model clusters.

Water sandwiched by a pair of aromatic rings in a proton-conducting metal–organic framework

Water–benzene interactions are investigated in a metal–organic framework incorporating strong H-bond nets which contribute to distinct proton conduction.

Serine and Cysteine π-Interactions in Nature: A Comparison of the Frequency, Structure, and Stability of Contacts Involving Oxygen and Sulfur

Despite many DNA–protein π-interactions in high-resolution crystal structures, only four X–H···π or X···π interactions were found between serine (Ser) or cysteine (Cys) and DNA nucleobase π-systems in over 100 DNA–protein complexes (where X = O for Ser and X = S for Cys). Nevertheless, 126 non-covalent contacts occur between Ser or Cys and the aromatic amino acids in many binding arrangements within proteins. Furthermore, Ser and Cys protein–protein π-interactions occur with similar frequencies and strengths. Most importantly, due to the great stability that can be provided to biological macromolecules (up to –20 kJ mol–1 for neutral π-systems or –40 kJ mol–1 for cationic π-systems), Ser and Cys π-interactions should be considered when analyzing protein stability and function.

What are preferred water-aromatic interactions in proteins and crystal structures of small molecules?

The distribution of water molecules around aromatic rings in the protein structures and crystal structures of small molecules shows quite a small number of the strongest OH-π interactions, a larger number of parallel interactions, and the largest number of the weakest CH-O interactions.

Parallel water/aromatic interactions of non-coordinated and coordinated water

The parallel interactions of non-coordinated and coordinated water molecules with an aromatic ring were studied by analyzing data in the Cambridge structural database (CSD) and by using quantum chemical calculations. The CSD data show that water/aromatic contacts prefer parallel to OH/π interactions, which indicates the importance of parallel interactions. The results reveal the influence of water coordination to a metal ion; the interactions of aqua complexes are stronger. Coordinated water molecules prefer a parallel-down orientation in which one OH bond is parallel to the aromatic ring, whereas the other OH bond points to the plane of the ring. The interactions of aqua complexes with parallel-down water/benzene orientation are as strong as the much better known OH/π orientations. The strongest calculated interaction energy is -14.89 kcal mol(-1) . The large number of parallel contacts in crystal structures and the quite strong interactions indicate the importance of parallel orientation in water/benzene interactions.

C–H/O interactions of nucleic bases with a water molecule: a crystallographic and quantum chemical study

The C–H/O interactions of nucleic bases are substantially stronger than the C–H/O interactions of benzene and pyridine. These results can be very important for molecular recognition of DNA and RNA.