Search

We describe an approach to online analysis of micro-rheology data using a multi-scale time-correlation method. The method is particularly suited to process high-volume data streams and compress the relevant information in real time. Using this, we can obtain complex moduli of visco-elastic media without suffering from the high-frequency artefacts that are associated with the truncation errors in the most widely used versions of micro-rheology. Moreover, the present approach obviates the need to choose the time interval for data acquisition beforehand. We test our approach first on an artificial data set and then on experimental data obtained both for an optically trapped colloidal probe in water and a similar probe in poly-ethylene glycol solutions at various concentrations. In all cases, we obtain good agreement with the bulk rheology data in the region of overlap. We compare our method with the conventional Kramers-Kronig transform approach and find that the two methods agree over most of the frequency regime. For the same data set, the present approach is superior to Kramers-Kronig at high frequencies and can be made to perform at least comparable at low frequencies., Comment: 16 pages, 9 figures

We report a Molecular Dynamics study of homogenous bubble nucleation in a Lennard-Jones fluid. The rate of bubble nucleation is estimated using forward-flux sampling (FFS). We find that cavitation starts with compact bubbles rather than with ramified structures as had been suggested by Shen and Debenedetti (J. Chem. Phys. 111:3581, 1999). Our estimate of the bubble-nucleation rate is higher than predicted on the basis of Classical Nucleation Theory (CNT). Our simulations show that local temperature fluctuations correlate strongly with subsequent bubble formation - this mechanism is not taken into account in CNT.

We study the kinetics of the liquid-to-crystal transformation and of gel formation in colloidal suspensions of oppositely charged particles. We analyse, by means of both computer simulations and experiments, the evolution of a fluid quenched to a state point of the phase diagram where the most stable state is either a homogeneous crystalline solid or a solid phase in contact with a dilute gas. On the one hand, at high temperatures and high packing fractions, close to an ordered-solid/disordered-solid coexistence line, we find that the fluid-to-crystal pathway does not follow the minimum free energy route. On the other hand, a quench to a state point far from the ordered-crystal/disordered-crystal coexistence border is followed by a fluid-to-solid transition through the minimum free energy pathway. At low temperatures and packing fractions we observe that the system undergoes a gas-liquid spinodal decomposition that, at some point, arrests giving rise to a gel-like structure. Both our simulations and experiments suggest that increasing the interaction range favors crystallization over vitrification in gel-like structures.

We present a method for computing stationary distributions for activated processes in equilibrium and non-equilibrium systems using Forward Flux Sampling (FFS). In this method, the stationary distributions are obtained directly from the rate constant calculations for the forward and backward reactions; there is no need to perform separate calculations for the stationary distribution and the rate constant. We apply the method to the non-equilibrium rare event problem proposed by Maier and Stein, to nucleation in a 2-dimensional Ising system, and to the flipping of a genetic switch.

We report a computer-simulation study of the equilibrium phase diagram of a three-dimensional system of particles with a repulsive step potential. The phase diagram is obtained using free-energy calculations. At low temperatures, we observe a number of distinct crystal phases. We show that at certain values of the potential parameters the system exhibits the water-like thermodynamic anomalies: density anomaly and diffusion anomaly. The anomalies disappear with increasing the repulsive step width: their locations move to the region inside the crystalline phase., Comment: 6 pages, 5 figures

We report the phase diagram of interpenetrating Hertzian spheres. The Hertz potential is purely repulsive, bounded at zero separation, and decreases monotonically as a power law with exponent 5/2, vanishing at the overlapping threshold. This simple functional describes the elastic interaction of weakly deformable bodies and, therefore, it is a reliable physical model of soft macromolecules, like star polymers and globular micelles. Using thermodynamic integration and extensive Monte Carlo simulations, we computed accurate free energies of the fluid phase and a large number of crystal structures. For this, we defined a general primitive unit cell that allows for the simulation of any lattice. We found multiple re-entrant melting and first-order transitions between crystals with cubic, trigonal, tetragonal, and hexagonal symmetries., Comment: The inset in Fig. 4 in the previous version had incorrect ordinate values. This has been corrected here. This mistake does not affect the remaining content of the article

We compute rates and pathways for nucleation in a sheared two dimensional Ising model with Metropolis spin flip dynamics, using Forward Flux Sampling (FFS). We find a peak in the nucleation rate at intermediate shear rate. We analyse the origin of this peak using modified shear algorithms and committor analysis. We find that the peak arises from an interplay between three shear-mediated effects: shear-enhanced cluster growth, cluster coalescence and cluster breakup. Our results show that complex nucleation behaviour can be found even in a simple driven model system. This work also demonstrates the use of FFS for simulating rare events, including nucleation, in nonequilibrium systems., Comment: 11 pages, 13 figures

We conduct a numerical study of the dynamical behavior of a system of three-dimensional crosses, particles that consist of three mutually perpendicular line segments rigidly joined at their midpoints. In an earlier study [W. van Ketel et al., Phys. Rev. Lett. 94, 135703 (2005)] we showed that this model has the structural properties of an ideal gas, yet the dynamical properties of a strong glass former. In the present paper we report an extensive study of the dynamical heterogeneities that appear in this system in the regime where glassy behavior sets in. On the one hand, we find that the propensity of a particle to diffuse is determined by the structure of its local environment. The local density around mobile particles is significantly less than the average density, but there is little clustering of mobile particles, and the clusters observed tend to be small. On the other hand, dynamical susceptibility results indicate that a large dynamical length scale develops even at moderate densities. This suggests that propensity and other mobility measures are an incomplete measure of dynamical length scales in this system., Comment: 11 pages, 7 figures

We report a computer-simulation study of the equilibrium phase diagram of a three-dimensional system of particles with a repulsive step potential. Using free-energy calculations, we have determined the equilibrium phase diagram of this system. At low temperatures, we observe a number of distinct crystal phases. However, under certain conditions the system undergoes a glass transition in a regime where the liquid appears thermodynamically stable. We argue that the appearance of this amorphous low-temperature phase can be understood by viewing this one-component system as a pseudo-binary mixture., Comment: 4 pages, 4 figures

We report a study of the phase behavior of multiple-occupancy crystals through simulation. We argue that in order to reproduce the equilibrium behavior of such crystals it is essential to treat the number of lattice sites as a constraining thermodynamic variable. The resulting free-energy calculations thus differ considerably from schemes used for single-occupancy lattices. Using our approach, we obtain the phase diagram and the bulk modulus for a generalized exponential model that forms cluster crystals at high densities. We compare the simulation results with existing theoretical predictions. We also identify two types of density fluctuations that can lead to two sound modes and evaluate the corresponding elastic constants., Comment: 4 pages, 3 figures

We report a numerical study of the rate of crystal nucleation in a binary suspension of oppositely charged colloids. Two different crystal structures compete in the thermodynamic conditions under study. We find that the crystal phase that nucleates is metastable and, more surprisingly, its nucleation free energy barrier is not the lowest one. This implies that, during nucleation, there is insufficient time for sub-critical nuclei to relax to their lowest free-energy structure. Such behavior is in direct contradiction with the common assumption that the phase that crystallizes most readily is the one with the lowest free-energy barrier for nucleation. The phenomenon that we describe should be relevant for crystallization experiments where competing solid structures are not connected by an easy transformation., Comment: 4 pages 3 figures. Accepted Phys. Rev. Lett

The comprehensive computer simulation study of the phase diagram of the repulsive step potential system in three dimensions is represented. We show that the system with a simple purely repulsive isotropic potential demonstrates a number of unusual features. The maxima and minima on the melting curve are found for some regions of potential parameters. It is shown that the phase diagram in $\rho-T$ plane includes two isostructural crystalline parts separated by the disordered phase which is amorphous at low enough temperatures. Phase diagram in the (P-T) plane shows that the transition to the amorphous state occurs approximately along the extrapolated spinodals. Structural FCC-BCC phase transition is found at high densities., Comment: 4 pages, 5 figures

We perform Monte Carlo simulations on the hard-core attractive Yukawa system to test the Optimized Baxter Model that was introduced in [P.Prinsen and T. Odijk, J. Chem. Phys. 121, p.6525 (2004)] to study a fluid phase of spherical particles interacting through a short-range pair potential. We compare the chemical potentials and pressures from the simulations with analytical predictions from the Optimized Baxter Model. We show that the model is accurate to within 10 percent over a range of volume fractions from 0.1 to 0.4, interaction strengths up to three times the thermal energy and interaction ranges from 6 to 20 % of the particle diameter, and performs even better in most cases. We furthermore establish the consistency of the model by showing that the thermodynamic properties of the Yukawa fluid computed via simulations may be understood on the basis of one similarity variable, the stickiness parameter defined within the Optimized Baxter Model. Finally we show that the Optimized Baxter Model works significantly better than an often used, naive method determining the stickiness parameter by equating the respective second virial coefficients based on the attractive Yukawa and Baxter potentials., Comment: 11 pages, 8 figures

We analyse the efficiency of several simulation methods which we have recently proposed for calculating rate constants for rare events in stochastic dynamical systems, in or out of equilibrium. We derive analytical expressions for the computational cost of using these methods, and for the statistical error in the final estimate of the rate constant, for a given computational cost. These expressions can be used to determine which method to use for a given problem, to optimize the choice of parameters, and to evaluate the significance of the results obtained. We apply the expressions to the two-dimensional non-equilibrium rare event problem proposed by Maier and Stein. For this problem, our analysis gives accurate quantitative predictions for the computational efficiency of the three methods., Comment: 19 pages, 13 figures

We present three algorithms for calculating rate constants and sampling transition paths for rare events in simulations with stochastic dynamics. The methods do not require a priori knowledge of the phase space density and are suitable for equilibrium or non-equilibrium systems in stationary state. All the methods use a series of interfaces in phase space, between the initial and final states, to generate transition paths as chains of connected partial paths, in a ratchet-like manner. No assumptions are made about the distribution of paths at the interfaces. The three methods differ in the way that the transition path ensemble is generated. We apply the algorithms to kinetic Monte Carlo simulations of a genetic switch and to Langevin Dynamics simulations of intermittently driven polymer translocation through a pore. We find that the three methods are all of comparable efficiency, and that all the methods are much more efficient than brute force simulation., Comment: 18 pages, 11 figures

We report a series of Lattice-Boltzmann simulations of the sedimentation velocity of charged disks. In these simulations, we explicitly account for the hydrodynamic and electrostatic forces on disks and on their electrical double layer. By comparing our results with those for spheres with equal surface and charge, we can clarify the effect of the particle shape on the sedimentation process. We find that disks and spheres exhibit a different dependence of the sedimentation velocity on the Debye screening length. An analysis of the behavior of highly charged disks (beyond the scope of the linearized Poisson-Boltzmann equation) shows that, in that regime, the charge dependence of the sedimentation velocity of disks and spheres is similar. This suggests that, at high charge, the effective hydrodynamic shape of the disks becomes more spherical.

The evolutionary reason for the increase in gene length from archaea to prokaryotes to eukaryotes observed in large scale genome sequencing efforts has been unclear. We propose here that the increasing complexity of protein-protein interactions has driven the selection of longer proteins, as longer proteins are more able to distinguish among a larger number of distinct interactions due to their greater average surface area. Annotated protein sequences available from the SWISS-PROT database were analyzed for thirteen eukaryotes, eight bacteria, and two archaea species. The number of subcellular locations to which each protein is associated is used as a measure of the number of interactions to which a protein participates. Two databases of yeast protein-protein interactions were used as another measure of the number of interactions to which each \emph{S. cerevisiae} protein participates. Protein length is shown to correlate with both number of subcellular locations to which a protein is associated and number of interactions as measured by yeast two-hybrid experiments. Protein length is also shown to correlate with the probability that the protein is encoded by an essential gene. Interestingly, average protein length and number of subcellular locations are not significantly different between all human proteins and protein targets of known, marketed drugs. Increased protein length appears to be a significant mechanism by which the increasing complexity of protein-protein interaction networks is accommodated within the natural evolution of species. Consideration of protein length may be a valuable tool in drug design, one that predicts different strategies for inhibiting interactions in aberrant and normal pathways., Comment: 13 pages, 5 figures, 2 tables, to appear in Physica A

The turn-over frequency of the catalytic oxidation of CO at RuO2(110) was calculated as function of temperature and partial pressures using ab initio statistical mechanics. The underlying energetics of the gas-phase molecules, dissociation, adsorption, surface diffusion, surface chemical reactions, and desorption were obtained by all-electron density-functional theory. The resulting CO2 formation rate [in the full (T, p_CO, p_O2)-space], the movies displaying the atomic motion and reactions over times scales from picoseconds to seconds, and the statistical analyses provide insights into the concerted actions ruling heterogeneous catalysis and open thermodynamic systems in general., Comment: 4 pages including 3 figures, Related publications can be found at http://www.fhi-berlin.mpg.de/th/paper.html

The structure of the Baxter adhesive hard sphere fluid is examined using computer simulation. The radial distribution function (which exhibits unusual discontinuities due to the particle adhesion) and static structure factor are calculated with high accuracy over a range of conditions and compared with the predictions of Percus--Yevick theory. We comment on rigidity in percolating clusters and discuss the role of the model in the context of experiments on colloidal systems with short-range attractive forces., Comment: 14 pages, 7 figures. (For proceedings of "Structural arrest in colloidal systems with short-range attractive forces", Messina, December 2003)

We study the effect of shear flow on homogeneous crystal nucleation, using Brownian Dynamics simulations in combination with an umbrella sampling like technique. The symmetry breaking due to shear results in anisotropic radial distribution functions. The homogeneous shear rate suppresses crystal nucleation and leads to an increase of the size of the critical nucleus. These observations can be described by a simple, phenomenological extension of classical nucleation theory. In addition, we find that nuclei have a preferential orientation with respect to the direction of shear. On average the longest dimension of a nucleus is along the vorticity direction, while the shortest dimension is preferably perpendicular to that and slightly tilted with respect to the gradient direction., Comment: 10 pages, 8 figures, Submitted to J. Phys.: Condens. Matter

We use Brownian Dynamics simulations in combination with the umbrella sampling technique to study the effect of shear flow on homogeneous crystal nucleation. We find that a homogeneous shear rate leads to a significant suppression of the crystal nucleation rate and to an increase of the size of the critical nucleus. A simple, phenomenological extension of classical nucleation theory accounts for these observations. The orientation of the crystal nucleus is tilted with respect to the shear direction., Comment: 4 pages, 3 figures, Submitted to Phys. Rev. Lett

We present a robust scheme for solving the electrokinetic equations. This goal is achieved by combining the lattice-Boltzmann method (LB) with a discrete solution of the convection-diffusion equation for the different charged and neutral species that compose the fluid. The method is based on identifying the elementary fluxes between nodes, which ensures the absence of spurious fluxes in equilibrium. We show how the model is suitable to study electro-osmotic flows. As an illustration, we show that, by introducing appropriate dynamic rules in the presence of solid interfaces, we can compute the sedimentation velocity (and hence the sedimentation potential) of a charged sphere. Our approach does not assume linearization of the Poisson-Boltzmann equation and allows us for a wide variation of the Peclet number., Comment: 24 pages, 7 figures

We propose that the first stage of membrane fusion need not be the formation of a single stalk. Instead, we consider a scenario for defect-induced membrane fusion that proceeds cooperatively via multiple stalk formation. The defects (stalks or pores) attract each other via membrane-mediated capillary interactions that result in a condensation transition of the defects. The resulting dense phase of stalks corresponds to the so-called fusion intermediate.

We report density-functional based molecular dynamics simulations, that show that, with increasing pressure, liquid carbon undergoes a gradual transformation from a liquid with local three-fold coordination to a 'diamond-like' liquid. We demonstrate that this unusual structural change is well reproduced by an empirical bond order potential with isotropic long range interactions, supplemented by torsional terms. In contrast, state-of-the-art short-range bond-order potentials do not reproduce this diamond structure. This suggests that a correct description of long-range interactions is crucial for a unified description of the solid and liquid phases of carbon., Comment: 4 pages, 5 figures

The melting of a polydisperse hard disk system is investigated by Monte Carlo simulations in the semigrand canonical ensemble. This is done in the context of possible continuous melting by a dislocation unbinding mechanism, as an extension of the 2D hard disk melting problem. We find that while there is pronounced fractionation in polydispersity, the apparent density-polydispersity gap does not increase in width, contrary to 3D polydisperse hard spheres. The point where the Young's modulus is low enough for the dislocation unbinding to occur moves with the apparent melting point, but stays within the density gap, just like for the monodisperse hard disk system. Additionally, we find that throughout the accessible polydispersity range, the bound dislocation-pair concentration is high enough to affect the dislocation unbinding melting as predicted by Kosterlitz, Thouless, Halperin, Nelson and Young., Comment: 6 pages, 6 figures

We propose a model that can account for the experimentally observed phase behavior of DNA-nanoparticle assemblies (R. Jin et al., JACS 125, 1643 (2003); T. A. Taton et al., Science 289, 1757 (2000)). The binding of DNA-coated nano-particles by dissolved DNA linker can be described by exploiting an analogy with quantum particles obeying fractional statistics. In accordance with experimental findings, we predict that the phase-separation temperature of the nano-colloids increases with the DNA coverage of the colloidal surface. Upon the addition of salt, the demixing temperature increases logarithmically with the salt concentration. Our analysis suggests an experimental strategy to map microscopic DNA sequences onto the macroscopic phase behavior of the DNA-nanoparticle solutions. Such an approach should enhance the efficiency of methods to detect (single) mutations in specific DNA sequences.

We compute the equilibrium concentration of stacking faults and point defects in polydisperse hard-sphere crystals. We find that, while the concentration of stacking faults remains similar to that of monodisperse hard sphere crystals, the concentration of vacancies decreases by about a factor two. Most strikingly, the concentration of interstitials in the maximally polydisperse crystal may be some six orders of magnitude larger than in a monodisperse crystal. We show that this dramatic increase in interstitial concentration is due to the increased probability of finding small particles and that the small-particle tail of the particle size distribution is crucial for the interstitial concentration in a colloidal crystal., Comment: 6 pages, 4 figures

We investigate the phase behavior of tetrapods, hard non-convex bodies formed by 4 rods connected under tetrahedral angles. We predict that, depending on the relative lengths of the rods these particles can form a uniaxial nematic phase, and more surprisingly they can exhibit a cubatic phase, a special case of the biaxial nematic phase. These predictions may be experimentally testable, as experimental realizations of tetrapods have recently become available., Comment: 8 pages ReVTeX 4, including 3 EPS figures

We report a numerical calculation of the elastic constants of the fcc and hcp crystal phases of monodisperse hard-sphere colloids. Surprisingly, some of these elastic constants are very different (up to 20%), even though the free energy, pressure and bulk compressibility of the two crystal structures are very nearly equal. As a consequence, a moderate deformation of a hard-sphere crystal may make the hcp phase more stable than the fcc phase. This finding has implications for the design of patterned templates to grow colloidal hcp crystals. We also find that, below close packing, there is a small, but significant, difference between the distances between hexagonal layers (c/a ratios) of fcc and hcp crystals., Comment: 4 pages, 4 figures, accepted for publication in Physical Review Letters

Using a combination of Monte Carlo techniques, we locate the liquid--vapor critical point of adhesive hard spheres. We find that the critical point lies deep inside the gel region of the phase diagram. The (reduced) critical temperature and density are $\tau_c=0.1133\pm0.0005$ and $\rho_c=0.508\pm0.01$. We compare these results with the available theoretical predictions. Using a finite-size scaling analysis, we verify that the critical behavior of the adhesive hard sphere model is consistent with that of the 3D Ising universality class., Comment: 4 pages, 3 figures

In this paper, we report dynamic Monte Carlo simulations of melting and crystallization in a single-chain system. Their free energy barriers are calculated by the umbrella sampling method and can be described well by a simple expression DeltaF = n Deltaf+sigma (N-n)^(2/3), where n is the amount of molten bonds, Deltaf is the free energy change of each molten bond from a crystalline state, N is the chain length, and sigma is the surface free energy of crystallite. We found that, together with the expression Delta F = n Delta f+ sigma (N-n)^(1/2) for molecular nucleation, the molecular-weight dependent properties of the free-energy barriers for polymer primary and secondary nucleation, in particular, the molecular segregation during crystal growth, can be reproduced. Then for the mechanism of polymer crystallization, we suggested a quantitative model of intramolecular nucleation, as a direct development from the previous qualitative description of molecular nucleation model., Comment: 7 pages 11 figures

We report Monte Carlo simulations of a lattice-polymer model that can account for both polymer crystallization and liquid-liquid demixing in solutions of semiflexible homopolymers. In our model, neighboring polymer segments can have isotropic interactions that affect demixing, and anisotropic interactions that are responsible for freezing. However, our simulations show that the isotropic interactions also have a noticeable effect on the freezing curve, as do the anisotropic interactions on demixing. As the relative strength of the isotropic interactions is reduced, the liquid-liquid demixing transition disappears below the freezing curve. A simple, extended Flory-Huggins theory accounts quite well for the phase behavior observed in the simulations., Comment: Revtex, 7 pages, the content accepted by J. Chem. Phys

We report Monte Carlo simulations of the melting of a single-polymer crystallite. We find that, unlike most atomic and molecular crystals, such crystallites can be heated appreciably above their melting temperature before they transform to the disordered "coil" state. The surface of the superheated crystallite is found to be disordered. The thickness of the disordered layer increases with superheating. However, the order-disorder transition is not gradual but sudden. Free-energy calculations reveal the presence of a large free-energy barrier to melting., Comment: AMS-Latex, 4 pages with 5 figures, submitted to Phys. Rev. Lett

We present a new simulation scheme based on the Lattice-Boltzmann method to simulate the dynamics of charged colloids in an electrolyte. In our model we describe the electrostatics on the level of a Poisson-Boltzmann equation and the hydrodynamics of the fluid by the linearized Navier-Stokes equations. We verify our simulation scheme by means of a Chapman-Enskog expansion. Our method is applied to the calculation of the reduced sedimentation velocity U/U_0 for a cubic array of charged spheres in an electrolyte. We show that we recover the analytical solution first derived by Booth (F. Booth, J. Chem. Phys. 22, 1956 (1954)) for a weakly charged, isolated sphere in an unbounded electrolyte. The present method makes it possible to go beyond the Booth theory, and we discuss the dependence of the sedimentation velocity on the charge of the spheres. Finally we compare our results to experimental data., Comment: 18 pages, 5 figures, to appear in Phys. Rev. E

We report numerical calculations of the concentration of interstitials in hard-sphere crystals. We find that, in a three-dimensional fcc hard-sphere crystal at the melting point, the concentration of interstitials is 2 * 10^-8. This is some three orders of magnitude lower than the concentration of vacancies. A simple, analytical estimate yields a value that is in fair agreement with the numerical results., Comment: 12 pages, 2 figures; Submitted to J. Phys. Chem. B

We report a numerical simulation of the phase diagram of a simple model for membrane proteins constrained to move in a plane. In analogy with the corresponding three dimensional models, the liquid-gas transition becomes metastable as the range of attraction decreases. Spontaneous crystallization happens much more readily in the two dimensional models rather than in their three dimensional counterparts., Comment: Revtex, 9 pages, 7 figures

We propose an extension of the law of corresponding states that can be applied to systems - such as colloidal suspensions - that have widely different ranges of attractive interactions. We argue that, for such systems, the ``reduced'' second virial coefficient is a convenient parameter to quantify the effective range of attraction. This procedure allows us to give a simple definition of the effective range of attraction of potentials with different functional forms. The advantage of the present approach is that it allows us to estimate the relative location of the liquid-vapor and solid-fluid coexistence curves exclusively on basis of the knowledge of the pair-potential., Comment: REVTeX, 5 pages, 2 figures

We investigate the possibility that freely rotating cylinders with an aspect ratio $L/D=0.9$ exhibit a cubatic phase similar to the one found for a system of cut-spheres. We present theoretical arguments why a cubatic phase might occur in this particular system. Monte Carlo simulations do not confirm the existence of a cubatic phase for cylinders. However, they do reveal an unexpected phase behavior between the isotropic and crystalline phase., Comment: 24 pages, 12 figures, RevTex (Submitted to J. Chem. Phys.)

Temperature changes during the growth of lamellar polymer crystals give rise to steps on the surface of the crystals. It has recently been suggested that these steps could provide important insights into the mechanism of polymer crystallization. In particular, a characterization of the profiles of these steps might reveal the fixed-point attractor that underlies a recently proposed crystallization mechanism. Here we examine this hypothesis by performing simulations of such temperature jumps using the Sadler-Gilmer model. We find that for this model the step profiles do reveal the fixed-point attractor. However, for temperature decreases they also reflect the rounding of the crystal edge that occurs in this model and for temperature increases they also reflect the fluctuations in the thickness present in the crystal. We discuss the implications of these results for the interpretation of experimental step profiles., Comment: 8 pages, 7 figures, revtex

Recent work on the mechanism of polymer crystallization has led to a proposal for the mechanism of thickness selection which differs from those proposed by the surface nucleation theory of Lauritzen and Hoffman and the entropic barrier model of Sadler and Gilmer. This has motivated us to reexamine the model used by Sadler and Gilmer. We again find a fixed-point attractor which describes the dynamical convergence of the crystal thickness to a value just larger than the minimum stable thickness, l_min. This convergence arises from the combined effect of two constraints on the length of stems in a layer: it is unfavourable for a stem to be shorter than l_min and for a stem to overhang the edge of the previous layer. The relationship between this new mechanism and the explanation given by Sadler and Gilmer in terms of an entropic barrier is discussed. We also examine the behaviour of the Sadler-Gilmer model when an energetic contribution from chain folds is included., Comment: 15 pages, 13 figures, revtex

We consider the effect of intermolecular interactions on the optimal size-distribution of $N$ hard spheres that occupy a fixed total volume. When we minimize the free-energy of this system, within the Percus-Yevick approximation, we find that no solution exists beyond a quite low threshold ($\eta \thickapprox 0.260$). Monte Carlo simulations reveal that beyond this density, the size-distribution becomes bi-modal. Such distributions cannot be reproduced within the Percus-Yevick approximation. We present a theoretical argument that supports the occurrence of a non-monotonic size-distribution and emphasizing the importance of finite size effects., Comment: 21 pages, 6 figures, RevTex (Submitted to J. Chem. Phys.)

Based upon kinetic Monte Carlo simulations of crystallization in a simple polymer model we present a new picture of the mechanism by which the thickness of lamellar polymer crystals is constrained to a value close to the minimum thermodynamically stable thickness, l_{min}. The free energetic costs of the polymer extending beyond the edges of the previous crystalline layer and of a stem being shorter than l_{min} provide upper and lower constraints on the length of stems in a new layer. Their combined effect is to cause the crystal thickness to converge dynamically to a value close to l_{min} where growth with constant thickness then occurs. This description contrasts with those given by the two dominant theoretical approaches. However, at small supercoolings the rounding of the crystal profile does inhibit growth as suggested in Sadler and Gilmer's entropic barrier model., Comment: 12 pages, 13 figures, revtex