Molecular dynamics study on the relationships of modeling, structural and energy properties with sensitivity for RDX-based PBXs

In this paper,a primary model is established for MD(molecular dynamics) simulation for the PBXs(polymer-bonded explosives) with RDX(cyclotrimethylene trinitramine) as base explosive and PS as polymer binder.A series of results from the MD simulation are compared between two PBX models,which are represented by PBX1 and PBX2,respectively,including one PS molecular chain having 46 repeating units and two PS molecular chains with each having 23 repeating units.It has been found that their structural,interaction energy and mechanical properties are basically consistent between the two models.A systematic MD study for the PBX2 is performed under NPT conditions at five different temperatures,i.e.,195 K,245 K,295 K,345 K,and 395 K.We have found that with the temperature increase,the maximum bond length(L max) of RDX N N trigger bond increases,and the interaction energy(E N-N) between two N atoms of the N-N trigger bond and the cohesive energy density(CED) decrease.These phenomena agree with the experimental fact that the PBX becomes more sensitive as the temperature increases.Therefore,we propose to use the maximum bond length L max of the trigger bond of the easily decomposed and exploded component and the interaction energy E N-N of the two relevant atoms as theoretical criteria to judge or predict the relative degree of heat and impact sensitivity for the energetic composites such as PBXs and solid propellants.     

RELATIONSHIPS BETWEEN X—H STRETCHING OVERTONE FREQUENCIES AND BOND LENGTHS/BOND ANGLES

A relationship between the X-H (X = N, O, C, and so on) equilibrium bond length in a Morse oscillator and the X-H stretching overtone frequency shifts is obtained theoretically. We use the equation to discuss the empirical linear relationships that have been proposed for heterocyclics, alkanes and fluorinated benzenes. On the other hand, a unified relationship between the X-H bond angles and the experimental quantities (ω(?) and the coupling strength λ) is also presented for XH2, XH, and XH4 molecules or molecular fragments. Calculations of X-H bond angles for a number of molecules show that the results from our equations are in excellent agreement with the experimental values. Also we can extract the information of relative magnitude of bond coupling force field.     

SYNTHESIS, CHARACTERIZATION OF IMIDAZOLATE-BRIDGED HETEROPOLY- NUCLEAR COMPLEXES AND THEIR ACTIVITY IN DISMUTATION OF O~(-)_(2)

A series of imidazolate-bridged heteropolynuclear complexes containing Cu2+ or Zn2+ were synthesized and characterized by reflectance spectroscopy, NMR and X-ray diffraction analysis. The bonding nature and the stability of imidazolate bridges in the complexes were studied by ESR spectroscopy, and the catalytic activity of the complexes in dismutation of O-2 was determined by NBT method. Results obtained indicate that the central Cu with N4 and N2O2 square planar or N4O square pyramidal coordination in which there is a weak bond H2O or ClO-4 on axial position, has a comparatively higher activity, but that with N5 square pyramidal having a strong bond axial ligand has almost no activity. Thus the results imply a possible formation of Cu-O-2 intermediate adduct in the catalytic process by Cu, Zn-SOD.     

Structure of cis-[Pt(NH_3)(2-picoline)]~(2+) and DNA adduct and its bonding characteristics

Several methods including molecular mechanics, molecular dynamics, ONIOM that combines quantum chemistry with molecular mechanics and standard quantum chemistry are used to study the configuration and electron structures of an adduct of the DMA segment d(ATACATG*G*TACATA)-d(TATGTACCATGTAT) with cis-[Pt(NH3)(2-Picoline)]2+. The investigation shows that the configuration optimized by ONIOM is similar to that determined by NMR. Strong chemical bonds between Pt of the complex and two N7s of neighboring guanines in the DNA duplex and hydrogen bond between the NH3of the complex and O6 of a nearby guanine have a large impact on the configuration of the adduct. Chemical bonds, the aforementioned hydrogen bond, and the interaction between a methyl of the complex and a methyl of the base in close proximity are critical for the complex to specifically recognize DNA.     

Insights into the Decomposition and Polymerization of Cage Energetic Crystal Octanitrocubane under High Pressure: a DFT-D Study

Dispersion-corrected DFT calculations have been performed to study the crystal structure, initial decomposition mechanism, polymerization mechanism, electronic structure and absorption properties of the most powerful CHNO energetic compound octanitrocubane(ONC) under hydrostatic pressure of 0~200 GPa. Our results show that the lattice parameter c is sensitive to van der Waals interactions and the structure is the stiffest in the a direction. ONC decomposes by the breaking of N–O bond in nitro group at 195 GPa. At 200 GPa, ONC decomposes through the breaking of N–O bond in nitro group and polymerizes through the forming of a new C–O covalent bond between a carbon atom in the cage skeleton and one oxygen atom of nitro group in another ONC molecule. ONC becomes more and more sensitive under high pressure and transforms into a metallic phase in the pressure range of 175~200 GPa. ONC also has higher optical activity and wider absorption range under compression.     

Ni-mediated C–N activation of amides and derived catalytic transformations

Amide, as a ubiquitous functional group, is essential in various aspects of chemistry and biology. Although the history of studying amide is rich and fruitful, the synthetic application of amide is very limited due to the inertness of amide C–N bond. Recently,significant advances have been achieved towards the nickel-mediated C–N activation of amides. This approach allows a facile generation of acyl-nickel intermediates, and a number of unique transformations have been designed and realized based on the amide C–N bond activation. Focused on the catalytic transformation, this review summarizes and categorizes the recent advances on the synthetic applications of Ni-mediated C–N bond activation of amides.     

Formation Mechanism and Stability Study of Konjac Glucomannan Helical Structure

The formation mechanism and stability of konjac glucomannan (KGM) helical structure were investigated by molecular dynamic simulation and experimental method. The results indicate that the molecular conformation of KGM is a non-typical helical structure. In detail, helical structure of KGM is mainly sustained by acetyl group, whose size and stability are affected by the molecular polymerization degree of KGM. In vacuum among the non-bonding interactions, electrostatic force is the greatest factor affecting its helical structure, but in water solution, hydrogen bond affects the helical arrangement greatly. To our interest, temperature exhibits a reversible destroying effect to some extent; the helical structure will disappear completely and present a ruleless clew-like arrangement till 341 K. This work suggests that the method of combining molecular dynamic simulation and experiment tools can be effective in the study of KGM helical structure.     

Co-existing Intermolecular Halogen Bonding and Hydrogen Bonding in the Compound Trans-5,10-bis（1-bromodifluoro-acetyl-1-ethoxycarbonyl-methylidene）thianthrene

Trans-5,10-bis(1-bromodifluoroacetyl-l-ethoxycarbonyl-methylidene)thianthrene (1b) was prepared from the reaction of BrCF2COC(N2)CO2Et with thianthrene. X-ray single crystal diffraction analysis showed that the intermolecular halogen bonding and hydrogen bonding coexisted in this compound. The bromine atom acted as an electron acceptor in the halogen bond and an electron donor in the hydrogen bond. It is the first example that the bromine atom acted as such a dual role in the hydrogen and halogen bond.     

Adsorption of HCN on Ni/Pt(111) Bimetallic Surfaces Investigated with Density Functional Theory Method

We applied periodic density-functional theory to investigate the adsorption of HCN on x Ni@Pt(111) bimetallic surfaces(x = 1~4). The results have been compared with those obtained on pure Ni(111) and Pt(111) surfaces. For all bimetallic surfaces,HCN is preferentially tilted with the CN bond parallel to the surface,and adsorption energies increase with an increasing number of layer Ni atoms on the surface. The adsorption energies of HCN on all bimetallic surfaces are larger than that on the Pt(111) surface,whereas the adsorption energies of HCN on 3Ni@Pt(111) and 4Ni@Pt(111) are larger than that on the Ni(111) surface,indicating that the introduction of Ni to the Pt catalyst could increase the activity of bimetallic catalyst in the hydrogenation reaction for nitriles. Larger adsorption energy of HCN leads to a longer C–N bond length and a smaller CN vibrational frequency. The analysis of Bader charge and vibrational frequencies showed obvious weakening of the adsorbed C–N bond and an indication of sp2 hybridization of both carbon and nitrogen atoms.     

Pentazole anion cyclo-N_5~-: a rising star in nitrogen chemistry and energetic materials

正Nitrogen and carbon are common elements in nature. The development of nitrogen chemistry, however, is obviously left behind by that of carbon chemistry. Prof. Christ [1] attributed this reason to the stability difference between nitrogen-nitrogen bond and carbon-carbon bond. From N≡N,N=N to N–N bond, the bond dissociation energy increases significantly and the stability becomes worse rapidly, which is in contrast with stable C–C bond and unstable C≡C bond in carbon chemistry. Thus, lots of stable "polycarbon"     

SYNTHESIS OF 1-SUBSTITUTED-BENZENESULFO-NYL-3-(4'''',6''''-DISUBSTITUTED-s-TRIAZIN.2''''-YL)-1,3,2-DIAZAPHOSPHOLIDINES AND THEIR STRUCTURAL ANALYSES

The reaction of N-substituted-benzenesulfonyl-N'-(4,6-disubstituted-s-triazin-2-yl)-ethylene diamines with tri(dialkylamine)phosphine has been studied and the cyclic products have been obtained by one-pot reaction. It was found from ~1H NMR and ~(13)C NMR analyses that the single bond rotation between the nitrogen atom 3 in the 1,3,2-diazaphospholidines and the carbon atom 2 in the s-triazine was hindered at room temperature, which was further confirmed by X-ray diffraction results. The most stable conformation of the phosphorous heterocycle is envelope-form. From the bond length and the hybridization of the nitrogen atom outside the ring, it was deduced that there was a certain contribution of dπ-pπ bond between P and N atoms.     

Halogen Bonding： An AIM Analysis of the Weak Interactions

A series of complexes formed between halogen-containing molecules and ammonia have been investigated by means of the atoms in molecules （AIM） approach to gain a deeper insight into halogen bonding. The existence of the halogen bond critical points （XBCP） and the values of the electron density （Pb） and Laplacian of electron density （V2pb） at the XBCP reveal the closed-shell interactions in these complexes. Integrated atomic properties such as charge, energy, polarization moment, volume of the halogen bond donor atoms, and the corresponding changes （△） upon complexation have been calculated. The present calculations have demonstrated that the halogen bond represents different AIM properties as compared to the well-documented hydrogen bond. Both the electron density and the Laplacian of electron density at the XBCP have been shown to correlate well with the interaction energy, which indicates that the topological parameters at the XBCP can be treated as a good measure of the halogen bond strength In addition, an excellent linear relationship between the interatomic distance d（X…N） and the logarithm of Pb has been established.     

Conversion of N2O to N2 on MgO （001） Surface with Vacancy： A DFT Study

The adsorption and decomposition of NzO at regular and defect sites of MgO (001) surface have been studied using cluster models embedded in a large array of point charges (PCs) by DFT/B3LYP method. The results indicate that the MgO (001)surface with oxygen vacancies exhibits high catalytic reactivity toward N2O adsorptive-decomposition. It is different from the regular MgO surface or the surface with magnesium vacancies.Much elongation of O—N bond of N2O after adsorption at oxy-gen vacancy site with O end down shows that O—N bond has been broken with concurrent production of N2, leaving a regu-lar site instead of the original oxygen vacancy site (F center ).The MgO (001) surface with magnesium vacancies hardly ex-hibits catalytic reactivity. It can be concluded that N2O dissoci-ation likely occurs at oxygen vacancy sites of MgO (001) sur-face, which is consistent with the generally accepted viewpoint in the experiments. The potential energy surface (PES) reflects that the dissociation process of N2O does not virtually need to surmount a given energy barrier.     

~(11)B NMR AND QUANTUM CHEMISTRY STUDIES OF BORATRANE COMPOUNDS

The ~(11)B chemical shifts of boratrane compounds were measured. Quantum chemistry calculations of some of the compounds wore carried out using the CNDO/2 method. Both experimental and calculated results confirmed the existence of the B←N dative bonds in boratranes. The exporimental results indicate that the ~(11)B chemical shifts are influenced by the size of the rings, the characteristics of substituents and steric factors. The calculated results indicate that the strengths of the B←N dative bonds in boratranes are about half of that of a purely "covalent two-electron bond" and therefore, the dative bonds B←N are stronger than those of Si←N in the analogous silatranes, and the boratranes are more stable than silatranes.     

EFFECTS OF SECONDARY STRUCTURE ON ELASTIC BEHAVIOR OF PROTEIN-LIKE CHAINS

The elastic behavior of protein-like chains was investigated by using the Pruned-Enriched-Rosenbluth Method (PERM).Three typical protein-like chains such as all-α,all-β,and α+β(α/β) proteins were studied in our modified orientation dependent monomer-monomer interaction (ODI) model.We calculated the ratio of /N and shape factor <δ*> of protein-like chains in the process of elongation.In the meantime,we discussed the average energy per bond ＜U＞/N,average contact energy per bond ＜U＞c/N,average helical energy per bond ＜U＞h/N and average sheet energy per bond ＜U＞b/N.Three maps of contact formation,α-helix formation,β-sheet formation were depicted.All the results educe a view that the helix structure is the most stable structure,while the other two structures are easy to be destroyed.Besides,the average Helmholtz free energy per bond ＜A＞/Nis was presented.The force f obtained from the free energy was also discussed.It was shown that the chain extended itself spontaneously first.The force was studied in the process of elongation.Lastly,the energy contribution to elastic force fu was calculated too.It was noted that fu for all-β chains increased first,and then decreased with x0 increasing.     

Theoretical Studies of Water’s and Methanol’s Effects on Alcoholysis of N-Benzyl-3-oxo-β-sultam

The mechanisms about the water’s and methanol’s effects on the alcoholysis of N-benzyl-3-oxo-β-sultam together with their differences have been studied by using density func- tional theory at the B3LYP/6-31G* level. The results, in comparison with a previous study on the relative reaction without the assistance of water and methanol, show that the added water or methanol can remarkably reduce the energy barrier of alcoholysis reaction of N-benzyl-3-oxo- β-sultam and the most favorite pathway is the breaking of C–N bond instead of S–N. It is also found that the reaction energy barrier of methanol-assisted alcoholysis is a little higher than that of the water-assisted one.     

Theoretical Study of NO Dimer Adsorption and Dissociation on the CuCr2O4 （100） Surface

Theoretical simulation of the adsorption and dissociation of two NO molecules at the Cu^2＋, Cr^3＋ and bridge Cr^3＋ sites （b-Cr^3＋） on the normal spinel CuCr2O4 （100） surface has been carried out by density functional theory calculations. The results show that the formed N-down and O-down NO dimers are negatively charged. The formation of stable O-down dimers on the surface leads to the great elongation of N-O bond, which contributes to the NO reduction. The transition-state calculations indicate that the decomposition of O-down NO dimer at the b-Cr^3＋ site is most favorable and N2O is the major reduction product.     

Theoretical Study on Dihydrogen Bonds of NH3BH3 with Several Small Molecules

The dihydrogen bonds B-H...H-X （X= the complexes of NH3BH3 with HF, HCl, F, Cl, Br, C, O, N） in the dimer （NH3BH3）2 and HBr, H2CO, H20, and CH3OH were theoretically studied. The results show that formation of the dihydrogen bond leads to elongation and stretch frequency red shift of the BH and XH bonds, except that in the H2CO system, the CH bond blue shifts. For （NH3BH3）2 and the complexes of the halogenides, red shifts of the XH bonds are caused by the intermolecular hyperconjugation σ（BH）→σ^＊ （XH）. For the system of H2CO, a blue shift of the CH bond is caused by a decrease of the intramolecular hyperconjugation n（O→σ^＊ （CH）. In the other two systems, the red shift of OH bond is a secondary effect of the stronger traditional red-shifted H-bonds N-H... O. In all these systems, red shifts of the BH bonds are caused by two factors： negative repolarization and negative rehybridization of the BH bond, and decrease of occupancy on σ（BH） caused by the intermolecular hyperconjugation σ（BH）→σ^＊ （XH）.     

Theoretical Studies on the Hydrogen Bonds of Different Position Action Mechanisms of Thymine with Uracil

In this research, the hydrogen bonds Y···H-X(X = C, N; Y = N, O) of thymine and uracil have been theoretically studied. The results show that hydrogen bond leads to bond length elongation and stretches the frequency red-shift of N-H···Y. Meanwhile, the C-H···O bonds shorten and stretch the frequency blue-shift. They all belong to traditional hydrogen bonds. The intermolecular charge transfer caused by the intermolecular hyperconjugation ρ*(N–H) →n(Y) and intramolecular charge redistribution by intramolecular hyperconjugation ρ(C-H)→ρ*(C-N) play important roles in the formation of hydrogen bonds. According to the judgment standards proposed by Bader and Popelier, these hydrogen bonds have typical electron density topological properties. Electrostatic surface potential(ESP) is a useful physicochemical property of a molecule that provides insights into inter- and intramolecular associations, as well as the prediction of likely sites of electrophilic and nucleophilic metabolic attack.