Structure of liquid nitromethane: Comparison of simulation and diffraction studies

Unique O═N...O Pnicogen Interactions in Nitromethane Dimers: Evidence Using Matrix Isolation Infrared Spectroscopy and Computational Methodology

Geometries of nitromethane homodimers have been revisited using ab initio and density functional methodologies, following their formation within cryogenic matrices, confirmed using infrared spectroscopy. In contrast to the claim that the intermolecular interactions are due to dispersion forces or very weak hydrogen bonds, in the present work, concrete evidence for the prevalence of O═N...O pnicogen bonding has been presented. The pnicogen bonds have been found to be primarily responsible for the characteristic geometry of a homodimer. The formation of a nitromethane dimer with pnicogen bonding stabilization is evidenced using matrix isolation infrared spectroscopy with Ne and Ar matrices and computations. The interactions within homodimers have been characterized using quantum theory of atoms in molecules, natural bond orbital, energy decomposition, electrostatic potential mapping, and noncovalent interaction analyses. The larger intermolecular separations in liquid nitromethane indicated by previous molecular dynamics-based studies alongside the supposed invariance of infrared spectra in the gas phase and liquid phase could have led to the assumption of a lack of intermolecular interactions. However, the prevalence of hydrogen and pnicogen bonds across these larger than usual separations is affirmed by the geometries presented here.

Intermolecular Interaction in Methylene Halide (CH2F2, CH2Cl2, CH2Br2 and CH2I2) Dimers

The intermolecular interaction in difluoromethane, dichloromethane, dibromomethane, and diiodomethane dimers has been investigated using high level quantum chemical methods. The potential energy curve of intermolecular interaction along the C⋯C bond distance obtained using the coupled-cluster theory with singles, doubles, and perturbative triples excitations CCSD(T) were compared with values given by the same method, but applying the local (LCCSD(T)) and the explicitly correlated (CCSD(T)-F12) approximations. The accuracy of other theoretical methods-Hartree-Fock (HF), second order Møller-Plesset perturbation (MP2), and dispersion corrected DFT theory-were also presented. In the case of MP2 level, the canonical and the local-correlation cases combined with the density-fitting technique (DF-LMP2)theories were considered, while for the dispersion-corrected DFT, the empirically-corrected BLYP-D and the M06-2Xexchange-correlation functionals were applied. In all cases, the aug-cc-pVTZ basis set was used, and the results were corrected for the basis set superposition error (BSSE) using the counterpoise method. For each molecular system, several dimer geometries were found, and their mutual orientations were compared with the nearest neighbor orientations obtained in recent neutron scattering studies. The nature of the intermolecular interaction energy was discussed.

Structural Rearrangement of Energetic Materials under an External Electric Field: A Case Study of Nitromethane

As a significant stimulus, the external electric field (EEF) can change the decomposition mechanism and energy release of energetic materials (EMs). Hence, understanding the response of EMs to an EEF is greatly meaningful for their safe usage. Herein, the structural arrangement, a crucial factor in the impact sensitivity and detonation performance of EMs, under the EEF ranging from 0.0 to 0.5 V/Å was investigated via molecular dynamics simulation. Nitromethane (NM) was taken as a case study due to the simple structure. The simulation results show that there exists a critical EEF strength between 0.2 and 0.3 V/Å, which can induce the transition of NM molecules from relatively disordered distribution to solidlike ordered and compacted arrangement with a large density. In this ordered structure, NM dipoles are aligned in a head-to-tail pattern parallel to the EEF direction because of the favored dipole-dipole interactions and weak C-H···O hydrogen bonds. As the EEF strength is enhanced, the potential energy and cohesive energy density of the NM system gradually decrease and increase, respectively, indicative of high thermodynamics stability of ordered arrangement. The results reported here also shed light on the potential of the EEF to induce the nucleation and crystallization to explore new polymorphs of EMs.

Neutron Scattering of Aromatic and Aliphatic Liquids

Organic solvents, such as cyclohexane, cyclohexene, methylcyclohexane, benzene and toluene, are widely used as both reagents and solvents in industrial processes. Despite the ubiquity of these liquids, the local structures that govern the chemical properties have not been studied extensively. Herein, we report neutron diffraction measurements on liquid cyclohexane, cyclohexene, methylcyclohexane, benzene and toluene at 298 K to obtain a detailed description of the local structure in these compounds. The radial distribution functions of the centres of the molecules, as well as the partial distribution functions for the double bond for cyclohexene and methyl group for methylcyclohexane and toluene have been calculated. Additionally, probability density functions and angular radial distribution functions were extracted to provide a full description of the local structure within the chosen liquids. Structural motifs are discussed and compared for all liquids, referring specifically to the functional group and aromaticity present in the different liquids.

Theoretical gas to liquid shift of (15)N isotropic nuclear magnetic shielding in nitromethane using ab initio molecular dynamics and GIAO/GIPAW calculations

Chemical shift requires the knowledge of both the sample and a reference magnetic shielding. In few cases as nitrogen (15N), the standard experimental reference corresponds to its liquid phase. Theoretical estimate of NMR magnetic shielding parameters of compounds in their liquid phase is then mandatory but usually replaced by an easily-get gas phase value, forbidding direct comparisons with experiments. We propose here to combine ab initio molecular dynamic simulations with the calculations of magnetic shielding using GIAO approach on extracted cluster's structures from MD. Using several computational strategies, we manage to accurately calculate 15N magnetic shielding of nitromethane in its liquid phase. Theoretical comparison between liquid and gas phase allows us to extrapolate an experimental value for the 15N magnetic shielding of nitromethane in gas phase between -121.8 and -120.8 ppm.

Exploding Nitromethane in Silico, in Real Time

Nitromethane (NM) is widely applied in chemical technology as a solvent for extraction, cleaning, and chemical synthesis. NM was considered safe for a long time, until a railroad tanker car exploded in 1958. We investigate the detonation kinetics and explosion reaction mechanisms in a variety of systems consisting of NM, molecular oxygen, and water vapor. Reactive molecular dynamics allows us to simulate reactions in time-domain, as they occur in real life. High polarity of the NM molecule is shown to play a key role, driving the first exothermic step of the reaction. Rapid temperature and pressure growth stimulate the subsequent reaction steps. Oxygen is important for faster oxidation, whereas its optimal concentration is in agreement with the proposed reaction mechanism. Addition of water (50 mol %) inhibits detonation; however, water does not prevent detonation entirely. The reported results provide important insights for improving applications of NM and preserving the safety of industrial processes.

Self-enhanced catalytic activities of functionalized graphene sheets in the combustion of nitromethane: molecular dynamic simulations by molecular reactive force field

Functionalized graphene sheet (FGS) is a promising additive that enhances fuel/propellant combustion, and the determination of its mechanism has attracted much interest. In the present study, a series of molecular dynamic simulations based on a reactive force field (ReaxFF) are performed to explore the catalytic activity (CA) of FGS in the thermal decay of nitromethane (NM, CH3NO2). FGSs and pristine graphene sheets (GSs) are oxidized in hot NM liquid to increase their functionalities and subsequently show self-enhanced CAs during the decay. The CAs result from the interatomic exchanges between the functional groups on the sheets and the NM liquid, i.e., mainly between H and O atoms. CA is dependent on the density of NM, functionalities of sheets, and temperature. The GSs and FGSs that originally exhibit different functionalities tend to possess similar functionalities and consequently similar CAs as temperature increases. Other carbon materials and their oxides can accelerate combustion of other fuels/propellants similar to NM, provided that they can be dispersed and their key reaction steps in combustion are similar to NM.

Formation of supramolecular structures in organic solvents

The unusual interactions in organic liquids such as methane derivatives, arenes, and alkanes by the infrared (IR) method were revealed. The transformations of molecular shapes, arising from nonclassical hydrogen and dihydrogen bonds, as well as water factor provide the existence of supramolecular structures in organic fluids. The interpretation of the obtained results in terms of the quantum-chemical calculations has been suggested.

A gold cyano complex in nitromethane: MD simulation and X-ray diffraction

The solvation structure around the dicyanoaurate(I) anion (Au(CN)2-) in a dilute nitromethane (CH3NO2) solution is presented from X-ray diffraction measurements and molecular dynamics simulation (NVT ensemble, 460 nitromethane molecules at room temperature). The simulations are based on a new solute-solvent force-field fitted to a training set of quantum-chemically derived interaction energies. Radial distribution functions from experiment and simulation are in good agreement. The solvation structure has been further elucidated from MD data. Several shells can be identified. We obtain a solvation number of 13-17 nitromethane molecules with a strong preference to be oriented with their methyl groups towards the solute.

MOLECULAR DYNAMIC SIMULATIONS OF SOLID NITROMETHANE UNDER HIGH PRESSURES

We report ReaxFF molecular dynamic simulations of structure change of crystalline nitromethane and the formation of hydrogen bond under high pressure. Under high pressure, the angles between C–N bonds and X, Y and Z axes have changed. Through the calculation of g(r) of O and H atoms, we found a new peak near 1.6 Å, which indicates the formation of the hydrogen bond between O and H atoms. We calculated the distribution of the angles of the C–N bonds orientations, the distribution of the dihedral angle of CNOO , and the charge distribution of nitromethane molecules under various pressures, and made a comparison between low and high pressures. The effects of hydrogen bonding in high explosive materials are discussed.

Solution structure of NaNO3 in water: diffraction and molecular dynamics simulation study

The structure of a series of aqueous sodium nitrate solutions (1.9-7.6 M) was studied using a combination of experimental and theoretical methods. The results obtained from diffraction (X-ray, neutron) and molecular dynamics simulation have been compared and the capabilities and limitations of the methods in describing solution structure are discussed. For the solutions studied, diffraction methods were found to perform very well in description of hydration spheres of the sodium ion but do not yield detailed structural information on the anion's hydration structure. Molecular dynamics simulations proved to be a suitable tool in the detailed interpretation of the hydration sphere of ions, ion pair formation, and bulk structure of solutions.

Complete structural characterization of metallacyclic complexes in solution-phase using simultaneously X-ray diffraction and molecular dynamics simulation

Wide-angle X-ray diffraction and molecular dynamics simulation has been used to perform complete structural characterization of nitromethane solution of a 16-membered gold(I) ring. The joint application of these two methods was an adequate tool to describe not only the structure of the complex but also the solvation properties of the complex in nitromethane and the effect of the solvation on the bulk structure. It has been found that a relatively diffuse slightly distorted solvation shell is formed around the complex, following the shape of the molecule. Nitromethane molecules in the solvation sphere are distributed randomly; no special orientation can be detected. The interaction energy of the complex with nitromethane molecules is attractive. In bulk, besides the antiparallel orientation of the nitromethane molecules, T-shape orientation and long-range order in antidipole orientation can also be detected.