Researcher profile

T. Schilling

T. Schilling contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
12works
0followers
3topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

12 published item(s)

preprint2016arXiv

Monolayers of hard rods on planar substrates: I. Equilibrium

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.

preprint2013arXiv

Depletion interaction effects on the tunneling conductivity of nanorod suspensions

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.

preprint2013arXiv

Mixtures of anisotropic and spherical colloids: Phase behavior, confinement, percolation phenomena and kinetics

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.

preprint2012arXiv

Description of hard sphere crystals and crystal-fluid interfaces: a critical comparison between density functional approaches and a phase field crystal model

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.

preprint2012arXiv

Tunneling conductivity in composites of attractive colloids

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 $σ$ 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 $ϕ$ of the conducting fillers does not determine the behavior of $σ$ at larger concentrations, in striking contrast to what is predicted by percolation theory. In particular, we find that at low $ϕ$ the conductivity is governed almost entirely by the stickiness of the attraction, while at larger $ϕ$ values $σ$ 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 $ϕ$ 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 $ϕ$ is varied. Our predictions could be tested in experiments by changing the potential profile with different depletants in polymer nanocomposites.

preprint2011arXiv

2D versus 3D Freezing of a Lennard-Jones Fluid in a Slit Pore

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.

preprint2011arXiv

How Close to Two Dimensions Does a Lennard-Jones System Need to Be to Produce a Hexatic Phase?

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.

preprint2010arXiv

Free energies, vacancy concentrations and density distribution anisotropies in hard--sphere crystals: A combined density functional and simulation study

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.

preprint2010arXiv

Precursor-mediated crystallization process in suspensions of hard spheres

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.

preprint2009arXiv

Hard sphere fluids in annular wedges: density distributions and depletion potentials

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.

preprint2009arXiv

Osmotic compression of droplets of hard rods: A computer simulation study

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.

preprint2009arXiv

Solid-solid phase transition in hard ellipsoids

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.