arxiv: v1 [cond-mat.soft] 15 Nov 2011

Transcription

1 Wha consiues a simple liquid? Tond S. Ingebigsen, Thomas B. Schøde, and Jeppe C. Dye DNRF cene Glass and Time, IMFUFA, Depamen of Sciences, Roskilde Univesiy, Posbox 260, DK-4000 Roskilde, Denmak (Daed: Novembe 16, 2011) axiv: v1 [cond-ma.sof] 15 Nov 2011 Simple liquids ae adiionally defined as many-body sysems of classical paicles ineacing via adially symmeic pai poenials. We sugges ha a simple liquid should be defined insead by he popey of having song coelaions beween viial and poenial enegy equilibium flucuaions in he NV T ensemble. Thee is consideable ovelap beween he wo definiions, bu also some noable diffeences. In he new definiion, fo insance, simpliciy is no a popey of he inemolecula poenial, because a liquid is usually only songly coelaing in pa of is phase diagam. Moeove, no all aomic liquids ae simple and no all simple liquids ae aomic (i.e., wih adially symmeic pai poenials). The main pa of he pape moivaes he new definiion of a simple liquid by showing fom compue simulaions ha a liquid is songly coelaing if and only if is inemolecula ineacions may be ignoed beyond he fis coodinaion shell (FCS). This is demonsaed by NV T simulaions of sucue and dynamics of seveal model liquids wih a shifed-foces cuoff placed a he fis minimum of he adial disibuion funcion, delimiing he FCS. The liquids sudied ae: Invese powe-law sysems ( n pai poenials wih n = 4, 6, 12, 18), Lennad-Jones (LJ) models (he sandad LJ model, hee genealized Kob-Andesen binay LJ mixues, he Wahnsom binay LJ mixue), he Buckingham model, he Dzuguov model, he LJ Gaussian model, he Gaussian coe model, he Hansen-McDonald molen sal model, he Wahnsom OTP model, he asymmeic dumbbell model, and he SPC/E wae model. The final pa of he pape summaizes seveal popeies of songly coelaing liquids, emphasizing ha hese ae simple han liquids in geneal. Simple liquids as defined hee as hose wih song viial/poenialenegy coelaions may be chaaceized (1) chemically by he fac ha he liquid s popeies ae fully deemined by ineacions fom he molecules of he FCS, (2) physically by he fac ha hee ae isomophs in he phase diagam, i.e., cuves along which a numbe of popeies like excess enopy, sucue, and dynamics, ec ae invaian in educed unis, and (3) mahemaically by he fac ha houghou he phase diagam he educed-coodinae consan-poenial enegy hypesufaces define a one-paamee family of compac Riemannian manifolds. No igoous poof is given ha he chemical chaaceizaion follows fom he song coelaion popey, bu we show ha he FCS chaaceizaion is consisen wih he exisence of isomophs in he phase diagam of all songly coelaing liquids. I. INTRODUCTION Going back o Plao, classificaion o caegoizaion is he episemological pocess ha goups objecs based on simila popeies [1]. Having pimaily biological examples in mind, Aisole held ha caegoies ae discee eniies chaaceized by popeies shaed by hei membes [2]. Locke in 1690 disinguished beween he nominal and he eal essence of an objec [3]. The nominal essence comes fom expeience and epesens he objec s appeaance, wheeas he eal essence epesens he objec s deepe, consiuing feaues. Fo insance, he eal essence of a maeial hing is is aomic consiuion, because his is he causal basis of all he hing s obsevable popeies [4]. A scienific classificaion is paiculaly useful if i eflecs he eal essence of he objecs in quesion by idenifying hei undelying common feaues, fom which he moe obvious and easily obsevable nominal popeies follow. Having in mind Locke s concep of he eal essence, we ague below fo a new definiion of he class of simple liquids. The guiding quesion is: When is a liquid simple? The concep of a simple liquid is old and emains cenal, as evidenced by he 2003 book ile Basic conceps fo simple and complex liquids [5] o he eview of nonsimple liquids eniled Theoy of complicaed liquids fom 2007 [6]. Geneaions of liquid sae heoiss wee inoduced o his exciing opic by sudying Hansen and McDonald s classic exbook Theoy of simple liquids [7]. This book fis appeaed in 1976 following a decade of specacula pogess in he heoy of liquids, caalyzed by some of he vey fis scienific compue simulaions. Physiciss love simple sysems. This eflecs he fundamenal paadigm of physics ha in ode o capue a given phenomenon, simple is bee. Mos classificaions in physics ae clea cu, fo example he classificaion of elemenay paicles ino bayons and lepons, wheeas classificaions in ohe sciences usually have a wealh of bodeline cases. Due o he divesiy of molecules i is easonable o expec a definiion of simple liquids o be of he lae ype. In Ref. 7 a simple liquid is defined as a classical sysem of appoximaely spheical, nonpola molecules ineacing via pai poenials. This definiion has much in common wih he chemiss nonassociaed liquids [8], bu hee ae some significan diffeences. Chemiss geneally egad a liquid as simple even if i consiss of lage elongaed molecules, as long as hese ae wihou inenal degees of feedom and ineac puely via van de Waals

2 2 foces. Many physiciss would pobably disagee. Thus i is fa fom ivial o ask: Wha chaaceizes a simple liquid? Moe accuaely: Given a classical many-body sysem wih poenial enegy as a funcion of he bodies cenes of masses and hei spaial oienaions, is i possible o give a quaniaive cieion of how simple he sysem is? If yes, is simpliciy encoded uniquely in he poenial enegy funcion o can he degee of simpliciy vay houghou he phase diagam? Recen woks idenified and descibed he popeies of wha we have, fo lack of a bee name, emed songly coelaing liquids [9 20]. In hese liquids, he viial W and he poenial enegy U coelae songly in hei consan-volume hemal-equilibium flucuaions. Recall ha he aveage viial W gives he conibuion o he pessue fom inemolecula ineacions added o he ideal-gas em Nk B T deiving fom momenum anspo via paicle moion (N is he numbe of paicles, k B Bolzmann s consan, and T he empeaue): pv = Nk B T + W. (1) Hee p is he pessue and V he volume. The em songly coelaing liquid efes o he case when he W U coelaion coefficien in he NV T ensemble is lage han 0.9 [13]. If angula backes denoe an NV T ensemble aveage, he coelaion coefficien R is defined by R = W U ( W )2 ( U) 2. (2) An example of a songly coelaing liquid is he sandad Lennad-Jones liquid a ypical condensed-phase sae poins (even a zeo pessue), i.e., no fa fom he solid-liquid coexisence line. Many ohe sysems, including some simple molecula models, have been shown o be songly coelaing; we efe he eade o papes ha deive and documen he seveal simple popeies of songly coelaing liquids [9 19], ecenly biefly eviewed in Ref.20. These popeies ae summaized in Sec. IV A afe he pesenaion of he simulaion esuls. The pesen wok is moivaed by developmens iniiaed by ecen findings by Behie and Tajus [21, 22]. These auhos showed ha fo he viscous Kob-Andesen binay Lennad-Jones mixue [23, 24] he dynamics is no epoduced popely by cuing he poenials a hei minima accoding o he well-known Weeks-Chandle- Andesen (WCA) ecipe [25]. Subsequenly, he ole of he cuoff was sudied in wo papes showing ha placing a shifed-foces cuoff a he fis minimum of he pai coelaion funcions, gives good esuls fo Lennad- Jones ype sysems [26, 27]. This applies no only a modeae densiies, bu also a vey high densiies. Applying he same cuoff o wae does no wok popely [28]. Wae is an example of a non-songly coelaing liquid wih R 0 a ambien condiions (his is a consequence of wae s densiy maximum). These findings led us o wonde whehe i is a geneal popey of songly coelaing liquids ha he inemolecula ineacions may be ignoed beyond he FCS wihou compomising simulaion accuacy o any significan exen. The main pa of his pape shows ha, indeed, using such an FCS cuoff gives accuae simulaion esuls if he liquid is songly coelaing. The pape pesens esuls obained fom compue simulaions of 17 diffeen sysems, some of which ae songly coelaing and some no. We focus on he ole of he FCS in deemining liquid sucue and dynamics. Sucue is pobed by he adial disibuion funcion (RDF), dynamics by he incoheen inemediae scaeing funcion (ISF) a he waveveco coesponding o he fis peak of he RDF. The numeical evidence pesened is clea. By vaying he cuoff of he inemolecula foces we find ha in ode o ge accuae simulaion esuls i is enough o ake ino accoun he ineacions wihin he FCS if and only if he liquid is songly coelaing. In ohe wods, fo songly coelaing liquids ineacions beyond he FCS ae unimpoan, and his applies only fo hese liquids. This is an empiical finding fo which hee ae a pesen no compelling agumens. We ague biefly in he final pa of he pape, howeve, ha he finding is consisen wih known popeies of songly coelaing liquids. Secion II pesens he esuls fom molecula dynamics simulaions of he 17 diffeen sysems. Secion III summaizes he simulaion esuls. Secion IV gives an oveview of he many simple popeies songly coelaing liquids have, moivaing he suggesion ha a liquid is o be defined as simple wheneve i is songly coelaing a he sae poin in quesion. Secion V gives a few concluding emaks. II. MOLECULAR DYNAMICS SIMULATIONS OF ATOMIC AND MOLECULAR LIQUIDS In a compue simulaion he ineacions, which usually exend in pinciple o infiniy, ae uncaed a some disance c beyond which hey ae ignoed [26]. To avoid a disconinuiy in he foce, which can seveely affec he simulaion esuls [26, 29], he simulaions epoed below all use poenials uncaed such ha he foce goes coninuously o zeo a c. This is done by applying a so-called shifed foces (SF) cuoff [29] whee, if he pai poenial is v() and he pai foce is f() = v (), he oal shifed foce is given by f SF () = { f() f( c ) if < c, 0 if > c. (3) This coesponds o using he following pai poenial below c : v SF () = v() v ( c )( c ) v( c ). Using a SF cuoff gives moe accuae esuls and bee numeical sabiliy han using he sandad shifed-poenial

3 3 (SP) cuoff [26]. Ineesingly, his is so despie he fac ha a SF cuoff does no have he coec pai foce fo any, wheeas he foce is coec below c fo he SP cuoff. Appaenly, avoiding disconinuiy of he foce a c is moe impoan han mainaining he coec foce. I was ecenly discussed [14] why adding a linea em o he pai poenial affecs neihe sucue no dynamics o any significan exen. The eason is ha, when one neaes-neighbo disance deceases, ohes incease in such a way ha hei sum is viually consan. This agumen is exac in one dimension and holds o a good appoximaion in 3D consan-volume ensemble simulaions[14] (in consan-pessue ensemble simulaions he volume flucuaes and he agumen no longe applies). Coulomb ineacions have also been eaed by he SF cuoff pocedue. Alhough he Coulomb ineacion is long-anged and condiionally convegen, an SF cuoff gives esuls close o hose of he sandad, moe involved Ewald summaion mehod when c is sufficienly lage [30, 31]. All simulaions wee pefomed in he NV T ensemble wih peiodic bounday condiions and he Nose-Hoove algoihm [32, 33] using he RUMD [34] molecula dynamics package, which is opimized fo simulaions, in paicula, of small sysems (a few housand eniies). Fo he molecula models sudied below bond lenghs wee held fixed using he ime-symmeical cenal diffeence algoihm [35 37], which is a cenal-diffeence disceizaion of he consained equaions of moion peseving ime evesibiliy. The effec on sucue and dynamics of vaying he pai-poenial cuoff c was ecenly invesigaed fo he single-componen Lennad-Jones liquid (SCLJ) and he Kob-Andesen binay LJ (KABLJ) mixue [26, 27]. Fo boh sysems i was found ha if a SF cuoff is applied insead of he commonly used SP cuoff, he sandad cuoff c = 2.5σ can be deceased o σ while sill giving he coec physics. The value c = σ is close o he fis minimum of he RDF, implying ha all neaesneighbo FCS ineacions ae accouned fo. Deceasing he cuoff fuhe quickly affecs he simulaions esuls, howeve, an effec ha is ponounced fo he dynamics in he viscous egime [21, 22]. A. Fou invese-powe-law (IPL) fluids We conside fis sysems wih 100% coelaion beween viial and poenial enegy equilibium flucuaions in he NV T ensemble. I follows fom he definiion of he viial W = 1/3 i i U [29] ha he necessay and sufficien condiion fo W o coelae pefecly wih U is ha he poenial enegy is a so-called Eule homogeneous funcion of he paicle coodinaes i. This is clealy he case fo sysems wih invese powe-law (IPL) pai poenials (u() n ), bu noe ha poenials wih non-ivial angula dependence may also be Eule homogeneous. We simulaed single-componen IPL pai-poenial sysems wih exponens n = 18, 12, 6, 4 a densiy ρ = 5. Each sysem was simulaed a wo diffeen empeaues. The fou IPL poenials ae shown in Fig. 1 he simulaed sysems ange fom vey hash epulsive (n = 18) o quie sof and long anged (n = 4). Each sysem was simulaed wih fou diffeen SF cuoffs. The ole of he cuoff is invesigaed by choosing hee diffeen cuoffs: one placed a he fis minimum he RDF (ed colo), one coesponding o he half heigh of he RDF fom is minimum o is maximum (blue colo), and one placed o he igh of he RDF fis minimum displaced in disance o he igh he same amoun as he diffeence beween he fis and he second cuoff (geen colo). u() IPL, n = 18 IPL, n = 12 IPL, n = 6 IPL, n = 4 IPL, n = 3 Invese powe-law: u() = ε (σ/) n 5.0 FIG. 1: Invese powe-law (IPL) pai poenials wih exponens n = 18, 12, 6, 4 fom lef o igh, especively: black, ed, geen, blue. In he following we invesigae fo seveal sysems whehe i is possible o choose an FCS cuoff, i.e., a cuoff a he fis minimum of he pai coelaion funcion, and sill ge he coec physics. We sa by sudying seveal songly coelaing aomic liquids. Then daa ae pesened fo some aomic liquids ha ae no songly coelaing. Finally daa ae given fo wo molecula liquids ha ae songly coelaing as well as a wae model (ha is no). Deails of he models sudied, numbe of paicles, ec, ae given in Appendix A. The RDFs g() ae shown in Fig. 2. The simulaions wih a SF cuoff a he fis minimum of he RDF efeed o as FCS-cuoff simulaions give a faihful epesenaion of he sucue. The inses quanify as funcions of he cuoff he deviaions in RDF beween esuls fo an FCS cuoff and ue esuls by inegaing he numeical diffeence in pai coelaion funcion. Clealy, deviaions incease shaply when he cuoff enes he FCS (blue cosses).

4 4 g() 5.0 T = 0.7: c = T = 0.7: c = 50 T = 0.3: c = T = 0.3: c = 59 g () - g() d T = 0.7 T = c IPL 18 (ρ = 5) c Figue 3 shows he incoheen inemediae scaeing funcions (ISFs) evaluaed a he waveveco coesponding o he fis maximum of he saic sucue faco fo he low-empeaue sae poins of each of he fou IPL sysems. A good epesenaion of he dynamics is obained fo all IPL sysems when he FCS cuoff is used. g() g() g() = 5.300) = 0 = 59) = T = : c = T = : c = 48 T = : c = T = : c = 59 g () - g() d T = T = c IPL 12 (ρ = 5) = 5.300) = 0 = 59) = T = 0.70: c = T = 0.70: c = 65 T = 0.15: c = T = 0.15: c = 70 g () - g() d T = 0.70 T = c IPL 6 (ρ = 5) R( c = 5.300) = 0 (c) = 70) = T = 0.7: c = T = 0.7: c = 59 T = 0.1: c = T = 0.1: c = 80 g () - g() d c c 0 T = 0.7 T = c IPL 4 (ρ = 5) = 5.300) = 0 (d) = 80) = FIG. 2: RDFs fo single-componen IPL fluids wih exponens n = 18, 12, 6, 4, each simulaed a wo empeaues a densiy ρ = 5. The black and oange cuves show efeence simulaion esuls wih lage cuoffs epesening he ue IPL behavio, he ed and geen dos give esuls fom simulaions wih an FCS cuoff (maked by he veical ed dashed lines). The inses quanify he deviaions in he RDF fom he efeence RDF as funcions of he cuoff; a boh sae poins deviaions incease damaically when he cuoff enes he FCS (blue cosses). In each subfigue he viial/poenial-enegy coelaion coefficien R is given fo he simulaed cuoff (his quaniy is exacly uniy fo IPL sysems wih infinie cuoff) ( q = 6.64, ) ( q = 6.72, ) ( q = 6.72, ) IPL 18 (ρ = 5, T = 0.30) c = c = c = 59 (min) c = = 5.300) = 0 = 59) = IPL 12 (ρ = 5, T = 0) c = c = c = 59 (min) c = = 5.300) = 0 = 59) = IPL 6 (ρ = 5, T = 0.15) (c) c = c = c = 70 (min) c = = 5.300) = 0 = 70) =

5 5 ( q = 6.61, ) IPL 4 (ρ = 5, T = 0.10) (d) c = c = c = 80 (min) c = = 5.000) = 0 = 80) = FIG. 3: Incoheen inemediae scaeing funcions (ISFs) fo he IPL fluids a he lowes empeaue sae poins of Fig. 2. The black cuves give esuls fo a lage cuoff, he ed cosses fo FCS cuoffs (maked by he veical ed dashed lines in Fig. 2). n = 18, T = 0.30; n = 12, T = 0; (c) n = 6, T = 0.15; (d) n = 4, T = u() - AA Genealized (12, 04)-KABLJ AA Genealized (12, 06)-KABLJ AA Genealized (12, 10)-KABLJ Genealized Lennad-Jones: u() = ε/(m-n) [ n (σ/) m - m (σ/) n ] FIG. 4: The AA-paicle genealized Kob-Andesen (KABLJ) pai poenials wih fixed epulsive exponen 12 and hee diffeen aacive exponens n = 4, 6, 10. The model paamees ae: ɛ AA = 0, ɛ AB = 0, and ɛ BB = 0, σ AA = 1.12, σ AB = 0.90, σ BB = 0.99, m A = m B = 1. B. Five Lennad-Jones liquids Nex, we conside he mos sudied poenial in he hisoy of compue simulaions, he Lennad-Jones (LJ) pai poenial, [( σ ) 12 ( σ ) 6 ] u LJ () = 4ɛ. (4) Hee σ and ɛ, define, especively, he lengh and enegy scale (dimensionless unis defined by σ = ɛ = 1 ae used below). This poenial does no have 100% viial/poenial-enegy coelaions, bu sill has quie song coelaions wih coelaion coefficiens R > 0.9 in he condensed-fluid pa of he phase diagam (and, ineesingly, also in he cysalline phase [14]). We sudied he single-componen LJ (SCLJ) liquid, hee genealized 80/20 Kob-Andesen binay LJ (KABLJ) mixues wih epulsive exponen 12 and aacive exponens n = 4, 6, 10, and finally he Wahnsom binay LJ (WABLJ) liquid (Fig. 4 and Appendix A gives model deails). The influence of an SF cuoff on simulaion accuacy was invesigaed ecenly by Toxvaed and Dye [26, 27] fo he SCLJ liquid and he sandad KABLJ mixue (n = 6), bu fo compleeness we include esuls fo hese sysems hee as well. The ole of he cuoff is invesigaed by choosing hee diffeen cuoffs: one placed a he fis minimum he RDF (ed colo in Figs. 5-9), one coesponding o he half heigh of he RDF fom is minimum o is maximum (blue colo in Figs. 5-9), and one placed o he igh of he RDF fis minimum displaced in disance o he igh he same amoun as he diffeence beween he fis and he second cuoff (geen colo in Figs. 5-9). In Fig. 5 RDFs ae shown fo he SCLJ liquid a hee diffeen sae poins, anging fom high o low densiies. The ed coloed cicles show esuls fom simulaions wih an FCS cuoff (maked by he veical ed dashed line), he black cuves show he coesponding simulaions wih a lage cuoff (efeence sysem). The inses quanify he deviaions in he simulaed RDF fom he efeence RDF as a funcion of he cuoff. The efeence RDF of Fig. 5 is clealy epesened well using an FCS cuoff, while choosing he cuoff inside he FCS esuls in significan deeioaion of he ageemen. A low densiies [Figs. 5, and (c)], deviaions occu beween FCS cuoff simulaions and he efeence sysem. As menioned, he SCLJ liquid is songly coelaing in lage pas of is phase diagam, bu as densiy is loweed, he coelaions decease gadually and he liquid is no longe songly coelaing a he wo sae poins and (c), which have coelaion coefficiens R = 9 and R = 0, especively. These simulaions sugges ha only when a liquid is songly coelaing is i possible o ignoe ineacions beyond he FCS.

6 6 g() c = c = 38 g () - g() d c SCLJ (ρ = 5, T = 0.70) = 5.000) = 0.96 = 38) = ( q = 6.788, ) SCLJ (ρ = 5, T = 0.70) c = c = c = 38 (min) c = = 5.000) = 0.96 = 38) = g() 2.5 c = c = g () - g() d c SCLJ (ρ = 0, T = 0) = 6.000) = 9 = 1.663) = ( q = 6.38, ) SCLJ (ρ = 0, T = 0) c = c = 67 c = (min) c = = 6.000) = 9 = 1.663) = g() 2.5 c = c = g () - g() d c SCLJ (ρ = 5, T = 1.13) = 6.000) = 0 = 1.645) = 9 (c) 5.0 FIG. 5: RDFs fo he single-componen Lennad-Jones (SCLJ) liquid a hee diffeen sae poins: ρ = 5, T = 0.70 (R = 0.96); ρ = 0, T = 0 (R = 9); (c) ρ = 5, T = 1.13 (R = 0). The black cuves show efeence simulaions wih lage cuoffs, he ed dos show esuls fom simulaions wih an FCS cuoff (maked by he veical ed dashed lines). The inses quanify he deviaion in RDF fom he efeence RDF as funcions of he cuoff. A all hee sae poins deviaions incease significanly when he cuoff enes he FCS (blue cosses in he inses). Only fo sae poin, which is songly coelaing (i.e., R > 0.9), does an FCS cuoff lead o accuae esuls. Nex, we invesigae he SCLJ dynamics a he hee sae poins of Fig. 5. The dynamics is sudied via he incoheen inemediae scaeing funcion (ISF). The ISFs ae shown in Fig. 6; a all sae poins he dynamics is epesened well using an FCS cuoff. ( q = 6.42, ) SCLJ (ρ = 5, T = 1.13) (c) c = c = 40 c = (min) c = = 6.000) = 0 = 1.645) = FIG. 6: ISFs fo he SCLJ liquid a he sae poins of Fig. 5. The black cuves give he efeence cuoff esuls, he ed cuves he FCS cuoff esuls, he blue cuves esuls fo a cuoff a he half-heigh owads he fis maximum of he RDF, he geen cuves esuls fo a cuoff o he igh of he minimum a he same diffeence in cuoff disance as beween he ed and he blue cuves. ρ = 5, T = 0.70 (R = 0.96); ρ = 0, T = 1.15 (R = 9); (c) ρ = 5, T = 1.13 (R = 0). We poceed o invesigae mixues of wo diffeen (A and B) paicles ineacing wih LJ ype poenials. The cuoff used fo all hee ineacions (AA, AB, BB) ae placed a he same disance given by he AA-paicle ineacion. In Figs. 7 efeence and FCS cuoff esuls ae shown fo he AA-paicle RDFs of genealized KABLJ mixues wih epulsive exponen 12 and aacive exponens n = 4, 6, 10. Fo all invesigaed sae poins an FCS cuoff gives accuae esuls.

7 7 g AA () 5.0 c = c = g () - g() d c (12, 4)-KABLJ (ρ = 1.20, T = 0) = 4.700) = 4 = 1.433) = 0.95 A ( q = 7.23, ) (12, 4)-KABLJ (ρ = 1.20, T = 0) c = c = c = (min) c = = 4.700) = 4 = 1.433) = g AA () g AA () 5.0 c = c = g () - g() d c (12, 6)-KABLJ (ρ = 1.20, T = 0) = 4.700) = 0.95 = 1.429) = c = c = g () - g() d c (12, 10)-KABLJ (ρ = 1.20, T = 0) = 4.700) = 0.97 = 1.426) = 0.96 FIG. 7: RDFs fo genealized KABLJ mixues wih epulsive exponen 12 and diffeen aacive exponens n. The black cuves give he efeence cuoff esuls, he ed cuves he FCS cuoff esuls. The inses quanify he deviaion in RDF fom he efeence RDF as funcions of he cuoff. n = 4, ρ = 1.20, T = 0 (R = 4); n = 6, ρ = 1.20, T = 0 (R = 0.95); (c) n = 10, ρ = 1.20, T = 0 (R = 0.97). The A-paicle ISFs fo he sae poins of Fig. 7 ae shown in Fig. 8. Fo he KABLJ mixue also he dynamics placing he cuoff inside he FCS (blue cuves) fails o epoduce he dynamics popely, wheeas he dynamics is well appoximaed using an FCS cuoff (ed). Sligh deviaions ae noed fo he ed cuves, an issue consideed in Appendix B ha discusses alenaives fo delimiing he FCS. Simila esuls apply fo he B paicles (esuls no shown). (c) A ( q = 7.25, ) A ( q = 7.21, ) (12, 6)-KABLJ (ρ = 1.20, T = 0) c = c = c = (min) c = = 4.700) = 0.95 = 1.429) = (c) (12, 10)-KABLJ (ρ = 1.20, T = 0) c = c = c = (min) c = = 4.700) = 0.97 = 1.426) = FIG. 8: ISFs fo genealized KABLJ mixues wih epulsive exponen 12 and diffeen aacive exponens n. The ed and black cuves give, especively, esuls fo FCS cuoffs and lage efeence cuoffs. n = 4, ρ = 1.20, T = 0 (R = 4); n = 6, ρ = 1.20, T = 0 (R = 0.95); (c) n = 10, ρ = 1.20, T = 0 (R = 0.97). We close he invesigaion of aomic LJ liquids by showing a few esuls fo he Wahnsom-BLJ (WABLJ) mixue [38]. The lage (A)-paicle RDF and ISF ae shown, especively, in Figs. 9 and. Fo his songly coelaing LJ sysem sucue and dynamics ae also well epoduced using an FCS cuoff.

8 8 g AA () A ( q = 7.5, ) 5.0 c = c = g () - g() d c Wahnsöm-BLJ (ρ = 1.30, T = 0) = 4.500) = 0.99 = 1.427) = 0.98 Wahnsöm-BLJ (ρ = 1.30, T = 0) c = c = c = (min) c = = 4.500) = 0.99 = 1.427) = FIG. 9: The effec on sucue and dynamics of vaying he cuoff fo he Wahnsom (WABLJ) mixue. The ed and black cuves give, especively, esuls fo an FCS cuoff and a lage efeence cuoff. AA-paicle RDF a ρ = 1.30 and T = 0 (R = 0.99). The inse quanifies he deviaion in RDF fom he efeence RDF as a funcion of he cuoff. The A-paicle ISF a he same sae poin. In summay, fo all LJ ype sysems we find ha when hee ae song viial poenial-enegy coelaions (R > 0.9), an FCS cuoff gives accuae esuls fo boh sucue and dynamics. In he pas of he phase diagam whee coelaions ae weake, FCS cuoffs give pooe esuls. C. Buckingham liquid Nex, we conside he single-componen Buckingham liquid (SCB). The Buckingham poenial [39, 40] is simila o he LJ poenial, bu does no have an IPL epulsive em; insead he poenial s sho-disance behavio follows a seep exponenial. Consequenly, he Buckingham poenial does no divege a = 0. The paamees of he Buckingham poenial (Appendix A) wee chosen such ha he LJ poenial is well appoximaed in he lowe-enegy epulsive pa of he poenial (Fig. 10). u() - u() Buckingham: (12, 6)-Lennad Jones Buckingham u() = ε( 6/(α-6) exp[α(1 - / m )] - α/(α - 6) ( m /) 6 ) 5.0 FIG. 10: The Buckingham pai poenial (ed) and he LJ pai poenial (black). The paamees of he Buckingham poenial wee chosen such ha he LJ poenial is well appoximaed in he epulsive egion. Figues 11 and show, especively, he RDF and ISF fo he SCB liquid. g() ( q = 7.10, ) c = c = g () - g() d c Single-componen Buckingham (ρ = 0, T = 0) = 4.900) = 0.99 = 1.475) = c = c = c = (min) c = Single-componen Buckingham (ρ = 0, T = 0) = 4.900) = 0.99 = 1.475) = FIG. 11: The effec on sucue and dynamics of vaying he cuoff fo he single-componen Buckingham (SCB) liquid. The ed and black cuves give, especively, esuls fo an FCS cuoff and a lage efeence cuoff. RDF a ρ = 0 and T = 0 (R = 0.99). The inse quanifies he deviaion in RDF fom he efeence RDF as a funcion of he cuoff. ISF a he same sae poin. D. Dzuguov liquid This secion consides he Dzuguov (DZ) liquid [41]. Figue 12 shows he DZ pai poenial, which was oig-

9 9 inally suggesed as a model poenial fo which cysallizaion is impeded by enegeically punishing paicle sepaaions coesponding o he nex-neaes neighbo disance of an face-ceneed cubic laice. A sho disances he DZ pai poenial appoximaes he LJ poenial. ( q = 6.74, ) Dzuguov (ρ = 5, T = 0.75) c = c = (min) c = = 1.940) = 0.71 = 1.602) = 0.94 (12, 6)-Lennad-Jones (shifed) Dzuguov u() Dzuguov: u() = A( -n - B) exp[c / ( - a )] + B exp[d / ( - b)] FIG. 13: The effec on sucue and dynamics of vaying he cuoff fo he Dzuguov (DZ) liquid. The ed and black cuves give, especively, esuls fo an FCS cuoff and a lage efeence cuoff. RDF a ρ = 5 and T = 0.75 (R = 0.71). The inse quanifies he deviaion in RDF fom he efeence RDF as a funcion of he cuoff. ISF a he same sae poin. - FIG. 12: The Dzuguov (DZ) pai poenial [41] (blue). Also shown is he LJ pai poenial (black cuve). The DZ poenial appoximaes he LJ poenial aound he fis minimum, bu has a maximum a lage disances. E. Lennad-Jones Gaussian liquid The Lennad-Jones Gaussian liquid [42] (LJG) is a non-songly coelaing liquid wih he wo-minimum pai poenial shown in Fig. 14. The paamees of LJG model (Appendix A) ae such ha he LJG poenial minimum does no coincide wih ha of he SCLJ sysem [43]. Figues 13 and show, especively, he RDF and ISF of he DZ sysem. Fo his sysem he use of an FCS cuoff leads o poo esuls. This is no supising given he fac ha using an FCS cuoff emoves he maximum of he DZ poenial. Howeve, we noe ha he poo esuls of using an FCS cuoff is consisen wih he faily weak viial/poenial-enegy coelaions (R = 0.71), which suggess sudying ohe non-songly coelaing liquids. u() - - Lennad-Jones Gaussian Lennad-Jones Gaussian: u() = u 12,6 + εε 0 exp[-(- 0 ) 2 /2σ 0 2 ] g() c = c = g () - g() d c Dzuguov (ρ = 5, T = 0.75) = 1.940) = 0.71 = 1.602) = FIG. 14: The Lennad-Jones Gaussian (LJG) pai poenial [42]. Two disinc minima ae pesen. Resuls fom simulaing sucue and dynamics of he LJG liquid ae shown in Figs. 15 and, especively. The FCS cuoff does no give he coec RDF. Deviaions in he dynamics ae faily small, likely due o he fac ha he dynamics is fas (compae, e.g., he ime scale of he decay o ha of he DZ liquid in Fig. 13). Noe ha hee he FCS cuoff also emoves he second minimum, as was he case fo he DZ liquid.

10 10 g() 2.5 c = c = g () - g() d c Lennad-Jones Gaussian (ρ = 5, T = 0) = 5.000) = 4 = 1.242) = u() Gaussian coe model Gaussian coe model: u() = ε exp[-(/σ) 2 ] 2.5 FIG. 16: The Gaussian coe model (GCM) pai poenial. The foce a = 0 is zeo, implying ha a high densiies he paicles pefe o be siuaed on op of each ohe. This feaue, combined wih he model s geneally sof epulsion, esuls in seveal iniguing phenomena [46]. ( q = 8.7, ) 1 0 = 5.000) = 4 = 1.242) = 0.97 c = c = c = (min) c = 38 Lennad-Jones Gaussian (ρ = 5, T = 0) FIG. 15: The effec on sucue and dynamics of he cuoff fo he LJG liquid. The ed and black cuves give, especively, esuls fo an FCS cuoff and a lage efeence cuoff. RDF a ρ = 5 and T = 0 (R = 4). The inse quanifies he deviaion in RDF fom he efeence RDF (black cuve) as a funcion of he cuoff. ISFs a he same sae poin. Figue 17 shows he RDF and ISF fo he GCM liquid. The GCM cysallizes when an FCS cuoff is used. Fo his eason, obviously, an FCS cuoff is no able o epoduce sucue and dynamics of he efeence sysem. Noe, howeve, ha cysallizaion does no occu when he cuoff disance is chosen in he neighbohood of he FCS cuoff (see he inse). g() c = c = Cysal g () - g() d c Gaussian coe (ρ = 5, T = 1) = 5.000) = = 1.679) = Cysal F. Gaussian coe model The Gaussian coe model (GCM) [44, 45], which is no songly coelaing, is chaaceized by having a finie pai-poenial enegy a zeo sepaaion (Fig. 16). The high-densiy egime of he GCM model (ρ > ) has ecenly eceived aenion as a single-componen model glass fome [46] because i is no pone o cysallizaion and shows he chaaceisic feaues of glass-foming liquids (lage viscosiy, wo-sep elaxaion, ec). ( q = 6.25, ) c = c = 17 c = (min) c = Gaussian coe model (ρ = 5, T = 1) = 5.000) = = 1.679) = FIG. 17: The effec on sucue and dynamics of vaying he cuoff fo he GCM liquid. The ed and black cuves give, especively, esuls fo an FCS cuoff and a lage efeence cuoff. RDF a ρ = 5 and T = 1 (R = 0.77). The inse quanifies he deviaion in RDF fom he efeence RDF as a funcion of he cuoff. The ed cuve epesens a cysallized sae. ISF a he same sae poins.

11 11 G. The Hansen-McDonald molen sal model The final aomic sysem sudied is he so-called singlychaged molen sal model poposed by Hansen and Mc- Donald [47]. In Fig. 18 we see ha he sucue is no epesened well by he use of an FCS cuoff. Ineesingly, he dynamics is well epoduced using his cuoff even bee, in fac, han fo a lage cuoff (Fig. 18, geen cuve). g AA () A ( q = 4.10, ) c = c = g () - g() d c A singly-chaged molen sal (ρ = 0.37, T = 18) = 9.500) = 0.15 = 2.602) = A singly-chaged molen sal (ρ = 0.37, T = 18) c = c = 55 c = (min) c = = 9.500) = 0.15 = 2.602) = FIG. 18: The effec on sucue and dynamics of vaying he cuoff fo Hansen-McDonald singly-chaged molen sal model. Because of he compeing ineacions (Coulomb and n = 9 epulsive IPL) his model is no songly coelaing. The ed and black cuves give, especively, esuls fo an FCS cuoff and a lage efeence cuoff. AA-paicle RDF a ρ = 0.37 and T = 18 (R = 0.15). The inse quanifies he deviaion in RDF fom he efeence RDF as a funcion of he cuoff. A-paicle ISF a he same sae poin. Figues 19 and show he RDF and ISF of he OTP model. Boh quaniies ae well appoximaed using an FCS cuoff, alhough sligh deviaions ae noed fo he ISF (ed cuve, see Appendix B fo consideaions concening eliminaion of his behavio). The OTP model is a bode-line songly coelaing liquid (R = 0.91). g() ( q = 7.06, ) 5.0 c = c = g () - g() d c Lewis-Wahnsöm OTP (ρ = 0.33, T = 0.70) = 4.500) = 0.91 = 1.469) = 0.91 Lewis-Wahnsöm OTP (ρ = 0.33, T = 0.70) c = c = c = (min) c = = 4.500) = 0.91 = 1.469) = FIG. 19: The effec on sucue and dynamics of vaying he cuoff fo he Wahnsom OTP model. The ed and black cuves give, especively, esuls fo an FCS cuoff and a lage efeence cuoff. RDF a ρ = 0.33 and T = 0.70 (R = 0.91). The inse quanifies he deviaion in RDF fom he efeence RDF as a funcion of he cuoff. The spikes deive fom he bonds. ISFs a he same sae poin. H. Two songly coelaing molecula model liquids We finish he pesenaion of he numeical esuls by giving daa fo hee molecula model liquids. In his subsecion daa ae given fo wo songly coelaing molecula liquid models, he Lewis-Wahnsom OTP model[48, 49] and he asymmeic dumbbell model[19], models which epesen a molecule by hee and wo igidly bonded LJ sphees, especively. The nex subsecion gives daa fo a igid wae model. Figues 20 and show he coesponding figues fo he lage (A) paicle of he asymmeic dumbbell model a a viscous sae poin. The use of an FCS cuoff gives accuae esuls fo boh sucue and dynamics. The FCS cuoff was hee placed a he second minimum of he AA RDF, because he AA RDF has hee a lowe value han a he fis minimum. If he cuoff is placed a he fis minimum, clea deviaions ae noed (daa no shown).

12 12 g AA () 5.0 c = 00 c = 56 g () - g() d c The asymmeic dumbbell (ρ = 0.93, T = 6) = 00) = 0.96 = 56) = 0.96 O ( q = 9.23, ) = 4.980) = 8 = 65) = 8 Rigid SPC/E Wae (ρ = 0, T = 0) c = c = c = 65 (min) c = A ( q = 7.20, ) The asymmeic dumbbell (ρ = 0.93, T = 6) c = 00 c = c = 56 (min) c = = 00) = 0.96 = 56) = FIG. 20: The effec on sucue and dynamics of vaying he cuoff fo he asymmeic dumbbell model. The ed and black cuves give, especively, esuls fo an FCS cuoff and a lage efeence cuoff. Noe ha in his case he FCS cuoff is defined by using he second minimum (because he fis minimum is no he absolue minimum). RDF a ρ = 0.93 and T = 6 (R = 0.96). The inse quanifies he deviaion in RDF fom he efeence RDF (black cuve) as a funcion of he cuoff. A-paicle ISFs a he same sae poin. Fo he lowes cuoff he sysem cysallized. I. Rigid SPC/E wae We conside finally he igid SPC/E wae model [50] (Fig. 21). This model is no songly coelaing a ambien condiions (a fac eflecing wae s well-known densiy maximum [14]). The sucue of he SPC/E wae model is no well epesened using an FCS cuoff. The FCS cuoff dynamics on he ohe hand shows only sligh deviaions fom ha of he efeence cuve (black). g OO () c = c = 65 g () - g() d Rigid SPC/E Wae c (ρ = 0, T = 0) = 4.980) = 8 = 65) = 8 FIG. 21: The effec on sucue and dynamics of vaying he cuoff fo he igid SPC/E wae model [50]. The ed and black cuves give, especively, esuls fo an FCS and a lage efeence cuoff. Oxygen-oxygen RDF a ρ = 0 and T = 0 (R = 8). The inse quanifies he deviaion in RDF fom he efeence RDF as a funcion of he cuoff. Oxygen ISF a he same sae poin. III. SUMMARIZING THE SIMULATION RESULTS The pevious secion showed ha sucue and dynamics ae well appoximaed in simulaions using an FCS cuoff fo he following aomic and molecula sysems: The invese powe-law (IPL) sysems (n = 18, 12, 6, 4), he single-componen Lennad-Jones (SCLJ) liquid a densiy ρ = 5, he genealized Kob-Andesen binay Lennad- Jones (KABLJ) mixues, he Wahnsom binay Lennad-Jones (WABLJ) mixue, he single-componen Buckingham (SCB) liquid, he Lewis-Wahnsom OTP model, he asymmeic dumbbell model. These sysems ae all songly coelaing [13 17]. Thus fo songly coelaing liquids i is enough o know he inemolecula ineacions wihin he FCS in ode o accuaely simulae sucue and dynamics. The simulaions fuhe showed ha fo all of he following aomic and molecula sysems The Dzuguov (DZ) liquid, he Lennad-Jones Gaussian (LJG) liquid, he Gaussian coe model (GCM), he Hansen-McDonald molen sal model,

13 13 he igid SPC/E wae model, sucue and/o dynamics ae no popely epoduced in simulaions when an FCS cuoff is used fo all hese sysems lage cuoffs ae needed fo faihful epoducion of he sysem s physics. None of he lae liquids ae songly coelaing. An ineesing limiing case is he SCLJ liquid, which a low densiies ( a ypical empeaues) sops being songly coelaing concuenly, simulaions using an FCS cuoff become poo (Figs. 5 and (c)). In conclusion, a shifed-foces FCS cuoff leads o accuae esuls if and only if he liquid in quesion is songly coelaing a he simulaed sae poin. We know of no excepions o his ule. By inducion his suggess ha songly coelaing liquids ae chaaceized by he popey ha inemolecula ineacions beyond he FCS can be safely ignoed in simulaions. IV. DISCUSSION: THE REAL ESSENCE OF LIQUID SIMPLICITY As menioned in he inoducion a definiion of simple liquids is mos useful if i idenifies hei eal essence in Locke s undesanding of he em [3], he undelying fundamenal chaaceisic fom which hese liquids simple feaues, hei nominal essences, follow. We sugges below ha he class of simple liquids is o be idenified wih he class of songly coelaing liquids (Sec. IV C). This is moivaed by fis summaizing biefly he many simple popeies of songly coelaing liquids, hen showing ha his class of liquids can be chaaceized fom hee diffeen pespecives: mahemaically, physically, and chemically (Sec. IV B). This gives ise o hee vey diffeen chaaceizaions, indicaing ha he class of songly coelaing liquids is fundamenal and moivaing he suggesion ha he eal essence of liquid simpliciy is he exisence of song coelaions of viial/poenial-enegy equilibium N V T flucuaions. Finally, by connecing o he chemiss concep of nonassociaed liquids we discuss biefly which eal-wold liquids ae simple (Sec. IV D) and give examples of complex liquid popeies (Sec. IV E). A. Songly coelaing liquids An example of a songly coelaing liquid is he pooype simple liquid Agon a supeciical sae poins, whee R = 0.96 fo empeaues beween 200 and 660K and densiies beween 20 and 30 mol/l [9, 14]. All popeies of songly coelaing liquids follow fom he exisence of isomophs in hei phase diagam (see below). Some of hei simple popeies, howeve, wee idenified befoe isomophs wee defined in 2009 [16], fo insance ha all eigh fundamenal hemoviscoelasic esponse funcions ae given in ems of jus one, i.e., he dynamic Pigogine-Defay aio is close o uniy [12], aging may be descibed by adding meely one exa paamee [14, 15], powe-law densiy scaling [51] is obeyed o a good appoximaion, i.e., fo vaying densiy and empeaue he elaxaion ime is a funcion of ρ γ /T [19]. An isomoph is an equivalence class of sae poins in a sysem s phase diagam. Two sae poins (ρ 1, T 1 ) and (ρ 2, T 2 ) ae defined o be isomophic [16] if he following holds: Wheneve one micoconfiguaion of sae poin (1) and one of sae poin (2) have he same educed coodinaes (i.e., ρ 1/3 1 (1) i = ρ 1/3 2 (2) i fo all paicles i), hese wo micoconfiguaions have popoional configuaional Bolzmann facos, e U((1) 1,...,(1) N )/k BT 1 = C 12 e U((2) 1,...,(2) N )/k BT 2. (5) Fo mos sysems he isomoph concep is appoximae jus as he W U coelaions ae no pefec hus we do no equie Eq. (5) o be igoously obeyed fo all micoconfiguaions, bu only o a good appoximaion fo all physically elevan micoconfiguaions. By his is mean micoconfiguaions ha ae no a pioi unimpoan fo he physics. Being an equivalence class, an isomoph defines a coninuous cuve of sae poins in he liquid s phase diagam. Only liquids fo which he poenial enegy is an Eule homogeneous funcion, e.g., sysems wih IPL pai poenials, have exac isomophs. An IPL fluid wih exponen n has exac isomophs chaaceized by ρ γ /T = Cons. whee γ = n/3. Appendix A of Ref. 16 showed ha a liquid is songly coelaing if and only if i has isomophs o a good appoximaion. This was confimed by Ref. 16 and 17, which showed ha Lennad-Jones ype aomic liquids have good isomophs, as well as Ref. 52, which showed ha he songly coelaing Lewis-Wahnsom OTP and asymmeic dumbbell models have good isomophs. Equaion (5) has many consequences. These wee deived and discussed in deail in he oiginal isomoph pape fom 2009 (Ref. 16), o which he eade is efeed. Basically, sucue and dynamics a wo isomophic sae poins ae idenical in educed unis. In moe deail, quaniies ha ae invaian along an isomoph include (bu ae no limied o): 1. The excess enopy, i.e., he enopy in excess of he ideal gas enopy a he same densiy and empeaue his is he configuaional conibuion o he enopy (a em ha is negaive because a liquid is always moe odeed han an ideal gas a he same densiy and empeaue).

14 14 2. All N-body enopy ems, whee we emind he eade ha he excess enopy can be expanded in a seies of wo-body, hee-body, ec, ems; each em is invaian along an isomoph [16]. 3. The isochoic hea capaciy. 4. The sucue in educed unis (defined by i ρ 1/3 i fo all paicles i). No only he adial disibuion funcion in educed unis, bu all higheode disibuion funcions, ae isomoph invaian. 5. The Newonian NVE and Nosé-Hoove NVT equaions of moion in educed unis; likewise Bownian dynamics. 6. All auocoelaion funcions in educed unis. 7. All aveage elaxaion imes in educed unis. 8. Reduced anspo coefficiens like he diffusion coefficien, viscosiy, ec. Isomophs have he fuhe ineesing popey ha hee is no elaxaion fo an insananeous change of empeaue and densiy jumping fom an equilibaed sae poin o a diffeen sae poin ha is isomophic wih he iniial sae poin. The absence of elaxaion deives fom he fac ha he Bolzmann pobabiliies of scaled micoconfiguaions ae idenical. Simulaions have shown ha such isomoph jumps wok quie well fo he KABLJ liquid [16], fo he asymmeic dumbbell, and fo he Lewis-Wahnsom OTP molecula models [52]. Moeove, i has been shown ha he effecive empeaue of a glass pepaed by a empeaue-densiy jump fom an equilibium sae of a songly coelaing liquid depends only on he final densiy [11]; his povides ye anohe example of a simple feaue of hese liquids ha does no apply geneally. Some fuhe pedicions fo he class of songly coelaing liquids deiving fom he exisence of isomophs ae: The solid-liquid coexisence cuve is an isomoph [16, 17]. This implies invaiance along he coexisence cuve of he educed sucue faco, he educed viscosiy, he educed diffusion consan, ec, as well as pessue invaiance of he meling enopy, and a educed-uni Lindemann meling cieion [16]. Collapse of he wo-ode-paamee maps of Debenedei e al. [53 57] o one-dimensional cuves [16]. Isochonal supeposiion [58], i.e., he fac ha when pessue and empeaue ae vaied, he aveage elaxaion ime deemines he enie elaxaion specum [16]. The above lised popeies of songly coelaing liquids all eflec one o he ohe kind of simpliciy of songly coelaing liquids. A final, ecenly esablished simple popey is a hemodynamic sepaaion ideniy: fo all songly coelaing liquids, if s is he excess enopy pe paicle, he empeaue as a funcion of s and densiy ρ facoizes as follows [59] T = f(s)g(ρ). (6) This is a saighfowad consequence of he fac ha he cuves in he phase diagam of consan excess enopy coincide wih he cuves of consan isochoic hea capaciy, a popey of all songly coelaing liquids because boh quaniies ae isomoph invaians (poins 1 and 3 above). Equaion (6) implies 1) ha he densiyscaling exponen is a funcion of densiy only, 2) a configuaion-space vesion of he Güneisen equaion of sae, 3) ha fo sysems wih pai ineacions he nonzeo ems in he expansion g(ρ) = n C nρ n/3 ae pecisely hose appeaing in he pai poenial expanded as v() = n v n n, and 4) ha he isomophs of LJ liquids in paicula, he solid-liquid coexisence cuve ae given by (Aρ 4 Bρ 2 )/T = Cons. [60]. B. Mahemaical, physical, and chemical chaaceizaions of songly coelaing liquids A a given sae poin, if he aveage poenial enegy is denoed by U, he consan-poenial-enegy hypesuface is defined by Ω = {( 1,..., N ) R 3N U( 1,..., N ) = U }. This is a compac, Riemannian (3N 1)- dimensional diffeeniable manifold. Each sae poin has is own such hypesuface. In his way, a family of manifolds is defined houghou he phase diagam. In Appendix A of Ref. 16 i was shown ha he educed-uni consan-poenial-enegy manifold is invaian along a songly coelaing liquid s isomophs, and invaiance cuves in he phase diagam fo hese manifolds exis only fo songly coelaing liquids. Thus, fo such liquids hese manifolds consiue a one-paamee family of manifolds, no wo-paamee families as expeced fom he fac ha he phase diagam is wo-dimensional. This povides a mahemaical chaaceizaion of he class of songly coelaing liquids. The physical chaaceizaion of his class was discussed aleady: he exisence of isomophs in he phase diagam. A liquid is songly coelaing if and only if i has isomophs o a good appoximaion; his is shown in Appendix A of Ref. 16. The poof uilizes ha a liquid is songly coelaing if and only if is consan-viial hypesufaces in he 3N dimensional configuaion space ae almos idenical o is consan-poenial-enegy hypesufaces. The chemical chaaceizaion of songly coelaing liquids is he novel popey documened in he pesen pape: A liquid is songly coelaing a a given sae

15 15 poin if and only if he liquid s sucue and dynamics ae accuaely calculaed by simulaions ha ignoe ineacions beyond he fis coodinaion shell. This is an empiical finding fo which we, in conas o he mahemaical and physical chaaceizaions, have a pesen no compelling agumen. How can one undesand he FCS chaaceizaion of songly coelaing liquids? Noe fis ha he popey of insignificance of ineacions beyond he FCS is an isomoph invaian: If a liquid has good isomophs and if an FCS cuoff woks well a one sae poin, FCS cuoffs mus wok well fo all is isomophic sae poins. Thus he chemical chaaceizaion of songly coelaing liquids is nicely consisen wih he fac ha hese liquids have isomophs. Noe fuhe ha i has been shown fo he Lennad- Jones liquid ha almos all of he flucuaions in viial and poenial enegy come fom inepaicle sepaaions wihin he FCS [14]. The new FCS chaaceizaion of songly coelaing liquids shows ha hese liquids ae chaaceized by having a well-defined FCS; mos likely i is he exisence of a well-defined FCS which implies he almos cancellaion of he linea em of he shifedfoce poenial. The fac ha ineacions beyond he FCS may be ignoed shows ha ineacions ae effecively sho anged, which means ha he sucue is dominaed by wha may be emed packing effecs. C. Defining he class of simple liquids Secion IV A lised songly coelaing liquids many simple popeies. Secion IV B showed ha his liquid class may be chaaceized fom hee quie diffeen poins of view. Taken ogehe i appeas ha he class of songly coelaing liquids is fundamenal. Since he popeies of songly coelaing liquids ae geneally simple han hose of liquids in geneal, we now popose he definiion Simple liquids = Songly coelaing liquids This is he basic message of he pesen pape, which implies a quanificaion of he degee of simpliciy via he numbe R of Eq. (2), he NV T ensemble equilibium viial/poenial-enegy coelaion coefficien. Befoe poceeding o sugges ha he above be supplemened by an analyiciy equiemen of he poenial, we noe ha compaed o he commonly used definiion of simple liquids as hose wih adially symmeic pai ineacions, hee ae some noable diffeences: 1. Simpliciy is no an on/off popey, bu is quanified by a coninuously vaiable. 2. The degee of simpliciy is no a unique popey of a liquid s inemolecula ineacions, bu geneally vaies houghou he phase diagam. Indicaions ae ha mos songly coelaing liquids loose his popey a low densiies. 3. No all aomic liquids (i.e., wih adially symmeic pai ineacions) have simple egions in hei phase diagam (compae he above examples of he Dzuguov, Lennad-Jones Gaussian, Gaussian coe, and molen sal models); 4. No all simple liquids ae aomic (compae he Wahnsom OTP and he asymmeic dumbbell models). Accoding o he above poposed definiion of liquid simpliciy liquids, he case whee he poenial enegy is an Eule homogeneous funcion of he paicle posiions (R = 1) se he gold sandad fo simpliciy. This is consisen wih he many simple popeies of hese liquids. Usually one has in mind IPL sysems wih pai ineacions scaling wih disance as n, bu R = 1 sysems include also some wih angula dependences in he poenial enegy, as long as hese scale wih disance he same way as all ohe ineacions. Due o he absence of aacions, IPL fluids have no liquid-gas phase ansiion. In his sense i may seem sange o claim ha IPL fluids ae he simples liquids. Howeve, we have found ha moe ealisic songly coelaing liquids cease o be so when he liquid-vapo coexisence line is appoached, showing ha his phase ansiion canno be undesood in he famewok of simple liquids. This conass wih he liquid-solid phase ansiion, whee fo insance he fac ha he coexisence line fo simple liquids is an isomoph confimed by simulaion fo he LJ liquid [17] explains seveal peviously noed egulaiies [16]. Is he had-sphee fluid simple? One may define a configuaional viial funcion fo his sysem, bu i is no obvious how o define a poenial enegy funcion ha is diffeen fom zeo. Thus hee is no meaningful coelaion coefficien R fo had-sphee fluids. On he ohe hand, he had-sphee liquid may be egaded as he n limi of an IPL liquid, and i is well known ha fo insance he had-sphee adial disibuion funcion is close o ha of, e.g., an 20 IPL liquid if he effecive had-sphee adius of he IPL sysem is chosen popely. This would indicae ha had-sphee liquids ae simple, consisen wih he pevailing poin of view. Anohe ineesing case is ha of he Weeks-Chandle-Andesen (WCA) vesion of he LJ liquid, which cus off all aacions by puing he foce equal o zeo beyond he poenial enegy minimum. This liquid is songly coelaing [61] indeed moe so han he SCLJ liquid bu despie his we found in simulaions ha he WCALJ liquid has somewha pooe isomophs han he LJ liquid. We d like o sugges ha he had-sphee liquid and he WCALJ liquid should be boh excluded fom he class of simple liquids on he gounds ha hei poenials ae no analyic. Fo sysems ineacing wih pai poenials, fo insance, i makes good sense o add he exa equiemen ha he pai poenial is an analyical funcion of he invese pai disance, i.e., ha an expansion exiss of he fom