Showing 1–50 of 53 results for author: Luijten, E

Phase Separation and Ripening in a Viscoelastic Gel

Authors: Tine Curk, Erik Luijten

Abstract: The process of phase separation in elastic solids and viscous fluids is of fundamental importance to the stability and function of soft materials. We explore the dynamics of phase separation and domain growth in a viscoelastic material such as a polymer gel. Using analytical theory and Monte Carlo simulations we report a new domain growth regime, in which the domain size increases algebraically wi… ▽ More The process of phase separation in elastic solids and viscous fluids is of fundamental importance to the stability and function of soft materials. We explore the dynamics of phase separation and domain growth in a viscoelastic material such as a polymer gel. Using analytical theory and Monte Carlo simulations we report a new domain growth regime, in which the domain size increases algebraically with a ripening exponent $α$ that depends on the viscoelastic properties of the material. For a prototypical Maxwell material, we obtain $α=1$, which is markedly different from the well-known Ostwald ripening process with $α=1/3$. We generalize our theory to systems with arbitrary power-law relaxation behavior and discuss our findings in the context of the long-term stability of materials as well as recent experimental results on phase separation in cross-linked networks and cytoskeleton. △ Less

Submitted 23 March, 2023; v1 submitted 21 January, 2022; originally announced January 2022.

Comments: 12 pages, 5 figures

Accelerated simulation method for charge regulation effects

Authors: Tine Curk, Jiaxing Yuan, Erik Luijten

Abstract: The net charge of solvated entities, ranging from polyelectrolytes and biomolecules to charged nanoparticles and membranes, depends on the local dissociation equilibrium of individual ionizable groups. Incorporation of this phenomenon, \emph{charge regulation}, in theoretical and computational models requires dynamic, configuration-dependent recalculation of surface charges and is therefore typica… ▽ More The net charge of solvated entities, ranging from polyelectrolytes and biomolecules to charged nanoparticles and membranes, depends on the local dissociation equilibrium of individual ionizable groups. Incorporation of this phenomenon, \emph{charge regulation}, in theoretical and computational models requires dynamic, configuration-dependent recalculation of surface charges and is therefore typically approximated by assuming constant net charge on particles. Various computational methods exist that address this. We present an alternative, particularly efficient charge regulation Monte Carlo method (CR-MC), which explicitly models the redistribution of individual charges and accurately samples the correct grand-canonical charge distribution. In addition, we provide an open-source implementation in the LAMMPS molecular dynamics (MD) simulation package, resulting in a hybrid MD/CR-MC simulation method. This implementation is designed to handle a wide range of implicit-solvent systems that model discreet ionizable groups or surface sites. The computational cost of the method scales linearly with the number of ionizable groups, thereby allowing accurate simulations of systems containing thousands of individual ionizable sites. By matter of illustration, we use the CR-MC method to quantify the effects of charge regulation on the nature of the polyelectrolyte coil--globule transition and on the effective interaction between oppositely charged nanoparticles. △ Less

Submitted 21 January, 2022; v1 submitted 15 September, 2021; originally announced September 2021.

Journal ref: J. Chem. Phys. (2022)

Charge-regulation effects in nanoparticle self-assembly

Authors: Tine Curk, Erik Luijten

Abstract: Nanoparticles in solution acquire charge through dissociation or association of surface groups. Thus, a proper description of their electrostatic interactions requires the use of charge-regulating boundary conditions rather than the commonly employed constant-charge approximation. We implement a hybrid Monte Carlo/Molecular Dynamics scheme that dynamically adjusts the charges of individual surface… ▽ More Nanoparticles in solution acquire charge through dissociation or association of surface groups. Thus, a proper description of their electrostatic interactions requires the use of charge-regulating boundary conditions rather than the commonly employed constant-charge approximation. We implement a hybrid Monte Carlo/Molecular Dynamics scheme that dynamically adjusts the charges of individual surface groups of objects while evolving their trajectories. Charge-regulation effects are shown to qualitatively change self-assembled structures due to global charge redistribution, stabilizing asymmetric constructs. We delineate under which conditions the conventional constant-charge approximation may be employed and clarify the interplay between charge regulation and dielectric polarization. △ Less

Submitted 9 February, 2021; v1 submitted 17 September, 2020; originally announced September 2020.

Comments: 5 pages

Journal ref: Phys. Rev. Lett. 126, 138003 (2021)

Asymmetric electrolytes near structured dielectric interfaces

Authors: Huanxin Wu, Honghao Li, Francisco J. Solis, Monica Olvera de la Cruz, Erik Luijten

Abstract: The ion distribution of electrolytes near interfaces with dielectric contrast has important consequences for electrochemical processes and many other applications. To date, most studies of such systems have focused on geometrically simple interfaces, for which dielectric effects are analytically solvable or computationally tractable. However, all real surfaces display nontrivial structure at the n… ▽ More The ion distribution of electrolytes near interfaces with dielectric contrast has important consequences for electrochemical processes and many other applications. To date, most studies of such systems have focused on geometrically simple interfaces, for which dielectric effects are analytically solvable or computationally tractable. However, all real surfaces display nontrivial structure at the nanoscale and have, in particular, nonuniform local curvature. Using a recently developed, highly efficient computational method, we investigate the effect of surface geometry on ion distribution and interface polarization. We consider an asymmetric 2:1 electrolyte bounded by a sinusoidally deformed solid surface. We demonstrate that even when the surface is neutral, the electrolyte acquires a nonuniform ion density profile near the surface. This profile is asymmetric and leads to an effective charging of the surface. We furthermore show that the induced charge is modulated by the local curvature. The effective charge is opposite in sign to the multivalent ions and is larger in concave regions of the surface. △ Less

Submitted 16 September, 2018; originally announced September 2018.

Comments: J. Chem. Phys. (accepted)

Dielectric modulation of ion transport near interfaces

Authors: Hanne S. Antila, Erik Luijten

Abstract: Ion mobility and ionic conductance in nanodevices are known to deviate from bulk behavior, a phenomenon often attributed to surface effects. We demonstrate that dielectric mismatch between the electrolyte and the surface can qualitatively alter ionic transport in a counterintuitive manner. Instead of following the polarization-induced modulation of the concentration profile, mobility is enhanced o… ▽ More Ion mobility and ionic conductance in nanodevices are known to deviate from bulk behavior, a phenomenon often attributed to surface effects. We demonstrate that dielectric mismatch between the electrolyte and the surface can qualitatively alter ionic transport in a counterintuitive manner. Instead of following the polarization-induced modulation of the concentration profile, mobility is enhanced or reduced by changes in the ionic atmosphere near the interface and affected by a polarization force parallel to the surface. In addition to revealing this mechanism, we explore the effect of salt concentration and electrostatic coupling. △ Less

Submitted 13 September, 2018; originally announced September 2018.

Journal ref: Phys. Rev. Lett. 120, 135501 (2018)

Accurate and Efficient Numerical Simulation of Dielectrically Anisotropic Particles

Authors: Huanxin Wu, Erik Luijten

Abstract: A variety of electrostatic phenomena, including the structure of electric double layers and the aggregation of charged colloids and proteins, are affected by nonuniform electric permittivity. These effects are frequently ignored in analytical and computational studies, and particularly difficult to handle in situations where multiple dielectric contrasts are present, such as in colloids that are h… ▽ More A variety of electrostatic phenomena, including the structure of electric double layers and the aggregation of charged colloids and proteins, are affected by nonuniform electric permittivity. These effects are frequently ignored in analytical and computational studies, and particularly difficult to handle in situations where multiple dielectric contrasts are present, such as in colloids that are heterogeneous in permittivity. We present an extension to the Iterative Dielectric Solver developed by Barros and Luijten [Phys. Rev. Lett. 113, 017801 (2014)] that makes it possible to accurately compute the polarization of anisotropic particles with multiple dielectric contrasts. This efficient boundary-element method-based approach is applicable to geometries that are not amenable to other solvers, opening the possibility of studying collective phenomena of dielectrically anisotropic particles. We provide insight into the underlying physical reasons for this efficiency. △ Less

Submitted 6 September, 2018; originally announced September 2018.

Comments: J. Chem. Phys. (accepted)

The Effective Temperature Concept Tested in an Active Colloid Mixture

Authors: Ming Han, Jing Yan, Steve Granick, Erik Luijten

Abstract: Thermal energy agitates all matter and its competition with ordering tendencies is one of the most fundamental organizing principles in the physical world. Thus, it is natural to enquire if an effective temperature could result when external energy input enhances agitation. Potentially this could extend the insights of statistical thermodynamics to nonequilibrium systems, but despite proposals tha… ▽ More Thermal energy agitates all matter and its competition with ordering tendencies is one of the most fundamental organizing principles in the physical world. Thus, it is natural to enquire if an effective temperature could result when external energy input enhances agitation. Potentially this could extend the insights of statistical thermodynamics to nonequilibrium systems, but despite proposals that the effective temperature concept may apply to synthetic active matter, biological motors, granular materials and turbulent fluids, its predictive value remains unclear. Here, combining computer simulations and imaging experiments, we design a two-component system of driven Janus colloids such that collisions produced by external energy sources play the role of temperature, and in this system we demonstrate quantitative agreement with hallmarks of statistical thermodynamics for binary phase behavior: the archetypal phase diagram with equilibrium critical exponents, Gaussian displacement distributions, fluctuation-dissipation relations, and capillarity. These quantitative analogies to equilibrium expectations, observed in this decidedly nonequilibrium system, constitute an existence proof from which to compare future theories of nonequilibrium, but limitations of this concept are also highlighted. △ Less

Submitted 16 November, 2016; v1 submitted 13 November, 2016; originally announced November 2016.

Comments: Article unchanged; only corrected misspelling in author name

Journal ref: Proc. Natl. Acad. Sci. 114, 7513-7518 (2017)

Electric double layer of anisotropic dielectric colloids under electric fields

Authors: Ming Han, Huanxin Wu, Erik Luijten

Abstract: Anisotropic colloidal particles constitute an important class of building blocks for self-assembly directed by electrical fields. The aggregation of these building blocks is driven by induced dipole moments, which arise from an interplay between dielectric effects and the electric double layer. For particles that are anisotropic in shape, charge distribution, and dielectric properties, calculation… ▽ More Anisotropic colloidal particles constitute an important class of building blocks for self-assembly directed by electrical fields. The aggregation of these building blocks is driven by induced dipole moments, which arise from an interplay between dielectric effects and the electric double layer. For particles that are anisotropic in shape, charge distribution, and dielectric properties, calculation of the electric double layer requires coupling of the ionic dynamics to a Poisson solver. We apply recently proposed methods to solve this problem for experimentally employed colloids in static and time-dependent electric fields. This allows us to predict the effects of field strength and frequency on the colloidal properties. △ Less

Submitted 10 February, 2016; originally announced February 2016.

Comments: Proceedings of Geilo School 2015 - Cooperative particles: Patchy colloids, active matter and nanofluids

Journal ref: Eur. Phys. J. Special Topics 225, 685-698 (2016)

Orientationally Glassy Crystals of Janus Spheres

Authors: Shan Jiang, Jing Yan, Jonathan K. Whitmer, Stephen M. Anthony, Erik Luijten, Steve Granick

Abstract: Colloidal Janus spheres in water (one hemisphere attractive and the other repulsive) assemble into two-dimensional hexagonal crystals with orientational order controlled by anisotropic interactions. We exploit the decoupled translational and rotational order to quantify the orientational dynamics. Via imaging experiments and Monte Carlo simulations we demonstrate that the correlations in the orien… ▽ More Colloidal Janus spheres in water (one hemisphere attractive and the other repulsive) assemble into two-dimensional hexagonal crystals with orientational order controlled by anisotropic interactions. We exploit the decoupled translational and rotational order to quantify the orientational dynamics. Via imaging experiments and Monte Carlo simulations we demonstrate that the correlations in the orientation of individual Janus spheres exhibit glasslike dynamics that can be controlled via the ionic strength. Thus, these colloidal building blocks provide a particularly suitable model glass system for elucidating nontrivial dynamics arising from directional interactions, not captured by the consideration of just translational order. △ Less

Submitted 10 June, 2014; originally announced June 2014.

Journal ref: Phys. Rev. Lett. 112, 218301 (2014)

Dielectric effects in the self-assembly of binary colloidal aggregates

Authors: Kipton Barros, Erik Luijten

Abstract: Electrostatic interactions play an important role in numerous self-assembly phenomena, including colloidal aggregation. Although colloids typically have a dielectric constant that differs from the surrounding solvent, the effective interactions that arise from inhomogeneous polarization charge distributions are generally neglected in theoretical and computational studies. We introduce an efficient… ▽ More Electrostatic interactions play an important role in numerous self-assembly phenomena, including colloidal aggregation. Although colloids typically have a dielectric constant that differs from the surrounding solvent, the effective interactions that arise from inhomogeneous polarization charge distributions are generally neglected in theoretical and computational studies. We introduce an efficient technique to resolve polarization charges in dynamical dielectric geometries, and demonstrate that dielectric effects \emph{qualitatively} alter the predicted self-assembled structures, with surprising colloidal strings arising from many-body effects. △ Less

Submitted 6 June, 2014; originally announced June 2014.

Comments: Accepted for publication in Physical Review Letters

Journal ref: Phys. Rev. Lett. 113, 017801 (2014)

Efficient and accurate simulation of dynamic dielectric objects

Authors: Kipton Barros, Daniel Sinkovits, Erik Luijten

Abstract: Electrostatic interactions between dielectric objects are complex and of a many-body nature, owing to induced surface bound charge. We present a collection of techniques to simulate dynamical dielectric objects. We calculate the surface bound charge from a matrix equation using the Generalized Minimal Residue method (GMRES). Empirically, we find that GMRES converges very quickly. Indeed, our detai… ▽ More Electrostatic interactions between dielectric objects are complex and of a many-body nature, owing to induced surface bound charge. We present a collection of techniques to simulate dynamical dielectric objects. We calculate the surface bound charge from a matrix equation using the Generalized Minimal Residue method (GMRES). Empirically, we find that GMRES converges very quickly. Indeed, our detailed analysis suggests that the relevant matrix has a very compact spectrum for all non-degenerate dielectric geometries. Each GMRES iteration can be evaluated using a fast Ewald solver with cost that scales linearly or near-linearly in the number of surface charge elements. We analyze several previously proposed methods for calculating the bound charge, and show that our approach compares favorably. △ Less

Submitted 7 January, 2014; originally announced January 2014.

Comments: 16 pages, 2 figures

Journal ref: J. Chem. Phys. 140, 064903 (2014)

Monte Carlo cluster algorithm for fluid phase transitions in highly size-asymmetrical binary mixtures

Authors: Douglas J. Ashton, Jiwen Liu, Erik Luijten, Nigel B. Wilding

Abstract: Highly size-asymmetrical fluid mixtures arise in a variety of physical contexts, notably in suspensions of colloidal particles to which much smaller particles have been added in the form of polymers or nanoparticles. Conventional schemes for simulating models of such systems are hamstrung by the difficulty of relaxing the large species in the presence of the small one. Here we describe how the rej… ▽ More Highly size-asymmetrical fluid mixtures arise in a variety of physical contexts, notably in suspensions of colloidal particles to which much smaller particles have been added in the form of polymers or nanoparticles. Conventional schemes for simulating models of such systems are hamstrung by the difficulty of relaxing the large species in the presence of the small one. Here we describe how the rejection-free geometrical cluster algorithm (GCA) of Liu and Luijten [Phys. Rev. Lett 92, 035504 (2004)] can be embedded within a restricted Gibbs ensemble to facilitate efficient and accurate studies of fluid phase behavior of highly size-asymmetrical mixtures. After providing a detailed description of the algorithm, we summarize the bespoke analysis techniques of Ashton et al. [J. Chem. Phys. 132, 074111 (2010)] that permit accurate estimates of coexisting densities and critical-point parameters. We apply our methods to study the liquid--vapor phase diagram of a particular mixture of Lennard-Jones particles having a 10:1 size ratio. As the reservoir volume fraction of small particles is increased in the range 0--5%, the critical temperature decreases by approximately 50%, while the critical density drops by some 30%. These trends imply that in our system, adding small particles decreases the net attraction between large particles, a situation that contrasts with hard-sphere mixtures where an attractive depletion force occurs. △ Less

Submitted 21 July, 2010; originally announced July 2010.

Comments: 11 pages, 10 figures

Journal ref: J. Chem. Phys. 133, 194102 (2010)

Polyelectrolyte Condensation Induced by Linear Cations

Authors: Camilo Guáqueta, Erik Luijten

Abstract: We examine the role of the condensing agent in the formation of polyelectrolyte bundles, via grand-canonical Monte Carlo simulations. Following recent experiments we use linear, rigid divalent cations of various lengths to induce condensation. Our results clarify and explain the experimental results for short cations. For longer cations we observe novel condensation behavior owing to alignment o… ▽ More We examine the role of the condensing agent in the formation of polyelectrolyte bundles, via grand-canonical Monte Carlo simulations. Following recent experiments we use linear, rigid divalent cations of various lengths to induce condensation. Our results clarify and explain the experimental results for short cations. For longer cations we observe novel condensation behavior owing to alignment of the cations. We also study the role of the polyelectrolyte surface charge density, and find a nonmonotonic variation in bundle stability. This nonmonotonicity captures two trends that have been observed in separate experiments. △ Less

Submitted 21 September, 2007; originally announced September 2007.

Comments: To appear in Physical Review Letters

Journal ref: Phys. Rev. Lett. 99, 138302 (2007)

Salt-induced collapse and reexpansion of highly charged flexible polyelectrolytes

Authors: Pai-Yi Hsiao, Erik Luijten

Abstract: We study the salt-dependent conformations of dilute flexible polyelectrolytes in solution via computer simulations. Low concentrations of multivalent salt induce the known conformational collapse of individual polyelectrolyte chains, but as the salt concentration is increased further this is followed by a reexpansion. We explicitly demonstrate that multivalent counterions can overcompensate the… ▽ More We study the salt-dependent conformations of dilute flexible polyelectrolytes in solution via computer simulations. Low concentrations of multivalent salt induce the known conformational collapse of individual polyelectrolyte chains, but as the salt concentration is increased further this is followed by a reexpansion. We explicitly demonstrate that multivalent counterions can overcompensate the bare charge of the chain in the reexpansion regime. Both the degree of reexpansion and the occurrence of overcharging sensitively depend on ion size. Our findings are relevant for a wide range of salt-induced complexation phenomena. △ Less

Submitted 5 October, 2006; originally announced October 2006.

Journal ref: Phys. Rev. Lett. 97, 148301 (2006)

Simulation of phase transitions in highly asymmetric fluid mixtures

Authors: Jiwen Liu, Nigel B. Wilding, Erik Luijten

Abstract: We present a novel method for the accurate numerical determination of the phase behavior of fluid mixtures having large particle size asymmetries. By incorporating the recently developed geometric cluster algorithm within a restricted Gibbs ensemble, we are able to probe directly the density and concentration fluctuations that drive phase transitions, but that are inaccessible to conventional si… ▽ More We present a novel method for the accurate numerical determination of the phase behavior of fluid mixtures having large particle size asymmetries. By incorporating the recently developed geometric cluster algorithm within a restricted Gibbs ensemble, we are able to probe directly the density and concentration fluctuations that drive phase transitions, but that are inaccessible to conventional simulation algorithms. We develop a finite-size scaling theory that relates these density fluctuations to those of the grand-canonical ensemble, thereby enabling accurate location of critical points and coexistence curves of multicomponent fluids. Several illustrative examples are presented. △ Less

Submitted 10 August, 2006; originally announced August 2006.

Comments: Accepted for publication in Physical Review Letters

Journal ref: Phys. Rev. Lett. 97, 115705 (2006)

Quantum criticality, lines of fixed points, and phase separation in doped two-dimensional quantum dimer models

Authors: Stefanos Papanikolaou, Erik Luijten, Eduardo Fradkin

Abstract: We study phase diagrams of a class of doped quantum dimer models on the square lattice with ground-state wave functions whose amplitudes have the form of the Gibbs weights of a classical doped dimer model. In this dimer model, parallel neighboring dimers have attractive interactions, whereas neighboring holes either do not interact or have a repulsive interaction. We investigate the behavior of… ▽ More We study phase diagrams of a class of doped quantum dimer models on the square lattice with ground-state wave functions whose amplitudes have the form of the Gibbs weights of a classical doped dimer model. In this dimer model, parallel neighboring dimers have attractive interactions, whereas neighboring holes either do not interact or have a repulsive interaction. We investigate the behavior of this system via analytic methods and by Monte Carlo simulations. At zero doping, we confirm the existence of a Kosterlitz-Thouless transition from a quantum critical phase to a columnar phase. At low hole densities we find a dimer-hole liquid phase and a columnar phase, separated by a phase boundary which is a line of critical points with varying exponents. We demonstrate that this line ends at a multicritical point where the transition becomes first order and the system phase separates. The first-order transition coexistence curve is shown to become unstable with respect to more complex inhomogeneous phases in the presence of direct hole-hole interactions. We also use a variational approach to determine the spectrum of low-lying density fluctuations in the dimer-hole fluid phase. △ Less

Submitted 24 October, 2007; v1 submitted 12 July, 2006; originally announced July 2006.

Comments: New and vastly expanded version of the original letter; 31 pages, 8 sections, 2 appendices, 19 figures, 86 references, new abstract. Final published version on Phys. Rev. B

Journal ref: Phys. Rev. B 76, 134514 (2007)

Confinement-driven translocation of a flexible polymer

Authors: Angelo Cacciuto, Erik Luijten

Abstract: We consider the escape of a flexible, self-avoiding polymer chain out of a confined geometry. By means of simulations, we demonstrate that the translocation time can be described by a simple scaling law that exhibits a nonlinear dependence on the degree of polymerization and that is sensitive to the nature of the confining geometry. These results contradict earlier predictions but are in agreeme… ▽ More We consider the escape of a flexible, self-avoiding polymer chain out of a confined geometry. By means of simulations, we demonstrate that the translocation time can be described by a simple scaling law that exhibits a nonlinear dependence on the degree of polymerization and that is sensitive to the nature of the confining geometry. These results contradict earlier predictions but are in agreement with recently confirmed geometry-dependent expressions for the free energy of confinement. △ Less

Submitted 29 May, 2006; originally announced May 2006.

Comments: Accepted for publication in Phys. Rev. Lett

Journal ref: Phys. Rev. Lett. 96, 238104 (2006)

Colloidal stabilization via nanoparticle haloing

Authors: Jiwen Liu, Erik Luijten

Abstract: We present a detailed numerical study of effective interactions between micron-sized silica spheres, induced by highly charged zirconia nanoparticles. It is demonstrated that the effective interactions are consistent with a recently discovered mechanism for colloidal stabilization. In accordance with the experimental observations, small nanoparticle concentrations induce an effective repulsion t… ▽ More We present a detailed numerical study of effective interactions between micron-sized silica spheres, induced by highly charged zirconia nanoparticles. It is demonstrated that the effective interactions are consistent with a recently discovered mechanism for colloidal stabilization. In accordance with the experimental observations, small nanoparticle concentrations induce an effective repulsion that counteracts the intrinsic van der Waals attraction between the colloids and thus stabilizes the suspension. At higher nanoparticle concentrations an attractive potential is recovered, resulting in reentrant gelation. Monte Carlo simulations of this highly size-asymmetric mixture are made possible by means of a geometric cluster Monte Carlo algorithm. A comparison is made to results obtained from the Ornstein-Zernike equations with the hypernetted-chain closure. △ Less

Submitted 12 August, 2005; originally announced August 2005.

Journal ref: Phys. Rev. E 72, 061401 (2005)

Critical polymer-polymer phase separation in ternary solutions

Authors: Lei Guo, Erik Luijten

Abstract: We study polymer-polymer phase separation in a common good solvent by means of Monte Carlo simulations of the bond-fluctuation model. Below a critical, chain-length dependent concentration, no phase separation occurs. For higher concentrations, the critical demixing temperature scales nonlinearly with the total monomer concentration, with a power law relatively close to a renormalization-group p… ▽ More We study polymer-polymer phase separation in a common good solvent by means of Monte Carlo simulations of the bond-fluctuation model. Below a critical, chain-length dependent concentration, no phase separation occurs. For higher concentrations, the critical demixing temperature scales nonlinearly with the total monomer concentration, with a power law relatively close to a renormalization-group prediction based on "blob" scaling arguments. We point out that earlier simulations and experiments have tested this power-law dependence at concentrations outside the validity regime of the scaling arguments. The critical amplitudes of the order parameter and the zero-angle scattering intensity also exhibit chain-length dependences that differ from the conventional predictions but are in excellent agreement with the renormalization-group results. In addition, we characterize the variation of the average coil shape upon phase separation. △ Less

Submitted 26 July, 2005; originally announced July 2005.

Comments: To appear in J. Chem. Phys

Journal ref: J. Chem. Phys. 123, 074907 (2005)

Structure and stability of self-assembled actin-lysozyme complexes in salty water

Authors: Lori K. Sanders, Camilo Guaqueta, Thomas E. Angelini, Jae-Wook Lee, Scott C. Slimmer, Erik Luijten, Gerard C. L. Wong

Abstract: Interactions between actin, an anionic polyelectrolyte, and lysozyme, a cationic globular protein, have been examined using a combination of synchrotron small-angle x-ray scattering and molecular dynamics simulations. Lysozyme initially bridges pairs of actin filaments, which relax into hexagonally-coordinated columnar complexes comprised of actin held together by incommensurate one-dimensional… ▽ More Interactions between actin, an anionic polyelectrolyte, and lysozyme, a cationic globular protein, have been examined using a combination of synchrotron small-angle x-ray scattering and molecular dynamics simulations. Lysozyme initially bridges pairs of actin filaments, which relax into hexagonally-coordinated columnar complexes comprised of actin held together by incommensurate one-dimensional close-packed arrays of lysozyme macroions. These complexes are found to be stable even in the presence of significant concentrations of monovalent salt, which is quantitatively explained from a redistribution of salt between the condensed and the aqueous phases. △ Less

Submitted 25 July, 2005; originally announced July 2005.

Comments: Accepted for publication in Phys. Rev. Lett

Journal ref: Phys. Rev. Lett. 95, 108302 (2005)

Screening in Ionic Systems: Simulations for the Lebowitz Length

Authors: Young C. Kim, Erik Luijten, Michael E. Fisher

Abstract: Simulations of the Lebowitz length, $ξ_{\text{L}}(T,ρ)$, are reported for t he restricted primitive model hard-core (diameter $a$) 1:1 electrolyte for densi ties $ρ\lesssim 4ρ_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 $ξ_{\text{L}}$. On extrapolation to the bulk limit for $T\gtrsim 10T_c$ t… ▽ More Simulations of the Lebowitz length, $ξ_{\text{L}}(T,ρ)$, are reported for t he restricted primitive model hard-core (diameter $a$) 1:1 electrolyte for densi ties $ρ\lesssim 4ρ_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 $ξ_{\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 $ρ> {1/10}ρ_c$ (with $ρ_c a^3 \simeq 0.08$). At highe r densities $ξ_{\text{L}}$ rises above the Debye length, $ξ_{\text{D}} \prop to \sqrt{T/ρ}$, by 10-30% (upto $ρ\simeq 1.3ρ_c$); the variation is portrayed fairly well by generalized Debye-Hückel theory (Lee and Fisher, 19 96). On approaching criticality at fixed $ρ$ or fixed $T$, $ξ_{\text{L}}(T, ρ)$ remains finite with $ξ_{\text{L}}^c \simeq 0.30 a \simeq 1.3 ξ_{\text {D}}^c$ but displays a weak entropy-like singularity. △ Less

Submitted 14 July, 2005; originally announced July 2005.

Comments: 4 pages 5 figures

Journal ref: Phys. Rev. Lett. 95, 145701 (2005)

Geometric Cluster Algorithm for Interacting Fluids

Authors: Erik Luijten, Jiwen Liu

Abstract: We discuss a new Monte Carlo algorithm for the simulation of complex fluids. This algorithm employs geometric operations to identify clusters of particles that can be moved in a rejection-free way. It is demonstrated that this geometric cluster algorithm (GCA) constitutes the continuum generalization of the Swendsen-Wang and Wolff cluster algorithms for spin systems. Because of its nonlocal natu… ▽ More We discuss a new Monte Carlo algorithm for the simulation of complex fluids. This algorithm employs geometric operations to identify clusters of particles that can be moved in a rejection-free way. It is demonstrated that this geometric cluster algorithm (GCA) constitutes the continuum generalization of the Swendsen-Wang and Wolff cluster algorithms for spin systems. Because of its nonlocal nature, it is particularly well suited for the simulation of fluid systems containing particles of widely varying sizes. The efficiency improvement with respect to conventional simulation algorithms is a rapidly growing function of the size asymmetry between the constituents of the system. We study the cluster-size distribution for a Lennard-Jones fluid as a function of density and temperature and provide a comparison between the generalized GCA and the hard-core GCA for a size-asymmetric mixture with Yukawa-type couplings. △ Less

Submitted 17 March, 2005; originally announced March 2005.

Comments: To appear in "Computer Simulation Studies in Condensed-Matter Physics XVII". Edited by D.P. Landau, S.P. Lewis and H.B. Schuettler. Springer, Heidelberg, 2005

Generalized Geometric Cluster Algorithm for Fluid Simulation

Authors: Jiwen Liu, Erik Luijten

Abstract: We present a detailed description of the generalized geometric cluster algorithm for the efficient simulation of continuum fluids. The connection with well-known cluster algorithms for lattice spin models is discussed, and an explicit full cluster decomposition is derived for a particle configuration in a fluid. We investigate a number of basic properties of the geometric cluster algorithm, incl… ▽ More We present a detailed description of the generalized geometric cluster algorithm for the efficient simulation of continuum fluids. The connection with well-known cluster algorithms for lattice spin models is discussed, and an explicit full cluster decomposition is derived for a particle configuration in a fluid. We investigate a number of basic properties of the geometric cluster algorithm, including the dependence of the cluster-size distribution on density and temperature. Practical aspects of its implementation and possible extensions are discussed. The capabilities and efficiency of our approach are illustrated by means of two example studies. △ Less

Submitted 17 March, 2005; originally announced March 2005.

Comments: Accepted for publication in Phys. Rev. E. Follow-up to cond-mat/0412744

Journal ref: Phys. Rev. E 71, 066701 (2005)

Rejection-free Geometric Cluster Algorithm for Complex Fluids

Authors: Jiwen Liu, Erik Luijten

Abstract: We present a novel, generally applicable Monte Carlo algorithm for the simulation of fluid systems. Geometric transformations are used to identify clusters of particles in such a manner that every cluster move is accepted, irrespective of the nature of the pair interactions. The rejection-free and non-local nature of the algorithm make it particularly suitable for the efficient simulation of com… ▽ More We present a novel, generally applicable Monte Carlo algorithm for the simulation of fluid systems. Geometric transformations are used to identify clusters of particles in such a manner that every cluster move is accepted, irrespective of the nature of the pair interactions. The rejection-free and non-local nature of the algorithm make it particularly suitable for the efficient simulation of complex fluids with components of widely varying size, such as colloidal mixtures. Compared to conventional simulation algorithms, typical efficiency improvements amount to several orders of magnitude. △ Less

Submitted 29 December, 2004; originally announced December 2004.

Journal ref: Phys. Rev. Lett. 92, 035504 (2004)

Stabilization of colloidal suspensions by means of highly-charged nanoparticles

Authors: Jiwen Liu, Erik Luijten

Abstract: We employ a novel Monte Carlo simulation scheme to elucidate the stabilization of neutral colloidal microspheres by means of highly-charged nanoparticles [V. Tohver et al., Proc. Natl. Acad. Sci. U.S.A. 98, 8950 (2001)]. In accordance with the experimental observations, we find that small nanoparticle concentrations induce an effective repulsion that prevents gelation caused by the intrinsic van… ▽ More We employ a novel Monte Carlo simulation scheme to elucidate the stabilization of neutral colloidal microspheres by means of highly-charged nanoparticles [V. Tohver et al., Proc. Natl. Acad. Sci. U.S.A. 98, 8950 (2001)]. In accordance with the experimental observations, we find that small nanoparticle concentrations induce an effective repulsion that prevents gelation caused by the intrinsic van der Waals attraction between colloids. Higher nanoparticle concentrations induce an attractive potential which is, however, qualitatively different from the regular depletion attraction. We also show how colloid-nanoparticle size asymmetry and nanoparticle charge can be used to manipulate the effective interactions. △ Less

Submitted 22 December, 2004; v1 submitted 10 November, 2004; originally announced November 2004.

Comments: Accepted for publication in Physical Review Letters. See also S. Karanikas and A.A. Louis, cond-mat/0411279. Updated to synchronize with published version

Journal ref: Phys. Rev. Lett. 93, 247802 (2004)

Reversible Gel Formation of Triblock Copolymers Studied by Molecular Dynamics Simulation

Authors: Lei Guo, Erik Luijten

Abstract: Molecular dynamics simulations have been employed to study the formation of a physical (thermoreversible) gel by amphiphilic A-B-A triblock copolymers in aqueous solution. In order to mimic the structure of hydrogel-forming polypeptides employed in experiments [W.A. Petka et al., Science 281, 389 (1998)], the endblocks of the polymer chains are modeled as hydrophobic rods representing the alpha-… ▽ More Molecular dynamics simulations have been employed to study the formation of a physical (thermoreversible) gel by amphiphilic A-B-A triblock copolymers in aqueous solution. In order to mimic the structure of hydrogel-forming polypeptides employed in experiments [W.A. Petka et al., Science 281, 389 (1998)], the endblocks of the polymer chains are modeled as hydrophobic rods representing the alpha-helical part of the polypeptides whereas the central B-block is hydrophilic and semi-flexible. We have determined structural properties, such as the hydrophobic cluster-size distribution function, the geometric percolation point and pair correlation functions, and related these to the dynamical properties of the system. Upon decrease of the temperature, a network structure is formed in which bundles of endblocks act as network junctions. Both at short and medium distances an increased ordering is observed, as characterized by the pair correlation function. Micelle formation and the corresponding onset of geometric percolation induce a strong change in dynamical quantities, e.g., in the diffusion constant and the viscosity, and causes the system to deviate from the Stokes-Einstein relation. The dynamical properties show a temperature dependence that is strongly reminiscent of the behavior of glass-forming liquids. The appearance of a plateau in the stress autocorrelation function suggests that the system starts to exhibit a solid-like response to applied stress once the network structure has been formed, although the actual sol-gel transition occurs only at a considerably lower temperature. △ Less

Submitted 13 October, 2004; originally announced October 2004.

Comments: Submitted to J. Polymer Science Part B: Polymer Physics

Journal ref: J. Polym. Sci. Part B: Polym. Phys. 43, 959 (2005)

Precise Simulation of Near-critical Fluid Coexistence

Authors: Young C. Kim, Michael E. Fisher, Erik Luijten

Abstract: We present a novel method to derive liquid-gas coexisting densities, $ρ^{\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 = ρ- <ρ>_{L}$ are used to generate recursively an unbiased universal finite-size scaling function. Monte Carlo data for a hard-… ▽ More We present a novel method to derive liquid-gas coexisting densities, $ρ^{\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 = ρ- <ρ>_{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 $ρ^{\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_μ = -0.04_{4}$ and $0.2_{6}$ for the two models, respectively. △ Less

Submitted 23 June, 2003; v1 submitted 1 April, 2003; originally announced April 2003.

Comments: accepted for publication in Phys. Rev. Lett

Journal ref: Phys. Rev. Lett. vol. 91, 065701 (2003)

Dynamical Properties of the Slithering Snake Algorithm: A numerical test of the activated reptation hypothesis

Authors: L. Mattioni, J. P. Wittmer, J. Baschnagel, J. -L. Barrat, E. Luijten

Abstract: The correlations in the motion of reptating polymers in their melt are investigated by means of kinetic Monte Carlo simulations of the three dimensional slithering snake version of the bond-fluctuation model. Surprisingly, the slithering snake dynamics becomes inconsistent with classical reptation predictions at high chain overlap (either chain length $N$ or volume fraction $φ$) where the rela… ▽ More The correlations in the motion of reptating polymers in their melt are investigated by means of kinetic Monte Carlo simulations of the three dimensional slithering snake version of the bond-fluctuation model. Surprisingly, the slithering snake dynamics becomes inconsistent with classical reptation predictions at high chain overlap (either chain length $N$ or volume fraction $φ$) where the relaxation times increase much faster than expected. This is due to the anomalous curvilinear diffusion in a finite time window whose upper bound $τ_+$ is set by the chain end density $φ/N$. Density fluctuations created by passing chain ends allow a reference polymer to break out of the local cage of immobile obstacles created by neighboring chains. The dynamics of dense solutions of snakes at $t \ll τ_+$ is identical to that of a benchmark system where all but one chain are frozen. We demonstrate that it is the slow creeping of a chain out of its correlation hole which causes the subdiffusive dynamical regime. Our results are in good qualitative agreement with the activated reptation scheme proposed recently by Semenov and Rubinstein [Eur. Phys. J. B, {\bf 1} (1998) 87]. Additionally, we briefly comment on the relevance of local relaxation pathways within a slithering snake scheme. Our preliminary results suggest that a judicious choice of the ratio of local to slithering snake moves is crucial to equilibrate a melt of long chains efficiently. △ Less

Submitted 22 March, 2003; v1 submitted 17 December, 2002; originally announced December 2002.

Comments: 24 pages, 18 figures, submitted to EPJ E

Journal ref: Eur. Phys. J. E 10, 369 (2003).

Crossover critical behavior in the three-dimensional Ising model

Authors: Young C. Kim, Mikhail A. Anisimov, Jan V. Sengers, Erik Luijten

Abstract: 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 fluctua… ▽ More 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 △ Less

Submitted 8 March, 2002; originally announced March 2002.

Comments: 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

Journal ref: J. Stat. Phys. Vol. 110, 591-609 (2003)

The boundary between long-range and short-range critical behavior

Authors: Erik Luijten, Henk W. J. Blöte

Abstract: We investigate phase transitions of two-dimensional Ising models with power-law interactions, using an efficient Monte Carlo algorithm. For slow decay, the transition is of the mean-field type; for fast decay, it belongs to the short-range Ising universality class. We focus on the intermediate range, where the critical exponents depend continuously on the power law. We find that the boundary wit… ▽ More We investigate phase transitions of two-dimensional Ising models with power-law interactions, using an efficient Monte Carlo algorithm. For slow decay, the transition is of the mean-field type; for fast decay, it belongs to the short-range Ising universality class. We focus on the intermediate range, where the critical exponents depend continuously on the power law. We find that the boundary with short-range critical behavior occurs for interactions depending on distance r as r^{-15/4}. This answers a long-standing controversy between mutually conflicting renormalization-group analyses. △ Less

Submitted 27 December, 2001; originally announced December 2001.

Comments: See http://ariadne.mse.uiuc.edu/ for related papers

Journal ref: Phys. Rev. Lett. 89, 025703 (2002).

Universality class of criticality in the restricted primitive model electrolyte

Authors: Erik Luijten, Michael E. Fisher, Athanassios Z. Panagiotopoulos

Abstract: The 1:1 equisized hard-sphere electrolyte or restricted primitive model has been simulated via grand-canonical fine-discretization Monte Carlo. Newly devised unbiased finite-size extrapolation methods using temperature-density, (T, rho), loci of inflections, Q = <m^2>^2/<m^4> maxima, canonical and C_V criticality, yield estimates of (T_c, rho_c) to +- (0.04, 3)%. Extrapolated exponents and Q-rat… ▽ More The 1:1 equisized hard-sphere electrolyte or restricted primitive model has been simulated via grand-canonical fine-discretization Monte Carlo. Newly devised unbiased finite-size extrapolation methods using temperature-density, (T, rho), loci of inflections, Q = <m^2>^2/<m^4> maxima, canonical and C_V criticality, yield estimates of (T_c, rho_c) to +- (0.04, 3)%. Extrapolated exponents and Q-ratio are (gamma, nu, Q_c) = [1.24(3), 0.63(3); 0.624(2)] which support Ising (n = 1) behavior with (1.23_9, 0.630_3; 0.623_6), but exclude classical, XY (n = 2), SAW (n = 0), and n = 1 criticality with potentials phi(r)>Phi/r^{4.9} when r \to \infty. △ Less

Submitted 20 December, 2001; originally announced December 2001.

Journal ref: Phys. Rev. Lett. 88, 185701 (2002)

On quantum effects near the liquid-vapor transition in helium

Authors: Martin H. Müser, Erik Luijten

Abstract: The liquid-vapor transition in He-3 and He-4 is investigated by means of path-integral molecular dynamics and the quantum virial expansion. Both methods are applied to the critical isobar and the critical isochore. While previous path-integral simulations have mainly considered the lambda transition and superfluid regime in He-4, we focus on the vicinity of the critical point and obtain good agr… ▽ More The liquid-vapor transition in He-3 and He-4 is investigated by means of path-integral molecular dynamics and the quantum virial expansion. Both methods are applied to the critical isobar and the critical isochore. While previous path-integral simulations have mainly considered the lambda transition and superfluid regime in He-4, we focus on the vicinity of the critical point and obtain good agreement with experimental results for the molar volume and the internal energy down to subcritical temperatures. We find that an effective classical potential that properly describes the two-particle radial distribution function exhibits a strong temperature dependence near the critical temperature. This contrasts with the behavior of essentially classical systems like xenon, where the effective potential is independent of temperature. It is conjectured that, owing to this difference in behavior between classical and quantum-mechanical systems, the crossover behavior observed for helium in the vicinity of the critical point differs qualitatively from that of other simple liquids. △ Less

Submitted 14 May, 2001; originally announced May 2001.

Journal ref: J. Chem. Phys. 116, 1621 (2002)

Criticality in one dimension with inverse square-law potentials

Authors: Erik Luijten, Holger Messingfeld

Abstract: It is demonstrated that the scaled order parameter for ferromagnetic Ising and three-state Potts chains with inverse-square interactions exhibits a universal critical jump, in analogy with the superfluid density in helium films. Renormalization-group arguments are combined with numerical simulations of systems containing up to one million lattice sites to accurately determine the critical proper… ▽ More It is demonstrated that the scaled order parameter for ferromagnetic Ising and three-state Potts chains with inverse-square interactions exhibits a universal critical jump, in analogy with the superfluid density in helium films. Renormalization-group arguments are combined with numerical simulations of systems containing up to one million lattice sites to accurately determine the critical properties of these models. In strong contrast with earlier work, compelling quantitative evidence for the Kosterlitz--Thouless-like character of the phase transition is provided. △ Less

Submitted 10 April, 2001; originally announced April 2001.

Comments: To appear in Phys. Rev. Lett

Journal ref: Phys. Rev. Lett. 86, 5305 (2001)

The heat capacity of the restricted primitive model electrolyte

Authors: Erik Luijten, Michael E. Fisher, Athanassios Z. Panagiotopoulos

Abstract: The constant-volume heat capacity, C_V(T, rho), of the restricted primitive model (RPM) electrolyte is considered in the vicinity of its critical point. It is demonstrated that, despite claims, recent simulations for finite systems do not convincingly indicate the absence of a divergence in C_V(T, rho)--which would point to non-Ising-type criticality. The strong qualitative difference between C_… ▽ More The constant-volume heat capacity, C_V(T, rho), of the restricted primitive model (RPM) electrolyte is considered in the vicinity of its critical point. It is demonstrated that, despite claims, recent simulations for finite systems do not convincingly indicate the absence of a divergence in C_V(T, rho)--which would point to non-Ising-type criticality. The strong qualitative difference between C_V for the RPM and for a Lennard-Jones fluid is shown to result from the low critical density of the former. If one considers the theoretically preferable configurational heat-capacity density, C_V/V, the finite-size results for the two systems display qualitatively similar behavior on near-critical isotherms. △ Less

Submitted 21 November, 2000; originally announced November 2000.

Comments: 5 Pages, including 5 EPS figures. Also available as PDF file at http://pallas.umd.edu/~luijten/erikpubs.html

Journal ref: J. Chem. Phys. 114, 5468 (2001)

A Kac-potential treatment of nonintegrable interactions

Authors: Benjamin P. Vollmayr-Lee, Erik Luijten

Abstract: We consider d-dimensional systems with nonintegrable, algebraically decaying pairwise interactions. It is shown that, upon introduction of periodic boundary conditions and a long-distance cutoff in the interaction range, the bulk thermodynamics can be obtained rigorously by means of a Kac-potential treatment, leading to an exact, mean-field-like theory. This explains various numerical results re… ▽ More We consider d-dimensional systems with nonintegrable, algebraically decaying pairwise interactions. It is shown that, upon introduction of periodic boundary conditions and a long-distance cutoff in the interaction range, the bulk thermodynamics can be obtained rigorously by means of a Kac-potential treatment, leading to an exact, mean-field-like theory. This explains various numerical results recently obtained for finite systems in the context of ``nonextensive thermodynamics,'' and in passing exposes a strong regulator dependence not discussed in these studies. Our findings imply that, contrary to some claims, Boltzmann-Gibbs statistics are sufficient for a standard description of this class of nonintegrable interactions. △ Less

Submitted 3 September, 2000; originally announced September 2000.

Comments: RevTeX, 8 pages, 1 figure

Journal ref: Phys. Rev. E 63, 031108 (2001)

Crossover behavior in He3 and Xe near their liquid-vapor critical point

Authors: Erik Luijten, Horst Meyer

Abstract: We present a detailed discussion of the crossover from mean-field to Ising critical behavior upon approach of the critical point, both for He3 and Xe. By combining different sets of experimental data, we are able to cover an unusually large temperature range on either side of the critical temperature T_c. Below T_c, we thus can make an accurate comparison with a recent calculation for the crosso… ▽ More We present a detailed discussion of the crossover from mean-field to Ising critical behavior upon approach of the critical point, both for He3 and Xe. By combining different sets of experimental data, we are able to cover an unusually large temperature range on either side of the critical temperature T_c. Below T_c, we thus can make an accurate comparison with a recent calculation for the crossover of the coexistence curve. For the regime above T_c, an analysis of the compressibility demonstrates that the crossover regime in He3 is unexpectedly widened by a subtle interplay between quantum and critical fluctuations. △ Less

Submitted 29 June, 2000; originally announced June 2000.

Comments: Accepted for publication in Phys. Rev. E

Journal ref: Phys. Rev. E 62, 3257 (2000)

Comment on ``Scaling Laws for a System with Long-Range Interactions within Tsallis Statistics''

Authors: Benjamin P. Vollmayr-Lee, Erik Luijten

Abstract: In their recent Letter [Phys. Rev. Lett. 83, 4233 (1999)], Salazar and Toral (ST) study numerically a finite Ising chain with non-integrable interactions decaying like 1/r^(d+sigma) where -d <= sigma <= 0 (like ST, we discuss general dimensionality d). In particular, they explore a presumed connection between non-integrable interactions and Tsallis's non-extensive statistics. We point out that (… ▽ More In their recent Letter [Phys. Rev. Lett. 83, 4233 (1999)], Salazar and Toral (ST) study numerically a finite Ising chain with non-integrable interactions decaying like 1/r^(d+sigma) where -d <= sigma <= 0 (like ST, we discuss general dimensionality d). In particular, they explore a presumed connection between non-integrable interactions and Tsallis's non-extensive statistics. We point out that (i) non-integrable interactions provide no more motivation for Tsallis statistics than do integrable interactions, i.e., Gibbs statistics remain meaningful for the non-integrable case, and in fact provide a {\em complete and exact treatment}; and (ii) there are undesirable features of the method ST use to regulate the non-integrable interactions. △ Less

Submitted 25 May, 2000; originally announced May 2000.

Comments: Accepted for publication in Phys. Rev. Lett

Journal ref: Phys. Rev. Lett. 85, 470 (2000)

Test of renormalization predictions for universal finite-size scaling functions

Authors: Erik Luijten

Abstract: We calculate universal finite-size scaling functions for systems with an n-component order parameter and algebraically decaying interactions. Just as previously has been found for short-range interactions, this leads to a singular epsilon-expansion, where epsilon is the distance to the upper critical dimension. Subsequently, we check the results by numerical simulations of spin models in the sam… ▽ More We calculate universal finite-size scaling functions for systems with an n-component order parameter and algebraically decaying interactions. Just as previously has been found for short-range interactions, this leads to a singular epsilon-expansion, where epsilon is the distance to the upper critical dimension. Subsequently, we check the results by numerical simulations of spin models in the same universality class. Our systems offer the essential advantage that epsilon can be varied continuously, allowing an accurate examination of the region where epsilon is small. The numerical calculations turn out to be in striking disagreement with the predicted singularity. △ Less

Submitted 2 November, 1999; originally announced November 1999.

Comments: 6 pages, including 3 EPS figures. To appear in Phys. Rev. E. Also available as PDF file at http://www.cond-mat.physik.uni-mainz.de/~luijten/erikpubs.html

Journal ref: Phys. Rev. E 60, 7558 (1999)

Optimized energy calculation in lattice systems with long-range interactions

Authors: Michael Krech, Erik Luijten

Abstract: We discuss an efficient approach to the calculation of the internal energy in numerical simulations of spin systems with long-range interactions. Although, since the introduction of the Luijten-Blöte algorithm, Monte Carlo simulations of these systems no longer pose a fundamental problem, the energy calculation is still an O(N^2) problem for systems of size N. We show how this can be reduced to… ▽ More We discuss an efficient approach to the calculation of the internal energy in numerical simulations of spin systems with long-range interactions. Although, since the introduction of the Luijten-Blöte algorithm, Monte Carlo simulations of these systems no longer pose a fundamental problem, the energy calculation is still an O(N^2) problem for systems of size N. We show how this can be reduced to an O(N logN) problem, with a break-even point that is already reached for very small systems. This allows the study of a variety of, until now hardly accessible, physical aspects of these systems. In particular, we combine the optimized energy calculation with histogram interpolation methods to investigate the specific heat of the Ising model and the first-order regime of the three-state Potts model with long-range interactions. △ Less

Submitted 28 October, 1999; originally announced October 1999.

Comments: 10 pages, including 8 EPS figures. To appear in Phys. Rev. E. Also available as PDF file at http://www.cond-mat.physik.uni-mainz.de/~luijten/erikpubs.html

Journal ref: Phys. Rev. E 61, 2058 (2000)

Monte Carlo investigations of phase transitions: status and perspectives

Authors: Kurt Binder, Erik Luijten, Marcus Müller, Nigel B. Wilding, Henk W. J. Blöte

Abstract: Using the concept of finite-size scaling, Monte Carlo calculations of various models have become a very useful tool for the study of critical phenomena, with the system linear dimension as a variable. As an example, several recent studies of Ising models are discussed, as well as the extension to models of polymer mixtures and solutions. It is shown that using appropriate cluster algorithms, eve… ▽ More Using the concept of finite-size scaling, Monte Carlo calculations of various models have become a very useful tool for the study of critical phenomena, with the system linear dimension as a variable. As an example, several recent studies of Ising models are discussed, as well as the extension to models of polymer mixtures and solutions. It is shown that using appropriate cluster algorithms, even the scaling functions describing the crossover from the Ising universality class to the mean-field behavior with increasing interaction range can be described. Additionally, the issue of finite-size scaling in Ising models above the marginal dimension (d*=4) is discussed. △ Less

Submitted 19 August, 1999; originally announced August 1999.

Comments: 23 pages, including 14 PostScript figures. Presented at StatPhys-Taiwan, August 9-16, 1999. Also available as PDF file at http://www.cond-mat.physik.uni-mainz.de/~luijten/erikpubs.html

Journal ref: Physica A 281, 112 (2000)

Monte Carlo simulation of spin models with long-range interactions

Authors: Erik Luijten

Abstract: An efficient Monte Carlo algorithm for the simulation of spin models with long-range interactions is discussed. Its central feature is that the number of operations required to flip a spin is independent of the number of interactions between this spin and the other spins in the system. In addition, critical slowing down is strongly suppressed. In order to illustrate the range of applicability of… ▽ More An efficient Monte Carlo algorithm for the simulation of spin models with long-range interactions is discussed. Its central feature is that the number of operations required to flip a spin is independent of the number of interactions between this spin and the other spins in the system. In addition, critical slowing down is strongly suppressed. In order to illustrate the range of applicability of the algorithm, two specific examples are presented. First, some aspects of the Kosterlitz--Thouless transition in the one-dimensional Ising chain with inverse-square interactions are calculated. Secondly, the crossover from Ising-like to classical critical behavior in two-dimensional systems is studied for several different interaction profiles. △ Less

Submitted 15 June, 1999; originally announced June 1999.

Comments: 14 pages + 4 PostScript figures. To appear in Computer Simulation Studies in Condensed Matter Physics XII, edited by D. P. Landau, S. P. Lewis and H. B. Schuettler (Springer, Heidelberg, 1999). Also available as PDF file at http://www.cond-mat.physik.uni-mainz.de/~luijten/erikpubs.html

Do crossover functions depend on the shape of the interaction profile?

Authors: Erik Luijten, Kurt Binder

Abstract: We examine the crossover from classical to non-classical critical behaviour in two-dimensional systems with a one-component order parameter. Since the degree of universality of the corresponding crossover functions is still subject to debate, we try to induce non-universal effects by adding interactions with a second length scale. Although the crossover functions clearly depend on the range of t… ▽ More We examine the crossover from classical to non-classical critical behaviour in two-dimensional systems with a one-component order parameter. Since the degree of universality of the corresponding crossover functions is still subject to debate, we try to induce non-universal effects by adding interactions with a second length scale. Although the crossover functions clearly depend on the range of the interactions, they turn out to be remarkably robust against further variation of the interaction profile. In particular, we find that the earlier observed non-monotonic crossover of the effective susceptibility exponent occurs for several qualitatively different shapes of this profile. △ Less

Submitted 15 June, 1999; originally announced June 1999.

Comments: 7 pages + 4 PostScript figures. Accepted for publication in Europhysics Letters. Also available as PDF file at http://www.cond-mat.physik.uni-mainz.de/~luijten/erikpubs.html

Journal ref: Europhys. Lett. 47, 311 (1999)

Finite-size scaling above the upper critical dimension revisited: The case of the five-dimensional Ising model

Authors: Erik Luijten, Kurt Binder, Henk W. J. Blöte

Abstract: Monte Carlo results for the moments <M^k> of the magnetization distribution of the nearest-neighbor Ising ferromagnet in a L^d geometry, where L (4 \leq L \leq 22) is the linear dimension of a hypercubic lattice with periodic boundary conditions in d=5 dimensions, are analyzed in the critical region and compared to a recent theory of Chen and Dohm (CD) [X.S. Chen and V. Dohm, Int. J. Mod. Phys.… ▽ More Monte Carlo results for the moments <M^k> of the magnetization distribution of the nearest-neighbor Ising ferromagnet in a L^d geometry, where L (4 \leq L \leq 22) is the linear dimension of a hypercubic lattice with periodic boundary conditions in d=5 dimensions, are analyzed in the critical region and compared to a recent theory of Chen and Dohm (CD) [X.S. Chen and V. Dohm, Int. J. Mod. Phys. C (1998)]. We show that this finite-size scaling theory (formulated in terms of two scaling variables) can account for the longstanding discrepancies between Monte Carlo results and the so-called ``lowest-mode'' theory, which uses a single scaling variable tL^{d/2} where t=T/T_c-1 is the temperature distance from the critical temperature, only to a very limited extent. While the CD theory gives a somewhat improved description of corrections to the ``lowest-mode'' results (to which the CD theory can easily be reduced in the limit t \to 0, L \to \infty, tL^{d/2} fixed) for the fourth-order cumulant, discrepancies are found for the susceptibility (L^d <M^2>). Reasons for these problems are briefly discussed. △ Less

Submitted 7 January, 1999; originally announced January 1999.

Comments: 9 pages, 13 Encapsulated PostScript figures. To appear in Eur. Phys. J. B. Also available as PDF file at http://www.cond-mat.physik.uni-mainz.de/~luijten/erikpubs.html

Journal ref: Eur. Phys. J. B 9, 289 (1999)

Critical properties of the three-dimensional equivalent-neighbor model and crossover scaling in finite systems

Authors: Erik Luijten

Abstract: Accurate numerical results are presented for the three-dimensional equivalent-neighbor model on a cubic lattice, for twelve different interaction ranges (coordination number between 18 and 250). These results allow the determination of the range dependences of the critical temperature and various critical amplitudes, which are compared to renormalization-group predictions. In addition, the analy… ▽ More Accurate numerical results are presented for the three-dimensional equivalent-neighbor model on a cubic lattice, for twelve different interaction ranges (coordination number between 18 and 250). These results allow the determination of the range dependences of the critical temperature and various critical amplitudes, which are compared to renormalization-group predictions. In addition, the analysis yields an estimate for the interaction range at which the leading corrections to scaling vanish for the spin-1/2 model and confirms earlier conclusions that the leading Wegner correction must be negative for the three-dimensional (nearest-neighbor) Ising model. By complementing these results with Monte Carlo data for systems with coordination numbers as large as 52514, the full finite-size crossover curves between classical and Ising-like behavior are obtained as a function of a generalized Ginzburg parameter. Also the crossover function for the effective magnetic exponent is determined. △ Less

Submitted 23 February, 1999; v1 submitted 24 November, 1998; originally announced November 1998.

Comments: Corrected shift of critical temperature and some typos. To appear in Phys. Rev. E. 18 pages RevTeX, including 10 EPS figures. Also available as PDF file at http://www.cond-mat.physik.uni-mainz.de/~luijten/erikpubs.html

Report number: KOMA-98-25

Journal ref: Phys. Rev. E 59, 4997 (1999)

Shape of crossover between mean-field and asymptotic critical behavior in a three-dimensional Ising lattice

Authors: M. A. Anisimov, E. Luijten, V. A. Agayan, J. V. Sengers, K. Binder

Abstract: Recent numerical studies of the susceptibility of the three-dimensional Ising model with various interaction ranges have been analyzed with a crossover model based on renormalization-group matching theory. It is shown that the model yields an accurate description of the crossover function for the susceptibility. Recent numerical studies of the susceptibility of the three-dimensional Ising model with various interaction ranges have been analyzed with a crossover model based on renormalization-group matching theory. It is shown that the model yields an accurate description of the crossover function for the susceptibility. △ Less

Submitted 31 August, 1999; v1 submitted 20 October, 1998; originally announced October 1998.

Comments: 4 pages RevTeX + 3 PostScript figures. Uses epsf.sty and rotate.sty. Final version; accepted for publication in Physics Letters A

Journal ref: Phys. Lett. A 264, 63 (1999)

Nature of crossover from classical to Ising-like critical behavior

Authors: Erik Luijten, Kurt Binder

Abstract: We present an accurate numerical determination of the crossover from classical to Ising-like critical behavior upon approach of the critical point in three-dimensional systems. The possibility to vary the Ginzburg number in our simulations allows us to cover the entire crossover region. We employ these results to scrutinize several semi-phenomenological crossover scaling functions that are widel… ▽ More We present an accurate numerical determination of the crossover from classical to Ising-like critical behavior upon approach of the critical point in three-dimensional systems. The possibility to vary the Ginzburg number in our simulations allows us to cover the entire crossover region. We employ these results to scrutinize several semi-phenomenological crossover scaling functions that are widely used for the analysis of experimental results. In addition we present strong evidence that the exponent relations do not hold between effective exponents. △ Less

Submitted 31 July, 1998; originally announced July 1998.

Comments: 4 pages RevTeX 3.0/3.1, 4 Encapsulated PostScript figures. Uses epsf.sty. Also available as PDF file at http://www.cond-mat.physik.uni-mainz.de/~luijten/erikpubs.html

Report number: KOMA-98-12

Journal ref: Phys. Rev. E 58, R4060 (1998); 59, 7254 (1999)

Crossover scaling in two dimensions

Authors: Erik Luijten, Henk W. J. Blöte, Kurt Binder

Abstract: We determine the scaling functions describing the crossover from Ising-like critical behavior to classical critical behavior in two-dimensional systems with a variable interaction range. Since this crossover spans several decades in the reduced temperature as well as in the finite-size crossover variable, it has up to now largely evaded a satisfactory numerical determination. Using a new Monte C… ▽ More We determine the scaling functions describing the crossover from Ising-like critical behavior to classical critical behavior in two-dimensional systems with a variable interaction range. Since this crossover spans several decades in the reduced temperature as well as in the finite-size crossover variable, it has up to now largely evaded a satisfactory numerical determination. Using a new Monte Carlo method, we could obtain accurate results for sufficiently large interactions ranges. Our data cover the full crossover region both above and below the critical temperature and support the hypothesis that the crossover functions are universal. Also the so-called effective exponents are discussed and we show that these can vary nonmonotonically in the crossover region. △ Less

Submitted 25 June, 1997; originally announced June 1997.

Comments: 24 pages RevTeX 3.0/3.1, including 22 PostScript figures. Uses epsf.sty

Journal ref: Phys. Rev. E 56, 6540 (1997)

Nonmonotonical crossover of the effective susceptibility exponent

Authors: Erik Luijten, Henk W. J. Blöte, Kurt Binder

Abstract: We have numerically determined the behavior of the magnetic susceptibility upon approach of the critical point in two-dimensional spin systems with an interaction range that was varied over nearly two orders of magnitude. The full crossover from classical to Ising-like critical behavior, spanning several decades in the reduced temperature, could be observed. Our results convincingly show that th… ▽ More We have numerically determined the behavior of the magnetic susceptibility upon approach of the critical point in two-dimensional spin systems with an interaction range that was varied over nearly two orders of magnitude. The full crossover from classical to Ising-like critical behavior, spanning several decades in the reduced temperature, could be observed. Our results convincingly show that the effective susceptibility exponent gamma_eff changes nonmonotonically from its classical to its Ising value when approaching the critical point in the ordered phase. In the disordered phase the behavior is monotonic. Furthermore the hypothesis that the crossover function is universal is supported. △ Less

Submitted 19 June, 1997; originally announced June 1997.

Comments: 4 pages RevTeX 3.0/3.1, 5 Encapsulated PostScript figures. Uses epsf.sty. Accepted for publication in Physical Review Letters. Also available as PostScript and PDF file at http://www.tn.tudelft.nl/tn/erikpubs.html

Journal ref: Phys. Rev. Lett. 79, 561 (1997)

Classical critical behavior of spin models with long-range interactions

Authors: Erik Luijten, Henk W. J. Blöte

Abstract: We present the results of extensive Monte Carlo simulations of Ising models with algebraically decaying ferromagnetic interactions in the regime where classical critical behavior is expected for these systems. We corroborate the values for the exponents predicted by renormalization theory for systems in one, two, and three dimensions and accurately observe the predicted logarithmic corrections a… ▽ More We present the results of extensive Monte Carlo simulations of Ising models with algebraically decaying ferromagnetic interactions in the regime where classical critical behavior is expected for these systems. We corroborate the values for the exponents predicted by renormalization theory for systems in one, two, and three dimensions and accurately observe the predicted logarithmic corrections at the upper critical dimension. We give both theoretical and numerical evidence that above the upper critical dimension the decay of the critical spin-spin correlation function in finite systems consists of two different regimes. For one-dimensional systems our estimates for the critical couplings are more than two orders of magnitude more accurate than existing estimates. In two and three dimensions we give, to our knowledge, the first results for the critical couplings. △ Less

Submitted 14 May, 1997; v1 submitted 2 April, 1997; originally announced April 1997.

Comments: 20 pages RevTeX, 9 PostScript figures. Uses epsf.sty. Corrected scaling analysis at the upper critical dimension, in particular Eqs. (13a) and (13b). Also available as PostScript and PDF file at http://www.tn.tudelft.nl/tn/erikpubs.html

Journal ref: Phys. Rev. B 56, 8945 (1997)

Universality and the five-dimensional Ising model

Authors: Henk W. J. Blöte, Erik Luijten

Abstract: We solve the long-standing discrepancy between Monte Carlo results and the renormalization prediction for the Binder cumulant of the five-dimensional Ising model. Our conclusions are based on accurate Monte Carlo data for systems with linear sizes up to L=22. A detailed analysis of the corrections to scaling allows the extrapolation of these results to L=\infinity. Our determination of the criti… ▽ More We solve the long-standing discrepancy between Monte Carlo results and the renormalization prediction for the Binder cumulant of the five-dimensional Ising model. Our conclusions are based on accurate Monte Carlo data for systems with linear sizes up to L=22. A detailed analysis of the corrections to scaling allows the extrapolation of these results to L=\infinity. Our determination of the critical point, K_c=0.1139150 (4), is more than an order of magnitude more accurate than previous estimates. △ Less

Submitted 19 February, 1997; originally announced February 1997.

Comments: 6 pages LaTeX, 1 PostScript figure. Uses cite.sty (included) and epsf.sty. Also available as PostScript and PDF file at http://www.tn.tudelft.nl/tn/erikpubs.html

Journal ref: Europhys. Lett. 38, 565 (1997)