Surface tension of electrolytes: Hydrophilic and hydrophobic ions near an interface

We calculate the ion distributions around an interface in fluid mixtures of highly polar and less polar fluids (water and oil) for two and three ion species. We take into account the solvation and image interactions between ions and solvent. We show that hydrophilic and hydrophobic ions tend to undergo a microphase separation at an interface, giving rise to an enlarged electric double layer. We also derive a general expression for the surface tension of electrolyte systems, which contains a negative electrostatic contribution proportional to the square root of the bulk salt density. The amplitude of this square-root term is small for hydrophilic ion pairs, but is much increased for hydrophilic and hydrophobic ion pairs. For three ion species including hydrophilic and hydrophobic ions, we calculate the ion distributions to explain those obtained by x-ray reflectivity measurements.Comment: 8 figure

Solvent-free coarse-grained lipid model for large-scale simulations

A coarse-grained molecular model, which consists of a spherical particle and an orientation vector, is proposed to simulate lipid membrane on a large length scale. The solvent is implicitly represented by an effective attractive interaction between particles. A bilayer structure is formed by orientation-dependent (tilt and bending) potentials. In this model, the membrane properties (bending rigidity, line tension of membrane edge, area compression modulus, lateral diffusion coefficient, and flip-flop rate) can be varied over broad ranges. The stability of the bilayer membrane is investigated via droplet-vesicle transition. The rupture of the bilayer and worm-like micelle formation can be induced by an increase in the spontaneous curvature of the monolayer membrane.Comment: 13 pages, 19 figure

Role of the first coordination shell in determining the equilibrium structure and dynamics of simple liquids

The traditional view that the physical properties of a simple liquid are determined primarily by its repulsive forces was recently challenged by Berthier and Tarjus, who showed that in some cases ignoring the attractions leads to large errors in the dynamics [L. Berthier and G. Tarjus, Phys. Rev. Lett. 103, 170601 (2009); J. Chem. Phys. 134, 214503 (2011)]. We present simulations of the standard Lennard-Jones liquid at several condensed-fluid state points, including a fairly low density state and a very high density state, as well as simulations of the Kob-Andersen binary Lennard-Jones mixture at several temperatures. By varying the range of the forces, results for the thermodynamics, dynamics, and structure show that the determining factor for getting the correct statics and dynamics is not whether or not the attractive forces {\it per se} are included in the simulations. What matters is whether or not interactions are included from all particles within the first coordination shell (FCS) - the attractive forces can thus be ignored, but only at extremely high densities. The recognition of the importance of a local shell in condensed fluids goes back to van der Waals; our results confirm this idea and thereby the basic picture of the old hole- and cell theories for simple condensed fluids

Description of the fluctuating colloid-polymer interface

To describe the full spectrum of surface fluctuations of the interface between phase-separated colloid-polymer mixtures from low scattering vector q (classical capillary wave theory) to high q (bulk-like fluctuations), one must take account of the interface's bending rigidity. We find that the bending rigidity is negative and that on approach to the critical point it vanishes proportionally to the interfacial tension. Both features are in agreement with Monte Carlo simulations.Comment: 5 pages, 3 figures, 1 table. Accepted for publication in Phys. Rev. Let

Capillary Waves in a Colloid-Polymer Interface

The structure and the statistical fluctuations of interfaces between coexisting phases in the Asakura-Oosawa (AO) model for a colloid--polymer mixture are analyzed by extensive Monte Carlo simulations. We make use of a recently developed grand canonical cluster move with an additional constraint stabilizing the existence of two interfaces in the (rectangular) box that is simulated. Choosing very large systems, of size LxLxD with L=60 and D=120, measured in units of the colloid radius, the spectrum of capillary wave-type interfacial excitations is analyzed in detail. The local position of the interface is defined in terms of a (local) Gibbs surface concept. For small wavevectors capillary wave theory is verified quantitatively, while for larger wavevectors pronounced deviations show up. For wavevectors that correspond to the typical distance between colloids in the colloid-rich phase, the interfacial fluctuations exhibit the same structure as observed in the bulk structure factor. When one analyzes the data in terms of the concept of a wavevector-dependent interfacial tension, a monotonous decrease of this quantity with increasing wavevector is found. Limitations of our analysis are critically discussed.Comment: 12 pages, 15 figure

Viscous coalescence of droplets: a Lattice Boltzmann study

The coalescence of two resting liquid droplets in a saturated vapor phase is investigated by Lattice Boltzmann simulations in two and three dimensions. We find that, in the viscous regime, the bridge radius obeys a t^{1/2}-scaling law in time with the characteristic time scale given by the viscous time. Our results differ significantly from the predictions of existing analytical theories of viscous coalescence as well as from experimental observations. While the underlying reason for these deviations is presently unknown, a simple scaling argument is given that describes our results well.Comment: 12 pages, 10 figures; as published in Phys. Fluid

Isotropic-nematic interfacial tension of hard and soft rods: application of advanced grand canonical biased sampling techniques

Coexistence between the isotropic and the nematic phase in suspensions of rods is studied using grand canonical Monte Carlo simulations with a bias on the nematic order parameter. The biasing scheme makes it possible to estimate the interfacial tension gamma in systems of hard and soft rods. For hard rods with L/D=15, we obtain gamma ~ 1.4 kB T/L^2, with L the rod length, D the rod diameter, T the temperature, and kB the Boltzmann constant. This estimate is in good agreement with theoretical predictions, and the order of magnitude is consistent with experiments.Comment: 10 pages, 10 figure

Institutional determinants of construction safety management strategies of contractors in Hong Kong

From an institutional perspective, organisations are not only a production system; but also a social and cultural system. The external institutional environment in which organisations are embedded plays a critical role in shaping organisational structures and practices. Organisations as active agencies respond strategically to the institutional pressures and expectations in a variety of ways, such as conformity, compromise and avoidance. Building on the conceptual framework established by Oliver (1991), this study explores how contractors respond to institutional demands in terms of compliance with construction health and safety requirements, involvement in voluntary safety initiatives and safety campaigns, and commitment to zero accidents in the context of the Hong Kong construction industry. A case study was conducted to collect empirical evidence of contractors' strategic responses. The results indicate that the contractor prefers to adopt compromise and avoidance strategies in the face of the client's stringent safety requirements and tight project progress demands. This study contributes to the knowledge of safety management by offering an institutional explanation of contractors' safety management strategies

Solid-Liquid Phase Diagrams for Binary Metallic Alloys: Adjustable Interatomic Potentials

We develop a new approach to determining LJ-EAM potentials for alloys and use these to determine the solid-liquid phase diagrams for binary metallic alloys using Kofke's Gibbs-Duhem integration technique combined with semigrand canonical Monte Carlo simulations. We demonstrate that it is possible to produce a wide-range of experimentally observed binary phase diagrams (with no intermetallic phases) by reference to the atomic sizes and cohesive energies of the two elemental materials. In some cases, it is useful to employ a single adjustable parameter to adjust the phase diagram (we provided a good choice for this free parameter). Next, we perform a systematic investigation of the effect of relative atomic sizes and cohesive energies of the elements on the binary phase diagrams. We then show that this approach leads to good agreement with several experimental binary phase diagrams. The main benefit of this approach is not the accurately reproduction of experimental phase diagrams, but rather to provide a method by which material properties can be continuously changed in simulations studies. This is one of the keys to the use of atomistic simulations to understand mechanisms and properties in a manner not available to experiment

Interfacial tension of the isotropic--nematic interface in suspensions of soft spherocylinders

The isotropic to nematic transition in a system of soft spherocylinders is studied by means of grand canonical Monte Carlo simulations. The probability distribution of the particle density is used to determine the coexistence density of the isotropic and the nematic phases. The distributions are also used to compute the interfacial tension of the isotropic--nematic interface, including an analysis of finite size effects. Our results confirm that the Onsager limit is not recovered until for very large elongation, exceeding at least L/D=40, with L the spherocylinder length and D the diameter. For smaller elongation, we find that the interfacial tension increases with increasing L/D, in agreement with theoretical predictions.Comment: 8 pages, 7 figures, and also 1 tabl