A B3LYP and MP2(full) theoretical investigation into the strength of the C–NO2 bond upon the formation of the intermolecular hydrogen-bonding interaction between HF and the nitro group of nitrotriazole or its methyl derivatives

The changes of bond dissociation energy (BDE) in the C–NO₂ bond and nitro group charge upon the formation of the intermolecular hydrogen-bonding interaction between HF and the nitro group of 14 kinds of nitrotriazoles or methyl derivatives were investigated using the B3LYP and MP2(full) methods with the 6-311++G**, 6-311++G(2df,2p) and aug-cc-pVTZ basis sets. The strength of the C–NO₂ bond was enhanced and the charge of nitro group turned more negative in complex in comparison with those in isolated nitrotriazole molecule. The increment of the C–NO₂ bond dissociation energies correlated well with the intermolecular H-bonding interaction energies. Electron density shifts analyses showed that the electron density shifted toward the C–NO₂ bond upon complex formation, leading to the strengthened C–NO₂ bond and the possibly reduced explosive sensitivity.

A comparative theoretical investigation into the strength of the trigger-bond in the Na+, Mg2+ and HF complexes involving the nitro group of R–NO2 (R = –CH3, –NH2 and –OCH3) or the C = C bond of (E)-O2N–CH = CH–NO2

A comparative theoretical investigation into the change in strength of the trigger-bond upon formation of the Na⁺, Mg²⁺ and HF complexes involving the nitro group of RNO₂ (R?=? –CH₃, –NH₂, –OCH₃) or the C?=?C bond of (E)-O₂N–CH?=?CH–NO₂ was carried out using the B3LYP and MP2(full) methods with the 6-311++G**, 6-311++G(2df,2p) and aug-cc-pVTZ basis sets. Except for the Mg²⁺?π system with (E)-O₂N–CH?=?CH–NO₂ (i.e., C₂H₂N₂O₄?Mg²⁺), the strength of the trigger-bond X–NO₂ (X?=?C, N or O) was enhanced upon complex formation. Furthermore, the increment of bond dissociation energy of the X–NO₂ bond in the Na⁺ complex was far greater than that in the corresponding HF system. Thus, the explosive sensitivity in the former might be lower than that in the latter. For C₂H₂N₂O₄?Mg²⁺, the explosive sensitivity might also be reduced. Therefore, it is possible that introducing cations into the structure of explosives might be more efficacious at reducing explosive sensitivity than the formation of an intermolecular hydrogen-bonded complex. AIM, NBO and electron density shifts analyses showed that the electron density shifted toward the X–NO₂ bond upon complex formation, leading to a strengthened X–NO₂ bond and possibly reduced explosive sensitivity.

A B3LYP and MP2(full) theoretical investigation into explosive sensitivity upon the formation of the molecule-cation interaction between the nitro group of 3,4-dinitropyrazole and H+, Li+, Na+, Be2+ or Mg2+

The explosive sensitivity upon the formation of molecule-cation interaction between the nitro group of 3,4-dinitropyrazole (DNP) and H⁺, Li⁺, Na⁺, Be²⁺ or Mg²⁺ has been investigated using the B3LYP and MP2(full) methods with the 6-311++G** and 6-311++G(2df,2p) basis sets. The bond dissociation energy (BDE) of the C3–N7 trigger bond has also been discussed for the DNP monomer and the corresponding complex. The interaction between the oxygen atom of nitro group and H⁺ in DNP…H⁺ is partly covalent in nature. The molecule-cation interaction and bond dissociation energy of the C3–N7 trigger bond follow the order of DNP…Be²⁺ > DNP…Mg²⁺ > DNP…Li⁺ > DNP…Na⁺. Except for DNP…H⁺, the increment of the trigger bond dissociation energy in comparison with the DNP monomer correlates well with the molecule-cation interaction energy, natural charge of the nitro group, electron density ρ _BCP(C3–N7), delocalization energy E ⁽²⁾ and NBO charge transfer. The analyses of atoms in molecules (AIM), natural bond orbital (NBO) and electron density shifts have shown that the electron density of the nitro group shifts toward the C3–N7 trigger bond upon the formation of the molecule-cation interaction. Thus, the trigger bond is strengthened and the sensitivity of DNP is reduced.     

Theoretical investigation on the structure and thermodynamic properties of the 2,4-dinitroimidazole complex with methanol

The structure and thermodynamic properties of the 2, 4-dinitroimidazole complex with methanol were investigated using the B3LYP and MP2(full) methods with the 6-31++G(2d,p) and 6-311++G(3df,2p) basis sets. Four types of hydrogen bonds [N–H?O, C–H?O, O–H?O (nitro oxygen) and O–H?π] were found. The hydrogen-bonded complex having the highest binding energy had a N–H?O hydrogen bond. Analyses of natural bond orbital (NBO) and atoms-in-molecules (AIM) revealed the nature of the intermolecular hydrogen-binding interaction. The changes in thermodynamic properties from monomers to complexes with temperatures ranging from 200.0 to 800.0 K were investigated using the statistical thermodynamic method. Hydrogen-bonded complexes of 2,4-dinitroimidazole with methanol are fostered by low temperatures.

A B3LYP and MP2(full) theoretical investigation into the cooperativity effect between dihydrogen-bonding and H–M∙∙∙π (M = Li, Na, K) interactions among HF, MH with the π-electron donor C2H2, C2H4 or C6H6

The DFT-B3LYP/6-311++G(3df,2p) and MP2(full)/6-311++G(3df,2p) calculations were carried out on the binary complex formed by HM (M?=?Li, Na, K) and HF or the π-electron donor (C₂H₂, C₂H₄, C₆H₆), as well as the ternary system FH???HM???C₂H₂/C₂H₄/C₆H₆. The cooperativity effect between the dihydrogen-bonding and H–M???π interactions was investigated. The result shows that the equilibrium distances R _H???H and R _M???π in the ternary complex decrease and both the H???H and H–M???π interactions are strengthened when compared to the corresponding binary complex. The cooperativity effect of the dihydrogen bond on the H–M???π interaction is more pronounced than that of the M???π bond on the H???H interaction. Furthermore, the values of cooperativity effect follow the order of FH???HNa???π?>?FH???HLi???π?>?FH???HK???π and FH???HM???C₆H₆?>?FH???HM???C₂H₄?>?FH???HM???C₂H₂. The nature of the cooperativity effect was revealed by the analyses of the charge of the hydrogen atoms in H???H moiety, atom in molecule (AIM) and electron density shifts methods.

Density functional theory investigation of electrophilic addition reaction of bromine to tricyclo[4.2.2.22,5]dodeca-1,5-diene1

The geometry and the electronic structure of tricyclo[4.2.2.2^2,5]dodeca-1,5-diene (TCDD) molecule were investigated by DFT/B3LYP and /B3PW91 methods using the 6-311G(d,p) and 6-311++G(d,p) basis sets. The double bonds of TCDD molecule are syn-pyramidalized. The structure of π-orbitals and their mutual interactions for TCDD molecule were investigated. Potential energy surface (PES) of the TCDD-Br₂ system was studied by B3LYP/6-311++G(d,p) method and the configurations [molecular charge-transfer (CT) complex, transition states (TS1 and TS2), intermediate (INT) and product (P)] corresponding to the stationary points (minima or saddle points) were determined. Initially, a molecular CT-complex forms between Br₂ and TCDD. With a barrier of 22.336 kcal mol^-1 the CT-complex can be activated to an intermediate (INT) with energy 15.154 kcal mol^-1 higher than that of the CT-complex. The intermediate (INT) then transforms easily (barrier 5.442 kcal mol^-1) into the final, N-type product. The total bromination is slightly exothermic. Accompanying the breaking of Br-Br bond, C1-Br, C5-Br and C2-C6 bonds are formed, and C1 = C2 and C5 = C6 double bonds transform into single bonds. The direction of the reaction is determined by the direction of intramolecular skeletal rearrangement that is realized by the formation of C2-C6 bond.

Penta- and heteropentadienyl ligands coordinated to beryllium

In this work we have performed a systematic study of new organometallic complexes containing penta- and heteropentadienyl (CH₂CHCHCHX, X?=?CH₂, O, NH, S) ligands coordinated to beryllium. Calculated complexes were studied using the density functional theory (PBE) in combination with the 6-311++G(3d,2p) basis set. The coordination number on the beryllium atom varies according to the type of ligand. Pentadienyl ligand shows hapticities η¹ and η⁵, while heteropentadienyl ligands display η¹ and η² hapticities. A Wiberg bond indices study was performed in order to get information about their bond orders.

Structural characterization of the (MeSH)4 potential energy surface

A random walk on the PES for (MeSH)₄ clusters produced 50 structural isomers held together by hydrogen-bonding networks according to calculations performed at the B3LYP/6–311++G** and MP2/6–311++G** levels. The geometric motifs observed are somewhat similar to those encountered for the methanol tetramer, but the interactions responsible for cluster stabilization are quite different in origin. Cluster stabilization is not related to the number of hydrogen bonds. Two distinct, well-defined types of hydrogen bonds scattered over a wide range of distances are predicted.

Theoretical insights into the structures and mechanical properties of HMX/NQ cocrystal explosives and their complexes,and the influence of molecular ratios on their bonding energies

Molecular dynamics (MD) methods were employed to study the binding energies and mechanical properties of selected crystal planes of 1,3,5,7-tetranitro-1,3,5,7-tetrazacyclooctane (HMX)/nitroguanidine (NQ) cocrystals at different molecular molar ratios. The densities and detonation velocities of the cocrystals at different molar ratios were estimated. The intermolecular interaction and bond dissociation energy (BDE) of the N–NO₂ bond in the HMX:NQ (1:1) complex were calculated using the B3LYP, MP2(full) and M06-2X methods with the 6-311++G(d,p) and 6-311++G(2df,2p) basis sets. The results indicated that the HMX/NQ cocrystal prefers cocrystalizing in a 1:1 molar ratio, and the cocrystallization is dominated by the (0 2 0) and (1 0 0) facets. The K, G, and E values of the ratio of 1:1 are smaller than those of the other ratios, and the 1:1 cocrystal has the best ductility. The N–NO₂ bond becomes stronger upon the formation of the intermolecular H-bonding interaction and the sensitivity of HMX decreases in the cocrystal. This sensitivity change in the HMX/NQ cocrystal originates not only from the formation of the intermolecular interaction but also from the increment of the BDE of N–NO₂ bond in comparison with isolated HMX. The HMX/NQ (1:1) cocrystal exhibits good detonation performance. Reduced density gradient (RDG) reveals the nature of cocrystallization. Analysis of the surface electrostatic potential further confirmed that the sensitivity decreases in complex (or cocrystal) in comparison with that in isolated HMX.

Open image in new window

Graphical Abstract Binding energies and mechanical properties of HMX/NQ cocrystals in different molecular molar ratios were studied using molecular dynamics methods. The origin of the sensitivity change in the HMX/NQ cocrystal originates from formation of intermolecular interactions and the bond dissociation energy increment of the N–NO₂ bond

     

Non-covalent interactions – QTAIM and NBO analysis

MP2(full)/6-311++G(3df,3pd) calculations were carried out on complexes linked through various non-covalent Lewis acid – Lewis base interactions. These are: hydrogen bond, dihydrogen bond, hydride bond and halogen bond. The quantum theory of ´atoms in molecules´ (QTAIM) as well as the natural bond orbitals (NBO) method were applied to analyze properties of these interactions. It was found that for the A-H…B hydrogen bond as well as for the A-X…B halogen bond (X designates halogen) the complex formation leads to the increase of s-character in the A-atom hybrid orbital aimed toward the H or X atom. In opposite, for the A…H-B hydride bond, where the H-atom possesses negative charge, the decrease of s-character in the B-atom orbital is observed. All these changes connected with the redistribution of the electron charge being the effect of the complex formation are in line with Bent´s rule. The numerous correlations between energetic, geometrical, NBO and QTAIM parameters were also found.

Modulating weak intramolecular interactions through the formation of beryllium bonds: complexes between squaric acid and BeH2

The electronic structure of the two most stable isomers of squaric acid and their complexes with BeH₂ were investigated at the B3LYP/6-311?+?G(3df,2p)// B3LYP/6-31?+?G(d,p) level of theory. Squaric acid forms rather strong beryllium bonds with BeH₂, with binding energies of the order of 60 kJ?mol^?1. The preferential sites for BeH₂ attachment are the carbonyl oxygen atoms, but the global minima of the potential energy surfaces of both EZ and ZZ isomers are extra-stabilized through the formation of a BeH···HO dihydrogen bond. More importantly, analysis of the electron density of these complexes shows the existence of significant cooperative effects between the beryllium bond and the dihydrogen bond, with both becoming significantly reinforced. The charge transfer involved in the formation of the beryllium bond induces a significant electron density redistribution within the squaric acid subunit, affecting not only the carbonyl group interacting with the BeH₂ moiety but significantly increasing the electron delocalization within the four membered ring. Accordingly the intrinsic properties of squaric acid become perturbed, as reflected in its ability to self-associate.

Modeling the effect of H-bonding interactions and molecular packing on the molecular structure of [Ag(ethylnicotinate)2]NO3 complex

The gas phase molecular structure of a single isolated molecule of [Ag(Etnic)₂NO₃];1 where Etnic = Ethylnicotinate was calculated using B3LYP method. The H-bonding interaction between 1 with one (complex 2) and two (complex 3) water molecules together with the dimeric formula [Ag(Etnic)₂NO₃]₂;4 and the tetrameric formula [Ag(Etnic)₂NO₃]₄;5 were calculated using the same level of theory to model the effect of intermolecular interactions and molecular packing on the molecular structure of the titled complex. The H-bond dissociation energies of complexes 2 and 3 were calculated to be in the range of 12.220–14.253 and 30.106–31.055 kcal?mol^?1, respectively, indicating the formation of relatively strong H-bonds between 1 and water molecules. The calculations predict bidentate nitrate ligand in the case of 1 and 2, leading to distorted tetrahedral geometry around the silver ion with longer Ag–O distances in case of 2 compared to 1, while 3 has a unidentate nitrate ligand leading to a distorted trigonal planar geometry. The packing of two [Ag(Etnic)₂NO₃] complex units; 4 does not affect the molecular geometry around Ag(I) ion compared to 1. In the case of 5, the two asymmetric units of the formula [Ag(Etnic)₂NO₃] differ in the bonding mode of the nitrate group, where the geometry around the silver ion is distorted tetrahedral in one unit and trigonal planar in the other. The calculations predicted almost no change in the charge densities at the different atomic sites except at the sites involved in the C–H?O interactions as well as at the coordinated nitrogen of the pyridine ring.

A theoretical investigation into the cooperativity effect between the H∙∙∙O and H∙∙∙F– interactions and electrostatic potential upon 1:2 (F–:N-(Hydroxymethyl)acetamide) ternary-system formation

The cooperativity effects between the O/N–H???F^– anionic hydrogen-bonding and O/N–H???O hydrogen-bonding interactions and electrostatic potentials in the 1:2 (F^–:N-(Hydroxymethyl)acetamide (signed as “ha”)) ternary systems are investigated at the B3LYP/6-311++G** and MP2/6-311++G** levels. A comparison of the cooperativity effect in the “F^–???ha???ha” and “FH???ha^–???ha” systems is also carried out. The result shows that the increase of the H???O interaction energy in the O–H???O–H, N–H???O–H or N–H???O?=?C link is more notable than that in the O–H???O?=?C contact upon ternary-system formation. The cooperativity effect is found in the complex formed by the O/N–H???F^– and O/N–H???O interactions, while the anti-cooperativity effect is present in the system with only the O/N–H???F^– H-bond or the “FH???ha^–???ha” complex by the N^–???H–F contact. Atoms in molecules (AIM) analysis and shift of electron density confirm the existence of cooperativity. The most negative surface electrostatic potential (V _S,min ) correlates well with the interaction energy E′ _{int.(ha???F–)} and synergetic energy E _syn., respectively. The relationship between the change of V _S,min (i.e., ΔV _S,min ) and E _syn. is also found.

Low-energy conformers of pamidronate and their intramolecular hydrogen bonds: a DFT and QTAIM study

Extensive DFT and ab initio calculations were performed to characterize the conformational space of pamidronate, a typical pharmaceutical for bone diseases. Mono-, di- and tri-protic states of molecule, relevant for physiological pH range, were investigated for both canonical and zwitterionic tautomers. Semiempirical PM6 method were used for prescreening of the single bond rotamers followed by geometry optimizations at the B3LYP/6-31++G(d,p) and B3LYP/6-311++G(d,p) levels. For numerous identified low energy conformers the final electronic energies were determined at the MP2/6-311++G(2df,2p) level and corrected for thermal effects at B3LYP level. Solvation effects were also considered via the COSMO and C-PCM implicit models. Reasonable agreement was found between bond lengths and angle values in comparison with X-ray crystal structures. Relative equilibrium populations of different conformers were determined from molecular partition functions and the role of electronic, vibrational and rotational degrees of freedom on the stability of conformers were analyzed. For no level of theory is a zwitterionic structure stable in the gas-phase while solvation makes them available depending on the protonation state. Geometrically identified intramolecular hydrogen bonds were analyzed by QTAIM approach. All conformers exhibit strong inter-phosphonate hydrogen bonds and in most of them the alkyl-amine side chain is folded on the P-C-P backbone for further hydrogen bond formation.

Relationship between structure and entropy contributions in an anthraquinone mercapto derivative

The structural and thermodynamic properties of an anthraquinone derivative were studied by means of quantum-chemical calculations. Conformational analysis using ab initio and density functional theory methods revealed 14 low-energy conformers. In order to discuss similarities and differences in entropy of the conformers, the rotational and vibrational contributions to entropy were correlated with changes in conformer structure. The component of the moment of inertia perpendicular to the molecular plane gives significant input to ΔS _rot , whereas the largest contributions to the ΔS _vib have vibrations associated with the τ _S1C20 coordinate.

Synthetic and quantum chemical study on the regioselective addition of amines to methyl maleamate

Synthetic and theoretical studies were performed to gain insight into the regioselectivity in the mechanism of aspartyl-isoaspartyl formation, modeled by additions of ammonia and primary amines to methyl maleamate. Reactions between maleamate and aliphatic, araliphatic amines or O-methyl acetimidate lead to the formation of N-substituted isoasparaginates. The size of the amine and the activating effect of the amide and ester group on the double bond are the determining factors of the site of addition. The formation of both isomers was observed only in the case of ammonia addition. The regioselectivity was predicted on the basis of the charge distribution for low-energy methyl maleamate conformers, calculated at the B3LYP/6-311++G(2df,2pd)//B3LYP/6-31+G(d) level, both in gas phase and in methanol. The methyl isoasparaginate over methyl asparaginate product ratio was computed based on the free energy Boltzmann distribution of their conformers. The calculated 2 : 1 ratio is in agreement with the experimental regioselectivity of the addition of nitrogen nucleophiles.

A new interaction mechanism of LiNH2 with MgH2: magnesium bond

Quantum chemical calculations were performed for LiNH₂–HMgX (X?=?H, F, Cl, Br, CH₃, OH, and NH₂) complexes to propose a new interaction mechanism between them. This theoretical survey showed that the complexes are stabilized through the combinative interaction of magnesium and lithium bonds. The binding energies are in the range of 63.2–66.5 kcal mol^?1, i.e., much larger than that of the lithium bond. Upon complexation, both Mg–H and Li–N bonds are lengthened. Substituents increase Mg-H bond elongation and at the same time decrease Li-N bond elongation. These cyclic complexes were characterized with the presence of a ring critical point and natural population analysis charges.

DFT studies of the adsorption and dissociation of H2O on the Al13 cluster: origins of this reactivity and the mechanism for H2 release

A theoretical study of the chemisorption and dissociation pathways of water on the Al₁₃ cluster was performed using the hybrid density functional B3LYP method with the 6-311+G(d, p) basis set. The activation energies, reaction enthalpies, and Gibbs free energy of activation for the reaction were determined. Calculations revealed that the H₂O molecule is easily adsorbed onto the Al₁₃ surface, forming adlayers. The dissociation of the first H₂O molecule from the bimolecular H₂O structure via the Grotthuss mechanism is the most kinetically favorable among the five potential pathways for O–H bond breaking. The elimination of H₂ in the reaction of an H₂O molecule with a hydrogen atom on the Al cluster via the Eley–Rideal mechanism has a lower activation barrier than the elimination of H₂ in the reaction of two adsorbed H atoms or the reaction of OH and H. Following the adsorption and dissociation of H₂O, the structure of Al₁₃ is distorted to varying degrees.

Easy methods to study the smart energetic TNT/CL-20 co-crystal

2,4,6,8,10,12-hexanitro-2,4,6,8,10,12-hexaazaisowurtzitane (CL-20) is a high-energy nitramine explosive with high mechanical sensitivity. 2,4,6-trinitrotoluene (TNT) is insensitive but by no means a high performance explosive. To reveal the significant importance and smart-material functionality of the energetic-energetic co-crystals, the stability, mechanical and explosive properties TNT/CL-20 co-crystal, TNT crystal and CL-20 crystal were studied. Non-hydrogen bonded non-covalent interactions govern the structures of energetic-energetic co-crystals. However, it is very difficult to accurately calculate the non-covalent intermolecular interaction energies. In this paper, the local conformation and the intricate non-covalent interactions were effectively mapped and analyzed from the electron density (ρ) and its derivatives. The results show that the two components TNT and CL-20 are connected mainly by nitro–aromatic interactions, and nitro–nitro interactions. The steric interactions in TNT/CL-20 could not be confronted with the attractive interactions. Moreover, the scatter graph of TNT crystal reveals the reason why TNT is brittle. The detailed electrostatic potential analysis predicted that the detonation velocities (D) and impact sensitivity for the compounds both increase in the sequence of CL-20 > TNT/CL-20 co-crystal > TNT. Additionally, TNT/CL-20 co-crystal has better malleability than its pure components. This demonstrates the capacity and the feasibility of realizing explosive smart materials by co-crystallization, even if strong hydrogen bonding schemes are generally lacking in energetic materials.