Lennard-Jones potentials liquid structure simulation studies

Recently, many experiments have been performed on the structure and dynamics of liquids in porous glasses [175-190]. These studies are difficult to interpret because of the inhomogeneity of the sample. Simulations of water in a cylindrical cavity inside a block of hydrophilic Vycor glass have recently been performed [24,191,192] to facilitate the analysis of experimental results. Water molecules interact with Vycor atoms, using an empirical potential model which consists of (12-6) Lennard-Jones and Coulomb interactions. All atoms in the Vycor block are immobile. For details see Ref. 191. We have simulated samples at room temperature, which are filled with water to between 19 and 96 percent of the maximum possible amount. Because of the hydrophilicity of the glass, water molecules cover the surface already in nearly empty pores no molecules are found in the pore center in this case, although the density distribution is rather wide. When the amount of water increases, the center of the pore fills. Only in the case of 96 percent filling, a continuous aqueous phase without a cavity in the center of the pore is observed. 

Recently, detailed molecular pictures of the interfacial structure on the time and distance scales of the ion-crossing event, as well as of ion transfer dynamics, have been provided by Benjamin s molecular dynamics computer simulations [71, 75, 128, 136]. The system studied [71, 75, 136] included 343 water molecules and 108 1,2-dichloroethane molecules, which were separately equilibrated in two liquid slabs, and then brought into contact to form a box about 4 nm long and of cross-section 2.17 nmx2.17 nm. In a previous study [128], the dynamics of ion transfer were studied in a system including 256 polar and 256 nonpolar diatomic molecules. Solvent-solvent and ion-solvent interactions were described with standard potential functions, comprising coulombic and Lennard-Jones 6-12 pairwise potentials for electrostatic and nonbonded interactions, respectively. While in the first study [128] the intramolecular bond vibration of both polar and nonpolar solvent molecules was modeled as a harmonic oscillator, the next studies [71,75,136] used a more advanced model [137] for water and a four-atom model, with a united atom for each of two... 

---