Small Hard Spheres Research Articles

Temperature, Pressure, and Length-Scale Dependence of Solvation in Water-like Solvents. II. Large Solvophovic Solutes.

We use molecular simulation to determine solvation free energies, isochoric solvation energies and entropies, isobaric solvation enthalpies and entropies, partial molecular volumes, and isothermal density derivatives of the solvation free energy as a function of temperature and pressure for hard-sphere solutes with diameters ranging from 4 to 36 Å in TIP4P/2005 and Jagla water-like solvents exhibiting unusual thermodynamics. An important piece of our discussion focuses on the nanometer-sized solutes, for which simulation results are found to be accounted for by the most basic classical thermodynamic treatment contemplating bulk and interfacial contributions to the solvation free energy. Thus, since water's liquid-vapor surface tension is only special inasmuch as it takes unusually large values, solvent's water-like unusual thermodynamics manifests through a term proportional to the pressure in the solvation free energy. As a result, such solvent's unusual thermodynamics is found to be relevant to the temperature and pressure dependence of the isochoric solvation energy and entropy as well as to the isothermal density derivative of the solvation free energy. This sharply contrasts with the findings of the first part of this series indicating that the solvation free energy of small hard spheres responds to temperature and pressure changes as solvent's density does, with such a contrasting picture embodying a "pressure-density dichotomy." As for the length-scale dependence, we find the zero nominal pressure and the solvent's temperature of the maximum density as singular conditions for cavity surface-area size scaling of large solutes to occur for all solvation quantities. We finally argue that the overall study undertaken in this series suggests that water's unusual thermodynamics may be relevant to the thermodynamic stability of clusters of solvophobic units in the temperature-pressure plane. Some comments on the role of solute-solvent attractive interactions are also depicted.

Effect of particle exchange on the glass transition of binary hard spheres

We investigate the replica theory of the liquid-glass transition for a binary mixture of large and small additive hard spheres. We consider two different ansätze for this problem: the frozen glass ansatz (FGA) in which the exchange of large and small particles in a glass state is prohibited, and the exchange glass ansatz (EGA), in which it is allowed. We calculate the dynamical and thermodynamical glass transition points with the two ansätze. We show that the dynamical transition density of the FGA is lower than that of the EGA, while the thermodynamical transition density of the FGA is higher than that of the EGA. We discuss the algorithmic implications of these results for the density-dependence of the relaxation time of supercooled liquids. We particularly emphasize the difference between the standard Monte Carlo and swap Monte Carlo algorithms. Furthermore, we discuss the importance of particle exchange for estimating the configurational entropy.

Localization switching of a large object in a crowded cavity: A rigid/soft object prefers surface/inner positioning.

For living cells in the real world, a large organelle is commonly positioned in the inner region away from membranes, such as the nucleus of eukaryotic cells, the nucleolus of nuclei, mitochondria, chloroplast, Golgi body, etc. It contradicts the expectation by the current depletion-force theory in that the larger particle should be excluded from the inner cell space onto cell boundaries in a crowding media. Here we simply model a sizable organelle as a soft-boundary large particle allowing crowders, which are smaller hard spheres in the model, to intrude across its boundary. The results of Monte Carlo simulation indicate that the preferential location of the larger particle switches from the periphery into the inner region of the cavity by increasing its softness. An integral equation theory is further developed to account for the structural features of the model, and the theoretical predictions are found consistent with our simulation results.

Assessing the Exact Muffin-Tin Orbitals method for the Bain path of metals

We scrutinise the muffin-tin approximation and the screening within the framework of the Exact Muffin-Tin Orbitals method in the case of cubic and tetragonal crystal symmetries. Systematic total energy calculations are carried out for the Bain path including the body-centred cubic and face-centred cubic structures for a set of simple and transition metals. The present converged results in terms of potential sphere radius (S) and hard sphere radius (b) are in good agreement with previous theoretical calculations. We demonstrate that for all structures considered here, potential sphere radii around and slightly larger than the average Wigner–Seitz radius (w) yield accurate total energy results whereas S values smaller than w give large errors. It is shown that for converged total energies hard spheres with radii b = 0.7–0.8w should be used for an efficient screening within real space clusters consisting typically of 70–90 lattice sites. The less efficient convergence of the total energy in the case of small hard spheres is ascribed to the delocalisation of the screened spherical waves, which leads to inaccurate interstitial overlap matrix. The above conclusions are not significantly affected by the volume of the system.

Equilibrium Grain Boundary Segregation and Clustering of Impurities in Colloidal Polycrystalline Monolayers.

We investigate the segregation of impurities to grain boundaries in colloidal polycrystalline monolayers using video microscopy. A model colloidal alloy is prepared by embedding large spherical impurities in a polycrystalline monolayer of small host colloidal hard spheres, which stops grain growth at a finite grain size. The size ratio between the impurities and the host particles determines whether they behave as interstitial or substitutional impurities in the bulk crystal, akin to those in real alloys. We find that the partitioning of impurities between the grains and the grain boundaries is in very good agreement with the Langmuir-McLean adsorption model for equilibrium grain boundary segregation. This enables the direct measurement of the free energy of adsorption for the two types of impurities. Near saturation, we characterize the spatial distribution of the adsorbed impurities and find that it strongly depends on their interstitial or substitutional nature. This is because the relative importance of clustering and mixing due to nonadditivity is determined by geometrical constraints imposed by the crystalline host lattice.

Generic phase diagram of binary superlattices

Emergence of a large variety of self-assembled superlattices is a dramatic recent trend in the fields of nanoparticle and colloidal sciences. Motivated by this development, we propose a model that combines simplicity with a remarkably rich phase behavior applicable to a wide range of such self-assembled systems. Those systems include nanoparticle and colloidal assemblies driven by DNA-mediated interactions, electrostatics, and possibly, controlled drying. In our model, a binary system of large and small hard spheres (L and S, respectively) interacts via selective short-range ("sticky") attraction. In its simplest version, this binary sticky sphere model features attraction only between S and L particles. We show that, in the limit when this attraction is sufficiently strong compared with [Formula: see text], the problem becomes purely geometrical: the thermodynamically preferred state should maximize the number of [Formula: see text] contacts. A general procedure for constructing the phase diagram as a function of system composition f and particle size ratio r is outlined. In this way, the global phase behavior can be calculated very efficiently for a given set of plausible candidate phases. Furthermore, the geometric nature of the problem enables us to generate those candidate phases through a well-defined and intuitive construction. We calculate the phase diagrams for both 2D and 3D systems and compare the results with existing experiments. Most of the 3D superlattices observed to date are featured in our phase diagram, whereas several more are predicted for future discovery.

Effective stiffness and formation of secondary structures in a protein-like model.

We use Wang-Landau and replica exchange techniques to study the effect of an increasing stiffness on the formation of secondary structures in protein-like systems. Two possible models are considered. In both models, a polymer chain is formed by tethered beads where non-consecutive backbone beads attract each other via a square-well potential representing the tendency of the chain to fold. In addition, smaller hard spheres are attached to each non-terminal backbone bead along the direction normal to the chain to mimic the steric hindrance of side chains in real proteins. The two models, however, differ in the way bending rigidity is enforced. In the first model, partial overlap between consecutive beads is allowed. This reduces the possible bending angle between consecutive bonds thus producing an effective entropic stiffness that competes with a short-range attraction, and leads to the formation of secondary structures characteristic of proteins. We discuss the low-temperature phase diagram as a function of increasing interpenetration and find a transition from a planar, beta-like structure, to helical shape. In the second model, an energetic stiffness is explicitly introduced by imposing an infinitely large energy penalty for bending above a critical angle between consecutive bonds, and no penalty below it. The low-temperature phase of this model does not show any sign of protein-like secondary structures. At intermediate temperatures, however, where the chain is still in the coil conformation but stiffness is significant, we find the two models to predict a quite similar dependence of the persistence length as a function of the stiffness. This behaviour is rationalized in terms of a simple geometrical mapping between the two models. Finally, we discuss the effect of shrinking side chains to zero and find the above mapping to still hold true.

Effective electrostatic interactions among charged thermo-responsive microgels immersed in a simple electrolyte.

This work explores the nature and thermodynamic behavior of the effective electrostatic interactions among charged microgels immersed in a simple electrolyte, taking special interest in the effects due to the thermally induced variation of the microgel size while the remaining parameters (microgel charge and concentration, plus the amount of added salt) are kept constant. To this end, the rigorous approach obtained from applying the precise methodology of the dressed ion theory to the proper definition of the effective direct correlation functions, which emerge from tracing-out the degrees of freedom of the microscopic ions, is employed to provide an exact description of the parameters characterizing such interactions: screening length, effective permittivity, and renormalized charges. A model solution with three components is assumed: large permeable anionic spheres for the microgels, plus small charged hard spheres of equal size for the monovalent cations and anions. The two-body correlations among the components of this model suspension, used as the input for the determination of the effective interaction parameters, are here calculated by using the hyper-netted chain approximation. It is then found that at finite microgel concentrations the values of these parameters change as the microgel size increases, even though the ionic strength of the supporting electrolyte and the bare charge of the microgels remain fixed during this process. The variation of the screening length, as well as that of the effective permittivity, is rather small, but still interesting in view of the fact that the corresponding Debye length stays constant. The renormalized charges, in contrast, increase markedly as the microgels swell. The ratio of the renormalized charge to the corresponding analytic result obtained in the context of an extended linear response theory allows us to introduce an effective charge that accounts for the non-linear effects induced by the short-ranged association of microions to the microgels. The behavior of these effective charges as a function of the amount of added salt and the macroion charge, size, and concentration reveals the interplay among all these system parameters.

Theory and calculations of colloidal depletion interaction

Colloidal suspension is composed of particles with sizes between 1 nm and 1 m, suspended in liquid phase. The interaction between the particles consists of a hard core repulsive interaction and other kinds of repulsive and attractive interacions. Hard interaction forbids the particles from occupying the same places, resulting in a depletion effect. When big colloid particles are immersed in a colloid of small particles, each big particle has a depletion layer where the small particles cannot enter due to the hard interaction. The depletion layers of two big particles overlap when they are close enough so that extra free volume of the small particles increases and therefore the entropy of the small particles increase, thus an effective interaction between big particles is induced. This effective interaction is the so-called depletion interaction. In this review the concepts and an intuitive explanation of depletion interaction of colloidal suspensions are presented. The numerical calculation methods, including the acceptance ratio method, Wang-Landau-type method, and density functional theory method, are briefly reviewed. Several useful analytic approximations are presented. Stating from the depletion interaction between two flat plates, the Derjaguin approximation is introduced through the Asakura- Oosawa model. With this approximation, the approximate formulas of depletion interaction between two hard spheres, between a hard sphere and a hard wall, and between a hard sphere and curved hard walls in a small hard sphere colloid are derived. The depletion interaction between two hard spheres in a thin rod colloid and a thin disk colloid are also derived in the Derjaguin approximation.

Indentation of a rigid sphere into an elastic substrate with surface tension and adhesion.

The surface tension of compliant materials such as gels provides resistance to deformation in addition to and sometimes surpassing that owing to elasticity. This paper studies how surface tension changes the contact mechanics of a small hard sphere indenting a soft elastic substrate. Previous studies have examined the special case where the external load is zero, so contact is driven by adhesion alone. Here, we tackle the much more complicated problem where, in addition to adhesion, deformation is driven by an indentation force. We present an exact solution based on small strain theory. The relation between indentation force (displacement) and contact radius is found to depend on a single dimensionless parameter: ω=σ(μR)-2/3((9π/4)Wad)-1/3, where σ and μ are the surface tension and shear modulus of the substrate, R is the sphere radius and Wad is the interfacial work of adhesion. Our theory reduces to the Johnson-Kendall-Roberts (JKR) theory and Young-Dupre equation in the limits of small and large ω, respectively, and compares well with existing experimental data. Our results show that, although surface tension can significantly affect the indentation force, the magnitude of the pull-off load in the partial wetting liquid-like limit is reduced only by one-third compared with the JKR limit and the pull-off behaviour is completely determined by ω.

Rigidity of a spherical capsule switches the localization of encapsulated particles between inner and peripheral regions under crowding condition: Simple model on cellular architecture

We have investigated the inhomogeneous interior of confined spherical cavities as capsules containing encapsulated binary hard sphere mixtures for different compositions and cavity wall rigidity. Such a greatly simplified model manifests the effects of macromolecular crowding arising from excluded volume interactions in a tiny cell or a cellular nucleus. By fixing the number of large particles, the level of crowding is adjusted by changing the amount of small hard spheres in the cavity. For a rigid cavity, large spheres tend to pack in liquid-like order apart from the surface to the center of the cavity as the crowding level is increased. Whereas, for a soft cavity, larger spheres tend to blend with small spheres in the peripheral region at near the boundary of the cavity, and are susceptible to be depleted from the interior of the cavity as the cavity becomes more crowded. These results may help future elucidation of the thermodynamic pathways to stabilize the inhomogeneous structure of mixtures confined in cavities, such as the derepression of genome materials around the interior rim of the nucleus in a cancerous cell.

Monte Carlo methods for estimating depletion potentials in highly size-asymmetrical hard sphere mixtures

We investigate Monte Carlo simulation strategies for determining the effective ("depletion") potential between a pair of hard spheres immersed in a dense sea of much smaller hard spheres. Two routes to the depletion potential are considered. The first is based on estimates of the insertion probability of one big sphere in the presence of the other; we describe and compare three such methods. The second route exploits collective (cluster) updating to sample the depletion potential as a function of the separation of the big particles; we describe two such methods. For both routes, we find that the sampling efficiency at high densities of small particles can be enhanced considerably by exploiting "geometrical shortcuts" that focus the computational effort on a subset of small particles. All the methods we describe are readily extendable to particles interacting via arbitrary potentials.

Thermodynamically consistent closure approximation for hard spheres systems

We present a new closure relation that is an extension of the Percus-Yevick approximation. In the proposed closure, we introduce an additional term and a mixing coefficient that can be determined by imposing a condition of thermodynamic self-consistency. Moreover, the mixing coefficient is calculated analytically within a linear approximation. In the case of a monodisperse system of hard spheres, we compare the results of our model to well-established thermodynamic expressions and also to the structural properties of fairly known closure approximations. In the second case, and using an equivalent scheme, the new closure relation is extended to the depletion potential between two large hard spheres immersed in a liquid of small hard spheres. In both cases, the results of our model are in good agreement with numerical simulations performed at intermediate concentrations.

Predicting cavity formation free energy: how far is the Gaussian approximation valid?

We examine the range of validity of the Gaussian model for various water-like liquids whose intermolecular potentials differ from SPC/E water, to provide insight into the temperature dependence of the hydrophobic effect for small hard sphere solutes. We find that low compressibility liquids that have more close-packed network structures show much larger deviations from Gaussian fluctuations for low or zero occupancies relative to more compressible fluids with more open networks. Water appears to be a unique molecular fluid in possessing equilibrium density fluctuations that are faithfully described by the Gaussian theory. We ascribe this success to the fact, shown here, that the orientational correlations near a small hard sphere solute involve remarkably little reorganization from the bulk, which is a consequence of water's low solvent reorganization enthalpy and entropy.

Depletion potentials in highly size-asymmetric binary hard-sphere mixtures: Comparison of simulation results with theory

We report a detailed study, using state-of-the-art simulation and theoretical methods, of the effective (depletion) potential between a pair of big hard spheres immersed in a reservoir of much smaller hard spheres, the size disparity being measured by the ratio of diameters q ≡ σ(s)/σ(b). Small particles are treated grand canonically, their influence being parameterized in terms of their packing fraction in the reservoir η(s)(r). Two Monte Carlo simulation schemes--the geometrical cluster algorithm, and staged particle insertion--are deployed to obtain accurate depletion potentials for a number of combinations of q ≤ 0.1 and η(s)(r). After applying corrections for simulation finite-size effects, the depletion potentials are compared with the prediction of new density functional theory (DFT) calculations based on the insertion trick using the Rosenfeld functional and several subsequent modifications. While agreement between the DFT and simulation is generally good, significant discrepancies are evident at the largest reservoir packing fraction accessible to our simulation methods, namely, η(s)(r) = 0.35. These discrepancies are, however, small compared to those between simulation and the much poorer predictions of the Derjaguin approximation at this η(s)(r). The recently proposed morphometric approximation performs better than Derjaguin but is somewhat poorer than DFT for the size ratios and small-sphere packing fractions that we consider. The effective potentials from simulation, DFT, and the morphometric approximation were used to compute the second virial coefficient B(2) as a function of η(s)(r). Comparison of the results enables an assessment of the extent to which DFT can be expected to correctly predict the propensity toward fluid-fluid phase separation in additive binary hard-sphere mixtures with q ≤ 0.1. In all, the new simulation results provide a fully quantitative benchmark for assessing the relative accuracy of theoretical approaches for calculating depletion potentials in highly size-asymmetric mixtures.

Self-Assembly of a Colloidal Interstitial Solid with Tunable Sublattice Doping

We determine the phase diagram of a binary mixture of small and large hard spheres with a size ratio of 0.3 using free-energy calculations in Monte Carlo simulations. We find a stable binary fluid phase, a pure face-centered-cubic (fcc) crystal phase of the small spheres, and binary crystal structures with LS and LS(6) stoichiometries. Surprisingly, we demonstrate theoretically and experimentally the stability of a novel interstitial solid solution in binary hard-sphere mixtures, which is constructed by filling the octahedral holes of an fcc crystal of large spheres with small spheres. We find that the fraction of octahedral holes filled with a small sphere can be completely tuned from 0 to 1. Additionally, we study the hopping of the small spheres between neighboring octahedral holes, and interestingly, we find that the diffusion increases upon increasing the density of small spheres.

Fatigue behavior of foreign object damaged 7075 heat treated aluminum alloy coated with PVD WC/C

The effect of a physically vapor deposited (PVD) WC/C coating on the fatigue behavior of as produced and foreign object damaged (FOD) solution heat treated and aged 7075 aluminum alloy was studied. Coated and uncoated samples were tested under rotating bending to determine the fatigue strengths between 104 and 106 cycles in both damaged and smooth condition. FOD was produced with single shots of small hard steel spheres impacting at 100 m/s in the minimum cross section. SEM was used to characterize the features of the fracture surfaces.

Binary fluid with attractions near a planar wall

It is well known that a mixture of big and small hard spheres next to a planar wall will exhibit segregation based on their size difference. The larger spheres will tend to locate next to the substrate because the overall system entropy loss per unit area is less. In the present study we determine the role of attraction between the small particles and the wall to displace the larger particles. Both fluids density-functional theory and discontinuous molecular dynamics simulations demonstrate that at a certain attractive potential, which is on the order of the thermal energy, the large particles can indeed be dislodged from the surface layer so the small particles are now the major surface component. Exploration of a range of parameters, including relative sphere size and concentration, as well as attractions between the small spheres in the bulk, shows that the phenomenon is quite robust.

Entropic potential field formed for a linear-motor protein near a filament: Statistical-mechanical analyses using simple models

We report a new progress in elucidating the mechanism of the unidirectional movement of a linear-motor protein (e.g., myosin) along a filament (e.g., F-actin). The basic concept emphasized here is that a potential field is entropically formed for the protein on the filament immersed in solvent due to the effect of the translational displacement of solvent molecules. The entropic potential field is strongly dependent on geometric features of the protein and the filament, their overall shapes as well as details of the polyatomic structures. The features and the corresponding field are judiciously adjusted by the binding of adenosine triphosphate (ATP) to the protein, hydrolysis of ATP into adenosine diphosphate (ADP)+Pi, and release of Pi and ADP. As the first step, we propose the following physical picture: The potential field formed along the filament for the protein without the binding of ATP or ADP+Pi to it is largely different from that for the protein with the binding, and the directed movement is realized by repeated switches from one of the fields to the other. To illustrate the picture, we analyze the spatial distribution of the entropic potential between a large solute and a large body using the three-dimensional integral equation theory. The solute is modeled as a large hard sphere. Two model filaments are considered as the body: model 1 is a set of one-dimensionally connected large hard spheres and model 2 is a double helical structure formed by two sets of connected large hard spheres. The solute and the filament are immersed in small hard spheres forming the solvent. The major findings are as follows. The solute is strongly confined within a narrow space in contact with the filament. Within the space there are locations with sharply deep local potential minima along the filament, and the distance between two adjacent locations is equal to the diameter of the large spheres constituting the filament. The potential minima form a ringlike domain in model 1 while they form a pointlike one in model 2. We then examine the effects of geometric features of the solute on the amplitudes and asymmetry of the entropic potential field acting on the solute along the filament. A large aspherical solute with a cleft near the solute-filament interface, which mimics the myosin motor domain, is considered in the examination. Thus, the two fields in our physical picture described above are qualitatively reproduced. The factors to be taken into account in further studies are also discussed.

Phase Transitions in a Mixture of Hard-Sphere Dipoles and Neutral Hard Spheres

Using the reference hypernetted chain (RHNC) integral equation theory and a rigorous stability analysis method, we investigate the phase behavior of a mixture of hard-sphere dipoles and neutral hard spheres based on the correlations of the homogeneous isotropic phase. Lowering the temperature down to the points where the RHNC equations fail to have a solution, several fluctuations strongly increase. At low densities our results indicate the onset of chain formation, which is similar with the pure DHS system. At high densities, the results indicate the appearance of isotropic-to-ferroelectric transitions (small neutral hard spheres concentrations) and demixing transitions (large neutral hard spheres concentrations).