BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Gerhard Kahl

Gerhard Kahl appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
cond-mat.softcond-mat.stat-mechphysics.comp-phcond-mat.dis-nnphysics.chem-ph

Catalog footprint

What is connected

28works
5topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

28 published item(s)

preprint2021arXiv

Deformable hard particles particles confined in a disordered porous matrix

With suitably designed Monte Carlo simulations we have investigated the properties of mobile, impenetrable, yet deformable particles that are immersed into a porous matrix, the latter one realized via a frozen configuration of spherical particles. By virtue of a model put forward by Batista and Miller [Phys. Rev. Lett. {\bf 105}, 088305 (2010)] the fluid particles can change under the impact of their surrounding (i.e., either other fluid particles or the matrix) their shape within the class of ellipsoids of revolution; such a change in shape is related to an energy change which is fed into suitably defined selection rules in the deformation "moves" of the Monte Carlo simulations. This concept represents a simple, yet powerful model of realistic, deformable molecules with complex internal structures (such as dendrimers or polymers). For the evaluation of the properties of the system we have used the well-known quenched-annealed protocol (with its characteristic double average prescription) and have analysed the simulation data in terms of static properties (radial distribution function and aspect ratio distribution of the ellipsoids) and dynamic features (notably the mean squared displacement). Our data provide evidence that the degree of deformability of the fluid particles has a distinct impact on the aforementioned properties of the system.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Elasticity in crystals with high density of local defects : insights from ultra-soft colloids

In complex crystals close to melting or at finite temperatures, different types of defects are ubiquitous and their role becomes relevant in the mechanical response of these solids. Conventional elasticity theory fails to provide a microscopic basis to include and account for the motion of point-defects in an otherwise ordered crystalline structure. We study the elastic properties of a point-defect rich crystal within a first-principles theoretical framework derived from microscopic equations of motion. This framework allows us to make specific predictions pertaining to the mechanical properties which we can validate through deformation experiments performed in Molecular Dynamics simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Inverse patchy colloids: synthesis, modeling and self-organization

Inverse patchy colloids are nano- to micro-scale particles with a surface divided into differently charged regions. This class of colloids combines directional, selective bonding with a relatively simple particle design: owing to the competitive interplay between the orientation-dependent attraction and repulsion -- induced by the interactions between like/oppositely charged areas -- experimentally accessible surface patterns are complex enough to favor the stabilization of specific structures of interest. Most important, the behavior of heterogeneously charged units can be ideally controlled by means of external parameters, such as the pH and the salt concentration. We present a concise review about this class of systems, spanning the range from the synthesis of model inverse patchy particles to their self-assembly, covering their coarse-grained modeling and the related numerical/analytical treatments.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

On the yielding of a defect-rich model crystal under shear: insights from molecular dynamics simulations

Point defects in real crystals at finite temperatures are inevitable. Their dynamics severely influence the mechanical properties of crystals under shear giving rise to nonlinear effects such as ductility. Therefore, it is crucial to explore the interplay of the equilibrium point-defect diffusion timescales and shear-induced timescales to understand the plastic deformation of crystals. Using extensive nonequilibrium molecular dynamics simulations, we present a study on the yielding behavior of cluster crystals (CC), an archetypal model for defect-rich crystals where the crystalline structure is characterized by multiple occupancies (cluster) of particles at a lattice site. In equilibrium, particles diffuse via site-to-site hopping keeping the crystalline structure intact. We consider the CCs at a fixed density and different temperatures where it remains in the FCC structure, and the diffusion timescales of particles vary depending on the temperature. We choose the range of shear rates, which covers timescales higher and much lower than the equilibrium diffusion timescales at high temperatures. For the considered range of shear rates and temperatures, both the macro and microscopic responses of CCs to the shear suggest that the yielding scenario remains independent of the diffusion of particles. It involves the plastic deformation of the underlying crystalline structure. The averaged local bond order parameters and local angle measurements demonstrate the structural changes and cooperative movement of the center of masses of the clusters close to the yield point. A comparison with the soft-sphere (SS) FCC crystal reveals the similarities in the yielding behavior of both systems. Nonetheless, diffusion of particles influences certain features, such as a less prominent increase in the local bond order parameters and local angles close to the yield point.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Hydrodynamic correlations of viscoelastic fluids by multiparticle collision dynamics simulations

The emergent fluctuating hydrodynamics of a viscoelastic fluid modeled by the multiparticle collision dynamics (MPC) approach is studied. The fluid is composed of flexible, Gaussian phantom polymers, which interact by local momentum-conserving stochastic MPC collisions. For comparison, the analytical solution of the linearized Navier-Stokes equation is calculated, where viscoelasticity is taken into account by a time-dependent shear relaxation modulus. The fluid properties are characterized by the transverse velocity autocorrelation function in Fourier space as well as in real space. Various polymer lengths are considered---from dumbbells to (near-)continuous polymers. Viscoelasticity affects the fluid properties and leads to strong correlations, which overall decay exponentially in Fourier space. In real space, the center-of-mass velocity autocorrelation function of individual polymers exhibits a long-time tail independent of polymer length, which decays as $t^{-3/2}$, similar to a Newtonian fluid, in the asymptotic limit $t \to \infty$. Moreover, for long polymers an additional power-law decay appears at time scales shorter than the longest polymer relaxation time with the same time dependence, but negative correlations, and the polymer length dependence $L^{-1/2}$. Good agreement is found between the analytical and simulation results.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

On the degeneracy of ordered ground state configurations of the aspherical Gaussian core model

We provide rigorous evidence that the ordered ground state configurations of a system of parallel oriented, ellipsoidal particles, interacting via a Gaussian interaction (termed in literature as Gaussian core nematics) {\it must} be infinitely degenerate: we have demonstrated that these configurations originate from the related ground state configuration of the corresponding symmetric Gaussian core system via a suitable stretching operation of this lattice in combination with an arbitrary rotation. These findings explain related observations in former investigations, which then remained unexplained. Our conclusions have far reaching consequences for the search of ground state configurations of other nematic particles.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Reliable computational prediction of supramolecular ordering of complex molecules under electrochemical conditions

We propose a computationally lean, two-stage approach that reliably predicts self-assembly behavior of complex charged molecules on a metallic surfaces under electrochemical conditions. Stage one uses ab initio simulations to provide reference data for the energies (evaluated for archetypical configurations) to fit the parameters of a conceptually much simpler and computationally less expensive model of the molecules: classical, spherical particles, representing the respective atomic entities, a soft but perfectly conductive wall potential represents the metallic surface. Stage two feeds the energies that emerge from this classical model into highly efficient and reliable optimization techniques to identify via energy minimization the ordered ground state configurations of the molecules. We demonstrate the power of our approach by successfully reproducing, on a semi-quantitative level, the intricate supramolecular ordering observed experimentally for PQP$^+$ and ClO$_4^-$ molecules at an Au(111)-electrolyte interface, including the formation of open-porous, self-hosts--guest, and stratified bilayer phases as a function of the electric field at the solid--liquid interface. We also discuss the role of the perchlorate ions in the self-assembly process, whose positions could not be identified in the related experimental investigations.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Equilibrium structures of anisometric, quadrupolar particles confined to a monolayer

We investigate the structural properties of a two-dimensional system of ellipsoidal particles carrying a linear quadrupole moment in their center. These particles represent a simple model for a variety of uncharged, non-polar conjugated organic molecules. Using optimization tools based on ideas of Evolutionary Algorithms, we first examine the ground state structures as we vary the aspect ratio of the particles and the pressure. Interestingly, we find, besides the intuitively expected T-like configurations, a variety of complex structures, characterized with up to three different particle orientations. In an effort to explore the impact of thermal fluctuations, we perform constant-pressure Molecular Dynamics simulations within a range of rather low temperatures. We observe that ground state structures formed by particles with a large aspect ratio are in particular suited to withstand fluctuations up to rather high temperatures. Our comprehensive investigations allow for a deeper understanding of molecular or colloidal monolayer arrangements under the influence of a typical electrostatic interaction on a coarse-grained level.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Generalized inverse patchy colloid model

We generalize the inverse patchy colloid model that was originally developed for heterogeneously charged particles with two identical polar patches and an oppositely charged equator to a model that can have a considerably richer surface pattern. Based on a Debye-Hueckel framework, we propose a coarse-grained description of the effective pair interactions that is applicable to particles with an arbitrary patch decoration. We demonstrate the versatility of this approach by applying it to models with (i) two differently charged and/or sized patches, and (ii) three, possibly different patches.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Self-assembly of DNA-functionalized colloids

Colloidal particles grafted with single-stranded DNA (ssDNA) chains can self-assemble into a number of different crystalline structures, where hybridization of the ssDNA chains creates links between colloids stabilizing their structure. Depending on the geometry and the size of the particles, the grafting density of the ssDNA chains, and the length and choice of DNA sequences, a number of different crystalline structures can be fabricated. However, understanding how these factors contribute synergistically to the self-assembly process of DNA-functionalized nano- or micro-sized particles remains an intensive field of research. Moreover, the fabrication of long-range structures due to kinetic bottlenecks in the self-assembly are additional challenges. Here, we discuss the most recent advances from theory and experiment with particular focus put on recent simulation studies.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Antinematic local order in dendrimer liquids

We use monomer-resolved numerical simulations to study the positional and orientational structure of a dense dendrimer solution, focusing on the effects of dendrimers' prolate shape and deformability on the short-range order. Our results provide unambiguous evidence that the nearest-neighbor shell of a tagged particle consists of a mixture of crossed, side-by-side, side-to-end, and end-to-end pair configurations, imposing antinematic rather than nematic order observed in undeformable rodlike particles. This packing pattern persists even at densities where particle overlap becomes sizable. We demonstrate that the antinematic arrangement is compatible with the A15 crystal lattice reported in several dendrimer compounds.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Crystal phases of soft spheres systems in a slab geometry

We have identified the ground state configurations of soft particles (interacting via inverse power potentials) confined between two hard, impenetrable walls. To this end we have used a highly reliable optimization scheme at {\it vanishing} temperature while varying the wall separation over a representative range. Apart from the expected layered triangular and square structures (which are compatible with the three dimensional bulk fcc lattice), we have identified a cascade of highly complex intermediate structures. Taking benefit of the general scaling properties of inverse power potentials, we could identify -- for a given softness value -- one single master curve which relates the energy to the wall separation, irrespective of the density of the system. Via extensive Monte Carlo simulations, we have performed closer investigations of these intermediate structures at {\it finite} temperature: we could provide evidence to which extent these particle arrangements remain stable over a relatively large temperature range.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Explicit spatial description of fluid inclusions in porous matrices in terms of an inhomogeneous integral equation

We study the fluid inclusion of both Lennard-Jones particles and particles with competing interaction ranges --short range attractive and long range repulsive (SALR)-- in a disordered porous medium constructed as a controlled pore glass in two dimensions. With the aid of a full two-dimensional Ornstein-Zernike approach, complemented by a Replica Ornstein-Zernike integral equation, we explicitly obtain the spatial density distribution of the fluid adsorbed in the porous matrix and a good approximation for the average fluid-matrix correlations. The results illustrate the remarkable differences between the adsorbed Lennard-Jones (LJ) and SALR systems. In the latter instance, particles tend to aggregate in clusters which occupy pockets and bays in the porous structure, whereas the LJ fluid uniformly wets the porous walls. A comparison with Molecular Dynamics simulations shows that the two-dimensional Ornstein-Zernike approach with a Hypernetted Chain closure together with a sensible approximation for the fluid-fluid correlations can provide an accurate picture of the spatial distribution of adsorbed fluids for a given configuration of porous material.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Inverse patchy colloids with small patches: fluid structure and dynamical slowing down

Inverse Patchy Colloids (IPCs) differ from conventional patchy particles because their patches repel (rather than attract) each other and attract (rather than repel) the part of the colloidal surface that is free of patches. These particular features occur, .e.g., in heterogeneously charged colloidal systems. Here we consider overall neutral IPCs carrying two, relatively small, polar patches. Previous studies of the same model under planar confinement have evidenced the formation of branched, disordered aggregates composed of ring-like structures. We investigate here the bulk behavior of the system via molecular dynamics simulations, focusing on both the structure and the dynamics of the fluid phase in a wide region of the phase diagram. Additionally, the simulation results for the static observables are compared to the Associative Percus Yevick solution of an integral equation approach based on the multi-density Ornstein-Zernike theory. A good agreement between theoretical and numerical quantities is observed even in the region of the phase diagram where the slowing down of the dynamics occurs.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Taking one charge off a two-dimensional Wigner crystal

A planar array of identical charges at vanishing temperature forms a Wigner crystal with hexagonal symmetry. We take off one (reference) charge in a perpendicular direction, hold it fixed, and search for the ground state of the whole system. The planar projection of the reference charge should then evolve from a six-fold coordination (center of a hexagon) for small distances to a three-fold arrangement (center of a triangle), at large distances $d$ from the plane. The aim of this paper is to describe the corresponding non-trivial lattice transformation. For that purpose, two numerical methods (direct energy minimization and Monte Carlo simulations), together with an analytical treatment, are presented. Our results indicate that the $d=0$ and $d\to\infty$ limiting cases extend for finite values of $d$ from the respective starting points into two sequences of stable states, with intersecting energies at some value $d_t$; beyond this value the branches continue as metastable states.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Theoretical and numerical investigations of inverse patchy colloids in the fluid phase

We investigate the structural and thermodynamic properties of a new class of patchy colloids, referred to as inverse patchy colloids (IPCs) in their fluid phase via both theoretical methods and simulations. IPCs are nano- or micro- meter sized particles with differently charged surface regions. We extend conventional integral equation schemes to this particular class of systems: our approach is based on the so-called multi-density Ornstein-Zernike equation, supplemented by the associative Percus-Yevick approximation (APY). To validate the accuracy of our framework, we compare the obtained results with data extracted from $NpT$ and $NVT$ Monte Carlo simulations. In addition, other theoretical approaches are used to calculate the properties of the system: the reference hypernetted-chain (RHNC) method and the Barker-Henderson thermodynamic perturbation theory. Both APY and RHNC frameworks provide accurate predictions for the pair distribution functions: APY results are in slightly better agreement with MC data, in particular at lower temperatures where the RHNC solution does not converge.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Computing Gibbs free energy differences by interface pinning

We propose an approach for computing the Gibbs free energy difference between phases of a material. The method is based on the determination of the average force acting on interfaces that separate the two phases of interest. This force, which depends on the Gibbs free energy difference between the phases, is computed by applying an external harmonic field that couples to a parameter which specifies the two phases. Validated first for the Lennard-Jones model, we demonstrate the flexibility, efficiency and practical applicability of this approach by computing the melting temperatures of sodium, magnesium, aluminum and silicon at ambient pressure using density functional theory. Excellent agreement with experiment is found for all four elements, except for silicon, for which the melting temperature is, in agreement with previous simulations, seriously underestimated.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Wigner crystals for a planar, equimolar binary mixture of classical, charged particles

We have investigated the ground state configurations of an equimolar, binary mixture of classical charged particles (with nominal charges $Q_1$ and $Q_2$) that condensate on a neutralizing plane. Using efficient Ewald summation techniques for the calculation of the ground state energies, we have identified the energetically most favourable ordered particle arrangements with the help of a highly reliable optimization tool based on ideas of evolutionary algorithms. Over a large range of charge ratios, $q = Q_2 / Q_1$, we identify six non-trivial ground states, some of which show a remarkable and unexpected structural complexity. For $0.59 \lesssim q < 1$ the system undergoes a phase separation where the two charge species populate in a hexagonal arrangement spatially separated areas.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Predicting Patchy Particle Crystals: Variable Box Shape Simulations and Evolutionary Algorithms

We consider several patchy particle models that have been proposed in literature and we investigate their candidate crystal structures in a systematic way. We compare two different algorithms for predicting crystal structures: (i) an approach based on Monte Carlo simulations in the isobaric-isothermal ensemble and (ii) an optimization technique based on ideas of evolutionary algorithms. We show that the two methods are equally successful and provide consistent results on crystalline phases of patchy particle systems.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Clustering, conductor-insulator transition and phase separation of an ultrasoft model of electrolytes

We investigate the clustering and phase separation of a model of ultrasoft, oppositely charged macroions by a combination of Monte Carlo and Molecular Dynamics simulations. Static and dynamic diagnostics, including the dielectric permittivity and the electric conductivity of the model, show that ion pairing induces a sharp conductor-insulator transition at low temperatures and densities, which impacts the separation into dilute and concentrated phases below a critical temperature. Preliminary evidence is presented for a possible tricritical nature of the phase diagram of the model.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Dynamic arrest of colloids in porous environments: disentangling crowding and confinement

Using numerical simulations we study the slow dynamics of a colloidal hard-sphere fluid adsorbed in a matrix of disordered hard-sphere obstacles. We calculate separately the contributions to the single-particle dynamic correlation functions due to free and trapped particles. The separation is based on a Delaunay tessellation to partition the space accessible to the centres of fluid particles into percolating and disconnected voids. We find that the trapping of particles into disconnected voids of the matrix is responsible for the appearance of a nonzero long-time plateau in the single-particle intermediate scattering functions of the full fluid. The subdiffusive exponent $z$, obtained from the logarithmic derivative of the mean-squared displacement, is observed to be essentially unaffected by the motion of trapped particles: close to the percolation transition, we determined $z \simeq 0.5$ for both the full fluid and the particles moving in the percolating void. Notably, the same value of $z$ is found in single-file diffusion and is also predicted by mode-coupling theory along the diffusion-localisation line. We also reveal subtle effects of dynamic heterogeneity in both the free and the trapped component of the fluid particles, and discuss microscopic mechanisms that contribute to this phenomenon.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Ultrasoft primitive model of polyionic solutions: structure, aggregation, and dynamics

We introduce an ultrasoft core model of interpenetrating polycations and polyanions with continuous Gaussian charge distributions, to investigate polyelectrolyte aggregation in dilute and semi-dilute, salt-free solutions. The model is studied by a combination of approximate theories (random phase approximation and hypernetted chain theory) and numerical simulations. The calculated pair structure, thermodynamics, phase diagram and polyion dynamics of the symmetric version of the model (the "ultrasoft restricted primitive model" or UPRM) differ from the corresponding properties of the widely studied "restricted primitive model" (RPM) where ions have hard cores. At sufficiently low temperatures and densities, oppositely charged polyions form weakly interacting, polarizable neutral pairs. The clustering probabilities, dielectric behavior and electrical conductivity point to a line of sharp conductor-insulator transitions in the density-temperature plane. At very low temperatures the conductor-insulator transition line terminates near the top of a first order coexistence curve separating a high-density, liquid phase from a low-density, vapor phase. The simulation data hint at a tricritical behavior, reminiscent of that observed of the two-dimensional Coulomb Gas, which contrasts with the Ising criticality of its three-dimensional counterpart, the RPM.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Hopping and microscopic dynamics of ultrasoft particles in cluster crystals

We have investigated the slow dynamics of ultrasoft particles in crystalline cluster phases, where point particles interact through the generalized exponential potential u(r) = ε\exp[-(r/σ)^n], focusing on the cluster fcc phase of this model with n=4. In an effort to elucidate how the mechanisms of mass transport depend on the microscopic dynamics and in order to mimic a realistic scenario in a related experiment we have performed molecular dynamics, Brownian dynamics, and Monte Carlo simulations. In molecular dynamics simulations the diffusion of particles proceeds through long-range jumps, guided by strong correlations in the momentum direction. In Monte Carlo and Brownian dynamics simulations jump events are short-ranged, reflecting the purely configurational nature of the dynamics. In contrast to what was found in models of glass-forming liquids, the effect of Newtonian and stochastic microscopic dynamics on the long-time relaxation cannot be accounted for by a temperature-independent rescaling of the time units. From the obvious qualitative discrepancies in the short time behavior between the three simulation methods and the non-trivial difference in jump length distributions, long time relaxation, and dynamic heterogeneity, we learn that a more complex modeling of the dynamics in realistic systems of ultrasoft colloids is required.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Impact of random obstacles on the dynamics of a dense colloidal fluid

Using molecular dynamics simulations we study the slow dynamics of a colloidal fluid annealed within a matrix of obstacles quenched from an equilibrated colloidal fluid. We choose all particles to be of the same size and to interact as hard spheres, thus retaining all features of the porous confinement while limiting the control parameters to the packing fraction of the matrix, Φm, and that of the fluid, Φf. We conduct detailed investigations on several dynamic properties, including the tagged-particle and collective intermediate scattering functions, the mean-squared displacement, and the van Hove function. We show the confining obstacles to profoundly impact the relaxation pattern of various quantifiers pertinent to the fluid. Varying the type of quantifier (tagged-particle or collective) as well as Φm and Φf, we unveil both discontinuous and continuous arrest scenarios. Furthermore, we discover subdiffusive behavior and demonstrate its close connection to the matrix structure. Our findings partly confirm the various predictions of a recent extension of mode-coupling theory to the quenched-annealed protocol.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Monomer-resolved simulations of cluster-forming dendrimers

We present results of monomer-resolved Monte Carlo simulations for a system of amphiphilic dendrimers of second generation. Our investigations validate a coarse-grained level description based on the zero-density limit effective pair-interactions for low and intermediate densities, which predicted the formation of stable, finite aggregates in the fluid phase. Indeed, we find that these systems form a homogeneous fluid for low densities, which on increasing the density spontaneously transforms into a fluid of clusters of dendrimers. Although these clusters are roughly spherical in nature for intermediate densities, also more complex structures are detected for the highest densities considered.

0 citations0 reviewsNo code linkedOpen access
preprint2006arXiv

Liquid-vapor transition of systems with mean field universality class

We have considered a system where the interaction, v(r) = v_IS(r) + xi^2 v_MF(r), is given as a linear combination of two potentials, each of which being characterized with a well-defined critical behavior: for v_IS(r) we have chosen the potential of the restricted primitive model which is known to belong to the Ising 3D (IS) universality class, while for v_MF(r) we have considered a long-range interaction in the Kac-limit, displaying mean field (MF) behavior. We study the performance of two theoretical approaches and of computer simulations in the critical region for this particular system and give a detailed comparison between theories and simulation of the critical region and the location of the critical point. Both, theory and simulation give evidence that the system belongs to the MF universality class for any positive value of xi and that it shows only non-classical behavior for xi=0. While in this limiting case theoretical approaches are known to fail, we find good agreement for the critical properties between the theoretical approaches and the simulations for xi^2 larger than 0.05.

0 citations0 reviewsNo code linkedOpen access
preprint2006arXiv

Phase behaviour of a symmetrical binary fluid mixture

We have investigated the phase behaviour of a symmetrical binary fluid mixture for the situation where the chemical potentials $μ_1$ and $μ_2$ of the two species differ. Attention is focused on the set of interparticle interaction strengths for which, when $μ_1=μ_2$, the phase diagram exhibits both a liquid-vapor critical point and a tricritical point. The corresponding phase behaviour for the case $μ_1\neμ_2$ is investigated via integral-equation theory calculations within the mean spherical approximation (MSA), and grand canonical Monte Carlo (GCMC) simulations. We find that two possible subtypes of phase behaviour can occur, these being distinguished by the relationship between the critical lines in the full phase diagram in the space of temperature, density, and concentration. We present the detailed form of the phase diagram for both subtypes and compare with the results from GCMC simulations, finding good overall agreement. The scenario via which one subtype evolves into the other, is also studied, revealing interesting features.

0 citations0 reviewsNo code linkedOpen access
preprint2005arXiv

Thermodynamically self-consistent liquid state theories for systems with bounded potentials

The mean spherical approximation (MSA) can be solved semi-analytically for the Gaussian core model (GCM) and yields - rather surprisingly - exactly the same expressions for the energy and the virial equations. Taking advantage of this semi-analytical framework, we apply the concept of the self-consistent Ornstein-Zernike approximation (SCOZA) to the GCM: a state-dependent function K is introduced in the MSA closure relation which is determined to enforce thermodynamic consistency between the compressibility route and either the virial or energy route. Utilizing standard thermodynamic relations this leads to two different differential equations for the function K that have to be solved numerically. Generalizing our concept we propose an integro-differential-equation based formulation of the SCOZA which, although requiring a fully numerical solution, has the advantage that it is no longer restricted to the availability of an analytic solution for a particular system. Rather it can be used for an arbitrary potential and even in combination with other closure relations, such as a modification of the hypernetted chain approximation.

0 citations0 reviewsNo code linkedOpen access