BZPEERReading, trust and collaboration for researchDiscover
Workspaces
Network
Opportunities
Account
Researcher profile
Yu. D. Fomin contributes to research discovery and scholarly infrastructure.
Trust snapshot
Actions
Identity and collaboration
Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.
Log in to claimDirect collaboration
Claim this author entity first to unlock direct invitations.
Research graph
Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.
BZPEER is loading the nearby papers, people, topics and institutions for this page.
Published work
We report an ab-initio simulation of liquid Indium in a wide range of pressures and temperatures. We calculate equation of state, thermal expansion and compressibility coefficients. The structure of the system is analyzed by radial distribution functions and structure factors. The results are compared with available experimental data.
We investigate a behavior of a core-softened system in a thin slit pore (the width of the pore is equil to three diameters of the particles). In previous studies it was shown that strongly confined systems form crystalline phases which consist of several triangular or square layers. These phases can be also considered as cuts of FCC or HCP structures. We show that the behavior of the core-softened system is more complex. We observe also a quasicrystalline phase. Moreover, the phase with two triangular layers appears at lower densities than the one with two square layers which is in contrast to the behavior of the systems studied before. These results demonstrate that the phase behavior of strongly confined systems can be even more complex than it was supposed before.
In the present paper we study the thermodynamic properties of superionic conductor $SrCl_2$ at high temperatures by means of molecular dynamics method. Firstly, we calculate the melting line. Then we compute the equations of state and the response functions (heat capacity, thermal expansion coefficient, etc) at the temperatures up to the melting. We show that the response functions show maxima or minima at the temperatures well above the temperature of transition into the conductive state, and therefore are not related to this transition.
Phase diagram of a two-dimensional system with a potential which stabilizes Kagome lattice is calculated. It is shown that this system demonstrate a set of crystalline and the regions of stability of these phases are calculated. The scenarios of melting of triangular and square crystals of the system are determined.
Investigation of the longitudinal and transverse excitations in liquids is of great importance for understanding the fundamentals of the liquid state of matter. One of the important questions is the temperature and density dependence of the frequency of the excitations. In our recent works it was shown that while in simple liquids the frequency of longitudinal excitations increases when the temperature is increased isochorically, in water the frequency can anomalously decrease with the temperature increase. In the present manuscript we study the dispersion curves of longitudinal and transverse excitations of water and liquid silicon modelled by Stillinger-Weber potential. We show that both substances demonstrate the anomaly of the dispersion curves, but it the case of water it is more pronounced.
The random disorder can drastically change the melting scenario of two-dimensional systems and has to be taken into account in the interpretation of the experimental results. We present the results of the molecular dynamics simulations of the two dimensional system with the core-softened potentials with two repulsive steps in which a small fraction of the particles is pinned, inducing quenched disorder. It is shown that without the quenched disorder the system with small repulsive shoulder, which is close in the shape to the soft disks, melts in accordance with the melting scenario proposed in Refs. [30,31,33] (first-order liquid-hexatic and continuous hexatic-solid transitions). Random pinning widens the hexatic phase, but leaves the melting scenario unchanged. For the system with larger repulsive step at high densities the conventional first-order transition takes place without random pinning. However, in the presence of disorder the single first-order transition transforms into two transitions, one of them (solid-hexatic) is the continuous KTHNY-like transition, while the hexatic to isotropic liquid transition occurs as the first order transition.
In experiments the two-dimensional systems are realized mainly on solid substrates which introduce quenched disorder due to some inherent defects. The defects of substrates influence the melting scenario of the systems and have to be taken into account in the interpretation of the experimental results. We present the results of the molecular dynamics simulations of the two dimensional system with the core-softened potential in which a small fraction of the particles is pinned, inducing quenched disorder.The potentials of this type are widely used for the qualitative description of the systems with the water-like anomalies. In our previous publications it was shown that the system demonstrates an anomalous melting scenario: at low densities the system melts through two continuous transition in accordance with the Kosterlitz-Thouless-Halperin-Nelson-Young (KTHNY) theory with the intermediate hexatic phase, while at high densities the conventional first order melting transition takes place. We find that the well-known disorder-induced widening of the hexatic phase occurs at low densities, while at high density part of the phase diagram random pinning transforms the first-order melting into two transitions: the continuous KTHNY-like solid-hexatic transition and first-order hexatic-isotropic liquid transition.
Supercritical state has been viewed as an intermediate state between gases and liquids with largely unknown physical properties. Here, we address the important ability of supercritical fluids to sustain collective excitations. We directly study propagating modes on the basis of correlation functions calculated from extensive molecular dynamics simulations, and find that the supercritical system sustains propagating solid-like transverse modes below the Frenkel line but becomes devoid of transverse modes above the line where it supports longitudinal modes only. Important thermodynamic implications of this finding are discussed. We directly detect positive sound dispersion (PSD) below the Frenkel line where transverse modes are operative, quantitatively explain its magnitude on the basis of transverse and longitudinal velocities. PSD disappears above the Frenkel line which therefore demarcates the supercritical phase diagram into two areas where PSD does and does not operate.
It is well known that some liquids can demonstrate anomalous behavior. Interestingly, this be- havior can be qualitatively reproduced with simple core-softened isotropic pair-potential systems. Although anomalous properties of liquids usually take place at low and moderate temperatures it was recently recognized that many important phenomena can appear in supercritical fuids. However, no studies of supercritical behavior of core-softened fuids is reported. This paper reports a study of dynamical crossover in supercritical core-softened systems. The crossover line is calculated from three different criteria and good agreement between them is observed. It is found that the behavior of the dynamical crossover line of core-softened systems is quite complex due to its quasi-binary nature.
Dynamical crossover in water is studied by means of computer simulation. The crossover temperature is calculated from the behavior of velocity autocorrelation functions. The results are compared with experimental data. It is shown that the qualitative behavior of the dynamical crossover line is similar to the melting curve behavior. Importantly, the crossover line belongs to experimentally achievable $(P,T)$ region which stimulates the experimental investigation in this field.
It is well known that confining a liquid into a pore strongly alters the liquid behavior. Investigations of the effect of confinement are of great importance for many scientific and technological applications. Here we present a molecular dynamics study of the behavior of cyclohexane confined in carbon slit pores. The local structure and orientational ordering of cyclohexane molecules are investigated. It is shown that the system freezes with decreasing the pore width, and the freezing temperature of nanoconfined cyclohenae is higher than the bulk one.
We present the molecular simulation study of the behavior of water and sodium chloride solution confined in lizardite asbestos nanotube which is a typical example of hydrophilic confinement. The local structure, orientational and dynamic properties are studied. It is shown that the diffusion coefficient drops about two orders of magnitude comparing to the bulk case, and water in lizardite asbestos tubes experiences vitrification rather then crystallization upon cooling in accordance with the results for some other hydrophilic confinements. The behavior of sodium chloride solutions also considered and the formation of double layer is observed. It is shower that both sodium and chlorine have larger diffusion coefficients then water.
We consider the phase behavior of two-dimensional ($2D$)system of particles with an isotropic core-softened potential introduced in our previous publications. As one can expect from the qualitative consideration for the three dimensional case, the system demonstrates a reentrant-melting transition at low densities along with waterlike anomalies in the fluid phase near the melting maximum. In contrast with the three dimensional case, in two dimensions melting is a continuous two-stage transition in the low density part of the phase diagram with an intermediate hexatic phase corresponding to the Kosterlitz-Thouless-Halperin-Nelson-Young (KTHNY) scenario. At the same time, at high densities the system melts through one first-order transition. We also show, that the order of the region of anomalous diffusion and the regions of density and structural anomalies are inverted in comparison with the $3D$ case and have silicalike sequence.
Recent advance in understanding the supercritical state posits the existence of a new line above the critical point separating two physically distinct states of matter: rigid liquid and non-rigid gas-like fluid. The location of this line, the Frenkel line, remains unknown for important real systems. Here, we map the Frenkel line on the phase diagram of supercritical iron using molecular dynamics simulations. On the basis of our data, we propose a general recipe to locate the Frenkel line for any system, the recipe that importantly does not involve system-specific detailed calculations and relies on the knowledge of the melting line only. We further discuss the relationship between the Frenkel line and the metal-insulator transition in supercritical liquid metals. Our results enable predicting the state of supercritical iron in several conditions of interest. In particular, we predict that liquid iron in the Jupiter core is in the "rigid liquid" state and is highly conducting. We finally analyse the evolution of iron conductivity in the core of smaller planets such as Earth and Venus as well as exoplanets: as planets cool off, the supercritical core undergoes the transition to the rigid-liquid conducting state at the Frenkel line.
In the present paper, using a molecular dynamics simulation, we study a nature of melting of a two-dimensional ($2D$) system of classical particles interacting through a purely repulsive isotropic core-softened potential which is used for the qualitative description of the anomalous behavior of water and some other liquids. We show that the melting scenario drastically depends on the potential softness and changes with increasing the width of the smooth repulsive shoulder. While at small width of the repulsive shoulder the melting transition exhibits what appears to be weakly first-order behavior, at larger values of the width a reentrant-melting transition occurs upon compression for not too high pressures, and in the low density part of the $2D$ phase diagram melting is a continuous two-stage transition, with an intermediate hexatic phase in accordance with the Kosterlitz-Thouless-Halperin-Nelson-Young (KTHNY) scenario. On the other hand, at high density part of the phase diagram one first-order transition takes place. These results may be useful for the qualitative understanding the behavior of water confined between two hydrophobic plates.
Some fluids are known to behave anomalously. The so-called structural anomaly which means that the fluid becomes less structures under isothermal compression is among the most frequently discussed ones. Several methods for quantifying the degree of structural order are described in the literature and are used for calculating the regions of structural anomalies. It is implied that all of the structural order determinations yield qualitatively identical results. However, no explicit comparison was made. This paper presents such a comparison for the first time. the results of some definitions are shown to contradict the intuitive notion of a fluid. On the basis of this comparison we show that the structural anomaly can be most reliably determined from the behavior of the excess entropy.
In the framework of the density functional theory of freezing proposed in our previous works, we calculate the phase diagram of two-dimensional system of particles interacting through the repulsive shoulder potential. This potential consists of the hard core and repulsive shoulder of the larger radius. It is shown that at low densities the system melts through the continuous transition in accordance with the Kosterlitz-Thouless-Halperin-Nelson-Young (KTHNY) scenario, while at high densities the conventional first order transition takes place.
In the present paper we propose the van der Waals-like model, which allows a purely analytical study of fluid properties including the equation of state, phase behavior and supercritical fluctuations. We take a square-well system as an example and calculate its liquid - gas transition line and supercritical fluctuations. Employing this model allows us to calculate not only the thermodynamic response functions (isothermal compressibility $β_T$, isobaric heat capacity $C_P$, density fluctuations $ζ_T$, and thermal expansion coefficient $α_T$), but also the correlation length in the fluid $ξ$. It is shown that the bunch of extrema widens rapidly upon departure from the critical point. It seems that the Widom line defined in this way cannot be considered as a real boundary that divides the supercritical region into the gaslike and liquidlike regions. As it has been shown recently, the new dynamic line on the phase diagram in the supercritical region, namely the Frenkel line, can be used for this purpose.
We present the molecular dynamics study of benzene molecules confined into the single wall carbon nanotube. The local structure and orientational ordering of benzene molecules are investigated. It is found that the molecules mostly group in the middle distance from the axe of the tube to the wall. The molecules located in the vicinity of the wall demonstrate some deviation from planar shape. There is a tilted orientational ordering of the molecules which depends on the location of the molecule. It is shown that the diffusion coefficient of the benzene molecules is very small at the conditions we report here.
Behavior of supercritical fluids attracts a lot of attention nowadays. It is important both from the point of view of fundamental science and technological applications. However, up to now the progress in the field is rather moderate. In this article we report a computational study of supercritical carbon dioxide which is one of the most important fluids for chemical industry. We study the response functions of $CO_2$ in supercritical regime and calculate the locations of their maxima (Widom lines). We also make preliminary calculations of the line of crossover of microscopic dynamics of particles (Frenkel line). The conclusions on the Frenkel line location can be applied to study of the atmosphere of Venus.
Recently, we have proposed a new dynamic line on the phase diagram in the supercritical region. Crossing this line corresponds to the radical changes of the fluid properties. Here, we focus on the dynamics of model Lennard-Jones and Soft-Sphere fluids. We show that the change of the dynamics from the liquid-like to gas-like can be established on the basis of the velocity autocorrelation function and mean-square displacement. Using the rigorous criterion, we show that the crossover of particle dynamics and key liquid properties occurs at the same line. We further show that positive sound dispersion disappears in the vicinity of this line in both kinds of systems. The dynamic line bears no relationship to the existence of the critical point. We find that the region of existence of liquid-like dynamics narrows with the increase of the exponent of the repulsive part of inter-particle potential.
The state of matter above the critical point is terra incognita, and is loosely discussed as a physically homogeneous flowing state where no differences can be made between a liquid and a gas and where properties undergo no marked or distinct changes with pressure and temperature. In particular, the structure of supercritical state is currently viewed to be the same everywhere on the phase diagram, and to change only gradually and in a featureless way while moving along any temperature and pressure path above the critical point. Here, we demonstrate that this is not the case, but that there is a well-defined structural crossover instead. Evidenced by the qualitative changes of distribution functions of interatomic distances and angles, the crossover demarcates liquid-like and gas-like configurations and the presence of medium-range structural correlations. Importantly, the discovered structural crossover is closely related to both dynamic and thermodynamic crossovers operating in the supercritical state, providing new unexpected fundamental interlinks between the supercritical structure, dynamics and thermodynamics.
We report a molecular dynamics study of transport coefficients and infinite frequency shear mod- ulus of liquid iron at high temperatures and high pressures. We observe a simultaneous rise of both shear viscosity and diffusion coefficient along the melting line and estimate if liquid iron can vitrify under Earth-core conditions. We show that in frames of the model studied in our work iron demonstrates a moderate increase of viscosity along the melting line. It is also demonstrated that in the limit of high temperatures and high pressures the liquid iron behaves similar to the soft spheres system with exponent n=4.6.
In this paper we present a simulation study of density, structural and diffusion anomalies in core-softened system introduced in our previous publications. It is well-known, that with appropriate parametrization, core-softened systems are remarkable model liquids that exhibit anomalous properties observed in tetrahedral liquids such as silica and water. It is widely believed that core-softened potentials demonstrate the water-like sequence of anomalies. We show that with increasing the depth of the attractive part of the potential the order of the region of anomalous diffusion and the regions of density and structural anomalies is inverted and have the silica-like sequence. We also show that the slope of the Widom line is negative like in water.
The present article presents a molecular dynamics study of several anomalies of core-softened systems. It is well known that many core-softened liquids demonstrate diffusion anomaly. Usual intuition relates the diffusion coefficient to shear viscosity via Stockes-Einstein relation. However, it can break down at low temperature. In this respect it is important to see if viscosity also demonstrates anomalous behavior.
In the present article we carry out a molecular dynamics study of the core-softened system and show that the existence of the water-like anomalies in this system depends on the trajectory in $P-ρ-T$ space along which the behavior of the system is studied. For example, diffusion and structural anomalies are visible along isotherms as a function of density, but disappears along the isochores and isobars as a function of temperature. On the other hand, the diffusion anomaly may be seen along adiabats as a function of temperature, density and pressure. It should be noted that it may be no signature of a particular anomaly along a particular trajectory, but the anomalous region for that particular anomaly can be defined when all possible trajectories in the same space are examined (for example, signature of diffusion anomaly is evident through the crossing of different isochors. However, there is no signature of diffusion anomaly along a particular isochor). We also analyze the applicability of the Rosenfeld entropy scaling relations to this system in the regions with the water-like anomalies. It is shown that the validity of the Rosenfeld scaling relation for the diffusion coefficient also depends on the trajectory in the $P-ρ-T$ space along which the kinetic coefficients and the excess entropy are calculated.
In this paper we present a simulation study of water-like anomalies in core-softened system introduced in our previous publications. We investigate the anomalous regions for a system with the same functional form of the potential but with different parameters and show that the order of the region of anomalous diffusion and the region of density anomaly is inverted with increasing the width of the repulsive shoulder.
It is generally agreed that the supercritical region of a liquid consists of one single state (supercritical fluid). On the other hand, we show here that liquids in this region exist in two qualitatively different states: "rigid" and "non-rigid" liquid. Rigid to non-rigid transition corresponds to the condition τ ~ τ0, where τis liquid relaxation time and τ0 is the minimal period of transverse quasi-harmonic waves. This condition defines a new dynamic line on the phase diagram, and corresponds to the loss of shear stiffness of a liquid at all available frequencies, and consequently to the qualitative change of many important liquid properties. We analyze the dynamic line theoretically as well as in real and model liquids, and show that the transition corresponds to the disappearance of high-frequency sound, qualitative changes of diffusion and viscous flow, increase of particle thermal speed to half of the speed of sound and reduction of the constant volume specific heat to 2kB per particle. In contrast to the Widom line that exists near the critical point only, the new dynamic line is universal: it separates two liquid states at arbitrarily high pressure and temperature, and exists in systems where liquid - gas transition and the critical point are absent overall.
We report a detailed simulation study of the phase behavior of core softened system with attractive well. Different repulsive shoulder widthes and attractive well depthes are considered which allows to monitor the influence of repulsive and attractive forces on the phase diagram of the system. Thermodynamic anomalies in the systems are also studied. It is shown that the diffusion anomaly is stabilized by small attraction.
Molecular dynamics computer simulation has been used to compute the self-diffusion coefficient, and shear viscosity of soft-sphere fluids, in which the particles interact through the soft-sphere or inverse power pair potential. The calculations have been made along the melting line in a wide range of pressures and temperatures. The validity of scaling relations for thermodynamic parameters and kinetic coefficients was checked. It was shown that the Stokes-Einstein relationship is obeyed if the Barker diameter is used as a characteristic length scale. It was also shown that the viscosity is non-monotonic along the isochores as predicted by Ya. Rosenfeld. It was shown that the viscosity is strongly growing along the melting line, however, this increase does not stimulate the glass transition because the relaxation time is decreasing.
We show that the existence of the water-like anomalies in kinetic coefficients in the core-softened systems depends on the trajectory in $ρ-T$ plane along which the kinetic coefficients are calculated. In particular, it is shown that the diffusion anomaly does exist along the isotherms, but disappears along the isochores. We analyze the applicability of the Rosenfeld entropy scaling relations to the systems with the core-softened potentials demonstrating the water-like anomalies. It is shown that the validity of the of Rosenfeld scaling relation for the diffusion coefficient also depends on the trajectory in the $ρ-T$ plane along which the kinetic coefficients and the excess entropy are calculated. In particular, it is valid along isochors, but it breaks down along isotherms.
In the present article we carry out a molecular dynamics study of the core-softened system and show that the existence of the water-like anomalies in this system depends on the trajectory in $P-ρ-T$ space along which the behavior of the system is studied. For example, diffusion and structural anomalies are visible along isotherms, but disappears along the isochores and isobars, while density anomaly exists along isochors. We analyze the applicability of the Rosenfeld entropy scaling relations to this system in the regions with the water-like anomalies. It is shown that the validity of the of Rosenfeld scaling relation for the diffusion coefficient also depends on the trajectory in the $P-ρ-T$ space along which the kinetic coefficients and the excess entropy are calculated.
This articles presents a simulation study of the applicability of the Rosenfeld entropy scaling to the systems which can not be approximated by effective hard spheres. Three systems are studied: Herzian spheres, Gauss Core Model and soft repulsive shoulder potential. These systems demonstrate the diffusion anomalies at low temperatures: the diffusion increases with increasing density or pressure. It is shown that for the first two systems which belong to the class of bounded potentials the Rosenfeld scaling formula is valid only in the infinite temperature limit where there are no anomalies. For the soft repulsive shoulder the scaling formula is valid already at sufficiently low temperatures, however, out of the anomaly range.
We report a computer-simulation study of the equilibrium phase diagram of a three-dimensional system of particles with a repulsive step potential. The phase diagram is obtained using free-energy calculations. At low temperatures, we observe a number of distinct crystal phases. We show that at certain values of the potential parameters the system exhibits the water-like thermodynamic anomalies: density anomaly and diffusion anomaly. The anomalies disappear with increasing the repulsive step width: their locations move to the region inside the crystalline phase.