Hierarchically Doped Plasmonic Nanocrystal Metamaterials.

Assembling plasmonic nanocrystals in regular superlattices can produce effective optical properties not found in homogeneous materials. However, the range of these metamaterial properties is limited when a single nanocrystal composition is selected for the constituent meta-atoms. Here, we show how continuously varying doping at two length scales, the atomic and nanocrystal scales, enables tuning of both the frequency and bandwidth of the collective plasmon resonance in nanocrystal-based metasurfaces, while these features are inextricably linked in single-component superlattices. Varying the mixing ratio of indium tin oxide nanocrystals with different dopant concentrations, we use large-scale simulations to predict the emergence of a broad infrared spectral region with near-zero permittivity. Experimentally, tunable reflectance and absorption bands are observed, owing to in- and out-of-plane collective resonances. These spectral features and the predicted strong near-field enhancement establish this multiscale doping strategy as a powerful new approach to designing metamaterials for optical applications.

Plasmonic Response of Complex Nanoparticle Assemblies.

Optical properties of nanoparticle assemblies reflect distinctive characteristics of their building blocks and spatial organization, giving rise to emergent phenomena. Integrated experimental and computational studies have established design principles connecting the structure to properties for assembled clusters and superlattices. However, conventional electromagnetic simulations are too computationally expensive to treat more complex assemblies. Here we establish a fast, materials agnostic method to simulate the optical response of large nanoparticle assemblies incorporating both structural and compositional complexity. This many-bodied, mutual polarization method resolves limitations of established approaches, achieving rapid, accurate convergence for configurations including thousands of nanoparticles, with some overlapping. We demonstrate these capabilities by reproducing experimental trends and uncovering far- and near-field mechanisms governing the optical response of plasmonic semiconductor nanocrystal assemblies including structurally complex gel networks and compositionally complex mixed binary superlattices. This broadly applicable framework will facilitate the design of complex, hierarchically structured, and dynamic assemblies for desired optical characteristics.

Characterizing Protein-Protein Interactions and Viscosity of a Monoclonal Antibody from Low to High Concentration Using Small-Angle X-ray Scattering and Molecular Dynamics Simulations.

Understanding protein-protein interactions and formation of reversible oligomers (clusters) in concentrated monoclonal antibody (mAb) solutions is necessary for designing stable, low viscosity (η) concentrated formulations for processing and subcutaneous injection. Here we characterize the strength (K) of short-range anisotropic attractions (SRA) for 75-200 mg/mL mAb2 solutions at different pH and cosolute conditions by analyzing structure factors (Seff(q)) from small-angle X-ray scattering (SAXS) using coarse-grained molecular dynamics simulations. Best fit simulations additionally provide cluster size distributions, fractal dimensions, cluster occluded volume, and mAb coordination numbers. These equilibrium properties are utilized in a model to account for increases in viscosity caused by occluded volume in the clusters (packing effects) and dissipation of stress across lubricated fractal clusters. Seff(q) is highly sensitive to K at 75 mg/mL where mAbs can mutually align to form SRA contacts but becomes less sensitive at 200 mg/mL as steric repulsion due to packing becomes dominant. In contrast, η at 200 mg/mL is highly sensitive to SRA and the average cluster size from SAXS/simulation, which is observed to track the cluster relaxation time from shear thinning. By analyzing the distribution of sub-bead hot spots on the 3D mAb surface, we identify a strongly attractive hydrophobic patch in the complementarity determining region (CDR) at pH 4.5 that contributes to the high K and consequently large cluster sizes and high η. Adding NaCl screens electrostatic interactions and increases the impact of hydrophobic attraction on cluster size and raises η, whereas nonspecific binding of Arg attenuates all SRA, reducing η. The hydrophobic patch is absent at higher pH values, leading to smaller K, smaller clusters, and lower η. This work constitutes a first attempt to use SAXS and CG modeling to link both structural and rheological properties of concentrated mAb solutions to the energetics of specific hydrophobic patches on mAb surfaces. As such, our work opens an avenue for future research, including the possibility of designing coarse-grained models with physically meaningful interacting hot spots.

Subclass Effects on Self-Association and Viscosity of Monoclonal Antibodies at High Concentrations.

The effects of a subclass of monoclonal antibodies (mAbs) on protein-protein interactions, formation of reversible oligomers (clusters), and viscosity (η) are not well understood at high concentrations. Herein, we quantify a short-range anisotropic attraction between the complementarity-determining region (CDR) and CH3 domains (KCDR-CH3) for vedolizumab IgG1, IgG2, or IgG4 subclasses by fitting small-angle X-ray scattering (SAXS) structure factor Seff(q) data with an extensive library of 12-bead coarse-grained (CG) molecular dynamics simulations. The KCDR-CH3 bead attraction strength was isolated from the strength of long-range electrostatic repulsion for the full mAb, which was determined from the theoretical net charge and a scaling parameter ψ to account for solvent accessibility and ion pairing. At low ionic strength (IS), the strongest short-range attraction (KCDR-CH3) and consequently the largest clusters and highest η were observed with IgG1, the subclass with the most positively charged CH3 domain. Furthermore, the trend in KCDR-CH3 with the subclass followed the electrostatic interaction energy between the CDR and CH3 regions calculated with the BioLuminate software using the 3D mAb structure and molecular interaction potentials. Whereas the equilibrium cluster size distributions and fractal dimensions were determined from fits of SAXS with the MD simulations, the degree of cluster rigidity under flow was estimated from the experimental η with a phenomenological model. For the systems with the largest clusters, especially IgG1, the inefficient packing of mAbs in the clusters played the largest role in increasing η, whereas for other systems, the relative contribution from stress produced by the clusters was more significant. The ability to relate η to short-range attraction from SAXS measurements at high concentrations and to theoretical characterization of electrostatic patches on the 3D surface is not only of fundamental interest but also of practical value for mAb discovery, processing, formulation, and subcutaneous delivery.

Correction: Magnetic field enabled in situ control over the structure and dynamics of colloids interacting via SALR potentials.

Correction for 'Magnetic field enabled in situ control over the structure and dynamics of colloids interacting via SALR potentials' by Hashir M. Gauri et al., Soft Matter, 2023, 19, 4439-4448, https://doi.org/10.1039/D3SM00354J.

Magnetic field enabled in situ control over the structure and dynamics of colloids interacting via SALR potentials.

Colloidal suspensions are an ideal model for studying crystallization, nucleation, and glass transition mechanisms, due to the precise control of interparticle interactions by changing the shape, charge, or volume fraction of particles. However, these tuning parameters offer insufficient active control over interparticle interactions and reconfigurability of assembled structures. Dynamic control over the interparticle interactions can be obtained through the application of external magnetic fields that are contactless and chemically inert. In this work, we demonstrate the dual nature of magnetic nanoparticle dispersions to program interactions between suspended nonmagnetic microspheres using an external magnetic field. The nanoparticle dispersion simultaneously behaves as a continuous magnetic medium at the microscale and a discrete medium composed of individual particles at the nanoscale. This enables control over a depletion attractive potential and the introduction of a magnetic repulsive potential, allowing a reversible transition of colloidal structures within a rich phase diagram by applying an external magnetic field. Active control over competing interactions allows us to create a model system encompassing a range of states, from large fractal clusters to low-density Wigner glass states. Monitoring the dynamics of colloidal particles reveals dynamic heterogeneity and a marked slowdown associated with approaching the Wigner glass state.

Multiscale modeling of solute diffusion in triblock copolymer membranes.

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.

Modular mixing in plasmonic metal oxide nanocrystal gels with thermoreversible links.

Gelation offers a powerful strategy to assemble plasmonic nanocrystal networks incorporating both the distinctive optical properties of constituent building blocks and customizable collective properties. Beyond what a single-component assembly can offer, the characteristics of nanocrystal networks can be tuned in a broader range when two or more components are intimately combined. Here, we demonstrate mixed nanocrystal gel networks using thermoresponsive metal-terpyridine links that enable rapid gel assembly and disassembly with thermal cycling. Plasmonic indium oxide nanocrystals with different sizes, doping concentrations, and shapes are reliably intermixed in linked gel assemblies, exhibiting collective infrared absorption that reflects the contributions of each component while also deviating systematically from a linear combination of the spectra for single-component gels. We extend a many-bodied, mutual polarization method to simulate the optical response of mixed nanocrystal gels, reproducing the experimental trends with no free parameters and revealing that spectral deviations originate from cross-coupling between nanocrystals with distinct plasmonic properties. Our thermoreversible linking strategy directs the assembly of mixed nanocrystal gels with continuously tunable far- and near-field optical properties that are distinct from those of the building blocks or mixed close-packed structures.

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 the design of reconfigurable materials and fueled assemblies.

Colloidal Nanocrystal Gels from Thermodynamic Principles.

Gels assembled from solvent-dispersed nanocrystals are of interest for functional materials because they promise the opportunity to retain distinctive properties of individual nanocrystals combined with tunable, structure-dependent collective behavior. By incorporating stimuli-responsive components, these materials could also be dynamically reconfigured between structurally distinct states. However, nanocrystal gels have so far been formed mostly through irreversible aggregation, which has limited the realization of these possibilities. Meanwhile, gelation strategies for larger colloidal microparticles have been developed using reversible physical or chemical interactions. These approaches have enabled the experimental navigation of theoretically predicted phase diagrams, helping to establish an understanding of how thermodynamic behavior can guide gel formation in these materials. However, the translation of these principles to the nanoscale poses both practical and fundamental challenges. The molecules guiding assembly can no longer be safely assumed to be vanishingly small compared to the particles nor large compared to the solvent.In this Account, we discuss recent progress toward the assembly of tunable nanocrystal gels using two strategies guided by equilibrium considerations: (1) reversible chemical bonding between functionalized nanocrystals and difunctional linker molecules and (2) nonspecific, polymer-induced depletion attractions. The effective nanocrystal attractions, mediated in both approaches by a secondary molecule, compete against stabilizing repulsions to promote reversible assembly. The structure and properties of the nanocrystal gels are controlled microscopically by the design of the secondary molecule and macroscopically by its concentration. This mode of control is compelling because it largely decouples nanocrystal synthesis and functionalization from the design of interactions that drive assembly. Statistical thermodynamic theory and computer simulation have been applied to simple models that describe the bonding motifs in these assembling systems, furnish predictions for conditions under which gelation is likely to occur, and suggest strategies for tuning and disassembling the gel networks. Insights from these models have guided experimental realizations of reversible gels with optical properties in the infrared range that are sensitive to the gel structure. This process avoids time-consuming and costly trial-and-error experimental investigations to accelerate the development of nanocrystal gel assemblies.These advances highlight the need to better understand interactions between nanocrystals, how interactions give rise to gel structure, and properties that emerge. Such an understanding could suggest new approaches for creating stimuli-responsive and dissipative assembled materials whose properties are tunable on demand through directed reconfiguration of the underlying gel microstructure. It may also make nanocrystal gels amenable to computationally guided design using inverse methods to rapidly optimize experimental parameters for targeted functionalities.

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.

Dynamics of equilibrium-linked colloidal networks.

Colloids that attractively bond to only a few neighbors (e.g., patchy particles) can form equilibrium gels with distinctive dynamic properties that are stable in time. Here, we use a coarse-grained model to explore the dynamics of linked networks of patchy colloids whose average valence is macroscopically, rather than microscopically, constrained. Simulation results for the model show dynamic hallmarks of equilibrium gel formation and establish that the colloid-colloid bond persistence time controls the characteristic slow relaxation of the self-intermediate scattering function. The model features re-entrant network formation without phase separation as a function of linker concentration, centered at the stoichiometric ratio of linker ends to nanoparticle surface bonding sites. Departures from stoichiometry result in linker-starved or linker-saturated networks with reduced connectivity and shorter characteristic relaxation times with lower activation energies. Underlying the re-entrant trends, dynamic properties vary monotonically with the number of effective network bonds per colloid, a quantity that can be predicted using Wertheim's thermodynamic perturbation theory. These behaviors suggest macroscopic in situ strategies for tuning the dynamic response of colloidal networks.

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.

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 the 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.

Gelation of plasmonic metal oxide nanocrystals by polymer-induced depletion attractions.

Gelation of colloidal nanocrystals emerged as a strategy to preserve inherent nanoscale properties in multiscale architectures. However, available gelation methods to directly form self-supported nanocrystal networks struggle to reliably control nanoscale optical phenomena such as photoluminescence and localized surface plasmon resonance (LSPR) across nanocrystal systems due to processing variabilities. Here, we report on an alternative gelation method based on physical internanocrystal interactions: short-range depletion attractions balanced by long-range electrostatic repulsions. The latter are established by removing the native organic ligands that passivate tin-doped indium oxide (ITO) nanocrystals while the former are introduced by mixing with small PEG chains. As we incorporate increasing concentrations of PEG, we observe a reentrant phase behavior featuring two favorable gelation windows; the first arises from bridging effects while the second is attributed to depletion attractions according to phase behavior predicted by our unified theoretical model. Our assembled nanocrystals remain discrete within the gel network, based on X-ray scattering and high-resolution transmission electron microscopy. The infrared optical response of the gels is reflective of both the nanocrystal building blocks and the network architecture, being characteristic of ITO nanocrystals' LSPR with coupling interactions between neighboring nanocrystals.

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 nanocrystals' surfaces, and are therefore not easily transferable. Here, we report a general gelation strategy via nonspecific 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.

Coarse-Grained Molecular Dynamics Simulations for Understanding the Impact of Short-Range Anisotropic Attractions on Structure and Viscosity of Concentrated Monoclonal Antibody Solutions.

Understanding protein-protein interactions in concentrated therapeutic monoclonal antibody (mAb) solutions is desirable for improved drug discovery, processing, and administration. Here, we deduce both the net protein charge and the magnitude and geometry of short-ranged, anisotropic attractions of a mAb across multiple concentrations and cosolute conditions by comparing structure factors S(q) obtained from small-angle X-ray scattering experiments with those from molecular dynamics (MD) simulations. The simulations, which utilize coarse-grained 12-bead models exhibiting a uniform van der Waals attraction, uniform electrostatic repulsion, and short-range attractions between specific beads, are versatile enough to fit S(q) of a wide range of protein concentrations and ionic strength with the same charge on each bead and a single anisotropic short-range attraction strength. Cluster size distributions (CSDs) obtained from best fit simulations reveal that the experimental structure is consistent with small reversible oligomers in even low viscosity systems and help quantify the impact of these clusters on viscosity. The ability to systematically use experimental S(q) data together with MD simulations to discriminate between different possible protein-protein interactions, as well as to predict viscosities from protein CSDs, is beneficial for designing mAbs and developing formulation strategies that avoid high viscosities and aggregation at high concentration.

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.

Influence of pore morphology on the diffusion of water in triblock copolymer membranes.

Understanding the transport properties of water in self-assembled block copolymer morphologies is important for furthering the use of such materials as water-purifying membranes. In this study, we used coarse-grained dissipative particle dynamics simulations to clarify the influence of pore morphology on the self-diffusion of water in linear-triblock-copolymer membranes. We considered representative lamellar, cylindrical, and gyroid morphologies and present results for both the global and local diffusivities of water in the pores. Our results suggest that the diffusivity of water in the confined, polymer-coated pores differs from that in the unconfined bulk. Explicitly, in confinement, the mobility of water is reduced by the hydrodynamic friction arising from the hydrophilic blocks coating the pore walls. We demonstrate that in lamella and cylindrical morphologies, the latter effects can be rendered as a universal function of the pore size relative to the brush height of the hydrophilic blocks.

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 a 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.