Towards a quantification of disorder in materials: distinguishing equilibrium and glassy sphere packings

Robust algorithm to generate a diverse class of dense disordered and ordered sphere packings via linear programming

We have formulated the problem of generating dense packings of nonoverlapping, nontiling nonspherical particles within an adaptive fundamental cell subject to periodic boundary conditions as an optimization problem called the adaptive-shrinking cell (ASC) formulation [S. Torquato and Y. Jiao, Phys. Rev. E 80, 041104 (2009)]. Because the objective function and impenetrability constraints can be exactly linearized for sphere packings with a size distribution in d-dimensional Euclidean space R(d), it is most suitable and natural to solve the corresponding ASC optimization problem using sequential-linear-programming (SLP) techniques. We implement an SLP solution to produce robustly a wide spectrum of jammed sphere packings in R(d) for d=2, 3, 4, 5, and 6 with a diversity of disorder and densities up to the respective maximal densities. A novel feature of this deterministic algorithm is that it can produce a broad range of inherent structures (locally maximally dense and mechanically stable packings), besides the usual disordered ones (such as the maximally random jammed state), with very small computational cost compared to that of the best known packing algorithms by tuning the radius of the influence sphere. For example, in three dimensions, we show that it can produce with high probability a variety of strictly jammed packings with a packing density anywhere in the wide range [0.6, 0.7408...], where π/√18 = 0.7408... corresponds to the density of the densest packing. We also apply the algorithm to generate various disordered packings as well as the maximally dense packings for d=2, 4, 5, and 6. Our jammed sphere packings are characterized and compared to the corresponding packings generated by the well-known Lubachevsky-Stillinger (LS) molecular-dynamics packing algorithm. Compared to the LS procedure, our SLP protocol is able to ensure that the final packings are truly jammed, produces disordered jammed packings with anomalously low densities, and is appreciably more robust and computationally faster at generating maximally dense packings, especially as the space dimension increases.

Solidification and ordering during directional drying of a colloidal dispersion

During drying, colloidal dispersions undergo processes such as solidification, cracking, and the draining of interstitial pores. Here we show that the solidification of polystyrene and silica dispersions, during directional drying, occurs in two separate stages. These correspond to the initial ordering and subsequent aggregation of the colloidal particles. Transitions between these stages are observed as changes in transparency and color that propagate as distinct fronts along the drying layer. The dynamics of these fronts are shown to arise from a balance between compressive capillary forces and the electrostatic and van der Waals forces described by DLVO theory. This suggests a simple method by which the maximum interparticle repulsion between particles can be measured through the optical inspection of the dynamics of a drying dispersion, under a microscope.

Dense packing in the monodisperse hard-sphere system: a numerical study

We report a numerical study of the close packing of monodisperse hard spheres. The close packings of hard spheres are produced by the Lubachesky-Stillinger (LS) compression algorithm and span the range from the disordered states to the ordered states. We provide quantitative evidence for the claim that the density and structural order of the arrested close packing can be determined by the compression rate, i.e., with slower rates producing denser and more ordered structures. Through deeply analyzing the structure of the resulting arrested close packings, a transition region has been identified in the plane of density and reciprocal compression rate, in between what have been historically thought of as amorphous and crystalline packings. We also find clear system size dependences in studying the structural properties of the packings from the disordered ones to the ordered ones. These detailed investigations, on the structure of the arrested close packings, may provide a link between the glassy states and the crystalline states in the hard spheres.

Pressure-energy correlations in liquids. IV. "Isomorphs" in liquid phase diagrams

This paper is the fourth in a series devoted to identifying and explaining the properties of strongly correlating liquids, i.e., liquids where virial and potential energy correlate better than 90% in their thermal equilibrium fluctuations in the NVT ensemble. For such liquids we here introduce the concept of "isomorphic" curves in the phase diagram. A number of thermodynamic, static, and dynamic isomorph invariants are identified. These include the excess entropy, the isochoric specific heat, reduced-unit static and dynamic correlation functions, as well as reduced-unit transport coefficients. The dynamic invariants apply for both Newtonian and Brownian dynamics. It is shown that after a jump between isomorphic state points the system is instantaneously in thermal equilibrium; consequences of this for generic aging experiments are discussed. Selected isomorph predictions are validated by computer simulations of the Kob-Andersen binary Lennard-Jones mixture, which is a strongly correlating liquid. The final section of the paper relates the isomorph concept to phenomenological melting rules, Rosenfeld's excess entropy scaling, Young and Andersen's approximate scaling principle, and the two-order parameter maps of Debenedetti and co-workers. This section also shows how the existence of isomorphs implies an "isomorph filter" for theories for the non-Arrhenius temperature dependence of viscous liquids' relaxation time, and it explains isochronal superposition for strongly correlating viscous liquids.

Geometrical ambiguity of pair statistics: point configurations

Point configurations have been widely used as model systems in condensed-matter physics, materials science, and biology. Statistical descriptors, such as the n -body distribution function g(n), are usually employed to characterize point configurations, among which the most extensively used is the pair distribution function g(2). An intriguing inverse problem of practical importance that has been receiving considerable attention is the degree to which a point configuration can be reconstructed from the pair distribution function of a target configuration. Although it is known that the pair-distance information contained in g(2) is, in general, insufficient to uniquely determine a point configuration, this concept does not seem to be widely appreciated and general claims of uniqueness of the reconstructions using pair information have been made based on numerical studies. In this paper, we present the idea of the distance space called the D space. The pair distances of a specific point configuration are then represented by a single point in the D space. We derive the conditions on the pair distances that can be associated with a point configuration, which are equivalent to the realizability conditions of the pair distribution function g(2). Moreover, we derive the conditions on the pair distances that can be assembled into distinct configurations, i.e., with structural degeneracy. These conditions define a bounded region in the D space. By explicitly constructing a variety of degenerate point configurations using the D space, we show that pair information is indeed insufficient to uniquely determine the configuration in general. We also discuss several important problems in statistical physics based on the D space, including the reconstruction of atomic structures from experimentally obtained g(2) and a recently proposed "decorrelation" principle. The degenerate configurations have relevance to open questions involving the famous traveling salesman problem.

A tetrahedral entropy for water

We introduce the space-dependent correlation function C (Q)(r) and time-dependent autocorrelation function C (Q)(t) of the local tetrahedral order parameter Q identical with Q(r,t). By using computer simulations of 512 waterlike particles interacting through the transferable interaction potential with five points (TIP5 potential), we investigate C (Q)(r) in a broad region of the phase diagram. We find that at low temperatures C (Q)(t) exhibits a two-step time-dependent decay similar to the self-intermediate scattering function and that the corresponding correlation time tau(Q) displays a dynamic cross-over from non-Arrhenius behavior for T > T (W) to Arrhenius behavior for T < T (W), where T (W) denotes the Widom temperature where the correlation length has a maximum as T is decreased along a constant-pressure path. We define a tetrahedral entropy S (Q) associated with the local tetrahedral order of water molecules and find that it produces a major contribution to the specific heat maximum at the Widom line. Finally, we show that tau(Q) can be extracted from S (Q) by using an analog of the Adam-Gibbs relation.

Using compressibility factor as a predictor of confined hard-sphere fluid dynamics

We study the correlations between the diffusivity (or viscosity) and the compressibility factor of bulk hard-sphere fluid as predicted by the ultralocal limit of the barrier hopping theory. Our specific aim is to determine if these correlations observed in the bulk equilibrium hard-sphere fluid can be used to predict the self-diffusivity of fluid confined between a slit-pore or a rectangular channel. In this work, we consider a single-component and a binary mixture of hard spheres. To represent confining walls, we use purely reflecting hard walls and interacting square-well walls. Our results clearly show that the correspondence between the diffusivity and the compressibility factor can be used along with the knowledge of the confined fluid's compressibility factor to predict its diffusivity with quantitative accuracy. Our analysis also suggests that a simple measure, the average fluid density, can be an accurate predictor of confined fluid diffusivity for very tight confinements ( approximately 2-3 particle diameters wide) at low to intermediate density conditions. Together, these results provide further support for the idea that one can use robust connections between thermodynamic and dynamic quantities to predict dynamics of confined fluids from their thermodynamics.

The characteristic crystallographic element norm: A descriptor of local structure in atomistic and particulate systems

We introduce the characteristic crystallographic element (CCE) norm as a powerful descriptor of local structure in atomistic and particulate systems. The CCE-norm is sensitive both to radial and orientational deviations from perfect local order. Unlike other measures of local order, the CCE-norm decreases monotonically with increasing order, is zero for a perfectly ordered environment, and is strictly discriminating among different, competing crystal structures in imperfectly ordered systems. The CCE-norm descriptor can be used as a sensitive, quantitative measure to detect and track changes in local order in atomistic and general particulate systems. In a specific example we show the ability of the CCE-norm to monitor the onset and evolution of order in an initially amorphous, densely packed assembly of hard-sphere chains generated through extensive Monte Carlo simulations [Phys. Rev. Lett. 100, 050602 (2008)].

Dense sphere packings from optimized correlation functions

Elementary smooth functions (beyond contact) are employed to construct pair correlation functions that mimic jammed disordered sphere packings. Using the g_{2} -invariant optimization method of Torquato and Stillinger [J. Phys. Chem. B 106, 8354 (2002)], parameters in these functions are optimized under necessary realizability conditions to maximize the packing fraction varphi and average number of contacts per sphere Z . A pair correlation function that incorporates the salient features of a disordered packing and that is smooth beyond contact is shown to permit a varphi of 0.6850: this value represents a 45% reduction in the difference between the maximum for congruent hard spheres in three dimensions, pi/sqrt[18] approximately 0.7405 and 0.64, the approximate fraction associated with maximally random jammed packings in three dimensions. We show that, surprisingly, the continued addition of elementary functions consisting of smooth sinusoids decaying as r;{-4} permits packing fractions approaching pi/sqrt[18] . A translational order metric is used to discriminate between degrees of order in the packings presented. We find that to achieve higher packing fractions, the degree of order must increase, which is consistent with the results of a previous study [Torquato, Phys. Rev. Lett. 84, 2064 (2000)].

Anomalous structure and dynamics of the Gaussian-core fluid

It is known that there are thermodynamic states for which the Gaussian-core fluid displays anomalous properties such as expansion upon isobaric cooling (density anomaly) and increased single-particle mobility upon isothermal compression (self-diffusivity anomaly). Here, we investigate how temperature and density affect its short-range translational structural order, as characterized by the two-body excess entropy. We find that there is a wide range of conditions for which the short-range translational order of the Gaussian-core fluid decreases upon isothermal compression (structural order anomaly). As we show, the origin of the structural anomaly is qualitatively similar to that of other anomalous fluids (e.g., water or colloids with short-range attractions) and is connected to how compression affects static correlations at different length scales. Interestingly, we find that the self-diffusivity of the Gaussian-core fluid obeys a scaling relationship with the two-body excess entropy that is very similar to the one observed for a variety of simple liquids. One consequence of this relationship is that the state points for which structural, self-diffusivity, and density anomalies of the Gaussian-core fluid occur appear as cascading regions on the temperature-density plane; a phenomenon observed earlier for models of waterlike fluids. There are, however, key differences between the anomalies of Gaussian-core and waterlike fluids, and we discuss how those can be qualitatively understood by considering the respective interparticle potentials of these models. Finally, we note that the self-diffusivity of the Gaussian-core fluid obeys different scaling laws depending on whether the two-body or total excess entropy is considered. This finding, which deserves more comprehensive future study, appears to underscore the significance of higher-body correlations for the behavior of fluids with bounded interactions.

Anomalous waterlike behavior in spherically-symmetric water models optimized with the relative entropy

Recent efforts have attempted to understand many of liquid water's anomalous properties in terms of effective spherically-symmetric pairwise molecular interactions entailing two characteristic length scales (so-called "core-softened" potentials). In this work, we examine the extent to which such simple descriptions of water are representative of the true underlying interactions by extracting coarse-grained potential functions that are optimized to reproduce the behavior of an all-atom model. To perform this optimization, we use a novel procedure based upon minimizing the relative entropy, a quantity that measures the extent to which a coarse-grained configurational ensemble overlaps with a reference all-atom one. We show that the optimized spherically-symmetric water models exhibit notable variations with the state conditions at which they were optimized, reflecting in particular the shifting accessibility of networked hydrogen bonding interactions. Moreover, we find that water's density and diffusivity anomalies are only reproduced when the effective coarse-grained potentials are allowed to vary with state. Our results therefore suggest that no state-independent spherically-symmetric potential can fully capture the interactions responsible for water's unique behavior; rather, the particular way in which the effective interactions vary with temperature and density contributes significantly to anomalous properties.

Relation of water anomalies to the excess entropy

Using the five-site transferable intermolecular potential (TIP5P) we perform molecular dynamics simulations to investigate the relationship between the excess entropy and anomalies of water. We find that the two-body excess entropy is an ideal quantity to predict the regions of structural, dynamic, and thermodynamic anomalies of water in its pressure-temperature and density-temperature phase diagrams. From the excess entropy we can also predict the location of the Widom line, associated with the hypothesized liquid-liquid critical point in supercooled water.

Structural and thermodynamic properties of different phases of supercooled liquid water

Computer simulation results are reported for a realistic polarizable potential model of water in the supercooled region. Three states, corresponding to the low density amorphous ice, high density amorphous ice, and very high density amorphous ice phases are chosen for the analyses. These states are located close to the liquid-liquid coexistence lines already shown to exist for the considered model. Thermodynamic and structural quantities are calculated, in order to characterize the properties of the three phases. The results point out the increasing relevance of the interstitial neighbors, which clearly appear in going from the low to the very high density amorphous phases. The interstitial neighbors are found to be, at the same time, also distant neighbors along the hydrogen bonded network of the molecules. The role of these interstitial neighbors has been discussed in connection with the interpretation of recent neutron scattering measurements. The structural properties of the systems are characterized by looking at the angular distribution of neighboring molecules, volume and face area distribution of the Voronoi polyhedra, and order parameters. The cumulative analysis of all the corresponding results confirms the assumption that a close similarity between the structural arrangement of molecules in the three explored amorphous phases and that of the ice polymorphs I(h), III, and VI exists.

Shear-rate-dependent structural order and viscosity of a fluid with short-range attractions

We study a model short-range attractive fluid under shear. For this system, the strength of interparticle attractions strongly influences the equilibrium structural order. We find that shear monotonically decreases structural order regardless of the strength of the attractions. There is a strong correlation between shear-rate-dependent viscosity and a structural order metric, suggesting a structurally based constitutive equation. This correlation also holds for the Lennard-Jones fluid.

A computational investigation of thermodynamics, structure, dynamics and solvation behavior in modified water models

We investigate the properties of geometrically modified water models by performing molecular dynamics simulations of perturbations of the extended simple point charge (SPC/E) model of water over a wide range of temperatures at 1 bar. The geometric modification consists of altering the H-O-H angle in SPC/E. The dipole moment is held constant by altering the O-H bond length, while the electrostatic charges are left unchanged. We find that a H-O-H angle of at least 100 degrees is necessary for the appearance of density anomalies and of solubility extrema with respect to temperature for small apolar solutes. We observe the occurrence of two incompatible types of structural order in these models: Tetrahedral, with waterlike translational order for bent models with H-O-H angles in excess of 100 degrees ; and linear, with Lennard-Jones-like orientationally averaged translational order for smaller H-O-H angles. Increasing the H-O-H angle causes the density to increase, while at the same time shifting waterlike anomalies to progressively higher temperatures. For bent models with H-O-H angle greater than SPC/E's, we observe arrest of translational motion at 300 K (115 degrees) and 330 K (120 degrees).

Waterlike hierarchy of anomalies in a continuous spherical shouldered potential

We investigate by molecular dynamics simulations a continuous isotropic core-softened potential with attractive well in three dimensions, introduced by Franzese [J. Mol. Liq. 136, 267 (2007)], that displays liquid-liquid coexistence with a critical point and waterlike density anomaly. Besides the thermodynamic anomalies, here we find diffusion and structural anomalies. The anomalies, not observed in the discrete version of this model, occur with the same hierarchy that characterizes water. We discuss the differences in the anomalous behavior of the continuous and the discrete model in the framework of the excess entropy, calculated within the pair correlation approximation.

Use of the cage formation probability for obtaining approximate phase diagrams

In this work, we introduce the idea of cage formation probability, defined by considering the angular space needed by a particle in order to leave a cage given an average distance to its neighbors. Considering extreme fluctuations, two phases appear as a function of the number of neighbors and their distances to a central one: Solid and fluid. This allows us to construct an approximated phase diagram based on a geometrical approach. As an example, we apply this probability concept to hard disks in two dimensions and hard spheres in three dimensions. The results are compared with numerical simulations using a Monte Carlo method.

Structural anomalies of fluids: origins in second and higher coordination shells

Compressing or cooling a fluid typically enhances its static interparticle correlations. However, there are notable exceptions. Isothermal compression can reduce the translational order of fluids that exhibit anomalous waterlike trends in their thermodynamic and transport properties, while isochoric cooling (or strengthening of attractive interactions) can have a similar effect on fluids of particles with short-range attractions. Recent simulation studies by Yan [Phys. Rev. E 76, 051201 (2007)] on the former type of system and Krekelberg [J. Chem. Phys. 127, 044502 (2007)] on the latter provide examples where such structural anomalies can be related to specific changes in second and more distant coordination shells of the radial distribution function. Here, we confirm the generality of this microscopic picture through analysis, via molecular simulation and integral equation theory, of coordination shell contributions to the two-body excess entropy for several related model fluids which incorporate different levels of molecular resolution. The results suggest that integral equation theory can be an effective and computationally inexpensive tool for assessing, based on the pair potential alone, whether new model systems are good candidates for exhibiting structural (and hence thermodynamic and transport) anomalies.

Tuning density profiles and mobility of inhomogeneous fluids

Density profiles are the most common measure of inhomogeneous structure in confined fluids, but their connection to transport coefficients is poorly understood. We explore via simulation how tuning particle-wall interactions to flatten or enhance the particle layering of a model confined fluid impacts its self-diffusivity, viscosity, and entropy. Interestingly, interactions that eliminate particle layering significantly reduce confined fluid mobility, whereas those that enhance layering can have the opposite effect. Excess entropy helps to understand and predict these trends.

Influence of adhesion and friction on the geometry of packings of spherical particles

We study the effect of both adhesion and friction on the geometry of monosized packings of spheres by means of discrete element simulations. We use elastic properties that are characteristic of materials typically used for particulate processing (Young's modulus in the range 20-200 GPa). The geometrical features, both global and local, of the packings are studied using a variety of approaches in order to investigate their ability to quantify the effect of adhesion and/or friction. We show that both adhesion and friction interaction decrease the packing fraction. The very localized ordering that adhesion triggers is particularly investigated by use of the radial distribution function, the ordering parameter Q(6) , and four triclinic cells that allow a description of the microstructure at the local level. We show that the probability of occurrence of these triclinic cells is approximately proportional to their degree of freedom when neither adhesion nor friction plays a role. We find that the introduction of adhesive interactions increases the probability of occurrence of those cells that have the lowest degree of freedom.

Structure of the first- and second-neighbor shells of simulated water: quantitative relation to translational and orientational order

We perform molecular dynamics simulations of water using the five-site transferable interaction potential (TIP5P) model to quantify structural order in both the first shell (defined by four nearest neighbors) and second shell (defined by twelve next-nearest neighbors) of a central water molecule. We find that the anomalous decrease of orientational order upon compression occurs in both shells, but the anomalous decrease of translational order upon compression occurs mainly in the second shell. The decreases of translational order and orientational order upon compression (called the "structural anomaly") are thus correlated only in the second shell. Our findings quantitatively confirm the qualitative idea that the thermodynamic, structural, and hence dynamic anomalies of water are related to changes upon compression in the second shell.

How short-range attractions impact the structural order, self-diffusivity, and viscosity of a fluid

We present molecular simulation data for viscosity, self-diffusivity, and the local structural ordering of (i) a hard-sphere fluid and (ii) a square-well fluid with short-range attractions. The latter fluid exhibits a region of dynamic anomalies in its phase diagram, where its mobility increases upon isochoric cooling, which is found to be a subset of a larger region of structural anomalies, in which its pair correlations strengthen upon isochoric heating. This "cascade of anomalies" qualitatively resembles that found in recent simulations of liquid water. The results for the hard-sphere and square-well systems also show that the breakdown of the Stokes-Einstein relation upon supercooling occurs for conditions where viscosity and self-diffusivity develop different couplings to the degree of pairwise structural ordering of the liquid. We discuss how these couplings reflect dynamic heterogeneities. Finally, we note that the simulation data suggest how repulsive and attractive glasses may generally be characterized by two distinct levels of short-range structural order.

Entropy, local order, and the freezing transition in Morse liquids

The behavior of the excess entropy of Morse and Lennard-Jones liquids is examined as a function of temperature, density, and the structural order metrics. The dominant pair correlation contribution to the excess entropy is estimated from simulation data for the radial distribution function. The pair correlation entropy (S2) of these simple liquids is shown to have a threshold value of (-3.5+/-0.3)kB at freezing. Moreover, S2 shows a T(-2/5) temperature dependence. The temperature dependence of the pair correlation entropy as well as the behavior at freezing closely correspond to earlier predictions, based on density functional theory, for the excess entropy of repulsive inverse power and Yukawa potentials [Rosenfeld, Phys. Rev. E 62, 7524 (2000)]. The correlation between the pair correlation entropy and the local translational and bond orientational order parameters is examined, and, in the case of the bond orientational order, is shown to be sensitive to the definition of the nearest neighbors. The order map between translational and bond orientational order for Morse liquids and solids shows a very similar pattern to that seen in Lennard-Jones-type systems.

Does confining the hard-sphere fluid between hard walls change its average properties?

We use grand canonical transition-matrix Monte Carlo and discontinuous molecular dynamics simulations to generate precise thermodynamic and kinetic data for the equilibrium hard-sphere fluid confined between smooth hard walls. These simulations show that the pronounced inhomogeneous structuring of the fluid normal to the confining walls, often the primary focus of density functional theory studies, has a negligible effect on many of its average properties over a surprisingly broad range of conditions. We present one consequence of this insensitivity to confinement: a simple analytical equation relating the average density of the confined fluid to that of the bulk fluid with equal activity. Nontrivial implications of confinement for average fluid properties do emerge in this system, but only when the fluid is both (i) dense and (ii) confined to a gap smaller than approximately three particle diameters. For this limited set of conditions, we find that "in-phase" oscillatory deviations in excess entropy and self-diffusivity (relative to the behavior of the bulk fluid at the same average density) occur as a function of gap size. These paired thermodynamic/kinetic deviations from bulk behavior appear to reflect the geometric packing frustration that arises when the confined space cannot naturally accommodate an integer number of particle layers.

Underconstrained jammed packings of nonspherical hard particles: ellipses and ellipsoids

Continuing on recent computational and experimental work on jammed packings of hard ellipsoids [Donev, Science 303, 990 (2004)] we consider jamming in packings of smooth strictly convex nonspherical hard particles. We explain why an isocounting conjecture, which states that for large disordered jammed packings the average contact number per particle is twice the number of degrees of freedom per particle (Z[over]=2d{f}) , does not apply to nonspherical particles. We develop first- and second-order conditions for jamming and demonstrate that packings of nonspherical particles can be jammed even though they are underconstrained (hypoconstrained, Z[over]<2d{f}). We apply an algorithm using these conditions to computer-generated hypoconstrained ellipsoid and ellipse packings and demonstrate that our algorithm does produce jammed packings, even close to the sphere point. We also consider packings that are nearly jammed and draw connections to packings of deformable (but stiff) particles. Finally, we consider the jamming conditions for nearly spherical particles and explain quantitatively the behavior we observe in the vicinity of the sphere point.

Entropy, diffusivity, and structural order in liquids with waterlike anomalies

The excess entropy, defined as the difference between the entropies of the liquid and the ideal gas under identical density and temperature conditions, is studied as a function of density and temperature for liquid silica and a two-scale ramp potential, both of which are known to possess waterlike liquid state anomalies. The excess entropy for both systems is evaluated using a fairly accurate pair correlation approximation. The connection between the excess entropy and the density and diffusional anomalies is demonstrated. Using the pair correlation approximation to the excess entropy, it can be shown that if the energetically favorable local geometries in the low and high density limits have different symmetries, then a structurally anomalous regime can be defined in terms of orientational and translational order parameters, as in the case of silica and the two-scale ramp system but not for the one-scale ramp liquid. Within the category of liquids with waterlike anomalies, we show that the relationship between the macroscopic entropy and internal energy is sufficient to distinguish between those with local anisotropy and consequent open packings at low densities and those with isotropic interactions but multiple length scales. Since it is straightforward to evaluate the pair correlation entropy and internal energy from simulations or experimental data, such plots should provide a convenient means to diagnose the existence as well as type of anomalous behavior in a range of liquids, including ionic and intermetallic melts and complex fluids with ultrasoft repulsions.

Excess-entropy-based anomalies for a waterlike fluid

Many thermodynamic and dynamic properties of water display unusual behavior at low enough temperatures. In a recent study, Yan et al. [Phys. Rev. Lett. 95, 130604 (2005)] identified a spherically symmetric two-scale potential that displays many of the same anomalous properties as water. More specifically, for select parametrizations of the potential, one finds that the regions where isothermal compression anomalously (i) decreases the fluid's structural order, (ii) increases its translational self-diffusivity, and (iii) increases its entropy form nested domes in the temperature-density plane. These property relationships are similar to those found for more realistic models of water. In this work, the authors provide evidence that suggests that the anomalous regions specified above can all be linked through knowledge of the excess entropy. Specifically, the authors show how entropy scaling relationships developed by Rosenfeld [Phys. Rev. A 15, 2545 (1977)] can be used to describe the region of diffusivity anomalies and to predict the state conditions for which anomalous viscosity and thermal conductivity behavior might be found.

Supercooling of aqueous NaCl and KCl solutions under acoustic levitation

The supercooling capability of aqueous NaCl and KCl solutions is investigated at containerless state by using acoustic levitation method. The supercooling of water is obviously enhanced by the alkali metal ions and increases linearly with the augmentation of concentrations. Furthermore, the supercooling depends on the nature of ions and is 2-3 K larger for NaCl solution than that for KCl solution in the present concentration range: Molecular dynamics simulations are performed to reveal the intrinsic correlation between supercoolability and microstructure. The translational and orientational order parameters are applied to quantitatively demonstrate the effect of ionic concentration on the hydrogen-bond network and ice melting point. The disrupted hydrogen-bond structure determines essentially the concentration dependence of supercooling. On the other hand, the introduced acoustic pressure suppresses the increase of supercooling by promoting the growth and coalescence of microbubbles, the effective nucleation catalysts, in water. However, the dissolved ions can weaken this effect, and moreover the degree varies with the ion type. This results in the different supercoolability for NaCl and KCl solutions under the acoustic levitation conditions.

Structural anomalies for a three dimensional isotropic core-softened potential

Using molecular dynamics simulations we investigate the structure of a system of particles interacting through a continuous core-softened interparticle potential. We found for the translational order parameter t a local maximum at a density rho(t-max) and a local minimum at rho(t-min)>rho(t-max). Between rho(t-max) and rho(t-min), the t parameter anomalously decreases upon increasing pressure. For the orientational order parameter Q(6) a maximum was observed at a density rho(t-max)<rho(Qmax)<rho(t-min). For densities between rho(Qmax) and rho(t-min), both the translational (t) and orientational (Q(6)) order parameters have anomalous behavior. We know that this system also exhibits density and diffusion anomalies. We found that the region in the pressure-temperature phase diagram of the structural anomaly englobes the region of the diffusion anomaly that is larger than the region limited by the temperature of maximum density. This cascade of anomalies (structural, dynamic, and thermodynamic) for our model has the same hierarchy as that observed for the simple point charge/extended water.

How does column packing microstructure affect column efficiency in liquid chromatography?

Full three-dimensional computer simulations of the fluid flow and dispersion characteristics of model nonporous chromatographic packings are reported. Interstitial porosity and packing defects are varied in an attempt to understand the chromatographic consequences of the packing microstructure. The tracer zone dispersion is calculated in the form of plate height as a function of fluid velocity for seven model particle packs where particles are selectively removed from the packs in clusters of varying size and topology. In an attempt to examine the consequences of loose but random packs, the velocities and zone dispersion of seven defect-free packs are simulated over the range 0.36< or =epsilon< or =0.50, where epsilon is the interstitial porosity. The results indicate that defect-free loose packings can give good chromatographic efficiency but the efficiency can vary depending on subtle details of the pack. When the defect population increases, the zone dispersion increases accordingly. For a particle pack where 6% of the particles are removed from an epsilon=0.36 pack, approximately 33% of the column efficiency is lost. These results show that it is far more important in column packing to prevent defect sites leading to inhomogeneous packing rather than obtaining the highest density pack with the smallest interstitial void volume.

Attraction-driven disorder in a hard-core colloidal monolayer

Monte Carlo simulation techniques were employed to explore the effect of short-range attraction on the orientational ordering in a two-dimensional assembly of monodisperse spherical particles. We find that if the range of square-well attraction is approximately 15% of the particle diameter, the dense attractive fluid shows the same ordering behavior as the same density fluid composed of purely repulsive hard spheres. Fluids with an attraction range larger than 15% show an enhanced tendency to crystallization, while disorder occurs for fluids with an attractive range shorter than 15% of the particle diameter. A possible link with the existence of "repulsive" and "attractive" states in dense colloidal systems is discussed.

Family of tunable spherically symmetric potentials that span the range from hard spheres to waterlike behavior

We investigate the equation of state, diffusion coefficient, and structural order of a family of spherically symmetric potentials consisting of a hard core and a linear repulsive ramp. This generic potential has two characteristic length scales: the hard and soft core diameters. The family of potentials is generated by varying their ratio, lambda. We find negative thermal expansion (thermodynamic anomaly) and an increase of the diffusion coefficient upon isothermal compression (dynamic anomaly) for 0< or =lambda<6/7. As in water, the regions where these anomalies occur are nested domes in the (T, p) or (T, P) planes , with the thermodynamic anomaly dome contained entirely within the dynamic anomaly dome. We calculate translational and orientational order parameters (t and Q6), and project equilibrium state points onto the (t, Q6) plane, or order map. The order map evolves from waterlike behavior to hard-sphere-like behavior upon varying lambda between 4/7 and 6/7. Thus, we traverse the range of liquid behavior encompassed by hard spheres (lanbda=1) and waterlike (lambda approximately 4/7) with a family of tunable spherically symmetric potentials by simply varying the ratio of hard to soft-core diameters. Although dynamic and thermodynamic anomalies occur almost across the entire range 0< or=lambda< or=1, waterlike structural anomalies (i.e., decrease in both t and Q6 upon compression and strictly correlated T and Q6 in the anomalous region) occur only around lambda=4/7. Waterlike anomalies in structure, dynamics and thermodynamics arise solely due to the existence of two length scales, with their ratio lambda being the single control parameter, orientation-dependent interactions being absent by design.

Entropy-based measure of structural order in water

We analyze the nature of the structural order established in liquid TIP4P water in the framework provided by the multiparticle correlation expansion of the statistical entropy. Different regimes are mapped onto the phase diagram of the model upon resolving the pair entropy into its translational and orientational components. These parameters are used to quantify the relative amounts of positional and angular order in a given thermodynamic state, thus allowing a structurally unbiased definition of low-density and high-density water. As a result, the structurally anomalous region--within which both types of order are simultaneously disrupted by an increase of pressure at constant temperature--is clearly identified through extensive molecular-dynamics simulations.

Coarse-grained strategy for modeling protein stability in concentrated solutions. II: phase behavior

We use highly efficient transition-matrix Monte Carlo simulations to determine equilibrium unfolding curves and fluid phase boundaries for solutions of coarse-grained globular proteins. The model we analyze derives the intrinsic stability of the native state and protein-protein interactions from basic information about protein sequence using heteropolymer collapse theory. It predicts that solutions of low hydrophobicity proteins generally exhibit a single liquid phase near their midpoint temperatures for unfolding, while solutions of proteins with high sequence hydrophobicity display the type of temperature-inverted, liquid-liquid transition associated with aggregation processes of proteins and other amphiphilic molecules. The phase transition occurring in solutions of the most hydrophobic protein we study extends below the unfolding curve, creating an immiscibility gap between a dilute, mostly native phase and a concentrated, mostly denatured phase. The results are qualitatively consistent with the solution behavior of hemoglobin (HbA) and its sickle variant (HbS), and they suggest that a liquid-liquid transition resulting in significant protein denaturation should generally be expected on the phase diagram of high-hydrophobicity protein solutions. The concentration fluctuations associated with this transition could be a driving force for the nonnative aggregation that can occur below the midpoint temperature.