### Article

## Universality of hydrogen bond distributions in liquid and supercritical water

Monte Carlo and molecular dynamics computer simulations using the rigid TIP4P and the flexible BJH intermolecular H2O potentials were carried out for 50 states of supercritical water characterizing a very wide range of thermodynamic conditions, 573 ≤ T ≤ 1273 K; 0.02 ≤ rho ≤ 1.67 g/cm3; 10 ≤ P ≤ 10,000 MPa. Good agreement with available experimental data of the simulated thermodynamic and structural properties give confidence to the quantitative statistical analysis of intermolecular hydrogen bonding under the conditions studied. Energetic, geometric, and angular characteristics of supercritical H-bonds and their distributions at a given temperature remain almost invariant over the entire density range studied from dilute gas-like (~ 0.03 g·cm(-3)) to highly compressed liquid-like (~ 1.5 g·cm(-3)) states. The increase of temperature from ambient to supercritical affects the characteristics of H-bonding in water much more dramatically than the changes in density along any supercritical isotherm. Compared to H-bonds in liquid water under ambient conditions, the H-bonds at 773 K are, on average, 10% weaker, 5% longer, and less linear. Both above and below the H-bonding percolation threshold the fractions of H2O molecules involved in a certain number of H-bonds in liquid and supercritical water closely follow the universal binomial distribution as a function of the average number of H-bonds per one water molecule in the system, as predicted by the independent bond theory. This universal distribution remains intact even when dynamic criteria of H-bonding lifetimes are additionally applied.

The quantum nuclear effects are studied in water using centroid molecular dynamics (CMD) method. The aim is the calibration of CMD implementation in LAMMPS. The calculated intramolecular energy, atoms gyration radii and radial distribution functions are shown in comparison with previous works. The work is assumed to be the step toward to solution of the discrepancy between the simulation results and the experimental data of liquid n-alkane properties in our previous works.

Molecular dynamics study of shear viscosity behavior of liquid aluminum is performed. The embedded atom method potential is used at the simulation of isobaric cooling. The viscosity is calculated using the Green–Kubo formula. The stress autocorrelation functions are obtained in the range 300–1200 K. The calculated kinematic viscosity is in agreement with the experimental data for the temperatures above melting temperature. The steep change of the shear viscosity is found below 650 K which we associate with the glass transition and is in a good agreement with the temperature which is obtained using the calorimetric criterion Kolotova et al (2015 J. Non-Cryst. Solids 429 98). The viscosity coefficient can not be calculated using the direct atomistic simulations below that temperature

It is well recognized that excessive ice accumulation at low-temperature conditions can cause significant damage to civil infrastructure. The passive anti-icing surfaces provide a promising solution to suppress ice nucleation and enhance ice removal. However, despite extensive efforts, it remains a challenge to design anti-icing surfaces with low ice adhesion. Using all-atom molecular dynamics (MD) simulations, we show that surfaces with single-walled carbon nanotube array (CNTA) significantly reduce ice adhesion due to the extremely low solid areal fraction. It was found that the CNTA surface exhibits up to a 45% decrease in the ice adhesion strength in comparison with the atomically smooth graphene surface. The details of the ice detachment from the CNTA surface were examined for different water-carbon interaction energies and temperatures of the ice cube. Remarkably, the results of MD simulations demonstrate that the ice detaching strength depends linearly on the ratio of the ice-surface interaction energy and the ice temperature. These results open the possibility for designing novel robust surfaces with low ice adhesion for passive anti-icing applications.

Molecular-level knowledge of the thermodynamic, structural, and transport properties of water confined by interfaces and nanopores of various materials is crucial for quantitative understanding and prediction of many natural and technological processes, including carbon sequestration, water desalination, nuclear waste storage, cement chemistry, fuel cell technology, etc. Computational molecular modeling is capable to significantly complement the experimental investigations of such systems by providing invaluable atomic-scale information leading to improved understanding of the specific effects of the substrate structure and composition on the structure, dynamics and reactivity of interfacial and nano-confined aqueous solutions. This paper offers a brief overview of recent efforts to quantify some of these effects for individual H2O molecules and hydrated ions confined at the interfaces and in nanopores of several typical hydrophilic and hydrophobic materials. The first molecular layer of aqueous solution at all substrates is often highly ordered, indicating reduced translational and orientational mobility of the H2O molecules. This ordering cannot be simply described as “ice-like”, but rather resembles the behavior of supercooled water or amorphous ice, although with very significant substrate-specific variations.

In this work, we perform coexistence simulations of methane hydrates for pressures up to 5000 bar for different water models. We calculate the kinetic stability boundary of the superheated metastable sI structure and analyze the effects of the heating rate, system size and cage occupancy. We also report molecular dynamics simulation of several possible structuresfor the new hydrogen hydrate clathrate. We show the strength of molecular simulation as a supplement tool for the analysis of experimental data. © 2015 by Nova Science Publishers, Inc. All rights reserved.

Cooling of tokamak boundary plasma owing to radiation of non-fully stripped lithium ions is considered as a promising way for protection of plasma facing elements (PFE) in tokamak. It may be effectively realized when the main part of lithium ions are involved in the closed circuit of migration between plasma and PFE surface. Such an approach may be implemented with the use of lithium device whose hot (500-600 °C) area to be effected by plasma serves as a Li-emitter and the cold part (∼180 °C) as a Li-collector in the shadow. Capillary-pore system (CPS) provides the returning of collected and condensed lithium to emitting zone by capillary forces. The main goals of the last T-11M lithium experiments were investigating Li ions transport in the tokamak scrape of layer (SOL) and their collecting by different kinds of limiters. The design of devices based on lithium CPS with different ratio of emitting/collecting area is the main subject of this paper. © 2015 The Authors.

The dynamics of a two-component Davydov-Scott (DS) soliton with a small mismatch of the initial location or velocity of the high-frequency (HF) component was investigated within the framework of the Zakharov-type system of two coupled equations for the HF and low-frequency (LF) fields. In this system, the HF field is described by the linear Schrödinger equation with the potential generated by the LF component varying in time and space. The LF component in this system is described by the Korteweg-de Vries equation with a term of quadratic influence of the HF field on the LF field. The frequency of the DS soliton`s component oscillation was found analytically using the balance equation. The perturbed DS soliton was shown to be stable. The analytical results were confirmed by numerical simulations.

Radiation conditions are described for various space regions, radiation-induced effects in spacecraft materials and equipment components are considered and information on theoretical, computational, and experimental methods for studying radiation effects are presented. The peculiarities of radiation effects on nanostructures and some problems related to modeling and radiation testing of such structures are considered.

Let G be a semisimple algebraic group whose decomposition into the product of simple components does not contain simple groups of type A, and P⊆G be a parabolic subgroup. Extending the results of Popov [7], we enumerate all triples (G, P, n) such that (a) there exists an open G-orbit on the multiple flag variety G/P × G/P × . . . × G/P (n factors), (b) the number of G-orbits on the multiple flag variety is finite.

I give the explicit formula for the (set-theoretical) system of Resultants of m+1 homogeneous polynomials in n+1 variables