BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Tanja Schilling

Tanja Schilling appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
cond-mat.softcond-mat.stat-mechphysics.comp-phcond-mat.mtrl-sciComputer Visioncond-mat.dis-nncond-mat.mes-hallmath-phmath.MPphysics.atm-clusquant-ph

Catalog footprint

What is connected

34works
11topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

34 published item(s)

preprint2022arXiv

Hard sphere crystal nucleation rates: Reconciliation of simulation and experiment

Over the past two decades, a large number of studies addressed the topic of crystal nucleation in suspensions of hard spheres. The shared result of all these efforts is that, at low super-saturations, experimentally observed nucleation rates and numerically computed ones differ by more than ten orders of magnitude. We present precise simulation results of crystal nucleation rate densities in the meta-stable hard sphere liquid. To compare these rate densities to experimentally measured ones, we propose an interpretation of the experimental data as a combination of nucleation and crystal growth processes (rather than purely the nucleation process). This interpretation may resolve the long standing dispute about the differing rates.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

The Interplay between Memory and Potentials of Mean Force: A Discussion on the Structure of Equations of Motion for Coarse Grained Observables

The underdamped, non-linear, generalized Langevin equation is widely used to model coarse-grained dynamics of soft and biological materials. By means of a projection operator formalism, we show under which approximations this equation can be obtained from the Hamiltonian dynamics of the underlying microscopic system and in which cases it makes sense to introduce a potential of mean force. We discuss shortcomings of previous derivations presented in the literature and demonstrate the implications of our derivation for the structure of memory terms and their connection to generalized fluctuation-dissipation relations. We show, in particular, that the widely used, simple structure which contains a potential of mean force, a memory term which is linear in the observable, and a fluctuating force which is related to the memory term by a fluctuation-dissipation relation, is neither exact nor can it, in general, be derived as a controlled approximation to the exact dynamics.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Coarse-Grained Modelling Out of Equilibrium

Active matter, responsive ("smart") materials and materials under time-dependent load are systems out of thermal equilibrium. To construct coarse-grained models for such systems, one needs to integrate out a distribution of microstates that evolves in time. This is a challenging task. As a preparation to the topic, we recall equilibrium coarse-graining methods, both theoretical and numerical, such as projection operator formalisms, united atom simulations, numerical reconstruction of memory and Markov State Modelling. Then we review recent developments in theoretical approaches to the non-equilibrium coarse-graining problem, in particular, time-dependent projection operator formalisms, dynamic density functional theory and power functional theory, as well as numerical schemes to contruct explicitly time-dependent memory kernels.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Evaluation of memory effects at phase transitions and during relaxation processes

We propose to describe the dynamics of phase transitions in terms of a non-stationary Generalized Langevin Equation for the order parameter. By construction, this equation is non-local in time, i.e.~it involves memory effects whose intensity is governed by a memory kernel. In general, it is a hard task to determine the physical origin and the extent of the memory effects based on the underlying microscopic equations of motion. Therefore we propose to relate the extent of the memory kernel to quantities that are experimentally observed such as the induction time and the duration of the phase transformation process. Using a simple kinematic model, we show that the extent of the memory kernel is positively correlated with the duration of the transition and of the same order of magnitude, while the distribution of induction times does not have an effect. This observation is tested at the example of several model systems, for which we have run Molecular Dynamics simulations: a modified Potts-model, a dipole gas, an anharmonic spring in a bath and a nucleation problem. All these cases are shown to be consistent with the simple theoretical model.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Shape, connectedness percolation and electrical conductivity of clusters in suspensions of hard platelets

Using Monte Carlo simulations, we investigate how geometric percolation and electrical conductivity in suspensions of hard conducting platelets are affected by the addition of platelets and their degree of spontaneous alignment. For aspect ratios $10, 25$ and $50$, we consistently observe a monotonically decreasing percolation threshold as a function of volume fraction. In the nematic phase, the distribution of particles inside the percolating clusters becomes less spherically symmetric and the aspect ratio of the clusters increases. However, the clusters are also anisotropically shaped in the isotropic phase, although their aspect ratio remains constant as a function of volume fraction. Mapping the percolating clusters of platelets to linear resistor networks, and assigning unit conductance to all connections, we find a constant conductivity both across the isotropic-nematic transition and in the respective stable phases. This behaviour is consistent with the other observed topological properties of the networks. On the contrary, using an anisotropic conductance model, the network conductivity decreases with increasing volume fraction in the isotropic, and further diminishes at the onset of the nematic. Hence, our observations consistently suggest that unlike for rod-like fillers, the network structures that arise from platelet suspensions are neither very sensitive to the particle aspect ratio nor to alignment, thus rendering platelets less versatile fillers for dispersion in conductive composites.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Structural localization in the Classical and Quantum Fermi-Pasta-Ulam Model

We study the statistics and short-times dynamics of the classical and the quantum Fermi-Pasta-Ulam chain in thermal equilibrium. We analyze the distributions of single-particle configurations by integrating out the rest of the system. At low temperatures we observe a systematic increase in the mobility of the chain when transitioning from classical to quantum mechanics due to zero-point energy effects. We analyze the consequences of the quantum dispersion on the dynamics at short times of configurational correlation functions.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Continuum percolation expressed in terms of density distributions

We present a new approach to derive the connectivity properties of pairwise interacting n-body systems in thermal equilibrium. We formulate an integral equation that relates the pair connectedness to the distribution of nearest neighbors. For one-dimensional systems with nearest-neighbor interactions, the nearest-neighbor distribution is, in turn, related to the pair correlation function g through a simple integral equation. As a consequence, for those systems, we arrive at an integral equation relating g to the pair connectedness, which is readily solved even analytically if g is specified analytically. We demonstrate the procedure for a variety of pair-potentials including fully penetrable spheres as well as impenetrable spheres, the only two systems for which analytical results for the pair connectedness exist. However, the approach is not limited to nearest-neighbor interactions in one dimension. Hence, we also outline the treatment of external fields and long-ranged interactions, and we illustrate how the formalism can applied to higher-dimensional systems using the three-dimensional ideal gas as an example.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Unusual geometric percolation of hard nanorods in the uniaxial nematic liquid crystalline phase

We investigate by means of continuum percolation theory and Monte Carlo simulations how spontaneous uniaxial symmetry breaking affects geometric percolation in dispersions of hard rod-like particles. If the particle aspect ratio exceeds about twenty, percolation in the nematic phase can be lost upon adding particles to the dispersion. This contrasts with percolation in the isotropic phase, where a minimum particle loading is always required to obtain system-spanning clusters. For sufficiently short rods, percolation in the uniaxial nematic mimics that of the isotropic phase, where the addition of particles always aids percolation. For aspect ratios between twenty and infinity, but not including infinity, we find re-entrance behavior: percolation in the low-density nematic may be lost upon increasing the amount of nanofillers but can be re-gained by the addition of even more particles to the suspension. Our simulation results for aspect ratios of 5, 10, 20, 50 and 100 strongly support our theoretical predictions, with almost quantitative agreement. We show that a new closure of the connectedness Ornstein-Zernike equation, inspired by Scaled Particle Theory, is more accurate than the Lee-Parsons closure that effectively describes the impact of many-body direct contacts.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Clearing out a maze: The hungry random walker and its anomalous diffusion

We study chemotaxis in a porous medium using as a model a biased ("hungry") random walk on a percolating cluster. In close resemblance to the 1980s arcade game Pac-Man, the hungry random walker consumes food, which is initially distributed in the maze, and biases its movement towards food-filled sites. We observe that, on the percolating cluster, the mean-squared displacement of the pacman process shows anomalous dynamics, which follow a power law with a dynamical exponent different from both that of a self avoiding random walk as well as that of an unbiased random walk. The change in dynamics with the propensity to move towards food is well described by a dynamical exponent that depends continuously on this propensity, and results in slower differential growth when compared to the unbiased random walk.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The Early Crystal Nucleation Process in Hard Spheres shows Synchronised Ordering and Densification

We investigate the early part of the crystal nucleation process in the hard sphere fluid using data produced by computer simulation. We find that hexagonal order manifests continuously in the overcompressed liquid, beginning approximately one diffusion time before the appearance of the first `solid-like' particle of the nucleating cluster, and that a collective influx of particles towards the nucleation site occurs simultaneously to the ordering process: the density increases leading to nucleation are generated by the same individual particle displacements as the increases in order. We rule out the presence of qualitative differences in the early nucleation process between medium and low overcompressions, and also provide evidence against any separation of translational and orientational order on the relevant lengthscales.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Connectivity Percolation in Suspensions of Attractive Square-Well Spherocylinders

We have studied the connectivity percolation transition in suspensions of attractive square-well spherocylinders by means of Monte Carlo simulation and connectednes percolation theory. In the 1980s the percolation threshold of slender fibres has been predicted to scale as the fibres' inverse aspect ratio (Phys.~Rev.~B {\bf 30}, 3933 (1984)). The main finding of our study is that the attractive spherocylinder system reaches this inverse scaling regime at much lower aspect ratios than found in suspensions of hard spherocylinders. We explain this difference by showing that third virial corrections of the pair connectedness functions, which are responsible for the deviation from the scaling regime, are less important for attractive potentials than for hard particles.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Dissipation in a Crystallization Process

We discuss the crystallization process from the supersaturated melt in terms of its non-equilibrium properties. In particular, we quantify the amount of heat that is produced irreversibly when a suspension of hard spheres crystallizes. This amount of heat can be interpreted as arising from the resistance of the system against undergoing phase transition. We identify an intrinsic compression rate that separates a quasi-static regime from a regime of rapid crystallization. In the former the disspated heat grows linearly in the compression rate. In the latter the system crystallizes more easily, because new relaxation channels are opened, at the cost of forming a higher fraction of non-equilibrium crystal structures. In analogy to a shear-thinning fluid, the system shows a decreased resistance when it is driven rapidly.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Non-universal Voronoi cell shapes in amorphous ellipsoid packings

In particulate systems with short-range interactions, such as granular matter or simple fluids, local structure plays a pivotal role in determining the macroscopic physical properties. Here, we analyse local structure metrics derived from the Voronoi diagram of configurations of oblate ellipsoids, for various aspect ratios $α$ and global volume fractions $ϕ_g$. We focus on jammed static configurations of frictional ellipsoids, obtained by tomographic imaging and by discrete element method simulations. In particular, we consider the local packing fraction $ϕ_l$, defined as the particle's volume divided by its Voronoi cell volume. We find that the probability $P(ϕ_l)$ for a Voronoi cell to have a given local packing fraction shows the same scaling behaviour as function of $ϕ_g$ as observed for random sphere packs. Surprisingly, this scaling behaviour is further found to be independent of the particle aspect ratio. By contrast, the typical Voronoi cell shape, quantified by the Minkowski tensor anisotropy index $β=β_0^{2,0}$, points towards a significant difference between random packings of spheres and those of oblate ellipsoids. While the average cell shape $β$ of all cells with a given value of $ϕ_l$ is very similar in dense and loose jammed sphere packings, the structure of dense and loose ellipsoid packings differs substantially such that this does not hold true. This non-universality has implications for our understanding of jamming of aspherical particles.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Percolation in suspensions of hard nanoparticles: From spheres to needles

We investigate geometric percolation and scaling relations in suspensions of nanorods, covering the entire range of aspect ratios from spheres to extremely slender needles. A new version of connectedness percolation theory is introduced and tested against specialized Monte Carlo simulations. The theory accurately predicts percolation thresholds for aspect ratios as low as 10. The percolation threshold for rod-like particles of aspect ratios below 1000 deviates significantly from the inverse aspect ratio scaling prediction, thought to be valid in the limit of infinitely slender rods and often used as a rule of thumb for nano-fibers in composite materials. Hence, most fibers that are currently used as fillers in composite materials cannot be regarded as practically infinitely slender for the purposes of percolation theory. Comparing percolation thresholds of hard rods and new benchmark results for ideal rods, we find that (i) for large aspect ratios, they differ by a factor that is inversely proportional to the connectivity distance between the hard cores, and (ii) they approach the slender rod limit differently.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Percolation in suspensions of polydisperse hard rods : quasi-universality and finite-size effects

We present a study of connectivity percolation in suspensions of hard spherocylinders by means of Monte Carlo simulation and connectedness percolation theory. We focus attention on polydispersity in the length, the diameter and the connectedness criterion, and invoke bimodal, Gaussian and Weibull distributions for these. The main finding from our simulations is that the percolation threshold shows quasi universal behaviour, i.e., to a good approximation it depends only on certain cumulants of the full size and connectivity distribution. Our connectedness percolation theory hinges on a Lee-Parsons type of closure recently put forward that improves upon the often-used second virial approximation [ArXiv e-prints, May 2015, 1505.07660]. The theory predicts exact universality. Theory and simulation agree quantitatively for aspect ratios in excess of 20, if we include the connectivity range in our definition of the aspect ratio of the particles. We further discuss the mechanism of cluster growth that, remarkably, differs between systems that are polydisperse in length and in width, and exhibits non-universal aspects.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Crystal growth from a supersaturated melt: relaxation of the solid-liquid dynamic stiffness

We discuss the growth process of a crystalline phase out of a metastable over-compressed liquid that is brought into contact with a crystalline substrate. The process is modeled by means of molecular dynamics. The particles interact via the Lennard-Jones potential and their motion is locally thermalized by Langevin dynamics. We characterize the relaxation process of the solid-liquid interface, showing that the growth speed is maximal for liquid densities above the solid coexistence density, and that the structural properties of the interface rapidly converge to equilibrium-like properties. In particular, we show that the off-equilibrium dynamic stiffness can be extracted using capillary wave theory arguments, even if the growth front moves fast compared to the typical diffusion time of the compressed liquid, and that the dynamic stiffness converges to the equilibrium stiffness in times much shorter than the diffusion time.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Crystal nucleation mechanism in melts of short polymer chains under quiescent conditions and under shear flow

We present a molecular dynamics simulation study of crystal nucleation from undercooled melts of n-alkanes, and we identify the molecular mechanism of homogeneous crystal nucleation under quiescent conditions and under shear flow. We compare results for n-eicosane(C20) and n-pentacontahectane(C150), i.e. one system below the entanglement length and one above. Under quiescent conditions, we observe that entanglement does not have an effect on the nucleation mechanism. For both chain lengths, the chains first align and then straighten locally. Then the local density increases and finally positional ordering sets in. At low shear rates the nucleation mechanism is the same as under quiescent conditions, while at high shear rates the chains align and straighten at the same time. We report on the effects of shear rate and temperature on the nucleation rates and estimate the critical shear rates, beyond which the nucleation rates increase with the shear rate. We show that the viscosity of the system is not affected by the crystalline nuclei.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Solid phase properties and crystallization in simple model systems

We review theoretical and simulational approaches to the description of equilibrium bulk crystal and interface properties as well as to the nonequilibrium processes of homogeneous and heterogeneous crystal nucleation for the simple model systems of hard spheres and Lennard-Jones particles. For the equilibrium properties of bulk and interfaces, density functional theories employing fundamental measure functionals prove to be a precise and versatile tool, as exemplified with a closer analysis of the hard sphere crystalliquid interface. A detailed understanding of the dynamic process of nucleation in these model systems nevertheless still relies on simulational approaches. We review bulk nucleation and nucleation at structured walls and examine in closer detail the influence of walls with variable strength on nucleation in the Lennard-Jones fluid. We find that a planar crystalline substrate induces the growth of a crystalline film for a large range of lattice spacings and interaction potentials. Only a strongly incommensurate substrate and a very weakly attractive substrate potential lead to crystal growth with a non-zero contact angle.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Strong Effect of Weak Charging in Suspensions of Anisotropic Colloids

Suspensions of hard colloidal particles frequently serve as model systems in studies on fundamental aspects of phase transitions. But often colloidal particles that are considered as ``hard'' are in fact weakly charged. If the colloids are spherical, weak charging has a only a weak effect on the structural properties of the suspension, which can be easily corrected for. However, this does not hold for anisotropic particles. We introduce a model for the interaction potential between charged ellipsoids of revolution (spheroids) based on the Derjaguin approximation of Debye--Hückel Theory and present a computer simulation study on aspects of the system's structural properties and phase behaviour. In line with previous experimental observations, we find that even a weak surface charge has a strong impact on the correlation functions. A likewise strong impact is seen on the phase behaviour, in particular, we find stable cubatic order in suspensions of oblate ellipsoids.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Crystallization in Glassy Suspensions of Hard Ellipsoids

We have carried out computer simulations of overcompressed suspensions of hard monodisperse ellipsoids and observed their crystallization dynamics. The system was compressed very rapidly in order to reach the regime of slow, glass-like dynamics. We find that, although particle dynamics become sub-diffusive and the intermediate scattering function clearly develops a shoulder, crystallization proceeds via the usual scenario: nucleation and growth for small supersaturations, spinodal decomposition for large supersaturations. In particular, we compared the mobility of the particles in the regions where crystallization set in with the mobility in the rest of the system. We did not find any signature in the dynamics of the melt that pointed towards the imminent crystallization events.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Crystallization mechanism in melts of short n-alkane chains

We study crystallization in a model system for eicosane (C20) by means of molecular dynamics simulation and we identify the microscopic mechanisms of homogeneous crystal nucleation and growth. For the nucleation process, we observe that chains first align and then straighten. Then the local density increases and finally the monomer units become ordered positionally. The subsequent crystal growth process is characterized by a sliding-in motion of the chains. Chains preferably attach to the crystalline cluster with one end and then move along the stems of already crystallized chains towards their final position. This process is cooperative, i.e. neighboring chains tend to get attached in clusters rather than independently

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Euclidean Upgrade from a Minimal Number of Segments

In this paper, we propose an algebraic approach to upgrade a projective reconstruction to a Euclidean one, and aim at computing the rectifying homography from a minimal number of 9 segments of known length. Constraints are derived from these segments which yield a set of polynomial equations that we solve by means of Gröbner bases. We explain how a solver for such a system of equations can be constructed from simplified template data. Moreover, we present experiments that demonstrate that the given problem can be solved in this way.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Influence of random pinning on the crystallization process in suspensions of hard spheres

We discuss crystal formation in supersaturated suspensions of monodisperse hard spheres with a concentration of hard spheres randomly pinned in space and time. The pinning procedure introduces an external length scale and an external time scale that restrict the accessible number of configurations and ultimately the number of pathways leading to crystallization. We observe a significant drop in the nucleation rate density at a characteristic pinning concentration that can be directly related to the structure of the critical nucleus and the dynamics of its formation in the unpinned system.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Solvent hydrodynamics speed up crystal nucleation in suspensions of hard spheres

We present a computer simulation study on the crystal nucleation process in suspensions of hard spheres, fully taking into account the solvent hydrodynamics. If the dynamics of collodial crystallization were purely diffusive, the crystal nucleation rate densities would drop as the inverse of the solvent viscosity. However, we observe that the nucleation rate densities do not scale in this way, but are enhanced at high viscosities. This effect might explain the large discrepancy between the nuclation rate densities obtained by simulation and experiment that have reported in the literature so far.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Connectivity percolation in suspensions of hard platelets

We present a study on connectivity percolation in suspensions of hard platelets by means of Monte Carlo simulation. We interpret our results using a contact-volume argument based on an effective single--particle cell model. It is commonly assumed that the percolation threshold of anisotropic objects scales as their inverse aspect ratio. While this rule has been shown to hold for rod-like particles, we find that for hard plate-like particles the percolation threshold is non-monotonic in the aspect ratio. It exhibits a shallow minimum at intermediate aspect ratios and then saturates to a constant value. This effect is caused by the isotropic-nematic transition pre-empting the percolation transition. Hence the common strategy to use highly anisotropic, conductive particles as fillers in composite materials in order to produce conduction at low filler concentration is expected to fail for plate-like fillers such as graphene and graphite nanoplatelets.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Free Energies by Thermodynamic Integration Relative to an Exact Solution, Used to Find the Handedness-Switching Salt Concentration for DNA

Sets of free energy differences are useful for finding the equilibria of chemical reactions, while absolute free energies have little physical meaning. However finding the relative free energy between two macrostates by subtraction of their absolute free energies is a valuable strategy in certain important cases. We present calculations of absolute free energies of biomolecules, using a combination of the well-known Einstein Molecule method (for treating the solute) with a conceptually related method of recent genesis for computing free energies of liquids (to treat the solvent and counterions). The approach is based on thermodynamic integration from a detailed atomistic model to one which is simplified but analytically solvable, thereby giving the absolute free energy as that of the tractable model plus a correction term found numerically. An example calculation giving the free energy with respect to salt concentration for the B- and Z-isomers of all-atom duplex DNA in explicit solvent and counterions is presented. The coexistence salt concentration is found with unprecedented accuracy.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Hydrodynamic interactions of colloidal spheres under shear flow

Particles that are immersed in a fluid exchange momentum via the fluid, hence their Brownian motion is correlated. By means of multiparticle-collision dynamics simulations we study the interactions between two colloidal beads in a sheared fluid suspension. Recently, this topic has been addressed in experiments on colloidal particles trapped by optical tweezers in a microfluidic device [PRL, 103, 230602 (2009)] and theoretically by means of a Langevin model [Eur. Phys. J E, 33, 313 (2010)]. Although we neglect the rotational degrees of freedom of the colloids, and employ a very simple coupling between the colloids and the flow field, we can reproduce the experimental data and partly explain why it differs from theory.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

On the influence of a patterned substrate on crystallization in suspensions of hard spheres

We present a computer simulation study on crystal nucleation and growth in supersaturated suspensions of mono-disperse hard spheres induced by a triangular lattice substrate. The main result is that compressed substrates are wet by the crystalline phase (the crystalline phase directly appears without any induction time), while for stretched substrates we observe heterogeneous nucleation. The shapes of the nucleated crystallites fluctuate strongly. In the case of homogeneous nucleation amorphous precursors have been observed (Phys. Rev. Lett. {\bf 105}(2):025701 (2010)). For heterogeneous nucleation we do not find such precursors. The fluid is directly transformed into highly ordered crystallites.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

A method to compute absolute free energies or enthalpies of fluids

We propose a new method to compute the free energy or enthalpy of fluids or disordered solids by computer simulation . The main idea is to construct a reference system by freezing one representative configuration, and then carry out a thermodynamic integration. We present a strategy and an algorithm which allows to sample the thermodynamic integration path even in the case of liquids, despite the fact that the particles can diffuse freely through the system. The method is described in detail and illustrated with applications to hard sphere fluids and solids with mobile defects.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

A General Algorithm for Sampling Rare Events in Non-Equilibrium and Non-Stationary Systems

Although many computational methods for rare event sampling exist, this type of calculation is not usually practical for general nonequilibrium conditions, with macroscopically irreversible dynamics and away from both stationary and metastable states. A novel method for calculating the time-series of the probability of a rare event is presented which is designed for these conditions. The method is validated for the cases of the Glauber-Ising model under time-varying shear flow, the Kawasaki-Ising model after a quench into the region between nucleation dominated and spinodal decomposition dominated phase change dynamics, and the parallel open asymmetric exclusion process (p-o ASEP). The method requires a subdivision of the phase space of the system: it is benchmarked and found to scale well for increasingly fine subdivisions, meaning that it can be applied without detailed foreknowledge of the physically important reaction pathways.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Crystallization in suspensions of hard spheres: A Monte Carlo and Molecular Dynamics simulation study

The crystallization of a metastable melt is one of the most important non equilibrium phenomena in condensed matter physics, and hard sphere colloidal model systems have been used for several decades to investigate this process by experimental observation and computer simulation. Nevertheless, there is still an unexplained discrepancy between simulation data and experimental nucleation rate densities. In this paper we examine the nucleation process in hard spheres using molecular dynamics and Monte Carlo simulation. We show that the crystallization process is mediated by precursors of low orientational bond-order and that our simulation data fairly match the experimental data sets.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Indium-Gallium Segregation in CuIn$_{x}$Ga$_{1-x}$Se$_2$: An ab initio based Monte Carlo Study

Thin-film solar cells with CuIn$_x$Ga$_{1-x}$Se$_2$ (CIGS) absorber are still far below their efficiency limit, although lab cells reach already 19.9%. One important aspect is the homogeneity of the alloy. Large-scale simulations combining Monte Carlo and density functional calculations show that two phases coexist in thermal equilibrium below room temperature. Only at higher temperatures, CIGS becomes more and more a homogeneous alloy. A larger degree of inhomogeneity for Ga-rich CIGS persists over a wide temperature range, which may contribute to the low observed efficiency of Ga-rich CIGS solar cells.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Computing absolute free energies of disordered structures by molecular simulation

We present a Monte Carlo simulation technique by which the free energy of disordered systems can be computed directly. It is based on thermodynamic integration. The central idea is to construct an analytically solvable reference system from a configuration which is representative for the state of interest. The method can be applied to lattice models (e.g., the Ising model) as well as off-lattice molecular models. We focus mainly on the more challenging off-lattice case. We propose a Monte Carlo algorithm, by which the thermodynamic integration path can be sampled efficiently. At the examples of the hard sphere liquid and a hard disk solid with a defect we discuss several properties of the approach.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Fluctuating Interfaces in Liquid Crystals

We review and compare recent work on the properties of fluctuating interfaces between nematic and isotropic liquid-crystalline phases. Molecular dynamics and Monte Carlo simulations have been carried out for systems of ellipsoids and hard rods with aspect ratio 15:1, and the fluctuation spectrum of interface positions (the capillary wave spectrum) has been analyzed. In addition, the capillary wave spectrum has been calculated analytically within the Landau-de Gennes theory. The theory predicts that the interfacial fluctuations can be described in terms of a wave vector dependent interfacial tension, which is anisotropic at small wavelengths (stiff director regime) and becomes isotropic at large wavelengths (flexible director regime). After determining the elastic constants in the nematic phase, theory and simulation can be compared quantitatively. We obtain good agreement for the stiff director regime. The crossover to the flexible director regime is expected at wavelengths of the order of several thousand particle diameters, which was not accessible to our simulations.

0 citations0 reviewsNo code linkedOpen access