Siddharth Dasgupta

BENZENE FORMS HYDROGEN BONDS WITH WATER

Fully rotationally resolved spectra of three isotopic species of 1:1 clusters of benzene with water (H2O, D2O, and HDO) were fit to yield moments of inertia that demonstrate unambiguously that water is positioned above the benzene plane in nearly free internal rotation with both hydrogen atoms pointing toward the π cloud. Ab initio calculations (MP2 level of electron correlation and 6-31 G** basis set with basis set superposition error corrections) predict a binding energy De ≳ 1.78 kilocalories per mole. In both the experimental and theoretical structures, water is situated nearly 1 angstrom within the van der Waals contacts of the monomers, a clear manifestation of hydrogen bond formation in this simple model of aqueous-π electron interactions.

Initiation Mechanisms and Kinetics of Pyrolysis and Combustion of JP-10 Hydrocarbon Jet Fuel

In order to investigate the initiation mechanisms and kinetics associated with the pyrolysis of JP-10 (exo-tricyclo[5.2.1.0(2,6)]decane), a single-component hydrocarbon jet fuel, we carried out molecular dynamics (MD) simulations employing the ReaxFF reactive force field. We found that the primary decomposition reactions involve either (1) dissociation of ethylene from JP-10, resulting in the formation of a C(8) hydrocarbon intermediate, or (2) the production of two C(5) hydrocarbons. ReaxFF MD leads to good agreement with experiment for the product distribution as a function of temperature. On the basis of the rate of consumption of JP-10, we calculate an activation energy of 58.4 kcal/mol for the thermal decomposition of this material, which is consistent with a strain-facilitated C-C bond cleavage mechanism in JP-10. This compares well with the experimental value of 62.4 kcal/mol. In addition, we carried out ReaxFF MD studies of the reactive events responsible for oxidation of JP-10. Here we found overall agreement between the thermodynamic energies obtained from ReaxFF and quantum-mechanical calculations, illustrating the usefulness of ReaxFF for studying oxidation of hydrocarbons. The agreement of these results with available experimental observations demonstrates that ReaxFF can provide useful insights into the complicated thermal decomposition and oxidation processes of important hydrocarbon fuels.

Carbon Cluster Formation during Thermal Decomposition of Octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazocine and 1,3,5-Triamino-2,4,6-trinitrobenzene High Explosives from ReaxFF Reactive Molecular Dynamics Simulations

We report molecular dynamics (MD) simulations using the first-principles-based ReaxFF reactive force field to study the thermal decomposition of 1,3,5-triamino-2,4,6-trinitrobenzene (TATB) and octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazocine (HMX) at various densities and temperatures. TATB is known to produce a large amount (15-30%) of high-molecular-weight carbon clusters, whereas detonation of nitramines such as HMX and RDX (1,3,5-trinitroperhydro-1,3,5-triazine) generate predominantly low-molecular-weight products. In agreement with experimental observation, these simulations predict that TATB decomposition quickly (by 30 ps) initiates the formation of large carbonaceous clusters (more than 4000 amu, or approximately 15-30% of the total system mass), and HMX decomposition leads almost exclusively to small-molecule products. We find that HMX decomposes readily on this time scale at lower temperatures, for which the decomposition rate of TATB is about an order of magnitude slower. Analyzing the ReaxFF MD results leads to the detailed atomistic structure of this carbon-rich phase of TATB and allows characterization of the kinetics and chemistry related to this phase and their dependence on system density and temperature. The carbon-rich phase formed from TATB contains mainly polyaromatic rings with large oxygen content, leading to graphitic regions. We use these results to describe the initial reaction steps of thermal decomposition of HMX and TATB in terms of the rates for forming primary and secondary products, allowing comparison to experimentally derived models. These studies show that MD using the ReaxFF reactive force field provides detailed atomistic information that explains such macroscopic observations as the dramatic difference in carbon cluster formation between TATB and HMX. This shows that ReaxFF MD captures the fundamental differences in the mechanisms of such systems and illustrates how the ReaxFF may be applied to model complex chemical phenomena in energetic materials. The studies here illustrate this for modestly sized systems and modest periods; however, ReaxFF calculations of reactive processes have already been reported on systems with approximately 10(6) atoms. Thus, with suitable computational facilities, one can study the atomistic level chemical processes in complex systems under extreme conditions.

Hessian-biased force fields from combining theory and experiment

We describe a new approach of combining experimental and theoretical information to develop accurate valence force fields. We combine the Hessian from ab initio calculations with the structural and spectroscopic data from experiment to generate a new Hessian that is used for extracting force fields. Emphasis is placed on obtaining accurate geometries and accurate frequencies. This technique is illustrated by calculations on formaldehyde, formate anion, thioformaldehyde, carbonyl chloride, and carbonyl fluoride.

Explanation of the Colossal Detonation Sensitivity of Silicon Pentaerythritol Tetranitrate (Si-PETN) Explosive

DFT calculations have identified the novel rearrangement shown here for decomposition of the Si derivative of the PETN explosive [PentaErythritol TetraNitrate (PETN), C(CH(2)ONO(2))(4)] that explains the very dramatic increase in sensitivity observed experimentally. The critical difference is that Si-PETN allows a favorable five-coordinate transition state in which the new Si-O and C-O bonds form simultaneously, leading to a transition state barrier of 33 kcal/mol (it is 80 kcal/mol for PETN) and much lower than the normal O-NO(2) bond fission observed in other energetic materials (approximately 40 kcal/mol). In addition this new mechanism is very exothermic (45 kcal/mol) leading to a large net energy release at the very early stages of Si-PETN decomposition that facilitates a rapid temperature increase and expansion of the reaction zone.

A detailed model for the decomposition of nitramines: RDX and HMX

A unified decomposition scheme for two very important cyclic nitramines used as primary explosives – RDX and HMX – has been constructed using ab initio Density Functional Theory (DFT) calculations. Molecular parameters such as vibrational frequencies and moments of inertia corresponding to the computed potential energy profile of unimolecular decomposition of these nitramines were then used to obtain the thermochemistry of all identified species and reaction rate constants of each individual channel. These primary decomposition reactions were then combined with: (i) important secondary reactions of the key reactive radical intermediates, such as CH2NNO2 (Methylene Nitramine MN), CH2N, NO, NO2, OH, etc.; (ii) existing nitramie reaction networks [33]. We have developed an improved mechanism for the detailed chemistry of nitramines which can be applied to combustion and detonation phenomena of this class of energetic materials.

Strategies for Multiscale Modeling and Simulation of Organic Materials: Polymers and Biopolymers

Abstract Advances in theory and methods are making it practical to consider fully first principles (de novo) predictions of structures, properties and processes for organic materials. However, despite the progress there remains an enormous challenge in bridging the vast range of distances and time scales between de novo atomistic simulations and the quantitative continuum models for the macroscopic systems essential in industrial design and operations. Recent advances relevant to such developments include: quantum chemistry including continuum solvation and force field embedding, de novo force fields to describe phase transitions, molecular dynamics (MD) including continuum solvent, non equilibrium MD for rheology and thermal conductivity and mesoscale simulations. To provide some flavor for the opportunities we will illustrate some of the progress and challenges by summarizing some recent developments in methods and their applications to polymers and biopolymers. Four different topics will be covered: (1) hierarchical modeling approach applied to modeling olfactory receptors, (2) stabilization of leucine zipper coils by introduction of trifluoroleucine, (3) modeling response of polymers sensors for electronic nose, and (4) diffusion of gases in amorphous polymers.

The Hessian biased singular value decomposition method for optimization and analysis of force fields

We present methodology (HBFF/SVD) for optimizing the form and parameters of force fields (FF) for molecular dynamics simulations through utilizing information about properties such as the geometry, Hessian, polarizability, stress (crystals), and elastic constants (crystals). This method is based on singular value decomposition (SVD) of the Jacobian describing the partial derivatives in various properties with respect to FF parameters. HBFF/SVD is effective for optimizing the parameters for accurate FFs of organic, inorganic, and transition metal compounds. In addition it provides information on the validity of the functional form of the FF for describing the properties of interest. This method is illustrated by application to organic molecules (CH2O, C2H4, C4H6, C6H8, C6H6, and naphthalene) and inorganic molecules (Cl2CrO2 and Cl2MoO2).

Electronic Structures of Group 9 Metallocorroles with Axial Ammines

The electronic structures of metallocorroles (tpfc)M(NH(3))(2) and (tfc)M(NH(3))(2) (tpfc is the trianion of 5,10,15-(tris)pentafluorophenylcorrole, tfc is the trianion of 5,10,15-trifluorocorrole, and M = Co, Rh, Ir) have been computed using first principles quantum mechanics [B3LYP flavor of Density Functional Theory (DFT) with Poisson-Boltzmann continuum solvation]. The geometry was optimized for both the neutral systems (formal M(III) oxidation state) and the one-electron oxidized systems (formally M(IV)). As expected, the M(III) systems have a closed shell d(6) configuration; for all three metals, the one-electron oxidation was calculated to occur from a ligand-based orbital (highest occupied molecular orbital (HOMO) of B(1) symmetry). The ground state of the formal M(IV) system has M(III)-Cπ character, indicating that the metal remains d(6), with the hole in the corrole π system. As a result the calculated M(IV/III) reduction potentials are quite similar (0.64, 0.67, and 0.56 V vs SCE for M = Ir, Rh and Co, respectively), whereas the differences would have been large for purely metal-based oxidations. Vertically excited states with substantial metal character are well separated from the ground state in one-electron-oxidized cobalt (0.27 eV) and rhodium (0.24 eV) corroles, but become closer in energy in the iridium (0.15 eV) analogues. The exact splittings depend on the chosen functional and basis set combination and vary by ~0.1 eV.

Thermal Decomposition of Energetic Materials by ReaxFF Reactive Molecular Dynamics

We report the study of thermal decomposition of 1,3,5‐trinitrohexahydro‐s‐triazine (RDX) bonded with polyurethane (Estane) and of the bulk hydrazine by molecular dynamics (MD) simulations equipped with the reactive force field (ReaxFF). For the polymer binder explosive, the simulation results show that the thermal decomposition of RDX is affected by the presence of the polymer binder Estane. Generally, with addition of Estane the decomposition of RDX slows down. Final products including N2, H2O, CO, CO2 and intermediates NO2, NO and HONO have been identified from the thermal decomposition processes. For the bulk hydrazine, it is found that with the increase of temperature, its decomposition increases and more N2 and H2 are generated, but NH3 molecules are consumed much faster at higher temperatures. This simulation work provides us an approach to quickly test the response of various energetic materials to thermal conditions.