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

Thomas M. Truskett

Thomas M. Truskett 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.mtrl-scicond-mat.stat-mechcond-mat.mes-hallphysics.chem-phphysics.app-phphysics.comp-phphysics.flu-dyn

Catalog footprint

What is connected

36works
8topics
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

36 published item(s)

preprint2022arXiv

Geometric model of crack-templated networks for transparent conductive films

Crack-templated networks, metallic frameworks fabricated from crack patterns in sacrificial thin films, can exhibit high optical transmittance, high electric conductivity, and a host of other properties attractive for applications. Despite advances in preparing, characterizing, and analyzing optoelectronic performance of cracked template networks, limited efforts have focused on predicting how their disordered structures help determine their electrical and optical properties and explain their interrelationships. We introduce a geometric modeling approach for crack-templated networks and use simulation to compute their wavelength- and incident angle-dependent optical transmittance and sheet resistivity. We explore how these properties relate to one another and to those of metallic meshes with periodically ordered aperture arrays. We consider implications of the results for optoelectronic applications, compare figure-of-merit predictions to experimental data, and highlight an opportunity to extend the modeling approach using inverse methods.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Machine learning-assisted design of material properties

Designing functional materials requires a deep search through multidimensional spaces for system parameters that yield desirable material properties. For cases where conventional parameter sweeps or trial-and-error sampling are impractical, inverse methods that frame design as a constrained optimization problem present an attractive alternative. However, even efficient algorithms require time- and resource-intensive characterization of material properties many times during optimization, imposing a design bottleneck. Approaches that incorporate machine learning can help address this limitation and accelerate the discovery of materials with targeted properties. In this article, we review how to leverage machine learning to reduce dimensionality in order to effectively explore design space, accelerate property evaluation, and generate unconventional material structures with optimal properties. We also discuss promising future directions, including integration of machine learning into multiple stages of a design algorithm and interpretation of machine learning models to understand how design parameters relate to material properties.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Treating random sequential addition via the replica method

While many physical processes are non-equilibrium in nature, the theory and modeling of such phenomena lag behind theoretical treatments of equilibrium systems. The diversity of powerful theoretical tools available to describe equilibrium systems has inspired strategies that map non-equilibrium systems onto equivalent equilibrium analogs so that interrogation with standard statistical mechanical approaches is possible. In this work, we revisit the mapping from the non-equilibrium random sequential addition process onto an equilibrium multi-component mixture via the replica method, allowing for theoretical predictions of non-equilibrium structural quantities. We validate the above approach by comparing the theoretical predictions to numerical simulations of random sequential addition.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Assembling Inorganic Nanocrystal Gels

Inorganic nanocrystal gels retain distinct properties of individual nanocrystals while offering tunable, network structure-dependent characteristics. We review different mechanisms for assembling gels from colloidal nanocrystals including (1) controlled destabilization, (2) direct bridging, (3) depletion, as well as linking mediated by (4) coordination bonding or (5) dynamic covalent bonding, and we highlight how each impacts gel properties. These approaches use nanocrystal surface chemistry or the addition of small molecules to mediate inter-nanocrystal attractions. Each method offers advantages in terms of gel stability, reversibility, or tunability and presents new opportunities for design of reconfigurable materials and fueled assemblies.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Effects of linker flexibility on phase behavior and structure of linked colloidal gels

Colloidal nanocrystal gels can be assembled using a difunctional "linker" molecule to mediate bonding between nanocrystals. The conditions for gelation and the structure of the gel are controlled macroscopically by the linker concentration and microscopically by the linker's molecular characteristics. Here, we demonstrate using a toy model for a colloid-linker mixture that linker flexibility plays a key role in determining both phase behavior and structure of the mixture. We fix the linker length and systematically vary its bending stiffness to span the flexible, semiflexible, and rigid regimes. At fixed linker concentration, flexible-linker and rigid-linker mixtures phase separate at low colloid volume fractions in agreement with predictions of first-order thermodynamic perturbation theory, but the semiflexible-linker mixtures do not. We correlate and attribute this qualitatively different behavior to undesirable "loop" linking motifs that are predicted to be more prevalent for linkers with end-to-end distances commensurate with the locations of chemical bonding sites on the colloids. Linker flexibility also influences the spacing between linked colloids, suggesting strategies to design gels with desired phase behavior, structure, and by extension, structure-dependent properties.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Inverse methods for design of soft materials

Functional soft materials, comprising colloidal and molecular building blocks that self-organize into complex structures as a result of their tunable interactions, enable a wide array of technological applications. Inverse methods provide systematic means for navigating their inherently high-dimensional design spaces to create materials with targeted properties. While multiple physically motivated inverse strategies have been successfully implemented in silico, their translation to guiding experimental materials discovery has thus far been limited to a handful of proof-of-concept studies. In this Perspective, we discuss recent advances in inverse methods for design of soft materials that address two challenges: (1) methodological limitations that prevent such approaches from satisfying design constraints and (2) computational challenges that limit the size and complexity of systems that can be addressed. Strategies that leverage machine learning have proven particularly effective, including methods to discover order parameters that characterize complex structural motifs and schemes to efficiently compute macroscopic properties from the underlying structure. We also highlight promising opportunities to improve the experimental realizability of materials designed computationally, including discovery of materials with functionality at multiple thermodynamic states, design of externally directed assembly protocols that are simple to implement in experiments, and strategies to improve the accuracy and computational efficiency of experimentally relevant models.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Stability of force-driven shear flows in nonequilibrium molecular simulations with periodic boundaries

We analyze the hydrodynamic stability of force-driven parallel shear flows in nonequilibrium molecular simulations with three-dimensional periodic boundary conditions. We show that flows simulated in this way can be linearly unstable, and we derive an expression for the critical Reynolds number as a function of the geometric aspect ratio of the simulation domain. Approximate periodic extensions of Couette and Poiseuille flows are unstable at Reynolds numbers two orders of magnitude smaller than their aperiodic equivalents because the periodic boundaries impose fundamentally different constraints on the flow. This instability has important implications for simulating shear rheology and for designing nonequilibrium simulation methods that are compatible with periodic boundary conditions.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Universal Gelation of Metal Oxide Nanocrystals via Depletion Attractions

Nanocrystal gelation provides a powerful framework to translate nanoscale properties into bulk materials and to engineer emergent properties through the assembled microstructure. However, many established gelation strategies rely on chemical reactions and specific interactions, e.g., stabilizing ligands or ions on the surface of the nanocrystals, and are therefore not easily transferrable. Here, we report a general gelation strategy via non-specific and purely entropic depletion attractions applied to three types of metal oxide nanocrystals. The gelation thresholds of two compositionally distinct spherical nanocrystals agree quantitatively, demonstrating the adaptability of the approach for different chemistries. Consistent with theoretical phase behavior predictions, nanocrystal cubes form gels at a lower polymer concentration than nanocrystal spheres, allowing shape to serve as a handle to control gelation. These results suggest that the fundamental underpinnings of depletion-driven assembly, traditionally associated with larger colloidal particles, are also applicable at the nanoscale.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Wertheim's thermodynamic perturbation theory with double-bond association and its application to colloid-linker mixtures

We extend Wertheim's thermodynamic perturbation theory to derive the association free energy of a multicomponent mixture for which double bonds can form between any two pairs of the molecules' arbitrary number of bonding sites. This generalization reduces in limiting cases to prior theories that restrict double bonding to at most one pair of sites per molecule. We apply the new theory to an associating mixture of colloidal particles ("colloids") and flexible chain molecules ("linkers"). The linkers have two functional end groups, each of which may bond to one of several sites on the colloids. Due to their flexibility, a significant fraction of linkers can "loop" with both ends bonding to sites on the same colloid instead of bridging sites on different colloids. We use the theory to show that the fraction of linkers in loops depends sensitively on the linker end-to-end distance relative to the colloid bonding-site distance, which suggests strategies for mitigating the loop formation that may otherwise hinder linker-mediated colloidal assembly.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Assembly of nothing: Equilibrium fluids with designed structured porosity

Controlled micro- to meso-scale porosity is a common materials design goal with possible applications ranging from molecular gas adsorption to particle size selective permeability or solubility. Here, we use inverse methods of statistical mechanics to design an isotropic pair interaction that, in the absence of an external field, assembles particles into an inhomogeneous {\em fluid} matrix surrounding pores of prescribed size ordered in a lattice morphology. The pore size can be tuned via modification of temperature or particle concentration. Moreover, modulating density reveals a rich series of microphase-separated morphologies including pore- or particle-based lattices, pore- or particle-based columns, and bicontinuous or lamellar structures. Sensitivity of pore assembly to the form of the designed interaction potential is explored.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Consequences of minimizing pair correlations in fluids for dynamics, thermodynamics, and structure

Liquid-state theory, computer simulation, and numerical optimization are used to investigate the extent to which positional correlations of a hard-sphere fluid--as characterized by the radial distribution function and the two-particle excess entropy--can be suppressed via the introduction of auxiliary pair interactions. The corresponding effects of such interactions on total excess entropy, density fluctuations, and single-particle dynamics are explored. Iso-g processes, whereby hard-sphere-fluid pair structure at a given density is preserved at higher densities via the introduction of a density-dependent, soft repulsive contribution to the pair potential, are considered. Such processes eventually terminate at a singular density, resulting in a state that--while incompressible and hyperuniform--remains unjammed and exhibits fluid-like dynamic properties. The extent to which static pair correlations can be suppressed to maximize pair disorder in a fluid with hard cores, determined via direct functional maximization of two-body excess entropy, is also considered. Systems approaching a state of maximized two-body entropy display a progressively growing bandwidth of suppressed density fluctuations, pointing to a relation between "stealthiness" and maximal pair disorder in materials.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Fluids with competing interactions: I. Decoding the structure factor to detect and characterize self-limited clustering

We use liquid state theory and computer simulations to gain insights into the shape of the structure factor for fluids of particles interacting via a combination of short-range attractions and long-range repulsions. Such systems can reversibly morph between homogeneous phases and states comprising compact self-limiting clusters. We first highlight trends with respect to the presence and location of the intermediate-range order (IRO) pre-peak in the structure factor, which is commonly associated with clustering, for wide ranges of the tunable parameters that control interparticle interactions (e.g., Debye screening length). Next, for approximately 100 different cluster phases at various conditions (where aggregates range in size from six to sixty monomers), we quantitatively relate the shape of the structure factor to physical characteristics including intercluster distance and cluster size. We also test two previously postulated criteria for identifying the emergence of clustered phases that are based on IRO peak-height and -width, respectively. We find that the criterion based on peak-width, which encodes the IRO thermal correlation length, is more robust across a wide range of conditions and interaction strengths but nonetheless approximate. Ultimately, we recommend a hybrid heuristic drawing on both pre-peak height and width for positively identifying the emergence of clustered states.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Fluids with competing interactions: II. Validating a free energy model for equilibrium cluster size

Using computer simulations, we validate a simple free energy model that can be analytically solved to predict the equilibrium size of self-limiting clusters of particles in the fluid state governed by a combination of short-range attractive and long-range repulsive pair potentials. The model is a semi-empirical adaptation and extension of the canonical free energy-based result due to Groenewold and Kegel [J. Phys. Chem. B, 105 (2001)], where we use new computer simulation data to systematically improve the cluster-size scalings with respect to the strengths of the competing interactions driving aggregation. We find that one can adapt a classical nucleation like theory for small energetically-frustrated aggregates provided one appropriately accounts for a size-dependent, microscopic energy penalty of interface formation, which requires new scaling arguments. This framework is verified in part by considering the extensive scaling of intracluster bonding, where we uncover a superlinear scaling regime distinct from (and located between) the known regimes for small and large aggregates. We validate our model based on comparisons against approximately 100 different simulated systems comprising compact spherical aggregates with characteristic (terminal) sizes between six and sixty monomers, which correspond to wide ranges in experimentally-controllable parameters.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Graphoepitaxy for translational and orientational ordering of monolayers of rectangular nanoparticles

The combinations of particle aspect ratio and enthalpic-barrier templates that lead to translational and orientational ordering of monolayers of rectangular particles are determined using Monte Carlo simulations and density functional theory. For sufficiently high enthalpic barriers, we find that only specific combinations of particle sizes and template spacings lead to ordered arrays. The pattern multiplication factor provided by the template extends to approximately ten times the smallest dimension of the particle.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Inverse Design for Self Assembly via On-the-Fly Optimization

Inverse methods of statistical mechanics have facilitated the discovery of pair potentials that stabilize a wide variety of targeted lattices at zero temperature. However, such methods are complicated by the need to compare, within the optimization framework, the energy of the desired lattice to all possibly relevant competing structures, which are not generally known in advance. Furthermore, ground-state stability does not guarantee that the target will readily assemble from the fluid upon cooling from higher temperature. Here, we introduce a molecular dynamics simulation-based, optimization design strategy that iteratively and systematically refines the pair interaction according to the fluid and crystalline structural ensembles encountered during the assembly process. We successfully apply this probabilistic, machine-learning approach to the design of repulsive, isotropic pair potentials that assemble into honeycomb, kagome, square, rectangular, truncated square and truncated hexagonal lattices.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

On the Formation of Equilibrium Gels via a Macroscopic Bond Limitation

Restricting the number of attractive physical "bonds" that can form between particles in a fluid suppresses the usual demixing phase transition to very low particle concentrations, allowing for the formation of open, percolated, and homogeneous states, aptly called equilibrium or "empty" gels. Most demonstrations of this concept have directly limited the microscopic particle valence via anisotropic (patchy) attractions; however, an alternative macroscopic valence limitation would be desirable for greater experimental tunability and responsiveness. One possibility, explored in this paper, is to employ primary particles with attractions mediated via a secondary species of linking particles. In such a system, the linker-to-primary particle ratio serves as a macroscopic control parameter for the average microscopic valence. We show that the phase behavior of such a system predicted by Wertheim's first order perturbation theory is consistent with equilibrium gel formation: the primary particle concentrations corresponding to the two-phase demixing transition are significantly suppressed at both low and high linker-to-primary particle ratios. Extensive molecular dynamics simulations validate these theoretical predictions but also reveal the presence of loops of bonded particles, which are neglected in the theory. Such loops cause densification and inhibit percolation, and hence the range of viable empty gel state conditions is somewhat reduced relative to the Wertheim theory predictions.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Equilibrium cluster fluids: Pair interactions via inverse design

Inverse methods of statistical mechanics are becoming productive tools in the design of materials with specific microstructures or properties. While initial studies have focused on solid-state design targets (e.g, assembly of colloidal superlattices), one can alternatively design fluid states with desired morphologies. This work addresses the latter and demonstrates how a simple iterative Boltzmann inversion strategy can be used to determine the isotropic pair potential that reproduces the radial distribution function of a fluid of amorphous clusters with prescribed size. The inverse designed pair potential of this "ideal" cluster fluid, with its broad attractive well and narrow repulsive barrier at larger separations, is qualitatively different from the so-called SALR form most commonly associated with equilibrium cluster formation in colloids, which features short-range attractive (SA) and long-range repulsive (LR) contributions. These differences reflect alternative mechanisms for promoting cluster formation with an isotropic pair potential, and they in turn produce structured fluids with qualitatively different static and dynamic properties. Specifically, equilibrium simulations show that the amorphous clusters resulting from the inverse designed potentials display more uniformity in size and shape, and they also show greater spatial and temporal resolution than those resulting from SALR interactions.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Local structure-mobility relationships of confined fluids reverse upon supercooling

We examine the structural and dynamic properties of confined binary hard-sphere mixtures designed to mimic realizable colloidal thin films. Using computer simulations, governed by either Newtonian or overdamped Langevin dynamics, together with other techniques including a Fokker-Planck equation-based method, we measure the position-dependent and average diffusivities of particles along structurally isotropic and inhomogeneous dimensions of the fluids. At moderate packing fractions, local single-particle diffusivities normal to the direction of confinement are higher in regions of high total packing fraction; however, these trends are reversed as the film is supercooled at denser average packings. Auxiliary short-time measurements of particle displacements mirror data obtained for experimental supercooled colloidal systems. We find that average dynamics can be approximately predicted based on the distribution of available space for particle insertion across orders of magnitude in diffusivity regardless of the governing microscopic dynamics.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Origin and Detection of Microstructural Clustering in Fluids with Spatial-Range Competitive Interactions

Fluids with competing short-range attractions and long-range repulsions mimic dispersions of charge-stabilized colloids that can display equilibrium structures with intermediate range order (IRO), including particle clusters. Using simulations and analytical theory, we demonstrate how to detect cluster formation in such systems from the static structure factor and elucidate links to macrophase separation in purely attractive reference fluids. We find that clusters emerge when the thermal correlation length encoded in the IRO peak of the structure factor exceeds the characteristic lengthscale of interparticle repulsions. We also identify qualitative differences between the dynamics of systems that form amorphous versus micro-crystalline clusters.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Tuning structure and mobility of solvation shells surrounding tracer additives

Molecular dynamics simulations and a stochastic Fokker-Planck equation based approach are used to illuminate how position-dependent solvent mobility near one or more tracer particle(s) is affected when tracer-solvent interactions are rationally modified to affect corresponding solvation structure. For tracers in a dense hard-sphere fluid, we compare two types of tracer-solvent interactions: (1) a hard-sphere-like interaction; and (2) a soft repulsion extending beyond the hard core designed via statistical mechanical theory to enhance tracer mobility at infinite dilution by suppressing coordination-shell structure (Carmer et al., Soft Matter 8 (2011)). For the latter case, we show that the mobility of surrounding solvent particles is also increased by addition of the soft repulsive interaction, which helps to rationalize the mechanism underlying the tracer's enhanced diffusivity. However, if multiple tracer surfaces are in closer proximity (as at higher tracer concentrations), similar interactions that disrupt local solvation structure instead suppress the position-dependent solvent dynamics.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Dimensionality and design of isotropic interactions that stabilize honeycomb, square, simple cubic, and diamond lattices

We use inverse methods of statistical mechanics and computer simulations to investigate whether an isotropic interaction designed to stabilize a given two-dimensional (2D) lattice will also favor an analogous three-dimensional (3D) structure, and vice versa. Specifically, we determine the 3D ordered lattices favored by isotropic potentials optimized to exhibit stable 2D honeycomb (or square) periodic structures, as well as the 2D ordered structures favored by isotropic interactions designed to stabilize 3D diamond (or simple cubic) lattices. We find a remarkable `transferability' of isotropic potentials designed to stabilize analogous morphologies in 2D and 3D, irrespective of the exact interaction form, and we discuss the basis of this cross-dimensional behavior. Our results suggest that the discovery of interactions that drive assembly into certain 3D periodic structures of interest can be assisted by less computationally intensive optimizations targeting the analogous 2D lattices.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Graphoepitaxy for Pattern Multiplication of Nanoparticle Monolayers

We compute the free energy minimizing structures of particle monolayers in the presence of enthalpic barriers of a finite height \b{eta}Vext using classical density functional theory and Monte Carlo simulations. We show that a periodic square template with dimensions up to at least ten times the particle diameter disrupts the formation of the entropically favored hexagonally close-packed 2D lattice in favor of a square lattice. The results illustrate how graphoepitaxy can successfully order nanoparticulate films into desired patterns many times smaller than those of the prepatterned template.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Position-dependent and pair diffusivity profiles from steady-state solutions of color reaction-counterdiffusion problems

We show how the steady-state solution of the Smoluchowski (Fokker-Planck) equation for a color reaction-counterdiffusion problem, together with equilibrium trajectory information (e.g., from molecular simulations or confocal microscopy experiments), can be used to determine position-dependent or pair diffusivity profiles for condensed phases directly without the need for Bayesian inference, likelihood maximization, or the analysis of time-dependent data.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Predicting the structure of fluids with piecewise constant interactions: Comparing the accuracy of five efficient integral equation theories

We use molecular dynamics simulations to test integral equation theory predictions for the structure of fluids of spherical particles with eight different piecewise-constant pair interaction forms comprising a hard core and a combination of two shoulders and/or wells. Since model pair potentials like these are of interest for discretized or coarse-grained representations of effective interactions in complex fluids (e.g., for computationally intensive inverse optimization problems), we focus here on assessing how accurately their properties can be predicted by analytical or simple numerical closures including Percus-Yevick, hypernetted chain, reference hypernetted chain, first-order mean spherical approximation, and a modified first-order mean spherical approximation. To make quantitative comparisons between the predicted and simulated radial distribution functions, we introduce a cumulative structural error metric. For equilibrium fluid state points of these models, we find that the reference hypernetted chain closure is the most accurate of the tested approximations as characterized by this metric or related thermodynamic quantities.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

User-friendly, theory-based web applet for rapidly predicting structure and thermodynamics of complex fluids

Based on a recently introduced analytical strategy [Hollingshead et al., J. Chem. Phys. 139, 161102 (2013)], we present a web applet that can quickly and semi-quantitatively estimate the equilibrium radial distribution function and related thermodynamic properties of a fluid from knowledge of its pair interaction. We present a detailed description of the applet's features and intended workflow, followed by a description of how the applet can be used to illustrate two (of many possible) concepts of interest for introductory statistical mechanics courses: the transition from ideal gas-like behavior to correlated-liquid behavior with increasing density and the tradeoff between dominant length scales with changing temperature in a system with ramp-shaped repulsions. The latter type of interaction qualitatively captures distinctive thermodynamic properties of liquid water because its energetic bias toward locally open structures mimics that of water's hydrogen-bond network.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Fine discretization of pair interactions and an approximate analytical strategy for predicting equilibrium behavior of complex fluids

We study whether fine discretization (i.e., terracing) of continuous pair interactions, when used in combination with first-order mean-spherical approximation theory, can lead to a simple and general analytical strategy for predicting the equilibrium structure and thermodynamics of complex fluids. Specifically, we implement a version of this approach to predict how screened electrostatic repulsions, solute-mediated depletion attractions, or ramp-shaped repulsions modify the radial distribution function and the potential energy of reference hard-sphere fluids, and we compare the predictions to exact results from molecular simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Phase Behavior of Materials with Isotropic Interactions Designed by Inverse Strategies to Favor Diamond and Simple Cubic Lattice Ground States

We use molecular simulation to construct equilibrium phase diagrams for two recently introduced model materials with isotropic, soft-repulsive pair interactions designed to favor diamond and simple cubic lattice ground states, respectively, over a wide range of densities [Jain et al., Soft Matter 9 14 (2013)]. We employ free energy based Monte Carlo simulation techniques to precisely trace the inter-crystal and fluid-crystal coexistence curves. We find that both model materials display rich polymorphic phase behavior featuring stable crystals corresponding to the target ground-state structures, as well as a variety of other crystalline (e.g., hexagonal and body-centered cubic) phases and multiple reentrant melting transitions.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Predicting how nanoconfinement changes the relaxation time of a supercooled liquid

The properties of nanoconfined fluids can be strikingly different from those of bulk liquids. A basic unanswered question is whether the equilibrium and dynamic consequences of confinement are related to each other in a simple way. We study this question by simulation of a liquid comprising asymmetric dumbbell-shaped molecules, which can be deeply supercooled without crystallizing. We find that the dimensionless structural relaxation times - spanning six decades as a function of temperature, density, and degree of confinement - collapse when plotted versus excess entropy. The data also collapse when plotted versus excess isochoric heat capacity, a behaviour that follows from the existence of isomorphs in the bulk and confined states.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Enhancing tracer diffusivity by tuning interparticle interactions and coordination shell structure

This study uses a combination of stochastic optimization, statistical mechanical theory, and molecular simulation to test the extent to which the long-time dynamics of a single tracer particle can be enhanced by rationally modifying its interactions--and hence static correlations--with the other particles of a dense fluid. Specifically, a simulated annealing strategy is introduced that, when coupled with test-particle calculations from an accurate density functional theory, finds interactions that maximize either the tracer's partial molar excess entropy or a related pair-correlation measure (i.e., two quantities known to correlate with tracer diffusivity in other contexts). The optimized interactions have soft, Yukawa-like repulsions, which extend beyond the hard-sphere interaction and disrupt the coordination-shell cage structure surrounding the tracer. Molecular and Brownian dynamics simulations find that tracers with these additional soft repulsions can diffuse more than three times faster than bare hard spheres in a moderately supercooled fluid, despite the fact that the former appear considerably larger than the latter by conventional definitions of particle size.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Generalizing Rosenfeld's excess-entropy scaling to predict long-time diffusivity in dense fluids of Brownian particles: From hard to ultrasoft interactions

Computer simulations are used to test whether a recently introduced generalization of Rosenfeld's excess-entropy scaling method for estimating transport coefficients in systems obeying molecular dynamics can be extended to predict long-time diffusivities in fluids of particles undergoing Brownian dynamics in the absence of interparticle hydrodynamic forces. Model fluids with inverse-power-law, Gaussian-core, and Hertzian pair interactions are considered. Within the generalized Rosenfeld scaling method, long-time diffusivities of ultrasoft Gaussian-core and Hertzian particle fluids, which display anomalous trends with increasing density, are predicted (to within 20%) based on knowledge of interparticle interactions, excess entropy, and scaling behavior of simpler inverse-power-law fluids.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Impact of surface roughness on diffusion of confined fluids

Using event-driven molecular dynamics simulations, we quantify how the self diffusivity of confined hard-sphere fluids depends on the nature of the confining boundaries. We explore systems with featureless confining boundaries that treat particle-boundary collisions in different ways and also various types of physically (i.e., geometrically) rough boundaries. We show that, for moderately dense fluids, the ratio of the self diffusivity of a rough wall system to that of an appropriate smooth-wall reference system is a linear function of the reciprocal wall separation, with the slope depending on the nature of the roughness. We also discuss some simple practical ways to use this information to predict confined hard-sphere fluid behavior in different rough-wall systems.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Implications of the effective one-component analysis of pair correlations in colloidal fluids with polydispersity

Partial pair-correlation functions of colloidal suspensions with continuous polydispersity can be challenging to characterize from optical microscopy or computer simulation data due to inadequate sampling. As a result, it is common to adopt an effective one-component description of the structure that ignores the differences between particle types. Unfortunately, whether this kind of simplified description preserves or averages out information important for understanding the behavior of the fluid depends on the degree of polydispersity and can be difficult to assess, especially when the corresponding multicomponent description of the pair correlations is unavailable for comparison. Here, we present a computer simulation study that examines the implications of adopting an effective one-component structural description of a polydisperse fluid. The square-well model that we investigate mimics key aspects of the experimental behavior of suspended colloids with short-range, polymer-mediated attractions. To characterize the partial pair-correlation functions and thermodynamic excess entropy of this system, we introduce a Monte Carlo sampling strategy appropriate for fluids with a large number of pseudo-components. The data from our simulations at high particle concentrations, as well as exact theoretical results for dilute systems, show how qualitatively different trends between structural order and particle attractions emerge from the multicomponent and effective one-component treatments, even with systems characterized by moderate polydispersity. We examine consequences of these differences for excess-entropy based scalings of shear viscosity, and we discuss how use of the multicomponent treatment reveals similarities between the corresponding dynamic scaling behaviors of attractive colloids and liquid water that the effective one-component analysis does not capture.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Mapping between long-time molecular and Brownian dynamics

We use computer simulations to test a simple idea for mapping between long-time self diffusivities obtained from molecular and Brownian dynamics. The strategy we explore is motivated by the behavior of fluids comprising particles that interact via inverse-power-law pair potentials, which serve as good reference models for dense atomic or colloidal materials. Based on our simulation data, we present an empirical expression that semi-quantitatively describes the "atomic" to "colloidal" diffusivity mapping for inverse-power-law fluids, but also for model complex fluids with considerably softer (star-polymer, Gaussian-core, or Hertzian) interactions. As we show, the anomalous structural and dynamic properties of these latter ultrasoft systems pose problems for other strategies designed to relate Newtonian and Brownian dynamics of hard-sphere-like particles.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Structural signatures of mobility on intermediate time scales in a supercooled fluid

We use computer simulations to explore the manner in which the particle displacements on intermediate time scales in supercooled fluids correlate to their dynamic structural environment. The fluid we study, a binary mixture of hard spheres, exhibits classic signatures of dynamic heterogeneity, including a bifurcated single-particle displacement distribution (i.e., subpopulations of immobile and mobile particles). We find that immobile particles, during the course of their displacements, exhibit stronger average pair correlations to their neighbors than mobile particles, but not necessarily higher average coordination numbers. We discuss how the correlation between structure and single-particle dynamics depends on observation time.

0 citations0 reviewsNo code linkedOpen access