Vasily V. Pisarev

The nitrite admixtures in cement and concrete are used as corrosion inhibitors for steel reinforcement and also as anti-freezing agents. The characterization of the protective properties should account for the decrease in the concentration of free NO₂⁻ ions in the pores of cement concretes due to their adsorption. Here we applied the classical molecular dynamics computer simulation approach to quantitatively study the molecular scale mechanisms of nitrite adsorption from NaNO₂ aqueous solution on a portlandite surface. We used a new parameterization to model the hydrated NO₂⁻ ions in combination with the recently upgraded ClayFF force field (ClayFF-MOH) for the structure of portlandite. The new NO₂⁻ parameterization makes it possible to reproduce the properties of hydrated NO₂⁻ ions in good agreement with experimental data. In addition, the ClayFF-MOH model improves the description of the portlandite structure by explicitly taking into account the bending of Ca-O-H angles in the crystal and on its surface. The simulations showed that despite the formation of a well-structured water layer on the portlandite (001) crystal surface, NO₂⁻ ions can be strongly adsorbed. The nitrite adsorption is primarily due to the formation of hydrogen bonds between the structural hydroxyls on the portlandite surface and both the nitrogen and oxygen atoms of the NO₂⁻ ions. Due to that, the ions do not form surface adsorption complexes with a single well-defined structure but can assume various local coordinations. However, in all cases, the adsorbed ions did not show significant surface diffusional mobility. Moreover, we demonstrated that the nitrite ions can be adsorbed both near the previously-adsorbed hydrated Na⁺ ions as surface ion pairs, but also separately from the cations.

The behavior of brucite over wide ranges of temperatures and pressures is of great interest for fundamental geochemistry and geophysics. Brucite layers and their octahedral Mg(OH)₆ structural units constitute an important structural part of layered dense magnesium hydrous silicates (DMHS), which play a major role in mineral equilibria controlling water balance in the subduction zones of the upper mantle. The ClayFF force field was originally developed for atomistic computer simulations of clays and other layered minerals and their hydrated interfaces. The crystallographic parameters of brucite at 25 °C and 1 bar were used, among several others, to develop the original ClayFF parametrization. Its new recent modification, ClayFF-MOH, can more accurately account for the bending of Mg–O–H angles in the brucite structure, and it was used here to test the applicability of this simple classical model over very wide ranges of temperature and pressure well beyond the range of its original implementation (up to 600 °C and 15 GPa). The pressure and temperature dependencies of brucite crystallographic parameters, the compressibility of the crystal lattice, the coefficients of thermal expansion, and the vibrational spectra were calculated in a series of classical molecular dynamics simulations using the ClayFF-MOH model and compared with a diverse set of available experimental data, including X-ray diffractometry, neutron scattering, IR and Raman spectroscopy. These new results demonstrated that ClayFF-MOH, as simple and approximate as it is, can be quite accurate in predicting many mineral properties at subduction zone conditions, which greatly expands the area of its applicability.

The use of a new modification of the ClayFF force field for molecular dynamics modeling of ettringite crystals and interaction of aqueous NaCl and Na₂SO₄ solutions with their surface shows that the possibility of explicitly taking into account metal–O–H interactions in the system leads to the formation of stronger hydrogen bonds in the crystal structure and on the surface and to greater localization of the atoms of both the crystalline phase and solutions in the subsurface zone. The relative fractions of the inner- and outer-spheric adsorption of Na⁺, Cl^–, and (SO₄)^2– ions and the mobility of H₂O molecules in both solutions also change. The calculated crystal lattice parameters and the density of ettringite remained almost unchanged between the old and new versions of the ClayFF force field, but the accuracy of reproduction of the elastic characteristics of the crystal markedly increased.

In this paper, we simulate the nucleation and growth of crystalline nuclei in a molybdenum film cooled at different rates confined between two amorphous walls. We also compare the results for the wall-confined and wall-free systems. We apply the same methodology as in the work (Kirova and Pisarev 2019 J. Cryst. Growth 528 125266) which is based on reconstructing the probability density function for the largest crystalline nucleus in the system. The size of the nucleus and the asphericity parameter are considered as the reaction coordinates. We demonstrate that in both the free and confined systems there are two mechanisms of crystal growth: the attachment of atoms to the biggest crystal from the amorphous phase and the merging of the biggest crystal cluster with small ones (coalescence). We show that the attachment mechanism is dominant in the melt cooled down at a slower rate, and the mechanism gradually shifts to coalescence as cooling rate increases. We also observe the formation of long-lived crystal clusters and demonstrate that amorphous walls do not affect their geometric characteristics. However, system confined between walls demonstrates higher glass-forming ability.

Computational predictions of the high-pressure viscosity of hydrocarbon mixtures could help to accelerate the development of fuels and lubricants with improved performance. In this study, we use molecular dynamics simulations to study the viscosity and density of methylcyclohexane, 1-methylnaphthalene, and their binary mixtures at 323 K and pressures of up to 500 MPa. The simulation results are in excellent agreement with previous experiments available up to 100 MPa for both pure compounds (200 MPa for 1-methylnaphthalene) and the binary mixtures. For 1-methylnaphthalene, the viscosity initially increases slower-than-exponential with pressure before it reaches an inflection point and then increases faster-than-exponential. The inflection point (300 MPa) occurs at a pressure slightly below the one at which 1-methylnaphthalene is expected to enter the supercooled phase (400 MPa). For methylcyclohexane, the increase in viscosity with pressure is slower-than-exponential over the entire pressure range studied. The binary mixtures show intermediate pressure–viscosity responses between the two pure cases. The applicability of equations commonly used to describe the pressure dependence of viscosity, as well as the viscosity of binary mixtures, is evaluated against the computational predictions.

In this paper, the predictive power of molecular dynamics methods is demonstrated for the cases of model paraffinic and aromatic lubricant liquids at pressures up to 400 MPa. The shear viscosity and self-diffusion coefficients are calculated for 2,2,4-trimethylpentane (C₈H₁₈) at 298 K and 1,1-diphenylethane (C₁₄H₁₄) at 333 K. Three force fields with different levels of accuracy are compared by the ability to predict the experimental data. The Stokes–Einstein correlation between viscosity and self-diffusion is demonstrated for both compounds.