Liquid-vapor coexistence is calculated via molecular dynamics for a variety of parallelepiped shaped molecules. Models are constructed as an array of tangential hard spheres interacting with an attractive square-well potential. Each shape is formed by varying the number of spheres in their three sides.
View Article and Find Full Text PDFIntegral equations of uniform fluids have been considered unable to predict any characteristic feature of the fluid-solid phase transition, including the shoulder that arises in the second peak of the fluid-phase radial distribution function, RDF, of hard-core systems obtained by computer simulations, at fluid densities very close to the structural two-step phase transition. This reasoning is based on the results of traditional integral approximations, like Percus-Yevick, PY, which does not show such a shoulder in hard-core systems, neither in two nor three dimensions. In this work, we present results of three Ansätze, based on the PY theory, that were proposed to remedy the lack of PY analytical solutions in two dimensions.
View Article and Find Full Text PDFA model of a room temperature ionic liquid can be represented as an ion attached to an aliphatic chain mixed with a counter ion. The simple model used in this work is based on a short rigid tangent square well chain with an ion, represented by a hard sphere interacting with a Yukawa potential at the head of the chain, mixed with a counter ion represented as well by a hard sphere interacting with a Yukawa potential of the opposite sign. The length of the chain and the depth of the intermolecular forces are investigated in order to understand which of these factors are responsible for the lowering of the critical temperature.
View Article and Find Full Text PDFMolecular dynamics simulations are performed to clarify the reasons for the disagreement found in a previous publication [G. A. Chapela, F.
View Article and Find Full Text PDFA vibrating version of patchy particles in two dimensions is introduced to study self-assembly of kagome lattices, disordered networks of looping structures, and linear arrays. Discontinuous molecular dynamics simulations in the canonical ensemble are used to characterize the molecular architectures and thermodynamic conditions that result in each of those morphologies, as well as the time evolution of lattice formation. Several versions of the new model are tested and analysed in terms of their ability to produce kagome lattices.
View Article and Find Full Text PDFThe phase behavior of a two-dimensional square-well model of width 1.5σ, with emphasis on the low-temperature and/or high-density region, is studied using Monte Carlo simulation in the canonical and isothermal-isobaric ensembles, and discontinuous molecular-dynamics simulation in the canonical ensemble. Several properties, such as equations of state, Binder cumulant, order parameters, and correlation functions, were computed.
View Article and Find Full Text PDFLiquid-vapor coexistence and interfacial properties of short lineal rigid vibrating chains with three tangent monomers in two and three dimensions are calculated. The effect of the range and position of a long ranged square well attractive potential is studied. Orthobaric densities, vapor pressures, surface tensions, and interfacial widths are reported.
View Article and Find Full Text PDFThe effect of flexibility on liquid-vapor and interfacial properties of tangent linear vibrating square well chains is studied. Surface tension, orthobaric densities, vapor pressures, and interfacial thicknesses are reported and analyzed using corresponding states principles. Discontinuous molecular dynamics simulations in two and three dimensions are performed on rigid tangent linear vibrating square well chains of different lengths.
View Article and Find Full Text PDFThe liquid-vapor phase diagrams of equal size diameter σ binary mixtures of screened potentials have been reported for several ranges of interaction using Monte Carlo simulation methods [J. B. Caballero, A.
View Article and Find Full Text PDFLiquid-vapor coexistence and interfacial properties of square wells in two dimensions are calculated. Orthobaric densities, vapor pressures, surface tensions, and interfacial thicknesses are reported. Results are presented for a series of potential widths λ* = 1.
View Article and Find Full Text PDFA four-site rigid water model is presented, whose parameters are fitted to reproduce the experimental static dielectric constant at 298 K, the maximum density of liquid water and the equation of state at low pressures. The model has a positive charge on each of the three atomic nuclei and a negative charge located at the bisector of the HOH bending angle. This charge distribution allows increasing the molecular dipole moment relative to four-site models with only three charges and improves the liquid dielectric constant at different temperatures.
View Article and Find Full Text PDFLiquid-vapor interfacial properties of square well chains are calculated. Surface tension, orthobaric densities, and vapor pressures are reported. Spinodal decomposition with a discontinuous molecular dynamics simulation program is used to obtain the results which are compared to previously published data for orthobaric densities and vapor pressures.
View Article and Find Full Text PDFMolecular aggregates are formed by heteronuclear vibrating square-well dumbbells. In a recent article [G. A.
View Article and Find Full Text PDFDiscrete perturbation theory (DPT) is a powerful tool to study systems interacting with potentials that are continuous but can be approximated by a piecewise continuous function composed of horizontal segments. The main goal of this work is to analyze the effect of several variables to improve the representation of continuous potentials in order to take advantage of DPT. The main DPT parameters chosen for the purpose are the starting location and size of the horizontal segments used to divide the full range of the potential and its maximum reach.
View Article and Find Full Text PDFSurface tensions and liquid-vapor orthobaric densities are calculated for a wide variety of vibrating square well dumbbells using discontinuous molecular dynamics simulations. The size of the vibration well, the elongation or bond distance of the two particles of the dumbbell, the asymmetry in size (and interaction range) of the two particles, and the depth of the interaction well are the variables whose effects are systematically evaluated in this work. Extensive molecular dynamics simulations were carried out and the orthobaric liquid-vapor densities are compared with those obtained previously by other authors using different methods of simulation for rigid and vibrating square well dumbbells.
View Article and Find Full Text PDFThe liquid-vapor phase equilibria and surface tension of the TIP4P/2005 water model is obtained by using the Ewald summation method to determine the long range Lennard-Jones and electrostatic interactions. The method is implemented in a straightforward manner into standard simulation programs. The computational cost of using Ewald sums in dispersion interactions of water is estimated in direct simulation of interfaces.
View Article and Find Full Text PDFThe crystal formation of NaCl in water is studied by extensive molecular dynamics simulations. Ionic solutions at room temperature and various concentrations are studied using the SPC/E and TIP4P/2005 water models and seven force fields of NaCl. Most force fields of pure NaCl fail to reproduce the experimental density of the crystal, and in solution some favor dissociation at saturated conditions, while others favor crystal formation at low concentration.
View Article and Find Full Text PDFThe Wolf method for the calculation of electrostatic interactions is applied in a liquid phase and at the liquid-vapor interface of water and its results are compared with those from the Ewald sums method. Molecular dynamics simulations are performed to calculate the radial distribution functions at room temperature. The interface simulations are used to obtain the coexisting densities and surface tension along the coexistence curve.
View Article and Find Full Text PDFMolecular dynamics simulations of pure water at the liquid-vapor interface are performed using direct simulation of interfaces in a liquid slab geometry. The effect of intramolecular flexibility on coexisting densities and surface tension is analyzed. The dipole moment profile across the liquid-vapor interface shows different values for the liquid and vapor phases.
View Article and Find Full Text PDF