Tuning energetic properties through co-crystallisation - a high-pressure experimental and computational study of nitrotriazolone: 4,4'-bipyridine.

We report the preparation of a co-crystal formed between the energetic molecule 3-nitro-1,2,4-triazol-5-one (NTO) and 4,4'-bipyridine (BIPY), that has been structurally characterised by high-pressure single crystal and neutron powder diffraction data up to 5.93 GPa. No phase transitions or proton transfer were observed up to this pressure. At higher pressures the crystal quality degraded and the X-ray diffraction patterns showed severe twinning, with the appearance of multiple crystalline domains. Computational modelling indicates that the colour changes observed on application of pressure can be attributed to compression of the unit cell that cause heightened band dispersion and band gap narrowing that coincides with a shortening of the BIPY π⋯π stacking distance. Modelling also suggests that the application of pressure induces proton migration along an N-H⋯N intermolecular hydrogen bond. Impact-sensitivity measurements show that the co-crystal is less sensitive to initiation than NTO, whereas computational modelling suggests that the impact sensitivities of NTO and the co-crystal are broadly similar.

Accounting for super-, plateau- and mesa-rate burning by lead and copper-based ballistic modifiers in double-base propellants: a computational study.

We present a first-principles computational study to understand the action of lead and copper-based ballistic modifiers in the combustion of double-base propellants (DBPs). We show that lead oxide clusters are easily broken down upon addition of small amounts of carbon and the resulting graphitic matrix, dispersed with weakly bound and exposed Pb sites, acts as a Lewis acid to bind small molecule Lewis bases such as NO2 and CH2O that form in the combustion flame. This accounts for super-rate burning, where the fuel burn rate is enhanced. We also show how carbon availability accounts for the plateau- and mesa-rate burning effects, where the fuel burn rate is suppressed. In contrast, cluster integrity on binding carbon to copper oxide is retained, and interaction with NO2 and CH2O is essentially negligible. Carbon binds more strongly to copper oxide, however, and we therefore propose that when carbon levels start to fall this results in the lead oxide clusters being starved of carbon, which leads to plateau and mesa burning. Taken together, the calculations support a general model that accounts for the super-, plateau- and mesa-rate ballistic modifier burning effects.

High-pressure structural studies and pressure-induced sensitisation of 3,4,5-trinitro-1H-pyrazole.

Herein we report the first high-pressure study of the energetic material 3,4,5-trinitro-1H-pyrazole (3,4,5-TNP) using neutron powder diffraction and single-crystal X-ray diffraction. A new high-pressure phase, termed Form II, was first identified through a substantial change in the neutron powder diffraction patterns recorded over the range 4.6-5.3 GPa, and was characterised further by compression of a single crystal to 5.3 GPa in a diamond-anvil cell using X-ray diffraction. 3,4,5-TNP was found to be sensitive to initiation under pressure, as demonstrated by its unexpected and violent decomposition at elevated pressures in successive powder diffraction experiments. Initiation coincided with the sluggish phase transition from Form I to Form II. Using a vibrational up-pumping model, its increased sensitivity under pressure can be explained by pressure-induced mode hardening. These findings have potential implications for the safe handling of 3,4,5-TNP, on the basis that shock- or pressure-loading may lead to significantly increased sensitivity to initiation.

Is the impact sensitivity of RDX polymorph dependent?

Impact-sensitivity predictions based on the vibrational up-pumping model show a strong polymorph dependency for RDX and highlight that one of the high-pressure forms, which forms during shock-wave experiments, is appreciably more susceptible to mechanical initiation. The origin of the predicted impact sensitivity variation can be attributed to vibrational mode hardening by pressure and to differences in the molecular conformation of RDX in the four polymorphs studied. These polymorphs present different distributions of molecular vibrations within their respective up-pumping windows, which leads to their varying ability to up-pump and trap the vibrational energy that arises from mechanical insult.

Towards understanding the catalytic properties of lead-based ballistic modifiers in double base propellants.

Lead-based ballistic modifiers are additives in double base propellants (DBPs) which render the burn rate insensitive to changes in pressure within a defined pressure range, thus imparting greater control of combustion conditions. In-coming European legislation will soon ban the use of lead in propellant formulations, however, and few suitable candidate replacement materials are currently available. In an effort to understand better the unique properties offered by lead-based modifiers, we present a first-principles computational study on Pb, PbO, PbO2, SnO2 and Bi2O3, all of which have been investigated experimentally as ballistic modifier materials. Our study demonstrates that various quantifiable properties exist for the lead-based materials. Overall, they have narrower electronic band gaps, lower surface energies and lower surface work functions than the lead-free systems, indicating a greater propensity to form stable chemical surfaces with higher catalytic activity. We also show that of the set, only Pb and α-PbO can support the formation of a weakly bound layer of amorphous carbon, a key experimental observable in the burning of DBPs.

Experimental and simulation study of the high-pressure behavior of squalane and poly-α-olefins.

The equation of state, dynamical properties, and molecular-scale structure of squalane and mixtures of poly-α-olefins at room temperature are studied with a combination of state-of-the-art, high-pressure experiments and molecular-dynamics simulations. Diamond-anvil cell experiments indicate that both materials are non-hydrostatic media at pressures above ∼1 GPa. The equation of state does not exhibit any sign of a first-order phase transition. High-pressure x-ray diffraction experiments on squalane show that there are no Bragg peaks, and hence, the apparent solidification occurs without crystallization. These observations are complemented by a survey of the equation of state and dynamical properties using simulations. The results show that molecular diffusion is essentially arrested above about 1 GPa, which supports the hypothesis that the samples are kinetically trapped in metastable amorphous-solid states. The shear viscosity becomes extremely large at very high pressures, and the coefficient governing its increase from ambient pressure is in good agreement with the available literature data. Finally, simulated radial distribution functions are used to explore the evolution of the molecular-scale structure with increasing pressure. Subtle changes in the short-range real-space correlations are related to a collapse of the molecular conformations with increasing pressure, while the evolution of the static structure factor shows excellent correlation with the available x-ray diffraction data. These results are of indirect relevance to oil-based lubricants, as the pressures involved are comparable to those found in engines, and hence, the ability of lubricating thin films to act as load-bearing media can be linked to the solidification phenomena studied in this work.

Triglycerides as Novel Phase-Change Materials: A Review and Assessment of Their Thermal Properties.

Latent Heat Storage (LHS) with Phase-Change Materials (PCMs) represents a high energy density storage technology which could be applied in a variety of applications such as waste heat recovery and integration of renewable energy technologies in energy systems. To increase the sustainability of these storage solutions, PCMs have to be developed with particular regard to bio-origin and biodegradability. Triglycerides represent an interesting class of esters as the main constituents of animal and vegetable fats, with attractive thermal properties. In order to be used as PCMs, the thermal behaviour of triglycerides has to be fully understood, as in some cases they have been reported to show polymorphism and supercooling. This study assesses the suitability of triglycerides as PCMs by reviewing the literature published so far on their behaviour and properties. In particular, melting points, enthalpies of fusion, polymorphism, thermal conductivities, heat capacities and thermal cycling stabilities are considered, with a focus on LHS and thermal energy storage applications. In addition, the efforts conducted regarding modelling and the prediction of melting points and enthalpies based on chemical structures are summarized and assessed.

A Prolific Solvate Former, Galunisertib, under the Pressure of Crystal Structure Prediction, Produces Ten Diverse Polymorphs.

The solid form screening of galunisertib produced many solvates, prompting an extensive investigation into possible risks to the development of the favored monohydrate form. Inspired by crystal structure prediction, the search for neat polymorphs was expanded to an unusual range of experiments, including melt crystallization under pressure, to work around solvate formation and the thermal instability of the molecule. Ten polymorphs of galunisertib were found; however, the structure predicted to be the most stable has yet to be obtained. We present the crystal structures of all ten unsolvated polymorphs of galunisertib, showing how state-of-the-art characterization methods can be combined with emerging computational modeling techniques to produce a complete structure landscape and assess the risk of late-appearing, more stable polymorphs. The exceptional conformational polymorphism of this prolific solvate former invites further development of methods, computational and experimental, that are applicable to larger, flexible molecules with complex solid form landscapes.

Vibrationally induced metallisation of the energetic azide α-NaN3.

As initiation of an energetic material requires rupture of a covalent bond, and therefore population of antibonding electronic states, consideration of the electronic band gap has dominated initiation mechanisms for solid state materials. Most prominent are models based on metallisation, where static mechanical perturbation leads to closing of the electronic band gap. This work explores an alternative mechanism for the dynamic metallisation of a model energetic material, where vibrational excitation resulting from mechanical impact is found to induce transient metallisation of α-NaN3. The normal coordinates associated with bending the azido anion close the electronic band gap, facilitating the formation of highly reactive species important for initiation of energetic materials. The DFT simulated vibrational spectrum of α-NaN3 exhibits excellent reproduction of the experimental low-temperature inelastic neutron scattering spectrum (INS).

Structure and properties of nitrocellulose: approaching 200 years of research.

This review brings together almost 200 years of fragmented research on the structure of nitrocellulose to give an overview that covers production to application in composite materials. As a mouldable plastic, energetic rocket propellant and biomolecular binding membrane, nitrocellulose still finds widespread practical application today despite the inception of synthetic plastics. The influence of different cellulose source materials affects the structure and properties of nitrocellulose in ways that are not fully understood, and so this review brings together relatively recent developments in the understanding of cellulose nanostructures to highlight where the gaps in understanding now reside. The influence of nitration conditions on the material properties of nitrocellulose is described, together with the proposed mechanisms and equilibria associated with these synthetic routes. The reported crystal structures of nitrocellulose are also reviewed, and the confirmed structural features are separated from those yet to be proven. We also consider how nitrocellulose interacts with other compounds, to help explain the distinct properties of its composite materials. This review points to further work that is required to obtain well founded structural models of nitrocellulose, while highlighting opportunities to control and direct its structure to improve its material properties.

Temperature dependent solid-state proton migration in dimethylurea-oxalic acid complexes.

The phenomenon of solid-state proton migration within molecular complexes containing short hydrogen bonds is investigated in two dimethylurea-oxalic acid complexes. Extensive characterisation by both X-ray and neutron diffraction shows that proton migration along the hydrogen bond can be induced in these complexes as a function of temperature. This emphasises the subtle features of the hydrogen bond potential well in such short hydrogen bonded complexes, both intrinsically and in the effect of the local crystalline environment. Based on these findings, the synthesis and analysis of a series of solid-state molecular complexes is shown to be a potential route to designing materials with tuneable proton migration effects.

High-pressure reversibility in a plastically flexible coordination polymer crystal.

Single crystals which exhibit mechanical flexibility are promising materials for advanced technological applications. Before such materials can be used, a detailed understanding of the mechanisms of bending is needed. Using single crystal X-ray diffraction and microfocus Raman spectroscopy, we study in atomic detail the high-pressure response of the plastically flexible coordination polymer [Zn(µ-Cl)2(3,5-dichloropyridine)2]n (1). Contradictory to three-point bending, quasi-hydrostatic compression of (1) is completely reversible, even following compression to over 9 GPa. A structural phase transition is observed at ca. 5 GPa. DFT calculations show this transition to result from the pressure-induced softening of low-frequency vibrations. This phase transition is not observed during three-point-bending. Microfocus synchrotron X-ray diffraction revealed that bending yields significant mosaicity, as opposed to compression. Hence, our studies indicate of overall disparate mechanical responses of bulk flexibility and quasi-hydrostatic compression within the same crystal lattice. We suspect this to be a general feature of plastically bendable materials.

Predicting the impact sensitivity of a polymorphic high explosive: the curious case of FOX-7.

The impact sensitivity (IS) of FOX-7 polymorphs is predicted by phonon up-pumping to decrease as layers of FOX-7 molecules flatten. Experimental validation proved anomalous owing to a phase transition during testing, raising questions regarding impact sensitivity measurement and highlighting the need for models to predict IS of polymorphic energetic materials.

Towards Computational Screening for New Energetic Molecules: Calculation of Heat of Formation and Determination of Bond Strengths by Local Mode Analysis.

The reliable determination of gas-phase and solid-state heats of formation are important considerations in energetic materials research. Herein, the ability of PM7 to calculate the gas-phase heats of formation for CNHO-only and inorganic compounds has been critically evaluated, and for the former, comparisons drawn with isodesmic equations and atom equivalence methods. Routes to obtain solid-state heats of formation for a range of single-component molecular solids, salts, and co-crystals were also evaluated. Finally, local vibrational mode analysis has been used to calculate bond length/force constant curves for seven different chemical bonds occurring in CHNO-containing molecules, which allow for rapid identification of the weakest bond, opening up great potential to rationalise decomposition pathways. Both metrics are important tools in rationalising the design of new energetic materials through computational screening processes.

Synthesis and characterization of the mixed-ligand coordination polymer Cu3Cl(N4C-NO2)2.

Reduction of copper(ii) chloride using sodium ascorbate in the presence of pure sodium 5-nitro-tetrazolate (NaNT) forms copper(i) 5-nitrotetrazolate - a known initiatory explosive (DBX-1) - and the novel mixed-ligand copper(i) chloride 5-nitrotetrazolate coordination polymer Cu3Cl(N4C-NO2)2, as well as mixtures of both. The reaction is controlled by the presence of seed crystals and transition metal compounds other than CuCl2. Cu3Cl(N4C-NO2)2 is obtained as a wine-red, air stable, water-insoluble, crystalline and highly sensitive explosive material with a greater crystal density, lower thermal stability and a higher sensitivity toward hydrolysis and shock than DBX-1. Efforts to obtain the stable and pure starting material are improved by crystallisation of NaNT as a tetrahydrate. Cu3Cl(N4C-NO2)2 and Na(H2O)4(NT) were characterised by single crystal and powder XRD, IR spectroscopy, magnetic and thermal measaurements, elemental analysis, particle size measurements, mass spectrometry, and by drop weight testing.

Temperature- and pressure-induced proton transfer in the 1:1 adduct formed between squaric acid and 4,4'-bipyridine.

We have applied a combination of spectroscopic and diffraction methods to study the adduct formed between squaric acid and bypridine, which has been postulated to exhibit proton transfer associated with a single-crystal to single-crystal phase transition at ca. 450 K. A combination of X-ray single-crystal and very-high flux powder neutron diffraction data confirmed that a proton does transfer from the acid to the base in the high-temperature form. Powder X-ray diffraction measurements demonstrated that the transition was reversible but that a significant kinetic energy barrier must be overcome to revert to the original structure. Computational modeling is consistent with these results. Modeling also revealed that, while the proton transfer event would be strongly discouraged in the gas phase, it occurs in the solid state due to the increase in charge state of the molecular ions and their arrangement inside the lattice. The color change is attributed to a narrowing of the squaric acid to bipyridine charge-transfer energy gap. Finally, evidence for the possible existence of two further phases at high pressure is also presented.

The crystal structure of beta-RDX-an elusive form of an explosive revealed.

The crystal structure of the highly metastable beta-form of RDX shows that the molecules adopt different conformations compared to the alpha-form and that, contrary to previous reports, the beta-form obtained at ambient pressure is not the same form as that obtained at elevated temperatures and pressures.

A new polymorph of N,N'-dimethylurea characterized by X-ray diffraction and first-principles lattice dynamics calculations.

In this study we present a combined crystallographic and computational study of a new polymorph of N,N'-dimethylurea (DMU) with P2(1)2(1)2 space group symmetry, along with a revised theoretical study of the previously known phase in its corrected space group (Fdd2). X-ray diffraction studies show crystal structures that are very similar, differing only in the relative orientation of the hydrogen-bonded molecular chains that are common to both phases. The vibrational spectra were obtained from B3LYP hybrid functional lattice dynamics calculations and compared with the experimental data for the known phase. The free-energy difference between the forms is derived from the Gamma-point optical mode frequencies, and amounts to less than 1 kJ mol(-1) across the temperature range of interest. The electronic densities-of-states of both phases are also computed, yielding only marginal differences in valence and conduction band compositions and band gap widths. Taken together, the results highlight the small but important differences separating the two crystal lattices.

An acetonitrile-solvated cocrystal of piroxicam and succinic acid with co-existing zwitterionic and non-ionized piroxicam molecules.

This work reports a new acetonitrile (ACN)-solvated cocrystal of piroxicam (PRX) and succinic acid (SA), 2C15H13N3O4S·0.5C4H6O4·C2H3N or PRX:SA:ACN (4:1:2), which adopts the triclinic space group P-1. The outcome of crystallization from ACN solution can be controlled by varying only the PRX:SA ratio, with a higher PRX:SA ratio in solution unexpectedly favouring a lower stoichiometric ratio in the solid product. In the new solvate, zwitterionic (Z) and non-ionized (NI) PRX molecules co-exist in the asymmetric unit. In contrast, the nonsolvated PRX-SA cocrystal contains only NI-type PRX molecules. The ACN molecule entrapped in PRX-SA·ACN does not form any hydrogen bonds with the surrounding molecules. In the solvated cocrystal, Z-type molecules form dimers linked by intermolecular N-H...O hydrogen bonds, whereas every pair of NI-type molecules is linked to SA via N-H...O and O-H...N hydrogen bonds. Thermogravimetry and differential scanning calorimetry suggest that thermal desolvation of the solvate sample occurs at 148 °C, and is followed by recrystallization, presumably of a multicomponent PRX-SA structure. Vibrational spectra (IR and Raman spectroscopy) of PRX-SA·ACN and PRX-SA are also used to demonstrate the ability of spectroscopic techniques to distinguish between NI- and Z-type PRX molecules in the solid state. Hence, vibrational spectroscopy can be used to distinguish the PRX-SA cocrystal and its ACN solvate.

A novel crystal form of metacetamol: the first example of a hydrated form.

We report the crystal structure and crystallization conditions of a first hydrated form of metacetamol (a hemihydrate), C8H9NO2·0.5H2O. It crystallizes from metacetamol-saturated 1:1 (v/v) water-ethanol solutions in a monoclinic structure (space group P21/n) and contains eight metacetamol and four water molecules per unit cell. The conformations of the molecules are the same as in polymorph II of metacetamol, which ensures the formation of hydrogen-bonded dimers and R22(16) ring motifs in its crystal structure similar to those in polymorph II. Unlike in form II, however, these dimers in the hemihydrate are connected through water molecules into infinite hydrogen-bonded molecular chains. Different chains are linked to each other by metacetamol-water and metacetamol-metacetamol hydrogen bonds, the latter type being also present in polymorph I. The overall noncovalent network of the hemihydrate is well developed and several types of hydrogen bonds are responsible for its formation.