Search

We present a study of pedestrian motion along a corridor in a non-panic regime (e.g., schools, hospitals or airports). Such situations have been discussed so far within the Social Force Model (SFM). We suggest to enrich this model by interactions based on the velocity of the particles and some randomness, both of which we introduce using the ideas of the Vicsek Model (VM). This new model allows to introduce fluctuations for a given average speed and geometry, and considering that the alignment interactions are modulated by an external control parameter (the noise $\eta$) allows to introduce phase transitions between ordered and disordered states. We have compared simulations of pedestrian motion along a corridor using (a) the VM with two boundary conditions (periodic and bouncing back) and with or without desired direction of motion, (b) the SFM, and (c) the new model SFM+VM. The study of steady-state configurations in the VM with confined geometry shows the expected bands perpendicular to the motion direction, while in the SFM and SFM+VM particles order in stripes of a given width $w$ along the direction of motion. The results in the SFM+VM case show that $w(t)\simeq t^\alpha$ has a diffusive-like behavior at low noise $\eta$ (dynamic exponent $\alpha \approx 1/2$), while it is sub-diffusive at high values of external noise ($\alpha \approx 1/4$). We observe the order-disorder transition in the VM with both boundary conditions, but the application of a desired direction condition inhibits the existence of disorder as expected. For the SFM+VM case we find a susceptibility maximum which increases with system size as a function of noise strength indicative of a order-disorder transition in the whole range of densities and speeds studied. From our results we conclude that the new SFM+VM model is a well-suited model to describe non-panic evacuation with diverse degrees of disorder., Comment: 8 pages, 7 figures

The persistence length of macromolecules is one of their basic characteristics, describing their intrinsic local stiffness. However, it is difficult to extract this length from physical properties of the polymers, different recipes may give answers that disagree with each other. Monte Carlo simulations are used to elucidate this problem, giving a comparative discussion of two lattice models, the self-avoiding walk model extended by a bond bending energy, and bottle-brush polymers described by the bond fluctuation model. The conditions are discussed under which a description of such macromolecules by Kratky-Porod worm-like chains holds, and the question to what extent the persistence length depends on external conditions (such as solvent quality) is considered. The scattering function of semiflexible polymers is discussed in detail, a comparison to various analytic treatments is given, and an outlook to experimental work is presented., Comment: 18 pages, 17 figures, to be published in Polymer Science Ser. C (2013)

Using the pruned-enriched Rosenbluth Monte Carlo algorithm, the scattering functions of semiflexible macromolecules in dilute solution under good solvent conditions are estimated both in $d=2$ and $d=3$ dimensions, considering also the effect of stretching forces. Using self-avoiding walks of up to $N = 25600$ steps on the square and simple cubic lattices, variable chain stiffness is modeled by introducing an energy penalty $\epsilon_b$ for chain bending; varying $q_b=\exp (- \epsilon_b/k_BT)$ from $q_b=1$ (completely flexible chains) to $q_b = 0.005$, the persistence length can be varied over two orders of magnitude. For unstretched semiflexible chains we test the applicability of the Kratky-Porod worm-like chain model to describe the scattering function, and discuss methods for extracting persistence length estimates from scattering. While in $d=2$ the direct crossover from rod-like chains to self-avoiding walks invalidates the Kratky-Porod description, it holds in $d=3$ for stiff chains if the number of Kuhn segments $n_K$ does not exceed a limiting value $n^*_K$ (which depends on the persistence length). For stretched chains, the Pincus blob size enters as a further characteristic length scale. The anisotropy of the scattering is well described by the modified Debye function, if the actual observed chain extension $$ (end-to-end distance in the direction of the force) as well as the corresponding longitudinal and transverse linear dimensions $ - ^2$, $$ are used., Comment: 23 pages, 16 figures, 1 table

By large-scale Monte Carlo simulations of semiflexible polymers in $d=2$ dimensions the applicability of the Kratky-Porod model is tested. This model is widely used as "standard model" for describing conformations and force versus extension curves of stiff polymers. It is shown that semiflexible polymers in $d=2$ show a crossover from hard rods to self-avoiding walks, the intermediate Gaussian regime (implied by the Kratky-Porod model) is completely absent. Hence the latter can also describe force versus extension curves only if the contour length is only a few times larger than the persistence length. Consequences for experiments on biopolymers at interfaces are briefly discussed., Comment: 6 pages, 4 figures

By Molecular Dynamics simulation of a coarse-grained bead-spring type model for a cylindrical molecular brush with a backbone chain of $N_b$ effective monomers to which with grafting density $\sigma$ side chains with $N$ effective monomers are tethered, several characteristic length scales are studied for variable solvent quality. Side chain lengths are in the range $5 \le N \le 40$, backbone chain lengths are in the range $50 \le N_b \le 200$, and we perform a comparison to results for the bond fluctuation model on the simple cubic lattice (for which much longer chains are accessible, $N_b \le 1027$, and which corresponds to an athermal, very good, solvent). We obtain linear dimensions of side chains and the backbone chain and discuss their $N$-dependence in terms of power laws and the associated effective exponents. We show that even at the Theta point the side chains are considerably stretched, their linear dimension depending on the solvent quality only weakly. Effective persistence lengths are extracted both from the orientational correlations and from the backbone end-to-end distance; it is shown that different measures of the persistence length (which would all agree for Gaussian chains) are not mutually consistent with each other, and depend distinctly both on $N_b$ and the solvent quality. A brief discussion of pertinent experiments is given., Comment: 30 pages, 13 figures, 1 table

Obtaining reliable estimates of the statistical properties of complex macromolecules by computer simulation is a task that requires high computational effort as well as the development of highly efficient simulation algorithms. We present here an algorithm combining local moves, the pivot algorithm, and an adjustable simulation lattice box for simulating dilute systems of bottle-brush polymers with a flexible backbone and flexible side chains under good solvent conditions. Applying this algorithm to the bond fluctuation model, very precise estimates of the mean square end-to-end distances and gyration radii of the backbone and side chains are obtained, and the conformational properties of such a complex macromolecule are studied. Varying the backbone length (from $N_b=67$ to $N_b=1027$), side chain length (from N=0 to N=24 or 48), the scaling predictions for the backbone behavior as well as the side chain behavior are checked. We are also able to give a direct comparison of the structure factor between experimental data and the simulation results., Comment: 9 pages, 10 figures

When two types of side chains (A,B) are densely grafted to a (stiff) backbone and the resulting bottle-brush polymer is in a solution under poor solvent conditions, an incompatibility between A and B leads to microphase separation in the resulting cylindrical brush. The possible types of ordering are reminiscent of the ordering of block copolymers in cylindrical confinement. Starting from this analogy, Leibler's theory of microphase separation in block copolymer melts is generalized to derive a description of the system in the weak segregation limit. Also molecular dynamics simulation results of a corresponding coarse-grained bead-spring model are presented. Using side chain lengths up to N = 50 effective monomers, the ratio of the Lennard-Jones energy parameter between unlike monomers $(\epsilon_{AB})$ and monomers of the same kind $(\epsilon _{AA} = \epsilon_{BB})$ is varied. Various correlation functions are analyzed to study the conditions when (local) Janus cylinder-type ordering and when (local) microphase separation in the direction along the cylinder axis occurs. Both the analytical theory and the simulations give evidence for short range order due to a tendency towards microphase separation in the axial direction, with a wavelength proportional to the side chain gyration radius, irrespective of temperature and grafting density, for a wide range of these parameters., Comment: 26 pages, 19 figures

When the local intrinsic stiffness of a polymer chain varies over a wide range, one can observe both a crossover from rigid-rod-like behavior to (almost) Gaussian random coils and a further crossover towards self-avoiding walks in good solvents. Using the pruned-enriched Rosenbluth method (PERM) to study self-avoiding walks of up to $N_b=50000$ steps and variable flexibility, the applicability of the Kratky-Porod model is tested. Evidence for non-exponential decay of the bond-orientational correlations $<\cos \theta (s) >$ for large distances $s$ along the chain contour is presented, irrespective of chain stiffness. For bottle-brush polymers on the other hand, where experimentally stiffness is varied via the length of side-chains, it is shown that these cylindrical brushes (with flexible backbones) are not described by the Kratky-Porod wormlike chain model, since their persistence length is (roughly) proportional to their cross-sectional radius, for all conditions of practical interest., Comment: 6 pages, 5 figures, to be published in Europhys. Lett. (2010)

The adsorption of a bottle-brush polymer end-grafted with one chain end of its backbone to a flat substrate surface is studied by Monte Carlo simulation of a coarse-grained model, that previously has been characterized in the bulk, assuming a dilute solution under good solvent conditions. Applying the bond fluctuation model on the simple cubic lattice, we vary the backbone chain length $N_b$ from $N_b=67$ to $N_b = 259$ effective monomeric units, the side chain length $N$ from N=6 to N=48, and the grafting density $\sigma=1$, i.e., parameters that correspond well to the experimentally accessible range. When the adsorption energy strength $\epsilon$ is varied, we find that the adsorption transition (which becomes well-defined in the limit $N_b \rightarrow \infty$, for arbitrary finite $N$) roughly occurs at the same value $\epsilon_c$ as for ordinary linear chains (N=0), at least within our statistical errors. Mean square end-to-end distances and gyration radii of the side chains are obtained, as well as the monomer density profile in the direction perpendicular to the adsorbing surface. We show that for longer side chains the adsorption of bottle-brushes is a two step process, the decrease of the perpendicular linear dimension of side chains with adsorption energy strength can even be non-monotonic. Also the behavior of the static structure factor $S(q)$ is analyzed, evidence for a quasi-two-dimensional scaling is presented, and consequences for the interpretation of experiments are discussed., Comment: 18 pages, 21 figures, to be published in J. Chem. Phys. (2010)

Many small proteins fold via a first-order "all-or-none" transition directly from an expanded coil to a compact native state. Here we study an analogous direct freezing transition from an expanded coil to a compact crystallite for a simple flexible homopolymer. Wang-Landau sampling is used to construct the 1D density of states for square-well chains of length 128. Analysis within both the micro-canonical and canonical ensembles shows that, for a chain with sufficiently short-range interactions, the usual polymer collapse transition is preempted by a direct freezing or "folding" transition. A 2D free-energy landscape, built via subsequent multi-canonical sampling, reveals a dominant folding pathway over a single free-energy barrier. This barrier separates a high entropy ensemble of unfolded states from a low entropy set of crystallite states and the transition proceeds via the formation of a transition-state folding nucleus. Despite the non-unique homopolymer ground state, the thermodynamics of this direct freezing transition are identical to the thermodynamics of two-state protein folding. The model chain satisfies the van't Hoff calorimetric criterion for two-state folding and an Arrhenius analysis of the folding/unfolding free energy barrier yields a Chevron plot characteristic of small proteins., Comment: 10 pages, 5 figures, prepared for Physics Procedia: 23rd Annual Workshop on Recent Developments in Computer Simulation Studies in Condensed Matter Physics, Athens, GA

One-dimensional flexible objects are abundant in physics, from polymers to vortex lines to defect lines and many more. These objects structure their environment and it is natural to assume that the influence these objects exert on their environment depends on the distance from the line-object. But how should this be defined? We argue here that there is an intrinsic length scale along the undulating line that is a measure of its "stiffness" (i.e., orientational persistence), which yields a natural way of defining the variation of physical properties normal to the undulating line. We exemplify how this normal variation can be determined from a computer simulation for the case of a so-called bottle-brush polymer, where side chains are grafted onto a flexible backbone., Comment: 4 pages, 4 figures

As an example for the interplay of structure, dynamics, and phase behavior of macromolecular systems, this article focuses on the problem of bottle-brush polymers with either rigid or flexible backbones. On a polymer with chain length $N_b$, side-chains with chain length $N$ are endgrafted with grafting density $\sigma$. Due to the multitude of characteristic length scales and the size of these polymers (typically these cylindrical macromolecules contain of the order of 10000 effective monomeric units) understanding of the structure is a challenge for experiment. But due to excessively large relaxation times (particularly under poor solvent conditions) such macromolecules also are a challenge for simulation studies. Simulation strategies to deal with this challenge, both using Monte Carlo and Molecular Dynamics Methods, will be briefly discussed, and typical results will be used to illustrate the insight that can be gained., Comment: 8 pages, 7 figures, to be published in the proceedings of "NIC Symposium 2010", eds. G. Muenster et al. (NIC, Juelich, 2010)

Extensive Monte Carlo results are presented for a lattice model of a bottle-brush polymer under good solvent or Theta solvent conditions. Varying the side chain length, backbone length, and the grafting density for a rigid straight backbone, both radial density profiles of monomers and side chain ends are obtained, as well as structure factors describing the scattering from a single side chain and from the total bottle-brush polymer. To describe the structure in the interior of a very long bottle-brush, a periodic boundary condition in the direction along the backbone is used, and to describe effects due to the finiteness of the backbone length, a second set of simulations with free ends of the backbone is performed. In the latter case, the inhomogeneity of the structure in the direction along the backbone is carefully investigated. We use these results to test various phenomenological models that have been proposed to interpret experimental scattering data for bottle-brush macromolecules. These models aim to extract information on the radial density profile of a bottle-brush from the total scattering via suitable convolution approximations. Possibilities to improve such models, guided by our simulation results, are discussed., Comment: 15 pages, 16 figures

Employing simplified models in computer simulation is on the one hand often enforced by computer time limitations but on the other hand it offers insights into the molecular properties determining a given physical phenomenon. We employ this strategy to the determination of the phase behaviour of quadrupolar fluids, where we study the influence of omitting angular degrees of freedom of molecules via an effective spherically symmetric potential obtained from a perturbative expansion. Comparing the liquid-vapor coexistence curve, vapor pressure at coexistence, interfacial tension between the coexisting phases, etc., as obtained from both the models with the full quadrupolar interactions and the (approximate) isotropic interactions, we find discrepancies in the critical region to be typically (such as in the case of carbon dioxide) of the order of 4%. However, when the Lennard-Jones parameters are rescaled such that critical temperatures and critical densities of both models coincide with the experimental results, almost perfect agreement between the above-mentioned properties of both models is obtained. This result justifies the use of isotropic quadrupolar potentials. We present also a detailed comparison of our simulations with a combined integral equation/density functional approach and show that the latter provides an accurate description except for the vicinity of the critical point., Comment: Phys. Rev. E, accepted

Monte Carlo simulations are presented for a coarse-grained model of real quadrupolar fluids. Molecules are represented by particles interacting with Lennard-Jones forces plus the thermally averaged quadrupole-quadrupole interaction. The properties discussed include the vapor-liquid coexistence curve, the vapor pressure along coexistence, and the surface tension. The full isotherms are also accessible over a wide range of temperatures and densities. It is shown that the critical parameters (critical temperature, density, and pressure) depend almost linearly on a quadrupolar parameter $q=Q^{*4} /T^*$, $Q^*$ is the reduced quadrupole moment of the molecule and $T^*$ the reduced temperature. The model can be applied to a variety of small quadrupolar molecules. We focus on carbon dioxide as a test case, but consider nitrogen and benzene, too. Experimental critical temperature, density and quadrupolar moment are sufficient to fix the parameters of the model. The resulting agreement with experiments is excellent and marks a significant improvement over approaches which neglect quadrupolar effects. The same coarse-grained model was also applied in the framework of Perturbation Theory (PT) in the Mean Spherical Approximation (MSA). As expected, the latter deviates from the Monte Carlo results in the critical region, but is reasonably accurate at lower temperatures., Comment: J. Chem. Phys. (2008), to appear

A polymer chain tethered to a surface may be compact or extended, adsorbed or desorbed, depending on interactions with the surface and the surrounding solvent. This leads to a rich phase diagram with a variety of transitions. To investigate these transitions we have performed Monte Carlo simulations of a bond-fluctuation model with Wang-Landau and umbrella sampling algorithms in a two-dimensional state space. The simulations' density of states results have been evaluated for interaction parameters spanning the range from good to poor solvent conditions and from repulsive to strongly attractive surfaces. In this work, we describe the simulation method and present results for the overall phase behavior and for some of the transitions. For adsorption in good solvent, we compare with Metropolis Monte Carlo data for the same model and find good agreement between the results. For the collapse transition, which occurs when the solvent quality changes from good to poor, we consider two situations corresponding to three-dimensional (hard surface) and two-dimensional (very attractive surface) chain conformations, respectively. For the hard surface, we compare tethered chains with free chains and find very similar behavior for both types of chains. For the very attractive surface, we find the two-dimensional chain collapse to be a two-step transition with the same sequence of transitions that is observed for three-dimensional chains: a coil-globule transition that changes the overall chain size is followed by a local rearrangement of chain segments., Comment: 17 pages, 12 figures, to appear in J. Chem. Phys

Scaling predictions and results from self-consistent field calculations for bottle-brush polymers with a rigid backbone and flexible side chains under good solvent conditions are summarized and their validity and applicability is assessed by a comparison with Monte Carlo simulations of a simple lattice model. It is shown that under typical conditions, as they are also present in experiments, only a rather weak stretching of the side chains is realized, and then the scaling predictions based on the extension of the Daoud-Cotton blob picture are not applicable. Also two-component bottle brush polymers are considered, where two types (A,B) of side chains are grafted, assuming that monomers of different kind repel each other. In this case, variable solvent quality is allowed for, such that for poor solvent conditions rather dense cylinder-like structures result. Theories predict ``Janus Cylinder''-type phase separation along the backbone in this case. The Monte Carlo simulations, using the pruned-enriched Rosenbluth method (PERM) then are restricted to rather short side chain length. Nevertheless, evidence is obtained that the phase separation between an A-rich part of the cylindrical molecule and a B-rich part can only occur locally. The correlation length of this microphase separation can be controlled by the solvent quality. This lack of a phase transition is interpreted by an analogy with models for ferromagnets in one space dimension., Comment: 27pages, 33 figures

A lattice model for a symmetrical copolymer "bottle brush" molecule, where two types (A,B) of flexible side chains are grafted with one chain end to a rigid backbone, is studied by a variant of the pruned-enriched Rosenbluth method (PERM), allowing for simultaneous growth of all side chains in the Monte Carlo sampling. Choosing repulsive binary interactions between unlike monomers and varying the solvent quality, it is found that phase separation into an $A$-rich part of the cylindrical molecule and a $B$-rich part can occur only locally. Long range order (in the direction of the backbone) does not occur, and hence the transition from the randomly mixed state of the bottle brush to the phase separated structure is strongly rounded, in contrast to corresponding mean field predictions. This lack of a phase transition can be understood from an analogy with spin models in one space dimension. We predict that the range of microphase separation along the bottle brush backbone can be controlled on the nanoscale by varying the solvent quality., Comment: 5 pages, 7 figures. Edited for typos