Kinetics of the γ–δ phase transition in energetic nitramine-octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazocine

The solid, secondary explosive nitramine-octahydro-1,3,5,7-tetranitro-1,3,5,7 or HMX has four different stable polymorphs which have different molecular conformations, crystalline structures, and densities, making structural phase transitions between these nontrivial. Previous studies of the kinetics of the β–δ HMX structural transition found this to happen by a nucleation and growth mechanism, where growth was governed by the heat of fusion, or melting, even though the phase transition temperature is more than 100 K below the melting point. A theory known as virtual melting could easily justify this since the large volume difference in the two phases creates a strain at their interface that can lower the melting point to the phase transition temperature through a relaxation of the elastic energy. To learn more about structural phase transitions in organic crystalline solids and virtual melting, here we use time-resolved X-ray diffraction to study another structural phase transition in HMX, γ–δ. Again, second order kinetics are observed which fit to the same nucleation and growth model associated with growth by melting even though the volume change in this transition is too small to lower the melting point by interfacial strain. To account for this, we present a more general model illustrating that melting over a very thin layer at the interface between the two phases reduces the total interfacial energy and is therefore thermodynamically favorable and can drive the structural phase transition in the absence of large volume changes. Our work supports the idea that virtual melting may be a more generally applicable mechanism for structural phase transitions in organic crystalline solids

Strain-induced phase transformation under compression in a diamond anvil cell: Simulations of a sample and gasket

Combined high pressure phase transformations (PTs) and plastic flow in a sample within a gasket compressed in diamond anvil cell (DAC) are studied for the first time using finite element method. The key point is that phase transformations are modelled as strain-induced, which involves a completely different kinetic description than for traditional pressure-induced PTs. The model takes into account, contact sliding with Coulomb and plastic friction at the boundaries between the sample, gasket, and anvil. A comprehensive computational study of the effects of the kinetic parameter, ratio of the yield strengths of high and low-pressure phases and the gasket, sample radius, and initial thickness on the PTs and plastic flow is performed. A new sliding mechanism at the contact line between the sample, gasket, and anvil called extrusion-based pseudoslip is revealed, which plays an important part in producing high pressure. Strain-controlled kinetics explains why experimentally determined phase transformation pressure and kinetics (concentration of high pressure phase vs. pressure) differ for different geometries and properties of the gasket and the sample: they provide different plastic strain, which was not measured. Utilization of the gasket changes radial plastic flow toward the center of a sample, which leads to high quasi-homogeneous pressure for some geometries. For transformation to a stronger high pressure phase, plastic strain and concentration of a high-pressure phase are also quasi-homogeneous. This allowed us to suggest a method of determining strain-controlled kinetics from experimentation, which is not possible for weaker and equal-strength high-pressure phases and cases without a gasket. Some experimental phenomena are reproduced and interpreted. Developed methods and obtained results represent essential progress toward the understanding of PTs under compression in the DAC. This will allow one optimal design of experiments and conditions for synthesis of new high pressure phases

Plastic flows and phase transformations in materials under compression in diamond anvil cell: Effect of contact sliding

Modeling of coupled plastic flows and strain-induced phase transformations (PTs) under high pressure in a diamond anvil cell is performed with the focus on the effect of the contact sliding between sample and anvils. Finite element software ABAQUS is utilized and a combination of Coulomb friction and plastic friction is considered. Results are obtained for PTs to weaker, equal-strength, and stronger high pressure phases, using different scaling parameters in a strain-controlled kinetic equation, and with various friction coefficients. Compared to the model with cohesion, artificial shear banding near the constant surface is eliminated. Sliding and the reduction in friction coefficient intensify radial plastic flow in the entire sample (excluding a narrow region near the contact surface) and a reduction in thickness. A reduction in the frictioncoefficient to 0.1 intensifies sliding and increases pressure in the central region. Increases in both plastic strain and pressure lead to intensification of strain-induced PT. The effect of self-locking of sliding is revealed. Multiple experimental phenomena are reproduced and interpreted. Thus, plastic flow and PT can be controlled by controlling friction

Lattice Instability during Solid-Solid Structural Transformations under a General Applied Stress Tensor: Example of Si I → Si II with Metallization

The density functional theory was employed to study the stress-strain behavior and elastic instabilities during the solid-solid phase transformation (PT) when subjected to a general stress tensor, as exemplified for semiconducting Si I and metallic Si II, where metallization precedes the PT, so stressed Si I can be a metal. The hydrostatic PT occurs at 76 GPa, while under uniaxial loading it is 11 GPa (3.7 GPa mean pressure), 21 times lower. The Si I → Si II PT is described by a critical value of the phase-field’s modified transformation work, and the PT criterion has only two parameters given six independent stress elements. Our findings reveal novel, more practical synthesis routes for new or known high-pressure phases under predictable nonhydrostatic loading, where competition of instabilities can serve for phase selection rather than free energy minima used for equilibrium processing

Coupled phase transformations and plastic flows under torsion at high pressure in rotational diamond anvil cell: Effect of contact sliding

A three-dimensional large-sliding contact model coupled with strain-induced phase transformations (PTs) and plastic flow in a disk-like sample under torsion at high pressure in rotational diamond anvil cell (RDAC) is formulated and studied. Coulomb and plastic friction are combined and take into account variable parameters due to PT. Results are obtained for weaker, equal-strength, and stronger high pressure phases, and for three values of the kinetic coefficient in a strain-controlled kinetic equation and friction coefficient. All drawbacks typical of problem with cohesion are overcome, including eliminating mesh-dependent shear band and artificial plastic zones. Contact sliding intensifies radial plastic flow, which leads to larger reduction in sample thickness. Larger plastic strain and increased pressure in the central region lead to intensification of PT. However, the effect of the reduction in the friction coefficient on PT kinetics is nonmonotonous. Sliding increases away from the center and with growing rotation and is weakly dependent on the kinetic coefficient. Also, cyclic back and forth torsion is studied and compared to unidirectional torsion. Multiple experimental phenomena, e.g., pressure self-multiplication effect, steps (plateaus) at pressure distribution, flow to the center of a sample, and oscillatory pressure distribution for weaker high-pressure phase, are reproduced and interpreted. Reverse PT in high pressure phase that flowed to the low pressure region is revealed. Possible misinterpretation of experimental PT pressure is found. Obtained results represent essential progress toward understanding of strain-induced PTs under compression and shear in RDAC and may be used for designing experiments for synthesis of new high pressure phases and reduction in PT pressure for known phases, as well as for determination of PT kinetics from experiments

Strain-induced phase transformations under compression, unloading, and reloading in a diamond anvil cell

Strain-induced phase transformations (PTs) in a sample under compression, unloading, and reloading in a diamond anvil cell are investigated in detail, by applying finite element method. In contrast to previous studies, the kinetic equation includes the pressure range in which both direct and reverse PTs occur simultaneously. Results are compared to the case when “no transformation” region in the pressure range exists instead, for various values of the kinetic parameters and ratios of the yield strengths of low and high pressure phases. Under unloading (which has never been studied before), surprising plastic flow and reverse PT are found, which were neglected in experiments and change interpretation of experimental results. They are caused both by heterogeneous stress redistribution and transformation-induced plasticity. After reloading, the reverse PT continues followed by intense direct PT. However, PT is less pronounced than after initial compression and geometry of transformed zone changes. In particular, a localized transformed band of a weaker high pressure phase does not reappear in comparison with the initial compression. A number of experimental phenomena are reproduced and interpreted