Computational method on highly efficient D-π-A-π-D-based different molecular acceptors for organic solar cells applications and non-linear optical behaviour.

Eight molecular structures (BT-A1 to BT-A8) with high-performance non-fullerene acceptor (NFA) were selected for organic solar cells (OSCs) and non-linear optical (NLO) applications. Their electronic, photovoltaic (PV) and optoelectronic properties were tuned by adding powerful electron-withdrawing groups to the acceptor (A) of the D-π-A-π-D structure. Using time-dependent density functional theory (TD-DFT) techniques, based on the laws of quantum chemical calculations, the absorption spectra, stability of the highest and lowest-energy molecular orbitals (HOMO/LUMOs), electron density, intramolecular charge transfer (ICT), transition density matrix (TDM), were examined. The binding energy (Eb) and density of states (DOS) were probed to realize the optoelectronic analysis of the structures BT-A1 to BT-A8. Noncovalent interactions (NCIs) based on a reduced density gradient (RDG) were used to describe the nature and strength of D-A interactions in the molecules BT-A1 to BT-A8. The new refined molecules BT-A1 to BT-A8 exhibited strong absorbance bands between 408-721 nm and high electron transfer contribution (ETC) ranges between 87-96 %, along with the smallest excitation energies (Ex) between 1.71-3.55 eV in the solvent dichloromethane. Dipolar moment strengths ranging from 0.38 to 4.72 Debye in both the excited and ground states have determined with good solubility properties of BT-A1 to BT-A8 in polar solvent. Highly effective charge mobilities and prevention of charge recombination have been demonstrated by the electron (0.18-0.41 eV) and hole RE values (0.13-0.89 eV) for the new compounds. Power conversion efficiencies (PCE) of BT-A1 to BT-A8 were nearly the same because of better outcomes compared to the molecules in the BT. Compared to poly[4.8-bis(5-(2-ethylhexyl)thiophen-2-yl)benzo[1,2-b: 4,5-b']dithiophene-2,6- diyl-alt-(4-2-ethylhexyl)-3-fluorothieno[3,4-b]thiophene-)-2-carboxylate-2-6-diyl)] (PTB7-Th), the open circuit voltages (Voc) of compounds BT-A1 to BT-A8 were ranged from 1.52 to 2.13 eV. The polarizability (α) and hyperpolarizability (ß) of the molecules BT-A1 to BT-A8 were used to determine the non-linear optical (NLO) properties. The results showed that BT-A2, BT-A6 and BT-A7 have good NLO activity. This computational analysis demonstrates the superiority of the molecules with NFA. Hence the compounds are advised for the use in production of high-performance OSCs and NLO activity.

Theoretical investigation on regulating photophysical properties and proton transfer behavior by electronegativity for near-infrared emitting styryl dyes.

Our main focus is to explore the atomic electronegativity-dependent photoinduced behavior of styryl derivatives (HBO, HBS, and HBSe). The results of structural parameter calculation by the DFT method show that the intramolecular hydrogen bonds of normal and tautomer form are strengthened and weakened, respectively, in an excited state (S1), which is conducive to the excited intramolecular proton transfer (ESIPT) process. The enhancement of excited hydrogen bond is beneficial to the ESIPT process from the aspects of infrared vibration frequency (IR), Mulliken's charge analysis, and density gradient reduction (RDG). Additionally, by determining the bond energy with the band critical point (BCP) parameter, we found that the lower the electronegativity of the atom, the larger the hydrogen bond strength at the excited state and the more likely ESIPT reaction occurs. Meanwhile, the intramolecular H-bonds O-H…N in HBO, HBS, and HBSe are enhanced with the weakened electron-withdrawing capacity of the atom (from O to S and Se). Subsequently, frontier molecular orbital (FMOs) and charge density difference (CDD) analyses essentially revealed that electron redistribution induces the ESIPT process. Low atomic electronegativity exhibits the high chemical activity of the excited state. Furthermore, to demonstrate the electronegativity-dependent ESIPT behavior of the system, we built potential energy curves (PECs) and located the transition states (TS) of proton transfer processes.

Theoretical design and evaluation of efficient small donor molecules for organic solar cells.

CONTEXT: The development of high-efficiency photovoltaic devices is the need of time with increasing demand for energy. Herein, we designed seven small molecule donors (SMDs) with A-π-D-π-A backbones containing various acceptor groups for high-efficiency organic solar cells (OSCs). Molecular engineering was performed by substituting the acceptor group in the synthesized compound (BPR) with another highly efficient acceptor group to improve the photoelectric performance of the molecule. METHOD: The photovoltaic, optoelectronic, and photophysical properties of the proposed compounds (BP1-BP7) were investigated in comparison to BPR using DFT and TD-DFT at MPW1PW91/6-311G(d,p) level of theory. All molecules we designed have red-shifted absorption spectra. The modification of the acceptor fragment of the BPR resulted in a reduced HOMO-LUMO energy gap; thus, the designed compounds (BP1-BP7) had improved optoelectronic responses as compared with the BPR molecule. Various key factors that are crucial for efficient SMDs such as exciton binding energy, frontier molecular orbitals (FMOs), absorption maximum (λmax), open circuit voltage (VOC), dipole moment (µ), excitation charge mobilities, and the transition density matrix of (BPR, BP1-BP7) have also been studied. Low reorganizational energy (holes and electrons) values provide high charge mobility, and all the designed compounds are efficient in this regard. Here, BP6 exhibits low excitation energy (1.66 eV), highest open circuit voltage (2.00 V), normalized VOC (77.23), and fill factor (0.931). Consequently, the superiority of the designed molecules advises experimenters to envision future developments in extremely effective OSC devices.

Crystal morphology prediction of CL-20 and 1,4-DNI co-crystal at different temperatures.

CONTEXT: The morphologies of hexanitrohexaazaisowurtzitane (CL-20) and 1,4-dinitroimidazole (1,4-DNI) co-crystal under vacuum or solvent at different temperatures were predicted. The CL-20/1,4-DNI co-crystal has six important growth crystal planes: (002), (011), (101), (11‒1), (110), (111). The areas of (002), (101), and (011) planes account for a relatively large proportion, which are important crystal planes that affect the crystal morphology. The crystal habits at different temperatures were simulated. The simulation results showed that the crystal plane attachment energy of CL-20 and 1,4-DNI co-crystal increases with the increase of temperature, indicating that the increase of temperature is conducive to the growth of crystal planes. The aspect ratio decreases with the increase of temperature and the morphology of co-crystal becomes more spherical at a higher temperature. The theoretical predictions are in good agreement with the experiment. The simulation results can provide guidance for the crystallization of CL-20/1,4-DNI to obtain a nearly spherical crystal morphology. METHODS: The CL-20/1,4-DNI unit cell structure was geometrically optimized by the COMPASS force field. The AE model was used to predict the morphology of CL-20/1,4-DNI under vacuum, resulting in the most morphologically important growth planes. Ethyl acetate was selected as the solvent. The interaction energy between the solvent and the crystal plane, and the attachment energies in solvent at 298 K, 320 K, 340 K, 360 K, and 380 K were predicted. The NVT ensemble is used in the molecular dynamics calculation process. The simulation step is 1 fs and the total simulation time is 500 ps. The Andersen thermostat is selected as the temperature control method. In the potential energy calculation, the atom-based and Ewald methods were selected to calculate the van der Waals force and the electrostatic interaction force, respectively.

Facile synthesis of triazine-based microporous organic network for high-efficient adsorption of flumequine and nadifloxacin: A comprehensive study on adsorption mechanisms and practical application potentials.

Flumequine (FLU) and nadifloxacin (NAD), as emerging contaminants, have received extensive attention recently. In this study, a triazine-based microporous organic network (TMON) was synthetized and developed as an excellent adsorbent for FLU and NAD. The adsorption behavior and influence factors were investigated in both single and binary systems. Insight into the adsorption mechanisms were conducted through experiments, models, and computational studies, from macro and micro perspectives including functional groups, adsorption sites, adsorption energy and frontier molecular orbital. The results showed that the maximum adsorption capacities of TMON for FLU and NAD are 325.27 and 302.28 mg/g under 30 °C higher than records reported before. TMON exhibits the better adaptability and anti-interference ability for influence factors, leading to the preferable application effect in kinds of real water samples. TMON also shows the application potentials for the adsorption of other quinolone antibiotics and CO2 capture. Hydrogen-bonding interaction played the most critical role compared to π-π stacking effect, π-π electron-donor-acceptor interaction, CH-π interaction, and hydrophobic interaction during the adsorption. TMON could be regarded as a promising environmental adsorbent for its large surface area, stable physical and chemical properties, excellent recyclability, and wide range of applications.

Decomposition mechanism of 1,3,5-trinitro-2,4,6-trinitroaminobenzene under thermal and shock stimuli using ReaxFF molecular dynamics simulations.

To obtain atomic-level insights into the decomposition behavior of 1,3,5-trinitro-2,4,6-trinitroaminobenzene (TNTNB) under different stimulations, this study applied reactive molecular dynamics simulations to illustrate the effects of thermal and shock stimuli on the TNTNB crystal. The results show that the initial decomposition of the TNTNB crystal under both thermal and shock stimuli starts with the breakage of the N-NO2 bond. However, the C6 ring in TNTNB undergoes structural rearrangement to form a C3-C5 bicyclic structure at a constant high temperature. Then, the C3 and C5 rings break in turn. The main final products of TNTNB under shock are N2, CO2, and H2O, while NO,  N2, H2O and CO are formed instead at 1 atm under a constant high temperature. Pressure is the main reason for this difference. High pressure promotes the complete oxidation of the reactants.

Nitrogen-rich polycyclic pentazolate salts as promising energetic materials: theoretical investigating.

Pentazolate (cyclo-N5-) salts are nitrogen-rich compounds with great development potential as energetic materials due to their full nitrogen anion. However, the densities of available N5- salts are generally low, which seriously lowers their performances. It is necessary to screen out cyclo-N5- salts with high density. To this end, eight new non-metallic cyclo-N5- salts based on fused heterocycle were designed. -NH2, -NO2, and -O- groups were introduced into the compounds to adjust and improve the detonation performance and impact sensitivity of cyclo-N5- salts. By theoretical calculations and Hirshfeld surface analyses, the densities, heats of formation, detonation performance, sensitivities, and crystal structures of the title compounds were predicted, and the contribution of hydrogen bond as well as π-π stacking to the stability of cyclo-N5- salt was revealed. The results indicate that the densities of title compounds are higher than 1.85 g cm‒3, and the sensitivities of these compounds are predicted to be lower than that of HMX. The detonation properties of a (D = 9.47 km s-1, P = 41.21 GPa) and d (D = 9.44 km s-1, P = 40.26 GPa) are better than those of HMX. These mean that using fused ring as a cation and introducing proper substituents are an effective method to improve cyclo-N5- salt's density and balance the detonation performance and sensitivity.

Atomic perspective revealing for combustion evolution of nitromethane/nano-aluminum hydride composite.

Adding aluminum hydride (AlH3) into energetic materials (EMs) can improve their combustion and energy performance effectively. However, the potential mechanism of AlH3 on EMs is still unclear. Based on the ReaxFF-lg method, the thermal decomposition of nitromethane/nano-aluminum hydride (NM/nano-AlH3) composites were studied. The addition of AlH3 reduces the energy barrier and increases the energy release during the decomposition of NM, accelerates the decomposition of NM. The main way of AlH3 oxidation involves the capture of O atoms from NM. The results show that AlH3 content and passivated layer affect the oxidation and hydrogen release of AlH3. The explosion of small particle size AlH3 leads to rapid oxidation and hydrogen release. The oxidation of large particle size AlH3 is dominated by the inward and outward diffusion of O and Al atoms. The products of NM/nano-AlH3 composites are H2O, CO2, N2 gases, and Al clusters. This work is expected to guide the application of AlH3 in EMs.

Encapsulating aluminum nanoparticles into carbon nanotubes for combustion: a molecular dynamics study.

Metal nanoparticles are easily deactivated by migration-aggregation in combustion. Encapsulated nanoparticles are one of the tools for coping with the stability challenges of metal nanoparticles. The self-assembly details of aluminum nanoparticles (ANPs) encapsulated into carbon nanotubes (CNTs) were demonstrated by molecular dynamics simulations. The simulation results show that ANPs can completely self-roll into CNTs to form a stable core-shell structure by inertial force and van der Waals force. Inside the tubes, ANPs move toward the cap at a velocity of 2.27 Å ps-1. However, it increases to 3.17 Å ps-1 when near the cap of CNTs. The initiation of the ANPs' oxidation and degradation can be effectively checked by coating CNTs. The diffusion of the Al atoms in the encapsulated ANPs occurred earlier than their oxidation in combustion, verified by using ReaxFF molecular dynamics simulations. The morphological evolutions of the nanostructures in the initial combustion of the encapsulated ANPs are predicted. The interplay between the encapsulated ANPs' responses and external stimuli is classified into core-shell separation, shell damage, and core-shell burst, which provides insights into the oxidation mechanism of encapsulated nanoparticles.

Coarse-grained force-field for large scale molecular dynamics simulations of polyacrylamide and polyacrylamide-gels based on quantum mechanics.

We developed a new coarse-grained (CG) molecular dynamics force field for polyacrylamide (PAM) polymer based on fitting to the quantum mechanics (QM) equation of state (EOS). In this method, all nonbond interactions between representative beads are parameterized using a series of QM-EOS, which significantly improves the accuracy in comparison to common CG methods derived from atomistic molecular dynamics. This CG force-field has both higher accuracy and improved computational efficiency with respect to the OPLS atomistic force field. The nonbond components of the EOS were obtained from cold-compression curves on PAM crystals with rigid chains, while the covalent terms that contribute to the EOS were obtained using relaxed chains. For describing PAM gels we developed water-PAM interaction parameters using the same method. We demonstrate that the new CG-PAM force field reproduces the EOS of PAM crystals, isolated PAM chains, and water-PAM systems, while successfully predicting such experimental quantities as density, specific heat capacity, thermal conductivity and melting point.

Theoretical research about nonmetallic energetic salts with pentazolate anion.

Significant progress has been made in the synthesis of nitrogen-rich high-energy salts by pairing pentazolate anion (cyclo-N5-) with different cations since cyclo-N5- was synthesized. It is difficult to screen out cyclo-N5- salts with high energy quickly and effectively in experiment, while theoretical research can realize this goal. Herein, nineteen high-energy salts, which were composed of tetrazole cation and cyclo-N5- anion, were designed. And their properties were studied via density functional theory and volume-based thermodynamic methods. The results indicate that most salts have high densities, low sensitivities, and good detonation properties. In particular, salt 14 (ρCalib = 1.802 g/cm3, ΔHf = 1058.4 kJ/mol, D = 9.38 km/s, P = 39.10 GPa, h50 = 44.92 cm) exhibits excellent detonation performance (approximating that of 2,4,6,8,10,12-hexanitro-2,4,6,8,10,12-hexaazaisowurtzitane (CL-20)) superior to 1,3,5,7-tetranitro-1,3,5,7-tetrazocine (HMX), and lower impact sensitivity than CL-20 or HMX. Hence salt 14 is regarded as promising candidates for high-performance energetic materials.

Reactive molecular dynamics simulation of thermal decomposition for nitromethane/nano-aluminum composites.

The thermal decomposition of pure nitromethane (NM) and NM/nano-aluminum (Al) composites was simulated by reactive molecular dynamics with ReaxFF-lg corrected force field parameters. The initial decomposition pathway of NM molecules in pure NM is C-N bond rupture. However, NM is decomposed early by the initial pathway of N-O bond rupture when it mixes with nano-Al because of the strong attraction of Al to O. The decomposition process of NM/nano-Al can be divided into three stages: adsorption, slow decomposition, and rapid decomposition. The addition of nano-Al particles decreases the energy barrier in decomposition, increases the released energy, and reduces the decomposition temperature of NM. Adding 3% Al to the explosive can make the detonation pressure 3.083% higher than that of pure system. Compared with pure NM, the energy barrier of 16% Al composite is 25.63 kcal/mol lower and the energy released is 22.99 kcal/mol more. There is an optimal amount of Al contents being added to the NM composite by which the largest total numbers of gaseous products (N2, H2O, and CO2) are released. The effect of Al additives on CO2 production is the most obvious. The maximum detonation pressure can be achieved by adding an appropriate amount of nano-Al, which is similar to the experimental results. Graphical abstract.

Thermal Decomposition Mechanism of 1,3,5,7-Tetranitro-1,3,5,7-tetrazocane Accelerated by Nano-Aluminum Hydride (AlH3): ReaxFF-Lg Molecular Dynamics Simulation.

ReaxFF-low-gradient reactive force field with CHONAl parameters is used to simulate thermal decomposition of 1,3,5,7-tetranitro-1,3,5,7-tetrazocane (HMX) and AlH3 composite. Perfect AlH3 and surface-passivated AlH3 particles were constructed to mix with HMX. The simulation results indicate HMX is adsorbed on the surface of particles to form O-Al and N-Al bonds. The decomposition of HMX and AlH3 composite is an exothermic reaction without energy barrier, but the decomposition of pure HMX needs to overcome the energy barrier of 133.57 kcal/mol. Active nano-AlH3 causes HMX to decompose rapidly at low temperature, and the primary decomposition pathway is the rupture of N-O and C-N bonds. Adiabatic simulation shows that the energy release and temperature increase of HMX/AlH3 is much larger than those of the HMX system. Surface-passivated AlH3 particles only affect the initial decomposition rate of HMX. In HMX and AlH3 composites, the strong attraction of Al in AlH3 to O and the activation of the intermediate reaction by H2 cause HMX to decompose rapidly. The final decomposition products of pure HMX are H2O, N2, and CO2, and those of HMX/AlH3 are H2O, N2, and Al-containing clusters dominated by C-Al. The final gas production shows that the specific impulse of HMX/AlH3 is larger than that of HMX.

Theoretical Study of Effects of Anchoring Groups on Photovoltaic Properties of a Triarylamine-Based p-Type Sensitizer.

The effects of anchoring groups on triarylamine-based p-type dyes were studied by substituting the strong electron-withdrawing carboxyl group with the weak electron-withdrawing pyridyl and the electron-rich catechol groups. Judged by the index t, the charge separation would be improved greatly when the carboxyl group of P4 is replaced by the pyridyl or catechol groups. Although carboxyl as an anchoring group lowers the HOMO energy and facilitates the hole injection in comparison with pyridyl and catechol groups, the weak electron-withdrawing pyridyl and the electron-rich catechol groups facilitate the charge separation. E g becomes narrow as the electron-withdrawing abilities of the anchoring groups decrease or as the conjugation extends. Both the extended π-spacers and the substitution of carboxyl with pyridyl and catechol groups promote the redshifts of adsorption wavelengths. The oscillator strengths for all dyes are over 2.00, indicating that all the dyes are able to harvest the sunlight strongly. The ΔG CR values of P4, DF4, and DZ4 are smaller than those of the other dyes. Also, these dyes have larger adsorption over infrared visible light, indicating that these dyes may be good candidates for p-type DSSCs.

Crystal morphology prediction of 2,2',4,4',6,6'-hexanitrostilbene (HNS) by molecular scale simulation.

The spiral growth model was applied to predict the crystal morphology of 2,2',4,4',6,6'-hexanitrostilbene (HNS). We selected solvents of N,N-dimethylformamide (DMF), N-methyl pyrrolidone (NMP), and nitric acid (NA) to control the crystal morphologies of HNS. Molecular dynamic simulations were used to relax the constructed interface model. The relative growth rate of important face is calculated by the spiral growth expression. The predicted crystal shapes are flaky in three solvents. Only (100), (001), and (011) faces are generated in DMF, NMP, and NA. The aspect ratios of the predicted HNS crystal morphologies in DMF, NMP, and NA are 23.00, 15.45, and 4.85, respectively. In addition, we analyzed the properties on each face using periodic bond chain, molecular arrangement, and roughness model. The excellent agreement between the predicted morphologies and the experimental images is clearly evident. These simulation results can provide guidance for the recrystallization of HNS. Graphical abstract.

Theoretical Investigation of Energetic Salts with Pentazolate Anion.

Energetic salts based on pentazolate anion (cyclo-N5-) have attracted much attention due to their high nitrogen contents. However, it is an enormous challenge to efficiently screen out an appropriate cation that can match well with cyclo-N5-. The vertical electron affinity (VEA) of the cations and vertical ionization potential (VIP) of the anions for 135 energetic salts and some cyclo-N5- salts were calculated by the density functional theory (DFT). The magnitudes of VEA and VIP, and their matchability were analyzed. The results based on the calculations at the B3LYP/6-311++G(d,p) and B3LYP/aug-cc-pVTZ levels indicate that there is an excellent compatibility between cyclo-N5- and cation when the difference between the VEA of cation and the VIP of cyclo-N5- anion is -2.8 to -1.0 eV. The densities of the salts were predicted by the DFT method. Relationship between the calculated density and the experimental density was established as ρExpt = 1.111ρcal - 0.06067 with a correlation coefficient of 0.905. This regression equation could be in turn used to calibrate the calculated density of the cyclo-N5- energetic salts accurately. This work provides a favorable way to explore the energetic salts with excellent performance based on cyclo-N5-.

Crystal Morphology Prediction and Anisotropic Evolution of 1,1-Diamino-2,2-dinitroethylene (FOX-7) by Temperature Tuning.

Temperature-induced morphological changes are one of the strategies for designing crystal shapes, but the role of temperature in enhancing or inhibiting crystal growth is not well understood yet. To meet the requirements of high density and low sensitivity, we need to control the crystal morphology of the energetic materials. We studied the crystal morphology of 1,1-diamino-2,2-dinitroethylene (FOX-7) in dimethyl sulfoxide/water mixed solvent by using the modified Hartman-Perdok theorem. Molecular dynamics simulations were used to determine the interaction of FOX-7 and solvents. The results showed that the crystal shape of FOX-7 is hexagonal, the (101) face is the largest exposed face and is adjacent to six crystal faces at 354 K. As the temperature goes down, the area of the (001) face is significantly reduced. The crystal morphology of FOX-7 at 324 K has a smaller aspect ratio of 4.72, and this temperature is suitable for tuning the morphology from slender hexagon into diamond. The prediction results are in remarkable agreement with the experiments. Moreover, we predicted the evolution path of FOX-7 morphology by Gibbs-Curie-Wulff theorem and explained the variation of crystal shape caused by different external conditions in the actual crystallization process.

Design of high performance p-type sensitizers with pyridinium derivatives as the acceptor by theoretical calculations.

Based on triphenylamine as an electron donor and thiophene as a π-linker, Series P and A p-type sensitizers were designed to investigate the effects of the different acceptors on the properties of the sensitizers. The optimized molecular structures, electronic and optical properties were investigated by density functional theory (DFT) and time-dependent DFT (TD-DFT). The results showed that the properties of the dyes can be tuned by the introduction of the different electron-withdrawing groups to the N atom in the pyridinium acceptor. Compared with the synthesized Series P dyes used in p-type sensitizers, the properties of Series A dyes, except for two dyes that cannot be used as p-type sensitizers, are improved by means of modifying pyridinium acceptors. Due to the suitable electron-withdrawing ability of the hexafluorodiacetylamino group in its acceptor, A6 has the narrowest energy gap (1.90 eV), the largest driving force of hole injection (ΔG inj, -0.68 eV), the high light harvesting efficiency (LHE, 0.9984) and the smallest internal reorganization energy (λ int, 5.05 kcal mol-1). Hence, A6 not only enhances electronic excitation, but also improves the reorganization energy. Importantly, A6 shows the largest red shift and the maximum integral values of the adsorption over the visible light, as well as the strongest adsorption energy (-74.80 kcal mol-1) on a NiO surface. Thus, A6 may be a promising sensitizer for the p-type dye-sensitized solar cells (DSSCs), and the acceptor of A6 may provide a new and easily accessible high performance acceptor for p-type sensitizers.

Reactive molecular dynamics simulation of the high-temperature pyrolysis of 2,2',2'',4,4',4'',6,6',6''-nonanitro-1,1':3',1''-terphenyl (NONA).

2,2',2'',4,4',4'',6,6',6''-Nonanitro-1,1':3',1''-terphenyl (NONA) is currently recognized as an excellent heat-resistant explosive. To improve the atomistic understanding of the thermal decomposition paths of NONA, we performed a series of reactive force field (ReaxFF) molecular dynamics simulations under extreme conditions of temperature and pressure. The results show that two distinct initial decomposition mechanisms are the homolytic cleavage of the C-NO2 bond and nitro-nitrite (NO2 → ONO) isomerization followed by NO fission. Bimolecular and fused ring compounds are found in the subsequent decomposition of NONA. The product identification analysis under finite time steps showed that the gaseous products are CO2, N2, and H2O. The amount of CO2 is energetically more favorable for the system at high temperature or low density. The carbon-containing clusters are a favorable growth pathway at low temperatures, and this process was further demonstrated by the analysis of diffusion coefficients. The increase of the crystal density accelerates the decomposition of NONA judged by the analysis of reaction kinetic parameters and activation barriers. In the endothermic and exothermic stages, a 20% increase in NONA density increases the activation energies by 3.24 and 0.48 kcal mol-1, respectively. The values of activation energies (49.34-49.82 kcal mol-1) agree with the experimental data in the initial decomposition stage.

Molecular Modeling on Morphology of 3,4-Bis(3-nitrofurazan-4-yl)furoxan Crystals in Dichloroethane or Benzene Mixture Solvents.

According to the experiments, DNTF crystallizes in benzene/methylbenzene (1:1), benzene/methylbenzene/ethanol (2:3:5), and sym-dichloroethane solvents into two similar crystal shapes, namely strip and tetrahedral. There is a possibility that solvent changes the crystal morphology. In order to explain this phenomenon, the DNTF growth interface model was constructed according to the actual solution environment. The interaction energy between the solvent phase and the DNTF crystal face was studied by means of molecular dynamics simulation. The crystal morphology of DNTF was predicted using the classical modified attachment energy model (MAE) in benzene, methylbenzene, benzene/methylbenzene (1:1), benzene/methylbenzene/ethanol (2:3:5), and sym-dichloroethane. The results show that the DNTF growths are mainly dominated by the (011), (001), (101), (110), (111), and (11[Formula: see text]) faces in vacuum. However, only a few faces will remain in the solvents, of which the (011) and (101) faces are exposed in benzene, methylbenzene, and benzene/methylbenzene (1:1), and only the (111) faces constitute the crystal shape of the DNTF in benzene/methylbenzene/ethanol (2:3:5) and sym-dichloroethane. The predicted results successfully explained the observed phenomena in the experiment. The simulation results can provide some guidance for the crystallization process of DNTF.