Search

Growth of hard--rod monolayers via deposition is studied in a lattice model using rods with discrete orientations and in a continuum model with hard spherocylinders. The lattice model is treated with kinetic Monte Carlo simulations and dynamic density functional theory while the continuum model is studied by dynamic Monte Carlo simulations equivalent to diffusive dynamics. The evolution of nematic order (excess of upright particles, "standing--up" transition) is an entropic effect and is mainly governed by the equilibrium solution, {rendering a continuous transition} (paper I, J. Chem. Phys. 145, 074902 (2016)). Strong non--equilibrium effects (e.g. a noticeable dependence on the ratio of rates for translational and rotational moves) are found for attractive substrate potentials favoring lying rods. Results from the lattice and the continuum models agree qualitatively if the relevant characteristic times for diffusion, relaxation of nematic order and deposition are matched properly. Applicability of these monolayer results to multilayer growth is discussed for a continuum--model realization in three dimensions where spherocylinders are deposited continuously onto a substrate via diffusion., Comment: 22 pages, 14 figures

We analyze the structure of the state space of chess by means of transition path sampling Monte Carlo simulation. Based on the typical number of moves required to transpose a given configuration of chess pieces into another, we conclude that the state space consists of several pockets between which transitions are rare. Skilled players explore an even smaller subset of positions that populate some of these pockets only very sparsely. These results suggest that the usual measures to estimate both, the size of the state space and the size of the tree of legal moves, are not unique indicators of the complexity of the game, but that topological considerations are equally important.

The equilibrium properties of hard rod monolayers are investigated in a lattice model (where position and orientation of a rod are restricted to discrete values) as well as in an off--lattice model featuring spherocylinders with continuous positional and orientational degrees of freedom. Both models are treated using density functional theory and Monte Carlo simulations. Upon increasing the density of rods in the monolayer, there is a continuous ordering of the rods along the monolayer normal ("standing up" transition). The continuous transition also persists in the case of an external potential which favors flat--lying rods in the monolayer. This behavior is found in both the lattice and the continuum model. For the lattice model, we find very good agreement between the results from the specific DFT used (lattice fundamental measure theory) and simulations. The properties of lattice fundamental measure theory are further illustrated by the phase diagrams of bulk hard rods in two and three dimensions., Comment: 15 pages

We study by simulation and theory how the addition of insulating spherical particles affects the conductivity of fluids of conducting rods, modeled by spherocylinders. The electrical connections are implemented as tunneling processes, leading to a more detailed and realistic description than a discontinuous percolation approach. We find that the spheres enhance the tunneling conductivity for a given concentration of rods and that the enhancement increases with rod concentration into the regime where the conducting network is well established. By reformulating the network of rods using a critical path analysis, we quantify the effect of depletion-induced attraction between the rods due to the spheres. Furthermore, we show that our conductivity data are quantitatively reproduced by an effective medium approximation, which explicitly relates the system tunneling conductance to the structure of the rod-sphere fluid., Comment: 8 pages, 8 figures

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.

In materials science the phase field crystal approach has become popular to model crystallization processes. Phase field crystal models are in essence Landau-Ginzburg-type models, which should be derivable from the underlying microscopic description of the system in question. We present a study on classical density functional theory in three stages of approximation leading to a specific phase field crystal model, and we discuss the limits of applicability of the models that result from these approximations. As a test system we have chosen the three--dimensional suspension of monodisperse hard spheres. The levels of density functional theory that we discuss are fundamental measure theory, a second-order Taylor expansion thereof, and a minimal phase-field crystal model. We have computed coexistence densities, vacancy concentrations in the crystalline phase, interfacial tensions and interfacial order parameter profiles, and we compare these quantities to simulation results. We also suggest a procedure to fit the free parameters of the phase field crystal model., Comment: 21 pages

In conductor-insulator nanocomposites in which conducting fillers are dispersed in an insulating matrix the electrical connectedness is established by interparticle tunneling or hopping processes. These systems are intrinsically non-percolative and a coherent description of the functional dependence of the conductivity $\sigma$ on the filler properties, and in particular of the conductor-insulator transition, requires going beyond the usual continuum percolation approach by relaxing the constraint of a fixed connectivity distance. In this article we consider dispersions of conducting spherical particles which are connected to all others by tunneling conductances and which are subjected to an effective attractive square well potential. We show that the conductor-insulator transition at low contents $\phi$ of the conducting fillers does not determine the behavior of $\sigma$ at larger concentrations, in striking contrast to what is predicted by percolation theory. In particular, we find that at low $\phi$ the conductivity is governed almost entirely by the stickiness of the attraction, while at larger $\phi$ values $\sigma$ depends mainly on the depth of the potential well. As a consequence, by varying the range and depth of the potential while keeping the stickiness fixed, composites with similar conductor-insulator transitions may display conductivity variations of several orders of magnitude at intermediate and large $\phi$ values. By using a recently developed effective medium theory and the critical path approximation we explain this behavior in terms of dominant tunneling processes which involve interparticle distances spanning different regions of the square-well fluid structure as $\phi$ is varied. Our predictions could be tested in experiments by changing the potential profile with different depletants in polymer nanocomposites., Comment: 10 pages, 8 figures

We report on a computer simulation study of a Lennard-Jones liquid confined in a narrow slit pore with tunable attractive walls. In order to investigate how freezing in this system occurs, we perform an analysis using different order parameters. Although some of the parameters indicate that the system goes through a hexatic phase, other parameters do not. This shows that to be certain whether a system has a hexatic phase, one needs to study not only a large system, but also several order parameters to check all necessary properties. We find that the Binder cumulant is the most reliable one to prove the existence of a hexatic phase. We observe an intermediate hexatic phase only in a monolayer of particles confined such that the fluctuations in the positions perpendicular to the walls are less then 0.15 particle diameters, i. e. if the system is practically perfectly 2d.

We perform a comparative study of the free energies and the density distributions in hard sphere crystals using Monte Carlo simulations and density functional theory (employing Fundamental Measure functionals). Using a recently introduced technique (Schilling and Schmid, J. Chem. Phys 131, 231102 (2009)) we obtain crystal free energies to a high precision. The free energies from Fundamental Measure theory are in good agreement with the simulation results and demonstrate the applicability of these functionals to the treatment of other problems involving crystallization. The agreement between FMT and simulations on the level of the free energies is also reflected in the density distributions around single lattice sites. Overall, the peak widths and anisotropy signs for different lattice directions agree, however, it is found that Fundamental Measure theory gives slightly narrower peaks with more anisotropy than seen in the simulations. Among the three types of Fundamental Measure functionals studied, only the White Bear II functional (Hansen-Goos and Roth, J. Phys.: Condens. Matter 18, 8413 (2006)) exhibits sensible results for the equilibrium vacancy concentration and a physical behavior of the chemical potential in crystals constrained by a fixed vacancy concentration., Comment: 17 pages, submitted to Phys. Rev. E

We report on a large scale computer simulation study of crystal nucleation in hard spheres. Through a combined analysis of real and reciprocal space data, a picture of a two-step crystallization process is supported: First dense, amorphous clusters form which then act as precursors for the nucleation of well-ordered crystallites. This kind of crystallization process has been previously observed in systems that interact via potentials that have an attractive as well as a repulsive part, most prominently in protein solutions. In this context the effect has been attributed to the presence of metastable fluid-fluid demixing. Our simulations, however, show that a purely repulsive system (that has no metastable fluid-fluid coexistence) crystallizes via the same mechanism., Comment: 4 figures

We present a computer simulation study of a (6,12)-Lennard-Jones fluid confined to a slit pore, formed by two uniform planes. These interact via a (3,9)-Lennard-Jones potential with the fluid particles. When the fluid approaches the liquid-to-solid transition we first observe layering parallel to the walls. In order to investigate the nature of the freezing transition, we performed a detailed analysis of the bond-orientational order parameter in the layers. Changing the distance between the slit wall, we studied howthe behaviour changes from a quasi-2d crystallization without hysteresis and an intermediate hexatic phase to a standard 3d nucleation., Comment: This paper has been withdrawn by the author due to more complete study presented in arXiv:1104.0611

We present a computer simulation study of the crystalline phases of hard ellipsoids of revolution. A previous study [Phys. Rev. E, \textbf{75}, 020402 (2007)] showed that for aspect ratios $a/b\ge 3$ the previously suggested stretched-fcc phase [Mol. Phys., \textbf{55}, 1171 (1985)] is unstable with respect to a simple monoclinic phase with two ellipsoids of different orientations per unit cell (SM2). In order to study the stability of these crystalline phases at different aspect ratios and as a function of density we have calculated their free energies by thermodynamic integration. The integration path was sampled by an expanded ensemble method in which the weights were adjusted by the Wang-Landau algorithm. We show that for aspect ratios $a/b\ge 2.0$ the SM2 structure is more stable than the stretched-fcc structure for all densities above solid-nematic coexistence. Between $a/b=1.55$ and $a/b=2.0$ our calculations reveal a solid-solid phase transition.

We analyze the density distribution and the adsorption of solvent hard spheres in an annular slit formed by two large solute spheres or a large solute and a wall at close distances by means of fundamental measure density functional theory, anisotropic integral equations and simulations. We find that the main features of the density distribution in the slit are described by an effective, two--dimensional system of disks in the vicinity of a central obstacle. For large solute--solvent size ratios, the resulting depletion force has a straightforward geometrical interpretation which gives a precise "colloidal" limit for the depletion interaction. For intermediate size ratios 5...10 and high solvent packing fractions larger than 0.4, the explicit density functional results show a deep attractive well for the depletion potential at solute contact, possibly indicating demixing in a binary mixture at low solute and high solvent packing fraction., Comment: 39 pages

We present a computer simulation study of a (6,12)-Lennard-Jones fluid confined to a slit pore, formed by two uniform planes. These interact via (3,9)-Lennard-Jones potential with the fluid particles. When the fluid approaches the liquid-to-solid transition we first observe layering parallel to the walls. In order to investigate the nature of the freezing transition we performed a detailed analysis of the bond-orientational order parameter in the layers. We found no signs of hexatic order which would indicate a melting scenario of the Kosterlitz-Thouless type. An analysis of the mean-square displacement shows that the particles can easily move between the layers, making the crystallization a 3d-like process. This is consistent with the fact that we observe a considerable hysteresis in the heating-freezing curves, showing that the crystallization transition proceeds as an activated process., Comment: submitter to J. Phys.: Condens. Matter

By means of computer simulations we study how droplets of hard, rod-like particles optimize their shape and internal structure under the influence of the osmotic compression caused by the presence of spherical particles that act as depletion agents. At sufficiently high osmotic pressures the rods that make up the drops spontaneously align to turn them into uniaxial nematic liquid crystalline droplets. The nematic droplets or "tactoids" that are formed this way are not spherical but elongated, resulting from the competition between the anisotropic surface tension and the elastic deformation of the director field. In agreement with recent theoretical predictions we find that sufficiently small tactoids have a uniform director field, whilst large ones are characterized by a bipolar director field. From the shape and director-field transformation of the droplets we are able to estimate the surface anchoring strength and an average of the elastic constants of the hard-rod nematic.

When non-adsorbing polymers are added to an isotropic suspension of rod-like colloids, the colloids effectively attract each other via depletion forces. We performed Monte Carlo simulations to study the phase diagram of such rod-polymer mixture. The colloidal rods were modelled as hard spherocylinders; the polymers were described as spheres of the same diameter as the rods. The polymers may overlap with no energy cost, while overlap of polymers and rods is forbidden. Large amounts of depletant cause phase separation of the mixture. We estimated the phase boundaries of isotropic-isotropic coexistence both, in the bulk and in confinement. To determine the phase boundaries we applied the grand canonical ensemble using successive umbrella sampling [J. Chem. Phys. 120, 10925 (2004)], and we performed a finite-size scaling analysis to estimate the location of the critical point. The results are compared with predictions of the free volume theory developed by Lekkerkerker and Stroobants [Nuovo Cimento D 16, 949 (1994)]. We also give estimates for the interfacial tension between the coexisting isotropic phases and analyse its power-law behaviour on approach of the critical point.

Above a certain density threshold, suspensions of rod-like colloidal particles form system-spanning networks. Using Monte Carlo simulations, we investigate how the depletion forces caused by spherical particles affect these networks in isotropic suspensions of rods. Although the depletion forces are strongly anisotropic and favor alignment of the rods, the percolation threshold of the rods decreases significantly. The relative size of the effect increases with the aspect ratio of the rods. The structural changes induced in the suspension by the depletant are characterized in detail and the system is compared to an ideal fluid of freely interpenetrable rods., Comment: 4 pages, 3 figures

We present a study of the isotropic-nematic interface in a system of hard spherocylinders. First we compare results from Monte Carlo simulations and Onsager density functional theory for the interfacial profiles of the orientational order parameter and the density. Those interfacial properties that are not affected by capillary waves are in good agreement, despite the fact that Onsager theory overestimates the coexistence densities. Then we show results of a Monte Carlo study of the capillary waves of the interface. In agreement with recent theoretical investigations (Eur.Phys.J. E {\bf 18} 407 (2005)) we find a strongly anistropic capillary wave spectrum. For the wave-numbers accessed in our simulations, the spectrum is quadratic, i.e.elasticity does not play a role. We conjecture that this effect is due to the strong bending rigidity of the director field in suspensions of spherocylinders., Comment: 8 pages, 10 figures

Coexistence between the isotropic and the nematic phase in suspensions of rods is studied using grand canonical Monte Carlo simulations with a bias on the nematic order parameter. The biasing scheme makes it possible to estimate the interfacial tension gamma in systems of hard and soft rods. For hard rods with L/D=15, we obtain gamma ~ 1.4 kB T/L^2, with L the rod length, D the rod diameter, T the temperature, and kB the Boltzmann constant. This estimate is in good agreement with theoretical predictions, and the order of magnitude is consistent with experiments., Comment: 10 pages, 10 figures

The isotropic to nematic transition in a system of soft spherocylinders is studied by means of grand canonical Monte Carlo simulations. The probability distribution of the particle density is used to determine the coexistence density of the isotropic and the nematic phases. The distributions are also used to compute the interfacial tension of the isotropic--nematic interface, including an analysis of finite size effects. Our results confirm that the Onsager limit is not recovered until for very large elongation, exceeding at least L/D=40, with L the spherocylinder length and D the diameter. For smaller elongation, we find that the interfacial tension increases with increasing L/D, in agreement with theoretical predictions., Comment: 8 pages, 7 figures, and also 1 table

The shell structure of superheavy nuclei is investigated within various parametrizations of relativistic and nonrelativistic nuclear mean field models. The heaviest known even-even nucleus 264Hs is used as a benchmark to estimate the predictive value of the models. From that starting point, doubly magic spherical nuclei are searched in the region Z=110-140 and N=134-298. They are found at (Z=114, N=184), (Z=120, N=172), or at (Z=126, N=184), depending on the parametrization., Comment: 16 pages RevTeX, 2 tables, 2 low resolution Gif figures (high resolution PostScript versions are available at http://www.th.physik.uni-frankfurt.de/~bender/nucl_struct_publications.html or at ftp://th.physik.uni-frankfurt.de/pub/bender ), submitted to Phys. Rev. C

Causality requires, that the (anti-) commutator of two interacting field operators vanishes for space-like coordinate differences. This implies, that the Fourier transform of the spectral function of this quantum field should vanish in the space-like domain. We find, that this requirement imposes some constraints on the use of resummed propagators in high temperature gauge theory., Comment: 15 pages, latex, 3 figures included. Major rewriting of paper. Available also at ftp://tpri6c.gsi.de/pub/phenning/h95comm/ or at http://www.gsi.de/~phenning/pahtft.html