Mean spherical approximation for the Lennard-Jones-like two Yukawa model: Comparison against Monte Carlo data

Monte Carlo simulation studies are performed for the Lennard-Jones like two Yukawa (LJ2Y) potential to show how properties of this model fluid depend on the replacement of the soft repulsion by the hard-core repulsion. Different distances for the positioning of hard core have been explored. We have found, that for temperatures that are slightly lower and slightly higher of the critical point temperature for the Lennard-Jones fluid, placing the hard core at distances that are shorter than zero-potential energy is well justified by thermodynamic properties that are practically the same as in original LJ2Y model without hard core. However, going to extreme conditions with the high temperature one should be careful since presence of the hard core provokes changes in the properties of the system. The later is extremely important when the mean spherical approximation (MSA) theory is applied to treat the Lennard-Jones-like fluid.Comment: 11 pages, 13 figure

Virial coefficients and vapor-liquid equilibria of the EXP6 and 2-Yukawa fluids

Virial coefficients B₂ through B₄ and the vapor-liquid equilibria for the EXP6 and 2-Yukawa (2Y) fluids have been determined using numerical integrations and Gibbs ensemble simulations, respectively. The chosen 2Y models have been recently determined as an appropriate reference fluid for the considered EXP6 models.Вiрiальнi коефiцiєнти в iд B₂ до B₄ i фазова рiвновага пара - рiдина у EXP6 та 2-Юкава (2Y) плинах розрахованi, вiдповiдно, з допомогою чисельного iнтегрування та на основi комп’ютерного експерименту з використанням ансамблю Гiбса. Вибранi 2Y модельнi системи нещодавно були запропонованi як базиснi для EXP6 плинiв, що розглядаються

Higher order glass-transition singularities in colloidal systems with attractive interactions

The transition from a liquid to a glass in colloidal suspensions of particles interacting through a hard core plus an attractive square-well potential is studied within the mode-coupling-theory framework. When the width of the attractive potential is much shorter than the hard-core diameter, a reentrant behavior of the liquid-glass line, and a glass-glass-transition line are found in the temperature-density plane of the model. For small well-width values, the glass-glass-transition line terminates in a third order bifurcation point, i.e. in a A_3 (cusp) singularity. On increasing the square-well width, the glass-glass line disappears, giving rise to a fourth order A_4 (swallow-tail) singularity at a critical well width. Close to the A_3 and A_4 singularities the decay of the density correlators shows stretching of huge dynamical windows, in particular logarithmic time dependence.Comment: 19 pages, 12 figures, Phys. Rev. E, in prin

Molecular Modeling and Simulation: Force Field Development, Evaporation Processes and Thermophysical Properties of Mixtures

To gain physical insight into the behavior of fluids on a microscopic level as well as to broaden the data base for thermophysical properties especially for mixtures, molecular modeling and simulation is utilized in this work. Various methods and applications are discussed, including a procedure for the development of new force field models. The evaporation of liquid nitrogen into a supercritical hydrogen atmosphere is presented as an example for large scale molecular dynamics simulation. System-size dependence and scaling behavior are discussed in the context of Kirkwood-Buff integration. Further, results for thermophysical mixture properties are presented, i.e. the Henry’s law constant of aqueous systems and diffusion coefficients of a ternary mixture

Percolation line, response functions, and Voronoi polyhedra analysis in supercritical water

The problem of a physical relevance (meaning) of percolation in supercritical fluids is addressed considering a primitive model of water. Two different criteria, physical and configurational, are used for the cluster definition in Monte Carlo simulations over a range of pressures to determine the percolation line and skewness, and a theoretical analytic equation of state is used to evaluate response functions. It is found that both criteria yield practically the same percolation line. However, unlike the findings for simple fluids, the loci of the response function extrema exhibit density/pressure dependence quite different from that of the percolation line. The only potential coincidence between the loci of the extrema of a thermodynamic property and a detectable structural change is found for the coefficient of isothermal compressibility and Voronoi neighbors distribution skewness maximum

Application of molecular simulations: Insight into liquid bridging and jetting phenomena

Molecular dynamics simulations have been performed on pure liquid water, aqueous solutions of sodium chloride, and polymer solutions exposed to a strong external electric field with the goal to gain molecular insight into the structural response to the field. Several simulation methodologies have been used to elucidate the molecular mechanisms of the processes leading to the formation of liquid bridges and jets (in the production of nanofibers). It is shown that in the established nanoscale structures, the molecules form a chain with their dipole moments oriented parallel to the applied field throughout the entire sample volume. The presence of ions may disturb this structure leading to its ultimate disintegration into droplets; the concentration dependence of the threshold field required to stabilize a liquid column has been determined. Conformational changes of the polymer in the jetting process have also been observed