Quantum nuclear effects in water using centroid molecular dynamics
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-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.
Molecular modelling is used to calculate transport properties and to study relaxation of liquid n-triacontane (C30H62). The problem is important in connection with the behavior of liquid isolators in a pre-breakdown state. Two all-atom models and a united-atom model are used. Shear viscosity is calculated using the Green–Kubo formula. The force fields are compared with each other using the following criteria: the required time for one molecular dynamics step, the compliance of the main physical and transport properties with experimental values. The problem of the system equilibration is considered. The united-atom potential is used to model the n-triacontane liquid with an initial directional orientation. The time of relaxation to the disordered state, when all molecules orientations are randomized, are obtained. The influence of the molecules orientations on the shear viscosity value and the shear viscosity relaxation are treated.
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.
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 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