Density functional theory study of the interaction between 3‐nitro‐1,2,4‐triazole‐5‐one and ammonia

Two fully optimized geometries of 3‐nitro‐1,2,4‐triazol‐5‐one (NTO)–NH₃ complexes have been obtained with the density function theory (DFT) method at the B3LYP/6‐311++G** level. The intermolecular interaction energy is calculated with zero point energy (ZPE) correction and basis set superposition error (BSSE) correction. The greatest corrected intermolecular interaction of the NTO–NH₃ complexes is ?37.58 kJ/mol. Electrons in complex systems transfer from NH₃ to NTO. The strong hydrogen bonds contribute to the interaction energies dominantly. Natural bond orbital (NBO) analysis is performed to reveal the origin of the interaction. Based on vibrational analysis, the changes of thermodynamic properties from the monomer to complexes with the temperature ranging from 200 K to 800 K have been obtained using the statistical thermodynamic method. It is found that two NTO–NH₃ complexes can be produced spontaneously from NTO and NH₃ at normal temperature. © 2005 Wiley Periodicals, Inc. Int J Quantum Chem, 2005     

Density Functional Theory Studies on Intermolecular Interactions of 4‐Amino‐3,5‐dinitropyrazole with NH3 and H2O

Density functional theory (DFT) method with 6‐311++G** basis set was applied to study intermolecular interactions of 4‐amino‐3,5‐dinitropyrazole (LLM‐116)/NH₃ and LLM‐116/H₂O supermolecules. Four optimized stable supermolecules were found on the potential energy surface. The intermolecular interaction energy was calculated with basis set superposition error (BSSE) correction and zero point energy (ZPE) correction. The greatest corrected intermolecular interaction energies of LLM‐116/NH₃ and LLM‐116/H₂O supermolecules are –42.75 and –19.09 kJ×mol^‐1 respectively, indicating that the intensity of interaction between LLM‐116 and NH₃ is stronger than that of LLM‐116/H₂O. The intermolecular interaction is an exothermic process accompanied by a decrease in the probability of supermolecules formation, and the interactions become weak as temperature increase. Natural bond orbital (NBO) analysis was performed to reveal the origin of interaction. The IR spectra were obtained and assigned by vibrational analysis. Based on vibrational analysis, the changes of thermodynamic properties from LLM‐116 to supermolecules with temperature ranging from 200.0 to 400.0 K were obtained using statistical thermodynamic method.     

Theoretical study on the interactions of halogen‐bonds and pnicogen‐bonds in phosphine derivatives with Br2, BrCl,and BrF

The MP2 ab initio quantum chemistry methods were utilized to study the halogen‐bond and pnicogen‐bond system formed between PH₂X (X = Br, CH₃, OH, CN, NO₂, CF₃) and BrY (Y = Br, Cl, F). Calculated results show that all substituent can form halogen‐bond complexes while part substituent can form pnicogen‐bond complexes. Traditional, chlorine‐shared and ion‐pair halogen‐bonds complexes have been found with the different substituent X and Y. The halogen‐bonds are stronger than the related pnicogen‐bonds. For halogen‐bonds, strongly electronegative substituents which are connected to the Lewis acid can strengthen the bonds and significantly influenced the structures and properties of the compounds. In contrast, the substituents which connected to the Lewis bases can produce opposite effects. The interaction energies of halogen‐bonds are 2.56 to 32.06 kcal·mol^?1; The strongest halogen‐bond was found in the complex of PH₂OH???BrF. The interaction energies of pnicogen‐bonds are in the range 1.20 to 2.28 kcal·mol^?1; the strongest pnicogen‐bond was found in PH₂Br???Br₂ complex. The charge transfer of lp(P) ? σ*(Br? Y), lp(F) ? σ*(Br? P), and lp(Br) ? σ*(X? P) play important roles in the formation of the halogen‐bonds and pnicogen‐bonds, which lead to polarization of the monomers. The polarization caused by the halogen‐bond is more obvious than that by the pnicogen‐bond, resulting in that some halogen‐bonds having little covalent character. The symmetry adapted perturbation theory (SAPT) energy decomposition analysis showes that the halogen‐bond and pnicogen‐bond interactions are predominantly electrostatic and dispersion, respectively.     

Theoretical Study on Hydrogen Bonding of Mono‐ and Dihydrated Complexes of 7‐(3′‐Pyridyl)indole in Excited States

The intermolecular hydrogen bonds of mono‐ and dihydrated complexes of 7‐(3′‐Pyridyl)indole (7‐3′PI) have been investigated using the time‐dependent density functional theory (TD‐DFT) method. The electrostatic potential analysis of monomer 7‐3′PI and 7‐(3′‐Pyridyl)indole‐water (7‐3′PI‐W) indicates that an intermolecular hydrogen bond between two waters can be formed for 7‐(3′‐Pyridyl)indole‐2water (7‐3′PI‐2W) complex. The calculated bond lengths of the intermolecular hydrogen bonds of 7‐3′PI‐W and 7‐3′PI‐2W in the S₁ state (the first excited singlet state) are all shortened compared to the ground state. By the analysis of bond length, charge population and infrared spectra, it is demonstrated that the intermolecular hydrogen bonds of 7‐3′PI‐W and 7‐3′PI‐2W are all strengthened upon electronic excitation to the S1 state. Moreover, the fluorescence of 7‐3′PI‐W and 7‐3′PI‐2W are all red‐shifted to larger wavelength compared to monomer 7‐3′PI. The red‐shift of fluorescence peak of 7‐3′PI‐W and 7‐3′PI‐2W should be attributed to the change of hydrogen bond interaction before and after photoexcitation. Therefore, it can be concluded that the intermolecular hydrogen bonding strengthening in the excited S₁ state induces the fluorescence weakening of 7‐3′PI.     

Polymorphism of 3‐(5‐phenyl‐1,3,4‐oxadiazol‐2‐yl)‐ and 3‐[5‐(pyridin‐4‐yl)‐1,3,4‐oxadiazol‐2‐yl]‐2H‐chromen‐2‐ones

This study of 3‐(5‐phenyl‐1,3,4‐oxadiazol‐2‐yl)‐2H‐chromen‐2‐one, C₁₇H₁₀N₂O₃, 1 , and 3‐[5‐(pyridin‐4‐yl)‐1,3,4‐oxadiazol‐2‐yl]‐2H‐chromen‐2‐one, C₁₆H₉N₃O₃, 2 , was performed on the assumption of the potential anticancer activity of the compounds. Three polymorphic structures for 1 and two polymorphic structures for 2 have been studied thoroughly. The strongest intermolecular interaction is stacking of the `head‐to‐head' type in all the studied crystals. The polymorphic structures of 1 differ with respect to the intermolecular interactions between stacked columns. Two of the polymorphs have a columnar or double columnar type of crystal organization, while the third polymorphic structure can be classified as columnar‐layered. The difference between the two structures of 2 is less pronounced. Both crystals can be considered as having very similar arrangements of neighbouring columns. The formation of polymorphic modifications is caused by a subtle balance of very weak intermolecular interactions and packing differences can be identified only using an analysis based on a study of the pairwise interaction energies.     

A full dimensional grid empowered simulation of the CO2 + CO2 processes

A recently introduced bond–bond formulation of the intermolecular interaction has been extended to six‐atom systems to the end of assembling a new potential energy surface (PES) and has been incorporated into a grid empowered simulator able to handle the modeling of the CO₂ + CO₂ processes. The proposed PES is full dimensional and accounts for the dependence of the intermolecular interaction on some basic physical properties of the colliding partners, including modulations induced by the monomer deformation. The used analytical formulation of the interaction involves a limited number of parameters, each having a clear physical meaning. Guess values for these parameters can also be obtained from analytical correlation formulae. Such estimates can then be fine tuned by exploiting experimental and theoretical information. The resulting PES well describes stretched and bent asymptotic CO₂ monomers as well as the CO₂–CO₂ interaction in the most and less stable configurations. On this potential massive quasiclassical elastic and inelastic detailed scattering trajectories have been integrated, by exploiting the innovative computational technologies of the grid. The efficiency of the approach used and the reliability of the estimates of the dynamical properties obtained in this way is such that we can now plan a systematic evaluation of the state specific rate coefficient matrix elements needed for space craft reentry modeling. Here, we present probabilities and cross sections useful to rationalize some typical mechanisms characterizing the vibrational transitions of the CO₂ + CO₂ system on the flexible monomer proposed PES. On such PES, the key dynamical outcomes are: (a) there is a strong energy interchange between symmetric stretching of the reactants and bending of the products (and viceversa) while asymmetric stretching is strongly adiabatic (b) reactant energy is more efficiently allocated (with respect to the rigid monomers PES) as product vibration when reactant stretching modes are excited while the contrary is true when the reactant bending mode is excited. © 2012 Wiley Periodicals, Inc.     

The Generalized Energy‐Based Fragmentation Approach with an Improved Fragmentation Scheme: Benchmark Results and Illustrative Applications

We propose an improved fragmentation scheme for the generalized energy‐based fragmentation (GEBF) approach, which improves the accuracy of the GEBF approach in total energy calculations and intermolecular interactions. The main modification is to introduce some two‐fragment‐centered primitive subsystems, which are neglected in the previous GEBF implementation. Numerical calculations demonstrate that the present GEBF approach can provide more accurate ground‐state energies and intermolecular interactions. The present GEBF approach with the M06‐2X functional and the cc‐pVTZ basis set are employed to investigate the structures and binding energies in two dimeric species, which are related to pseudopolymorphism of a phenyleneethynylene‐based π‐conjugated molecule. A comparison of the binding free energies in a dimeric species and its corresponding model without C? H???F contacts reveal that the substitution of fluorine atoms weakens the binding of monomers in the dimeric species formed by intermolecular O? H???O hydrogen bonds, but strengthens the binding in the dimer formed by the π–π stacking interaction. Therefore, the C? H???F contacts in these two dimeric species are demonstrated to play a less significant role.     

Ethyl 3‐[1‐(5,5‐dimethyl‐2‐oxo‐1,3,2‐dioxaphosphorin‐2‐yl)propan‐2‐ylidene]carbazate: a combined X‐ray and density functional theory (DFT) study

In the title compound, C₁₁H₂₁N₂O₅P, one of the two carbazate N atoms is involved in the C=N double bond and the H atom of the second N atom is engaged in an intramolecular hydrogen bond with an O atom from the dimethylphosphorin‐2‐yl group, which is in an uncommon cis position with respect to the carbamate group. The cohesion of the crystal structure is also reinforced by weak intermolecular hydrogen bonds. Density functional theory (DFT) calculations at the B3LYP/6‐311++g(2d,2p) level revealed the lowest energy structure to have a Z configuration at the C=N bond, which is consistent with the configuration found in the X‐ray crystal structure, as well as a less stable E counterpart which lies 2.0 kcal mol⁻¹ higher in potential energy. Correlations between the experimental and computational studies are discussed.     

Theoretical Study on the Intermolecular Interactions between 1,1‐Diamino‐2,2‐Dinitroethylene and H2O

The density functional method was applied to the study of 1,1‐diamino‐2,2‐dinitroethylene (Fox‐7)/H₂O dimer. All the possible dimers ( 1, 2 and 3 ), as well as the monomers, were fully optimized with the DFT method at the B3LYP/6‐311++G^** level. The basis set superposition errors (BSSE) are 4.62, 4.07 and 3.45 kJ/mol, and the zero point energy (ZPE) corrections for the interaction energies are 7.94, 5.66 and 6.40 kJ/mol for 1, 2 and 3 , respectively. Dimer 1 is the most stable, judged by binding energy. After BSSE and ZPE corrections, the greatest corrected intermolecular interaction energy of dimer 1 was predicted to be ?29.36 kJ/mol. The charge redistribution mainly occurs on the adjacent N–H··· O atoms and N–O··· H atoms between submolecules. The oxygen in the nitro group acts as a moderate hydrogen acceptor as compared to water oxygen. Based on the statistical thermodynamic method, the standard thermodynamic functions, heat capacities (C⁰_P), entropies (S⁰_T) and thermal corrections to enthalpy (H⁰_T), and the changes of thermodynamic properties on going from monomer to dimer over the temperature range 200.00‐700.00 K were predicted. It is energetically or thermodynamically favorable for Fox‐7 to bind with H₂O and to form dimer 1 at room temperature.     

Halogen‐bonded adduct of 1,2‐dibromo‐1,1,2,2‐tetrafluoroethane and 1,4‐diazabicyclo[2.2.2]octane

Halogen bonding is an intermolecular interaction capable of being used to direct extended structures. Typical halogen‐bonding systems involve a noncovalent interaction between a Lewis base, such as an amine, as an acceptor and a halogen atom of a halofluorocarbon as a donor. Vapour‐phase diffusion of 1,4‐diazabicyclo[2.2.2]octane (DABCO) with 1,2‐dibromotetrafluoroethane results in crystals of the 1:1 adduct, C₂Br₂F₄·C₆H₁₂N₂, which crystallizes as an infinite one‐dimensional polymeric structure linked by intermolecular N...Br halogen bonds [2.829 (3) Å], which are 0.57 Å shorter than the sum of the van der Waals radii.     

The effect of the intermolecular potential formulation on the state‐selected energy exchange rate coefficients in N2–N2 collisions

The rate coefficients for N₂–N₂ collision‐induced vibrational energy exchange (important for the enhancement of several modern innovative technologies) have been computed over a wide range of temperature. Potential energy surfaces based on different formulations of the intramolecular and intermolecular components of the interaction have been used to compute quasiclassically and semiclassically some vibrational to vibrational energy transfer rate coefficients. Related outcomes have been rationalized in terms of state‐to‐state probabilities and cross sections for quasi‐resonant transitions and deexcitations from the first excited vibrational level (for which experimental information are available). On this ground, it has been possible to spot critical differences on the vibrational energy exchange mechanisms supported by the different surfaces (mainly by their intermolecular components) in the low collision energy regime, though still effective for temperatures as high as 10,000 K. It was found, in particular, that the most recently proposed intermolecular potential becomes the most effective in promoting vibrational energy exchange near threshold temperatures and has a behavior opposite to the previously proposed one when varying the coupling of vibration with the other degrees of freedom. © 2014 Wiley Periodicals, Inc.     

Experimental and Theoretical Charge Density Study on Di‐2‐pyrazylamine (Hdpza) Molecule in Crystal

Electron density distribution of Di‐2‐pyrazylamine ( Hdpza ) is studied both by single‐crystal X‐ray diffraction method at 100K and theoretical calculation. Structural determination reveals that Hdpza molecules crystalize in a syn‐anti conformation with an intramolecular C? H?N hydrogen bond between two pyrazine rings and then gather together via two intermolecular N? H?N and C? H?N hydrogen interaction and π? π stacking interaction between pyrazine rings. Charge density analysis is made in terms of deformation density (Δπ), Laplacian distribution and topological analysis of total electron density based on multipole model and theoretical calculation. The agreement between experiment and theory is good. The topological properties at bond critical points of C? C and C? N bonds reveal a covalent bond character, and those of intermolecular interactions, such as hydrogen bonds and π? π stacking interactions, reveal a closed‐shell interaction. The potential energy curve of Hdpza molecule shows that the syn‐anti conformation is the most stable one (global minima) than the other two of syn‐syn and anti‐anti conformations.     

Ab initio Study of Electronic Structure and Properties of Crystalline of 1,5-Diamino-1,2,3,4-tetrazole

Density functional method was applied to study 1,5‐diamino‐1,2,3,4‐tetrazole (DAT, CH₄N₆) in both gaseous and bulk states. The banding and electronic structures of crystalline have been investigated at DFT‐B3LYP/ 6‐311G** level of theory. Relaxed crystal structure compares well with experimental data. The light fluctuation of the frontier orbital, which is mainly formed by atomic orbital of N(4) (heterocycle), is the most reactive part of the molecule, which is in good agreement with the experimental results. The energy gap is 9.035 eV, which indicates that DAT is an insulator. The distribution of electrostatic potential is uniform, indicating DAT is insensitive. The charge density of the intermolecular regions in the plane is not overlaid, indicating that the intermolecular interaction between the neighboring molecules along this direction in the bulk is very weak. The overlap populations of N(1)? N(2) bonds are much less than those of other bonds, therefore the N(1)? N(2) bonds first rupture by external stimuli.     

How the oxazole fragment influences the conformation of the tetraoxazocane ring in a cyclohexanespiro‐3′‐(1,2,4,5,7‐tetraoxazocane): single‐crystal X‐ray and theoretical study

Single crystals of (2S,5R)‐2‐isopropyl‐5‐methyl‐7‐(5‐methylisoxazol‐3‐yl)cyclohexanespiro‐3′‐(1,2,4,5,7‐tetraoxazocane), C₁₆H₂₆N₂O₅, have been studied via X‐ray diffraction. The tetraoxazocane ring adopts a boat–chair conformation in the crystalline state, which is due to intramolecular interactions. Conformational analysis of the tetraoxazocane fragment performed at the B3LYP/6‐31G(d,2p) level of theory showed that there are three minima on the potential energy surface, one of which corresponds to the conformation realized in the solid state, but not to a global minimum. Analysis of the geometry and the topological parameters of the electron density at the (3,?1) bond critical points (BCPs), and the charge transfer in the tetraoxazocane ring indicated that there are stereoelectronic effects in the O—C—O and N—C—O fragments. There is a two‐cross hyperconjugation in the N—C—O fragment between the lone electron pair of the N atom (lpN) and the antibonding orbital of a C—O bond (σ*_C—O) and vice versa between lpO and σ*_C—N. The oxazole substituent has a considerable effect on the geometry and the topological parameters of the electron density at the (3,?1) BCPs of the tetraoxazocane ring. The crystal structure is stabilized via intermolecular C—H…N and C—H…O hydrogen bonds, which is unambiguously confirmed with PIXEL calculations, a quantum theory of atoms in molecules (QTAIM) topological analysis of the electron density at the (3,?1) BCPs and a Hirshfeld analysis of the electrostatic potential. The molecules form zigzag chains in the crystal due to intermolecular C—H…N interactions being electrostatic in origin. The molecules are further stacked due to C—H…O hydrogen bonds. The dispersion component in the total stabilization energy of the crystal lattice is 68.09%.     

Hydrogen‐bonded structures of the isomeric 2‐, 3‐ and 4‐carbamoylpyridinium hydrogen chloranilates

In the three isomeric salts, all C₆H₇N₂O⁺·C₆HCl₂O₄⁻, of chloranilic acid (2,5‐dichloro‐3,6‐dihydroxy‐1,4‐benzoquinone) with 2‐, 3‐ and 4‐carbamoylpyridine, namely, 2‐carbamoylpyridinium hydrogen chloranilate (systematic name: 2‐carbamoylpyridinium 2,5‐dichloro‐4‐hydroxy‐3,6‐dioxocyclohexa‐1,4‐dienolate), (I), 3‐carbamoylpyridinium hydrogen chloranilate, (II), and 4‐carbamoylpyridinium hydrogen chloranilate, (III), acid–base interactions involving H‐atom transfer are observed. The shortest interactions between the cation and the anion in (I) and (II) are pyridinium N—H...(O,O) bifurcated hydrogen bonds, which act as the primary intermolecular interaction in each crystal structure. In (III), an amide N—H...(O,O) bifurcated hydrogen bond, which is much weaker than the bifurcated hydrogen bonds in (I) and (II), connects the cation and the anion.     

State‐specific multireference perturbation theory with improved virtual orbitals: Taming the ground state of F2, Be2, and N2

Adaptation of improved virtual orbitals (IVOs) in state‐specific multireference perturbation theory using Møller–Plesset multipartitioning of the Hamiltonian (IVO‐SSMRPT) is examined in which the IVO‐complete active space configuration interaction (CASCI) is used as an inexpensive alternative to the more involved CAS‐self‐consistent field (CASSCF) orbitals. Unlike the CASSCF approach, IVO‐CASCI does not bear tedious and costly iterations beyond those in the initial SCF calculation. The IVO‐SSMRPT is intruder‐free, and explicitly size‐extensive. In the present preliminary study, the IVO‐SSMRPT method which relies on a small reference space is applied to study potential energy surfaces (PES) of the ground state of challenging, multiconfigurational F₂, Be₂, and N₂ molecules. These systems provide a serious challenge to any ab initio methodology due to the presence of an intricate interplay of nondynamical and dynamical correlations to the entire PES. The quality of the computed PES has been judged by extracting spectroscopic parameters and vibrational levels. The reported results illustrate that the IVO‐SSMRPT method has a potential to yield accuracies as good as the CASSCF‐SSMRPT one with reduced computational labor. Even with small reference spaces, our estimates demonstrate a good agreement with the available experimental values, and some benchmark computations. The blend of accuracy and low computational cost of IVO‐SSMRPT should deserve future attention for the accurate treatment of electronic states of small to large molecular systems for which the wavefunction is characterized by various configurations. © 2015 Wiley Periodicals, Inc.     

Modeling the cyclopolymerization of diallyl ether and methyl α‐[(allyloxy)methyl]acrylate

The free‐radical cyclopolymerization of diallyl ether (1) and methyl α‐(allyloxymethyl)acrylate (2) has been modeled with the B3LYP/6‐31G* methodology, by making use of model compounds for the growing radicals. The cyclization of both monomers is exo, with activation barriers of 5.33 and 9.82 kcal/mol, respectively. To account for the polymerizabilities of these monomers, competing reactions have also been modeled. Although both monomers have a lower barrier for homopolymerization than for cyclization, cyclization dominates due to entropy. This explains the high cyclopolymerization vs. homopolymerization of monomer 2, although its monofunctional counterpart has been reported to homopolymerize well. It has also been shown that the degradative chain transfer by H‐abstraction from the allylic carbon is not effective with this monomer. Poor cyclopolymerization of the monomer 1 has been demonstrated by modeling the degradative chain transfer by H‐abstraction from the allylic carbon, which has been shown to compete very efficiently with polymerization reactions. Additionally, intermolecular propagation reaction has been shown to be facile due to cyclization, since the attacking monomer adopts a cyclic structure. © 2006 Wiley Periodicals, Inc. Int J Quantum Chem, 2007     

Folding‐Induced Modulation of Excited‐State Dynamics in an Oligophenylene–Ethynylene‐Tethered Spiral Perylene Bisimide Aggregate

The excited‐state photophysical behavior of a spiral perylene bisimide (PBI) folda‐octamer ( F8 ) tethered to an oligophenylene–ethynylene scaffold is comprehensively investigated. Solvent‐dependent UV/Vis and fluorescence studies reveal that the degree of folding in this foldamer is extremely sensitive to the solvent, thus giving rise to an extended conformation in CHCl₃ and a folded helical aggregate in methylcyclohexane (MCH). The exciton‐deactivation dynamics are largely governed by the supramolecular structure of F8 . Femtosecond transient absorption (TA) in the near‐infrared region indicates a photoinduced electron‐transfer process from the backbone to the PBI core in the extended conformation, whereas excitation power‐ and polarization‐dependent TA measurements combined with computational modeling showed that excitation energy transfer between the unit PBI chromophores is the major deactivation pathway in the folded counterpart.