Accurate simulation estimates of phase behaviour in ternary mixtures with prescribed composition

This paper describes an isobaric semi-grand canonical ensemble Monte Carlo scheme for the accurate study of phase behaviour in ternary fluid mixtures under the experimentally relevant conditions of prescribed pressure, temperature and overall composition. It is shown how to tune the relative chemical potentials of the individual components to target some requisite overall composition and how, in regions of phase coexistence, to extract accurate estimates for the compositions and phase fractions of individual coexisting phases. The method is illustrated by tracking a path through the composition space of a model ternary Lennard-Jones mixture.Comment: 6 pages, 3 figure

The fluid-fluid interface in a model colloid-polymer mixture: Application of grand canonical Monte Carlo to asymmetric binary mixtures

We present a Monte Carlo method to simulate asymmetric binary mixtures in the grand canonical ensemble. The method is used to study the colloid-polymer model of Asakura and Oosawa. We determine the phase diagram of the fluid-fluid unmixing transition and the interfacial tension, both at high polymer density and close to the critical point. We also present density profiles in the two-phase region. The results are compared to predictions of a recent density functional theory.Comment: 4 pages, 4 figure

Chimeric glutamate receptor subunits reveal the transmembrane domain is sufficient for NMDA receptor pore properties but some positive allosteric modulators require additional domains

NMDA receptors are ligand-gated ion channels that underlie transmission at excitatory synapses and play an important role in regulating synaptic strength and stability. Functional NMDA receptors require two copies of the GluN1 subunit coassembled with GluN2 (and/or GluN3) subunits into a heteromeric tetramer. A diverse array of allosteric modulators can upregulate or downregulate NMDA receptor activity. These modulators include both synthetic compounds and endogenous modulators, such as cis-unsaturated fatty acids, 24(S)-hydroxycholesterol, and various neurosteroids. To evaluate the structural requirements for the formation and allosteric modulation of NMDA receptor pores, we have replaced portions of the rat GluN1, GluN2A, and GluN2B subunits with homologous segments from the rat GluK2 kainate receptor subunit. Our results with these chimeric constructs show that the NMDA receptor transmembrane domain is sufficient to account for most pore properties, but that regulation by some allosteric modulators requires additional cytoplasmic or extracellular domains. SIGNIFICANCE STATEMENT Glutamate receptors mediate excitatory synaptic transmission by forming cation channels through the membrane that open upon glutamate binding. Although many compounds have been identified that regulate glutamate receptor activity, in most cases the detailed mechanisms that underlie modulation are poorly understood. To identify what parts of the receptor are essential for pore formation and sensitivity to allosteric modulators, we generated chimeric subunits that combined segments from NMDA and kainate receptors, subtypes with distinct pharmacological profiles. Surprisingly, our results identify separate domain requirements for allosteric potentiation of NMDA receptor pores by pregnenolone sulfate, 24(S)-hydroxycholesterol, and docosahexaenoic acid, three endogenous modulators derived from membrane constituents. Understanding where and how these compounds act on NMDA receptors should aid in designing better therapeutic agents

Development of a GIS-based seismic hazard screening tool

Geotechnical earthquake engineering hazards have consequences that are difficult to mitigate, especially for infrastructure systems with multiple and distributed components. In the last few decades significant progress has been made to provide more accurate and useful methods to evaluate hazards for complex systems. One of these advances involves the application of spatial analysis and geographic information systems (GIS), for not only presenting data as maps, but also providing more practical and usable solutions, such as calculating the hazard potential with spatial distribution. This thesis studies the evaluation of geotechnical earthquake engineering hazards within a GIS environment, using borehole-specific data and seismic ground motions. Existing methods and applications used to evaluate these hazards, as well as existing geotechnical database formats are presented and discussed. This research developed a GIS methodology to be used as a screening tool: to evaluate geotechnical earthquake engineering hazards from a database of borehole data and then display the results on a map --Abstract, page iii

Three-body interactions in complex fluids: virial coefficients from simulation finite-size effects

A simulation technique is described for quantifying the contribution of three-body interactions to the thermodynamical properties of coarse-grained representations of complex fluids. The method is based on comparing the third virial coefficient $B_3$ for a complex fluid with that of an approximate coarse-grained model described by a pair potential. To obtain $B_3$ we introduce a new technique which expresses its value in terms of the measured volume-dependent asymptote of a certain structural function. The strategy is applicable to both Molecular Dynamics and Monte Carlo simulation. Its utility is illustrated via measurements of three-body effects in models of star polymer and highly size-asymmetrical colloid-polymer mixtures.Comment: 13 pages, 8 figure

Liquid-gas coexistence and critical point shifts in size-disperse fluids

Specialized Monte Carlo simulations and the moment free energy (MFE) method are employed to study liquid-gas phase equilibria in size-disperse fluids. The investigation is made subject to the constraint of fixed polydispersity, i.e. the form of the `parent' density distribution $\rho^0(\sigma)$ of the particle diameters $\sigma$, is prescribed. This is the experimentally realistic scenario for e.g. colloidal dispersions. The simulations are used to obtain the cloud and shadow curve properties of a Lennard-Jones fluid having diameters distributed according to a Schulz form with a large (40%) degree of polydispersity. Good qualitative accord is found with the results from a MFE method study of a corresponding van der Waals model that incorporates size-dispersity both in the hard core reference and the attractive parts of the free energy. The results show that polydispersity engenders considerable broadening of the coexistence region between the cloud curves. The principal effect of fractionation in this region is a common overall scaling of the particle sizes and typical inter-particle distances, and we discuss why this effect is rather specific to systems with Schulz diameter distributions. Next, by studying a family of such systems with distributions of various widths, we estimate the dependence of the critical point parameters on $\delta$. In contrast to a previous theoretical prediction, size-dispersity is found to raise the critical temperature above its monodisperse value. Unusually for a polydisperse system, the critical point is found to lie at or very close to the extremum of the coexistence region in all cases. We outline an argument showing that such behaviour will occur whenever size polydispersity affects only the range, rather than the strength of the inter-particle interactions.Comment: 14 pages, 12 figure

Liquid-vapor interface of a polydisperse fluid

We report a Grand Canonical Monte Carlo simulation study of the liquid-vapor interface of a model fluid exhibiting polydispersity in terms of the particle size $\sigma$. The bulk density distribution, $\rho^0(\sigma)$, of the system is controlled by the imposed chemical potential distribution $\mu(\sigma)$. We choose the latter such that $\rho^0(\sigma)$ assumes a Schulz form with associated degree of polydispersity $\approx 14$. By introducing a smooth attractive wall, a planar liquid-vapor interface is formed for bulk state points within the region of liquid-vapor coexistence. Owing to fractionation, the pure liquid phase is enriched in large particles, with respect to the coexisting vapor. We investigate how the spatial non-uniformity of the density near the liquid-vapor interface affects the evolution of the local distribution of particle sizes between the limiting pure phase forms. We find (as previously predicted by density functional theory, Bellier-Castella {\em et al}, Phys. Rev. {\bf E65}, 021503 (2002)) a segregation of smaller particles to the interface. The magnitude of this effect is quantified for various $\sigma$ via measurements of the relative adsorption. Additionally, we consider the utility of various estimators for the interfacial width and highlight the difficulties of isolating the intrinsic contribution of polydispersity to this width.Comment: 9 pages, 10 Fig