Search

We report simulations on the homogeneous liquid-fcc nucleation of charged colloids for both low and high contact energy values. As a precursor for crystal formation, we observe increased local order at the position where the crystal will form, but no correlations with the local density. Thus, the nucleation is driven by order fluctuations rather than density fluctuations. Our results also show that the transition involves two stages in both cases, first a transition liquid-bcc, followed by a bcc-hcp/fcc transition. Both transitions have to overcome free energy barriers, so that a spherical bcc-like cluster is formed first, in which the final fcc-like structure is nucleated mainly at the surface of the crystallite. This means that the bcc-fcc phase transition is a heterogeneous nucleation, even though we start from a homogeneous bulk liquid. The height of the bcc-hcp/fcc free energy barrier strongly depends on the contact energies of the colloids. For low contact energy this barrier is low, so that the bcc-hcp/fcc transition happens spontaneously. For the higher contact energy, the second barrier is too high to be crossed spontaneously by the colloidal system. However, it was possible to ratchet the system over the second barrier and to transform the bcc nuclei into the stable hcp/fcc phase. The transitions are dominated by the first liquid-bcc transition and can be described by Classical Nucleation Theory using an effective surface tension., Comment: 9 pages, 6 figures

Nano-particle agglomeration plays an important role in processes such as spray drying and particle flame synthesis. These processes have in common that nano-particles collide at low concentrations and get irreversibly linked at the point of contact due to plastic deformation. In this paper, we investigate several models of irreversible connections, which require only local interactions between the colliding nano-particles and thus allow for scalable simulations. The models investigated here connect the particles upon collision by non-bonded strongly attractive interactions, bonded interactions or by binding agents placed at the point of contact. Models using spherically symmetric interactions form compact agglomerates and are therefore unsuitable to study agglomeration. In contrast, models that are either based on both central and angular potentials (type one) or on binding agents (type two) efficiently prevent restructuring of the agglomerates, and are therefore useful for modeling contacts formed by plastic deformation. Moreover, both types of models allow to control the rigidity and by that the degree of restructuring. The first type of model is computationally more efficient at low fractional dimensions of the aggregates, while the second gives easy access to local shear forces, which is important when breaking of agglomerates is to be considered. As example applications, we reproduce the well-known diffusion-limited agglomeration (DLA) and report results on soot aggregation., Comment: Updated due to missing references

Understanding dielectric spectra can reveal important information about the dynamics of solvents and solutes from the dipolar relaxation times down to electronic ones. In the late 1970s, Hubbard and Onsager predicted that adding salt ions to a polar solution would result in a reduced dielectric permittivity that arises from the unexpected tendency of solvent dipoles to align opposite to the applied field. So far, this effect has escaped an experimental verification, mainly because of the concomitant appearance of dielectric saturation from which the Hubbard-Onsager decrement cannot be easily separated. Here we develop a novel non-equilibrium molecular dynamics simulation approach to determine this decrement accurately for the first time. Using a thermodynamic consistent all-atom force field we show that for an aqueous solution containing sodium chloride around 4.8 Mol/l, this effect accounts for 12\% of the total dielectric permittivity. The dielectric decrement can be strikingly different if a less accurate force field for the ions is used. Using the widespread GROMOS parameters, we observe in fact an {\it increment} of the dielectric permittivity rather than a decrement. We can show that this increment is caused by ion pairing, introduced by a too low dispersion force, and clarify the microscopic connection between long-living ion pairs and the appearance of specific features in the dielectric spectrum of the solution.

With the help of a recently developed non-equilibrium approach, we investigate the ionic strength dependence of the Hubbard--Onsager dielectric decrement. We compute the depolarization of water molecules caused by the motion of ions in sodium chloride solutions from the dilute regime (0.035 M) up close to the saturation concentration (4.24 M), and find that the kinetic decrement displays a strong nonmonotonic behavior, in contrast to the prediction of available models. We introduce a phenomenological modification of the Hubbard--Onsager continuum theory, that takes into account the screening due to the ionic cloud at mean field level, and, is able to describe the kinetic decrement at high concentrations including the presence of a pronounced minimum.

Crystallization is often assumed to be a quasi-static process that is unaffected by details of particle transport other than the bulk diffusion coefficient. Therefore colloidal suspensions are frequently argued to be an ideal toy model for experimentally more difficult systems such as metal melts. In this letter, we want to challenge this assumption. To this aim, we have considered molecular dynamics simulations of the crystallization in a suspension of Yukawa-type colloids. In order to investigate the role of hydrodynamic interactions (HIs) mediated by the solvent, we modeled the solvent both implicitly and explicitly, using Langevin dynamics and the fluctuating Lattice Boltzmann method, respectively. Our simulations show a dramatic reduction of the crystal growth velocity due to HIs even at moderate hydrodynamic coupling. A detailed analysis shows that this slowdown is due to the wall-like properties of the crystal surface, which reduces the colloidal diffusion towards the crystal surface by hydrodynamic screening. Crystallization in suspensions therefore differs strongly from pure melts, making them less useful as a toy model than previously thought.

Forward flux sampling (FFS) provides a convenient and efficient way to simulate rare events in equilibrium or non-equilibrium systems. FFS ratchets the system from an initial state to a final state via a series of interfaces in phase space. The efficiency of FFS depends sensitively on the positions of the interfaces. We present two alternative methods for placing interfaces automatically and adaptively in their optimal locations, on-the-fly as an FFS simulation progresses, without prior knowledge or user intervention. These methods allow the FFS simulation to advance efficiently through bottlenecks in phase space by placing more interfaces where the probability of advancement is lower. The methods are demonstrated both for a single-particle test problem and for the crystallization of Yukawa particles. By removing the need for manual interface placement, our methods both facilitate the setting up of FFS simulations and improve their performance, especially for rare events which involve complex trajectories through phase space, with many bottlenecks., Comment: 13 pages, 9 figures

We report molecular dynamics simulations of the segregation of two overlapping chains in cylindrical confinement. We find that the entropic repulsion between the chains can be sufficiently strong to cause segregation on a time scale that is short compared to the one for diffusion. This result implies that entropic driving forces are sufficiently strong to cause rapid bacterial chromosome segregation., Comment: Minor changes. Added some references, corrected the labels in figure 6 and reformatted in two columns. Also added reference to published version in PRE

We study the link between three seeming-disparate cases of self-avoiding polymers: strongly overlapping multiple chains in dilute solution, chains under spherical confinement, and the onset of semi-dilute solutions. Our main result is that the free energy for overlapping $n$ chains is \emph{independent} of chain length and scales as $n^{9/4}$, slowly crossing over to $n^3$, as $n$ increases. For strongly confined polymers inside a spherical cavity, we show that rearranging the chains does not cost an additional free energy. Our results imply that, during cell cycle, \emph{global} reorganization of eukaryotic chromosomes in a large cell nucleus could be readily achieved., Comment: to appear in Physical Review Letters

We investigate the effective interaction between two like-charged rods in the regime of large electrostatic coupling parameters using both Molecular Dynamics simulation techniques and the recently introduced strong-coupling theory. We obtain attraction between the rods for elevated Manning parameters accompanied by an equilibrium surface-to-surface separation of the order of the Gouy-Chapman length. A continuous unbinding between the rods is predicted at a threshold Manning parameter of 2/3., Comment: 8 pages, 2 figures

We use Molecular Dynamics simulations to study attractive interactions and the underlying ionic correlations between parallel like-charged rods in the absence of additional salt. For a generic bulk system of rods we identify a reduction of short range repulsions as the origin of a negative osmotic coefficient. The counterions show signs of a weak three-dimensional order in the attractive regime only once the rod-imposed charge-inhomogeneities are divided out. We also treat the case of attraction between a single pair of rods for a few selected line charge densities and rod radii. Measurements of the individual contributions to the force between close rods are studied as a function of Bjerrum length. We find that even though the total force is always attractive at sufficiently high Bjerrum length, the electrostatic contribution can ultimately become repulsive. We also measure azimuthal and longitudinal correlation functions to answer the question how condensed ions are distributed with respect to each other and to the neighboring rod. For instance, we show that the prevalent image of mutually interlocked ions is qualitatively correct, even though modifications due to thermal fluctuations are usually strong., Comment: 14 pages, 14 figures, REVTeX4 style

In a previous paper a method was developed to subtract the interactions due to periodically replicated charges (or other long-range entities) in one spatial dimension. The method constitutes a generalized "electrostatic layer correction" (ELC) which adapts any standard 3D summation method to slab-like conditions. Here the implementation of the layer correction is considered in detail for the standard Ewald (EW3DLC) and the PPPM mesh Ewald (PPPMLC) methods. In particular this method offers a strong control on the accuracy and an improved computational complexity of O(N log N) for mesh-based implementations. We derive anisotropic Ewald error formulas and give some fundamental guidelines for optimization. A demonstration of the accuracy, error formulas and computation times for typical systems is also presented., Comment: 14 pages, 7 figures

We propose a new method to sum up electrostatic interactions in 2D slab geometries. It consists of a combination of two recently proposed methods, the 3D Ewald variant of Yeh and Berkowitz, J. Chem. Phys. 111 (1999) 3155, and the purely 2D method MMM2D by Arnold and Holm, to appear in Chem. Phys. Lett. 2002. The basic idea involves two steps. First we use a three dimensional summation method whose summation order is changed to sum up the interactions in a slab-wise fashion. Second we subtract the unwanted interactions with the replicated layers analytically. The resulting method has full control over the introduced errors. The time to evaluate the layer correction term scales linearly with the number of charges, so that the full method scales like an ordinary 3D Ewald method, with an almost linear scaling in a mesh based implementation. In this paper we will introduce the basic ideas, derive the layer correction term and numerically verify our analytical results., Comment: 10 pages, 7 figures

We present a new method, in the following called MMM2D, to accurately calculate the electrostatic energy and forces on charges being distributed in a two dimensional periodic array of finite thickness. It is not based on an Ewald summation method and as such does not require any fine-tuning of an Ewald parameter for convergence. We transform the Coulomb sum via a convergence factor into a series of fast decaying functions which can be easily evaluated. Rigorous error bounds for the energies and the forces are derived and numerically verified. Already for small systems our method is much faster than the traditional 2D-Ewald methods, but for large systems it is clearly superior because its time demand scales like O(N^{5/3}) with the number N of charges considered. Moreover it shows a rapid convergence, is very precise and easy to handle., Comment: 29 pages, 6 figures, needs elsart.cls (provided)