Search

A statistical mechanical theory is presented to predict the effects of macromolecular crowding on protein association equilibria, accounting for both excluded volume and attractive interactions between proteins and crowding molecules. Predicted binding free energies are in excellent agreement with simulation data over a wide range of crowder sizes and packing fraction. It is shown that attractive interactions between proteins and crowding agents counteract the stabilizing effects of excluded volume interactions. A critical attraction strength, for which there is no net effect of crowding, is almost independent of the crowder packing fraction., Comment: 5 pages, 5 figures, 2 tables, to appear in Physical Review Letters

Grand canonical simulations designed to resolve critical universality classes are reported for $z$:1 hard-core electrolyte models with diameter ratios $\lambda {=} a_+/a_- {\lesssim} 6$. For $z {=} 1$ Ising-type behavior prevails. Unbiased estimates of $T_c(\lambda)$ are within 1% of previous (biased) estimates but the critical densities are $\sim $5 % lower. Ising character is also established for the 2:1 and 3:1 equisized models, along with critical amplitudes and improved $T_c$ estimates. For $z {=} 3$, however, strong finite-size effects reduce the confidence level although classical and O$(n {\geq} 3)$ criticality are excluded., Comment: 4 pages, 3 figures

Precise calculations of the coexistence-curve diameters of a hard-core square-we ll (HCSW) fluid and the restricted primitive model (RPM) electrolyte exhibit mar ked deviations from rectilinear behavior. The HCSW diameter displays a $|t|^{1- alpha}$ singularity that sets in sharply for $|t|\equiv |T-T_c|/T_c\lesssim 10^{-3}$; this compares favorably with extensive data for ${SF}_6$, also reflec ted in C$_2$H$_6$, N$_2$, etc. By contrast, the curvature of the RPM diameter va ries slowly over a wide range $|t|\lesssim 0.1$; this behavior mirrors observati ons for liquid alkali metals, specifically Rb and Cs. Amplitudes for the leading singular terms can be estimated numerically but their values cannot be taken li terally., Comment: 9 pages and 4 figures

Simulations of the Lebowitz length, $\xi_{\text{L}}(T,\rho)$, are reported for t he restricted primitive model hard-core (diameter $a$) 1:1 electrolyte for densi ties $\rho\lesssim 4\rho_c$ and $T_c \lesssim T \lesssim 40T_c$. Finite-size eff ects are elucidated for the charge fluctuations in various subdomains that serve to evaluate $\xi_{\text{L}}$. On extrapolation to the bulk limit for $T\gtrsim 10T_c$ the low-density expansions (Bekiranov and Fisher, 1998) are seen to fail badly when $\rho > {1/10}\rho_c$ (with $\rho_c a^3 \simeq 0.08$). At highe r densities $\xi_{\text{L}}$ rises above the Debye length, $\xi_{\text{D}} \prop to \sqrt{T/\rho}$, by 10-30% (upto $\rho\simeq 1.3\rho_c$); the variation is portrayed fairly well by generalized Debye-H\"{u}ckel theory (Lee and Fisher, 19 96). On approaching criticality at fixed $\rho$ or fixed $T$, $\xi_{\text{L}}(T, \rho)$ remains finite with $\xi_{\text{L}}^c \simeq 0.30 a \simeq 1.3 \xi_{\text {D}}^c$ but displays a weak entropy-like singularity., Comment: 4 pages 5 figures

Individual processive molecular motors, of which conventional kinesin is the most studied quantitatively, move along polar molecular tracks and, by exerting a force ${\bm F} = (F_x,F_y,F_z)$ on a tether, drag cellular cargoes, {\em in vivo}, or spherical beads, {\em in vitro}, taking up to hundreds of nanometer-scale steps. From observations of velocities and the dispersion of displacements with time, under measured forces and controlled fuel supply (typically ATP), one may hope to obtain insight into the molecular motions undergone in the individual steps. In the simplest situation, the load force ${\bm F}$ may be regarded as a scalar resisting force, $F_x < 0$, acting parallel to the track: however, experiments, originally by Gittes {\em et al.} (1996), have imposed perpendicular (or vertical) loads, $F_z > 0$, while more recently Block and coworkers (2002, 2003) and Carter and Cross (2005) have studied {\em assisting} (or reverse) loads, $F_x > 0$, and also sideways (or transverse) loads $F_y \neq 0$.

A general method for estimating the Yang-Yang ratio, ${\cal R}_{\mu}$, and the coexistence-curve diameter of a model fluid via Monte Carlo simulations is presented on the basis of data for a hard-core square-well (HCSW) fluid and the restricted primitive model (RPM) electrolyte. The isothermal minima of $Q_{L}\equiv< m^{2}>^{2}_{L}/< m^{4}>_{L}$ are evaluated at $T_{c}$ in an $L\times L\times L$ box where $m = \rho - <\rho>_{L}$ is the density fluctuation. The ``complete'' finite-size scaling theory for the $Q_{\scriptsize min}^{\pm}(T_{c};L)$ incorporates pressure mixing in the scaling fields, thereby allowing for a Yang-Yang anomaly., Comment: Submitted for publication in Phys. Rev. E

In simulating continuum model fluids that undergo phase separation and criticality, significant gains in computational efficiency may be had by confining the particles to the sites of a lattice of sufficiently fine spacing, $a_{0}$ (relative to the particle size, say $a$). But a cardinal question, investigated here, then arises, namely: How does the choice of the lattice discretization parameter, $\zeta\equiv a/a_{0}$, affect the values of interesting parameters, specifically, critical temperature and density, $T_{\scriptsize c}$ and $\rho_{\scriptsize c}$? Indeed, for small $\zeta (\lesssim 4 $-$ 8)$ the underlying lattice can strongly influence the thermodynamic properties. A heuristic argument, essentially exact in $d=1$ and $d=2$ dimensions, indicates that for models with hard-core potentials, both $T_{\scriptsize c}(\zeta)$ and $\rho_{\scriptsize c}(\zeta)$ should converge to their continuum limits as $1/\zeta^{(d+1)/2}$ for $d\leq 3$ when $\zeta\to\infty$; but the behavior of the error is highly erratic for $d\geq 2$. For smoother interaction potentials, the convergence is faster. Exact results for $d=1$ models of van der Waals character confirm this; however, an optimal choice of $\zeta$ can improve the rate of convergence by a factor $1/\zeta$. For $d\geq 2$ models, the convergence of the {\em second virial coefficients} to their continuum limits likewise exhibit erratic behavior which is seen to transfer similarly to $T_{\scriptsize c}$ and $\rho_{\scriptsize c}$; but this can be used in various ways to enhance convergence and improve extrapolation to $\zeta = \infty$ as is illustrated using data for the restricted primitive model electrolyte., Comment: To appear in J. Phys. Chem. in honor of David Chandler

A novel algorithm is presented that yields precise estimates of coexisting liquid and gas densities, $\rho^{\pm}(T)$, from grand canonical Monte Carlo simulations of model fluids near criticality. The algorithm utilizes data for the isothermal minima of the moment ratio $Q_{L}(T;<\rho>_{L})$ $\equiv< m^{2}>_{L}^{2}/< m^{4}>_{L}$ in $L$$ \times$$ ...$$ \times$$ L$ boxes, where $m=\rho-<\rho>_{L}$. When $L$$ \to$$ \infty$ the minima, $Q_{\scriptsize m}^{\pm}(T;L)$, tend to zero while their locations, $\rho_{\scriptsize m}^{\pm}(T;L)$, approach $\rho^{+}(T)$ and $\rho^{-}(T)$. Finite-size scaling relates the ratio {\boldmath $\mathcal Y$}$ = $$(\rho_{\scriptsize m}^{+}-\rho_{\scriptsize m}^{-})/\Delta\rho_{\infty}(T)$ {\em universally} to ${1/2}(Q_{\scriptsize m}^{+}+Q_{\scriptsize m}^{-})$, where $\Delta\rho_{\infty}$$ = $$\rho^{+}(T)-\rho^{-}(T)$ is the desired width of the coexistence curve. Utilizing the exact limiting $(L$$ \to $$\infty)$ form, the corresponding scaling function can be generated in recursive steps by fitting overlapping data for three or more box sizes, $L_{1}$, $L_{2}$, $...$, $L_{n}$. Starting at a $T_{0}$ sufficiently far below $T_{\scriptsize c}$ and suitably choosing intervals $\Delta T_{j}$$ = $$T_{j+1}-T_{j}$$ > $0 yields $\Delta\rho_{\infty}(T_{j})$ and precisely locates $T_{\scriptsize c}$.

Grand canonical simulations at various levels, $\zeta=5$-20, of fine- lattice discretization are reported for the near-critical 1:1 hard-core electrolyte or RPM. With the aid of finite-size scaling analyses it is shown convincingly that, contrary to recent suggestions, the universal critical behavior is independent of $\zeta$ $(\grtsim 4)$; thus the continuum $(\zeta\to\infty)$ RPM exhibits Ising-type (as against classical, SAW, XY, etc.) criticality. A general consideration of lattice discretization provides effective extrapolation of the {\em intrinsically} erratic $\zeta$-dependence, yielding $(\Tc^ {\ast},\rhoc^{\ast})\simeq (0.0493_{3},0.075)$ for the $\zeta=\infty$ RPM., Comment: 4 pages including 4 figures

A critique is presented of the frequently used Bruce-Wilding (BW) mixed-field scaling method for estimating the critical points of nonsymmetric model fluids from grand canonical simulation data. An explicit, systematic technique for implementing this method is set out thereby revealing clearly a fortunate, close cancelation of contributions from the leading correction- to-scaling and thermal scaling functions that makes the method effective for Ising-type systems but which lacks a general theoretical basis. Pressure mixing is considered in this work which modifies the leading behavior of the critical density estimator while the critical temperature estimator maintains the leading behavior asserted by BW., Comment: Submitted to Journal of Physical Chemistry in honor of Hans Andersen

The vapor-liquid critical behavior of intrinsically asymmetric fluids is studied in finite systems of linear dimensions, $L$, focusing on periodic boundary conditions, as appropriate for simulations. The recently propounded ``complete'' thermodynamic $(L\to\infty)$ scaling theory incorporating pressure mixing in the scaling fields as well as corrections to scaling ${[arXiv:cond-mat/0212145]}$, is extended to finite $L$, initially in a grand canonical representation. The theory allows for a Yang-Yang anomaly in which, when $L\to\infty$, the second temperature derivative, $(d^{2}\mu_{\sigma}/dT^{2})$, of the chemical potential along the phase boundary, $\mu_{\sigma}(T)$, diverges when $T\to\Tc -$. The finite-size behavior of various special {\em critical loci} in the temperature-density or $(T,\rho)$ plane, in particular, the $k$-inflection susceptibility loci and the $Q$-maximal loci -- derived from $Q_{L}(T,<\rho>_{L}) \equiv < m^{2}>^{2}_{L}/< m^{4}>_{L}$ where $m \equiv \rho - <\rho>_{L}$ -- is carefully elucidated and shown to be of value in estimating $\Tc$ and $\rhoc$. Concrete illustrations are presented for the hard-core square-well fluid and for the restricted primitive model electrolyte including an estimate of the correlation exponent $\nu$ that confirms Ising-type character. The treatment is extended to the canonical representation where further complications appear., Comment: 23 pages in the two-column format (including 13 figures) This is Part II of the previous paper [arXiv:cond-mat/0212145]

We present a novel method to derive liquid-gas coexisting densities, $\rho^{\pm}(T)$, from grand canonical simulations (without knowledge of $\Tc$ or criticality class). The minima of $ Q_{L}\equiv< m^{2} >_{L}^{2}/< m^{4}>_{L}$ in an $L$$\times$$L \times$$L$ box with $m = \rho - <\rho>_{L}$ are used to generate recursively an unbiased universal finite-size scaling function. Monte Carlo data for a hard-core square-well fluid and for the restricted primitive model electrolyte yield $\rho^{\pm}$ to $\pm 1$-2% of $\rhoc$ down to 1 part in $10^4$-$10^3$ of $\Tc$ (and confirm well Ising character). Pressure mixing in the scaling fields is unequivocally revealed and indicates Yang-Yang ratios $R_{\mu} = -0.04_{4}$ and $0.2_{6}$ for the two models, respectively., Comment: accepted for publication in Phys. Rev. Lett

The thermodynamic behavior of a fluid near a vapor-liquid and, hence, asymmetric critical point is discussed within a general ``complete'' scaling theory incorporating pressure mixing in the nonlinear scaling fields as well as corrections to scaling. This theory allows for a Yang-Yang anomaly in which \mu_{\sigma}^{\prime\prime}(T), the second temperature derivative of the chemical potential along the phase boundary, diverges like the specific heat when T\to T_{\scriptsize c}; it also generates a leading singular term, |t|^{2\beta}, in the coexistence curve diameter, where t\equiv (T-T_{\scriptsize c}) /T_{\scriptsize c}. The behavior of various special loci, such as the critical isochore, the critical isotherm, the k-inflection loci, on which \chi^{(k)}\equiv \chi(\rho,T)/\rho^{k} (with \chi = \rho^{2} k_{\scriptsize B}TK_{T}) and C_{V}^{(k)}\equiv C_{V}(\rho,T)/\rho^{k} are maximal at fixed T, is carefully elucidated. These results are useful for analyzing simulations and experiments, since particular, nonuniversal values of k specify loci that approach the critical density most rapidly and reflect the pressure-mixing coefficient. Concrete illustrations are presented for the hard-core square-well fluid and for the restricted primitive model electrolyte. For comparison, a discussion of the classical (or Landau) theory is presented briefly and various interesting loci are determined explicitly and illustrated quantitatively for a van der Waals fluid., Comment: 21 pages in two-column format including 8 figures

The character of critical behavior in physical systems depends on the range of interactions. In the limit of infinite range of the interactions, systems will exhibit mean-field critical behavior, i.e., critical behavior not affected by fluctuations of the order parameter. If the interaction range is finite, the critical behavior asymptotically close to the critical point is determined by fluctuations and the actual critical behavior depends on the particular universality class. A variety of systems, including fluids and anisotropic ferromagnets, belongs to the three-dimensional Ising universality class. Recent numerical studies of Ising models with different interaction ranges have revealed a spectacular crossover between the asymptotic fluctuation-induced critical behavior and mean-field-type critical behavior. In this work, we compare these numerical results with a crossover Landau model based on renormalization-group matching. For this purpose we consider an application of the crossover Landau model to the three-dimensional Ising model without fitting to any adjustable parameters. The crossover behavior of the critical susceptibility and of the order parameter is analyzed over a broad range (ten orders) of the scaled distance to the critical temperature. The dependence of the coupling constant on the interaction range, governing the crossover critical behavior, is discussed, Comment: 10 pages in two-column format including 9 figures and 1 table. Submitted to J. Stat. Phys. in honor of M. E. Fisher's 70th birthday

Thermodynamics in the vicinity of a critical endpoint with nonclassical exponents $\alpha$, $\beta$, $\gamma$, $\delta$, $...$ is analyzed in terms of density variables (mole fractions, magnetizations, etc.). The shapes of the isothermal binodals or two-phase coexistence curves are found at andnear the endpoint for symmetric and nonsymmetric situations. The spectator- (or noncritical)-phase binodal at $T=T_{e}$ is characterized by an exponent $(\delta +1)/\delta$ $(\simeq 1.21)$ with leading corrections of relative order $1/\delta$ $(\simeq 0.21)$, $\theta_{4}/\beta\delta$ $(\simeq 0.34)$ and $1 -(\beta\delta)^{-1}$$(\simeq 0.36)$; in contrast to classical (van der Waals, mean field, $...$) theory, the critical endpoint binodal is singular with leading exponent $(1-\alpha)/\beta$ $(\simeq 2.73)$ and corrections which are elucidated; the remaining, $\lambda$-line binodals also display the `renormalized exponent,' $(1-\alpha)/\beta$ butwith more singular corrections. (The numerical values quoted here pertain to $(d=3)$-dimensional-fluid or Ising-type systems.), Comment: Revtex, 5 figures