-
Effects of Ligand vs. Linker on Phase Behavior and Mechanical Properties of Nanoparticle Gels
Authors:
Qizan Chen,
Dinesh Sundaravadivelu Devarajan,
Arash Nikoubashman,
Michael P. Howard,
Jeetain Mittal
Abstract:
Nanoparticle gels have attracted considerable attention due to their highly tunable properties. One strategy for producing nanoparticle gels involves utilizing strong local attractions between polymeric molecules, such as DNA hybridization or dynamic covalent chemistry, to form percolated nanoparticle networks. These molecules can be used in two distinct roles: as ``ligands'' with one end grafted…
▽ More
Nanoparticle gels have attracted considerable attention due to their highly tunable properties. One strategy for producing nanoparticle gels involves utilizing strong local attractions between polymeric molecules, such as DNA hybridization or dynamic covalent chemistry, to form percolated nanoparticle networks. These molecules can be used in two distinct roles: as ``ligands'' with one end grafted to a nanoparticle or as ``linkers'' with both ends free. Here, we explore how these roles shape the phase behavior and mechanical properties of gel-like nanoparticle assemblies using coarse-grained simulations. We systematically vary the interaction strength and bending stiffness of both ligands and linkers. We find that phase separation can be limited to low nanoparticle volume fractions by making the ligands rigid, consistent with previous studies on linked nanoparticle gels. At fixed interaction strength and volume fraction, both ligand- and linker-mediated nanoparticle assemblies show similar mechanical responses to variations in bending stiffness. However, a comparison between the two association schemes reveals that the linked nanoparticles form rigid percolated networks that are less stretchable than the ligand-grafted gels, despite exhibiting similar tensile strength. We attribute these differences between ligands and linkers to the distinct structural arrangement of nanoparticles. Our findings highlight the potential to utilize different association schemes to tune specific mechanical properties.
△ Less
Submitted 18 October, 2024;
originally announced October 2024.
-
relentless: Transparent, reproducible molecular dynamics simulations for optimization
Authors:
Adithya N Sreenivasan,
C. Levi Petix,
Zachary M. Sherman,
Michael P. Howard
Abstract:
relentless is an open-source Python package that enables the optimization of objective functions computed using molecular dynamics simulations. It has a high-level, extensible interface for model parametrization; setting up, running, and analyzing simulations natively in established software packages; and gradient-based optimization. We describe the design and implementation of relentless in the c…
▽ More
relentless is an open-source Python package that enables the optimization of objective functions computed using molecular dynamics simulations. It has a high-level, extensible interface for model parametrization; setting up, running, and analyzing simulations natively in established software packages; and gradient-based optimization. We describe the design and implementation of relentless in the context of relative entropy minimization, and we demonstrate its abilities to design pairwise interactions between particles that form targeted structures. relentless aims to streamline the development of computational materials design methodologies and promote the transparency and reproducibility of complex workflows integrating molecular dynamics simulations.
△ Less
Submitted 6 August, 2024;
originally announced August 2024.
-
Mesoscale simulations of diffusion and sedimentation in shape-anisotropic nanoparticle suspensions
Authors:
Yashraj M. Wani,
Penelope Grace Kovakas,
Arash Nikoubashman,
Michael P. Howard
Abstract:
We determine the long-time self-diffusion coefficient and sedimentation coefficient for suspensions of nanoparticles with anisotropic shapes (octahedra, cubes, tetrahedra, and spherocylinders) as a function of nanoparticle concentration using mesoscale simulations. We use a discrete particle model for the nanoparticles, and we account for solvent-mediated hydrodynamic interactions between nanopart…
▽ More
We determine the long-time self-diffusion coefficient and sedimentation coefficient for suspensions of nanoparticles with anisotropic shapes (octahedra, cubes, tetrahedra, and spherocylinders) as a function of nanoparticle concentration using mesoscale simulations. We use a discrete particle model for the nanoparticles, and we account for solvent-mediated hydrodynamic interactions between nanoparticles using the multiparticle collision dynamics method. Our simulations are compared to theoretical predictions and experimental data from existing literature, demonstrating good agreement in the majority of cases. Further, we find that the self-diffusion coefficient of the regular polyhedral shapes can be estimated from that of a sphere whose diameter is average of their inscribed and circumscribed sphere diameters.
△ Less
Submitted 5 March, 2024;
originally announced March 2024.
-
Multiscale modeling of solute diffusion in triblock copolymer membranes
Authors:
Anthony J. Cooper,
Michael P. Howard,
Sanket Kadulkar,
David Zhao,
Kris T. Delaney,
Venkat Ganesan,
Thomas M. Truskett,
Glenn H. Fredrickson
Abstract:
We develop a multiscale simulation model for diffusion of solutes through porous triblock copolymer membranes. The approach combines two techniques: self-consistent field theory (SCFT) to predict the structure of the self-assembled, solvated membrane and on-lattice kinetic Monte Carlo (kMC) simulations to model diffusion of solutes. Solvation is simulated in SCFT by constraining the glassy membran…
▽ More
We develop a multiscale simulation model for diffusion of solutes through porous triblock copolymer membranes. The approach combines two techniques: self-consistent field theory (SCFT) to predict the structure of the self-assembled, solvated membrane and on-lattice kinetic Monte Carlo (kMC) simulations to model diffusion of solutes. Solvation is simulated in SCFT by constraining the glassy membrane matrix while relaxing the brush-like membrane pore coating against the solvent. The kMC simulations capture the resulting solute spatial distribution and concentration-dependent local diffusivity in the polymer-coated pores; we parameterize the latter using particle-based simulations. We apply our approach to simulate solute diffusion through nonequilibrium morphologies of a model triblock copolymer, and we correlate diffusivity with structural descriptors of the morphologies. We also compare the model's predictions to alternative approaches based on simple lattice random walks and find our multiscale model to be more robust and systematic to parameterize. Our multiscale modeling approach is general and can be readily extended in the future to other chemistries, morphologies, and models for the local solute diffusivity and interactions with the membrane.
△ Less
Submitted 20 September, 2022;
originally announced September 2022.
-
Dynamic density functional theory for drying colloidal suspensions: Comparison of hard-sphere free-energy functionals
Authors:
Mayukh Kundu,
Michael P. Howard
Abstract:
Dynamic density functional theory (DDFT) is a promising approach for predicting the structural evolution of a drying suspension containing one or more types of colloidal particles. The assumed free-energy functional is a key component of DDFT that dictates the thermodynamics of the model and, in turn, the density flux due to a concentration gradient. In this work, we compare several commonly used…
▽ More
Dynamic density functional theory (DDFT) is a promising approach for predicting the structural evolution of a drying suspension containing one or more types of colloidal particles. The assumed free-energy functional is a key component of DDFT that dictates the thermodynamics of the model and, in turn, the density flux due to a concentration gradient. In this work, we compare several commonly used free-energy functionals for drying hard-sphere suspensions including local-density approximations based on the ideal-gas, virial, and Boublík-Mansoori-Carnahan-Starling-Leland (BMCSL) equations of state as well as a weighted-density approximation based on fundamental measure theory (FMT). To determine the accuracy of each functional, we model one- and two-component hard-sphere suspensions in a drying film with varied initial heights and compositions, and we compare the DDFT-predicted volume-fraction profiles to particle-based Brownian dynamics (BD) simulations. FMT accurately predicts the structure of the one-component suspensions even at high concentrations and when significant density gradients develop, but the virial and BMCSL equations of state provide reasonable approximations for smaller concentrations at a reduced computational cost. In the two-component suspensions, FMT and BMCSL are similar to each other but modestly overpredict the extent of stratification by size compared to BD simulations. This work provides helpful guidance for selecting thermodynamic models for soft materials in nonequilibrium processes such as solvent drying, solvent freezing, and sedimentation.
△ Less
Submitted 2 August, 2022;
originally announced August 2022.
-
Diffusion and sedimentation in colloidal suspensions using multiparticle collision dynamics with a discrete particle model
Authors:
Yashraj M. Wani,
Penelope Grace Kovakas,
Arash Nikoubashman,
Michael P. Howard
Abstract:
We study self-diffusion and sedimentation in colloidal suspensions of nearly-hard spheres using the multiparticle collision dynamics simulation method for the solvent with a discrete mesh model for the colloidal particles (MD+MPCD). We cover colloid volume fractions from 0.01 to 0.40 and compare the MD+MPCD simulations to Brownian dynamics simulations with free-draining hydrodynamics (BD) as well…
▽ More
We study self-diffusion and sedimentation in colloidal suspensions of nearly-hard spheres using the multiparticle collision dynamics simulation method for the solvent with a discrete mesh model for the colloidal particles (MD+MPCD). We cover colloid volume fractions from 0.01 to 0.40 and compare the MD+MPCD simulations to Brownian dynamics simulations with free-draining hydrodynamics (BD) as well as pairwise far-field hydrodynamics described using the Rotne--Prager--Yamakawa mobility tensor (BD+RPY). The dynamics in MD+MPCD suggest that the colloidal particles are only partially coupled to the solvent at short times. However, the long-time self-diffusion coefficient in MD+MPCD is comparable to that in BD and BD+RPY, and the sedimentation coefficient in MD+MPCD is in good agreement with that in BD+RPY, suggesting that MD+MPCD gives a reasonable description of the hydrodynamic interactions in colloidal suspensions. The discrete-particle MD+MPCD approach is convenient and readily extended to more complex shapes, and we determine the long-time self-diffusion coefficient in suspensions of nearly-hard cubes to demonstrate its generality.
△ Less
Submitted 12 October, 2021;
originally announced October 2021.
-
Effects of linker flexibility on phase behavior and structure of linked colloidal gels
Authors:
Michael P. Howard,
Zachary M. Sherman,
Adithya N Sreenivasan,
Stephanie A. Valenzuela,
Eric V. Anslyn,
Delia J. Milliron,
Thomas M. Truskett
Abstract:
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…
▽ More
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.
△ Less
Submitted 24 November, 2020;
originally announced November 2020.
-
Wertheim's thermodynamic perturbation theory with double-bond association and its application to colloid-linker mixtures
Authors:
Michael P. Howard,
Zachary M. Sherman,
Delia J. Milliron,
Thomas M. Truskett
Abstract:
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 associ…
▽ More
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.
△ Less
Submitted 10 December, 2020; v1 submitted 15 October, 2020;
originally announced October 2020.
-
Stratification of polymer mixtures in drying droplets: hydrodynamics and diffusion
Authors:
Michael P. Howard,
Arash Nikoubashman
Abstract:
We study the evaporation-induced stratification of a mixture of short and long polymer chains in a drying droplet using molecular simulations. We systematically investigate the effects of hydrodynamic interactions (HI) on this process by comparing hybrid simulations accounting for HI between polymers through the multiparticle collision dynamics technique with free-draining Langevin dynamics simula…
▽ More
We study the evaporation-induced stratification of a mixture of short and long polymer chains in a drying droplet using molecular simulations. We systematically investigate the effects of hydrodynamic interactions (HI) on this process by comparing hybrid simulations accounting for HI between polymers through the multiparticle collision dynamics technique with free-draining Langevin dynamics simulations neglecting the same. We find that the dried supraparticle morphologies are homogeneous when HI are included but are stratified in core--shell structures (with the short polymers forming the shell) when HI are neglected. The simulation methodology unambiguously attributes this difference to the treatment of the solvent in the two models. We rationalize the presence (or absence) of stratification by measuring phenomenological multicomponent diffusion coefficients for the polymer mixtures. The diffusion coefficients show the importance of not only solvent backflow but also HI between polymers in controlling the dried supraparticle morphology.
△ Less
Submitted 16 July, 2020; v1 submitted 19 May, 2020;
originally announced May 2020.
-
Inverse methods for design of soft materials
Authors:
Zachary M. Sherman,
Michael P. Howard,
Beth A. Lindquist,
Ryan B. Jadrich,
Thomas M. Truskett
Abstract:
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…
▽ More
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.
△ Less
Submitted 31 March, 2020;
originally announced April 2020.
-
Universal Gelation of Metal Oxide Nanocrystals via Depletion Attractions
Authors:
Camila A. Saez Cabezas,
Zachary M. Sherman,
Michael P. Howard,
Manuel N. Dominguez,
Shin Hum Cho,
Gary K. Ong,
Allison Green,
Thomas M. Truskett,
Delia J. Milliron
Abstract:
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, w…
▽ More
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.
△ Less
Submitted 25 March, 2020;
originally announced March 2020.
-
Stability of force-driven shear flows in nonequilibrium molecular simulations with periodic boundaries
Authors:
Michael P. Howard,
Antonia Statt,
Howard A. Stone,
Thomas M. Truskett
Abstract:
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…
▽ More
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.
△ Less
Submitted 13 May, 2020; v1 submitted 16 July, 2019;
originally announced July 2019.
-
Structure and phase behavior of polymer-linked colloidal gels
Authors:
Michael P. Howard,
Ryan B. Jadrich,
Beth A. Lindquist,
Fardin Khabaz,
Roger T. Bonnecaze,
Delia J. Milliron,
Thomas M. Truskett
Abstract:
Low-density "equilibrium" gels that consist of a percolated, kinetically arrested network of colloidal particles and are resilient to aging can be fabricated by restricting the number of effective bonds that form between the colloids. Valence-restricted patchy particles have long served as one archetypal example of such materials, but equilibrium gels can also be realized through a synthetically s…
▽ More
Low-density "equilibrium" gels that consist of a percolated, kinetically arrested network of colloidal particles and are resilient to aging can be fabricated by restricting the number of effective bonds that form between the colloids. Valence-restricted patchy particles have long served as one archetypal example of such materials, but equilibrium gels can also be realized through a synthetically simpler and scalable strategy that introduces a secondary linker, such as a small ditopic molecule, to mediate the bonds between the colloids. Here, we consider the case where the ditopic linker molecules are low-molecular-weight polymers and demonstrate using a model colloid-polymer mixture how macroscopic properties such as the phase behavior as well as the microstructure of the gel can be designed through the polymer molecular weight and concentration. The low-density window for equilibrium gel formation is favorably expanded using longer linkers, while necessarily increasing the spacing between all colloids. However, we show that blends of linkers with different sizes enable wider variation in microstructure for a given target phase behavior. Our computational study suggests a robust and tunable strategy for the experimental realization of equilibrium colloidal gels.
△ Less
Submitted 10 July, 2019;
originally announced July 2019.
-
The Role of Pressure in Inverse Design for Assembly
Authors:
Beth A. Lindquist,
Ryan B. Jadrich,
Michael P. Howard,
Thomas M. Truskett
Abstract:
Isotropic pairwise interactions that promote the self assembly of complex particle morphologies have been discovered by inverse design strategies derived from the molecular coarse-graining literature. While such approaches provide an avenue to reproduce structural correlations, thermodynamic quantities such as the pressure have typically not been considered in self-assembly applications. In this w…
▽ More
Isotropic pairwise interactions that promote the self assembly of complex particle morphologies have been discovered by inverse design strategies derived from the molecular coarse-graining literature. While such approaches provide an avenue to reproduce structural correlations, thermodynamic quantities such as the pressure have typically not been considered in self-assembly applications. In this work, we demonstrate that relative entropy optimization can be used to discover potentials that self-assemble into targeted cluster morphologies with a prescribed pressure when the iterative simulations are performed in the isothermal-isobaric ensemble. By tuning the pressure in the optimization, we generate a family of simple pair potentials that all self-assemble the same structure. Selecting an appropriate simulation ensemble to control the thermodynamic properties of interest is a general design strategy that could also be used to discover interaction potentials that self-assemble structures having, for example, a specified chemical potential.
△ Less
Submitted 4 June, 2019;
originally announced June 2019.
-
Crack formation and self-closing in shrinkable, granular packings
Authors:
H. Jeremy Cho,
Nancy B. Lu,
Michael P. Howard,
Rebekah A. Adams,
Sujit S. Datta
Abstract:
Many clays, soils, biological tissues, foods, and coatings are shrinkable, granular materials: they are composed of packed, hydrated grains that shrink when dried. In many cases, these packings crack during drying, critically hindering applications. However, while cracking has been widely studied for bulk gels and packings of non-shrinkable grains, little is known about how packings of shrinkable…
▽ More
Many clays, soils, biological tissues, foods, and coatings are shrinkable, granular materials: they are composed of packed, hydrated grains that shrink when dried. In many cases, these packings crack during drying, critically hindering applications. However, while cracking has been widely studied for bulk gels and packings of non-shrinkable grains, little is known about how packings of shrinkable grains crack. Here, we elucidate how grain shrinkage alters cracking during drying. Using experiments with model shrinkable hydrogel beads, we show that differential shrinkage can dramatically alter crack evolution during drying---in some cases, even causing cracks to spontaneously "self-close". In other cases, packings shrink without cracking or crack irreversibly. We developed both granular and continuum models to quantify the interplay between grain shrinkage, poromechanics, packing size, drying rate, capillarity, and substrate friction on cracking. Guided by the theory, we also found that cracking can be completely altered by varying the spatial profile of drying. Our work elucidates the rich physics underlying cracking in shrinkable, granular packings, and yields new strategies for controlling crack evolution.
△ Less
Submitted 15 May, 2019;
originally announced May 2019.
-
Modeling hydrodynamic interactions in soft materials with multiparticle collision dynamics
Authors:
Michael P. Howard,
Arash Nikoubashman,
Jeremy C. Palmer
Abstract:
Multiparticle collision dynamics (MPCD) is a flexible and robust mesoscale computational technique for simulating solvent-mediated hydrodynamic interactions in soft materials. Here, we provide a critical overview of the MPCD method and summarize its current strengths and limitations. The capabilities of the method are highlighted by reviewing its recent applications to simulate diverse phenomena,…
▽ More
Multiparticle collision dynamics (MPCD) is a flexible and robust mesoscale computational technique for simulating solvent-mediated hydrodynamic interactions in soft materials. Here, we provide a critical overview of the MPCD method and summarize its current strengths and limitations. The capabilities of the method are highlighted by reviewing its recent applications to simulate diverse phenomena, ranging from the flow of complex fluids and thermo-osmotic transport to bacterial swimming and active particle self-assembly. We also discuss outstanding challenges and emerging methodological developments that are expected to greatly expand the applicability of MPCD to other systems of technological importance.
△ Less
Submitted 27 February, 2019;
originally announced February 2019.
-
Quantized bounding volume hierarchies for neighbor search in molecular simulations on graphics processing units
Authors:
Michael P. Howard,
Antonia Statt,
Felix Madutsa,
Thomas M. Truskett,
Athanassios Z. Panagiotopoulos
Abstract:
We present an algorithm for neighbor search in molecular simulations on graphics processing units (GPUs) based on bounding volume hierarchies (BVHs). The BVH is compressed into a low-precision, quantized representation to increase the BVH traversal speed compared to a previous implementation. We find that neighbor search using the quantized BVH is roughly two to four times faster than current stat…
▽ More
We present an algorithm for neighbor search in molecular simulations on graphics processing units (GPUs) based on bounding volume hierarchies (BVHs). The BVH is compressed into a low-precision, quantized representation to increase the BVH traversal speed compared to a previous implementation. We find that neighbor search using the quantized BVH is roughly two to four times faster than current state-of-the-art methods using uniform grids (cell lists) for a suite of benchmarks for common molecular simulation models. Based on the benchmark results, we recommend using the BVH instead of a single cell list for neighbor list generation in molecular simulations on GPUs.
△ Less
Submitted 25 March, 2019; v1 submitted 23 January, 2019;
originally announced January 2019.
-
Cross-stream migration of a Brownian droplet in a polymer solution under Poiseuille flow
Authors:
Michael P. Howard,
Thomas M. Truskett,
Arash Nikoubashman
Abstract:
The migration of a Brownian fluid droplet in a parallel-plate microchannel was investigated using dissipative particle dynamics computer simulations. In a Newtonian solvent, the droplet migrated toward the channel walls due to inertial effects at the studied flow conditions, in agreement with theoretical predictions and recent simulations. However, the droplet focused onto the channel centerline w…
▽ More
The migration of a Brownian fluid droplet in a parallel-plate microchannel was investigated using dissipative particle dynamics computer simulations. In a Newtonian solvent, the droplet migrated toward the channel walls due to inertial effects at the studied flow conditions, in agreement with theoretical predictions and recent simulations. However, the droplet focused onto the channel centerline when polymer chains were added to the solvent. Focusing was typically enhanced for longer polymers and higher polymer concentrations with a nontrivial flow-rate dependence due to droplet and polymer deformability. Brownian motion caused the droplet position to fluctuate with a distribution that primarily depended on the balance between inertial lift forces pushing the droplet outward and elastic forces from the polymers driving it inward. The droplet shape was controlled by the local shear rate, and so its average shape depended on the droplet distribution.
△ Less
Submitted 17 December, 2018;
originally announced December 2018.
-
Unexpected secondary flows in reverse nonequilibrium shear flow simulations
Authors:
Antonia Statt,
Michael P. Howard,
Athanassios Z. Panagiotopoulos
Abstract:
We simulated two particle-based fluid models, namely multiparticle collision dynamics and dissipative particle dynamics, under shear using reverse nonequilibrium simulations (RNES). In cubic periodic simulation boxes, the expected shear flow profile for a Newtonian fluid developed, consistent with the fluid viscosities. However, unexpected secondary flows along the shear gradient formed when the s…
▽ More
We simulated two particle-based fluid models, namely multiparticle collision dynamics and dissipative particle dynamics, under shear using reverse nonequilibrium simulations (RNES). In cubic periodic simulation boxes, the expected shear flow profile for a Newtonian fluid developed, consistent with the fluid viscosities. However, unexpected secondary flows along the shear gradient formed when the simulation box was elongated in the flow direction. The standard shear flow profile was obtained when the simulation box was longer in the shear-gradient dimension than the flow dimension, while the secondary flows were always present when the flow dimension was at least 25% larger than the shear-gradient dimension. The secondary flows satisfy the boundary conditions imposed by the RNES and have a lower rate of viscous dissipation in the fluid than the corresponding unidirectional flows. This work highlights a previously unappreciated limitation of RNES for generating shear flow in simulation boxes that are elongated in the flow dimension, an important consideration when applying RNES to complex fluids like polymer solutions.
△ Less
Submitted 9 November, 2018;
originally announced November 2018.
-
Influence of hydrodynamic interactions on stratification in drying mixtures
Authors:
Antonia Statt,
Michael P. Howard,
Athanassios Z. Panagiotopoulos
Abstract:
Nonequilibrium molecular dynamics simulations are used to investigate the influence of hydrodynamic interactions on vertical segregation (stratification) in drying mixtures of long and short polymer chains. In agreement with previous computer simulations and theoretical modeling, the short polymers stratify on top of the long polymers at the top of the drying film when hydrodynamic interactions be…
▽ More
Nonequilibrium molecular dynamics simulations are used to investigate the influence of hydrodynamic interactions on vertical segregation (stratification) in drying mixtures of long and short polymer chains. In agreement with previous computer simulations and theoretical modeling, the short polymers stratify on top of the long polymers at the top of the drying film when hydrodynamic interactions between polymers are neglected. However, no stratification occurs at the same drying conditions when hydrodynamic interactions are incorporated through an explicit solvent model. Our analysis demonstrates that models lacking hydrodynamic interactions do not faithfully represent stratification in drying mixtures, in agreement with recent analysis of an idealized model for diffusiophoresis, and must be incorporated into such models in future.
△ Less
Submitted 30 March, 2018;
originally announced April 2018.
-
Efficient mesoscale hydrodynamics: multiparticle collision dynamics with massively parallel GPU acceleration
Authors:
Michael P. Howard,
Athanassios Z. Panagiotopoulos,
Arash Nikoubashman
Abstract:
We present an efficient open-source implementation of the multiparticle collision dynamics (MPCD) algorithm that scales to run on hundreds of graphics processing units (GPUs). We especially focus on optimizations for modern GPU architectures and communication patterns between multiple GPUs. We show that a mixed-precision computing model can improve performance compared to a fully double-precision…
▽ More
We present an efficient open-source implementation of the multiparticle collision dynamics (MPCD) algorithm that scales to run on hundreds of graphics processing units (GPUs). We especially focus on optimizations for modern GPU architectures and communication patterns between multiple GPUs. We show that a mixed-precision computing model can improve performance compared to a fully double-precision model while still providing good numerical accuracy. We report weak and strong scaling benchmarks of a reference MPCD solvent and a benchmark of a polymer solution with research-relevant interactions and system size. Our MPCD software enables simulations of mesoscale hydrodynamics at length and time scales that would be otherwise challenging or impossible to access.
△ Less
Submitted 6 March, 2018; v1 submitted 30 January, 2018;
originally announced January 2018.
-
Coupling of nanoparticle dynamics to polymer center-of-mass motion in semidilute polymer solutions
Authors:
Renjie Chen,
Ryan Poling-Skutvik,
Arash Nikoubashman,
Michael P. Howard,
Jacinta C. Conrad,
Jeremy C. Palmer
Abstract:
We investigate the dynamics of nanoparticles in semidilute polymer solutions when the nanoparticles are comparably sized to the polymer coils using explicit- and implicit-solvent simulation methods. The nanoparticle dynamics are subdiffusive on short time scales before transitioning to diffusive motion on long time scales. The long-time diffusivities scale according to theoretical predictions base…
▽ More
We investigate the dynamics of nanoparticles in semidilute polymer solutions when the nanoparticles are comparably sized to the polymer coils using explicit- and implicit-solvent simulation methods. The nanoparticle dynamics are subdiffusive on short time scales before transitioning to diffusive motion on long time scales. The long-time diffusivities scale according to theoretical predictions based on full dynamic coupling to the polymer segmental relaxations. In agreement with our recent experiments, however, we observe that the nanoparticle subdiffusive exponents are significantly larger than predicted by the coupling theory over a broad range of polymer concentrations. We attribute this discrepancy in the subdiffusive regime to the presence of an additional coupling mechanism between the nanoparticle dynamics and the polymer center-of-mass motion, which differs from the polymer relaxations that control the long-time diffusion. This coupling is retained even in the absence of many-body hydrodynamic interactions when the long-time dynamics of the colloids and polymers are matched.
△ Less
Submitted 9 February, 2018; v1 submitted 2 November, 2017;
originally announced November 2017.