Lennard-Jones fluid models simulations

An example drawn from Deitrick s work (Fig. 2) shows the chemical potential and the pressure of a Lennard-Jones fluid computed from molecular dynamics. The variance about the computed mean values is indicated in the figure by the small dots in the circles, which serve only to locate the dots. A test of the thermodynamic goodness of the molecular dynamics result is to compute the chemical potential from the simulated pressure by integrating the Gibbs-Duhem equation. The results of the test are also shown in Fig. 2. The point of the example is that accurate and affordable molecular simulations of thermodynamic, dynamic, and transport behavior of dense fluids can now be done. Currently, one can simulate realistic water, electrolytic solutions, and small polyatomic molecular fluids. Even some of the properties of micellar solutions and liquid crystals can be captured by idealized models [4, 5]. 

A commonly used model system in liquid crystal simulation is the Gay-Beme fluid. It can be regarded as a Lennard-Jones fluid generalised to ellipsoidal molecular cores. 

Wongkoblap et al.307 study Lennard-Jones fluids in finite pores, and compare their results with Grand canonical ensemble simulations of infinite pores. Slit pores of 3 finite layers of hexagonally arranged carbon atoms were constructed. They compare the efficiency of Gibbs ensemble simulations (where only the pore is modelled) with Canonical ensemble simulations where the pore is situated in a cubic cell with the bulk fluid, and find that while the results are mostly the same, the Gibbs ensemble method is more efficient. However, the meniscus is only able to be modelled in the canonical ensemble. 

A comparison of mesoscopic simulation methods with MD simulations has been performed by Denniston and Robbins.423 They study a binary mixture of simple Lennard-Jones fluids and map out the required parameters of the mesoscopic model from their MD simulation data. Their mapping scheme is more complete than those of previous workers because in addition to accounting for the interfacial order parameter and density profiles, they also consider the stress. Their mapping consists of using MD simulations to parameterise the popular mesoscale Lattice Boltzmann simulation technique and find that a... 

The adsorption of gas mixtures has been extensively studied. For example, Wendland et al. [64] applied the Bom—Green—Yvon approach using a coarse grained density to study the adsorption of subcritical Lennard-Jones fluids. In a subsequent paper, they tested their equations with simulated adsorption isotherms of several model mixtures [65]. They compared the adsorption of model gases with an equal molecular size but different adsorption potentials. They discussed the stmcture of the adsorbed phase, adsorption isotherms, and selectivity curves. Based on the vacancy solution theory [66], Nguyen and Do [67] developed a new technique for predicting the multicomponent adsorption equihbria of supercritical fluids in microporous carbons. They concluded that the degree of adsorption enhancement, due to the proximity of the pore... 

Subsequent applications of semigrand methods have been numerous, as species-identity changes have become a standard practice when simulating mixtures. We would fail in an attempt to mention all such uses, so instead we will sample some of the more interesting applications and extensions. Hautman and Klein [22] examined, by molecular dynamics a breathing Lennard-Jones fluid of fluctuating particle diameter the breathing modes are introduced to better model molecules that are treated as LJ atoms. Liu... 

In order to ascertain whether the 3-regime behavior observed in the experimental vibrational lifetimes is indeed a result of local density enhancements, Goodyear and Tucker [12] computed both vibrational lifetimes and local density enhancements from molecular dynamics simulation for a model solute-solvent SCF solution. These authors considered a diatomic solute in a 2-dimensional supercritical Lennard-Jones fluid of 1150 atoms (Fig. 1). In this model, each of the solute atoms was designated as a Lennard-Jones site, and the Lennard-Jones parameters between solute and solvent atoms were taken to be the same as those between solvent atoms. The vibrational lifetimes were computed using the standard, classical Landau-Teller expression [69,70,72,73,78], i.e. 

This section is devoted to studying the 2D Lennard-Jones model in order to serve as the basis in applying Steele s theory. In Section IVA the main studies about that model are summarized and commented on. In Section IVB, the most useful expressions for the equation of state of the model are given. In Section IVC we present results about the application of these equations, which are compared with other theoretical approaches to studying adsorption of 2D Lennard-Jones fluids onto perfectly flat surfaces. In Section FVD, the comparison with experimental results is made, including results for the adsorption isotherms, the spreading pressure, and the isosteric heat. Finally, in Section IVE we indicate briefly some details about the use of computer simulations to model the properties both of an isolated 2D Lennard-Jones system and of adsorbate-adsorbent systems. 

Nagayama et al. [57] carried out nonequilibrium molecular dynamic simulations to study the effect of interface wettability on the pressure driven flow of a Lennard-Jones fluid in a nanochannel. The velocity profile changed significantly depending on the wettability of the wall. The no-slip boundary condition breaks down for a hydrophobic wall. Siegel et al. [58] developed a two-dimensional computational model for fuel cells. 

In mixed MD/MC simulations, some of the atoms are moved by the MD method and some of the atoms are moved by the MC method. LaBerge et al. [48] demonstrated that this method rigorously converges to the same equilibrium state as either MC or canonical MD alone. Thus, it was shown that the intermption of the forces produced by the application of the MC moves does not incorrectly bias the evolution of the MD particles. This technique was applied by the above authors to a Lennard-Jones fluid. It was anticipated that this model would be superior to either MD or MC on its own, in systems where some particles are more efficiently sampled by MD (for instance solvent motions), while others are more efficiently sampled by MC (for instance highly correlated motions). 

To illustrate the relationship between the microscopic structure and experimentally accessible information, we compute pseudo-experimental solvation-force curves F h)/R [see Eq. (22)] as they would be determined in SEA experiments from computer-simulation data for T z [see Eqs. (93), (94), (97)]. Numerical values indicated by an asterisk are given in the customary dimensionless (i.e., reduced) units (see [33,75,78] for definitions in various model systems). Results are correlated with the microscopic structure of a thin film confined between plane parallel substrates separated by a distance = h. Here the focus is specifically on a simple fluid in which the interaction between a pair of film molecules is governed by the Lennard-Jones (12,6) potential [33,58,59,77,79-84]. A confined simple fluid serves as a suitable model for approximately spherical OMCTS molecules confined... 

We report here some results for a simple model of a one-component fluid interacting via a slightly modified Lennard-Jones potential, with angular-dependent associative forces. The model is considered in contact with the adsorbing surface. The principal aim of the simulation is to investigate the... 

Adsorption of hard sphere fluid mixtures in disordered hard sphere matrices has not been studied profoundly and the accuracy of the ROZ-type theory in the description of the structure and thermodynamics of simple mixtures is difficult to discuss. Adsorption of mixtures consisting of argon with ethane and methane in a matrix mimicking silica xerogel has been simulated by Kaminsky and Monson [42,43] in the framework of the Lennard-Jones model. A comparison with experimentally measured properties has also been performed. However, we are not aware of similar studies for simpler hard sphere mixtures, but the work from our laboratory has focused on a two-dimensional partly quenched model of hard discs [44]. That makes it impossible to judge the accuracy of theoretical approaches even for simple binary mixtures in disordered microporous media. 

The first MC (16) and MD (17) studies were used to simulate the properties of single particle fluids. Although the basic MC (11,12) and MD (12,13) methods have changed little since the earliest simulations, the systems simulated have continually increased in complexity. The ability to simulate complex interfacial systems has resulted partly from improvements in simulation algorithms (15,18) or in the interaction potentials used to model solid surfaces (19). The major reason, however, for this ability has resulted from the increasing sophistication of the interaction potentials used to model liquid-liquid interactions. These advances have involved the use of the following potentials Lennard-Jones 12-6 (20), Rowlinson (21), BNS... 

In the present work, we performed MC simulations at different operation conditions, constant fluid density and constant pressure, for calculating K2 to investigate the distribution behavior in the supercritical region. We selected C02, benzene, and graphitic slitpore as a model system by adopting the Lennard - Jones (LJ) potential function for intermolecular interactions. 

Figure 1.38. Molecular dynamics simulation of the density profiles for spherical molecules in a cylinder, mimicking SFg in controlled pore glass (CPG-10). Fluid-fluid and fluid-wall interaction modelled by Lennard-Jones interactions. Reference A. de Keizer. T. Michalski and G.H. Findenegg, Pure Appl. Chem. 63(1991) 1495.

Water Is a strongly three-dlmenslonally structured fluid (sec. 1.5.3c) with structure-originating Interactions reaching several molecular diameters. Considering this, simple models and/or simulations with a limited number of molecules are not really helpful. By "simple" we mean models in which water molecules are represented as point dipoles, point quadrupoles, or as molecules with Lennard-Jones Interactions plus an additional dipole, etc., and by "limited" less than, say 10 molecules, i.e. 10 molecules in each direction of a cubic box. Admittedly, for a number of simpler problems more embryonic models may suffice. For example, electrochemists often get away with a dipole Interpretation when focusing their attention solely on the Stern layer polarization. Helmholtz s equations for the jf-potential 3.9.9] is an illustration. 

---