Search

A generic but simple model is presented to evaluate the effect of the heterogeneous surface charge distribution of proteins and zwitterionic nanoparticles on their thermodynamic phase behavior. By considering surface charges as continuous "patches," the rich set of surface patterns that is embedded in proteins and charged patchy particles can readily be described. This model is used to study the fluid phase separation of charged particles where the screening length is of the same order of magnitude as the particle size. In particular, two types of charged particles are studied: dipolar fluids and proteinlike fluids. The former represents the simplest case of zwitterionic particles, whose charge distribution can be described by their dipole moment. The latter system corresponds to molecules/particles with complex surface charge arrangements such as those found in biomolecules. The results for both systems suggest a relation between the critical region, the strength of the interparticle interactions, and the arrangement of charged patches, where the critical temperature is strongly correlated to the magnitude of the dipole moment. Additionally, competition between attractive and repulsive charge-charge interactions seems to be related to the formation of fluctuating clusters in the dilute phase of dipolar fluids, as well as to the broadening of the binodal curve in protein-like fluids. Finally, a variety of self-assembled architectures are detected for dipolar fluids upon small changes to the charge distribution, providing the groundwork for studying the self-assembly of charged patchy particles. [ABSTRACT FROM AUTHOR]

Cubic colloids, sedimented on a surface and immersed in a solution of depletant molecules, were modeled with a family of shapes which smoothly varies from squares to circles. UsingWang-Landau simulations with expanded ensembles, we observe the formation of rhombic lattices, square lattices, hexagonal lattices, and a fluid phase. This systematic investigation includes locating transitions between all combinations of the three lattice structures upon changing the shape and transitions between the fluid and crystal upon changing the depletant concentration. The rhombic lattice deforms smoothly between square-like and hexagonal-like angles, depending on both the shape and the depletant concentration. Our results on the effect of the depletant concentration, depletant size, and colloid shape to influence the stability of the fluid and the lattice structures may help guide experimental studies with recently synthesized cubic colloids. [ABSTRACT FROM AUTHOR]

Here we review recent work by the authors to revisit the concept of extrapolating thermodynamic properties of classical systems using statistical mechanical principles. Specifically, we discuss how the combination of these principles with biased sampling techniques enables the prediction of free energy landscapes and other detailed information, such as structural properties, of the system in question. Remarkably accurate estimates of physical properties across a broad range of conditions have been achieved using this approach, greatly reducing the number of simulations needed to explore a given system's behaviour. While approximate, these extrapolations significantly amplify the amount of reasonably accurate information that can be extracted from simulations enabling a small set of them to feed data-intensive regression algorithms such as neural networks. Thus, this extrapolation methodology represents a useful tool for performing tasks such as high-throughput screening of physical properties, optimising force field parameters, exploring equilibrium phase behaviour, and enabling theory-guided data science for these systems. [ABSTRACT FROM AUTHOR]

Using flat-histogram Monte Carlo methods, we investigate the adsorptive behavior of the square-well fluid in two simple slit-pore-like models intended to capture fundamental characteristics of flexible adsorbent materials. Both models require as input thermodynamic information about the flexible adsorbent material itself. An important component of this work involves formulating the flexible pore models in the appropriate thermodynamic (statistical mechanical) ensembles, namely, the osmotic ensemble and a variant of the grand-canonical ensemble. Two-dimensional probability distributions, which are calculated using flat-histogram methods, provide the information necessary to determine adsorption thermodynamics. For example, we are able to determine precisely adsorption isotherms, (equilibrium) phase transition conditions, limits of stability, and free energies for a number of different flexible adsorbent materials, distinguishable as different inputs into the models. While the models used in this work are relatively simple from a geometric perspective, they yield non-trivial adsorptive behavior, including adsorption-desorption hysteresis solely due to material flexibility and so-called "breathing" of the adsorbent. The observed effects can in turn be tied to the inherent properties of the bare adsorbent. Some of the effects are expected on physical grounds while others arise from a subtle balance of thermodynamic and mechanical driving forces. In addition, the computational strategy presented here can be easily applied to more complex models for flexible adsorbents. [ABSTRACT FROM AUTHOR]

We have calculated the pair correlation functions for several binary mixtures composed of simple solutes in a Lennard-Jones solvent. In particular, we have studied the solute-solute pair correlation functions and their dependence on the total density, the solvent Lennard-Jones parameters, and on the solute-solute energy parameter. All the results were obtained from solving the Percus–Yevick equations, as well as from Monte Carlo simulations. The relevance of these results to the problem of hydrophobic interactions is also discussed. [ABSTRACT FROM AUTHOR]

Computer simulations are employed to obtain subcritical isotherms of small finite sized systems inside the coexistence region. For all temperatures considered, ranging from the triple point up to the critical point, the isotherms gradually developed a sequence of sharp discontinuities as the system size increased from ∼8 to ∼21 molecular diameters. For the smallest system sizes, and more so close to the critical point, the isotherms appeared smooth, resembling the continuous van der Waals loop obtained from extrapolation of an analytic equation of state outside the coexistence region. As the system size was increased, isotherms in the chemical potential-density plane developed first two, then four, and finally six discontinuities. Visual inspection of selected snapshots revealed that the observed discontinuities are related to structural transitions between droplets (on the vapor side) and bubbles (on the liquid side) of spherical, cylindrical, and tetragonal shapes. A capillary drop model was developed to qualitatively rationalize these observations. Analytic results were obtained and found to be in full agreement with the computer simulation results. The analysis shows that the shape of the subcritical isotherms is dictated by a single characteristic volume (or length scale), which depends on the surface tension, compressibility, and coexistence densities. For small reduced system volumes, the model predicts that a homogeneous fluid is stable across the whole coexistence region, thus explaining the continuous van der Waals isotherms observed in the simulations. When the liquid and vapor free energies are described by means of an accurate mean-field equation of state and surface tensions from simulation are employed, the capillary model is found to describe the simulated isotherms accurately, especially for large systems (i.e., larger than about 15 molecular diameters) at low temperature (lower than about 0.85 times the critical temperature). This implies that the Laplace pressure differences can be predicted for drops as small as five molecular diameters, and as few as about 500 molecules. The theoretical study also shows that the extrema or apparent spinodal points of the finite size loops are more closely related to (finite system size) bubble and dew points than to classical spinodals. Our results are of relevance to phase transitions in nanopores and show that first order corrections to nucleation energies in finite closed systems are power laws of the inverse volume. [ABSTRACT FROM AUTHOR]

We present a methodology based on grand-canonical transition-matrix Monte Carlo and finite-size scaling analysis to calculate surface tensions in binary mixtures. In particular, mixture transition-matrix Monte Carlo is first used to calculate apparent, system-size-dependent free-energy barriers separating coexisting fluid phases. Finite-size scaling is then used to extrapolate these values to the infinitely large system limit to determine the true thermodynamic surface tension. A key distinction of the methodology is that it yields the entire isothermal surface-tension curve for a binary mixture in a relatively small number of simulations. We demonstrate the utility of the method by calculating surface-tension curves for three binary Lennard-Jones mixtures. While we have only examined the surface tension of simple fluids in this work, the method is general and can be extended to molecular fluids as well as to determine interfacial tensions of liquid-liquid interfaces. [ABSTRACT FROM AUTHOR]

We present a method for directly locating density-driven phase transitions in multicomponent systems. Phase coexistence conditions are determined through manipulation of a total density probability distribution evaluated over a density range that includes both coexisting phases. Saturation quantities are determined through appropriate averaging of density-dependent mean values of a given property of interest. We discuss how to implement the method in both the grand-canonical and isothermal-isobaric semigrand ensembles. Calculations can be conducted using any of the recently introduced flat-histogram techniques. Here, we combine the general algorithm with a transition-matrix approach to produce an efficient self-adaptive technique for determining multicomponent phase equilibrium properties. To assess the performance of the new method, we generate phase diagrams for a number of binary and ternary Lennard-Jones mixtures. [ABSTRACT FROM AUTHOR]

We present a novel computational methodology for determining fluid-phase equilibria in binary mixtures. The method is based on a combination of highly efficient transition-matrix Monte Carlo and histogram reweighting. In particular, a directed grand-canonical transition-matrix Monte Carlo scheme is used to calculate the particle-number probability distribution, after which histogram reweighting is used as a postprocessing procedure to determine the conditions of phase equilibria. To validate the methodology, we have applied it to a number of model binary Lennard-Jones systems known to exhibit nontrivial fluid-phase behavior. Although we have focused on monatomic fluids in this work, the method presented here is general and can be easily extended to more complex molecular fluids. Finally, an important feature of this method is the capability to predict the entire fluid-phase diagram of a binary mixture at fixed temperature in a single simulation. [ABSTRACT FROM AUTHOR]

We present a kinetic theory of homogeneous bubble nucleation based on explicit calculation of the single-molecule evaporation and condensation rates as a function of the size of the vapor embryo. The surface condensation rate is calculated from the kinetic theory of gases, and the surface evaporation rate is related to the rate of escape of molecules from a potential well in the field established by the liquid-vapor interface. Equality of these rates corresponds naturally to the critical bubble. While the interface plays a crucial role in this respect, the kinetic nucleation theory does not invoke an explicit surface tension. The nucleation rate is derived from a population balance and depends only on the ratio of the evaporation to condensation rates. In contrast to classical theory, a nontrivial trend captured by the present theory is the increase in nucleation rate with decreasing temperature at fixed degree of metastability. Comparison with classical nucleation theory reveals markedly different supersaturation dependencies of the nucleation rate, while the predicted sizes of the critical bubble are in good agreement. [ABSTRACT FROM AUTHOR]

Results from a comparative study of recently introduced tail corrections are presented to calculate surface tension in systems where the dominant long-range interaction is dispersive. It is found that with standard tail corrections the surface tension decreases with increasing cutoff distance reaching a limiting value corresponding to the maximum cutoff possible.

Grand canonical transition-matrix Monte Carlo is used to investigate kinetic and thermodynamic stability limits of the liquid and vapor in the Lennard-Jones fluid. It is found that the properties of the supersaturated vapor at its limit of stability approach their saturation values with increasing system size.

Reports on a computational study of homogeneous liquid-vapor nucleation in the Lennard-Jones fluid. Use of the umbrella sampling Monte Carlo simulations to calculate free energy barriers to homogeneous liquid-vapor nucleation; Decrease in calculated free energy barriers with increased superheating.

The residual multiparticle entropy (RMPE) of two- and three-dimensional fluids changes sign near the freezing line, providing a quasiuniversal “one-phase” rule for the location of the liquid-solid transition. We present new simulation results for d-dimensional hard-sphere fluids (d=1-5) which show, however, that this freezing criterion fails in other spatial dimensions. The results also call into question the idea that a change in sign of the RMPE implies the emergence of a new kind of local structural order in the fluid. [ABSTRACT FROM AUTHOR]