Van Der Waals Complexes Research Articles

Interplay Between the N→C Dative Bond, Intramolecular Chalcogen Bond and π Conjugation in the Complexes Formed by Cyclo[18]carbon and C14 Polyyne with 1,2,5-Chalcogenadiazoles.

The N→C dative bond (DB), intramolecular chalcogen bond and π conjugation play important roles in determining the structures and properties of some molecular carbon materials and organic/polymeric photovoltaic materials. In this work, the interplay between the N→C DB, intramolecular chalcogen bond and π conjugation in the complexes formed by cyclo[18]carbon and C14 polyyne with 1,2,5-chalcogenadiazoles has been investigated in detail by using reliable quantum chemical calculations. This study has made four main findings. First, only the Te-containing complexes bound by N→C DBs are much more stable than their corresponding van der Waals (vdW) complexes. Second, in addition to through-bond π conjugations, through-space π conjugations also exist in some Se/Te-containing complexes bound by N→C DBs. Third, the cooperativity between intramolecular chalcogen bond, π conjugation between two monomers and N→C DB is not very strong and can be ignored. Fourth, compared to π conjugations, intramolecular Ch⋅⋅⋅C (Ch=O, S, Se, Te) chalcogen bonds play a secondary role in stabilizing the complexes bound by N→C DBs. These findings clearly indicate that the role of "conformational lock", popular in the field of organic optoelectronic materials, may have been greatly overestimated.

Improvements on the FANTASIO set-up and new spectroscopic results concerning the Ar-H2O van der Waals complex in the 2OH overtone region

We present a series of improvements of the FANTASIO experimental set-up that couples a supersonic expansion to a cavity ring-down spectrometer originally reported by M. Herman, et al., Mol. Phys. 105, 815 (2007). In particular, the implementation of a new pulsed slit valve of 80 mm in length is detailed. The design is based on a powerful piezo stack actuator. Schlieren imaging is used to verify the homogeneity of the molecular beam and the quality of the pulsed slit performance. The ability to form weakly bound van der Waals complexes is assessed using the cavity ring-down spectrometer and a new chirped-pulse Fourier transform microwave spectrometer with the 12–18 GHz operating range. The water hexamer and heptamer were observed using this microwave spectrometer. Following enhancements and an extended laser coverage of the CRDS in the 2OH region, we observed and analyzed seven sub-bands of the Ar-H216O van der Waals complex, along with one sub-band of the much less abundant Ar-H218O isotopologue.

Revisiting Numerical Solutions of Weakly Bound Noble Gases' Vibrational Energy Levels Modeled by the Improved Lennard-Jones Potential.

We revisit the numerical solutions of vibrational eigenstates of weakly bound homonuclear and heteronuclear noble gas pairs by applying a Fortran program based on the Numerov method. The harmonic, Lennard-Jones (LJ), Morse, Tang-Toennies (TT), and Improved Lennard-Jones (ILJ) potential models have been implemented to represent the potential energy curves (PECs). The obtained vibrational energies spectrum was tested on the experimental data and accurate ab initio calculations at CCSD(T)/CBS level. The vibrational eigenvalues and eigenfunctions can be reproduced accurately within the ILJ potential model. Moreover, considering the calculated lifetime of van der Waals (vdW) complexes, the implementation of ILJ rather than standard LJ potential model has a significant impact on the systems dynamics by providing more representative atomic trajectories when the function is incorporated in force fields for molecular dynamics (MD) simulations. Overall, the ILJ function is the best suited potential model for the representation of vibrational motions and the determination of vibrational energy levels of weakly bound systems, both at equilibrium and non-equilibrium conditions.

The role of solar photolysis in the atmospheric removal of methacrolein oxide and the methacrolein oxide-water van-der Waals complex in pristine environments.

Biogenic hydrocarbons are emitted into the Earth's atmosphere by terrestrial vegetation as by-products of photosynthesis. Isoprene is one such hydrocarbon and is the second most abundant volatile organic compound emitted into the atmosphere (after methane). Reaction with ozone represents an important atmospheric sink for isoprene removal, forming carbonyl oxides (Criegee intermediates) with extended conjugation. In this manuscript, we argue that the extended conjugation of these Criegee intermediates enables electronic excitation of these compounds, on timescales that are competitive with their slow unimolecular decay and bimolecular chemistry. We show that the complexation of methacrolein oxide with water enhances the absorption cross section of the otherwise dark S1 state, potentially revealing a new avenue for forming lower volatility compounds via tropospherically relevant photochemistry.

The Complex (Organic) Puzzle of the Formation of Hydrogen Cyanide and Isocyanide on Interstellar Ice Analogues.

Aiming to constrain the surface formation of HCN and HNC in the dense interstellar medium on ice-covered dust grains, we investigate the interaction of CN radicals with H2O and CO ices and their subsequent reactivity with H and H2. CN radicals can physisorb on both ices. However, on H2O ice, a hemibond formation is the most common binding mode, while on CO ice, the CN-CO van der Waals complex can form NCCO with a small energy barrier. We show low barrier or barrierless pathways to the formation of HCN and HNC for the reaction H + CN on both ices. Reactivity with H2 involves activation energy barriers to form HCN, which may be overcome by quantum tunneling, while HNC formation is unlikely. The formation of HCN and HNC competes with the formation of NH2CHO on H2O and HCOCN on CO.

Formation of Doubly and Triply Charged Fullerene Dimers in Superfluid Helium Nanodroplets.

Sequential ionization of fullerene cluster ions (C_{60})_{n}^{+} within multiply charged helium nanodroplets leads to the intriguing phenomenon of forming and stabilizing doubly and triply charged fullerene oligomers. While the formation of doubly charged dimers (C_{60})_{2}^{2+} has been predicted in earlier studies, the observation of even triply charged ones (C_{60})_{2}^{3+} is highly surprising. This remarkable resilience against Coulomb explosion is achieved through efficient cooling within the superfluid environment of helium nanodroplets and a sequential ionization scheme that populates covalently bound or physisorbed fullerene dimers. Calculations support the stability of four differently bonded (C_{60})_{2}^{2+} and (C_{60})_{2}^{3+} isomers and predict a low Coulomb barrier (<0.4 eV) preventing even dissociation of cold van der Waals complexes.

Vibrational Relaxation of D2O Induced by Collision with He: A Rigid Bender Close Coupling Study.

The vibrational relaxation of the first excited bending state of D2O induced by collision with He is studied at the close coupling level and using the Rigid Bender approximation. A new 4D potential energy surface is calculated and reported for this system. It is then used to determine the low-lying bound states of the D2O-He van der Waals complex and to perform scattering calculations. Collision rates are determined for pure rotational transitions as well as for rovibrational transitions within the first excited bending state. The results are compared with those obtained for the collision of D2O with other noble gases such as Ne and Ar. We also analyse the differences observed with respect to the H2O+He collisions and compare our results with experiment.

Theoretical Insights into the Vibrational Structure of Carbon Dioxide Rare-Gas Complexes.

Two new flexible-monomer two-body ab initio potential energy surfaces (PESs) for the neon and krypton van der Waals complexes with carbon dioxide were developed, extending our previous work on the Ar-CO2 molecule. The accuracy of the PESs was validated by their agreement with the vibrational spectrum of the rare-gas complexes. The intermolecular and intramolecular vibrational excitation energies were computed at the vibrational self-consistent field and vibrational configuration interaction levels of theory. Overall, the agreement between theory and experiment is excellent throughout the vibrational spectra. The observed slight splitting of the bending modes, resulting from their nondegeneracy in the complexes, is confirmed by our computations, and the results qualitatively agree with the experiment. The splitting increases with increasing polarizability of the rare-gas atom. Additionally, we explain a discrepancy in the mode assignment in the intermolecular region of the neon complex with our VCI character assignment.

Facilitating Electron Transfer by Resizing Cyclocarbon Acceptor from C18 to C16.

Recent advances in synthetic methods, combined with tip-induced on-surface chemistry, have enabled the formation of numerous cyclocarbon molecules. Here, we investigate computationally the experimentally studied C16 and C18 molecules as well as their van der Waals (vdW) complexes with several typical donor and acceptor molecules. Our results demonstrate a remarkable electron-withdrawing ability of cyclocarbon molecules. The vdW complexes of C16 and C18 exhibit a thermodynamically favorable photoinduced electron transfer (ET) from the donor partner to the cyclocarbons that occurs on a picosecond time scale. The lower reorganization energy of C16 compared to C18 leads to a significant acceleration of the ET reactions.

Mechanism of formation of p-benzylenephenol peroxide radical (p-PhC(O2•)HPhOH).

The reactions of radicals with O2 play the important role in the biological, medicinal, and industrial processes. The mechanism of this reaction is studied previously for the alkane, alkene, alkyne, phenols, and close-related radicals. According to these studies, the formation of intermediates in these reactions is predicted only for the aromatic radicals. Thus, the Van der Waals complexes of O2 with phenyl or benzyl radicals are predicted, as well as the π-π cluster for benzene. However, the possibility of the formation of such intermediate π-π clusters in the case of bisphenol radicals and the thermochemistry of its formation is not studied. Bisphenols are one of the main components of bio-oil, produced during pyrolysis of lignin-contained biomass. Synthetic bisphenols are used in polycarbonate plastics, epoxy resins, and thermal papers. Their mechanism of oxidation is important for the determination of fire safety of these materials, the possibility of using them as additives for fuels for the decreasing and the description of the ignition delays, as well as for the determination of its health risk assessment in medicine. The five DFT (M06-2X (i = 1), B3LYP (i = 2), wB97XD (i = 3), M08HX (i = 4), MN15 (i = 5)) approaches with 6-311 + + G(d,p) basis set are used for the determination of standard enthalpies of atomization (ΔraHo(Yi)) of considered compounds (molecules, radicals, and transition states). These values of ΔraHo(Yi) are corrected using the empirical linear calibration dependencies, reported previously. The different calibration dependencies are used for the hydrocarbons (including the aromatics and simple oxygenated derivatives) and for the peroxides. The corrected values of ΔraHo(Yi, CORR) are used according to Hess's law for the determination of ΔfHo(Yi, CORR). The most consistent values of ΔfHo(Y, MEAN) are derived from the coordination of the values of ΔfHo(Yi, CORR) using the intersection of their values of standard deviations (3SDi). These values of ΔfHo(Y, MEAN), as well as the B3LYP values of So(Y), which are accounting the frequency correction and internal rotations, as well as their temperature dependencies, are used for the determination of thermochemistry of considered reactions and of the calculation, within transition state theory (TST), of the values of high-pressure limits of the rate constant. The values of Ho(Yi), So(Yi), and Go(Yi) are calculated using the Gaussian 16w program. The temperature dependencies of thermochemical properties and the values of rate constants are determined using the ChemRate program (v.1.5). The optimized structures are visualized using the Chemcraft.

The Role of Bond Functions in Describing Intermolecular Electron Correlation for Van der Waals Dimers: A Study of (CH4)2 and Ne2.

We present a study of the intermolecular interactions in van der Waals complexes of methane and neon dimers within the framework of the CCSD method. This approach was implemented and applied to calculate and examine the behavior of the contracted two-particle reduced density matrix (2-RDM). It was demonstrated that the region near the minimum of the two-particle density matrix correlation part, corresponding to the primary bulk of the Coulomb hole contribution, exerts a significant influence on the dispersion interaction energetics of the studied systems. As a result, the bond functions approach was applied to improve the convergence performance for the intermolecular correlation energy results with respect to the size of the atomic basis. For this, substantial acceleration was achieved by introducing an auxiliary basis of bond functions centered on the minima of the 2-RDM. For both methane and neon dimers, this general conclusion was confirmed with a series of CCSD calculations for the 2-RDM and the correlation energies.

PES and transport properties of the He⋯HBr complex from kinetic theory and molecular dynamics simulations.

The ab initio intermolecular potential of the He⋯HBr van der Waals (vdW) complex was calculated at the CCSD(T)/a5zBF level of theory and expanded in terms of the orthogonal Legendre polynomials. The PES then was implemented to calculate the interaction viscosity (η12) and diffusion (D12) coefficients through classical Mason-Monchick approximation (MMA), quantum mechanical close-coupling (CC), and molecular dynamics (MD) simulations. Energy-dependent Senftleben-Beenakker (SB) cross-sections were calculated using rotationally averaged cross-sections, and Boltzmann averaging was used to reveal the temperature dependence of the SB cross-sections over the temperature range of T = 50-1000 K. The calculated transport properties from MMA are in close agreement with the CC results, especially for temperatures lower than T = 900 K. The ab initio potential data then were used to derive LJ (12,6) and Vashishta MD force fields, and the equilibrium MD simulation methods were implemented to extract η12 and D12 coefficients using Einstein formulas. It was found that the Vashishta 3-body interaction potential model shows better accuracy than the LJ (12,6) model in MD simulations of D12. For η12, however, both MD potential models are successful, and an average absolute deviation of lower than 1% was obtained when compared to the quantum mechanical CC method.

The temperature variation of the CH+ + H reaction rate coefficients: a puzzle finally understood?

CH+ was the first molecular ion identified in the interstellar medium and is found to be ubiquitous in interstellar clouds. However, its formation and destruction paths are not well understood, especially at low temperatures. A new theoretical approach based on the canonical variational transition state theory was used to study the H + CH+ reactive collisions. Rate coefficients for formation of C+ ions are calculated as a function of temperature. We considered the participation of a direct path and an indirect path in which the reactants should overcome an entropic barrier to form a van der Waals complex or pass through a CH2+ intermediate complex, respectively. We show that the contribution of both pathways to the formation of C+ has to be taken into account. The new reactive rate coefficients for the title reaction, complemented by reactive data for CH+/CH2+ in the H/H2/He mixture, have been used to simulate the corresponding kinetics experimentally measured using an Atomic Beam 22 Pole Trap apparatus at low temperature. A good agreement with the experimental findings was found at 50 K. At a lower temperature, the model overestimates the formation of C+. This shows that secondary reactions are not responsible for the weak C+ production in the experiments at such temperature. Then, we discuss the possible impact of non-adiabatic effects in the study of the H + CH+ reactive collisions and we found that such effects can be responsible for the decrease of the H + CH+ rate coefficients at low temperature. This study offers an explanation for the disagreement between H + CH+ theoretical and experimental rate coefficients which has been going on for 20 years and highlights the need for performing non-adiabatic studies for this simple chemical reaction.

Theoretical dynamics studies of the CH3 + HBr → CH4 + Br reaction: effects of isotope substitution and vibrational excitation of CH3.

The rate coefficient for two deuterium substituted isotopologues of reaction CH3 + HBr → CH4 + Br has been determined using the quasiclassical trajectory (QCT) method. We used the analytical potential energy surface (PES) fitted to high-level ab initio points in earlier work. The PES exhibits a pre-reaction van der Waals complex and a submerged potential barrier. The rate coefficients of the deuterated isotopologue reactions, similarly to the pure-protium isotopologue, show significant deviation from the Arrhenius law, namely, the activation energy is negative below about 600 K and positive above it: k[CH3 + DBr] = 1.35 × 10-11 exp(- 2472/T) + 5.85 × 10-13 exp(335/T) and k[CD3 + HBr] = 2.73 × 10-11 exp(- 2739/T) + 1.46 × 10-12 exp(363/T). The CH3 + DBr reaction is slower by a factor of 1.8, whereas CD3 + HBr isotopologue is faster by a factor of 1.4 compared to the HBr + CH3 system across a wide temperature range. The isotope effects are interpreted in terms of the properties of various regions of the PES. Quantum state-resolved simulations revealed that the reaction of CH3 with HBr becomes slower when any of the vibrational modes of the methyl radical is excited. This contradicts the assumption that vibrational excitation of methyl radicals enhances its reactivity, which is of historical importance: this assumption was used as an argument against the existence of negative activation energy in a decade-long controversy in the 1980s and 1990s.

Room temperature detection of the (H2)2 dimer.

The hydrogen dimer, (H2)2, is among the most weakly bound van der Waals complexes and a prototype species for first principles ab initio studies. The detection of the (H2)2 infrared absorption spectrum was reported more than sixty years ago at a temperature of 20 K. Due to the sharp decrease of the (H2)2 abundance with temperature, detection at room temperature was generally considered hardly achievable. Here we report the first room temperature detection of partly resolved rotational structures of (H2)2 by cavity ring down spectroscopy at sub-atmospheric pressures, in the region of the first overtone band of H2 near 1.2 μm. The quantitative analysis of the absorption features observed around ten allowed or forbidden transition frequencies of the monomer provides insight on the structure of this elusive species and a benchmark for future theoretical studies.

Primary processes in photophysics and photochemistry of a potential light-activated anti-cancer dirhodium complex.

Photophysics and photochemistry of a potential light-activated cytotoxic dirhodium complex [Rh2(µ-O2CCH3)2(bpy)(dppz)](O2CCH3)2, where bpy = 2,2'-bipyridine, dppz = dipyrido[3,2-a:2',3'-c]phenazine (Complex 1 or Rh2) in aqueous solutions was studied by means of stationary photolysis and time-resolved methods in time range from hundreds of femtoseconds to microseconds. According to the literature, Complex 1 demonstrates both oxygen-dependent (due to singlet oxygen formation) and oxygen-independent cytotoxicity. Photoexchange of an acetate ligand to a water molecule was the only observed photochemical reaction, which rate was increased by oxygen removal from solutions. Photoexcitation of Complex 1 results in the formation of the lowest triplet electronic excited state, which lifetime is less than 10ns. This time is too short for diffusion-controlled quenching of the triplet state by dissolved oxygen resulting in 1O2 formation. We proposed that singlet oxygen is produced by photoexcitation of weakly bound van der Waals complexes [Rh2…O2], which are formed in solutions. If this is true, no oxygen-independent light-induced cytotoxicity of Complex 1 exists. Residual cytotoxicity deaerated solutions are caused by the remaining [Rh2…O2] complexes.

Temperature-Dependent Line-Broadening Effects in CO2 Caused by Ar.

This computational study of line-broadening effects is based on an accurate, analytical representation of the intermonomer potential energy surface (PES) of the CO2 ⋅ Ar van der Waals (vdW) complex. The PES is employed to compute collisional broadening coefficients for rovibrational lines of CO2 perturbed by Ar. The semiclassical computations are performed using the modified Robert-Bonamy approach, including real and imaginary terms, and the exact trajectory model. The lines investigated are in the 10001←00011, 01101←00001, 00011←00001, and 00031←00001 vibrational bands and the computations are repeated at multiple temperatures. The computed results are in good agreement with the available experimental values, validating both the intermonomer PES developed and the methodology used. For lines in the 01101←00001 band of CO2 , temperature-dependent Ar-broadening coefficients are reported for the first time. The parameters presented should prove useful, among other applications, for the accurate experimental determination of CO2 and Ar abundances in planetary atmospheres.

Influence of Dispersion Interactions in the Diels-Alder Reaction of Cyclopentadiene, Anthracene and Tetracene with C60: A DFT Examination

Asbtract: The activation and thermodynamic parameters of the Diels-Alder reaction of cyclopentadiene, anthracene and tetracene with C60 fullerene are estimated at the DFT level using PBE0, B3LYP and CAM-B3LYP functionals and the effect of accounting for dispersion interactions in various semi-empirical Grimme models on these parameters. Structural characteristics of the reagents, products, pre-reaction van der Waals complex, and transition state are obtained. Taking into account dispersion corrections have practically no effect on the geometric characteristics of structures with the covalent type of the bond, as well as on the structure of the transition state. Still, it has a noticeable effect only on the structure of the pre-reaction complex, in which the interparticle interaction has a van der Waals character. At the same time, taking into account dispersion corrections plays a significant role in assessing the thermodynamic and activation parameters of the reaction. An analysis of the performed calculations showed that the combination of PBE0/def2-TZVP with the D4 model for dispersion interactions could be recommended for the Diels-Alder reaction considered above.

Potential Energy Surface and Bound States of Ne-Li2+(X2Σg+) van der Waals Complex Based on Ab Initio Calculations.

Theoretical studies of the potential energy surface and vibrational bound states calculations were performed for the ground state of the Ne-Li2+(X2Σg+) van der Waals (vdW) complex. The intermolecular interactions were investigated by using an accurate monoconfigurational RCCSD(T) method and large basis sets (aug-cc-pVnZ, n = T, Q, 5), extrapolated to the complete basis set (CBS) limit. In turn, the obtained raw data from RCCSD(T)/CBS(Q5) calculations were numerically interpolated using the Morse + vdW model and the Reproducing Kernel Hilbert Space (RKHS) polynomial method to generate analytic expressions for the 2D-PES. The RKHS interpolated PES was then used to assess the bound states of the Ne-Li2+(X2Σg+) system through nuclear quantum calculations. By studying the aspect of the potential energy surface, the analysis sheds light on the behavior of the Ne-Li2+(X2Σg+) complex and its interactions between repulsive and attractive forces with other particles. By examining the vibrational states and wave functions of the system, the researchers were able to gain a better understanding of the behavior of the Ne-Li2+(X2Σg+) complex. The calculated radial and angular distributions for all even and odd symmetries are discussed in detail. We observe that the radial distributions exhibit a more complicated nodal structure, representing stretching vibrational behavior in the neon atom along its radial coordinate. For the highest bound states, the situation is very different, and the energies surpass the angular barrier.

Interacting Quantum Atoms and Multipolar Electrostatic Study of XH···π Interactions.

The interaction energies of nine XH···π (X = C, N, and O) benzene-containing van der Waals complexes were analyzed, at the atomic and fragment levels, using QTAIM multipolar electrostatics and the energy partitioning method interacting quantum atoms/fragment (IQA/IQF). These descriptors were paired with the relative energy gradient method, which solidifies the connection between quantum mechanical properties and chemical interpretation. This combination provides a precise understanding, both qualitative and quantitative, of the nature of these interactions, which are ubiquitous in biochemical systems. The formation of the OH···π and NH···π systems is electrostatically driven, with the Qzz component of the quadrupole moment of the benzene carbons interacting with the charges of X and H in XH. There is the unexpectedly intramonomeric role of X-H (X = O, N) where its electrostatic energy helps the formation of the complex and its covalent energy thwarts it. However, the CH···π interaction is governed by exchange-correlation energies, thereby establishing a covalent character, as opposed to the literature's designation as a noncovalent interaction. Moreover, dispersion energy is relevant, statically and in absolute terms, but less relevant compared to other energy components in terms of the formation of the complex. Multipolar electrostatics are similar across all systems.