Identifying and Analyzing Intermolecular Bonding Interactions in van der Waals Molecules

The borazine dimer: the case of a dihydrogen bond competing with a classical hydrogen bond

Dimers of borazine were studied using matrix isolation infrared spectroscopy and ab initio quantum chemical calculations. Computations were performed at the MP2 and M06-2X levels of theory using the 6-311++G(d,p) and aug-cc-pVDZ basis sets for the various homodimers. At both levels of theory, an aligned stacked structure was found to be the global minimum, which was nearly isoenergetic to a parallel displaced structure. A T-shaped structure, where the N-H of one borazine pointed towards the N of the second borazine, was found to be a local minimum. In addition to these structures, a dihydrogen bonded structure, where the hydrogen attached to the nitrogen of borazine interacted with the hydrogen attached to the boron atom of another borazine, was also indicated. Experimentally, we observed the T-shaped dimer and the dihydrogen bonded dimer. This is one of the rare examples of experimental evidence for a dihydrogen bond, in a system other than in a metal hydride. These results for the borazine dimer were clearly different from the benzene dimer where the parallel displaced structure was found to be the global minimum followed by the T-shaped structure at the MP2/aug-cc-pVDZ level of theory. AIM, EDA and NBO analyses were carried out for all the structures to explore the nature of interactions.

Insight into π-hole interactions containing the inorganic heterocyclic compounds S2N2/SN2P2

Similar to σ-hole interactions, the π-hole interaction has attracted much attention in recent years. According to the positive electrostatic potentials above and below the surface of inorganic heterocyclic compounds S₂N₂ and three SN₂P₂ isomers (heterocyclic compounds 1-4), and the negative electrostatic potential outside the X atom of XH₃ (X = N, P, As), S₂N₂/SN₂P₂⋯XH₃ (X = N, P, As) complexes were constructed and optimized at the MP2/aug-cc-pVTZ level. The X atom of XH₃ (X = N, P, As) is almost perpendicular to the ring of the heterocyclic compounds. The π-hole interaction energy becomes greater as the trend goes from 1⋯XH₃ to 4⋯XH₃. These π-hole interactions are weak and belong to "closed-shell" noncovalent interactions. According to the energy decomposition analysis, of the three attractive terms, the dispersion energy contributes more than the electrostatic energy. The polarization effect also plays an important role in the formation of π-hole complexes, with the contrasting phenomena of decreasing electronic density in the π-hole region and increasing electric density outside the X atom of XH₃ (X = N, P, As). Graphical abstract Computed density difference plots for the complexes 3⋯NH ₃ (a ₁), 3⋯PH ₃ (b ₁), 3⋯AsH ₃ (c ₁) and electron density shifts for the complexes 3⋯NH ₃ (a ₂), 3⋯PH ₃ (b ₂),3⋯AsH ₃ (c ₂) on the 0.001 a.u. contour.

QTAIM and Stress Tensor Characterization of Intramolecular Interactions Along Dynamics Trajectories of a Light-Driven Rotary Molecular Motor

A quantum theory of atoms in molecules (QTAIM) and stress tensor analysis was applied to analyze intramolecular interactions influencing the photoisomerization dynamics of a light-driven rotary molecular motor. For selected nonadiabatic molecular dynamics trajectories characterized by markedly different S₁ state lifetimes, the electron densities were obtained using the ensemble density functional theory method. The analysis revealed that torsional motion of the molecular motor blades from the Franck-Condon point to the S₁ energy minimum and the S₁/S₀ conical intersection is controlled by two factors: greater numbers of intramolecular bonds before the hop-time and unusually strongly coupled bonds between the atoms of the rotor and the stator blades. This results in the effective stalling of the progress along the torsional path for an extended period of time. This finding suggests a possibility of chemical tuning of the speed of photoisomerization of molecular motors and related molecular switches by reshaping their molecular backbones to decrease or increase the degree of coupling and numbers of intramolecular bond critical points as revealed by the QTAIM/stress tensor analysis of the electron density. Additionally, the stress tensor scalar and vector analysis was found to provide new methods to follow the trajectories, and from this, new insight was gained into the behavior of the S₁ state in the vicinity of the conical intersection.

Experimental evidence for the blue-shifted hydrogen-bonded complexes of CHF3 with π-electron donors

Blue-shifted hydrogen-bonded complexes of fluoroform (CHF₃) with benzene (C₆H₆) and acetylene (C₂H₂) have been investigated using matrix isolation infrared spectroscopy and ab initio computations. For CHF₃-C₆H₆ complex, calculations performed at the B3LYP and MP2 levels of theory using 6-311++G (d,p) and aug-cc-pVDZ basis sets discerned two minima corresponding to a 1:1 hydrogen-bonded complex. The global minimum correlated to a structure, where the interaction is between the hydrogen of CHF₃ and the π-electrons of C₆H₆ and a weak local minimum was stabilized through H…F interaction. For the CHF₃-C₂H₂ complex, computation performed at MP2/aug-cc-pVDZ level of theory yielded two minima, corresponding to the cyclic C-H…π complex A (global) and a linear C-H…F (n-σ) complex B (local). Experimentally a blue-shift of 32.3cm^-1 and 7.7cm^-1 was observed in the ν₁ C-H stretching mode of CHF₃ sub-molecule in Ar matrix for the 1:1 C-H…π complexes of CHF₃ with C₆H₆ and C₂H₂ respectively. Natural bond orbital (NBO), Atoms-in-molecule (AIM) and energy decomposition (EDA) analyses were carried out to explain the blue-shifting and the nature of the interaction in these complexes.

Theoretical insights into the π-hole interactions in the complexes containing triphosphorus hydride (P3H3) and its derivatives

The π-hole of triphosphorus hydride (P₃H₃) and its derivatives Z₃X₃ (Z = P, As; X = H, F, Cl, Br) was discovered and analyzed. MP2/aug-cc-pVDZ calculations were performed on the π-hole interactions in the HCN...Z₃X₃ complexes and the mutual influence between π-hole interactions and the hydrogen bond in the HCN...HCN...Z₃X₃ and HCN...Z₃X₃...HCN complexes studied. The π-hole interaction belongs to the typical closed-shell noncovalent interaction. The linear relationship was found between the most positive electrostatic potential of the π-hole (V_S,max) and the interaction energy. Moreover, the V_S,max of the π-hole was also found to be linearly correlated to the electrostatic energy term, indicating the important contribution of the electrostatic energy term to the π-hole interaction. There is positive cooperativity between the π-hole interaction and the hydrogen bond in the termolecular complexes. The π-hole interaction has a greater influence on the hydrogen bond than vice versa. The mutual enhancing effect between the π-hole interaction and the hydrogen bond in the HCN...HCN...Z₃X₃ complexes is greater than that in the HCN...Z₃X₃...HCN complexes.

Cu@furfural imine-decorated halloysite as an efficient heterogeneous catalyst for promoting ultrasonic-assisted A3and KA2coupling reactions: a combination of experimental and computational study

A novel catalyst was synthesized by functionalization of HNTs with (3-chloropropyl) trimethoxysilan, thiosemicarbazide and furfural and incorporation of copper species.

Matrix isolation FTIR study of hydrogen-bonded complexes of methanol with heterocyclic organic compounds

Hydrogen bonded complexes of heterocyclic compounds with methanol were studied using matrix isolation FTIR spectroscopy and theoretical calculations.

A Quantum Theory of Atoms-in-Molecules Perspective and DFT Study of Two Natural Products: Trans-Communic Acid and Imbricatolic Acid

The topological features of the charge densities, ρ(r), and the chemical reactivity of two most biologically relevant and chemically interesting scaffold systems i.e. trans-communic acid and imbricatolic acid have been determined using density functional theory. To identify, characterize, and quantify efficiently, the non-covalent interactions of the atoms in the molecules have been investigated quantitatively using Bader's quantum theory of atoms-in-molecules (QTAIM) technique. The bond path is shown to persist for a range of weak H···H as well as C···H internuclear distances (in the range of 2.0–3.0 Å). These interactions exhibit all the hallmarks of a closed-shell weak interaction. To get insights into both systems, chemical reactivity descriptors, such as HOMO–LUMO, ionization potential, and chemical hardness, have been calculated and used to probe the relative stability and chemical reactivity. Some other useful information is also obtained with the help of several other electronic parameters, which are closely related to the chemical reactivity and reaction paths of the products investigated. Trans-communic acid seems to be chemically more sensitive when compared with imbricatolic acid due to its experimentally observed higher half-maximal inhibitory concentration (bioactivity parameter) value, which is in accordance with its higher chemical reactivity as theoretically predicted using density functional theory-based reactivity index. The quantum chemical calculations have also been performed in solution using different solvents, and the relative order of their structural and electronic properties as well as QTAIM-based parameters show patterns similar to those observed in gas phase only. This study further exemplifies the use and successful application of the bond path concept and the quantum theory of atoms-in-molecules.

The protonated 2-halogenated imidazolium cation as the noncovalent interaction donor: the σ-hole and π-hole interactions

The σ-hole and π-hole of the protonated 2-halogenated imidazolium cation (XC₃H₄N₂⁺; X = F, Cl, Br, I) were investigated and analyzed. The monomers of (CH₃)₃SiY(Y=F, Cl, Br, I), considered as the Lewis base, were combined with the σ-hole and π-hole of XC₃H₄N₂⁺ to form the σ-hole and π-hole interactions in the bimolecular complexes (CH₃)₃SiY · · · XC₃H₄N₂⁺ and (CH₃)₃SiY · · · C₃(X)H₄N₂⁺(X/Y=F, Cl, Br, I), respectively. For both the σ-hole and π-hole interactions, the equilibrium geometries of complexes show regular changes according to the sequence of heavy sequence of the noncovalent interaction acceptors and donors. The electrostatic energy is the main contribution in the formation of both kinds of interactions, it has linear relations with the V _S,max values of σ-hole and the V' _S,max values of π-hole. Both the σ-hole and π-hole interactions belong to the closed-shell and noncovalent interactions. The π-hole interactions are stronger than the σ-hole interactions. For the π-hole interactions, the contribution percents of the dispersion energies are somewhat greater than those of the σ-hole interactions, while it is contrary for the polarization energy. Graphical Abstract The protonated 2-halogenated imidazolium cation as the noncovalent interaction donor: the σ-hole and π-hole interactionsᅟ.

A QTAIM and stress tensor investigation of the torsion path of a light-driven fluorene molecular rotary motor

The utility of the QTAIM/stress tensor analysis method for characterizing the photoisomerization of light driven molecular rotary machines is investigated on the example of the torsion path in fluorene molecular motor. The scalar and vector descriptors of QTAIM/stress tensor reveal additional information on the bonding interactions between the rotating units of the motor, which cannot be obtained from the analysis of the ground and excited state potential energy surfaces. The topological features of the fluorene motor molecular graph display that, upon the photoexcitation a certain increase in the torsional stiffness of the rotating bond can be attributed to the increasing topological stability of the rotor carbon atom attached to the rotation axle. The established variations in the torsional stiffness of the rotating bond may cause transfer of certain fraction of the torsional energy to other internal degrees of freedom, such as the pyramidalization distortion. © 2016 Wiley Periodicals, Inc.

Does borazine-water behave like benzene-water? A matrix isolation infrared and ab initio study

Borazine is isoelectronic with benzene and is popularly referred to as inorganic benzene. The study of non-covalent interactions with borazine and comparison with its organic counterpart promises to show interesting similarities and differences. The motivation of the present study of the borazine-water interaction, for the first time, stems from such interesting possibilities. Hydrogen-bonded complexes of borazine and water were studied using matrix isolation infrared spectroscopy and quantum chemical calculations. Computations were performed at M06-2X and MP2 levels of theory using 6-311++G(d,p) and aug-cc-pVDZ basis sets. At both the levels of theory, the complex involving an N-H⋯O interaction, where the N-H of borazine serves as the proton donor to the oxygen of water was found to be the global minimum, in contrast to the benzene-water system, which showed an H-π interaction. The experimentally observed infrared spectra of the complexes corroborated well with our computations for the complex corresponding to the global minimum. In addition to the global minimum, our computations also located two local minima on the borazine-water potential energy surface. Of the two local minima, one corresponded to a structure where the water was the proton donor to the nitrogen of borazine, approaching the borazine ring from above the plane of the ring; a structure that resembled the global minimum in the benzene-water H-π complex. The second local minimum corresponded to an interaction of the oxygen of water with the boron of borazine, which can be termed as the boron bond. Clearly the borazine-water system presents a richer landscape than the benzene-water system.

DFT computations on the hydrogen bonding interactions between methacrylic acid-trimethylolpropane trimethacrylate copolymers and letrozole as drug delivery systems

The hydrogen bonding interactions between letrozole (Let) anticancer drug and three copolymers of methacrylic acid-trimethylolpropane trimethacrylate (M1–M3 as molecular imprinted polymers) were studied using density functional theory (DFT) at both B3LYP and B3PW91 levels. The binding energies were corrected for the basis set superposition error (BSSE) and zero-point vibrational energies (ZPVE) so that the most negative [Formula: see text] were measured for compounds 7 and 8 formed between M1 copolymer and endocyclic N1 and N2 atoms of drug, respectively. Also, among complexes 13–15 in which two copolymers were contributed in the formation of O–H[Formula: see text]N bonds with the drug, compound 13 (containing two M1 copolymers) showed the highest [Formula: see text] value. The interactions of all copolymers with drug were exergonic (spontaneous interaction) and exothermic. The QTAIM data supported the covalent character of the C–N, C–H, N–N, C–O, O–H and O–H[Formula: see text]N bonds, the intermediate nature of C[Formula: see text]N and C[Formula: see text]O bonds while the electrostatic character of C–H[Formula: see text]O, HC[Formula: see text]HC and CH[Formula: see text]N interactions. According to the [Formula: see text], [Formula: see text] and [Formula: see text] values, it was suggested that t complexes 7 and 8 (among two particles systems) as well as complex 13 (among three particles systems) can be the most promising drug delivery systems.

Non-covalent C-Cl…π interaction in acetylene-carbon tetrachloride adducts: Matrix isolation infrared and ab initio computational studies

Non-covalent halogen-bonding interactions between π cloud of acetylene (C2H2) and chlorine atom of carbon tetrachloride (CCl4) have been investigated using matrix isolation infrared spectroscopy and quantum chemical computations. The structure and the energies of the 1:1 C2H2-CCl4 adducts were computed at the B3LYP, MP2 and M05-2X levels of theory using 6-311++G(d,p) basis set. The computations indicated two minima for the 1:1 C2H2-CCl4 adducts; with the C-Cl…π adduct being the global minimum, where π cloud of C2H2 is the electron donor. The second minimum corresponded to a C-H…Cl adduct, in which C2H2 is the proton donor. The interaction energies for the adducts A and B were found to be nearly identical. Experimentally, both C-Cl…π and C-H…Cl adducts were generated in Ar and N2 matrixes and characterized using infrared spectroscopy. This is the first report on halogen bonded adduct, stabilized through C-Cl…π interaction being identified at low temperatures using matrix isolation infrared spectroscopy. Atoms in Molecules (AIM) and Natural Bond Orbital (NBO) analyses were performed to support the experimental results. The structures of 2:1 ((C2H2)2-CCl4) and 1:2 (C2H2-(CCl4)2) multimers and their identification in the low temperature matrixes were also discussed.

Comparative studies on group III σ-hole and π-hole interactions

The σ-hole of M2 H6 (M = Al, Ga, In) and π-hole of MH3 (M = Al, Ga, In) were discovered and analyzed, the bimolecular complexes M2 H6 ···NH3 and MH3 ···N2 P2 F4 (M = Al, Ga, In) were constructed to carry out comparative studies on the group III σ-hole interactions and π-hole interactions. The two types of interactions are all partial-covalent interactions; the π-hole interactions are stronger than σ-hole interactions. The electrostatic energy is the largest contribution for forming the σ-hole and π-hole interaction, the polarization energy is also an important factor to form the M···N interaction. The electrostatic energy contributions to the interaction energy of the σ-hole interactions are somewhat greater than those of the π-hole interactions. However, the polarization contributions for the π-hole interactions are somewhat greater than those for the σ-hole interactions.

Investigating molecular interactions and surface morphology of wax-doped asphaltenes

The nature and origin of bee-like microstructures (bees) in asphalt binders and their impact on asphalt oxidation have been the subject of extensive discussions in recent years.

PCl3–C6H6 heterodimers: evidence for P⋯π phosphorus bonding at low temperatures

A phosphorous trichloride (PCl₃)–benzene (C₆H₆) heterodimer was generated in a low temperature N₂ matrix and was characterized using infrared spectroscopy.

Hydrogen-Bonded Complexes of Phenylacetylene-Acetylene: Who is the Proton Donor?

Hydrogen-bonded complexes of C2H2 and phenylacetylene (PhAc) were studied using matrix isolation infrared spectroscopy and quantum chemical computations. Both C2H2 and PhAc, being potential proton donors, the question arises as to which of the two species would be the proton donor in the PhAc-C2H2 complex; a question that this work primarily addresses. The molecular structures, vibrational frequencies, and interaction energies of the PhAc-C2H2 complexes were calculated at the M06-2X and MP2 levels of theory, employing both 6-311++G(d,p) and aug-cc-pVDZ basis sets. At the M06-2X/aug-cc-pVDZ level, two nearly isoenergetic complexes (BSSE corrected) were indicated to be the global minima; one a C-H···π complex, where C2H2 served as a proton donor to the phenyl π-system in PhAc, and the other a C-H···π complex, where C2H2 served as a proton donor to the acetylene π-system in PhAc. Of the two, only the second complex was identified in the matrix, evidenced by a characteristic large shift in the ≡C-H stretch of C2H2. Experiments were also performed using PhAc deuterated at the acetylene hydrogen (PhAcD) to study the isotopic effects on the vibrational spectra of complexes. The isotopic studies further confirmed the structure of the complex trapped in the matrix, thereby presenting unambiguous evidence that C2H2 served as the proton donor to the acetylene π-system of PhAc. The theory of atoms-in-molecules (AIM), energy decomposition (EDA), and natural bond orbital (NBO) analysis were performed to understand the nature of the interactions involved in the complexes.