Identifying and Analyzing Intermolecular Bonding Interactions in van der Waals Molecules

Substituent effect on inter-ring interaction in paracyclophanes

The theoretical calculations, namely multipole-derived charge analysis, quantum theory of atom in molecules, and non-bonding interaction (NCI), were performed for [2.2]paracyclophanes, [2.2]paracyclophane-7,9-dienes, and [3.3]paracyclophanes optimized at B3LYP/6-311++G** level, including dispersion correction. The substituent effect of the electron donor N(Me)2 and electron acceptor NO2 group and the influence of the length of bridges joining the aromatic ring on aromatic ring interaction energy (AIE) and strain energy were discussed. The local and electrostatic character of the substituent effect in paracyclophanes was shown. The presence of the weak orbital through-space C···C interaction between the [3.3]paracyclophane ring and weak CH···O hydrogen bonds between the substituents in the different rings was shown.

3-D bond-paths of QTAIM and the stress tensor in neutral lithium clusters, Lim (m = 2-5), presented on the Ehrenfest force molecular graph

In this investigation we set out to understand the origins of non-nuclear attractors (NNAs) found for neutral lithium clusters Lim (m = 2-5) on the QTAIM molecular graph but not on the Ehrenfest force F(r) molecular graph. Therefore, we pursued the stress tensor σ(r) without using the dependency on the QTAIM partitioning, since previously σ(r) was only calculated within the QTAIM partitioning, to see if any indication of NNA character can be determined. Because the stress tensor σ(r) lacks an associated scalar- or vector-field as is the case for QTAIM and the Ehrenfest F(r) partitioning schemes respectively, a stress tensor σ(r) partitioning scheme cannot be constructed. Therefore, to overcome this difficulty we use next generation QTAIM, constructed from the most preferred directions of electronic charge density accumulation, to calculate the stress tensor σ(r) 3-D bond-paths on the Ehrenfest force F(r) molecular graph. Using next generation 3-D bond-paths within the Ehrenfest force F(r) partitioning, we can classify the degree of NNA character in the absence of NNAs. A much higher degree of NNA character is found to be present for the stress tensor σ(r) 3-D bond-paths than for the corresponding QTAIM or Ehrenfest force F(r) 3-D bond-paths. The stabilizing effect of the NNA is demonstrated by undertaking Li2 bond-path compression and stretching distortions sufficient to cause the annihilation of the NNA. The compression and stretching distortions also lead to a large increase in the 3-D bond-path asymmetry and persistent bond-path torsion respectively.

Strong proton-shared hydrogen bonding in a methyl imidazole⋯HCl complex: evidence from matrix isolation infrared spectroscopy and ab initio computations

Evidence for proton-shared hydrogen bonding is provided in a methyl imidazole⋯HCl complex using matrix isolation infrared spectroscopy and ab initio computations.

Prototypical cyclohexane dimers: spectroscopic evidence for σ stacking at low temperatures

In this work, the first unambiguous spectroscopic evidence for the existence of σ stacking interactions in cyclohexane dimers has been provided using matrix isolation infrared spectroscopy.

Stress Tensor Eigenvector Following with Next-Generation Quantum Theory of Atoms in Molecules: Excited State Photochemical Reaction Path from Benzene to Benzvalene

In this investigation, we considered both the scalar and 3-D vector-based measures of bonding using next generation quantum theory of atoms in molecules (QTAIM), constructed from the preferred direction of electronic charge density accumulation, to better understand the photochemical reaction associated with of the formation of benzvalene from benzene. The formation of benzvalene from benzene resulted in two additional C-C bonds forming compared with the benzene. The creation of the additional C-C bonds was explained in terms of an increasing the favorability of the reaction process by maximizing the bonding density. The topological instability of the benzvalene structure was determined using the scalar and vector-based measures to explain the short chemical half-life of benzvalene in terms of the competition between the formation of unstable new C-C bonding that also destabilizes nearest neighbor C-C bonds. The explosive character of benzvalene is indicated by the unusual tendency of the C-C bonds to rupture as easily as weak bonding. The topological instability of the short strong C-C bonds was explained by the existence of measures from conventional and next generation QTAIM that previously have only been observed in weak interactions; such measures included twisted 3-D bonding descriptors.

Enhancement of Atmospheric Nucleation by Highly Oxygenated Organic Molecules: A Density Functional Theory Study

New particle formation (NPF) by gas-particle conversion is the main source of atmospheric aerosols. Highly oxygenated organic molecules (HOMs) and sulfuric acid (SA) are important NPF participants. 2-Methylglyceric acid (MGA), a kind of HOMs, is a tracer of isoprene-derived secondary organic aerosols. The nucleation mechanisms of MGA with SA were studied using density functional theory and atmospheric cluster dynamics simulation in this study, along with that of MGA with methanesulfonic acid (MSA) as a comparison. Our theoretical works indicate that the (MGA)(SA) and (MGA)(MSA) clusters are the most stable ones in the (MGA) _i(SA) _j ( i = 1-2, j = 1-2) and (MGA) _i(MSA) _j ( i = 1-2, j = 1-2) clusters, respectively. Both the formation rates of (MGA)(SA) and (MGA)(MSA) clusters are quite large and could have significant contributions to NPF. The results imply that the homomolecular nucleation of MGA is unlikely to occur in the atmosphere, and MGA and SA can effectively contribute to heteromolecular nucleation mainly in the form of heterodimers. MSA exhibits properties similar to SA in its ability to form clusters with MGA but is slightly weaker than SA.

On the Uselessness of Bond Paths Linking Distant Atoms and on the Violation of the Concept of Privileged Exchange Channels

We refer to frequently used determinants suggesting dominant interactions between distant atoms in various dimers. First of all, we show, against the still-prevailling opinion, that, in general, bond paths have nothing in common with dominant intermolecular interactions and therefore they are useless in such cases. Quite the contrary, reliable information about dominant intermolecular interactions can be obtained by means of electrostatic potential maps, which very convincingly explain mutual orientation of molecules in a dimer. For the first time, numerous examples of interactions that violate both the concept of privileged exchange channels proposed by Pendás and his collaborators as well as inequalities obtained by Tognetti and Joubert for the β parameter related to secondary interactions are presented. The possible cause of this violation is suggested. We also show that the so-called counterintuitive bond paths result from quite natural behavior of the electron density gradient vector, i. e. searching for those areas of space that are characterized by large values of electron density or the most expanded its distributions.

Topological analysis and quantum mechanical structure of Hn+ clusters II. H5+ and H7+

Full typologocal and chemical features of H5+ and H7+ have been calculated at RHF/6-311G(3P) using AIM theory. The results have been compared with previous work on H3+.

Elusive hypervalent phosphorus⋯π interactions: evidence for paradigm transformation from hydrogen to phosphorus bonding at low temperatures

A paradigm transformation from hydrogen to phosphorus bonding is found to depend on the proton affinity of the interacting π-systems.

Comparison of pnicogen and tetrel bonds in complexes containing CX2 carbenes (X = F, Cl, Br, OH, OMe, NH2, and NMe2)

The similarities and differences of pnicogen and tetrel bonds formed by carbenes CX₂ with H₃AsO and H₃SiCN were investigated by carrying out ab initio calculations in association with topological analysis of electron density.

Computer-Aided Design of Molecularly Imprinted Polymers for Simultaneous Detection of Clenbuterol and Its Metabolites

Molecular imprinting technology (MIT) offers an effective technique for efficient separation and enrichment of specific analytes from complicated matrices and has been used for illicit veterinary drug detectionin recent years due to its high selectivity, good chemical stability, and simple preparation. The development of in silico-based approaches has enabled the simulation of molecularly imprinted polymers (MIPs) to facilitate the selection of imprinting conditions such as template, functional monomer, and the best suitable solvent. In this work, using density functional theory (DFT), the molecularly imprinted polymers of clenbuterol and its metabolites were designed by computer-aided at B3LYP/6-31 + G (d, p) level. Screening molecular imprinting components such as functional monomers, cross-linkers, and solvents has been achieved in the computational simulation considerations. The simulation results showed that methacrylic acid (MAA) is the best functional monomer; the optimal imprinting ratio for both clenbuterol (CLB) and its dummy template molecule of phenylephrine (PE) to functional monomer is 1:3, while the optimal imprinting ratio for the two dummy template molecules of CLB's metabolites is 1:5. Choosin gethyleneglycol dimethacrylate (EDGMA) as a crosslinker and aprotic solvents could increase the selectivity of the molecularly imprinted system. Atoms in Molecules (AIM) topology analysis was applied to investigate the template-monomer complexes bonding situation and helped to explain the nature of the reaction in the imprinting process. These theoretical predictions were also verified by the experimental results and found to be in good agreement with the computational results. The computer-simulated imprinting process compensates for the lack of clarity in the mechanism of the molecular imprinting process, and provides an important reference and direction for developing better recognition pattern towards CLB and its metabolite analytes in swine urine samples at the same time.

A theoretical study of energy transfer in Ar(1S) + SO2( X ̃ 1 A ') collisions: Cross sections and rate coefficients for vibrational transitions

Vibrational transitions, induced by collisions between rare-gas atoms and molecules, play a key role in many problems of interest in physics and chemistry. A theoretical investigation of the translation-to-vibration (T-V) energy transfer process in argon atom and sulfur dioxide molecule collisions is presented here. For such a purpose, the framework of the quasi-classical trajectory (QCT) methodology was followed over the range of translational energies 2 ≤ E_tr/kcal mol^-1 ≤ 100. A new realistic potential energy surface (PES) for the ArSO₂ system was developed using pairwise addition for the four-body energy term within the double many-body expansion. The topological features of the obtained function are compared with a previous one reported by Hippler et al. [J. Phys. Chem. 90, 6158 (1986)]. To test the accuracy of the PES, additional coupled cluster singles and doubles method with a perturbative contribution of connected triples calculations were carried out for the global minimum configuration. From dynamical calculations, the cross sections for the T-V excitation process indicate a barrier-type mechanism due to strong repulsive interactions between SO₂ molecules and the Ar atom. Corrections to zero-point energy leakage in QCT were carried out using vibrational energy quantum mechanical threshold of the complex and variations. Rate coefficients and cross sections are calculated for some vibrational transitions using pseudo-quantization approaches of the vibrational energy of products. Main attributes of the title molecular collision are discussed and compared with available information in the literature.

Bond paths between distant atoms do not necessarily indicate dominant interactions

The goal of the article is to revive discussion on the interpretation of bond paths linking distant atoms, particularly tracing weak interactions in dimers. According to the Pendás' concept of privileged exchange channel, a bond path is formed between this pair of competing atoms, which is associated with larger value of the exchange energy. We point out that, due to the short-range nature of the exchange energy, bond paths linking distant atoms clearly become doubtful indicators of dominant intermolecular interactions, particularly if some other characteristics (geometric, spectroscopic, based on electrostatic parameters, etc.) indicate other intermolecular interactions as dominant. Several such cases are thoroughly investigated. We show that electrostatic parameters are much more reliable indicators of dominant intermolecular interactions than bond paths. Then, we pay attention that the presence of ("unexpected", i.e., not necessarily indicating dominant intermolecular interactions) bond paths between pairs of atoms featuring highly expanded charge distributions can be easily explained by visual exploration of isodensity contour plots. As always pointing in the direction of the steepest increase, the gradient vector of the electron density favors areas of its high values gaining higher exchange energy, yet being blind to highly electron deficient areas nearby, which, however, can quite often be involved in dominant intermolecular interactions as strongly suggested by many other bonding analysis. We also suggest that an interatomic component of Hellmann-Feynman force would most likely be the most reliable indicator of attractive or repulsive character of individual interatomic interaction. © 2018 Wiley Periodicals, Inc.

Donor-Acceptor Complexes between 1,2,5-Chalcogenadiazoles (Te, Se, S) and the Pseudohalides CN- and XCN- (X=O, S, Se, Te)

Donor-acceptor (D-A) complexes between 3,4-dicyano-1,2,5-chalcogenadiazoles [chalcogen=Te (1 a), Se (1 b), S (1 c)] and the pseudohalides CN^- and XCN^- (X=O, S, Se, Te) were studied experimentally and theoretically. For 1 a, they were isolated as [K(18-crown-6)][1 a-CN] (2), [K(18-crown-6)][1 a-NCO] (3), [K(18-crown-6)][1 a-SCN] (4), [K(18-crown-6)][1 a-SeCN] (5), and [K][1 a-NCSe] (6) and characterized by X-ray diffraction (XRD), UV/Vis and NMR spectroscopy, and DFT and QTAIM calculations. For 1 b and 1 c, the complexes were not isolated due to unfavorable thermodynamics. In all isolated complexes, the D-A bonds, stabilized by negative hyperconjugation, were longer than the sum of the covalent radii and shorter than the sum of the van der Waals radii of the bonded atoms. In mixtures of 1 a, F^- , and SeCN^- , the complex [1 a-F]^- was selectively formed in accordance with thermodynamics. The reaction of 1 a with SeCN^- and the cyclic trimeric perfluoro-ortho-phenylene mercury afforded the complex [K(18-crown-6)][SCN]⋅(o-C₆ F₄ Hg)₃ , which was characterized by XRD.

Hydrogen bond docking preference in furans: OH⋯π vs. OH⋯O

The docking sites of hydrogen bonds in complexes formed between 2,2,2-trifluoroethanol (TFE), furan (Fu), and 2-methyl furan (MF) have been investigated. Using density functional theory (DFT) calculations, gas phase and matrix isolation FTIR spectroscopies, the strengths of OH⋯O and OH⋯π hydrogen bonds in the complexes were compared to find the docking preference. Calculations suggest that the hydrogen bond donor, TFE, is more likely to dock onto the oxygen atom of the aromatic furans ring, and consequently, the OH⋯O type hydrogen bond is relatively stronger than the OH⋯π type. The FTIR spectrum in the OH-stretching fundamental range obtained at room temperatures has been compared with that obtained at extremely low temperatures in the matrix. The fundamental and the red shifts of OH-stretching vibrations were observed in both FTIR spectra, confirming the formation of hydrogen bonded complexes. By assessing the ability of furan and MF to participate in the formation of OH⋯O hydrogen bond, the effect of ring methylation has been highlighted. From the calculated geometric and thermodynamic parameters as well as the frequency shift of the OH-stretching vibrations in complexes, TFE-MF is found to be more stable than TFE-Fu, which suggests that the strength of the OH⋯O hydrogen bond in TFE-MF originates from the high activity of the furan molecule caused by the methylation of the aromatic ring. The present study furthers the knowledge of docking preference in heteroaromatic molecules and is helpful to understand the nature of intermolecular interactions between hydrogen bond donors and acceptors, including both electron-deficient atoms and π cloud.

Hydration of 3-hydroxy-4,4-dimethylglutaric acid with dimethylamine complex and its atmospheric implications

Changes in temperature affects the distribution of isomers, which facilitates the understanding of new particle formation in the atmosphere.

Tetracyclo(9-methyl-2,7-carbazole) as a promising nanohoop for gas trapping: a multiscale study

H₂S, CS₂, NO₂, Br₂, HF, and C₂H₆ are the ideal adsorbates within the TCC host from their respective congeners.

Hydration of the methanesulfonate–ammonia/amine complex and its atmospheric implications

Methanesulfonate (MSA⁻), found in substantial concentrations in the atmosphere, is expected to enhance aerosol nucleation and the growth of nanoparticles, but the details of methanesulfonate clusters are poorly understood. In this study, MSA⁻ was chosen along with ammonia (NH₃) or three common amines and water (H₂O) to discuss the roles of ternary homogeneous nucleation and ion-induced nucleation in aerosol formation. We studied the structural characteristics and thermodynamics of the clusters using density functional theory at the PW91PW91/6-311++G(3df,3pd) level. The analysis of noncovalent interactions predicts that the amines can form more stable clusters with MSA⁻ than NH₃, in agreement with the results from structures and thermodynamics; however, the enhancement in stability for amines is not large enough to overcome the difference in the concentrations of NH₃ and amines under typical atmospheric conditions. In addition, the favorable free energies of formation for the (MSA⁻)(NH₃/amines)(H₂O)_n (n = 0–3) clusters at 298.15 K show that MSA⁻ could contribute to the aerosol nucleation process with binding NH₃/amines and H₂O up to n = 3. There are strong temperature and humidity dependences for the formation of complexes; higher humidity and temperature promote the formation of larger hydrates. Finally, for the (MSA⁻)(NH₃/amines)(H₂O)_n clusters, the evaporation rates were determined to further investigate the atmospheric implications.

Methanesulfonate (MSA⁻), found in substantial concentrations in the atmosphere, is expected to enhance aerosol nucleation and the growth of nanoparticles, but the details of methanesulfonate clusters are poorly understood.