Source author record

K. Binder

Explore this author’s publications, collaborators and research topics. This catalog record has not been claimed by its author.

ResearcherUnclaimed source record

cond-mat.soft cond-mat.stat-mech cond-mat.mtrl-sci cond-mat cond-mat.dis-nn cond-mat.mes-hall cond-mat.other cond-mat.str-el physics.chem-ph physics.flu-dyn

Catalog footprint

Research overview

20works

10topics

4close collaborators

Actions

Is this your work?

Open graph Browse papers

preprint2014arXiv

Magnetic order and electromagnon excitations in DyMnO3 studied by neutron scattering experiments

Magnetic order and excitations in multiferroic DyMnO3 were studied by neutron scattering experiments using a single crystal prepared with enriched 162Dy isotope. The ordering of Mn moments exhibits pronounced hysteresis arising from the interplay between Mn and Dy magnetism which possesses a strong impact on the ferroelectric polarization. The magnon dispersion resembles that reported for TbMnO3. We identify the excitations at the magnetic zone center and near the zone boundary in the b direction, which can possess electromagnon character. The lowest frequency of the zone-center magnons is in good agreement with a signal in a recent optical measurement so that this mode can be identified as the electromagnon coupled by the same Dzyaloshinski-Moriya interaction as the static multiferroic phase.

preprint2013arXiv

Mixtures of anisotropic and spherical colloids: Phase behavior, confinement, percolation phenomena and kinetics

Purely entropic systems such as suspensions of hard rods, platelets and spheres show rich phase behavior. Rods and platelets have successfully been used as models to predict the equilibrium properties of liquid crystals for several decades. Over the past years hard particle models have also been studied in the context of non-equilibrium statistical mechanics, in particular regarding the glass transition, jamming, sedimentation and crystallization. Recently suspensions of hard anisotropic particles also moved into the focus of materials scientists who work on conducting soft matter composites. An insulating polymer resin that is mixed with conductive filler particles becomes conductive when the filler percolates. In this context the mathematical topic of connectivity percolation finds an application in modern nano-technology. In this article, we briefly review recent work on the phase behavior, confinement effects, percolation transition and phase transition kinetics in hard particle models. In the first part, we discuss the effects that particle anisotropy and depletion have on the percolation transition. In the second part, we present results on the kinetics of the liquid-to-crystal transition in suspensions of spheres and of ellipsoids.

preprint2013arXiv

Star Polymers Confined in a Nanoslit: A Simulation Test of Scaling and Self-Consistent Field Theories

The free energy cost of confining a star polymer where $f$ flexible polymer chains containing $N$ monomeric units are tethered to a central unit in a slit with two parallel repulsive walls a distance $D$ apart is considered, for good solvent conditions. Also the parallel and perpendicular components of the gyration radius of the star polymer, and the monomer density profile across the slit are obtained. Theoretical descriptions via Flory theory and scaling treatments are outlined, and compared to numerical self-consistent field calculations (applying the Scheutjens-Fleer lattice theory) and to Molecular Dynamics results for a bead-spring model. It is shown that Flory theory and self-consistent field (SCF) theory yield the correct scaling of the parallel linear dimension of the star with $N$ , $f$ and $D$ , but cannot be used for estimating the free energy cost reliably. We demonstrate that the same problem occurs already for the confinement of chains in cylindrical tubes. We also briefly discuss the problem of a free or grafted star polymer interacting with a single wall, and show that the dependence of confining force on the functionality of the star is different for a star confined in a nanoslit and a star interacting with a single wall, which is due to the absence of a symmetry plane in the latter case.

preprint2013arXiv

The escape transition of a compressed star polymer: Self-consistent field predictions tested by simulation

The escape transition of a polymer "mushroom" (a flexible chain grafted to a flat non-adsorbing substrate surface in a good solvent) occurs when the polymer is compressed by a cylindrical piston of radius $R$ , that by far exceeds the chain gyration radius. At this transition, the chain conformation abruptly changes from a two-dimensional self-avoiding walk of blobs (of diameter $H$ , the height of the piston above the substrate) to a "flower conformation", i.e. stretched almost one-dimensional string of blobs (with end-to-end distance $\approx R$ ) and an "escaped" part of the chain, the "crown", outside the piston. The extension of this problem to the case of star polymers with $f$ arms is considered, assuming that the center of the star is grafted to the substrate. The question is considered whether under compression the arms escape all together, or whether there occurs an arm by arm escape under increasing compression. Both self-consistent field calculations and Molecular Dynamics simulations are found to favor the latter scenario.

preprint2012arXiv

Langevin Dynamics simulations of a 2-dimensional colloidal crystal under confinement and shear

Langevin Dynamics simulations are used to study the effect of shear on a two-dimensional colloidal crystal confined by structured parallel walls. When walls are sheared very slowly, only two or three crystalline layers next to the walls move along with them, while the inner layers of the crystal are only slightly tilted. At higher shear velocities, this inner part of the crystal breaks into several pieces with different orientations. The velocity profile across the slit is reminiscent of shear-banding in flowing soft materials, where liquid and solid regions coexist; the difference, however, is that in the latter case the solid regions are glassy while here they are crystalline. At even higher shear velocities, the effect of the shearing becomes smaller again. Also the effective temperature near the walls (deduced from the velocity distributions of the particles) decreases again when the wall velocity gets very large. When the walls are placed closer together, thereby introducing a misfit, a structure containing a soliton staircase arises in simulations without shear. Introducing shear increases the disorder in these systems until no solitons are visible any more. Instead, similar structures like in the case without misfit result. At high shear rates, configurations where the incommensurability of the crystalline structure is compensated by the creation of holes become relevant.

preprint2011arXiv

Microcanonical Determination of the Interface Tension of Flat and Curved Interfaces from Monte Carlo Simulations

The investigation of phase coexistence in systems with multi-component order parameters in finite systems is discussed, and as a generic example, Monte Carlo simulations of the two-dimensional q-state Potts model (q=30) on LxL square lattices (40<=L<=100) are presented. It is shown that the microcanonical ensemble is well-suited both to find the precise location of the first order phase transition and to obtain an accurate estimate for the interfacial free energy between coexisting ordered and disordered phases. For this purpose, a microcanonical version of the heatbath algorithm is implemented. The finite size behaviour of the loop in the curve describing the inverse temperature versus energy density is discussed, emphasizing that the extrema do not have the meaning of van der Waals-like "spinodal points" separating metastable from unstable states, but rather describe the onset of heterophase states: droplet/bubble evaporation/condensation transitions. Thus all parts of these loops, including the parts that correspond to a negative specific heat, describe phase coexistence in full thermal equilibrium. However, the estimates for the curvature-dependent interface tension of the droplets and bubbles suffer from unexpected and unexplained large finite size effects which need further study.

preprint2010arXiv

A slow process in confined polymer melts: layer exchange dynamics at a polymer solid interface

Employing Molecular Dynamics simulations of a chemically realistic model of 1,4-polybutadiene between graphite walls we show that the mass exchange between layers close to the walls is a slow process already in the melt state. For the glass transition of confined polymers this process competes with the slowing down due to packing effects and intramolecular rotation barriers.

preprint2010arXiv

Absorption/Expulsion of Oligomers and Linear Macromolecules in a Polymer Brush

The absorption of free linear chains in a polymer brush was studied with respect to chain size $L$ and compatibility $χ$ with the brush by means of Monte Carlo (MC) simulations and Density Functional Theory (DFT) / Self-Consistent Field Theory (SCFT) at both moderate, $σ_g = 0.25$ , and high, $σ_g = 1.00$ , grafting densities using a bead-spring model. Different concentrations of the free chains $0.0625 \le ϕ_o \le 0.375$ are examined. Contrary to the case of $χ= 0$ when all species are almost completely ejected by the polymer brush irrespective of their length $L$ , for $χ< 0$ we find that the degree of absorption (absorbed amount) $Γ(L)$ undergoes a sharp crossover from weak to strong ( $\approx 100%$ ) absorption, discriminating between oligomers, $1\le L\le 8$ , and longer chains. For a moderately dense brush, $σ_g = 0.25$ , the longer species, $L > 8$ , populate predominantly the deep inner part of the brush whereas in a dense brush $σ_g = 1.00$ they penetrate into the "fluffy" tail of the dense brush only. Gyration radius $R_g$ and end-to-end distance $R_e$ of absorbed chains thereby scale with length $L$ as free polymers in the bulk. Using both MC and DFT/SCFT methods for brushes of different chain length $32 \le N \le 256$ , we demonstrate the existence of unique {\em critical} value of compatibility $χ= χ^{c}<0$ . For $χ^{c}(ϕ_o)$ the energy of free chains attains the {\em same} value, irrespective of length $L$ whereas the entropy of free chain displays a pronounced minimum. At $χ^{c}$ all density profiles of absorbing chains with different $L$ intersect at the same distance from the grafting plane. The penetration/expulsion kinetics of free chains into the polymer brush after an instantaneous change in their compatibility $χ$ displays a rather rich behavior. We find three distinct regimes of penetration kinetics of free chains regarding the length $L$ : I ( $1\le L\le 8$ ), II ( $8 \le L \le N$ ), and III ( $L > N$ ), in which the time of absorption $τ$ grows with $L$ at a different rate. During the initial stages of penetration into the brush one observes a power-law increase of $Γ\propto t^α$ with power $α\propto -\ln ϕ_o$ whereby penetration of the free chains into the brush gets {\em slower} as their concentration rises.

preprint2010arXiv

Curvature Dependence of Surface Free Energy of Liquid Drops and Bubbles: A Simulation Study

We study the excess free energy due to phase coexistence of fluids by Monte Carlo simulations using successive umbrella sampling in finite LxLxL boxes with periodic boundary conditions. Both the vapor-liquid phase coexistence of a simple Lennard-Jones fluid and the coexistence between A-rich and B-rich phases of a symmetric binary (AB) Lennard-Jones mixture are studied, varying the density rho in the simple fluid or the relative concentration x_A of A in the binary mixture, respectively. The character of phase coexistence changes from a spherical droplet (or bubble) of the minority phase (near the coexistence curve) to a cylindrical droplet (or bubble) and finally (in the center of the miscibility gap) to a slab-like configuration of two parallel flat interfaces. Extending the analysis of M. Schrader, P. Virnau, and K. Binder [Phys. Rev. E 79, 061104 (2009)], we extract the surface free energy gamma (R) of both spherical and cylindrical droplets and bubbles in the vapor-liquid case, and present evidence that for R -> Infinity the leading order (Tolman) correction for droplets has sign opposite to the case of bubbles, consistent with the Tolman length being independent on the sign of curvature. For the symmetric binary mixture the expected non-existence of the Tolman length is confirmed. In all cases {and for a range of radii} R relevant for nucleation theory, gamma(R) deviates strongly from gamma (Infinity) which can be accounted for by a term of order gamma(Infinity)/gamma(R)-1 ~ 1/R^2. Our results for the simple Lennard-Jones fluid are also compared to results from density functional theory and we find qualitative agreement in the behavior of gamma(R) as well as in the sign and magnitude of the Tolman length.

preprint2010arXiv

Finite size scaling in Ising-like systems with quenched random fields: Evidence of hyperscaling violation

In systems belonging to the universality class of the random field Ising model, the standard hyperscaling relation between critical exponents does not hold, but is replaced by a modified hyperscaling relation. As a result, standard formulations of finite size scaling near critical points break down. In this work, the consequences of modified hyperscaling are analyzed in detail. The most striking outcome is that the free energy cost ΔF of interface formation at the critical point is no longer a universal constant, but instead increases as a power law with system size, ΔF proportional to $L^θ$ , with $θ$ the violation of hyperscaling critical exponent, and L the linear extension of the system. This modified behavior facilitates a number of new numerical approaches that can be used to locate critical points in random field systems from finite size simulation data. We test and confirm the new approaches on two random field systems in three dimensions, namely the random field Ising model, and the demixing transition in the Widom-Rowlinson fluid with quenched obstacles.

preprint2010arXiv

Monte Carlo simulations of the solid-liquid transition in hard spheres and colloid-polymer mixtures

Monte Carlo simulations at constant pressure are performed to study coexistence and interfacial properties of the liquid-solid transition in hard spheres and in colloid-polymer mixtures. The latter system is described as a one-component Asakura-Oosawa (AO) model where the polymer's degrees of freedom are incorporated via an attractive part in the effective potential for the colloid-colloid interactions. For the considered AO model, the polymer reservoir packing fraction is eta_p^r=0.1 and the colloid-polymer size ratio is q=sigma_p/σ=0.15 (with sigma_p and sigma the diameter of polymers and colloids, respectively). Inhomogeneous solid-liquid systems are prepared by placing the solid fcc phase in the middle of a rectangular simulation box creating two interfaces with the adjoined bulk liquid. By analyzing the growth of the crystalline region at various pressures and for different system sizes, the coexistence pressure p_co is obtained, yielding p_co=11.576 k_BT/sigma^3 for the hard sphere system and p_co=8.0 k_BT/sigma^3 for the AO model (with k_B the Boltzmann constant and T the temperature). Several order parameters are introduced to distinguish between solid and liquid phases and to describe the interfacial properties. From the capillary-wave broadening of the solid-liquid interface, the interfacial stiffness is obtained for the (100) crystalline plane, giving the values gamma=0.49 k_BT/sigma^2 for the hard-sphere system and gamma=0.95 k_BT/sigma^2 for the AO model.

preprint2009arXiv

Computer simulation studies of finite-size broadening of solid-liquid interfaces: From hard spheres to nickel

Using Molecular Dynamics (MD) and Monte Carlo (MC) simulations interfacial properties of crystal-fluid interfaces are investigated for the hard sphere system and the one-component metallic system Ni (the latter modeled by a potential of the embedded atom type). Different local order parameters are considered to obtain order parameter profiles for systems where the crystal phase is in coexistence with the fluid phase, separated by interfaces with (100) orientation of the crystal. From these profiles, the mean-squared interfacial width w^2 is extracted as a function of system size. We rationalize the prediction of capillary wave theory that w^2 diverges logarithmically with the lateral size of the system. We show that one can estimate the interfacial stiffness from the interfacial broadening, obtaining 0.5 k_B T/sigma^2 for hard spheres and 0.18 J/m^2 for Ni.

preprint2009arXiv

Hydrokinetic simulations of nanoscopic precursor films in rough channels

We report on simulations of capillary filling of high-wetting fluids in nano-channels with and without obstacles. We use atomistic (molecular dynamics) and hydrokinetic (lattice-Boltzmann) approaches which point out clear evidence of the formation of thin precursor films, moving ahead of the main capillary front. The dynamics of the precursor films is found to obey a square-root law as the main capillary front, z^2(t) ~ t, although with a larger prefactor, which we find to take the same value for the different geometries (2D-3D) under inspection. The two methods show a quantitative agreement which indicates that the formation and propagation of thin precursors can be handled at a mesoscopic/hydrokinetic level. This can be considered as a validation of the Lattice-Boltzmann (LB) method and opens the possibility of using hydrokinetic methods to explore space-time scales and complex geometries of direct experimental relevance. Then, LB approach is used to study the fluid behaviour in a nano-channel when the precursor film encounters a square obstacle. A complete parametric analysis is performed which suggests that thin-film precursors may have an important influence on the efficiency of nanochannel-coating strategies.

preprint2000arXiv

Defect fugacity, Spinwave Stiffness and T_c of the 2-d Planar Rotor Model

We obtain precise values for the fugacities of vortices in the 2-d planar rotor model from Monte Carlo simulations in the sector with {\em no} vortices. The bare spinwave stiffness is also calculated and shown to have significant anharmonicity. Using these as inputs in the KT recursion relations, we predict the temperature T_c = 0.925, using linearised equations, and $T_c = 0.899 \pm >.005$ using next higher order corrections, at which vortex unbinding commences in the unconstrained system. The latter value, being in excellent agreement with all recent determinations of T_c, demonstrates that our method 1) constitutes a stringent measure of the relevance of higher order terms in KT theory and 2) can be used to obtain transition temperatures in similar systems with modest computational effort.

preprint2000arXiv

Equation of state and critical behavior of polymer models: A quantitative comparison between Wertheim's thermodynamic perturbation theory and computer simulations

We present an application of Wertheim's Thermodynamic Perturbation Theory (TPT1) to a simple coarse grained model made of flexibly bonded Lennard-Jones monomers. We use both the Reference Hyper-Netted-Chain (RHNC) and Mean Spherical approximation (MSA) integral equation theories to describe the properties of the reference fluid. The equation of state, the density dependence of the excess chemical potential, and the critical points of the liquid--vapor transition are compared with simulation results and good agreement is found. The RHNC version is somewhat more accurate, while the MSA version has the advantage of being almost analytic. We analyze the scaling behavior of the critical point of chain fluids according to TPT1 and find it to reproduce the mean field exponents: The critical monomer density is predicted to vanish as $n^{-1/2}$ upon increasing the chain length $n$ while the critical temperature is predicted to reach an asymptotic finite temperature that is attained as $n^{-1/2}$ . The predicted asymptotic finite critical temperature obtained from the RHNC and MSA versions of TPT1 is found to be in good agreement with the $Θ$ point of our polymer model as obtained from the temperature dependence of the single chain conformations.

preprint1999arXiv

Cooling Rate Dependence and Dynamic Heterogeneity Below the Glass Transition in a Lennard-Jones Glass

We investigate a binary Lennard-Jones mixture with molecular dynamics simulations. We consider first a system cooled linearly in time with the cooling rate gamma. By varying gamma over almost four decades we study the influence of the cooling rate on the glass transition and on the resulting glass. We find for all investigated quantities a cooling rate dependence; with decreasing cooling rate the system falls out of equilibrium at decreasing temperatures, reaches lower enthalpies and obtains increasing local order. Next we study the dynamics of the melting process by investigating the most immobile and most mobile particles in the glass. We find that their spatial distribution is heterogeneous and that the immobile/mobile particles are surrounded by denser/less dense cages than an average particle.

preprint1998arXiv

Interfacial profiles between coexisting phases in thin films: Cahn Hilliard treatment versus capillary waves

We consider an interface between two demixed A and B phases, confined in a thin film between two antisymmetric walls, one of which prefers A and the other B. Above the wetting transition, the interface is stabilized in the center of the film. Based on a suitable extension of the Cahn-Hilliard gradient-square theory in combination with the capillary wave theory, we argue that the confinement influences the interfacial structure in two ways: It squeezes the intrinsic structure and cuts off the capillary wave spectrum. As a result, the interfacial width is proportional to the film thickness D in very thin films, and proportional to the square root of D in thicker films. These effects are then discussed in detail for the special case of an interface between demixed homopolymer phases. The width of the intrinsic profile is calculated analytically as a function of film thickness in the Cahn-Hilliard approximation (weak segregation limit) and in the Helfand theory (strong segregation limit), and numerically in the self-consistent field approximation. The results are confirmed by Monte Carlo simulations of a lattice model for a polymer mixture.

preprint1997arXiv

Diblock copolymers at a homopolymer-homopolymer-interface: a Monte Carlo simulation

The properties of diluted symmetric A-B diblock copolymers at the interface between A and B homopolymer phases are studied by means of Monte Carlo (MC) simulations of the bond fluctuation model. We calculate segment density profiles as well as orientational properties of segments, of A and B blocks, and of the whole chain. Our data support the picture of oriented ``dumbbells'', which consist of mildly perturbed A and B Gaussian coils. The results are compared to a self consistent field theory (SCFT) for single copolymer chains at a homopolymer interface. We also discuss the number of interaction contacts between monomers, which provide a measure for the ``active surface'' of copolymers or homopolymers close to the interface.

preprint1997arXiv

Interfaces in partly compatible polymer mixtures: A Monte Carlo simulation approach

The structure of polymer coils near interfaces between coexisting phases of symmetrical polymer mixtures (AB) is discussed, as well as the structure of symmetric diblock copolymers of the same chain length N adsorbed at the interface. The problem is studied by Monte Carlo simulations of the bond fluctuation model on the simple cubic lattice, using massively parallel computers (CRAY T3D). While homopolymer coils in the strong segregation limit are oriented parallel to the interface, the diblocks form ``dumbbells'' oriented perpendicular to the interface. However, in the dilute case (``mushroom regime'' rather than ``brush regime''), the diblocks are only weakly stretched. Distribution functions for monomers at the chain ends and in the center of the polymer are obtained, and a comparison to the self consistent field theory is made.

preprint1996arXiv

Finite-size scaling for near-critical continuum fluids at constant pressure

We consider the application of finite-size scaling methods to isothermal-isobaric (constant-NpT) simulations of pure continuum fluids. A finite-size scaling ansatz is made for the dependence of the relevant scaling operators on the particle number. To test the proposed scaling form, constant pressure simulations of the Lennard-Jones fluid at its liquid-vapour critical point are performed. The critical scaling operator distributions are obtained and their scaling with particle number found to be consistent with the proposed behaviour. The forms of these scaling distributions are shown to be identical to their Ising model counterparts. The relative merits of employing the constant-NpT and grand canonical (constant- $μ$ VT) ensembles for simulations of fluid criticality are also discussed.

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.

K. Binder

Research overview

Is this your work?

20 publications

Magnetic order and electromagnon excitations in DyMnO3 studied by neutron scattering experiments

Mixtures of anisotropic and spherical colloids: Phase behavior, confinement, percolation phenomena and kinetics

Star Polymers Confined in a Nanoslit: A Simulation Test of Scaling and Self-Consistent Field Theories

The escape transition of a compressed star polymer: Self-consistent field predictions tested by simulation

Langevin Dynamics simulations of a 2-dimensional colloidal crystal under confinement and shear

Microcanonical Determination of the Interface Tension of Flat and Curved Interfaces from Monte Carlo Simulations

A slow process in confined polymer melts: layer exchange dynamics at a polymer solid interface

Absorption/Expulsion of Oligomers and Linear Macromolecules in a Polymer Brush

Curvature Dependence of Surface Free Energy of Liquid Drops and Bubbles: A Simulation Study

Finite size scaling in Ising-like systems with quenched random fields: Evidence of hyperscaling violation

Monte Carlo simulations of the solid-liquid transition in hard spheres and colloid-polymer mixtures

Computer simulation studies of finite-size broadening of solid-liquid interfaces: From hard spheres to nickel

Hydrokinetic simulations of nanoscopic precursor films in rough channels

Defect fugacity, Spinwave Stiffness and T_c of the 2-d Planar Rotor Model

Equation of state and critical behavior of polymer models: A quantitative comparison between Wertheim's thermodynamic perturbation theory and computer simulations

Cooling Rate Dependence and Dynamic Heterogeneity Below the Glass Transition in a Lennard-Jones Glass

Interfacial profiles between coexisting phases in thin films: Cahn Hilliard treatment versus capillary waves

Diblock copolymers at a homopolymer-homopolymer-interface: a Monte Carlo simulation

Interfaces in partly compatible polymer mixtures: A Monte Carlo simulation approach

Finite-size scaling for near-critical continuum fluids at constant pressure

See the researcher in context

Building this map preview