Vapour-liquid equilibria for two centre Lennard-Jones diatomics and dipolar diatomics
|
|
- Damian Barrett
- 5 years ago
- Views:
Transcription
1 MOLECULAR PHYSICS, 1993, VOL. 80, No. 4, Vapour-liquid equilibria for two centre Lennard-Jones diatomics and dipolar diatomics By GIRIJA S. DUBEYt, SEAMUS F. O'SHEA Department of Chemistry, University of Lethbridge, Lethbridge, Alberta, Canada TIK 3M4 and PETER A. MONSON Department of Chemical Engineering, University of Massachusetts, Amherst, Massachusetts 01003, USA (Received 4 January 1993; accepted 20 April 1993) Gibbs ensemble Monte Carlo simulations have been used to calculate vapour-liquid equilibria of diatomic and dipolar diatomic fluids. Dipolar two-centre Lennard-Jones diatomics with elongations L*= L/a between 0.33 and 1-00 have been studied for various values of dipole moments 1/~ ~--- p.d/(r 1/2. The results for phase equilibrium are in good agreement with recent theoretical predictions based on cluster-expansion perturbation theory. 1. Introduction A major advance in the use of molecular simulations in the study of vapourliquid equilibria has been the use of the Gibbs ensemble [1,2] method. This permits the direct simulation of bulk vapour and liquid phases satisfying the conditions of phase equilibrium. Gibbs ensemble Monte Carlo simulations have now been performed for a number of potential models, including pure Lennard-Jones fluids and mixtures [3-5], Stockmayer fluids and their mixtures, square-well fluids and polyatomic fluids and mixtures [6-9], the symmetrical non-additive hard-sphere system [10], the hard-core two-yukawa fluid [11], the Lennard-Jones fluids with a quadrupole interaction [12,13], hard-core Yukawa fluids [14], the non-spherical Gay- Berne fluid [15], square-well diatomics [16], more realistic models for alkanes and water [17,18], and hard ellipsoids [19]. The purpose of the present work is to apply the Gibbs ensemble method to some interaction site models of diatomic molecules. The models are two-centre Lennard- Jones potentials with point charge dipoles. These models can be used to investigate the influence of molecular shape and polarity upon vapour-liquid equilibria. Although the dipolar diatomic systems studied do not correspond directly to real polar diatomic fluids which would, of course, have heteronuclear repulsive cores, the parameters considered correspond roughly to moderately polar fluids, such as methyl chloride. Our study also affords the opportunity to test the accuracy of some recent theoretical predictions [20,21] obtained using a cluster-expansion perturbation theory developed for interaction site fluids [22]. This theory addresses the calculation t Present address: Department of Chemistry, New York University, New York, N.Y , USA /93 $ Taylor & Francis Ltd.
2 _ (~/~v2) 998 G.S. Dubey et al. of the contribution to the free energy arising from long-ranged contributions to the site-site potential, such as dispersion or Coulombic interactions. At present it is the only theory of interaction site systems with Coulombic site-site interactions capable of treating fluid properties over the range of densities required for prediction of vapour-liquid coexistence properties. However, it has not yet been subjected to extensive testing by comparison with molecular simulation results. Details of model potentials and simulations computational details are discussed in section 2 and the results of the Gibbs-ensemble simulations for phase-equilibria and critical behaviour of our simulation results as well as theoretical results are presented in section 3 and section 4 summarizes our conclusions. 2. Model potentials and computational details The focus of the present paper is diatomic and dipolar-diatomic molecules interacting with site-site Lennard-Jones 12-6 potentials and site-site Coulombic potentials. The total intermolecular pair potentials is thus of the form 2 2 u(1,2) = Z Z{qiqj/rij + 4e[(~r/rij )12 - (a/rij)6]}' (1) i j where rij is the distance between site i on one molecule and sitej on the other, qi is the charge on site i and ql = --q2 for charge neutrality. We used the Ewald summation method [23] to treat the Coulombic site-site interactions in the simulation so that li~lj~l [~ erfc (alrij + nl) Ew = qiqj Irij + n[ "= = I,,1=o +l/(~l3) Z qiqj (4~2/k2) exp (-k2/4a 2) cos (k. r/j)] kr / ~ ~ + (2~/(3L3)) ~ iq;rel2 + V~xlf, (2) i=1 i vself ~ / ~ ] (3) ex. = q2 where d is the intramolecular charge separation, here the internuclear distance, and a is a convergence parameter which controls the relative weights of the real space and reciprocal space sums, rfj is the vector between sites i and j, n is a vector whose components are (nxl, nyl, nzl); nx, ny and nz are integers, L is the length of the simulation cell and k = 2nn/L. When a partial charge formulation of the Ewald sum expression is used, as it is here, the reciprocal space terms include a contribution arising from the interaction between partial charges within an individual molecule and these must be removed from the total energy [24]; -ex.vself provides that correction. Our implementation of the Gibbs ensemble is essentially the same as that described by Panagiotopoulos and coworkers [3,4] and the reader is referred to these references for a more detailed description of the technique. Our calculations were performed mostly with N = 364 molecules, although some calculations with N = 512 particles were made. For the N = 364 molecule system we initialized the
3 Vapour-liquid equilibria 999 simulation with 256 molecules at a liquid-like density in one box and 108 molecules at a lower, vapour-like density for the other box. For the 512 molecule system there were 256 molecules in each phase initially. Generally the simulations were started either with the molecules of each box on the sites of an a-nitrogen lattice, or in an equilibrated state from a lower temperature. A simulation cycle in the Gibbs ensemble Monte Carlo method comprises three distinct types of moves: individual molecule displacement within a box, transfer of volume from one box to the other and particle transfers from one box to the other. The relative numbers of each type of move are chosen to optimize convergence. In these calculations a cycle consisted of a single attempted displacement of each molecule, followed by one trial volume transfer, and a number of attempted particle transfers. In the individual molecule displacements, the particles are chosen and displaced randomly, both translationally and orientationally (within the boxes) following the well established Metropolis scheme for canonical NVT simulations. Particle transfers were achieved by creating a particle at a random position in one subsystem and annihilating a randomly chosen particle in other region. The success rate for particle transfers was maintained in the range about 1-3%. The convergence and reliability of the results are very sensitive to the transfer rate; too low a transfer rate prevents equilibrium between the boxes, and too high a rate destroys equilibration within the boxes. The calculation required for the insertion step permitted the evaluation of the chemical potential by Widom's method [25] as generalized to the Gibbs ensemble [3]. The Lennard-Jones potential was truncated at an intermolecular separation equal to half the box length, and the usual long tail corrections were applied to energy, pressure and chemical potential. The simulations consisted of a minimum cycles of which were equilibration and Ncycle = 5000 for determination of ensemble averages. 3. Results and discussion The phase equilibrium results from our Gibbs ensemble Monte Carlo simulations are summarized in tables 1-7. The results are tabulated in reduced units defined as L* = L/o, T* = kt/e; E* = E/Nr p* =/9o3; P* = po'3/r /z~ :/.td/(r /Z] r = #l/nr and p~ = I~/N~. To avoid confusion we have used #~ for dipole moment, and labelled the chemical potential #r and/z~ for the liquid and vapour phases, respectively. The uncertainty reported with each result is the standard deviation of averages calculated over 100-cycle blocks during the production phase of the simulation and are only a rough guide to the true uncertainties which are more likely to be larger in most cases. In some cases simulations were repeated, starting from different initial conditions or with a different number of attempted molecule transfers per cycle. We found that the results obtained were the same within their combined uncertainties. Nevertheless these uncertainties must be treated with some care since, at best, they represent a 'local' estimate of the variability; the coexistence curves show roughness which is somewhat larger than the uncertainties but which is unlikely to be physically correct. Generally, we see that the pressures and chemical potentials calculated for the coexisting phases agree quite well for a given temperature. There is some disagree-
4 1000 G.S. Dubey et al. Table 1. Summary of the results of the Monte Carlo simulation of the Gibbs ensemble for the 12-6 diatomic for bond length L* = The symbols are defined in the text. L* = /4i* = 0.0 Vapour phase Liquid phase N T* Ncycl e /~ P* E~ ~ Pl* /)1" El* /4* 364 2" " ' '73 4-0" I " '506 0" " ' ' " t " " ' ' " ment at the lowest temperatures studied, but we believe this is principally a consequence of significant uncertainty in the calculated values of the virial pressures of liquids at those temperatures. We have used these vapour-liquid equilibrium data to estimate the critical temperatures and densities. For this purpose we have assumed that the coexistence densities follow the law of rectilinear diameters [26] and a power law near the critical point, i.e. (Pl* + p'v)~2 = AT* + B. (4) pl* - p* = C 1 T* fl (5) where p* and ~* are the saturation vapour and liquid densities at temperature T*,/3 is the critical exponent (we have assumed that/3 = 1/3), and A, B, C are constants. The values of the critical constants are presented in tables 8 and 9. Table 2. Vapour-liquid coexistence data from the Gibbs ensemble Monte Carlo simulation for L* = 0"6. The symbols are defined in the text. L* = 0'6 /~ = 0'0 Vapour phase Liquid phase N T* Ncycl r /9~ Pv* E~ /~ f~* P1* El* /zl* 364 2" ' "37 0"396 0" "02 -I t-0' "03 4-0" "95-8"32 0'377 0'038-8"96-8"29 4-0" " ' " ' "11 0"340 0" "06 4-0' " "01 +0" " " ' "02 4-0" "03 0"004 -t " " " " " "79-8' " " " "096-2' '36-8'01 4-0" ' '13 4-0' ' " " "07 4-0"03 4-0" " "13 4-0"02
5 Vapour-liquid equilibria 1 O T* ~- 2.0 \ \ \ \ \ \ 1.5 \, \ 1.0 \ 0.0 0' I J 4 i ~.7 Figure 1. Vapour-liquid coexistence curves for bond length L* = , 0-6 and 1"00 (from top to bottom, respectively). The lines are the theoretical results of McGuigan et al. [21], the symbols are our Gibbs ensemble results. The 9 and "k represent the critical points T* predicted by theory and Gibbs simulation, respectively. 2.4~ <> 2.o~ ~. T* 1.8~ r~ I Io 1" <> t 0~4 0.5 Figure 2. Vapour-liquid coexistence curve for bond length L*= 0-67, where <> are the molecular dynamic results of Gupta [27], x are our Gibbs ensemble results and 9 and represent the critical points T* predicted by MD and Gibbs simulation, respectively.
6 1002 G.S. Dubey et al / 0.25 / P* 0.20 / 0.10 / "0q ,0 T* Figure 3. Vapour pressure curves for diatomic fluid with reduced bond length L* = l-0, 0.6 and (from left to right, respectively). The lines are the theoretical results of McGuigan et al. [21] and the symbols are our Gibbs ensemble results. The coexisting densities for the non-polar 12-6 diatomics are plotted against temperature in figures 1,2. For the case of L*= 0-67, results from molecular dynamics [27] calculations are also shown and the agreement is seen to be excellent. For the other three cases we have also plotted predictions from cluster perturbation theory [21,22]. The agreement is again very good, with some deterioration for the shortest bond length at the highest temperature. The agreement is especially good for L* = 1.0. The somewhat poorer agreement at lower values of L* can be understood at least in part in terms of the cluster diagrams at each order in the density contributing to the free energy which are neglected by the theory of Lupkowski and Monson [20]. These diagrams make diminishing contributions to the free energy as L* increases. Figure 3 shows the comparison of simulation and theoretical results for the vapour pressures. Very good agreement is again apparent. Table 8 shows a Table 3. Vapour-liquid coexistence data from the Gibbs ensemble Monte Carlo simulation for L* = The symbols are defined in the text. L* /z~ = 0"0 Vapour phase Liquid phase N 7"* Ncycl e /9~ P* E~' /~ /~* PI* El* 14" " I "07 +0" " "039 0" l -I " "061-1"53-7" '53-7" I " I-0"07 i0.02
7 Vapour-liquid equilibria 1003 Table 4. Vapour-liquid coexistence data from the Gibbs ensemble Monte Carlo simulation for L* = 1-0. The symbols are defined in the text. L* = 1.0 /~ = 0.0 Vapour phase Liquid phase N 7"* Ncycle P~ P* E~ /~ PI* el* E~* /~* : / t t comparison of the critical densities and temperatures estimated from our Gibbs ensemble results together with predictions from cluster-expansion perturbation theory, [22], the first-order perturbation theory of Fischer et al. [28] and from a recently developed semi-empirical equation of state due to Sowers and Sandier [29]. For the one value of L*, 0"3292, for which data are available from both theories and from simulation, the predictions from the cluster perturbation theory for the critical temperature represent an improvement over those from the first-order perturbation theory which does not include the influence of attractive forces upon the fluid structure at low and moderate density. The two levels of perturbation theory provide equally accurate preductions of the critical density, but the corrections to perturbation theory have a bigger effect on the critical temperature than on the critical density. The predictions from the semiempirical equation of state are also reliable for these systems. Table 5. Vapour-liquid coexistence data from Gibbs ensemble Monte Carlo simulation for dipolar diatomics with bond length, L* = 0.6 and/~ = The symbols are defined in the text. L* = 0-6 /~ = 2.0 Vapour phase Liquid phase N T* Ncycl e /gv* Pv* ~ /d~ ill* P~ ~ //'1" 364 2" "012 0"020-0"57-9' '007-11"96-8" "07 -t t ' '040-1" "392 0"047-10" " " "033 0"052-1"33-9' "86-9" ' " " " " " " " " "83-8"75 0"320 0" "73 +0"006 +0"003 +0" "
8 1004 G.S. Dubey et al od ~ D o [] 2.6 T* I i 1" p* Figure 4. Vapour-liquid coexistence curve for L*= 0.60 and /~ = 3-00, L*= 0.60 and /z~ = 2.00, and L* = 1.00 and/z~ = 3-00 (from top to bottom, respectively). The lines are the theoretical results of Lupkowski and Monson [22], the symbols are our Gibbs ensemble results. The 9 and ~r represent the critical points T'predicted by theory and Gibbs simulation, respectively. Figure 4 shows the coexisting densities against temperature for the dipolar diatomic systems considered. The predictions from the cluster perturbation theory are also shown [22]. For L* = 0.6 the agreement is best for the smaller of the two dipole moments considered. For the L* = 1.0 the agreement is qualitatively good but the theory consistently underestimates the liquid density. Figure 5 shows the corresponding results for the vapour pressures. Table 9 shows a comparison of the critical Table 6. Vapour-liquid coexistence data from Gibbs ensemble Monte Carlo simulation for dipolar diatomics with bond length, L* = 0-6 and/~ = The symbols are defined in the text. L* = 0.6 ~* = 3'0 Vapour phase Liquid phase N 7* Ncycl e /9~ Pv* E~' /~ Pl* Pl* El* Pt* ' " t " " "341 0" " " " "03 +0" " "78 0"305 0" " " ' " "006 +0" " '74-10"71 0" ' " " "06 4-0'01 +0" t "68 0" " ' "02+0" "
9 Vapour-liquid equilibria ~ ~ ~ "... I l 0.15 P* 0.10k / / i o/ u / / / /.: / 0.05! / ; /, /; / 1.o ~15 22o z~ a.0 a T" Figure 5. Vapour pressure curves for dipolar diatomics fluid with reduced bond length L* = 1.00 and M' = 3.00, L* = 0.60 and ~ = 2.00 and L* = 0-60 and ~ = 3.00 (from left to fight, respectively). The lines are the theoretical results of Lupkowski and Monson [221 and symbols are our Gibbs ensemble results. densities and temperatures for dipolar diatomic fluids estimated from our Gibbs ensemble results together with predictions from cluster perturbation theory. The theory tends to underestimate the critical density and overestimate the critical temperature but the overall agreement is quite good. 4. Conclusions In this paper we have applied the Gibbs ensemble to the study of vapour-liquid Table 7. Vapour-liquid coexistence data from Gibbs ensemble Monte Carlo simulation for dipolar diatomics with bond length, L* = 1-00, and p~ = The symbols are defined in the text. L*= 1-0 /h~ = 3"0 Vapour phase Liquid phase N T* Ncycl e /~ P* E~ P-~ g* Pl* Ei* /4" 364 1' "317 0" " :0-07 5: :0"05 5:0" " " "009 5:0"001 5:0-07 5:0-06 5: :0"009 5:0-08 5:0" ' "11 0"263 0" " : :0-05 5: "08 5: " " " :0"007 5:0"003 5: : " " " " " : :0"001
10 1006 G.S. Dubey et al. Table 8. Comparison of the Gibbs simulation estimates of the critical temperature T* and critical density p'for various bond lengths L* with those of duster perturbation theory 9 [21], molecular dynamics of [27], the first-order perturbation theory of Fischer et al. [28], and semiempirical equation of state theory of Sowers and Sandier [29]. Reference L* T* p* Simulation (this work) Simulation [27] Cluster perturbation theory [21] First-order perturbation theory [28] Equation of state [29] Simulation (this work) Cluster perturbation theory [21] Simulation (this work) Simulation [27] First-order perturbation theory [28] Equation of state [29] Simulation (this work) Cluster perturbation theory [21] equilibria of diatomic and dipolar diatomic fluids. The principal focus of the work has been to test the predictions of the cluster expansion perturbation theory [20-22] which incorporates the influence of the attractive forces upon the fluid structure which is important at moderate and low densities, and the influence of multipolar interactions which is important at all densities. For the case of 12-6 diatomics the simulation results are in very good agreement with predictions from cluster perturbation theory. These predictions are comparable in accuracy to those obtained from the optimized cluster theory for the atomic 12-6 potential [30]. The result for 12-6 diatomics with a point-charge dipole are also very good for the case of low dipole strength. For high dipole strength the results are still in quite good agreement with the theory although there is tendency for the theory of underestimate the saturated liquid density. Overall we believe that this comparison is an important demonstration that prediction of vapour-liquid equilibria for systems with both non-spherical molecular shape and polarity lies within the grasp of currently available theoretical methods. Table 9. Comparison of the Gibbs simulation estimates of the critical temperature T'and critical density p* for various bond lengths L* and dipolar strengths/~ with those from cluster perturbation theory. Reference L* p~ T* /9* Simulation (this work) 0"60 2" Cluster perturbation theory [22] '73 Simulation (this work) Cluster perturbation theory [22] ' "150 Simulation (this work) ' Cluster perturbation theory [22] ' 124
11 Vapour-liquid equilibria 1007 Research on this problem at the University of Lethbridge was supported by a grant from the Natural Sciences and Engineering Research Council of Canada. S. F. O'Shea and P. A. Monson are grateful to Professor J. S. Rowlinson and the staff of the Physical Chemistry Laboratory, Oxford University for their hospitality during the period when this project was initiated, and to Professor D. J. Tildesley for his help and encouragement with the Gibbs ensemble. References [1] GUBmNS, K. E., 1989, Molec. Simul., 2, 223. [2] PANAGtOrOPOULOS, A. Z., and STAPLE'rON, M., 1989, Fluid Phase Equil., 53, 133. [3] PANAGIOrOPOULOS, A. Z., 1987, Molec. Phys., 61, 813; Ibid., 62, 701. [4] PANAGIOTOPOULOS, A. Z., QUIRKE, N., STAPLErON, M., and TILDESLEV, D. J., 1988, Molec. Phys., 63, 527. [5] SMIT, B., WILLIAMS, C. P., HENDRIKS, E. M., and OE LEEUW, S. W., 1989, Molec. Phys., 68, 765. [6] DE LEEVW, S. W., S~IT, B., and WILLrAX4S, C. W., 1990, J. chem. Phys., 93, [7] VE6A, L., OE MIGUEL, E., RULL, L. F., JACKSON, G., and McLuRE, I. A., 1992, J. chem. Phys., 96, [8] PANAGIOTOPOULOS, A. Z., 1989, Molec. Phys., 62, 701. [9] SINOH, R. P., PITZER, K. S., DE PABLO, J. J., and PRAUSNITZ, J. M., 1990, J. chem. Phys., 92, [10] AMAR, J. G., 1989, Molec. Phys., 67, 739. [11] RtJDISILL, E. N., and CUMMINGS, P. T., 1989, Molec. Phys., 68, 629. [12] STAPLETON, M. R., T~LDESLEV, D. J., PANAGIOTOPOULOS, A. Z., and QUmKE, N., 1989, Molec. Simul., 2, 147. [13] SUIT, B., and WILLIAM, C. P., 1990, J. Phys.: Condens. Matter, 2, [14] SUIT, B., and FRENKEL, D., 1991, Molec. Phys., 74, 35. [15] DE MIGVEL, E., RULL, L. F., CHALAM, M. K., and GUaBINS, K. E., 1990, Molec. Phys., 71, [16] YETHIaAJ, A., and HALL, C. K., 1991, Molec. Phys., 72, 619. [17] DE PAaLO, J. J., and Pa_~USNIrZ, J. M., 1989, Fluid Phase Equil., 53, 177. [18] DE PAaLO, J. J., PRAUSNITZ, J. M., STRAUCn, H. J., and CUMMINGS, P. T., 1990, J. chem. Phys., 93, [19] DE MIGUEL, E. and ALLEN, M. P., 1992, Molec. Phys., 76, [20] LUPKOWSKI, M., and MONSON, P. A., 1987, J. chem. Phys., 87, [21] McGtJ~GAN, D. B., LUPKOWSKI, M., PAQUET, D., and Mo~soN, P. A., 1989, Molec. Phys., 67, 33. [22] LupKowsgI, M. and MONSON, P. A., 1989, Molec. Phys., 67, 53. [23] EWALD, P. P., 1921, Ann. Phys., 64, 253. [24] HEYES, D. M., 1983, CCP5 Quarterly, 8, 29. [25] WIDOM, B., 1963, J. chem. Phys., 39, [26] GUGGENHEIM, E. A. J., 1945, J. chem. Phys., 13, 253. [27] Gu~rA, S., 1988, J. phys. Chem., 92, [28] FISCHER, J., LUSTING, R., BREITENFELDER-MANSKE, H., and LEMMING, W., 1984, Molec. Phys., 52, 485. [29] SowEP.S, (3. M., and SANDLER, S. I., 1992, Molec. Phys., 77, 351. [30] StrNG, S. H., and CHANDLER, D., 1974, Phys. Rev. A, 9, 1688.
Phase transitions of quadrupolar fluids
Phase transitions of quadrupolar fluids Seamus F. O Shea Department of Chemistry, University of Lethbridge, Lethbridge, Alberta, Canada, T1K 3M4 Girija S. Dubey Brookhaven National Laboratory, Upton, New
More informationGibbs ensemble simulation of phase equilibrium in the hard core two-yukawa fluid model for the Lennard-Jones fluid
MOLECULAR PHYSICS, 1989, VOL. 68, No. 3, 629-635 Gibbs ensemble simulation of phase equilibrium in the hard core two-yukawa fluid model for the Lennard-Jones fluid by E. N. RUDISILL and P. T. CUMMINGS
More informationComputer simulation of vapor-liquid equilibria of linear dipolar fluids: Departures from the principle of corresponding states
Computer simulation of vapor-liquid equilibria of linear dipolar fluids: Departures from the principle of corresponding states Benito Garzón, Santiago Lago, a) and Carlos Vega Departamento de Química Física,
More informationCalculation of the chemical potential in the Gibbs ensemble
MOLECULAR PHYSICS, 1989, VOL. 68, NO. 4, 951-958 Calculation of the chemical potential in the Gibbs ensemble by B. SMIT Koninklijke/Shell-Laboratorium, Amsterdam, (Shell Research B.V.) Badhuisweg 3, 1031
More informationarxiv: v1 [physics.comp-ph] 22 Apr 2009
Comprehensive study of the vapour-liquid equilibria of the pure two-centre Lennard-Jones plus pointdipole fluid Jürgen Stoll, Jadran Vrabec, Hans Hasse arxiv:0904.3413v1 [physics.comp-ph] 22 Apr 2009 Institut
More informationJoy Billups Woller A THESIS. Presented to the Faculty of. The Graduate College in the University of Nebraska
DETERMINATION OF THE THERMODYNAMIC PROPERTIES OF FLUIDS BY GIBBS ENSEMBLE MONTE CARLO COMPUTER SIMULATIONS by Joy Billups Woller A THESIS Presented to the Faculty of The Graduate College in the University
More informationDIRECT EVALUATION OF VAPOUR-LIQUID EQUILIBRIA BY MOLECULAR DYNAMICS USING GIBBS-DUHEM INTEGRATION
Molecular Simulation, 1996, Vol. 17, pp. 21-39 Reprints available directly from the publisher Photocopying permitted by license only 0 1996 OPA (Overseas Publishers Association) Amsterdam B.V. Published
More informationMolecular dynamics simulation of the liquid vapor interface: The Lennard-Jones fluid
Molecular dynamics simulation of the liquid vapor interface: The Lennard-Jones fluid Matthias Mecke and Jochen Winkelmann a) Institut für Physikalische Chemie, Universität Halle-Wittenberg, Geusaer Str.,
More informationUB association bias algorithm applied to the simulation of hydrogen fluoride
Fluid Phase Equilibria 194 197 (2002) 249 256 UB association bias algorithm applied to the simulation of hydrogen fluoride Scott Wierzchowski, David A. Kofke Department of Chemical Engineering, University
More informationOn the Calculation of the Chemical Potential. Using the Particle Deletion Scheme
On the Calculation of the Chemical Potential Using the Particle Deletion Scheme Georgios C. Boulougouris,2, Ioannis G. Economou and Doros. Theodorou,3,* Molecular Modelling of Materials Laboratory, Institute
More informationEffect of a density-dependent potential on the phase behaviour of fluids
MOLECULAR PHYSICS, 1992, VOL. 77, No. 6, 1021-1031 Effect of a density-dependent potential on the phase behaviour of fluids By B. SMIT Koninklijke/Shell-Laboratorium, Amsterdam (Shell Research B.V.), PO
More informationMIHALY MEZEI. Department of Physiology and Biophysics, Mount Sinai School of Medicine, CUNY. New York, NY 10029
Molecular Simulation, 1992, Vol. 9, pp, 257-261 c Gordon and Breach THEORETICAL CALCULATION OF THE LIQUID-VAPOR COEXISTENCE CURVE OF WATER, CHLOROFORM AND METHANOL WITH THECAVITY-BIASED MONTE CARLO METHOD
More informationSupplemental Material for Temperature-sensitive colloidal phase behavior induced by critical Casimir forces
Supplemental Material for Temperature-sensitive colloidal phase behavior induced by critical Casimir forces Minh Triet Dang, 1 Ana Vila Verde, 2 Van Duc Nguyen, 1 Peter G. Bolhuis, 3 and Peter Schall 1
More informationCE 530 Molecular Simulation
1 CE 530 Molecular Simulation Lecture 14 Molecular Models David A. Kofke Department of Chemical Engineering SUNY Buffalo kofke@eng.buffalo.edu 2 Review Monte Carlo ensemble averaging, no dynamics easy
More information510 Subject Index. Hamiltonian 33, 86, 88, 89 Hamilton operator 34, 164, 166
Subject Index Ab-initio calculation 24, 122, 161. 165 Acentric factor 279, 338 Activity absolute 258, 295 coefficient 7 definition 7 Atom 23 Atomic units 93 Avogadro number 5, 92 Axilrod-Teller-forces
More informationAn explicit expression for finite-size corrections to the chemical potential
J. Phys.: Condens. Matter 1 (1989) 8659-8665. Printed in the UK An explicit expression for finite-size corrections to the chemical potential B Smitt and D Frenkelt t Koninklijke/Shell-Laboratorium, Amsterdam
More informationCE 530 Molecular Simulation
1 CE 530 Molecular Simulation Lecture 1 David A. Kofke Department of Chemical Engineering SUNY Buffalo kofke@eng.buffalo.edu 2 Time/s Multi-Scale Modeling Based on SDSC Blue Horizon (SP3) 1.728 Tflops
More informationMolecular simulation of adsorption from dilute solutions
Vol. 52 No. 3/2005 685 689 on-line at: www.actabp.pl Molecular simulation of adsorption from dilute solutions Werner Billes Rupert Tscheliessnig and Johann Fischer Institut für Verfahrens- und Energietechnik
More informationPhase Equilibria of binary mixtures by Molecular Simulation and PR-EOS: Methane + Xenon and Xenon + Ethane
International Journal of ChemTech Research CODEN( USA): IJCRGG ISSN : 0974-4290 Vol.5, No.6, pp 2975-2979, Oct-Dec 2013 Phase Equilibria of binary mixtures by Molecular Simulation and PR-EOS: Methane +
More informationAnalytic studies of the hard dumbell fluid III. Zero pole approximation for the hard dumbell Yukawa fluid
MOLECULAR PHYSICS, 1984, VOL. 53, No. 4, 849-863 Analytic studies of the hard dumbell fluid III. Zero pole approximation for the hard dumbell Yukawa fluid by P. T. CUMMINGS Department of Chemical Engineering,
More informationThermodynamics of Three-phase Equilibrium in Lennard Jones System with a Simplified Equation of State
23 Bulletin of Research Center for Computing and Multimedia Studies, Hosei University, 28 (2014) Thermodynamics of Three-phase Equilibrium in Lennard Jones System with a Simplified Equation of State Yosuke
More informationThermodynamics and structural properties of the dipolar Yukawa fluid. Citation Journal Of Chemical Physics, 1999, v. 111 n. 1, p.
Title Thermodynamics and structural properties of the dipolar Yukawa fluid Author(s) Szalai, I; Henderson, D; Boda, D; Chan, KY Citation Journal Of Chemical Physics, 1999, v. 111 n. 1, p. 337-344 Issued
More informationEvaluation of the CPY and PYX approximations for short ranged anisotropie potentials
MOLECULAR PHYSICS, 1983, VOL. 50, NO. 5, 1133-1140 Evaluation of the CPY and PYX approximations for short ranged anisotropie potentials by P. T. CUMMINGS t Departments of Mechanical Engineering and Chemistry,
More informationPhase equilibria and critical behavior of square-well fluids of variable width by Gibbs ensemble Monte Carlo simulation
Phase equilibria and critical behavior of square-well fluids of variable width by Gibbs ensemble Monte Carlo simulation Lourdes Vega, Enrique de Miguel, and Luis F. Ruil Departamento de fisica Atbmica
More informationPhase Equilibria and Molecular Solutions Jan G. Korvink and Evgenii Rudnyi IMTEK Albert Ludwig University Freiburg, Germany
Phase Equilibria and Molecular Solutions Jan G. Korvink and Evgenii Rudnyi IMTEK Albert Ludwig University Freiburg, Germany Preliminaries Learning Goals Phase Equilibria Phase diagrams and classical thermodynamics
More informationThermodynamic behaviour of mixtures containing CO 2. A molecular simulation study
Thermodynamic behaviour of mixtures containing. A molecular simulation study V. Lachet, C. Nieto-Draghi, B. Creton (IFPEN) Å. Ervik, G. Skaugen, Ø. Wilhelmsen, M. Hammer (SINTEF) Introduction quality issues
More informationThe condition of microscopic reversibility in Gibbs ensemble Monte Carlo simulations of phase equilibria
MOLECULAR PHYSICS, 1995, VOL. 85, NO. 3, 435-447 The condition of microscopic reversibility in Gibbs ensemble Monte Carlo simulations of phase equilibria By LUIS F. RULL, GEORGE JACKSON]" Departamento
More informationA study of orientational ordering in a fluid of dipolar Gay Berne molecules using density-functional theory
JOURNAL OF CHEMICAL PHYSICS VOLUME 116, NUMBER 20 22 MAY 2002 A study of orientational ordering in a fluid of dipolar Gay Berne molecules using density-functional theory Szabolcs Varga a) Department of
More informationEQUATION OF STATE DEVELOPMENT
EQUATION OF STATE DEVELOPMENT I. Nieuwoudt* & M du Rand Institute for Thermal Separation Technology, Department of Chemical Engineering, University of Stellenbosch, Private bag X1, Matieland, 760, South
More informationUnusual Entropy of Adsorbed Methane on Zeolite Templated Carbon. Supporting Information. Part 2: Statistical Mechanical Model
Unusual Entropy of Adsorbed Methane on Zeolite Templated Carbon Supporting Information Part 2: Statistical Mechanical Model Nicholas P. Stadie*, Maxwell Murialdo, Channing C. Ahn, and Brent Fultz W. M.
More informationThe Effect of Model Internal Flexibility Upon NEMD Simulations of Viscosity
Draft: September 29, 1999 The Effect of Model Internal Flexibility Upon NEMD Simulations of Viscosity N. G. Fuller 1 and R. L. Rowley 1,2 Abstract The influence of model flexibility upon simulated viscosity
More informationComputer simulations of soft repulsive spherocylinders
MOLECULAR PHYSICS, 2001, VOL. 99, NO. 20, 1719± 1726 Computer simulations of soft repulsive spherocylinders DAVID J. EARL, JAROSLAV ILNYTSKYI and MARK R. WILSON* Department of Chemistry, University of
More informationA New Intermolecular Potential Model for the n- Alkane Homologous Series
A New Intermolecular Potential Model for the n- Alkane Homologous Series Jeffrey R. Errington and Athanassios Z. Panagiotopoulos School of Chemical Engineering Cornell University Ithaca, NY 14853 and Institute
More informationReactive canonical Monte Carlo A new simulation technique for reacting or associating fluids
MOLECULAR PHYSICS, 1994, VOL. 81, No. 3, 717-733 Reactive canonical Monte Carlo A new simulation technique for reacting or associating fluids By J. KARL JOHNSON Complex Systems Theory Branch, Naval Research
More informationGibbs Ensemble Computer Simulations of Vapor Liquid Equilibrium of Hydrazine
CHME 498 Undergraduate Research Final Report New Mexico State University Gibbs Ensemble Computer Simulations of Vapor Liquid Equilibrium of Hydrazine Summer 2017 Gustav Barraza Faculty Adviser: Dr. Martha
More informationComputer simulation methods (2) Dr. Vania Calandrini
Computer simulation methods (2) Dr. Vania Calandrini in the previous lecture: time average versus ensemble average MC versus MD simulations equipartition theorem (=> computing T) virial theorem (=> computing
More informationRESEARCH NOTE. DOUGLAS HENDERSON? Department of Chemistry and Biochemistry, Brigham Young University, Provo, Utah USA
MOLECULAR PHYSICS, 1999, VOL. 96, No. 7, 1145-1149 RESEARCH NOTE A simple theory for the partial molar volumes of a binary mixture DOUGLAS HENDERSON? Department of Chemistry and Biochemistry, Brigham Young
More informationPORE SIZE DISTRIBUTION OF CARBON WITH DIFFERENT PROBE MOLECULES
PORE SIZE DISTRIBUTION OF CARBON WITH DIFFERENT PROBE MOLECULES Atichat Wongkoblap*, Worapot Intomya, Warangkhana Somrup, Sorod Charoensuk, Supunnee Junpirom and Chaiyot Tangsathitkulchai School of Chemical
More informationSurface tension of the two center Lennard-Jones plus quadrupole model fluid
Surface tension of the two center Lennard-Jones plus quadrupole model fluid Stephan Werth, Martin Horsch 1, Hans Hasse Laboratory of Engineering Thermodynamics, Department of Mechanical and Process Engineering,
More informationTransferable Potentials for Phase Equilibria. 1. United-Atom Description of n-alkanes
J. Phys. Chem. B 1998, 102, 2569-2577 2569 Transferable Potentials for Phase Equilibria. 1. United-Atom Description of n-alkanes Marcus G. Martin and J. Ilja Siepmann* Department of Chemistry, UniVersity
More informationHunter College of the CUNY. New York, NY. USA PLEASE SCROLL DOWN FOR ARTICLE
This article was downloaded by:[new York University] [New York University] On: 16 July 2007 Access Details: [subscription number 769426389] Publisher: Taylor & Francis Informa Ltd Registered in England
More informationIntermolecular Forces and Monte-Carlo Integration 열역학특수연구
Intermolecular Forces and Monte-Carlo Integration 열역학특수연구 2003.3.28 Source of the lecture note. J.M.Prausnitz and others, Molecular Thermodynamics of Fluid Phase Equiliria Atkins, Physical Chemistry Lecture
More informationCH.7 Fugacities in Liquid Mixtures: Models and Theories of Solutions
CH.7 Fugacities in Liquid Mixtures: Models and Theories of Solutions The aim of solution theory is to express the properties of liquid mixture in terms of intermolecular forces and liquid structure. The
More informationIntroduction Statistical Thermodynamics. Monday, January 6, 14
Introduction Statistical Thermodynamics 1 Molecular Simulations Molecular dynamics: solve equations of motion Monte Carlo: importance sampling r 1 r 2 r n MD MC r 1 r 2 2 r n 2 3 3 4 4 Questions How can
More informationMonte Carlo simulations of vapor-liquid-liquid equilibrium of ternary
1 Monte Carlo simulations of vapor-- equilibrium of ternary petrochemical mixtures Suren Moodley a, Kim Bolton b, and Deresh Ramjugernath a a Thermodynamics Research Unit, School of Chemical Engineering,
More informationIntroduction to molecular dynamics
1 Introduction to molecular dynamics Yves Lansac Université François Rabelais, Tours, France Visiting MSE, GIST for the summer Molecular Simulation 2 Molecular simulation is a computational experiment.
More informationAn Extended van der Waals Equation of State Based on Molecular Dynamics Simulation
J. Comput. Chem. Jpn., Vol. 8, o. 3, pp. 97 14 (9) c 9 Society of Computer Chemistry, Japan An Extended van der Waals Equation of State Based on Molecular Dynamics Simulation Yosuke KATAOKA* and Yuri YAMADA
More informationMolecular modelling and simulation of hydrogen bonding fluids
Molecular modelling and simulation of hydrogen bonding fluids Martin Horsch, Alptekin Celik and Hans Hasse Lehrstuhl für Thermodynamik (LTD) Technische Universität Kaiserslautern Rostock, 27 th March 13
More informationChem 4501 Introduction to Thermodynamics, 3 Credits Kinetics, and Statistical Mechanics
Chem 4501 Introduction to hermodynamics, 3 Credits Kinetics, and Statistical Mechanics Module Number 2 Active Learning Answers and Optional Problems/Solutions 1. McQuarrie and Simon, 2-6. Paraphrase: How
More informationNONADDITIVITY OF ORIENTATIONAL INTERACTIONS AND PERTURBATION THEORY FOR DIPOLE HARD SPHERES INTRODUCTION. G. B. Litinskii UDC :536.
Journal of Structural Chemistry. Vol. 51, No. 6, pp. 1081-1089, 010 Original Russian Text Copyright 010 by G. B. Litinskii NONADDITIVITY OF ORIENTATIONAL INTERACTIONS AND PERTURBATION THEORY FOR DIPOLE
More informationMayer-sampling Monte Carlo calculations of methanol virial. Katherine R. S. Shaul, Andrew J. Schultz, David A. Kofke*
Mayer-sampling Monte Carlo calculations of methanol virial coefficients Katherine R. S. Shaul, Andrew J. Schultz, David A. Kofke* Department of Chemical & Biological Engineering, University at Buffalo,
More informationCARBON 2004 Providence, Rhode Island. Adsorption of Flexible n-butane and n-hexane on Graphitized Thermal Carbon Black and in Slit Pores
CARBON Providence, Rhode Island Adsorption of Flexible n-butane and n-hexane on Graphitized Thermal Carbon Black and in Slit Pores D. D. Do* and H. D. Do, University of Queensland, St. Lucia, Qld 7, Australia
More informationInhomogeneous 2D Lennard Jones Fluid: Theory and Computer Simulation
Commun. Theor. Phys. 58 (212) 759 764 Vol. 58, No. 5, November 15, 212 Inhomogeneous 2D Lennard Jones Fluid: Theory and Computer Simulation R. Khordad Department of Physics, College of Sciences, Yasouj
More informationTitle Super- and subcritical hydration of Thermodynamics of hydration Author(s) Matubayasi, N; Nakahara, M Citation JOURNAL OF CHEMICAL PHYSICS (2000), 8109 Issue Date 2000-05-08 URL http://hdl.handle.net/2433/50350
More informationCritical temperature of in nitely long chains from Wertheim s perturbation theory
MOLECULAR PHYSICS, 2000, VOL. 98, NO. 17, 1295 ± 1308 Critical temperature of in nitely long chains from Wertheim s perturbation theory C. VEGA* and L. G. MACDOWELL Departamento de QuõÂ mica Fõ Â sica,
More informationINTERMOLECULAR FORCES
INTERMOLECULAR FORCES Their Origin and Determination By GEOFFREY C. MAITLAND Senior Research Scientist Schlumberger Cambridge Research, Cambridge MAURICE RIGBY Lecturer in the Department of Chemistry King's
More informationTHE DETAILED BALANCE ENERGY-SCALED DISPLACEMENT MONTE CARLO ALGORITHM
Molecular Simulation, 1987, Vol. 1, pp. 87-93 c Gordon and Breach Science Publishers S.A. THE DETAILED BALANCE ENERGY-SCALED DISPLACEMENT MONTE CARLO ALGORITHM M. MEZEI Department of Chemistry, Hunter
More informationEquations of State. Equations of State (EoS)
Equations of State (EoS) Equations of State From molecular considerations, identify which intermolecular interactions are significant (including estimating relative strengths of dipole moments, polarizability,
More informationWater structure near single and multi-layer nanoscopic hydrophobic plates of varying separation and interaction potentials
Bull. Mater. Sci., Vol. 31, No. 3, June 2008, pp. 525 532. Indian Academy of Sciences. Water structure near single and multi-layer nanoscopic hydrophobic plates of varying separation and interaction potentials
More informationMolecular modeling of hydrogen bonding fluids: vapor-liquid coexistence and interfacial properties
12 th HLRS Results and Review Workshop Molecular modeling of hydrogen bonding fluids: vapor-liquid coexistence and interfacial properties High Performance Computing Center Stuttgart (HLRS), October 8,
More informationIMPROVED METHOD FOR CALCULATING SURFACE TENSION AND APPLICATION TO WATER
IMPROVED METHOD FOR CALCULATING SURFACE TENSION AND APPLICATION TO WATER ABSTRACT Hong Peng 1, Anh V Nguyen 2 and Greg R Birkett* School of Chemical Engineering, The University of Queensland Brisbane,
More informationUnderstanding Molecular Simulation 2009 Monte Carlo and Molecular Dynamics in different ensembles. Srikanth Sastry
JNCASR August 20, 21 2009 Understanding Molecular Simulation 2009 Monte Carlo and Molecular Dynamics in different ensembles Srikanth Sastry Jawaharlal Nehru Centre for Advanced Scientific Research, Bangalore
More informationMonte Carlo study of quadrupolar Kihara fluids
MOLECULAR PHYSICS, 1992, VOL. 72, NO. 4, 881-895 Monte Carlo study of quadrupolar Kihara fluids By CARLOS VEGA Departamento de Quimica Fisica, Facultad de Ciencias Quimicas, Universidad Complutense de
More informationMonte Carlo Methods in Statistical Mechanics
Monte Carlo Methods in Statistical Mechanics Mario G. Del Pópolo Atomistic Simulation Centre School of Mathematics and Physics Queen s University Belfast Belfast Mario G. Del Pópolo Statistical Mechanics
More informationIntroduction. Monday, January 6, 14
Introduction 1 Introduction Why to use a simulation Some examples of questions we can address 2 Molecular Simulations Molecular dynamics: solve equations of motion Monte Carlo: importance sampling Calculate
More informationMolecular Simulation Study of the Vapor Liquid Interfacial Behavior of a Dimer-forming Associating Fluid
Molecular Simulation, Vol. 30 (6), 15 May 2004, pp. 343 351 Molecular Simulation Study of the Vapor Liquid Interfacial Behavior of a Dimer-forming Associating Fluid JAYANT K. SINGH and DAVID A. KOFKE*
More informationChain Length Effects of Linear Alkanes in Zeolite Ferrierite. 2. Molecular Simulations 1
3952 J. Phys. Chem. B 1998, 102, 3952-3958 Chain Length Effects of Linear Alkanes in Zeolite Ferrierite. 2. Molecular Simulations 1 Willy J. M. van Well, Xavier Cottin, Berend Smit, Jan H. C. van Hooff,
More informationImperfect Gases. NC State University
Chemistry 431 Lecture 3 Imperfect Gases NC State University The Compression Factor One way to represent the relationship between ideal and real gases is to plot the deviation from ideality as the gas is
More informationREPORT DOCUMENTATION PAGE
REPORT DOCUMENTATION PAGE Form Approved OMB NO. 0704-0188 Public Reporting burden for this collection of information is estimated to average 1 hour per response, including the time for reviewing instructions,
More informationMelting line of the Lennard-Jones system, infinite size, and full potential
THE JOURNAL OF CHEMICAL PHYSICS 127, 104504 2007 Melting line of the Lennard-Jones system, infinite size, and full potential Ethan A. Mastny a and Juan J. de Pablo b Chemical and Biological Engineering
More informationMultiple time step Monte Carlo
JOURNAL OF CHEMICAL PHYSICS VOLUME 117, NUMBER 18 8 NOVEMBER 2002 Multiple time step Monte Carlo Balázs Hetényi a) Department of Chemistry, Princeton University, Princeton, NJ 08544 and Department of Chemistry
More informationTitle Theory of solutions in the energy r of the molecular flexibility Author(s) Matubayasi, N; Nakahara, M Citation JOURNAL OF CHEMICAL PHYSICS (2003), 9702 Issue Date 2003-11-08 URL http://hdl.handle.net/2433/50354
More informationThermophysical Properties from the Equation of State of Mason and Co-workers ~
hlternational Journal of Tlwrmophysic.~ ~. I.'ol. IA', No. 4. 1997 Thermophysical Properties from the Equation of State of Mason and Co-workers ~ A. Boushehri 2" 3 and H. Eslami 2 The theory gives formuhts
More informationComputer simulation of vapor-liquid equilibria of linear quadrupolar fluids. Departures from the principle of corresponding states
Computer simulation of vapor-liquid equilibria of linear quadrupolar fluids. Departures from the principle of corresponding states Benito Garzck, Santiago Lago,a) and Carlos Vega Departamento de Quimica
More informationMonte Carlo methods for phase equilibria of uids 2 1. Introduction The phase behavior of uids and their mixtures is of central importance to many tech
TOPICAL REVIEW Monte Carlo methods for phase equilibria of uids Athanassios Z. Panagiotopoulos z Institute for Physical Science and Technology and Department of Chemical Engineering, University of Maryland,
More informationCritical point of electrolyte mixtures
Critical point of electrolyte mixtures THE JOURNAL OF CHEMICAL PHYSICS 123, 084903 2005 Antti-Pekka Hynninen and Marjolein Dijkstra Soft Condensed Matter Group, Debye Institute, Utrecht University, Princetonplein
More informationSchool of Chemical Engineering, Universidad del Valle, Ciudad Universitaria Melendez, Building 336, Apartado 25360, Cali, Colombia
Vapor-phase chemical equilibrium and combined chemical and vapor-liquid equilibrium for the ternary system ethylene + water + ethanol from reaction-ensemble and reactive Gibbs-ensemble molecular simulations
More informationMonte Carlo Calculations of Effective Surface Tension for Small Clusters
Monte Carlo Calculations of Effective Surface Tension for Small Clusters Barbara N. Hale Physics Department and Cloud and Aerosol Science Laboratory, University of Missouri- Rolla, Rolla, MO 65401, USA
More informationEdinburgh Research Explorer
Edinburgh Research Explorer Applications of Wang-Landau sampling to determine phase equilibria in complex fluids Citation for published version: Ganzenmueller, G & Camp, PJ 2007, 'Applications of Wang-Landau
More informationN-body simulations. Phys 750 Lecture 10
N-body simulations Phys 750 Lecture 10 N-body simulations Forward integration in time of strongly coupled classical particles or (rigid-body composite objects) dissipation Brownian motion Very general
More informationLennard-Jones as a model for argon and test of extended renormalization group calculations
JOURNAL OF CHEMICAL PHYSICS VOLUME 111, NUMBER 2 22 NOVEMBER 1999 Lennard-Jones as a model for argon and test of extended renormalization group calculations John A. White Department of Physics, American
More informationBinary Hard-Sphere Mixtures Within Spherical Pores
Journal of the Korean Physical Society, Vol. 35, No. 4, October 1999, pp. 350 354 Binary Hard-Sphere Mixtures Within Spherical Pores Soon-Chul Kim Department of Physics, Andong National University, Andong
More informationThe change in free energy on transferring an ion from a medium of low dielectric constantε1 to one of high dielectric constant ε2:
The Born Energy of an Ion The free energy density of an electric field E arising from a charge is ½(ε 0 ε E 2 ) per unit volume Integrating the energy density of an ion over all of space = Born energy:
More informationSystematic Coarse-Graining and Concurrent Multiresolution Simulation of Molecular Liquids
Systematic Coarse-Graining and Concurrent Multiresolution Simulation of Molecular Liquids Cameron F. Abrams Department of Chemical and Biological Engineering Drexel University Philadelphia, PA USA 9 June
More informationCE 530 Molecular Simulation
CE 530 Molecular Simulation Lecture 20 Phase Equilibria David A. Kofke Department of Chemical Engineering SUNY Buffalo kofke@eng.buffalo.edu 2 Thermodynamic Phase Equilibria Certain thermodynamic states
More informationVapor liquid equilibria of linear and branched alkanes from perturbation theory
JOURNAL OF CHEMICAL PHYSICS VOLUME 109, NUMBER 13 1 OCTOBER 1998 Vapor liquid equilibria of linear and branched alkanes from perturbation theory L. G. MacDowell and C. Vega Departamento de Química Física,
More informationThe broad topic of physical metallurgy provides a basis that links the structure of materials with their properties, focusing primarily on metals.
Physical Metallurgy The broad topic of physical metallurgy provides a basis that links the structure of materials with their properties, focusing primarily on metals. Crystal Binding In our discussions
More informationA MOLECULAR DYNAMICS SIMULATION OF WATER DROPLET IN CONTACT WITH A PLATINUM SURFACE
The 6th ASME-JSME Thermal Engineering Joint Conference March 16-2, 23 TED-AJ3-183 A MOLECULAR DYNAMICS SIMULATION OF WATER DROPLET IN CONTACT WITH A PLATINUM SURFACE Tatsuto KIMURA Department of Mechanical
More informationSimple model of sickle hemogloblin
Simple model of sickle hemogloblin THE JOURNAL OF CHEMICAL PHYSICS 125, 024902 2006 Andrey Shiryayev, a Xiaofei Li, and J. D. Gunton b Department of Physics, Lehigh University, Bethlehem, Pennsylvania
More informationChain length effects on aqueous alkane solubility near the solvent s critical point
Fluid Phase Equilibria 183 184 (2001) 289 294 Chain length effects on aqueous alkane solubility near the solvent s critical point Eric M. Yezdimer a,b,, Ariel A. Chialvo c,d, Peter T. Cummings b,e a Department
More informationChapter 2 Experimental sources of intermolecular potentials
Chapter 2 Experimental sources of intermolecular potentials 2.1 Overview thermodynamical properties: heat of vaporization (Trouton s rule) crystal structures ionic crystals rare gas solids physico-chemical
More informationAdsorption of Lennard-Jones Fluids in Carbon Slit Pores of a Finite Length. AComputer Simulation Study
1 Invited Contribution Adsorption of Lennard-Jones Fluids in Carbon Slit Pores of a Finite Length. AComputer Simulation Study A. Wongkoblap 1, S. Junpirom 2 and D.D. Do 1 * (1) Department of Chemical Engineering,
More informationGas-liquid phase separation in oppositely charged colloids: stability and interfacial tension
7 Gas-liquid phase separation in oppositely charged colloids: stability and interfacial tension We study the phase behaviour and the interfacial tension of the screened Coulomb (Yukawa) restricted primitive
More informationAuthor's personal copy
Adsorption (2011) 17: 265 271 DOI 10.1007/s10450-011-9325-7 Monte Carlo simulation of polymer adsorption Christopher J. Rasmussen Aleksey Vishnyakov Alexander V. Neimark Received: 15 May 2010 / Accepted:
More informationABSTRACT. Sarvin Moghaddam, Doctor of Philosophy, 2003
ABSTRACT Title of Dissertation: FINELY-DISCRETIZED LATTICE MODELS FOR THERMODYNAMIC PROPERTIES OF FLUIDS Sarvin Moghaddam, Doctor of Philosophy, 2003 Dissertation directed by: Professor Mikhail Anisimov
More informationSimulation Methods II
Simulation Methods II Maria Fyta Institute for Computational Physics Universität Stuttgart Summer Term 2018 SM II - contents First principles methods Hartree-Fock and beyond Density-funtional-theory Ab
More informationLocal molecular field theory for effective attractions between like charged objects in systems with strong Coulomb interactions
Local molecular field theory for effective attractions between like charged objects in systems with strong Coulomb interactions Yng-Gwei Chen and John D. Weeks Departments of Physics and Chemistry and
More informationTheory of infinite dilution chemical potential
Fluid Phase Equilibria Journal Volume 85, 141-151, 1993 141 Esam Z. Hamada and G.Ali Mansoorib, * a Department of Chemical Engineering, King Fahd University of Petroleum and Minerals, Dhahran 31261 (Saudi
More informationA new algorithm for Reverse Monte Carlo simulations
A new algorithm for Reverse Monte Carlo simulations Fernando Lus B. da Silva, Bo Svensson, Torbjörn Åkesson, and Bo Jönsson Citation: The Journal of Chemical Physics 109, 2624 (1998); doi: 10.1063/1.476861
More information