Evaluation of the CPY and PYX approximations for short ranged anisotropie potentials
|
|
- Georgia Burns
- 6 years ago
- Views:
Transcription
1 MOLECULAR PHYSICS, 1983, VOL. 50, NO. 5, Evaluation of the CPY and PYX approximations for short ranged anisotropie potentials by P. T. CUMMINGS t Departments of Mechanical Engineering and Chemistry, State University of New York at Stony Brook, Stony Brook, New York 11794, U.S.A. J. RAMS and C. G. GRAY Department of Physics, University of Guelph, Guelph, Ontario N1G 2W1, Canada and M. S. WERTHEIM Department of Mathematics, Hill Centre for the Mathematical Sciences, Rutgers University, Busch Campus, New Brunswick, New Jersey 08903, U.S.A. (Received 20June 1983 ; accepted 18 July 1983) In the first paper in this series, a new integral equation approximation --the Percus-Yevick excess (PYX) approximation--was introduced and contrasted with the corrected Percus-Yevick (CPY) approximation. In the present paper, these two approximations are evaluated for a molecular fluid, the molecules of which interact via a short ranged anisotropic potential prototypical of the overlap interactions encountered between homonuclear diatomic molecules. We conclude that the CPY and PYX approximations are quite accurate; interestingly, however, they unexpectedly yield very similar results leading us to conclude that, on balance, both approximations are equally good. 1. INTRODUCTION In the previous paper in this series [1, henceforth referred to as I ; equation (m.n) of I is denoted by (I.m.n)], the single super chain/linearized hypernetted chain (SSC/LHNC), single super chain f-expansion (SSCF) and Percus- Yevick (PY) approximations were evaluated for short ranged anisotropic potentials. This was done by solving these approximations numerically for the delta model (DM) fluid, the molecules of which interact via the DM potential given by (I.1.14). The molecular total correlation function h(12) (which depends on both orientations and positions of molecules 1 and 2 and is defined by (I.1.2)) predicted by the various approximations is resolved into a set of spherical harmonic coefficients h(11121 ; r) (see 3 and Appendix B of I) and then compared with Monte Carlo (MC) simulation results for the same quantities. The DM potential is prototypical of the short ranged anisotropic overlap t Present address : Department of Chemical Engineering, Thornton Hall, University of Virginia, Charlottesville, Virginia 22901, U.S.A. $ On leave from : Physics Department, Banaras Hindu University, Varanasi , India.
2 1134 P.T. Cummings et al. forces that characterize homonuclear diatomic molecules. On the basis of the results presented in I, it was concluded that SSC/LHNC is not a good approximation for short ranged anisotropic interactions. The SSCF approximation was found to be superior to SSC/LHNC, but nevertheless is unable to predict certain features in the MC results that are rendered qualitatively correctly by the PY approximation. In I, an additional new approximation--the Percus-Yevick excess (PYX) approximation--was introduced and contrasted to the corrected Percus-Yevick (CPY) approximation of Smith and Henderson [2]. The CPY and PYX approximation are expressed by approximate closure relations between h(12) and the molecular direct correlation function c(12) defined by the Ornstein- Zernike equation (I.1.3). If the molecular pair potential u(12) is divided into a reference part (ura(12)) and a perturbation part (Upert(12)), then the Mayer f-function f(12)=exp (-u(12)/kbt-1 (k B is Bohzmann's constant and T the absolute temperature) can be written as f(12) = fre~(12) + fp~rt(12) + f~t(12)fp~t(12). (1.1) Dividing h(12) and c(12) into reference and perturbation parts as well (i.e. h(12)=hret(12)+hpert(12) and similarly for c(12), where hre~(12 ) and Cref(12) are the total and direct correlation functions for a fluid interacting via the potential u~a(12)), the CPY and PYX approximations are respectively given by (see (I.A.17) and I.A.18)) Cpert(12 ) = fpert(12)[gref(12) + [1 + fre~(12)j[hpert(12) -- Cpert(12)] ] + f~(12)[hve,t(12)- %~t(12)], (1.2) %,n(12) = fp,rt(12)[gre,(12) +gra(12)[hpe~t(12)- %~t(12)]] +h,a(12)[hvc~t(12)-% rt(12)], (1.3) where g~ef(12)=h~(12)+l. In I it was suggested that renormalized PY approximations such as PYC and PYX (i.e. approximations that yield exact reference system quantities in the absence of the anisotropic perturbation) should yield superior results to PY by bringing the qualitative predictions of PY into better quantitative agreement with MC. In this paper, we report results for PYX and CPY that bear out this expectation. Perhaps a little surprisingly, the difference between CPY and PYX is not significant enough to definitively choose one as more accurate than the other; however, it will be clear from the results presented that both are superior overall to PY and SSCF. For the computations reported in 2, we choose a spherically symmetric reference system (so that qra(12)=qret(ra2) for q--u, f, h and c, rl~ being the distance between the centres of molecule 1 and 2). This simplifies the application of the closure (1.2) and (1.3) considerably. The spherically symmetric reference system chosen is the same as that used in 2 of I for SSCF : that is, the reference system Mayer f-function f~,,(rle ) is defined to be the unweighted angular average of f(12) as described by (I.2.10) and (I.2.11). Thus, in all the comparisons that follow, SSCF, PYX and CPY all share the same reference system, and the ability of each approximation to incorporate correctly the detailed effects of the anisotropic perturbation in the Mayer f-function is being evaluated.
3 CP Y and P YX approximations 1135 All the computational details for producing the results in the following section are to be found explicitly in I, or can be inferred from it. (The one noteworthy point is that products of quantities having orientational dependence in the CPY and PYX closures (1.2) and (1.3) are handled using the method applied in I to the similar product that arises in the PY approximation.) Thus, we do not discuss technical aspects of the computation here, and refer the interested reader to I for details. The state conditions employed in the paper are the same as those of I ; hence, the MC results used here are those described in detail in I. 2. RESULTS We shall consider g(r), the unweighted angular average of the molecular correlation function g(12)=h(12)+l, at the reduced density p*=poa=0.8, reduced temperature T*=kBT/ =0.719 and three values of 8, the strength of the anisotropic perturbation. (The parameters or, E and 8 are defined in (I.1.14).) In figure 1, g(r) for 8=0.1 is shown. At this small value of 8, CP, PYX and SSCF are all equally good, being indistinguishable on the scale of this figure ; the PY approximation is not good due to overestimation of the first peak, a point discussed in greater detail in I. For ~=0"25 (figure 2), we find that PYX and CPY are again indistinguishable and are very accurate, differing by at most 2 per cent from the MC results. The position of the principal peak in the SSCF g(r) is a little to the left of the MC results, as is the PY result. Thus, PYX and CPY are clearly superior at this value of 8. For 8 = 0.3 (figure 3), both SSCF and PY have incorrect peak positions and heights, while PYX and CPY are nearly indistinguishable and in good agreement with ] i I i g(r} 2 t i Z! l I = I r/o" 2 3 Figure 1. The angle average pair co-_relation function g(r) for the DM potential calculated at p*=0-8, T*=0.719 and 3=0.1 using the PY ( ), CPY (... ), PYX(... ) and SSCF(...) approximations. The circles are MC results. Note that on the scale of this graph CPY, PYX and SSCF are indistinguishable.
4 1136 P.T. Cummings et al. [ i I i 3 g(r) 2 0 ~ I, I 1 rio" 2 Figure 2. The angle average pair correlation function g(r) for the DM potential calculated at p*=0.8, T*=0.719 and The legend used is the same as that for figure 1. Note that on the scale of this graph CPY and PYX are indistinguishable. [ i ] i 2 - g(r), I r/o. 2 Figure 3. The angle average pair correlation function g(r) for the DM potential calculated at p*=0.8, T~=0.719 and ~=0-3. The legend used is the same as that for figure 1. MC with an error of at most 4 per cent. Clearly, however, even PYX and CPY are showing the initial signs of breakdown in accuracy. Turning now to the principal harmonic h(202; r) of h(12), we find in figure 4 (8--0.1) a mixed set of results. The CPY result obtains the first trough correctly, but overestimates the height of the first peak ; PYX performs
5 CP Y and P YX approximations r/d 2 3 Figure 4. The spherical harmonic coefficient h(202 ; r) for the DM potential calculated at p*=0-8, T*=0.719 and 8=0-1. The legend used is the same as that for figure h(202;r) -I I, I 1 r/o" 2 Figure 5. The spherical harmonic coefficient h(202 ; r) for the DM potential calculated at p*=0.8, T*=0.719 and 8=0.25. The legend used is the same as that for figure 1. better overall since it agrees with the first peak quite well and makes only a small underestimation of the depth of the first trough. Thus we would rate the approximations in descending order of accuracy as (PYX, CPY, SSCF, PY). For ~ =0-25, figure 5 shows that PYX is again superior, since PYX and CPY are equally good in predicting the first trough (both overestimating by about
6 1138 P.T. Cummings et al. i loi I i 0 I h(202;r) I : q C I z I p 1 r/o" 2 3 Figure 6. The spherical harmonic coefficient h(202 ; r) for the DM potential calculated at 0*=0-8, T*=0.719 and 8=0.3. The legend used is the same as that for figure h(220;r)! I ' I ' 10 J o 'o Lo I -10 a I i rio" 2 3 Figure 7. The spherical harmonic coefficient h(220 ; r) for the DM potential calculated at 0*=0.8, T*=0-719 and 8=0.25. The legend used is the same as that for figure. 15 per cent) while PYX is better for the first peak and beyond. Thus, the accuracy ordering is again (PYX, CPY, SSCF, PY). In figure 6, we display h(202 ; r) for 5=0.3, and here we find that CPY and PYX are overestimating the true result (although overall PYX is better) For this value of 5, SSCF appears to be best, so that the accuracy ordering is (SSCF, PYX, CPY, PY).
7 CP Y and P YX approximations 1139 We now consider two other harmonics (h(220 ; r) and h(224 ; r)) in detail for 3=0.25 only. In figure 7, we show the results for h(220 ; r). The SSCF result has the best magnitude for the first peak; however the position is not rendered well by SSCF and additionally the first trough is missed. CPY and PYX are roughly equivalent here, since they both have the peak and trough in their correct positions; in magnitudes, CPY is better than PYX for the first peak, while the reverse is true for the first trough. The PY approximation is qualitatively correct but the magnitudes are excessive. Thus we conclude that the accuracy ordering is (PYX = CPY, PY, SSCF). Although not shown, the 8=0-3 results yield the same accuracy ordering while at 3=0.1 we find (CPY, PYX, PY, SSCF). lfi I I ' 12 h(224;r) I, I, 1 r/d 2 Figure 8. The spherical harmonic coefficient h(224 ; r) for the DM potential calculated at p*=0-8, T*=0.719 and 8=0.25. The legend used is the same as that for figure 1. The h(224 ; r) harmonie for 3 = 0-25 is shown in figure 8. No approximation is very satisfactory for this harmonic coefficient and it is difficult to draw up an accuracy ordering for this harmonic coefficient; the same is true for 3=0-1 and 3=0.3, although we have not exhibited those results here. Also, we have not shown the h(222; r) harmonic coefficient. For 8=0.3, all the approximations are qualitative only and no accuracy ordering can be deduced. For 3=0.1, we find (CPY, PYX, PY, SSCF) and for 3=0-25 (CPY=PYX, PY, SSCF). 3. CONCLUSION As the results reported in 2 make clear, the CPY and PYX approximations represent definite improvements over PY and SSCF. For the centres distribution and the principal harmonic coefficient h(202 ; r) they are especially
8 1140 P.T. Cummings et al. good, while the accuracy deteriorates for the higher order harmonic coefficients considered. In choosing between CPY and PYX, we note that the differences between the approximations are small and neither one is consistently more accurate than the other when results are compared with MC. Consequently, we conclude that CPY and PYX are equally good. This is a little surprising in view of the somewhat richer structure of PYX that was noted in I ; however, as the example of the hypernetted chain and Percus-Yevick approximations for hard spheres illustrates [3], in liquid state physics approximations which in principle have a richer structure do not necessarily yield better results for some hamiltonian models. Overall, we have been encouraged by the success of the relatively simple CPY and PYX approximations; however, their utility for long-ranged interactions has yet to be determined and is presently under study. The authors gratefully acknowledge NSERC (Canada) for financial support, and the Computer Services Council of the University of Guelph for the allocation of computer resources. P.T.C. acknowledges the support of the National Science Foundation, under grant CHE for this research; M.S.W. acknowledges the support of the National Science Foundation under grants CHE and CHE REFERENCES [1] CUMMINCS, P. T., RAM, J., BARKER, R., GRAY, C. G., and WERTHEXM, M. S., 1983, Molec. Phys., 48, [2] SMITH, W. R., and HENDERSON, D., 1978, ~. chem. Phys., 69, 319. [3] See, for example, HANSEN, J.-P., and McDoNALD, I. R., 1976, Physics of Simple Fluids (Academic Press).
Analytic 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 informationInfluence of the Friedel oscillations on the structure of liquid rubidium
J. Phys. F: Met. Phys.. 12(1982)233-44 Printed in Great Britain Influence of the Friedel oscillations on the structure of liquid rubidium P T Cummingst and P A Egelstaff Department of Physics, University
More informationA New Uniform Phase Bridge Functional: Test and Its Application to Non-uniform Phase Fluid
Commun. Theor. Phys. (Beijing, China) 39 (2003) pp. 231 237 c International Academic Publishers Vol. 39, No. 2, February 15, 2003 A New Uniform Phase Bridge Functional: Test and Its Application to Non-uniform
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 informationUsing Fundamental Measure Theory to Treat the Correlation Function of the Inhomogeneous Hard-Sphere Fluid
Using Fundamental Measure Theory to Treat the Correlation Function of the Inhomogeneous Hard-Sphere Fluid Jeff B. Schulte, Patrick A. Kreitzberg, Chris V. Haglund, and David Roundy Department of Physics,
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 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 informationAuthor copy. Structure factors of binary liquid metal alloys within the square-well model. 1. Introduction. Central European Journal of Physics
Cent. Eur. J. Phys. 7(3) 2009 584-590 DOI: 10.2478/s11534-009-0064-2 Central European Journal of Physics Structure factors of binary liquid metal alloys within the square-well model Research Article Nikolay
More informationThree Semi-empirical Analytic Expressions for the Radial Distribution Function of Hard Spheres
Commun. Theor. Phys. (Beijing, China) 4 (2004) pp. 400 404 c International Academic Publishers Vol. 4, No. 3, March 5, 2004 Three Semi-empirical Analytic Expressions for the Radial Distribution Function
More informationCalculating thermodynamic properties from perturbation theory I. An analytic representation of square-well potential hard-sphere perturbation theory
Ž. Fluid Phase Equilibria 154 1999 1 1 Calculating thermodynamic properties from perturbation theory I. An analytic representation of square-well potential hard-sphere perturbation theory Bing-Jian Zhang
More informationMixtures, I. Hard Sphere Mixtures*
Proceedings of the Natioruil Academy of Scienccs Vol. 67, No. 4, pp. 1818-1823, December 1970 One- and Two-Fluid van der Waals Theories of Liquid Mixtures, I. Hard Sphere Mixtures* Douglas Henderson and
More informationMean spherical model-structure of liquid argon
Prami0a, Vol. 6, No 5, 1976, pp. 284-290. Printed in ndia. Mean spherical model-structure of liquid argon R V GOPALA RAO and T NAMMALVAR Department of Physical Chemistry, Jadavpur University, Calcutta
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 informationarxiv: v1 [cond-mat.stat-mech] 1 Oct 2009
Shell Structure of Confined Charges at Strong Coupling J. Wrighton, J. Dufty, M. Bonitz, 2 and H. Kählert 2 Department of Physics, University of Florida, Gainesville, FL 326 arxiv:9.76v [cond-mat.stat-mech]
More informationPhase Field Crystal (PFC) Model and Density Functional Theory (DFT) of Freezing
Phase Field Crystal (PFC) Model and Density Functional Theory (DFT) of Freezing Pyrite Project Meeting October 14 th 2010 Arvind Baskaran John Lowengrub Density functional Theory of Freezing [ Ramakrishnan
More informationPerturbation approach for equation of state for hard-sphere and Lennard Jones pure fluids
PRAMANA c Indian Academy of Sciences Vol. 76, No. 6 journal of June 2011 physics pp. 901 908 Perturbation approach for equation of state for hard-sphere and Lennard Jones pure fluids S B KHASARE and M
More informationA comparison between computer simulation and theoretical results for ionic solutions
MOLECULAR PHYSICS, 1987, VOL. 62, No. 2, 461-473 A comparison between computer simulation and theoretical results for ionic solutions by J. M. CAILLOL, D. LEVESQUE and J. J. WEIS Laboratoire de Physique
More informationVariational Approach to the Equilibrium Thermodynamic Properties of Simple Liquids. I*
THE JOURNAL OF CHEMICAL PHYSICS VOLUME 51, NUMBER 11 1 DECEMBER 1969 Variational Approach to the Equilibrium Thermodynamic Properties of Simple Liquids. I* G.Ali MANSOORi (1) AND Frank B. CANFIELD (2)
More informationUNIVERSITY OF CINCINNATI
UNIVERSITY OF CINCINNATI DATE: August 14, 2002 I, Manuel Valera, hereby submit this as part of the requirements for the degree of: DOCTORATE OF PHILOSOPHY (Ph.D.) in: Physics It is entitled: Density Functional
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 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 informationOn the local and nonlocal components of solvation thermodynamics and their relation to solvation shell models
JOURNAL OF CHEMICAL PHYSICS VOLUME 109, NUMBER 12 22 SEPTEMBER 1998 On the local and nonlocal components of solvation thermodynamics and their relation to solvation shell models Nobuyuki Matubayasi Institute
More informationOverlapping Associating Fluids with Directional Bonds in a Bulk and Near a Hard Wall: Monte Carlo Study
Journal of Statistical Physics, Vol. 100, Nos. 12, 2000 Overlapping Associating Fluids with Directional Bonds in a Bulk and Near a Hard Wall: Monte Carlo Study Douglas Henderson, 1 Ian Garcia, 2 Stefan
More informationProvided by the author(s) and NUI Galway in accordance with publisher policies. Please cite the published version when available. Title Structure of molecular liquids: cavity and bridge functions of the
More informationNote on the Perturbation Equation of State of Barker and Henderson*
fhe JOURNAL OF CHEMICAL PHYSICS VOLUME S1, NUMBER 12 15 DECEMBER 1969 Note on the Perturbation Equation of State of Barker Henderson* G.Ali MANSOORI, (1) Joe A. PROVINE, (2) AND Frank B. CANFIELD (3) School
More informationTHE JOURNAL OF CHEMICAL PHYSICS 126,
THE JOURNAL OF CHEMICAL PHYSICS 16 44503 007 Development of an equation of state for electrolyte solutions by combining the statistical associating fluid theory and the mean spherical approximation for
More informationDensity-functional theory for classical fluids and solids
Density-functional theory for classical fluids and solids C. Ebner Department of Physics, Ohio State University, Columbus, Ohio 43210 H. R. Krishnamurthy and Rahul Pandit Department of Physics, Indian
More informationA theory for calculating the number density distribution of. small particles on a flat wall from pressure between the two
theory for calculating the number density distribution of small particles on a flat wall from pressure between the two walls Kota Hashimoto a and Ken-ichi mano a a Department of Energy and Hydrocarbon
More informationThermodynamic Properties of Polymer Solutions
@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@?e?@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
More informationAuxiliary sites in the RISM approximation for molecular fluids
J. Phys. A: Math. Gen. 14 (1981) 1483-1512. Printed in Great Britain Auxiliary sites in the RISM approximation for molecular fluids P T Cummingsf, C G Gray and D E Sullivan Department of Physics, University
More informationThe inverse problem for simple classical liquids: a density functional approach
J. Phys.: Condens. Matter 9 (1997) L89 L98. Printed in the UK PII: S0953-8984(97)79859-X LETTER TO THE EDITOR The inverse problem for simple classical liquids: a density functional approach Yaakov Rosenfeld
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 informationScaled particle theory for hard sphere pairs. II. Numerical analysis
THE JOURNAL OF CHEMICAL PHYSICS 125, 204505 2006 Scaled particle theory for hard sphere pairs. II. Numerical analysis Swaroop Chatterjee and Pablo G. Debenedetti a Department of Chemical Engineering, Princeton
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 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 informationWeighted density functional theory of the solvophobic effect
PHYSICAL REVIEW E, VOLUME 64, 021512 Weighted density functional theory of the solvophobic effect Sean X. Sun Department of Chemistry, University of California, Berkeley, California 94720 Received 24 July
More informationExcess Entropy, Diffusion Coefficient, Viscosity Coefficient and Surface Tension of Liquid Simple Metals from Diffraction Data
Materials Transactions, Vol. 43, No. 1 (2002) pp. 67 to 72 c 2002 The Japan Institute of Metals EXPRESS REGULAR ARTICLE Excess Entropy, Diffusion Coefficient, Viscosity Coefficient and Surface Tension
More informationThe direct correlation functions and bridge functions for hard spheres near a large hard sphere
Title The direct correlation functions and bridge functions for hard spheres near a large hard sphere Author(s) Henderson, D; Chan, KY; Degrève, L Citation The Journal Of Chemical Physics, 1994, v. 101
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 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 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 informationarxiv:cond-mat/ v3 [cond-mat.stat-mech] 27 Mar 2001
Velocity correlations and the structure of nonequilibrium hard core fluids arxiv:cond-mat/8359v3 [cond-mat.stat-mech] 27 Mar 21 James F. Lutsko Physics Division, Starlab Rue Engelandstraat 555 B-118 Brussels,
More informationTreatment of Overlapping Divergences in the Photon Self-Energy Function*) Introduction
Supplement of the Progress of Theoretical Physics, Nos. 37 & 38, 1966 507 Treatment of Overlapping Divergences in the Photon Self-Energy Function*) R. L. MILLSt) and C. N. YANG State University of New
More informationStructure and phase behaviour of colloidal dispersions. Remco Tuinier
Structure and phase behaviour of colloidal dispersions Remco Tuinier Yesterday: Phase behaviour of fluids and colloidal dispersions Colloids are everywhere Hard sphere fluid at the base understanding fluids
More informationPair correlation function of short-ranged square-well fluids
THE JOURNAL OF CHEMICAL PHYSICS 122, 084510 2005 Pair correlation function of short-ranged square-well fluids J. Largo a and J. R. Solana b Departamento de Física Aplicada, Universidad de Cantabria, E-39005
More informationThe infrared properties of the energy spectrum in freely decaying isotropic turbulence
The infrared properties of the energy spectrum in freely decaying isotropic turbulence W.D. McComb and M.F. Linkmann SUPA, School of Physics and Astronomy, University of Edinburgh, UK arxiv:148.1287v1
More informationMonte Carlo Simulations for a Soft Sphere Fluid
Monte Carlo Simulations for a Soft Sphere Fluid Patrick Kreitzberg Advisor: Dr. David Roundy May 8 th, 2015 1 Contents 1 2 3 5 7 Abstract Introduction Methods Results and Discussion Conclusion Acknowledgments
More informationRotation in liquid 4 He: Lessons from a highly simplified model
JOURNAL OF CHEMICAL PHYSICS VOLUME 114, NUMBER 10 8 MARCH 2001 Rotation in liquid 4 He: Lessons from a highly simplified model Kevin K. Lehmann a) Department of Chemistry, Princeton University, Princeton,
More informationThe solution of the reference hypernetted-chain approximation for water-like models
MOLECULAR PHYSICS, 1988, VOL. 65, No. 5, 1105-1119 The solution of the reference hypernetted-chain approximation for water-like models by P. G. KUSALIKf and G. N. PATEY Department of Chemistry, University
More informationTheory of liquids and polymers Prof. Dr. Walter Schirmacher, WS 2010/11, Univ. Mainz
Theory of liquids and polymers Prof. Dr. Walter Schirmacher, WS /, Univ. Mainz Contents Introduction Structure of Liquids. Molecular distribution functions.................................. Scattering
More informationEquation of State of Dense Helium
Iowa State University From the SelectedWorks of Richard Alan Lesar October 31, 1988 Equation of State of Dense Helium Richard Alan Lesar, Los Alamos National Laboratory Available at: https://works.bepress.com/richard_lesar/27/
More informationSTRUCTURE OF LIQUID ALKALINE EARTH METALS AND METAL ALLOYS USING A NEW TRANSFERABLE LOCAL PSEUDOPOTENTIAL
Journal of Optoelectronics and Advanced Materials Vol. 5 No. 5 003 p. 181-191 STRUCTURE OF LIQUID ALKALINE EARTH METALS AND METAL ALLOYS USING A NEW TRANSFERABLE LOCAL PSEUDOPOTENTIAL H. Kes a S. S. Dalgic
More informationTheory of Interfacial Tension of Partially Miscible Liquids
Theory of Interfacial Tension of Partially Miscible Liquids M.-E. BOUDH-HIR and G.A. MANSOORI * University of Illinois at Chicago (M/C 063) Chicago, Illinois USA 60607-7052 Abstract The aim of this work
More informationPrepTest. October 1996 Logic Games Setups. Chapter Twenty: October 1996 Logic Games Setups
PrepTest October 1996 Logic Games Setups Chapter Twenty: October 1996 Logic Games Setups 305 PrepTest. October 1996 Game #1: 1. A 2. A 3. D 4. A 5. B This is an Advanced Linear: Unbalanced: Underfunded
More informationClose coupling results for inelastic collisions of NH3 and Ar. A stringent test of a spectroscopic potential
12 August 1994 ELSEVIER Chemical Physics Letters 226 ( 1994) 22-26 CHEMICAL PHYSICS LETTERS Close coupling results for inelastic collisions of NH3 and Ar. A stringent test of a spectroscopic potential
More informationCritical Properties of Isobaric Processes of Lennard-Jones Gases
Critical Properties of Isobaric Processes of Lennard-Jones Gases Akira Matsumoto Department of Material Sciences, College of Integrated Arts Sciences, Osaka Prefecture University, Sakai, Osaka, 599-8531,
More informationGeometry explains the large difference in the elastic properties of fcc and hcp crystals of hard spheres Sushko, N.; van der Schoot, P.P.A.M.
Geometry explains the large difference in the elastic properties of fcc and hcp crystals of hard spheres Sushko, N.; van der Schoot, P.P.A.M. Published in: Physical Review E DOI: 10.1103/PhysRevE.72.067104
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 informationTSTC Lectures: Theoretical & Computational Chemistry
TSTC Lectures: Theoretical & Computational Chemistry Rigoberto Hernandez May 2009, Lecture #2 Statistical Mechanics: Structure S = k B log W! Boltzmann's headstone at the Vienna Zentralfriedhof.!! Photo
More informationStructural properties of low-density liquid alkali metals
PRAMANA c Indian Academy of Sciences Vol. 65, No. 6 journal of December 2005 physics pp. 1085 1096 Structural properties of low-density liquid alkali metals A AKANDE 1, G A ADEBAYO 1,2 and O AKINLADE 2
More informationA Study of Integral Equations for Computing Radial Distribution Functions
Western Michigan University ScholarWorks at WMU Master's Theses Graduate College 12-1984 A Study of Integral Equations for Computing Radial Distribution Functions Robert C. Scherzer Western Michigan University
More informationDetermining Liquid Structure from the Tail of the Direct Correlation Function
Journal of Statistical Physics, Vol. 100, Nos. 12, 2000 Determining Liquid Structure from the Tail of the Direct Correlation Function Kirill Katsov 1 and John D. Weeks 1, 2 Received October 5, 1999; revised
More informationDEVIL CHAPTER 6 TEST REVIEW
IB PHYSICS Name: Period: Date: # Marks: 51 DEVIL PHYSICS BADDEST CLASS ON CAMPUS CHAPTER 6 TEST REVIEW 1. A cyclist rides around a circular track at a uniform speed. Which of the following correctly gives
More informationHowever, reliability analysis is not limited to calculation of the probability of failure.
Probabilistic Analysis probabilistic analysis methods, including the first and second-order reliability methods, Monte Carlo simulation, Importance sampling, Latin Hypercube sampling, and stochastic expansions
More informationMonte Carlo simulations of confined hard sphere fluids
Physics Procedia 3 (2010) 1493 1497 www.elsevier.com/locate/procedia Monte Carlo simulations of confined hard sphere fluids G.E. Schröder-Turk, S. Kuczera, D. Hörndlein, F. Krach, T. Hiester, K. Mecke
More informationStructural properties of low density liquid alkali metals
PRAMANA c Indian Academy of Sciences Vol. xx, No. x journal of xxxxxx 2005 physics pp. 1 12 Structural properties of low density liquid alkali metals A AKANDE 1, G A ADEBAYO 1,2 and O AKINLADE 2 1 The
More informationON ALGORITHMS FOR BROWNIAN DYNAMICS COMPUTER SIMULATIONS
COMPUTATIONAL METHODS IN SCIENCE AND TECHNOLOGY 4,35-42 (1998) ON ALGORITHMS FOR BROWNIAN DYNAMICS COMPUTER SIMULATIONS ARKADIUSZ C. BRAŃKA Institute of Molecular Physics, Polish Academy of Sciences, Smoluchowskiego
More informationSuggestions for Further Reading
Contents Preface viii 1 From Microscopic to Macroscopic Behavior 1 1.1 Introduction........................................ 1 1.2 Some Qualitative Observations............................. 2 1.3 Doing
More informationThe Liquid State ~~& R-E-S-O-N-A-N-C-E-I--Ju-n-e The Arrangement of Atoms.
The Liquid State 1. The Arrangement of Atoms K R Rao The liquid state of matter is of great practical importance. The arrangement of atoms in a liquid is more disordered than in a crystal, and can be studied
More informationProteins in solution: charge-tuning, cluster formation, liquid-liquid phase separation, and crystallization
HERCULES Specialized Course: Non-atomic resolution scattering in biology and soft matter Grenoble, September 14-19, 2014 Proteins in solution: charge-tuning, cluster formation, liquid-liquid phase separation,
More informationA THEORETICAL STUDY OF DILUTE AQUEOUS ELECTROLYTE SOLUTIONS PETER GERARD KUSALIK. B.Sc, The University of Lethbridge, 1981
A THEORETICAL STUDY OF DILUTE AQUEOUS ELECTROLYTE SOLUTIONS By PETER GERARD KUSALIK B.Sc, The University of Lethbridge, 1981 A THESIS SUBMITTED IN PARTIAL FULFILLMENT OF THE REQUIREMENTS FOR THE DEGREE
More informationDensity Functional Theory of the Interface between Solid and Superfluid Helium 4
Density Functional Theory of the Interface between Solid and Superfluid Helium 4 Frédéric Caupin and Tomoki Minoguchi Laboratoire de Physique Statistique de l Ecole Normale Supérieure associé aux Universités
More informationpk~ T g(r)= g g (r)rj". The exact coefficients up to first order are g, (r)=e(r 1), g, (r)=b(r 1)e(2 r)(8 6r+, 'r ), (2.5) where e is the Heaviside
PHYSICAL REVIEW A VOLUME 43, NUMBER 10 15 MAY 1991 Radial distribution function for hard spheres S. Bravo Yuste and A. Santos Departamento de Frsica, Universidad de Extremadura, 06071 Badajoz, Spain (Received
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 informationNON-EQUILIBRIUM MOLECULAR DYNAMICS CALCULATION OF THE
FluidPhaseEquilibria, 53 (1989) 191-198 Elsevier Science Publishers B.V., Amsterdam - Printed in The Netherlands 191 NON-EQUILIBRIUM MOLECULAR DYNAMICS CALCULATION OF THE TRANSPORT PROPERTIES OF CARBON
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 informationH = 1 2 τψ gψ4, (0.0.1)
32.1 Landau theory We have derived the Ginzburg-Landau Hamiltonian in Lecture 18. If we apply the mean field theory, that is, if we ignore fluctuations, then H = 1 2 τψ2 + 1 4 gψ4, (0.0.1) should be interpreted
More informationShear viscosity of liquid rubidium at the triple point
J. Phys. F: Met. Phys. 18 (1988) 1439-1447. Printed in the UK Shear viscosity of liquid rubidium at the triple point P T Cummingsi- and G P Morriss-l: i Department of Chemical Engineering, Thornton Hall,
More informationNUCLEATION IN FINITE SYSTEMS: THEORY AND COMPUTER SIMULATION*t. M. RAO and B. J. BERNE. Chemistry Dept, Columbia University, New York, U.S.A.
NUCLEATION IN FINITE SYSTEMS: THEORY AND COMPUTER SIMULATION*t M. RAO and B. J. BERNE Chemistry Dept, Columbia University, New York, U.S.A. Abstract. The free energy of formation of a droplet in a finite
More informationClose to close packing
Close to close packing Oliver Penrose Mathematics, Heriot-Watt University, Edinburgh EH14 4AS, UK and George Stell, Chemistry, SUNY, Stony Brook, NY, USA October 22, 2003 Journal of Statistical Physics,
More informationDisordered Hyperuniformity: Liquid-like Behaviour in Structural Solids, A New Phase of Matter?
Disordered Hyperuniformity: Liquid-like Behaviour in Structural Solids, A New Phase of Matter? Kabir Ramola Martin Fisher School of Physics, Brandeis University August 19, 2016 Kabir Ramola Disordered
More informationSTUDY OF INTEGRAL EQUATION THEORY AND HUMAN BETA DEFENSIN 28
STUDY OF INTEGRAL EQUATION THEORY AND HUMAN BETA DEFENSIN 28 CHEN WENWEN School of Chemical and Biomedical Engineering A thesis submitted to Nanyang Technological University in partial fulfillment of the
More informationHeadspace Raman Spectroscopy
ELECTRONICALLY REPRINTED FROM SEPTEMBER 2014 Molecular Spectroscopy Workbench Raman Spectroscopy We examine vapor-phase Raman spectroscopy through the acquisition of spectra from gas molecules confined
More informationRotational-translational relaxation effects in diatomic-gas flows
Rotational-translational relaxation effects in diatomic-gas flows V.V. Riabov Department of Computer Science, Rivier College, Nashua, New Hampshire 03060 USA 1 Introduction The problem of deriving the
More informationVapour-liquid equilibria for two centre Lennard-Jones diatomics and dipolar diatomics
MOLECULAR PHYSICS, 1993, VOL. 80, No. 4, 997-1007 Vapour-liquid equilibria for two centre Lennard-Jones diatomics and dipolar diatomics By GIRIJA S. DUBEYt, SEAMUS F. O'SHEA Department of Chemistry, University
More informationG : Statistical Mechanics
G5.651: Statistical Mechanics Notes for Lecture 8 I. DISTRIBUTION FUNCTIONS AND PERTURBATION THEORY A. General formulation Recall the expression for the configurational partition function: Z N = dr 1 dr
More informationarxiv: v1 [cond-mat.soft] 23 Nov 2018
Bulk structural informations from density functionals for patchy particles arxiv:8.9388v [cond-mat.soft] 3 Nov 8 Daniel Stopper,, Frank Hirschmann, Martin Oettel,, and Roland Roth Institute for Theoretical
More informationTwo Methods for Modelling the Propagation of Terahertz Radiation in a Layered Structure
Journal of Biological Physics 29: 141 148, 2003. 2003 Kluwer Academic Publishers. Printed in the Netherlands. 141 Two Methods for Modelling the Propagation of Terahertz Radiation in a Layered Structure
More informationarxiv:cond-mat/ v1 5 Aug 2002
Superfluidity in a Doped Helium Droplet E. W. Draeger and D. M. Ceperley Department of Physics and National Center for Supercomputing Applications, University of Illinois Urbana-Champaign, 68 Path Integral
More informationThe dielectric constant of polar hard dumb-bells A point charge model for chloromethane
MOLECULAR PHYSICS, 1982, VOL. 45, No. 5, 1099-1112 The dielectric constant of polar hard dumb-bells A point charge model for chloromethane by G. P. MORRISS School of Chemical Engineering, Cornell University,
More informationEXAM-3 PHYS 201 (Fall 2006), 10/31/06
EXAM-3 PHYS 201 (Fall 2006), 10/31/06 Name: Signature: Duration: 75 minutes Show all your work for full/partial credit! Include the correct units in your final answers for full credit! Unless otherwise
More informationStark effect of a rigid rotor
J. Phys. B: At. Mol. Phys. 17 (1984) 3535-3544. Printed in Great Britain Stark effect of a rigid rotor M Cohen, Tova Feldmann and S Kais Department of Physical Chemistry, The Hebrew University, Jerusalem
More informationPhotofragment angular momentum distributions in the molecular frame: Determination and interpretation
JOURNAL O CHEMICAL PHYSICS VOLUME 110, NUMBER 7 15 EBRUARY 1999 Photofragment angular momentum distributions in the molecular frame: Determination and interpretation T. Peter Rakitzis and Richard N. Zare
More informationPhysics 212: Statistical mechanics II Lecture XI
Physics 212: Statistical mechanics II Lecture XI The main result of the last lecture was a calculation of the averaged magnetization in mean-field theory in Fourier space when the spin at the origin is
More informationComputational Problem: Keplerian Orbits
Computational Problem: Keplerian Orbits April 10, 2006 1 Part 1 1.1 Problem For the case of an infinite central mass and an orbiting test mass, integrate a circular orbit and an eccentric orbit. Carry
More informationInteraction potentials for soft and hard ellipsoids
Interaction potentials for soft and hard ellipsoids R. Everaers and M.R. Ejtehadi Max-Planck-Institut für Polymerforschung, Postfach 348, D-550 Mainz, Germany January 7, 003 Using results from colloid
More informationDynamic force matching: Construction of dynamic coarse-grained models with realistic short time dynamics and accurate long time dynamics
for resubmission Dynamic force matching: Construction of dynamic coarse-grained models with realistic short time dynamics and accurate long time dynamics Aram Davtyan, 1 Gregory A. Voth, 1 2, a) and Hans
More informationMapping of charge polydispersity onto size polydispersity in colloidal suspensions
J. Phys.: Condens. Matter S (1991) 997-1002. Printed in the UK LETTER TO THE EDITOR Mapping of charge polydispersity onto size polydispersity in colloidal suspensions Hartmut Lowen, Jean-NoGI Row and Jean-Pierre
More informationPhase 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 informationBrownian Motion and Langevin Equations
1 Brownian Motion and Langevin Equations 1.1 Langevin Equation and the Fluctuation- Dissipation Theorem The theory of Brownian motion is perhaps the simplest approximate way to treat the dynamics of nonequilibrium
More information