Intermolecular interaction characteristics of the Boron-nitrogen rings,and with hexagonal boron nitride: A special molecular adsorption and weak interaction

With the discovery of cyclo[18]carbon, related theoretical research has become a very hot topic. Due to its unique molecular structure, many derivatives have been derived, especially its isoelectronic molecule cyclo B₉N₉. This work focused on the interaction between cyclo B₉N₉ molecules and explored the nature of their interaction. The interaction of the dimer of cyclo B₉N₉ was studied, quantum chemical calculation and wave function analysis were carried out, and compared with cyclo[18]carbon to find its commonness. It showed that the interaction between the cyclo B₉N₉ dimers is a weak electron generating stacking interaction. In addition, the adsorption of cyclo B₉N₉ was studied, and the interaction of graphene analogues hexagonal boron nitride (h-BN). Through calculation and analysis, it was found that the adsorption of cyclo B₉N₉ was much stronger than that of cyclo[18]carbon. In another sense, the π–π stacking interaction of boron nitrogen structure is stronger than that of carbon. To further validate this commonality, B₆N₆ was selected for the same calculation, and similar conclusions were obtained. It was also found that the more Boron nitrogen (BN), the stronger the stacking interaction.     

Ultrastrong Regulation Effect of the Electric Field on the All-Carboatomic Ring Cyclo[18]Carbon**

Cyclo[18]carbon has a very unique geometry and electronic structure. We found that an external electric field (EEF) has an ultrastrong regulation effect on various aspects of the cyclo[18]carbon: (1) The EEF makes the shape of the cyclo[18]carbon change from a circle to an oval, the elongation is particularly striking at a large EEF magnitude. (2) The EEF causes a huge polarization of distribution of in-plane π electrons, and strong EEF can even make some of the electrons detached from the carbon ring (3) EEF significantly lowers LUMO energy and reduces HOMO-LUMO gap (4) Large EEF leads to absorption band in the visible light range and thus makes the cyclo[18]carbon display color (5) Strong EEF causes a large number of new absorption peaks in IR spectrum. We also carefully analyzed how EEF deforms structure of the cyclo[18]carbon from the perspective of atomic forces and decomposition of energy variation, and the reason why the in-plane π electrons are much more polarizable by EEF than the out-of-plane π electrons is discussed. Moreover, we demonstrated that it is feasible to equivalently apply a strong EEF on the cyclo[18]carbon via a purely chemical and thus a more easily achieve way, namely introducing divalent alkaline earth metal cation.     

Model chemistry calculations of thiophene dimer interactions: origin of pi-stacking

The intermolecular interaction energies of thiophene dimers have been calculated by using an aromatic intermolecular interaction (AIMI) model (a model chemistry for the evaluation of intermolecular interactions between aromatic molecules). The CCSD(T) interaction energy at the basis set limit has been estimated from the MP2 interaction energy near the basis set limit and the CCSD(T) correction term obtained by using a medium-size basis set. The calculated interaction energies of the parallel and perpendicular thiophene dimers are -1.71 and -3.12 kcal/mol, respectively. The substantial attractive interaction in the thiophene dimer, even where the molecules are well separated, shows that the major source of attraction is not short-range interactions such as charge transfer but rather long-range interactions such as electrostatic and dispersion. The inclusion of electron correlation increases the attraction significantly. The dispersion interaction is found to be the major source of attraction in the thiophene dimer. The calculated total interaction energy of the thiophene dimer is highly orientation dependent. Although electrostatic interaction is substantially weaker than dispersion interaction, it is highly orientation dependent, and therefore electrostatic interaction play an important role in the orientation dependence of the total interaction energy. The large attractive interaction in the perpendicular dimer is the cause of the preference for the herringbone structure in the crystals of nonsubstituted oligothiophenes (alpha-terthienyls), and the steric repulsion between the beta-substituents is the cause of the pi-stacked structure in the crystals of some beta-substituted oligothiophenes.     

Origin of attraction and directionality of the pi/pi interaction: model chemistry calculations of benzene dimer interaction.

A model chemistry for the evaluation of intermolecular interaction between aromatic molecules (AIMI Model) has been developed. The CCSD(T) interaction energy at the basis set limit has been estimated from the MP2 interaction energy near the basis set limit and the CCSD(T) correction term obtained by using a medium size basis set. The calculated interaction energies of the parallel, T-shaped,and slipped-parallel benzene dimers are -1.48, -2.46, and -2.48 kcal/mol, respectively. The substantial attractive interaction in benzene dimer, even where the molecules are well separated, shows that the major source of attraction is not short-range interactions such as charge-transfer but long-range interactions such as electrostatic and dispersion. The inclusion of electron correlation increases attraction significantly. The dispersion interaction is found to be the major source of attraction in the benzene dimer. The orientation dependence of the dimer interaction is mainly controlled by long-range interactions. Although electrostatic interaction is considerably weaker than dispersion interaction, it is highly orientation dependent. Dispersion and electrostatic interactions are both important for the directionality of the benzene dimer interaction.     

Is Electrostatics Sufficient to Describe Hydrogen‐Bonding Interactions?

The stability and geometry of a hydrogen‐bonded dimer is traditionally attributed mainly to the central moiety A?H???B, and is often discussed only in terms of electrostatic interactions. The influence of substituents and of interactions other than electrostatic ones on the stability and geometry of hydrogen‐bonded complexes has seldom been addressed. An analysis of the interaction energy in the water dimer and several alcohol dimers—performed in the present work by using symmetry‐adapted perturbation theory—shows that the size and shape of substituents strongly influence the stabilization of hydrogen‐bonded complexes. The larger and bulkier the substituents are, the more important the attractive dispersion interaction is, which eventually becomes of the same magnitude as the total stabilization energy. Electrostatics alone are a poor predictor of the hydrogen‐bond stability trends in the sequence of dimers investigated, and in fact, dispersion interactions predict these trends better.     

A theoretical study of the cohesion of noble gases on graphite

The interactions of the noble gases with a graphene sheet are investigated theoretically. The short range repulsive interaction between the noble gas and each carbon atom is described using Hartree-Fock atomic densities and a local density functional theory with the exchange functional corrected for the finite range of the interaction by introducing a Rae-type correction depending on the effective number of electrons. The long range interactions are introduced as the sum of the Axilrod-Teller triple-dipole interaction plus the dipole-dipole and dipole-quadrupole dispersive attractions damped according to the theory of Jacobi and Csanak. The energy arising from the interactions between the permanent quadrupoles on the carbon atoms with the dipole they induce on the noble gas is negligible, being nonzero only on account of the atomistic structure of graphene. The mobile and delocalized nature of the graphene pi electrons causes the effective number of electrons to be around 500 rather than that of 12 appropriate for a system of entirely localized interactions with individual carbon atoms. Inclusion of the Axilrod-Teller term is required to obtain reliable predictions for the binding energies and equilibrium geometries. Absorption of a noble gas atom is predicted to occur at the site above the center of a six membered ring although this is preferred over two other sites by only about 5 meV. The methods presented for generating all the potentials can be applied to derive the interactions between any ion and carbon atom in the wall of a single-walled nanotube. Knowledge of these interactions is required to study the alkali halide nanocrystals encapsulated in single-walled carbon nanotubes of current interest.     

Achieving Adaptive Aromaticity in Cyclo[10]carbon by Screening Cyclo[n]carbon (n=8−24)

Discovery of species with adaptive aromaticity (being aromatic in both the lowest singlet and triplet states) is particularly challenging as cyclic species are generally aromatic either in the ground state or in the excited state only, according to Hückel's and Baird's rules. Inspired by the recent realization of cyclo[18]carbon, here we demonstrate that cyclo[10]carbon possesses adaptive aromaticity by screening cyclo[n]carbon (n=8?24), which is supported by nucleus‐independent chemical shift (NICS), anisotropy of the current‐induced density (ACID), π contribution of electron localization function (ELF_π) and electron density of delocalized bonds (EDDB) analyses. Further study reveals that the lowest triplet state of cyclo[10]carbon is formed by in‐plane ππ* excitation. Thus, the major contribution to the aromaticity from out‐of‐plane π molecular orbitals does not change significantly in the lowest singlet state. Our findings highlight a crucial role of out‐of‐plane π orbitals in maintaining aromaticity for both the lowest singlet and triplet states as well as the aromaticity dependence on the number of the carbon in cyclo[n]carbon.     

Reliable structures and energetics for two new delocalized pi...pi prototypes: cyanogen dimer and diacetylene dimer

Two new prototype delocalized pi[dot dot dot]pi complexes are introduced: the dimers of cyanogen, (N[triple bond]C-C[triple bond]N)(2), and diacetylene, (HC[triple bond]C-C[triple bond]CH)(2). These dimers have properties similar to larger delocalized pi...pi systems such as benzene dimer but are small enough that they can be probed in far greater detail with high accuracy electronic structure methods. Parallel-slipped and T-shaped structures of both cyanogen dimer and diacetylene dimer have been optimized with 15 different procedures. The effects of basis set size, theoretical method, counterpoise correction, and the rigid monomer approximation on the structure and energetics of each dimer have been examined. MP2 and CCSD(T) optimized geometries for all four dimer structures are reported, as well as estimates of the CCSD(T) complete basis set (CBS) interaction energy for every optimized geometry. The data reported here suggest that future optimizations of delocalized pi[dot dot dot]pi clusters should be carried out with basis sets of triple-zeta quality. Larger basis sets and the expensive counterpoise correction to the molecular geometry are not necessary. The rigid monomer approximation has very little effect on structure and energetics of these dimers and may be used without consequence. Due to a consistent cancellation of errors, optimization with the MP2 method leads to CCSD(T)/CBS interaction energies that are within 0.2 kcal mol(-1) of those for structures optimized with the CCSD(T) method. Future studies that aim to resolve structures separated by a few tenths of a kcal mol(-1) should consider the effects of optimization with the CCSD(T) method.     

Vibrational Spectra and Molecular Vibrational Behaviors of All-Carboatomic Rings,cyclo[18]carbon and Its Analogues

The vibrational spectra of cyclo[18]carbon and its analogues, cyclo[2n]carbon (n=3 to 15), were carefully simulated and characterized. The in-plane C−C stretching vibrations shows strong rigidity, while out-of-plane motions seem to be extremely flexible. The solvation effect can enhance signal strengths of the vibrational spectra, but does not evidently change the shape of the spectral curves. The infrared and Raman spectra of cyclo[2n]carbons are quite sensitive to ring size in the range of n=3 to 7, while the size only modestly affects peak positions and strengths for larger rings. Molecular dynamic trajectories show that the fluctuation period of the skeleton of cyclo[18]carbon is basically constant at different temperatures, and they are all about 300 fs. With increase of simulation temperature, the ring distortion due to thermal motion is notable and becomes much stronger. However, neither ring breaking nor isomerization in cyclo[18]carbon is observed during the simulations untill 298.15 K.     

Ab initio calculations of structures and interaction energies of toluene dimers including CCSD(T) level electron correlation correction

The intermolecular interaction energy of the toluene dimer has been calculated with the ARS-F model (a model chemistry for the evaluation of intermolecular interaction energy between ARomatic Systems using Feller's method), which was formerly called as the AIMI model III. The CCSD(T) (coupled cluster calculations with single and double substitutions with noniterative triple excitations) interaction energy at the basis set limit has been estimated from the second-order Moller-Plesset perturbation interaction energy at the basis set limit obtained by Feller's method and the CCSD(T) correction term obtained using a medium-size basis set. The cross (C(2)) dimer has the largest (most negative) interaction energy (-4.08 kcal/mol). The antiparallel (C(2h)) and parallel (C(S)) dimers (-3.77 and -3.41 kcal/mol, respectively) are slightly less stable. The dispersion interaction is found to be the major source of attraction in the toluene dimer. The dispersion interaction mainly determines the relative stability of the stacked three dimers. The electrostatic interaction of the stacked three dimers is repulsive. Although the T-shaped and slipped-parallel benzene dimers are nearly isoenergetic, the stacked toluene dimers are substantially more stable than the T-shaped toluene dimer (-2.62 kcal/mol). The large dispersion interaction in the stacked toluene dimers is the cause of their enhanced stability.     

Theoretical study on mesogenic core structures of nematic liquid crystalline compounds

The structures and intermolecular interaction energies of 10 dimers, included in the mesogenic core structures of typical liquid crystalline (LC) compounds, are obtained at the MP2/6-31G(d) level of theory. It is proved that the dispersion energy significantly contributes to the total interaction energy of these dimers. Even when bulky substituents are introduced into the core part, the interaction energy is still large. It is also revealed that when a long intermolecular distance is provided by a high steric repulsion originating from the linkage of two phenyl groups, the dispersion energy is significantly small. However, in this range of intermolecular distances, the electrostatic energy caused by a strong quadrupole-quadrupole attractive interaction plays a dominant role, and as a result, a rather stable dimer is formed. In all 10 dimers, the dispersion, electrostatic, and exchange-repulsion energies strongly depend on the geometrical orientation of the molecules. The calculated interaction energies of these dimers are also compared with the corresponding experimentally measured viscosities. The results suggest an explicit linear relationship between the interaction energies and viscosities.     

Exchange coupling across the cyanide bridge: structural and DFT interpretation of the magnetic properties of a binuclear chromium(III) complex

The reaction of [Cr(CN)6]3- with a mixture of trans-[Cr(cyclam)(OH)2]Cl, [Cr(cyclam)(OH)Cl]Cl and [Cr(cyclam)Cl2]Cl affords the cyanide bridged dimer, trans-[HO-Cr(cyclam)-NC-Cr(CN)5]-. The tetraphenylphosphonium salt of the anion crystallizes in space group P2(1)/n and shows a bent arrangement of the Cr1-CN-Cr2 unit with the Cr1-CN bond angle at 166.9 degrees and CN-Cr2 at 160.32 degrees . The Cr2-O bond, trans to the hexacyanide fragment, is very short at 1.902 A. Two dimers are held together by two hydrogen bonds connecting the Cr2-OH group of each dimer with one of the NH groups of the cyclam ligand of an adjacent molecule, leading to an almost linear configuration. These dimers of dimers get packed parallel to each other, generating layers separated by the tetraphenylphosphonium cations. Four of the cyanide groups of the anion are engaged in H-bonds with the four water molecules present in the structure or with a NH group of the macrocycle of an adjacent molecule. From magnetic susceptibility measurements, the dimer was found to exhibit antiferromagnetic interaction between the Cr(III) centers with J=-16 cm(-1)(H=-2JS(A)S(B)). Structural and magnetic parameters have been calculated by density functional theoretical methods at the B3LYP level. The exchange coupling constant, J, calculated for the dimer at the X-ray geometry is -23.2 cm(-1) which is in excellent agreement with the experimental value.     

PCILO Calculations on “true” van der Waals complexes

The PCILO method was applied to complexes where dispersion energy forms the dominant attractive contribution of the intermolecular interaction energy. With alkane dimers, especially with dimers of higher alkanes (C₄H₁₀)₂, (C₅H₁₂)₂, etc., and the ethylene dimer PCILO was successful as to stabilization energy, equilibrium distance, optimal dimer structure, and energy composition. The dispersion energy part is always underestimated, however. With (X₂)₂ dimers (X = N, F, Cl) PCILO failed, which might be explained by the presence of lone pairs in the subsystems.     

Origins of cooperative noncovalent host-guest chemistry in mixed valence complexes

The electronic effects resulting from noncovalent host-guest interactions between calix[6]arene and a ruthenium dimer, [Ru3O(OAc)6(CO)(ppy)]2-mu-pz (ppy=4-phenyl pyridine, pz=pyrazine), are presented. The noncovalent interaction is between the calix[6]arene and the ppy ligands of the dimer. The dimer can bind 2 equiv of calix[6]arene. The complex [Ru3O(OAc)6(CO)(ppy)]2-mu-pz forms a highly stable mixed valence ion with strong electronic coupling between the two Ru3 clusters. The strength of the electronic interaction is found to be moderated by calix[6]arene binding. Addition of calix[6]arene to the mixed valence ion causes the electronic coupling to decrease. The binding of calix[6]arene is found to be cooperative. The origins of cooperative binding are developed in terms of the potential energy surfaces associated with the symmetric and asymmetric mixed valence ion. In particular, it is found that symmetry breaking (through the binding of a single calix[6]arene) destabilizes the mixed valence state. Restoration of symmetry (through the binding of a second calix[6]arene) increases the stability of the mixed valence ion and provides an additional driving force for the binding of the second calix[6]arene.     

Assessing the Accuracy of SAPT(DFT) Interaction Energies by Comparison with Experimentally Derived Noble Gas Potentials and Molecular Crystal Lattice Energies

The density functional version of symmetry‐adapted perturbation theory, SAPT(DFT), is a computationally efficient method for calculating intermolecular interaction energies. We evaluate its accuracy by comparison with experimentally determined noble gas interaction potentials and sublimation enthalpies, most of which have not been previously calculated using this method. In order to compare the results with wavefunction methods, we also calculate these quantities using MP2 and, for noble gas dimers, using CCSD(T). For the crystal lattice energy calculations, we include corrections to the dispersion, electrostatic, and induction energies that account for the finite interaction distance cutoff and higher‐order induction contributions. Overall, the energy values extrapolated to the complete basis set limit show that SAPT(DFT) achieves significantly better agreement with experiment than MP2.     

Predicted stability and structure of (HXeCCH)n (n=2 or 4) clusters and of crystalline HXeCCH

Ab initio calculations predict the existence of the dimer and tetramer of HXeCCH. The interaction energies are -6.66 and -19.40 kcal mol-1 for the dimer and tetramer, respectively. For both complexes, larger blue shifts of the Xe-H stretching mode are found, while the Xe-C stretching modes are slightly redshifted. The stability and structure of HXeCCH crystals is predicted by density functional theory calculations with periodic boundary conditions. Strong electrostatic interactions are found between the monomers in the crystal. The results are first evidence for the existence of crystalline materials made of a novel class of noble gas molecules.     

Remarkable Size Effect on Photophysical and Nonlinear Optical Properties of All-Carboatomic Rings,Cyclo[18]carbon and Its Analogues

Inspired by recent experimental observation of molecular morphology and theoretical predictions of multiple properties of cyclo[18]carbon, we systematically studied the photophysical and nonlinear optical properties of cyclo[2N]carbons (N=3–15) allotropes through density functional theory. This work unveils the unusual optical properties of the sp-hybridized carbon rings with different sizes. The remarkable size dependence of the optical properties of these systems and their underlying nature are profoundly explored, and the relevance between aromaticity and optical properties are highlighted. The extrapolation curves fitted for energy level of frontier molecular orbitals, maximum absorption wavelength, and (hyper)polarizability of considered carbon rings are presented, which can be used to reliably predict corresponding properties for arbitrarily large carbon rings. The findings in this study will facilitate the exploration of potential application of cyclocarbons in the field of optical materials.     

Magnitude and orientation dependence of intermolecular interaction of perfluoropropane dimer studied by high-level ab initio calculations: comparison with propane dimer

Intermolecular interaction energies of 12 orientations of C(3)F(8) dimers were calculated with electron correlation correction by the second-order M?ller-Plesset perturbation method. The antiparallel C(2h) dimer has the largest interaction energy (-1.45 kcal/mol). Electron correlation correction increases the attraction considerably. Electrostatic energy is not large. Dispersion is mainly responsible for the attraction. Orientation dependence of the interaction energy of the C(3)F(8) dimer is substantially smaller than that of the C(3)H(8) dimer. The calculated interaction energy of the C(3)F(8) dimer at the potential minimum is 78% of that of the C(3)H(8) dimer (-1.85 kcal/mol), whereas the interaction energies of the CF(4) and C(2)F(6) dimers are larger than those of the CH(4) and C(2)H(6) dimers. The intermolecular separation in the C(3)F(8) dimer at the potential minimum is substantially larger than that in the C(3)H(8) dimer. The larger intermolecular separation due to the steric repulsion between fluorine atoms is the cause of the smaller interaction energy of the C(3)F(8) dimer at the potential minimum. The calculated intermolecular interaction energy potentials of the C(3)F(8) dimers using an all atom model OPLS-AA (OPLS all atom model) force field and a united atom model force field were compared with the ab initio calculations. Although the two force fields well reproduces the experimental vapor and liquid properties of perfluoroalkenes, the comparison shows that the united atom model underestimates the potential depth and orientation dependence of the interaction energy. The potentials obtained by the OPLS-AA force field are close to those obtained by the ab initio calculations.     

Intra‐ and Intermolecular Dispersion Interactions in [n]Cycloparaphenylenes: Do They Influence Their Structural and Electronic Properties?

Cycloparaphenylenes (CPPs) are nanosized structures with unique isolated and bulk properties, and are synthetic targets for the template‐driven bottom‐up synthesis of carbon nanotubes. Thus, a systematic understanding of the supramolecular order at the nanoscale is of utmost relevance for molecular engineering. In this study, it is found that intramolecular noncovalent (dispersion) interactions must be taken into account for obtaining accurate estimates of the structural and optoelectronic properties of [n]CPP compounds, and their influence as the number of repeat units increases from n=4 to n=12 is also analyzed, both in the gas phase and in solution. The supramolecular self‐assembly, for which both intra‐ and intermolecular noncovalent interactions are relevant, of [6]CPP is also investigated by calculating the binding energies of dimers taken along several crystal directions. These are also used to estimate the cohesive energy of the crystal, which is compared to the value obtained by means of dispersion‐corrected DFT calculations using periodic boundary conditions. The reasonable agreement between both computational strategies points towards a first estimate of the [6]CPP cohesive energy of around 50 kcal mol^?1.