1
|
Rivera-Rivera LY, Moore TC, Glotzer SC. Inverse design of triblock Janus spheres for self-assembly of complex structures in the crystallization slot via digital alchemy. SOFT MATTER 2023; 19:2726-2736. [PMID: 36974942 DOI: 10.1039/d2sm01593e] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/18/2023]
Abstract
The digital alchemy framework is an extended ensemble simulation technique that incorporates particle attributes as thermodynamic variables, enabling the inverse design of colloidal particles for desired behavior. Here, we extend the digital alchemy framework for the inverse design of patchy spheres that self-assemble into target crystal structures. To constrain the potentials to non-trivial solutions, we conduct digital alchemy simulations with constant second virial coefficient. We optimize the size, range, and strength of patchy interactions in model triblock Janus spheres to self-assemble the 2D kagome and snub square lattices and the 3D pyrochlore lattice, and demonstrate self-assembly of all three target structures with the designed models. The particles designed for the kagome and snub square lattices assemble into high quality clusters of their target structures, while competition from similar polymorphs lower the yield of the pyrochlore assemblies. We find that the alchemically designed potentials do not always match physical intuition, illustrating the ability of the method to find nontrivial solutions to the optimization problem. We identify a window of second virial coefficients that result in self-assembly of the target structures, analogous to the crystallization slot in protein crystallization.
Collapse
Affiliation(s)
| | - Timothy C Moore
- Department of Chemical Engineering, University of Michigan, Ann Arbor, MI, USA.
| | - Sharon C Glotzer
- Department of Chemical Engineering, University of Michigan, Ann Arbor, MI, USA.
- Biointerfaces Institute, University of Michigan, Ann Arbor, MI, USA
| |
Collapse
|
2
|
Sharma AK, Escobedo FA. Effect of particle anisotropy on the thermodynamics and kinetics of ordering transitions in hard faceted particles. J Chem Phys 2023; 158:044502. [PMID: 36725523 DOI: 10.1063/5.0135461] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/07/2023] Open
Abstract
Monte Carlo simulations were used to study the influence of particle aspect ratio on the kinetics and phase behavior of hard gyrobifastigia (GBF). First, the formation of a highly anisotropic nucleus shape in the isotropic-to-crystal transition in regular GBF is explained by the differences in interfacial free energies of various crystal planes and the nucleus geometry predicted by the Wulff construction. GBF-related shapes with various aspect ratios were then studied, mapping their equations of state, determining phase coexistence conditions via interfacial pinning, and computing nucleation free-energy barriers via umbrella sampling using suitable order parameters. Our simulations reveal a reduction of the kinetic barrier for isotropic-crystal transition upon an increase in aspect ratio, and that for highly oblate and prolate aspect ratios, an intermediate nematic phase is stabilized. Our results and observations also support two conjectures for the formation of the crystalline state from the isotropic phase: that low phase free energies at the ordering phase transition correlate with low transition barriers and that the emergence of a mesophase provides a steppingstone that expedites crystallization.
Collapse
Affiliation(s)
- Abhishek K Sharma
- R. F. Smith School of Chemical and Biomolecular Engineering, Cornell University, Ithaca, New York 14853, USA
| | - Fernando A Escobedo
- R. F. Smith School of Chemical and Biomolecular Engineering, Cornell University, Ithaca, New York 14853, USA
| |
Collapse
|
3
|
Escobedo FA. On the calculation of free energies over Hamiltonian and order parameters via perturbation and thermodynamic integration. J Chem Phys 2021; 155:114112. [PMID: 34551542 DOI: 10.1063/5.0061541] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
In this work, complementary formulas are presented to compute free-energy differences via perturbation (FEP) methods and thermodynamic integration (TI). These formulas are derived by selecting only the most statistically significant data from the information extractable from the simulated points involved. On the one hand, commonly used FEP techniques based on overlap sampling leverage the full information contained in the overlapping macrostate probability distributions. On the other hand, conventional TI methods only use information on the first moments of those distributions, as embodied by the first derivatives of the free energy. Since the accuracy of simulation data degrades considerably for high-order moments (for FEP) or free-energy derivatives (for TI), it is proposed to consider, consistently for both methods, data up to second-order moments/derivatives. This provides a compromise between the limiting strategies embodied by common FEP and TI and leads to simple, optimized expressions to evaluate free-energy differences. The proposed formulas are validated with an analytically solvable harmonic Hamiltonian (for assessing systematic errors), an atomistic system (for computing the potential of mean force with coordinate-dependent order parameters), and a binary-component coarse-grained model (for tracing a solid-liquid phase diagram in an ensemble sampled through alchemical transformations). It is shown that the proposed FEP and TI formulas are straightforward to implement, perform similarly well, and allow robust estimation of free-energy differences even when the spacing of successive points does not guarantee them to have proper overlapping in phase space.
Collapse
Affiliation(s)
- Fernando A Escobedo
- Robert Frederick Smith School of Chemistry and Biomolecular Engineering, Cornell University, Ithaca, New York 14853, USA
| |
Collapse
|
4
|
Lindquist BA. Inverse design of equilibrium cluster fluids applied to a physically informed model. J Chem Phys 2021; 154:174907. [PMID: 34241069 DOI: 10.1063/5.0048812] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Inverse design strategies have proven highly useful for the discovery of interaction potentials that prompt self-assembly of a variety of interesting structures. However, often the optimized particle interactions do not have a direct relationship to experimental systems. In this work, we show that Relative Entropy minimization is able to discover physically meaningful parameter sets for a model interaction built from depletion attraction and electrostatic repulsion that yield self-assembly of size-specific clusters. We then explore the sensitivity of the optimized interaction potentials with respect to deviations in the underlying physical quantities, showing that clustering behavior is largely preserved even as the optimized parameters are perturbed.
Collapse
Affiliation(s)
- Beth A Lindquist
- Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, USA
| |
Collapse
|
5
|
Moore TC, Anderson JA, Glotzer SC. Shape-driven entropic self-assembly of an open, reconfigurable, binary host-guest colloidal crystal. SOFT MATTER 2021; 17:2840-2848. [PMID: 33564812 DOI: 10.1039/d0sm02073g] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/12/2023]
Abstract
Entropically driven self-assembly of hard anisotropic particles, where particle shape gives rise to emergent valencies, provides a useful perspective for the design of nanoparticle and colloidal systems. Hard particles self-assemble into a rich variety of crystal structures, ranging in complexity from simple close-packed structures to structures with 432 particles in the unit cell. Entropic crystallization of open structures, however, is missing from this landscape. Here, we report the self-assembly of a two-dimensional binary mixture of hard particles into an open host-guest structure, where nonconvex, triangular host particles form a honeycomb lattice that encapsulates smaller guest particles. Notably, this open structure forms in the absence of enthalpic interactions by effectively splitting the structure into low- and high-entropy sublattices. This is the first such structure to be reported in a two-dimensional athermal system. We discuss the observed compartmentalization of entropy in this system, and show that the effect of the size of the guest particle on the stability of the structure gives rise to a reentrant phase behavior. This reentrance suggests the possibility for a reconfigurable colloidal material, and we provide a proof-of-concept by showing the assembly behavior while changing the size of the guest particles in situ. Our findings provide a strategy for designing open colloidal crystals, as well as binary systems that exhibit co-crystallization, which have been elusive thus far.
Collapse
Affiliation(s)
- Timothy C Moore
- Department of Chemical Engineering, University of Michigan, Ann Arbor, MI, USA.
| | - Joshua A Anderson
- Department of Chemical Engineering, University of Michigan, Ann Arbor, MI, USA.
| | - Sharon C Glotzer
- Department of Chemical Engineering, University of Michigan, Ann Arbor, MI, USA. and Biointerfaces Institute, University of Michigan, Ann Arbor, MI, USA
| |
Collapse
|
6
|
Ma Y, Aulicino JC, Ferguson AL. Inverse Design of Self-Assembling Diamond Photonic Lattices from Anisotropic Colloidal Clusters. J Phys Chem B 2021; 125:2398-2410. [PMID: 33647208 DOI: 10.1021/acs.jpcb.0c08723] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022]
Abstract
Colloidal nanoparticles with anisotropic interactions are promising building blocks for the fabrication of complex functional materials. A challenge in the self-assembly of colloidal particles is the rational design of geometry and chemistry to program the formation of a desired target structure. We report an inverse design procedure integrating Langevin dynamics simulations and evolutionary algorithms to engineer anisotropic patchy colloidal clusters to spontaneously assemble into a cubic diamond lattice possessing a complete photonic band gap. The combination of a tetrahedral cluster geometry and optimized placement of a single type of anisotropic interaction patch results in a colloidal building block predicted to assemble a cubic diamond lattice with more than 82% yield. This design represents an experimentally viable colloidal building block capable of high-fidelity assembly of a cubic diamond lattice.
Collapse
Affiliation(s)
- Yutao Ma
- Pritzker School of Molecular Engineering, University of Chicago, Chicago, Illinois 60637, United States
| | - Joseph C Aulicino
- Pritzker School of Molecular Engineering, University of Chicago, Chicago, Illinois 60637, United States
| | - Andrew L Ferguson
- Pritzker School of Molecular Engineering, University of Chicago, Chicago, Illinois 60637, United States
| |
Collapse
|
7
|
Quintela Matos I, Escobedo F. Congruent phase behavior of a binary compound crystal of colloidal spheres and dimpled cubes. J Chem Phys 2020; 153:214503. [DOI: 10.1063/5.0030174] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Isabela Quintela Matos
- School of Chemical and Biomolecular Engineering, Cornell University, Ithaca, New York 14853, USA
| | - Fernando Escobedo
- School of Chemical and Biomolecular Engineering, Cornell University, Ithaca, New York 14853, USA
| |
Collapse
|
8
|
Gutiérrez AB, Chapela GA. Effect of shape on liquid-vapor coexistence and surface properties of parallelepiped molecules. J Chem Phys 2020; 152:134501. [PMID: 32268761 DOI: 10.1063/1.5141352] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/03/2023] Open
Abstract
Liquid-vapor coexistence is calculated via molecular dynamics for a variety of parallelepiped shaped molecules. Models are constructed as an array of tangential hard spheres interacting with an attractive square-well potential. Each shape is formed by varying the number of spheres in their three sides. The initial density of the system is chosen close to the critical density of a SW fluid to obtain an equilibrated liquid-vapor coexistence curve by the process of spinodal decomposition. A pattern that relates the geometry of the molecular models and the existence or non-existence of a liquid-vapor orthobaric curve is shown.
Collapse
Affiliation(s)
- Anthony B Gutiérrez
- Departamento de Física, Universidad Autónoma Metropolitana-Iztapalapa, Av. San Rafael Atlixco 186, Col. Vicentina, 09340 México D.F., Mexico
| | - Gustavo A Chapela
- Departamento de Física, Universidad Autónoma Metropolitana-Iztapalapa, Av. San Rafael Atlixco 186, Col. Vicentina, 09340 México D.F., Mexico
| |
Collapse
|
9
|
Ma Y, Ferguson AL. Inverse design of self-assembling colloidal crystals with omnidirectional photonic bandgaps. SOFT MATTER 2019; 15:8808-8826. [PMID: 31603182 DOI: 10.1039/c9sm01500k] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/10/2023]
Abstract
Open colloidal lattices possessing omnidirectional photonic bandgaps in the visible or near-visible regime are attractive optical materials the realization of which has remained elusive. We report the use of an inverse design strategy termed landscape engineering that rationally sculpts the free energy self-assembly landscape using evolutionary algorithms to discover anisotropic patchy colloids capable of spontaneously assembling pyrochlore and cubic diamond lattices possessing complete photonic bandgaps. We validate the designs in computer simulations to demonstrate the defect-free formation of these lattices via a two-stage hierarchical assembly mechanism. Our approach demonstrates a principled strategy for the inverse design of self-assembling colloids for the bottom-up fabrication of desired crystal lattices.
Collapse
Affiliation(s)
- Yutao Ma
- Pritzker School of Molecular Engineering, University of Chicago, 5640 South Ellis Avenue, Chicago, IL 60637, USA.
| | - Andrew L Ferguson
- Pritzker School of Molecular Engineering, University of Chicago, 5640 South Ellis Avenue, Chicago, IL 60637, USA.
| |
Collapse
|
10
|
Lindquist BA, Jadrich RB, Howard MP, Truskett TM. The role of pressure in inverse design for assembly. J Chem Phys 2019; 151:104104. [DOI: 10.1063/1.5112766] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023] Open
Affiliation(s)
- Beth A. Lindquist
- McKetta Department of Chemical Engineering, University of Texas at Austin, Austin, Texas 78712, USA
- Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, USA
| | - Ryan B. Jadrich
- McKetta Department of Chemical Engineering, University of Texas at Austin, Austin, Texas 78712, USA
- Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, USA
- Center for Nonlinear Studies, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, USA
| | - Michael P. Howard
- McKetta Department of Chemical Engineering, University of Texas at Austin, Austin, Texas 78712, USA
| | - Thomas M. Truskett
- McKetta Department of Chemical Engineering, University of Texas at Austin, Austin, Texas 78712, USA
- Department of Physics, University of Texas at Austin, Austin, Texas 78712, USA
| |
Collapse
|
11
|
Chen D, Zhang G, Torquato S. Inverse Design of Colloidal Crystals via Optimized Patchy Interactions. J Phys Chem B 2018; 122:8462-8468. [DOI: 10.1021/acs.jpcb.8b05627] [Citation(s) in RCA: 22] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
|