arxiv: v2 [cond-mat.soft] 27 Jan 2012

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: May 23, 2018) axiv: v2 [cond-ma.sof] 27 Jan 2012 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 coelaion 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. Fo insance, in he new definiion simpliciy is no a popey of he inemolecula poenial because a liquid is usually only songly coelaing in pa of is phase diagam. Moeove, accoding o he new definiion no all simple liquids ae aomic (i.e., wih adially symmeic pai poenials) and no all aomic liquids ae simple. The main pa of he pape moivaes he new definiion of liquid simpliciy by pesening evidence 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 aomic and molecula model liquids wih a shifed-foces cuoff placed a he fis minimum of he adial disibuion funcion. The liquids sudied ae: invese powe-law sysems ( n pai poenials wih n = 18,6,4), Lennad- Jones (LJ) models (he sandad LJ model, wo genealized Kob-Andesen binay LJ mixues, he Wahnsöm binay LJ mixue), he Buckingham model, he Dzuguov model, he LJ Gaussian model, he Gaussian coe model, he Hansen-McDonald molen sal model, he Lewis-Wahnsöm OTP model, he asymmeic dumbbell model, and he igid SPC/E wae model. The final pa of he pape summaizes mos known popeies of songly coelaing liquids, showing ha hese ae simple han liquids in geneal. Simple liquids as defined hee may be chaaceized (1) chemically by he fac ha he liquid s popeies ae fully deemined by ineacions fom he molecules wihin he FCS, (2) physically by he fac ha hee ae isomophs in he phase diagam, i.e., cuves along which seveal popeies like excess enopy, sucue, and dynamics, ae invaian in educed unis, and (3) mahemaically by he fac ha he educed-coodinae consan-poenial enegy hypesufaces define a one-paamee family of compac Riemannian manifolds. No poof is given ha he chemical chaaceizaion follows fom he song coelaion popey, bu i is shown o be consisen wih he exisence of isomophs in songly coelaing liquids phase diagam. Finally, we noe ha he FCS chaaceizaion of simple liquids calls ino quesion he basis fo sandad peubaion heoy, accoding o which he epulsive and aacive foces play fundamenally diffeen oles fo he physics of 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 defined caegoies as 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, 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 eal essence, we ague below fo a new definiion of he class of simple liquids. Physiciss love simple sysems. This eflecs he fundamenal paadigm 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. The concep of a simple liquid is old, bu i emains cenal as evidenced by he 2003 book ile Basic conceps fo simple and complex liquids [5] o he eview of non-simple liquids eniled Theoy of complicaed liquids fom 2007 [6]. Geneaions of liquid-sae heoiss wee inoduced o he opic by sudying Hansen and Mc- Donald s exbook Theoy of simple liquids [7]. This book fis appeaed in 1976 following a peiod of specacula pogess in he heoy of liquids, caalyzed by some of he vey fis scienific compue simulaions. In Ref. 7 a simple liquid is defined as a classicalsysem of appoximaely spheical, nonpola molecules ineacing via pai poenials. This and closely elaed definiions of liquid simpliciy have been sandad fo many yeas [8 12]. In his definiion simple liquids have much in common wih he chemiss nonassociaed liquids [13], bu hee ae some diffeences. Chemiss geneally

2 2 egad a liquid as simple even if i consiss of elongaed molecules, as long as hese ae wihou inenal degees of feedom and ineac pimaily via van de Waals foces. Many physiciss would pobably disagee. Thus i is fa fom ivial o ask: Wha chaaceizes a simple liquid? Moe accuaely: Given a classical sysem of igid bodies wih poenial enegy as a funcion of he bodies cenes-of-masses and hei spaial oienaions, is i possible o give a quaniaive cieion fo how simple he sysem is? If yes, is simpliciy encoded uniquely in he poenial enegy funcion o may he degee of simpliciy vay houghou he phase diagam? Recen woks idenified and descibed he popeies of songly coelaing liquids [14 25]. In hese liquids, by definiion, he viial W and he poenial enegy U coelae songly in hei consan-volume hemalequilibium flucuaions. Recall ha he aveage viial W gives he conibuion o pessue fom inemoleculaineacions, addedohe ideal-gasemnk B T deiving fom momenum anspo via paicle moion(below p is he pessue, V he volume, N he numbe of paicles, k B Bolzmann s consan, and T he empeaue): pv = Nk B T + W. (1) The em songly coelaing liquid efes o he case when he WU coelaion coefficien in he NVT ensemble is lage han 0.9 [18]. If angula backes denoe an NVT 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, i.e., no fa fom he solid-liquid coexisence line. Many ohe sysems, including some 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 [14 25], eviewed biefly in Ref. 25. 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 [26, 27]. These auhos showed ha fo he viscous Kob-Andesen binay Lennad-Jones mixue [28, 29] he dynamics is no epoduced popely by cuing he poenials a hei minima accoding o he well-known Weeks-Chandle- Andesen (WCA) ecipe [30]. The ole of he cuoff was subsequenly sudied in wo papes, showing ha placingashifed-focescuoffa he fisminimum ofhe pai coelaion funcion gives good esuls fo Lennad-Jones ype sysems [31, 32]. This applies no only a modeae densiies, bu also a vey high densiies. Applying he same cuoff o wae does no wok vey well [33]. Wae is an example of a non-songly coelaing liquid wih R 0 a ambien condiions, a consequence of wae s densiy maximum [18]. These findings led us o speculae whehe i is a geneal popey of songly coelaing liquids ha he inemolecula ineacions may be ignoed beyond he FCS wihou compomising accuacy o any significan exen. The main pa of he pesen pape shows ha, indeed, using such an FCS cuoff gives accuae simulaion esuls if and only if he liquid is songly coelaing. The pape pesens esuls obained fom compue simulaions of 15 diffeen sysems, some of which ae songly coelaing. We invesigae he ole of he FCS in deemining liquid sucue and dynamics. Sucue is pobed by he adial disibuion funcion (RDF), dynamics by he incoheen o coheen inemediae scaeing funcion (ISF) a he waveveco of he saic sucue faco maximum. The numeical evidence 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. A pesen hee is no igoous agumen fo his empiical FCS popey, bu we show ha i is consisen wih known popeies of songly coelaing liquids. The FCS popey of songly coelaing liquids documened below emphasizes fuhe ha hese ae simple han liquids in geneal. A numbe of ohe simple popeies of songly coelaing liquids wee epoed peviously [14 22, 24, 25]. Alogehe, hese facs moivae ou new definiion of liquid simpliciy. Secion II pesens he esuls fom molecula dynamics simulaions and Sec. III summaizes he esuls. Secion IV gives an oveview of he many simple popeies of songly coelaing liquids, moivaing ou 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 inemolecula ineacions, which usually exend in pinciple o infiniy, ae uncaed a some cuoff disance c beyond which hey ae ignoed [31]. To avoid a disconinuiy in he foce, which can seveely affec he simulaion esuls [31, 34], he simulaions epoed below 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 [34 36] whee, if he pai poenial is v() and he pai foce is f() = v (), he shifed foce is given by

3 3 f SF () = { f() f( c ) if < c, 0 if > c. (3) he case fo sysems wih invese powe-law (IPL) pai poenials (v() n ), bu noe ha poenials wih non-ivial angula dependence may also be Eule homogeneous. This coesponds o using he following pai poenial below c : v SF () = v() v ( c )( c ) v( c ). Using an SF cuoff gives moe accuae esuls and bee numeical sabiliy han using he sandad shifed-poenial (SP) cuoff [31]. This is so despie he fac ha an SF cuoff does no have he coec pai foce fo any, wheeas he pai foce is coec below c fo an SP cuoff. Appaenly, avoiding disconinuiy of he foce a c is moe impoan han mainaining he coec foce. I was ecenly discussed why adding a linea em o he pai poenial affecs neihe sucue no dynamics o any significan exen [19]. 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 [19] (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-angedand condiionally convegen, when c is sufficienly lage an SF cuoff gives esuls close o hose of he sandad, much moe involved Ewald summaion mehod [37, 38]. All simulaions wee pefomed in he N V T ensemble wih peiodic bounday condiions and he Nose-Hoove algoihm [39 41]. We used he RUMD molecula dynamics package developed in-house, which is opimized fo simulaions of small sysems on sae-of-he-a GPU hadwae (a few housand eniies) [42]. Fo he molecula models bond lenghs wee held fixed using he imesymmeical cenal diffeence algoihm [43 45]. 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 songly coelaing aomic liquids. Then daa ae pesened fo aomic liquids ha ae no songly coelaing. Finally daa ae given fo wo songly coelaing molecula liquids and a wae model. Deails of he models sudied, numbe of paicles, ec, ae given in Appendix A. We simulaed single-componen IPL pai-poenial sysems wih exponens n = 18,6,4 a densiy ρ = 5. Each sysem was simulaed a wo empeaues. The simulaed sysems ange fom vey hash epulsive (n = 18) o quie sof and long anged (n = 4). Each sysem was simulaed wih fou SF cuoffs. The ole of he cuoff is invesigaed by choosing hee diffeen cuoffs: one placed a he fis minimum of he RDF, one coesponding o he half heigh of he RDF fom is fis maximum o he fis minimum, and one placed o he igh of he RDF fis minimum displaced he same amoun as he diffeence beween he fis and he second cuoff. The RDFs g() ae shown fo n = 18,6,4 in Fig. 1; n = 12 gives simila esuls (no shown). The simulaions wih an SF cuoff a he fis minimum of he RDF efeed o as FCS-cuoff simulaions give a faihful epesenaion of he sucue. The inses show he deviaions in RDF as funcions of he cuoff beween esuls fo an SF cuoff and ue esuls, quanified by inegaing he numeical diffeence in he pai coelaion funcion. Clealy, deviaions incease shaply when he cuoff enes he FCS (blue cosses). We simulaed also he n = 3 and n = 1 IPL fluids (esuls no shown), he lae is usually emed he onecomponen plasma. Fo boh sysems an FCS cuoff does no lead o he coec physics. This is mos likely because he poenials ae oo long anged fo such a sho cuoff o make good sense. In fac, boh ae so long anged ha hey do no have a pope hemodynamic limi fo which he exponen mus be lage han he dimension [46]. An indicaion ha an FCS cuoff woks pooly when he IPL exponen appoaches he dimension is seen fo he n = 4 simulaion, fo which he WU coelaion coefficien fo he FCS sysem sas o deviae significanly fom uniy. Moeove, almos invisible in he figue is he fac ha he n = 4 pai coelaion funcion s fis maximum deviaes slighly beween FCS and ue simulaions. A. Thee invese-powe-law (IPL) fluids We conside fis sysems wih 100% coelaion beween viial and poenial enegy equilibium flucuaions in he NVT ensemble. I follows fom he definiion of he viial W = 1/3 i i U [34] ha a 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 g() T = 0.7: c = T = 0.7: c = T = 0.3: c = T = 0.3: c = g () - g() d T = 0.7 T = c IPL 18 (ρ = 5) c = 5.300) = 0 = 1.559) =

4 4 g() T = 0.70: c = T = 0.70: c = T = 0.15: c = T = 0.15: c = g () - g() d T = 0.70 T = c IPL 6 (ρ = 5) c 0 R( c = 5.300) = 0 = 1.570) = ( q = 6.64, ) IPL 6 (ρ = 5, T = 0.15) c = c = c = (min) c = = 5.300) = 0 = 1.570) = g() T = 0.7: c = T = 0.7: c = T = 0.1: c = T = 0.1: c = g () - g() d T = 0.7 T = c c IPL 4 (ρ = 5) = 5.300) = 0 (c) = 1.580) = FIG. 1. RDFs fo single-componen IPL fluids wih exponens n = 18,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; 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.61, ) IPL 4 (ρ = 5, T = 0.10) (c) c = c = c = (min) c = = 5.000) = 0 = 1.580) = FIG. 2. Incoheen inemediae scaeing funcions (ISFs) fo he IPL fluids a he lowes-empeaue sae poins of Fig. 1. The black cuves give esuls fo a lage cuoff, he ed cosses fo an FCS cuoff (maked by he veical ed dashed lines in Fig. 1). n = 18, T = 0.30; n = 6, T = 0.15; (c) n = 4, T = B. Lennad-Jones ype liquids Nex, we conside he mos sudied poenial in he hisoy of compue simulaions, he Lennad-Jones (LJ) pai poenial, Figue 2 shows he incoheen inemediae scaeing funcions (ISFs) fo he low-empeaue sae poins of each of he hee IPL sysems evaluaed a he waveveco coesponding o he fis maximum of he saic sucue faco. A good epesenaion of he dynamics is obained fo all sysems when he FCS cuoff is used. ( q = 6.72, ) IPL 18 (ρ = 5, T = 0.30) c = c = c = (min) c = = 5.300) = 0 = 1.559) = [( σ ) 12 ( σ 6 ] v LJ () = 4ǫ. (4) ) Hee σ and ǫ, define, especively, he lengh and enegy scale of he ineacion (dimensionless unis defined by σ = ǫ = 1 ae used below). This poenial does no have 100% viial/poenial-enegy coelaion, bu has sill quie song coelaions wih coelaion coefficiens R > 0.9 in he condensed-fluid pa of he phase diagam (also in he cysalline phase [19]). We sudied he single-componen LJ (SCLJ) liquid, wo genealized 80/20 Kob-Andesen binay LJ (KABLJ) mixues wih epulsive exponen 12 and aacive exponens n = 4, 10, and he Wahnsom 50/50 binay LJ (WABLJ) mixue (Fig. 3 and Appendix A gives model deails). The influence of an SF cuoff on simulaion accuacy was invesigaed ecenly fo he SCLJ liquid and he sandad KABLJ mixue (n = 6) [31, 32], bu fo compleeness we include esuls fo he SCLJ sysem hee as well.

5 5 v() AA Genealized (12, 04)-KABLJ AA Genealized (12, 06)-KABLJ AA Genealized (12, 10)-KABLJ g() 4.0 c = c = g () - g() d c SCLJ (ρ = 5, T = 0) - Genealized Lennad-Jones: v() = ε/(m-n) [ n (σ/) m - m (σ/) n ] 4.0 FIG. 3. 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 = 1.50, and ǫ BB = 0, σ AA = 1.12, σ AB = 0.90, σ BB = 0.99, m A = m B = 1. The ole of he cuoff is again invesigaed by choosing hee diffeen cuoffs: one placed a he fis minimum of he RDF (ed colo in Figs. 4-8), one coesponding o he half heigh of he RDF fom maximum o he fis minumum (blue colo in Figs. 4-8), and one displaced o he igh of he minimum he same amoun as he diffeence beween he fis and he second cuoff (geen colo in Figs. 4-8). In Fig. 4 RDFs ae shown fo he SCLJ liquid a hee diffeen sae poins. 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 Figs. 4 and is clealy epesened well using an FCS cuoff, while choosing he cuoff inside he FCS esuls in significan deeioaion. A low densiies [Figs. 4(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 sae poin (c) whee R = 0. These simulaions sugges ha only when a liquid is songly coelaing is i possible o ignoe ineacions beyond he FCS. g() = 5.000) = 0.97 = 1.544) = c = c = g () - g() d c SCLJ (ρ = 5, T = 1.13) = 6.000) = 0 = 1.645) = 9 (c) FIG. 4. RDFs fo he single-componen Lennad-Jones (SCLJ) liquid a hee diffeen sae poins: ρ = 5, T = 0.70 (R = 0.96); ρ = 5, T = 0 (R = 0.97); (c) ρ = 5, T = 1.13 (R = 0). The black cuves show efeence simulaions wih lage cuoffs, he ed dos/cuve 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). Fo sae poins and, which ae songly coelaing (R > 0.9), an FCS cuoff leads o accuae esuls. This is no he case fo sae poin (c) ha is no songly coelaing. Nex, we invesigaed he SCLJ dynamics a he hee sae poins of Fig. 4. The dynamics is sudied via he ISF. The ISFs ae shown in Fig. 5; a all sae poins he dynamics is epesened well using an FCS cuoff. g() 4.0 c = c = g () - g() d c SCLJ (ρ = 5, T = 0.70) = 5.000) = 0.96 = 1.538) = ( q = 6.788, ) SCLJ (ρ = 5, T = 0.70) c = c = c = (min) c = = 5.000) = 0.96 = 1.538) =

6 6 Sysem ρ T R RDF max ISF max SCLJ SCLJ SCLJ SCLJ SCLJ TABLE I. Simulaion esuls fo five sae poins of he SCLJ liquid. Fo each sae poin is given densiy, empeaue, coelaion coefficien, maximum deviaion fom he ue RDF using an FCS cuoff, and maximum deviaion fom he ue incoheen sucue faco using an FCS cuoff. ( q = 6.788, ) ( q = 6.42, ) SCLJ (ρ = 5, T = 0) c = c = c = (min) c = = 5.000) = 0.97 = 1.544) = SCLJ (ρ = 5, T = 1.13) (c) c = c = 40 c = (min) c = = 6.000) = 0 = 1.645) = FIG. 5. ISFs fo he SCLJ liquid a he sae poins of Fig. 4. 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. ρ = 5, T = 0.70 (R = 0.96); ρ = 5, T = 0 (R = 0.97); (c) ρ = 5, T = 1.13 (R = 0). The n = 4 case is below ou definiion of a songly coelaing liquid R > 0.9 (a he simulaed sae poin his sysem has negaive viial), bu is sill pey well coelaing. Fo all invesigaed sae poins an FCS cuoff gives accuae esuls. We found he same using he sandad epulsive exponen n = 6 [32] (esuls no shown). g AA () g AA () c = c = g () - g() d c (12, 4)-KABLJ (ρ = 1.20, T = 0) = 4.700) = 4 = 1.433) = c = c = g () - g() d c (12, 10)-KABLJ (ρ = 1.20, T = 0) = 4.700) = 0.97 = 1.426) = 0.96 FIG. 6. RDFs fo genealized KABLJ mixues wih epulsive exponen 12 and aacive exponens n = 4,10. 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 = 10, ρ = 1.20, T = 0 (R = 0.97). Table I summaizes esuls fo five sae poins of he SCLJ liquid. The hid sae poin is slighly below he iple poin, i.e., in a slighly measable sae. The deviaions clealy incease as he W U coelaions decease. We poceed o invesigae mixues of wo diffeen paicles (A and B) ineacing wih LJ ype poenials. The cuoff used fo all hee ineacions (AA, AB, BB) ae placed a he same disance, efeing o σ AA. In Fig. 6 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, 10. The A-paicle ISFs fo he sae poins of Fig. 6 ae shown in Fig. 7. Fo he KABLJ mixue, placing he cuoff inside he FCS (blue cuves) fails o epoduce he dynamics popely, wheeas i 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 ae found fo he B paicles (esuls no shown).

7 7 A ( q = 7.23, ) (12, 4)-KABLJ (ρ = 1.20, T = 0) c = c = c = (min) c = = 4.700) = 4 = 1.433) = 0.95 In summay, fo all LJ ype sysems wheneve 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 WU coelaions ae weake, an FCS cuoff gives poo esuls. A ( q = 7.21, ) (12, 10)-KABLJ (ρ = 1.20, T = 0) c = c = c = (min) c = = 4.700) = 0.97 = 1.426) = FIG. 7. ISFs fo genealized KABLJ mixues wih epulsive exponen 12 and aacive exponens n = 4,10. The ed and black cuves give, especively, esuls fo FCS cuoffs and lage efeence cuoffs. n = 4, ρ = 1.20, T = 0 (R = 4); n = 10, ρ = 1.20, T = 0 (R = 0.97). We also simulaed he Wahnsöm binay LJ mixue [47], finding again ha wheneve R > 0.9 sucue and dynamics ae well epoduced using an FCS cuoff. We do no show hese esuls, bu show insead esuls fo he coheen scaeing funcion a one sae poin (Fig. 8). Again, he FCS cuoff (ed cosses) gives he coec dynamics wheeas educing he cuoff fuhe does no give pope esuls (blue cosses). F AA ( q = 7.48, ) Wahnsöm-BLJ (ρ = 1.30, T = 0) c = c = c = (min) c = = 4.500) = 0.99 = 1.427) = FIG. 8. Inemediae coheen scaeing funcion fo he Wahnsöm 50/50 binay LJ liquid a he waveveco coesponding o he saic sucue faco maximum. The ed and black cuves give, especively, esuls fo a FCS cuoff and lage efeence cuoff. C. Buckingham liquid Nex, we conside he single-componen Buckingham liquid (SCB). The Buckingham poenial [48, 49] 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 (Fig. 9). The paamees of he Buckingham poenial (Appendix A) wee chosen such ha he LJ poenial is well appoximaed in he epulsive egion (Fig. 9). v() - v() Buckingham: (12, 6)-Lennad Jones Buckingham v() = ε( 6/(α-6) exp[α(1 - / m )] - α/(α - 6) ( m /) 6 ) FIG. 9. 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 10 and show, especively, he RDF and ISF fo he SCB liquid. The FCS cuoff woks well. g() 4.0 c = c = g () - g() d c Single-componen Buckingham (ρ = 0, T = 0) = 4.900) = 0.99 = 1.475) =

8 8 ( q = 7.10, ) Single-componen Buckingham (ρ = 0, T = 0) c = c = c = (min) c = = 4.900) = 0.99 = 1.475) = FIG. 10. 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 Figue 11 shows he Dzuguov (DZ) pai poenial [50], which was oiginally suggesed as a model poenial impeding cysallizaion by enegeically punishing paicle sepaaions coesponding o he nex-neaes neighbo disance of a face-ceneed cubic laice. A sho disances he DZ pai poenial appoximaes he LJ poenial. (12, 6)-Lennad-Jones (shifed) Dzuguov g() F( q = 6.94, ) g () - g() d c = 1.940) = 0.71 = 1.602) = c = c = (min) c = c = c = Dzuguov (ρ = 0, T = 0.75) Dzuguov (ρ = 0, T = 0.75) = 1.940) = 0.71 = 1.602) = FIG. 12. The effec on sucue and dynamics of vaying he cuoff fo he Dzuguov (DZ) liquid a ρ = 0 and T = 0.75 (R = 0.71). The ed and black cuves give, especively, esuls fo an FCS cuoff and a lage efeence cuoff. RDF whee he inse quanifies he deviaion fom he efeence RDF as a funcion of he cuoff. The inemediae coheen scaeing funcion a he same sae poin, including hee esuls fo a cuoff wihin he FCS (blue cosses). Dzuguov: u() = A( -n - B) exp[c / ( - a )] + B exp[d / ( - b)] u() FIG. 11. The Dzuguov (DZ) pai poenial [50] (blue). Also shown is he LJ pai poenial(black cuve). The DZ poenial appoximaes he LJ poenial aound he fis minimum, bu has a local maximum a lage disances. E. Lennad-Jones Gaussian liquid Figues 12 and show, especively, he RDF and inemediae coheen scaeing funcion 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. Wha is impoan hee, howeve, is ha he poo FCS cuoff esuls coelae wih he faily weak viial/poenial-enegy coelaions (R = 0.71). The Lennad-Jones Gaussian liquid [51] (LJG) is a non-songly coelaing liquid wih he wo-minimum pai poenial shown in Fig. 13. The paamees of LJG model (Appendix A) ae such ha he second LJG poenial minimum does no coincide wih ha of he SCLJ sysem [52].

9 9 v() - - Lennad-Jones Gaussian Lennad-Jones Gaussian: v() = v LJ () + εε 0 exp[-(- 0 ) 2 /2σ 0 2 ] 4.0 FIG. 13. The Lennad-Jones Gaussian (LJG) pai poenial [51] consuced by adding a Gaussian o he pai LJ poenial. Two disinc minima ae pesen. Resuls fom simulaing sucue and dynamics of he LJG liquid ae shown in Figs. 14 and. The FCS cuoff does no give he coec RDF, wheeas deviaions inhedynamicsaefailysmall. NoehaheFCScuoff emoves he second minimum. g() ( q = 8.7, ) c = c = g () - g() d c Lennad-Jones Gaussian (ρ = 5, T = 0) = 5.000) = 4 = 1.242) = = 5.000) = 4 = 1.242) = 0.97 c = c = c = (min) c = 38 Lennad-Jones Gaussian (ρ = 5, T = 0) F. Gaussian coe model The Gaussian coe model (GCM) [53, 54], which is no songly coelaing, is defined by a Gaussian pai poenial and hus has a finie poenial enegy a zeo sepaaion. The high-densiy egime of he GCM model (ρ > 1.5) has ecenly eceived aenion as a singlecomponen model glass fome[55] because i is no pone o cysallizaion and shows he chaaceisic feaues of glass-foming liquids(lage viscosiy, wo-sep elaxaion, ec). Figue 15 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. g() ( q = 6.25, ) c = c = g () - g() d c Gaussian coe (ρ = 5, T = 1) = 5.000) = = 1.679) = Cysal c = c = 17 c = (min) c = Cysal Gaussian coe model (ρ = 5, T = 1) = 5.000) = = 1.679) = FIG. 15. 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 G. The Hansen-McDonald molen sal model FIG.14. Theeffeconsucueanddynamicsofhecuofffo he LJG liquid. The ed and black cuves give, especively, esuls fo an FCScuoff andalage 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. ISF a he same sae poin. The final aomic sysem sudied is he singly-chaged molen sal model of Hansen and McDonald [56]. In Fig. 16 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. 16).

10 10 g AA () 1.5 c = c = g () - g() d c A singly-chaged molen sal (ρ = 0.37, T = 18) = 9.500) = 0.15 = 2.602) = g() c = c = g () - g() d c Lewis-Wahnsöm OTP (ρ = 0.33, T = 0.70) = 4.500) = 0.91 = 1.469) = A ( q = 4.10, ) A singly-chaged molen sal (ρ = 0.37, T = 18) c = c = 55 c = (min) c = = 9.500) = 0.15 = 2.602) = FIG. 16. 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. ( q = 7.06, ) Lewis-Wahnsöm OTP (ρ = 0.33, T = 0.70) c = c = c = (min) c = = 4.500) = 0.91 = 1.469) = FIG. 17. The effec on sucue and dynamics of vaying he cuoff fo he Wahnsöm OTP model. The ed and black cuves give, especively, esuls fo an FCS cuoff and a lage efeence cuoff. RDF of he LJ paicles 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. ISF 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-Wahnsöm OTP [57, 58] and he asymmeic dumbbell [24] 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 17 and show he LJ paicle 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 his). The OTP model is a bode-line songly coelaing liquid (R = 0.91). Figues 18 and show 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 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). g AA () c = c = g () - g() d c The asymmeic dumbbell (ρ = 0.93, T = 6) = 4.000) = 0.96 = 1.556) =

11 11 A ( q = 7.20, ) The asymmeic dumbbell (ρ = 0.93, T = 6) c = c = c = (min) c = = 4.000) = 0.96 = 1.556) = 0.96 O ( q = 9.23, ) = 4.980) = 8 = 65) = 8 Rigid SPC/E Wae (ρ = 0, T = 4.00) c = c = c = 65 (min) c = FIG. 18. 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 (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 ISF a he same sae poin. FIG. 19. The effec on sucue and dynamics of vaying he cuoff fo he igid SPC/E wae model [59]. The ed and black cuves give, especively, esuls fo an FCS and a lage efeence cuoff. Oxygen-oxygen RDF a ρ = 0 and T = 4.00 (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: Invese powe-law (IPL) sysems (n = 18,6,4), I. Rigid SPC/E wae We conside finally he igid SPC/E wae model [59] (Fig. 19). This model is no songly coelaing a ambien condiions, a fac ha eflecs wae s well-known densiy maximum [19]. 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 = 4.00) 4.0 = 4.980) = 8 = 65) = 8 single-componen Lennad-Jones liquid a densiy ρ = 5, genealized Kob-Andesen binay Lennad-Jones mixues, Wahnsöm binay Lennad-Jones mixue, single-componen Buckingham liquid, Lewis-Wahnsöm OTP model, asymmeic dumbbell model. These sysems ae all songly coelaing [18 22]. 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 showed fuhe ha fo all of he following aomic and molecula sysems sucue and/o dynamics ae no popely epoduced when an FCS cuoff is used: Dzuguov (DZ) liquid, Lennad-Jones Gaussian (LJG) liquid, Gaussian coe model (GCM), Hansen-McDonald molen sal model,

12 12 igid SPC/E wae model. None of hese liquids ae songly coelaing. Fo all hese sysems lage cuoffs ae needed in ode o faihfully epoduce he sysem s physics. In conclusion, a shifed-foces FCS cuoff leads o accuae esuls if and only if he liquid is songly coelaing a he sae poin in quesion. We know of no excepions o his ule. This suggess ha songly coelaing liquids ae chaaceized by he popey ha inemolecula ineacions beyond he FCS can be ignoed. IV. 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 [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 he many simple popeies of songly coelaing liquids (Sec. IVA), hen showing ha his class of liquids can be chaaceized fom hee diffeen pespecives: mahemaically, physically, and chemically (Sec. IV B). This gives hee vey diffeen bu equivalen chaaceizaions, indicaing ha he class of songly coelaing liquids is fundamenal and fuhe moivaing he suggesion ha he eal essence of liquid simpliciy is he exisence of song coelaions of viial/poenialenegy equilibium N V T flucuaions. By connecing o he chemiss concep of non-associaed liquids we hen discuss which eal-wold liquids ae simple (Sec. IV D), discuss biefly liquids nea inefaces (Sec. IV E), and give examples of complex liquid popeies (Sec. IV F). Finally, Sec. IVG poins ou ha ou esuls call ino quesion adiional peubaion heoy, which is based on assuming quie diffeen oles of he aacive and he epulsive foces. A. Songly coelaing liquids and hei simple popeies Mos popeies of songly coelaing liquids follow fom he exisence of isomophs in hei phase diagam (see below). Some of hei popeies wee idenified befoe isomophs wee defined in 2009 [21], howeve, 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[17], agingmay be descibed by adding meely one exa paamee [19, 20], powe-law densiy scaling [60] is obeyed o a good appoximaion, i.e., fo vaying densiy and empeaue he elaxaion ime is a funcion of ρ γ /T [24]. 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 if he following holds [21]: 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 )/kbt1 = C 12 e U((2) 1,...,(2) N )/kbt2. (5) Fo mos sysems he isomoph concep is appoximae jus as WU coelaions ae aely pefec. Thus we do no equie Eq. (5) o be igoously obeyed fo all micoconfiguaions, bu only o a good appoximaion and only 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 v() n has exac isomophs chaaceized by ρ γ /T = Cons. whee γ = n/3. Appendix A of Ref. 21 showed ha a liquid is songly coelaing if and only if i has isomophs o a good appoximaion. This was confimed in Refs. 21 and 22, which showed ha Lennad-Jones ype aomic liquids have good isomophs. Refeence 61 showed ha he songly coelaing Lewis-Wahnsöm OTP and asymmeic dumbbell molecula models also have good isomophs. Equaion (5) has many consequences. These wee deived and discussed in deail in he oiginal isomoph pape fom 2009 (Ref. 21), o which he eade is efeed. Basically, sucue and dynamics a wo isomophic sae poins ae idenical in educed unis. 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 same densiy and empeaue). 2. All N-body enopy ems. Recall ha he excess enopy can be expanded in a seies of wo-body, hee-body, ec, ems. Each em is invaian along an isomoph [21]. 3. The isochoic hea capaciy.

13 13 4. The sucue in educed unis (defined by i ρ 1/3 i fo all paicles i). No only he adial disibuion funcion, bu all highe-ode disibuion funcions in educed unis 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 when jumping fom an equilibaed sae poin o a diffeen sae poin isomophic wih he iniial sae. The absence of elaxaion deives fom he fac ha he Bolzmann pobabiliies of scaled micoconfiguaions ae idenical. Such isomoph jumps wok quie well fo he KABLJ liquid [21], fo he asymmeic dumbbell, and fo he Lewis-Wahnsöm OTP molecula models [61]. Moeove, 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 [16]; his povides ye anohe example of a simple feaue of hese liquids. 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 [21, 22]. This implies invaiance along he coexisence cuve of he educed sucue faco, he educed viscosiy, he educed diffusion consan, ec, aswell aspessueinvaianceofhe meling enopy, and a educed-uni Lindemann meling cieion [21]. Collapse of he wo-ode-paamee maps of Debenedei e al. [62 66] o one-dimensional cuves [21]. Isochonal supeposiion [67], i.e., he fac ha when pessue and empeaue ae vaied, he aveage elaxaion ime deemines he enie elaxaion specum [21]. The above lised popeies of songly coelaing liquids all eflec simple feaues 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 [68] T = f(s)h(ρ). (6) Equaion (6) has a numbe of consequences [68], including he Guneisen equaion of sae and ha he isomophs of LJ liquids in paicula, he LJ solid-liquid coexisence cuve aegiven by (Aρ 4 Bρ 2 )/T = Cons. [69, 70]. B. Mahemaical, physical, and chemical chaaceizaions of songly coelaing liquids Aagivensaepoin, ifheaveagepoenialenegyis denoed by U, he consan-poenial-enegy hypesufaceisdefinedbyω = {( 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. 21 i was shown ha he educed-uni consan-poenial-enegy manifold is invaian along a songly coelaing liquid s isomophs, and ha invaiance cuves exis fo hese manifolds 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. The above povides a mahemaical chaaceizaion of he class of songly coelaing liquids. The physical chaaceizaion 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 (Appendix A of Ref. 21). The chemical chaaceizaion of songly coelaing liquids is he popey documened in he pesen pape: A liquid is songly coelaing a a given sae poin if and only if he liquid s sucue and dynamics ae accuaely calculaed by shifed-foces cuoff simulaions ha ignoe ineacions beyond he fis coodinaion shell. This is an empiical finding fo which we have a pesen no compelling agumens, bu can a pesen meely poin ou wo hings. Fis, he popey of insignificance of ineacions beyond he FCS is an isomoph invaian: If a liquidhasgoodisomophsandifanfcscuoffwokswell a one sae poin, FCS cuoffs mus wok well fo all is isomophic sae poins. Thus he chemical chaaceizaion of songly coelaing liquids is consisen wih he fac ha hese liquids have isomophs. Secondly, almos all of he flucuaions in viial and poenial enegy of he LJ liquid come fom inepaicle sepaaions wihin he FCS [19]. 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 shifed-foce 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.

14 14 C. Defining he class of simple liquids Secion IV A lised songly coelaing liquids many simple popeies. Secion IVB showed ha his liquid class may be chaaceized fom hee quie diffeen poins of view. Clealy 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 NVT ensemble equilibium viial/poenial-enegy coelaion coefficien. Compaed o he sandad definiion of simple liquids (as hose wih adially symmeic pai ineacions) hee ae some noable diffeences: 1. Simpliciy is quanified by a coninuously vaiable, i is no an on/off popey. 2. The degee of simpliciy geneally vaies houghou he phase diagam. Mos songly coelaing liquids loose his popey as he ciical poin is appoached; on he ohe hand many o mos nonsongly coelaing liquids ae expeced o become songly coelaing a vey high pessues [71]. 3. No all aomic liquids (i.e., wih adially symmeic pai ineacions) have simple egions in he low-pessue pa of he phase diagam (compae he Dzuguov, Lennad-Jones Gaussian, Gaussian coe, and molen sal models). 4. No all simple liquids ae aomic (compae he Wahnsöm OTP and he asymmeic dumbbell models). Accoding o he new definiion of liquid simpliciy he case whee he poenial enegy is an Eule homogeneous funcion of he paicle posiions (R = 1) ses he gold sandad fo simpliciy. This is consisen wih he many simple popeies of IPL liquids. 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, moe ealisic songly coelaing liquids like he LJ liquid cease o be so when he liquid-vapo coexisence line is appoached nea he ciical poin, 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 simulaions of he LJ liquid [22] explains seveal peviously noed egulaiies [21, 69]. 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 a a suiably chosen empeaue. 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 [72], bu in simulaions ha he WCALJ liquid has somewha pooe isomophs han he LJ liquid. I is possible 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 could make 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 v() = n v n n. Such an exa analyiciy equiemen would no exclude any songly coelaing liquids occuing in naue whee all poenials ae expeced o be analyic. D. Which liquids in he eal wold ae simple? Real-wold liquids may be classified accoding o he naue of he chemical bonds beween he molecules. Thee ae five ypes of bonds [73], lised below wih a few ypical examples (polymeic sysems may be added as a sepaae class): Van de Waals bonds (agon, oluene, buane,...); Meallic bonds (gold, aluminum, alloys,...); Hydogen bonds (wae, glyceol, ehanol,...); Ionic bonds(molen sodium chloide and poassium niae, oom-empeaue ionic liquids,...); Covalen bonds (silica and boae mels,...). Mos liquids involve elemens of moe han one ype of chemical bonds. Fo insance, van de Waals foces ae pesen in all liquids; he fis class consiss meely of hose liquids ha only have van de Waals foces. Anohe bodeline example is a dipola oganic liquid like di-buyl-phhalae, whee van de Waals as well as Coulomb foces ae pesen; a liquid like glyceol has also song dipola ineacions, i.e., an elemen of he ionic bonds ha defines class 4), ec. Based on compue simulaions and known popeies of liquids we believe ha mos o all van de Waals and meallic liquids ae songly coelaing [14, 17, 19], i.e., simple. Liquids ha ae no simple ae he hydogen-,

15 15 ionically, and covalenly bonding liquids. In hese cases he viial/poenial-enegy coelaions ae weakened by one o he ohe fom of compeing ineacions. Meals play a special ole as simple liquids, because hei ineaomic foces deive fom collecive ineacions beween ion coes and fee elecons [8]. The esuling ineacion is a non-diecional ineacion beween symmeic ion coes, i.e., hese sysems ae simple in he adiional sense. Peliminay compue simulaions of ous show ha meals ae songly coelaing [18], so meals ae simple also in he sense of his pape, as well. Howeve, no all isomoph invaians apply fo meals. Fo insance, he elecon gas can influence he collecive dynamics wihou any visible sucual and elaxaional counepa [74, 75], so isomoph invaiance mos likely beaks down fo hese (fas) collecive degees of feedom. I should be emphasized ha he above consideaions efe o ambien o modeae pessue condiions. I was ecenly suggesed ha if cysallizaion is avoided, all liquids become songly coelaing a sufficienly high pessue [71]. Thus, e.g., he molen silicaes of he eah s uppe manle ae pediced o be simple han molen silicaes a ambien pessue. E. Liquids nea inefaces Since he popey of being simple canno be ead of fom knowledge of he poenial alone, i is of inees o conside liquids unde moe geneal cicumsances, fo insance unde confinemen o geneally nea inefaces. Liquids nea inefaces show ich and complicaed behavio. Fo insance, a liquid confined o he nanoscale may change is dynamic popeies seveal odes of magniude compaed o he bulk sysem. Pedicing hese changes is an impoan challenge elevan fo biological sysems, engineeed devices, ec. Recenly, i was shown ha some liquids eain bulk liquid behavio in confinemen[76 79]. Moe specifically, i was shown ha Rosenfeld s excess enopy scaling in he bulk pesiss in confinemen and is o a good appoximaion independen of he wall-fluid ineacion sengh and he degee of confinemen. This was shown fo LJ and had-sphee liquids, suggesing he possibiliy of exending he concep of a simple liquid o apply beyond bulk sysems. Given ha ineacions in songly coelaing liquids ae limied in angeby he adiusofhe FCS, one mayspeculae ha hese liquids ae simple also by having he popey ha an exenal field a one poin affecs he liquid ove shoe disances han fo a geneal liquid. Clealy, moe wok is needed o claify he elevance and consequences of he pesen definiion of liquid simpliciy nea inefaces and in exenal fields [80, 81]. F. A noe on complex liquid behavio Liquids ha ae no simple in he above defined sense of he em ofen have complex popeies [64, 82 85]. Wae wih is coelaion coefficien close o zeo a ambien condiions is a pime example of a complex liquid. I is well known fo wae ha a ceain egion of sae poins in he densiy/empeaue phase diagam exhibis anomalous hemodynamic behavio in he sense ha isobaic heaing implies densificaion. Numeical evidence indicaes ha hese sae poins lie wihin a lage egion wih diffusion anomaly, i.e., an inceased diffusiviy upon isohemal compession [64], a egion ha in un lies wihin a lage egion of sucual anomaly chaaceized by deceasing ode upon isohemal compession [64]. Diffeen ode paamees exis fo chaaceizing he sucual ode of liquids, some of which elae puely o an inegal ove he RDF [82 84]. In his way i is possible o calculae he conibuion o sucual anomalies fom he diffeen coodinaion shells [82 84]. I has been shown [83, 84] ha he sucual anomaly of wae and waelike liquids is no a:... fis-shell effec. Rahe, hey eflec how sucuing in second and moe disan coodinaion shells esponds o changes in hemodynamic o sysem paamees. The full anomalous behavio of wae deives fom ineacions beyond he FCS [82, 84]. This is consisen wih he esuls pesened in his pape, since he sucue and dynamics of songly coelaing liquids ae given exclusively by he ineacions wihin he FCS. G. To which exen do he assumpions of sandad peubaion heoy hold? Thefinding hahe FCS playsacucialolefoalage class of sysems may be aken as a moden demonsaion of he classic van de Waals picue of liquids in he sense ha such liquids can be undesood in ems of packing effecs [86]. On he ohe hand, ou esuls call ino quesion he basis of adiional peubaion heoy, which is also usually aced back o van de Waals [87]. Peubaion heoy is based on he assumpion of eniely diffeen oles being played by he epulsive and he aacive foces [7, 30, 80, 86 89]: The epulsive foces define he sucue and educe he enopy compaed o ha of an ideal gas a same densiy and empeaue, he aacive foces educe he pessue and enegy compaed o ha of an ideal gas. Fom he findings of he pesen and a pevious pape [32] i is clea, howeve, ha his picue applies only when he FCS coincides oughly wih he pai poenial minimum, i.e., a faily low pessue. A vey high pessues he whole FCS is wihin he ange of he epulsive foces, hee he aacive foces play lile ole fo simple liquids. In geneal, wha is impoan fo a songly coelaing liquid is o ake ino accoun popely all foces fom paicles wihin he FCS and