Aqueous divalent metal-nitrate interactions: hydration versus ion pairing

Nitrate aqueous solutions, Mg(NO(3))(2), Ca(NO(3))(2), Sr(NO(3))(2), and Pb(NO(3))(2), are investigated using Raman spectroscopy and free energy profiles from molecular dynamics (MD) simulations. Analysis of the in-plane deformation, symmetric stretch, and asymmetric stretch vibrational modes of the nitrate ions reveal perturbation caused by the metal cations and hydrating water molecules. Results show that Pb(2+) has a strong tendency to form contact ion pairs with nitrate relative to Sr(2+), Ca(2+), and Mg(2+), and contact ion pair formation decreases with decreasing cation size and increasing cation charge density: Pb(2+) > Sr(2+) > Ca(2+) > Mg(2+). In the case of Mg(2+), the Mg(2+)-OH(2) intermolecular modes indicate strong hydration by water molecules and no contact ion pairing with nitrate. Free energy profiles provide evidence for the experimentally observed trend and clarification between solvent-separated, solvent-shared, and contact ion pairs, particularly for Mg(2+) relative to other cations.     

Coarse-grain modelling using an equation-of-state many-body potential: application to fluid mixtures at high temperature and high pressure

ABSTRACT

A many-body, coarse-grain model, termed the product gas mixture model, is presented that accurately describes the thermodynamic behaviour of molecular mixtures. The coarse-grain model is developed by first approximating the mixture as a van der Waals one-fluid, and subsequently applying an exponential-6 equation-of-state to describe the forces and energies between particles in the spirit of the many-body model pioneered by Pagonabarraga and Frenkel. Isothermal dissipative particle dynamics simulations are carried out at thermochemical states that occur during decomposition of a prototypical energetic material, RDX (1,3,5-trinitro-1,3,5-triazinane). The product gas mixture model performance is assessed by comparing to an explicit-molecule model and a hard-core coarse-grain model based on the exponential-6 pair potential. Overall, the many-body, coarse-grain model is shown to accurately capture the structural and thermodynamic properties for the wide variety of thermochemical states considered, while the hard-core coarse-grain model cannot. The many-body, coarse-grain model overcomes the issues of transferability, scaling consistency and unphysical ordered phase behaviour that often afflict coarse-grain models. While specific thermochemical conditions related to RDX decomposition are considered, the results are generally applicable to the thermodynamic behaviour of other fluid mixtures at both moderate and extreme conditions.     

Study of structural and electronic transport properties of Ce-doped LaMnO3

The structural and electronic transport properties of La_1−x Ce _x MnO₃ (x=0.0–1.0) have been studied. All the samples exhibit orthorhombic crystal symmetry and the unit cell volume decreases with Ce doping. They also make a metal-insulator transition (MIT) and transition temperature increases with increase in Ce concentration up to 50% doping. The system La_0.5Ce_0.5MnO₃ also exhibits MIT instead of charge-ordered state as observed in the hole doped systems of the same composition.     

Highly scalable discrete-particle simulations with novel coarse-graining: accessing the microscale

ABSTRACT

Simulating energetic materials with complex microstructure is a grand challenge, where until recently, an inherent gap in computational capabilities had existed in modelling grain-scale effects at the microscale. We have enabled a critical capability in modelling the multiscale nature of the energy release and propagation mechanisms in advanced energetic materials by implementing, in the widely used LAMMPS molecular dynamics (MD) package, several novel coarse-graining techniques that also treat chemical reactivity. Our innovative algorithmic developments rooted within the dissipative particle dynamics framework, along with performance optimisations and application of acceleration technologies, have enabled extensions in both the length and time scales far beyond those ever realised by atomistic reactive MD simulations. In this paper, we demonstrate these advances by modelling a shockwave propagating through a microstructured material and comparing performance with the state-of-the-art in atomistic reactive MD techniques. As a result of this work, unparalleled explorations in energetic materials research are now possible.     

A molecular dynamics study of alkaline earth metal-chloride complexation in aqueous solution

The relative stability of alkaline earth metals (M2+ = Mg2+, Ca2+, Sr2+, and Ba2+) and their chloride complexes in aqueous solution is examined through molecular dynamics simulations using a flexible SPC water model with an internally consistent set of metal ion force field parameters. For each metal-chloride ion pair in aqueous solution, the free energy profile was calculated via potential of mean force simulations. The simulations provide detailed thermodynamic information regarding the relative stability of the different types of metal-chloride pairs. The free energy profiles indicate that the preference for contact ion pair formation increases with ionic radius and is closely related to the metal hydration free energies. The water residence times within the first hydration shells are in agreement with residence times reported in other computational studies. Calculated association constants suggest an increase in metal-chloride complexation with increasing cation radii that is inconsistent with experimentally observed trends. Possible explanations for this discrepancy are discussed.     

Two connected inverse spectral transforms

We establish a transformation which connects the potentials of the one-dimensional Dirac and Klein-Gordon operators. This transformation links the solutions of the nonlinear evolution equations solvable by means of the two inverse spectral transforms which use the Dirac and Klein-Gordon direct and inverse spectral problems.