Theoretical Investigations on the Structure,Density, Thermodynamic and Performance Properties of Bis（2,2-dinitropropyl） formal

Density functional method was used to investigate the IR spectrum,heat of formation and thermal stability of a new energetic material bis(2,2-dinitropropyl) formal(BDNPF).The detonation velocity and pressure were evaluated by using the Kamlet-Jacobs equations based on the theoretical density and heat of formation.The bond dissociation energies for the weakest bonds were analyzed to investigate the thermal stability of the title compound.The results show that the C(1)-N(1) bond is predicted to be the trigger bond during pyrolysis.The crystal structure obtained by molecular mechanics belongs to the P21 space group,with the lattice parameters to be Z = 2,a = 11.5254,b = 6.2168,c = 9.5000  and ρ= 1.66 g/cm3.     

Computational studies on the crystal structure, thermodynamic properties, detonation performance, and pyrolysis mechanism of 2,4,6,8-tetranitro-1,3,5,7-tetraazacubane as a novel high energy density material

Studies have suggested that octanitrocubane (ONC) is one of the most powerful non-nuclear high energy density material (HEDM) currently known. 2,4,6,8-Tetranitro-1,3,5,7-tetraazacubane (TNTAC) studied in this work may also be a novel HEDM due to its high nitrogen content and crystal density. Density functional theory and molecular mechanics methods have been employed to study the crystal structure, IR spectrum, electronic structure, thermodynamic properties, gas-phase and condensed-phase heat of formation, detonation performance, and pyrolysis mechanism of TNTAC. The TNTAC has a predicted density of about 2.12 g/cm(3), and its detonation velocity (10.42 km/s) and detonation pressure (52.82 GPa) are higher than that of ONC. The crystalline packing is P2(1)2(1)2(1), and the corresponding cell parameters are Z = 4, a = 8.87 ?, b = 8.87 ?, and c = 11.47 ?. Both the density of states of the predicted crystal and the bond dissociation energy of the molecule in gas phase show that the cage C-N bond is the trigger bond during thermolysis. The activation energy of the pyrolysis initiation reaction obtained from the B3LYP/6-311++G(2df,2p) level is 125.98 kJ/mol, which indicates that TNTAC meets the thermal stability request as an exploitable HEDM.     

Theoretical Insights into Intermolecular Hydrogen-Bonding Strengthening in Fluorenone-Methanol Complexes Induced by Electronic Excitation and Bulk Solvent Effect

This work presents a theoretical insight into the variation of the site-specific intermolecular hydrogen-bonding （HB）, formed between C=O group of fluorenone （FN） and O-H groups of methanol （MeOL） molecules, induced by both the electronic excitation and the bulk solvent effect. Through the calculation of molecular ground- and excited-state properties, we not only demonstrate the characters of HB strengthening induced by electronic excitation and the bulk solvent effect but also reveal the underlying physical mechanism which leads to the HB variation. The strengthening of the intermolecular HB in electronically excited states and in liquid solution is characterized by the reduced HB bond-lengths and the red-shift IR spectra accompanied by the increasing intensities of IR absorption corresponding to the characteristic vibrational modes of the O-H and C--O stretching. The HB strengthening in the excited electronic states and in solution mainly arises from the charge redistribution of the FN molecule induced by the electronic excitation and bulk solvent instead of the intermolecular charge transfer. The charge redistribution of the solute molecule increases the partial dipole moment of FN molecule and the FN-MeOL intermolecular interaction, which subsequently leads to the HB strengthening. With the bulk solvent effect getting involved, the theoretical IR spectra of HBed FN-MeOL complexes agree much better with the experiments than those of gas-phase FN-MeOL dimer. All the calculations are carried out based on our developed analytical approaches for the first and second energy derivatives of excited electronic state within the time-dependent density functional theory.     

Stable Structure and Characteristic Spectroscopy of （CH3COCH3）n （n = 1-7） Clusters

The stable acetone molecule clusters have been studied by using Becke＇s threeparameter （B3LYP） density functional theory （DFT） with standard 6-31G（d） basis set. The calculated results show that the optimal structures of acetone clusters are cyclic and the cycles become larger and larger with the increase of cluster size. The strongest vibration peaks for neutral clusters are C=O stretching vibration. The C=O stretching peaks of cyclic acetone clusters split into double ones when n≥3, the frequencies are red-shifted and corresponding intensifies increase with the increase of cluster Size.     

Theory Study of the Adsorption of Hydrocyanic Acid onto Small Silver Clusters

Density functional theory calculations have been performed to study the interaction of small silver clusters, Ag2 ~Ag9, with HCN. The adsorption of HCN on-top site of the silver cluster, among various possible sites, is energetically preferred. The adsorption energies of HCN on the silver clusters reach a local maximum at n = 4, which is only about 0.450 eV, indicating that the adsorbed HCN molecule is weakly perturbed. The adsorbed C–N and C–H stretching frequencies are blue- and red-shifted compared with the values of free HCN, respectively.     

UV-Vis Spectrum and the Third-order Nonlinear Optical Properties of the Chiral Camphorderived β-diketonate Platinum Complexes

UV-Vis spectrum and the third-order nonlinear optical properties of the chiral camphor-derived β-diketonate have been studied at the B3LYP/6-31G* level. The results showed that the introduction of electron-drawing group-CF3 and-C3F7 on β-diketonate made the strongest absorption peak red-shift and the lowest energy absorption blue-shied. Introduction of-OC2H5 on the benzene or pyridine ring made the lowest energy absorption blue-shift. When the-C2H3 was introduced on the benzene or pyridine ring, the lowest energy absorption was red-shifted. Introduction of electron-donating group on β-diketonate can enlarge their nonlinear optical properties. On the contrary, the introduction of electron-drawing group dropped it down.     

Bonding of Hydroxyl and Epoxy Groups on Graphene： Insights from Density Functional Calculations

Density functional theory and GGA-PW91 exchange correlation function were performed to simulate the bonding behavior of hydroxyl and epoxy groups on the graphene surface. We compared the different binding energies for two epoxy groups, as well as one hydroxyl group and one epoxy group on all possible positions within a 6-fold ring, respectively. The calculated results suggest that two oxygen-containing groups always tend to bind with the neighboring carbon atoms at the opposite sides. Moreover, two hydroxyl groups on the meta position are unstable, and one of the hydroxyl groups easily migrates to the para position. In contrast to the disperse arrangement, the aggregation of multiply hydroxyl groups largely enhances the binding energy of every hydroxyl group. It is worth noting that the binding sites and hydrogen bonds play an important role in stability. Our work further points out the number of oxygen-containing groups and the location of oxide region largely influence the electronic properties of graphene oxide.     

Experimental and Theoretical Study of Reactions between Manganese Oxide Cluster Cations and Hydrogen Sulfide

Manganese oxide cluster cations Mnm180n＋ were prepared by laser ablation and reacted with hydrogen sulfide （H2S） in a fast flow reactor under thermal collision conditions. A time-of-flight mass spectrometer was used to detect the cluster distributions before and after the interactions with H2S. The experiments suggest that oxygen-for-sulfur （O/S） ex- change reaction to release water took place in the reactor for most of the manganese oxide cluster cations： MnmlSOn＋＋H2S→Mnm18On-1S＋＋H218O. Density functional theory cal- culations were performed for reaction mechanisms of Mn202＋＋H2S, Mn203＋＋H2S, and Mn204＋＋H2S. The computational results indicate these O/S exchange reactions are both thermodynamically and kinetically favorable, thus in good agreement with the experimental observations. The O/S exchange reactions identified in this gas-phase cluster study parallel similar behavior of related condensed phase reaction systems.     

The First-principlies Calculations of the Electronic Structures and Optical properties of Ⅱ-Ⅲ2-Ⅵ4 （Ⅱ=Zn,Cd; Ⅲ= In; Ⅵ=Se,Te）

The electronic structures and optical properties of II-III2-VI4 （II = Zn, Cd; III = In; VI = Se, Te） compounds are studied by the density functional theory （DFT） using the Vienna ab initio simulation package （VASP）. Geometrical optimization of the unit cell is in good agreement with the experimental data. Our calculations show that the valence band maximum （VBM） and conduction band minimum （CBM） are located at G resulting in a direct energy gap. The optical properties are analyzed, and the independent second harmonic generation （SHG） coefficients are determined. By an analysis of the band structure, we can get that SHG response of the system can be attributed to the transitions from the bands near the top of valence band that are derived from the Se/Te p states to the unoccupied bands contributed by the p states of In atoms.     

Theoretical Study on Nitroso-Substituted Derivatives of Azetidine as Potential High Energy Density Compounds

At the B3 PW91/6-311+G(d,p)//MP2/6-311+G(d,p) level, molecular densities, detonation velocities, and detonation pressures of nitroso substituted derivatives of azetidine with their thermal stabilities were investigated to look for high energy density compounds(HEDCs). It was found that the azetidine derivatives had high heat of formation(HOF) and large bond dissociation energy(BDE). Intramolecular hydrogen bonds were located in three molecules(1, 4, and 5), and the molecular stability were improved markedly as well. For 5 and 6, the detonation performances(D= 9.36 km/s and 10.80 km/s, P= 44.42 GPa and 60.70 GPa, respectively) meet requirements as high energy density compounds. This work may provide basic information for further study of title compounds.     

3,6-二肼基-1,2,4,5-四嗪的晶体结构及理论研究

通过缓慢蒸发溶剂法培养得到3, 6-二肼基-1, 2, 4, 5-四嗪(DHT)的单晶, 用X-射线单晶衍射仪进行了结构测定. 该晶体属于单斜晶系, P2(1)/c空间群, a=4.032 (4) ?, b=5.649 (6) ?, c=12.074 (14) ?. β=99.32°, Z=2, V=271.4(5) ?3. DHT分子中肼基N原子与四嗪环基本位于同一平面呈现轴对称结构, 分子中的大量氢键使之形成箭尾形排列的三维网络结构. 通过实验测得DHT的燃烧热为1787kJ?mol-1, 5s爆发点为454 K. 在DFT-B3LYP/6-311G*水平下对DHT的电子结构和自然键轨道进行了分析. 通过原子化能的方法计算得到DHT的标准生成热为1075 kJ?mol-1, 与实验值接近. 爆轰性能计算表明, DHT在密度为1.64 g?cm-3时, 爆速和爆压分别为9.27 km?s-1和36.02 GPa, 高于TNT和HMX.     

DFT studies on a high energy density cage compound 4-trinitroethyl-2,6,8,10,12-pentanitrohezaazaisowurtzitane

Polynitro cage compound 4-trinitroethyl-2,6,8,10,12-pentanitrohexaazaisowurtzitane has the same framework with but higher stability than CL-20 and is a potential new high energy density compound (HEDC). In this paper, the B3LYP/6-31G(d,p) method of density functional theory (DFT) has been used to study its heat of formation, IR spectrum, and thermodynamic properties. The stability of the compound was evaluated by the bond dissociation energies. The calculated results show that the first step of pyrolysis is the rupture of the N-NO(2) bond in the side chain and verify the experimental observation that the title compound has better stability than CL-20. The crystal structure obtained by molecular mechanics belongs to the P2(1)2(1)2(1) space group, with lattice parameters a = 12.59 ?, b = 10.52 ?, c = 12.89 ?, Z = 4, and ρ = 2.165 g·cm(-3). Both the detonation velocity of 9.767 km·s(-1) and the detonation pressure of 45.191 GPa estimated using the Kamlet-Jacobs equation are better than those of CL-20. Considering that this cage compound has a better detonation performance and stability than CL-20, it may be a superior HEDC.     

Quantum Chemical Studies on Structure and Detonation Performance of Bis(2,2-dinitropropyl ethylene)formal

Based on the full optimized molecular geometric structures via B3LYP/6-311+G(2d,p) method, a new gem-dinitro energetic plasticizer, bis(2,2-dinitropropyl ethylene)formal was investigated in order to search for high-performance energetic material. IR spectrum, heat of formation, and detonation performances were predicted. The bond dissociation energies and bond orders for the weakest bonds were analyzed to investigate the thermal stability of the title compound. The results show that the four N-NO₂ BDEs are nearly equal to the values of 164.38 kJ/mol, which shows that the title compound is a stable compound. The detonation velocity and pressure were evaluated by using Kamlet-Jacobs equations basedon the theoretical density and condensed HOF. The crystal structure obtained by molec-ular mechanics belongs to P2₁ space group, with lattice parameters Z=2, a=13.8017 ?, b=13.4072 ?, c=5.5635 ?.     

Supramolecular Self-assembly and Functionalization of Porphyrin-based Systems

Porphyrins are abundant in nature. They have been frequently employed as building blocks in the construction of nanoarchitectures and functional supramolecular systems. Recently, a series of novel porphyrin molecules including small molecules and polymers have been originally designed and synthesized with the aim of producing nanostructures with controllable-growth and materials with high-performance. Literature coverage is through 2004–2012. This review gives a full summary of related studies in our group.     

Frist-principles Band Structures Calculation of Tin-phthalocyanine

We adopt the density function theory with generalized approximation by the Beeke exchange plus Lee-Yang-Parr correlation functional to calculate the electronic first-principles band structure of tin-phthalocyanine （SnPc）. The intermolecular interaction related to transport behavior was analyzed from the F-point wave function as well as from the bandwidths and band gaps. From the calculated bandwidths of the frontier bands as well as the effective masses of the electron and hole, it can be concluded that the mobility of the electron is about two times larger than that of the hole. Furthermore, when several bands near the Fermi surface are taken into account, we find that the interband gaps within the unoccupied bands are generally smaller than those of the occupied bands, indicating that the electron can hop from one band to another which is much easier than the hole. This may happen through electron-phonon coupling for instance, thus effectively yielding an even larger mobility for the electron than for the hole. These facts indicate that in SnPc the electrons are the dominant carriers in transport, in contrast to most organic materials.     

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.     

Theoretical Investigation on Interaction between Guanine and Luteolin

The interacting patterns of the luteolin and guanine have been investigated by using the density functional theory B3LYP method with 6-31＋G＊ basis set. Eighteen stable structures for the luteolin-guanine complexes have been found respectively. The results indicate that the complexes are mainly stabilized by the hydrogen bonding interactions. Meanwhile, both the number and strength of hydrogen bond play important roles in determining the stability of the complexes which can form two or more hydrogen bonds. Theories of atoms in molecules and natural bond orbital have also been utilized to investigate the hydrogen bonds involved in all the systems. The interaction energies of all the complexes which were corrected by basis set superposition error are 6.04-56.94 kJ/mol. The calculation results indicate that there are strong hydrogen bonding interactions in the luteolin-guanine complexes. We compared the interaction between luteolin and four bases of DNA, and found luteolin-thymine was the strongest and luteolin-adenine was the weakest. The interaction between luteolin and DNA bases are all stronger than luteolin-water.     

Selective Sensing Characteristics of Ca Doped BeO Nano-sized Tube toward H20 and NH3

By means of density functional calculations, the structural and electronic properties of chemical modification of pristine and Ca-doped BeO nanotubes were investigated with NH3 and H20 molecules. It was found that the NH3 and H20 molecules can be adsorbed on the Be atom of the tube sidewall with the adsorption energies of about 36.1 and 39.0 kcal/mol, respectively. Density of states analysis shows that the electronic properties of the BeONT are slightly changed after the adsorption processes. Substitution of a Be atom in the tube surface with a Ca atom increases the adsorption energies by about 7.4 and 14.7 kcal/mol for NH3 and H20, respectively. Unlike the pristine tube, the electronic properties of Ca-doped BeONT are sensitive to NH3 and H20 molecules. Also, the Ca-doped tube is much more sensitive to H20 molecule than NH3 one.     

Theoretical Investigation on the Absence of Spore Photoproduct Analogue at Cytosine-Thymine Site

As a well known DNA photolesion product, the special UV induced pyrimidine dimmer called spore photoproduct （SP）, has aroused strong research interests. The SP formation was reported solely between two adjacent thymidine residues. It remains unclear in pervious experimental observations why there is an absence of the cytosine-derived SP-like photoprod- uct formation at the cytosine containing DNA strand, although the cytosine residue holds great similarity to thymine in terms of molecular structure. From a theoretical perspec- tive, we have explored this issue in this work by means of density functional theory at the B3LYP/6-311＋＋G（d,p） //B3LYP/6-31G（d,p） level for the DNA dinucleotide fragment, cy- tosine phosphate thymine （CpT）. Key factors blocking the formation of the SP-like product between two adjacent cytosine and thymidine residues are revealed. Instead of undergoing photochemical SP reaction, a photophysical deactivation pathway back to the ground state turns out to be favorable for the CpT dinucleotide fragment.