Inverse design of a pyrochlore lattice of DNA origami through model-driven experiments

Sophisticated statistical mechanics approaches and human intuition have demonstrated the possibility of self-assembling complex lattices or finite-size constructs. However, attempts so far have mostly only been successful in silico and often fail in experiment because of unpredicted traps associated with kinetic slowing down (gelation, glass transition) and competing ordered structures. Theoretical predictions also face the difficulty of encoding the desired interparticle interaction potential with the experimentally available nano- and micrometer-sized particles. To overcome these issues, we combine SAT assembly (a patchy-particle interaction design algorithm based on constrained optimization) with coarse-grained simulations of DNA nanotechnology to experimentally realize trap-free self-assembly pathways. We use this approach to assemble a pyrochlore three-dimensional lattice, coveted for its promise in the construction of optical metamaterials, and characterize it with small-angle x-ray scattering and scanning electron microscopy visualization.

Structure and elasticity of model disordered, polydisperse, and defect-free polymer networks

The elasticity of disordered and polydisperse polymer networks is a fundamental problem of soft matter physics that is still open. Here, we self-assemble polymer networks via simulations of a mixture of bivalent and tri- or tetravalent patchy particles, which result in an exponential strand length distribution analogous to that of experimental randomly cross-linked systems. After assembly, the network connectivity and topology are frozen and the resulting system is characterized. We find that the fractal structure of the network depends on the number density at which the assembly has been carried out, but that systems with the same mean valence and same assembly density have the same structural properties. Moreover, we compute the long-time limit of the mean-squared displacement, also known as the (squared) localization length, of the cross-links and of the middle monomers of the strands, showing that the dynamics of long strands is well described by the tube model. Finally, we find a relation connecting these two localization lengths at high density and connect the cross-link localization length to the shear modulus of the system.

The oxDNA Coarse-Grained Model as a Tool to Simulate DNA Origami

This chapter introduces how to run molecular dynamics simulations for DNA origami using the oxDNA coarse-grained model.

A simple solution to the problem of self-assembling cubic diamond crystals

The self-assembly of colloidal diamond (CD) crystals is considered as one of the most coveted goals of nanotechnology, both from the technological and fundamental points of view. For applications, colloidal diamond is a photonic crystal which can open new possibilities of manipulating light for information processing. From a fundamental point of view, its unique symmetry exacerbates a series of problems that are commonly faced during the self-assembly of target structures, such as the presence of kinetic traps and the formation of crystalline defects and alternative structures (polymorphs). Here we demonstrate that all these problems can be systematically addressed via SAT-assembly, a design framework that converts self-assembly into a Boolean satisfiability problem (SAT). Contrary to previous solutions (requiring four or more components), we prove that the assembly of the CD crystal only requires a binary mixture. Moreover, we use molecular dynamics simulations of a system composed by nearly a million nucleotides to test a DNA nanotechnology design that constitutes a promising candidate for experimental realization.

Designing Enhanced Entropy Binding in Single-Chain Nanoparticles

Single-chain nanoparticles (SCNPs) are a new class of bio- and soft-matter polymeric objects in which a fraction of the monomers are able to form equivalently intra- or interpolymer bonds. Here we numerically show that a fully entropic gas-liquid phase separation can take place in SCNP systems. Control over the discontinuous (first-order) change-from a phase of independent diluted (fully-bonded) polymers to a phase in which polymers entropically bind to each other to form a (fully-bonded) polymer network-can be achieved by a judicious design of the patterns of reactive monomers along the polymer chain. Such a sensitivity arises from a delicate balance between the distinct entropic contributions controlling the binding.

SAT-assembly: a new approach for designing self-assembling systems

We propose a general framework for solving inverse self-assembly problems, i.e. designing interactions between elementary units such that they assemble spontaneously into a predetermined structure. Our approach uses patchy particles as building blocks, where the different units bind at specific interaction sites (the patches), and we exploit the possibility of having mixtures with several components. The interaction rules between the patches is determined by transforming the combinatorial problem into a Boolean satisfiability problem (SAT) which searches for solutions where all bonds are formed in the target structure. Additional conditions, such as the non-satisfiability of competing structures (e.g. metastable states) can be imposed, allowing to effectively design the assembly path in order to avoid kinetic traps. We demonstrate this approach by designing and numerically simulating a cubic diamond structure from four particle species that assembles without competition from other polymorphs, including the hexagonal structure.

Hydrodynamic instability and flow reduction in polymer brush coated channels

A polymer brush is a passive medium. At equilibrium the knowledge of its chemical composition and thickness is enough for a full system characterization. However, when the brush is exposed to fluid flow it reveals a much more intriguing nature, in which filamentous protrusions and the way they interact among themselves and with the surrounding fluid are of outmost importance. Here we investigate such a rich behavior via numerical simulations. We focus on the brush hydrodynamic response at low Reynolds numbers, observing a significant fluid flow reduction inside a polymer-brush coated channel. We find that the reduction of the flow inside the channel is significantly larger than what would happen if the brush effect consisted only in reducing the effective channel width. This amplified reduction is understood as being due to the morphological instability of the brush-liquid interface which is shown to have an elastic origin: the mechanical stress acting on the brush due to the imposed flow is partially released by the interface modulation. In turn, this modulation dissipates more energy than a flat interface in the surrounding fluid, causing a reduction of flow velocity. Our results and interpretations provide an explanation for recent experimental measurements.

Sharing the Load: Stress Redistribution Governs Fracture of Polymer Double Networks

The stress response of polymer double networks depends not only on the properties of the constituent networks but also on the interactions arising between them. Here, we demonstrate, via coarse-grained simulations, that both their global stress response and their microscopic fracture mechanics are governed by load sharing through these internetwork interactions. By comparing our results with affine predictions, where stress redistribution is by definition homogeneous, we show that stress redistribution is highly inhomogeneous. In particular, the affine prediction overestimates the fraction of broken chains by almost an order of magnitude. Furthermore, homogeneous stress distribution predicts a single fracture process, while in our simulations, fracture of sacrificial chains takes place in two steps governed by load sharing within a network and between networks, respectively. Our results thus provide a detailed microscopic picture of how inhomogeneous stress redistribution after rupture of chains governs the fracture of polymer double networks.

Effect of Internal Architecture on the Assembly of Soft Particles at Fluid Interfaces

Monolayers of soft colloidal particles confined at fluid interfaces are at the core of a broad range of technological processes, from the stabilization of responsive foams and emulsions to advanced lithographic techniques. However, establishing a fundamental relation between their internal architecture, which is controlled during synthesis, and their structural and mechanical properties upon interfacial confinement remains an elusive task. To address this open issue, which defines the monolayer's properties, we synthesize core-shell microgels, whose soft core can be chemically degraded in a controlled fashion. This strategy allows us to obtain a series of particles ranging from analogues of standard batch-synthesized microgels to completely hollow ones after total core removal. Combined experimental and numerical results show that our hollow particles have a thin and deformable shell, leading to a temperature-responsive collapse of the internal cavity and a complete flattening after adsorption at a fluid interface. Mechanical characterization shows that a critical degree of core removal is required to obtain soft disk-like particles at an oil-water interface, which present a distinct response to compression. At low packing fractions, the mechanical response of the monolayer is dominated by the outer polymer chains forming a corona surrounding the particles within the interfacial plane, regardless of the presence of a core. By contrast, at high compression, the absence of a core enables the particles to deform in the direction orthogonal to the interface and to be continuously compressed without altering the monolayer structure. These findings show how fine, single-particle architectural control during synthesis can be engineered to determine the interfacial behavior of microgels, enabling one to link particle conformation with the resulting material properties.

OxDNA.org: a public webserver for coarse-grained simulations of DNA and RNA nanostructures

OxDNA and oxRNA are popular coarse-grained models used by the DNA/RNA nanotechnology community to prototype, analyze and rationalize designed DNA and RNA nanostructures. Here, we present oxDNA.org, a graphical web interface for running, visualizing and analyzing oxDNA and oxRNA molecular dynamics simulations on a GPU-enabled high performance computing server. OxDNA.org automatically generates simulation files, including a multi-step relaxation protocol for structures exported in non-physical states from DNA/RNA design tools. Once the simulation is complete, oxDNA.org provides an interactive visualization and analysis interface using the browser-based visualizer oxView to facilitate the understanding of simulation results for a user's specific structure. This online tool significantly lowers the entry barrier of integrating simulations in the nanostructure design pipeline for users who are not experts in the technical aspects of molecular simulation. The webserver is freely available at oxdna.org.

A Primer on the oxDNA Model of DNA: When to Use it, How to Simulate it and How to Interpret the Results

The oxDNA model of Deoxyribonucleic acid has been applied widely to systems in biology, biophysics and nanotechnology. It is currently available via two independent open source packages. Here we present a set of clearly documented exemplar simulations that simultaneously provide both an introduction to simulating the model, and a review of the model's fundamental properties. We outline how simulation results can be interpreted in terms of-and feed into our understanding of-less detailed models that operate at larger length scales, and provide guidance on whether simulating a system with oxDNA is worthwhile.

Effect of Chain Polydispersity on the Elasticity of Disordered Polymer Networks

Due to their unique structural and mechanical properties, randomly cross-linked polymer networks play an important role in many different fields, ranging from cellular biology to industrial processes. In order to elucidate how these properties are controlled by the physical details of the network (e.g., chain-length and end-to-end distributions), we generate disordered phantom networks with different cross-linker concentrations C and initial densities ρinit and evaluate their elastic properties. We find that the shear modulus computed at the same strand concentration for networks with the same C, which determines the number of chains and the chain-length distribution, depends strongly on the preparation protocol of the network, here controlled by ρinit. We rationalize this dependence by employing a generic stress-strain relation for polymer networks that does not rely on the specific form of the polymer end-to-end distance distribution. We find that the shear modulus of the networks is a nonmonotonic function of the density of elastically active strands, and that this behavior has a purely entropic origin. Our results show that if short chains are abundant, as it is always the case for randomly cross-linked polymer networks, the knowledge of the exact chain conformation distribution is essential for correctly predicting the elastic properties. Finally, we apply our theoretical approach to literature experimental data, qualitatively confirming our interpretations.

Charge affinity and solvent effects in numerical simulations of ionic microgels

Ionic microgel particles are intriguing systems in which the properties of thermo-responsive polymeric colloids are enriched by the presence of charged groups. In order to rationalize their properties and predict the behaviour of microgel suspensions, it is necessary to develop a coarse-graining strategy that starts from the accurate modelling of single particles. Here, we provide a numerical advancement of a recently-introduced model for charged co-polymerized microgels by improving the treatment of ionic groups in the polymer network. We investigate the thermoresponsive properties of the particles, in particular their swelling behaviour and structure, finding that, when charged groups are considered to be hydrophilic at all temperatures, highly charged microgels do not achieve a fully collapsed state, in favorable comparison to experiments. In addition, we explicitly include the solvent in the description and put forward a mapping between the solvophobic potential in the absence of the solvent and the monomer-solvent interactions in its presence, which is found to work very accurately for any charge fraction of the microgel. Our work paves the way for comparing single-particle properties and swelling behaviour of ionic microgels to experiments and to tackle the study of these charged soft particles at a liquid-liquid interface.

Shape control of soft patchy nanoparticles under confinement

Molecular building blocks undergoing a hierarchical assembly process form nano-scale objects which can further assemble into supramolecular structures. When the intermediate units have a limited valence in bonding, complex structures with tailored properties can be created. Here, we consider a composite, star-shaped particle made of f diblock copolymer chains uniformly grafted on a spherical colloid and investigate its first self-assembly stage both in the bulk and under lateral confinement. By means of numerical simulations, we show that, in the bulk, this system develops aggregates whose number and size depend on the temperature as well as on the relative ratio of solvophobic monomers. The emerging aggregates are referred to as patches and impart directionality in bonding to the complex particle. We further characterize how we can control, by changing the lateral confinement, the shape of the brush and the patch properties as a function of the distance between the confining walls. We find that the number of the patches can be determined by tuning the degree of confinement imposed on the particle. Finally, we employ a continuum mechanics model, known as the Liquid Drop Model, to gain insight into the elastic properties of the system. This theoretical approach allows to connect the patch properties to the elastic response of the composite particle.

Designer protein assemblies with tunable phase diagrams in living cells

Protein self-organization is a hallmark of biological systems. Although the physicochemical principles governing protein-protein interactions have long been known, the principles by which such nanoscale interactions generate diverse phenotypes of mesoscale assemblies, including phase-separated compartments, remain challenging to characterize. To illuminate such principles, we create a system of two proteins designed to interact and form mesh-like assemblies. We devise a new strategy to map high-resolution phase diagrams in living cells, which provide self-assembly signatures of this system. The structural modularity of the two protein components allows straightforward modification of their molecular properties, enabling us to characterize how interaction affinity impacts the phase diagram and material state of the assemblies in vivo. The phase diagrams and their dependence on interaction affinity were captured by theory and simulations, including out-of-equilibrium effects seen in growing cells. Finally, we find that cotranslational protein binding suffices to recruit a messenger RNA to the designed micron-scale structures.

Numerical insights on ionic microgels: structure and swelling behaviour

Recent progress has been made in the numerical modelling of neutral microgel particles with a realistic, disordered structure. In this work we extend this approach to the case of co-polymerised microgels where a thermoresponsive polymer is mixed with acidic groups. We compare the cases where counterions directly interact with microgel charges or are modelled implicitly through a Debye-Hückel description. We do so by performing extensive numerical simulations of single microgels across the volume phase transition (VPT) varying the temperature and the fraction of charged monomers. We find that the presence of charges considerably alters the microgel structure, quantified by the monomer density profiles and by the form factors of the microgels, particularly close to the VPT. We observe significant deviations between the implicit and explicit models, with the latter comparing more favourably to available experiments. In particular, we observe a shift of the VPT temperature to larger values as the amount of charged monomers increases. We also find that below the VPT the microgel-counterion complex is almost neutral, while it develops a net charge above the VPT. Interestingly, under these conditions the collapsed microgel still retains a large amount of counterions inside its structure. Since these interesting features cannot be captured by the implicit model, our results show that it is crucial to explicitly include the counterions in order to realistically model ionic thermoresponsive microgels.

Modeling Microgels with a Controlled Structure across the Volume Phase Transition

Thermoresponsive microgels are soft colloids that find widespread use as model systems for soft matter physics. Their complex internal architecture, made of a disordered and heterogeneous polymer network, has been so far a major challenge for computer simulations. In this work, we put forward a coarse-grained model of microgels whose structural properties are in quantitative agreement with results obtained with small-angle X-ray scattering experiments across a wide range of temperatures, encompassing the volume phase transition. These results bridge the gap between experiments and simulations of individual microgel particles, paving the way to theoretically address open questions about their bulk properties with unprecedented nano- and microscale resolution.

TacoxDNA: A user-friendly web server for simulations of complex DNA structures, from single strands to origami

Simulations of nucleic acids at different levels of structural details are increasingly used to complement and interpret experiments in different fields, from biophysics to medicine and materials science. However, the various structural models currently available for DNA and RNA and their accompanying suites of computational tools can be very rarely used in a synergistic fashion. The tacoxDNA webserver and standalone software package presented here are a step toward a long-sought interoperability of nucleic acids models. The webserver offers a simple interface for converting various common input formats of DNA structures and setting up molecular dynamics (MD) simulations. Users can, for instance, design DNA rings with different topologies, such as knots, with and without supercoiling, by simply providing an XYZ coordinate file of the DNA centre-line. More complex DNA geometries, as designable in the cadnano, CanDo and Tiamat tools, can also be converted to all-atom or oxDNA representations, which can then be used to run MD simulations. Though the latter are currently geared toward the native and LAMMPS oxDNA representations, the open-source package is designed to be further expandable. TacoxDNA is available at http://tacoxdna.sissa.it. © 2019 Wiley Periodicals, Inc.

q-Independent Slow Dynamics in Atomic and Molecular Systems

Investigating million-atom systems for very long simulation times, we demonstrate that the collective density-density correlation time (τ_{α}) in simulated supercooled water and silica becomes wave-vector independent (q^{0}) when the probing wavelength is several times larger than the interparticle distance. The q independence of the collective density-density correlation functions, a feature clearly observed in light-scattering studies of some soft-matter systems, is thus a genuine feature of many (but not all) slow-dynamics systems, either atomic, molecular, or colloidal. Indeed, we show that when the dynamics of the density fluctuations includes particle-type diffusion, as in the case of the Lennard-Jones binary-mixture model, the q^{0} regime does not set in and the relaxation time continues to scale as τ_{α}∼q^{-2} even at small q.

Microgels Adsorbed at Liquid-Liquid Interfaces: A Joint Numerical and Experimental Study

Soft particles display highly versatile properties with respect to hard colloids and even more so at fluid-fluid interfaces. In particular, microgels, consisting of a cross-linked polymer network, are able to deform and flatten upon adsorption at the interface due to the balance between surface tension and internal elasticity. Despite the existence of experimental results, a detailed theoretical understanding of this phenomenon is still lacking due to the absence of appropriate microscopic models. In this work, we propose an advanced modeling of microgels at a flat water/oil interface. The model builds on a realistic description of the internal polymeric architecture and single-particle properties of the microgel and is able to reproduce its experimentally observed shape at the interface. Complementing molecular dynamics simulations with in situ cryo-electron microscopy experiments and atomic force microscopy imaging after Langmuir-Blodgett deposition, we compare the morphology of the microgels for different values of the cross-linking ratios. Our model allows for a systematic microscopic investigation of soft particles at fluid interfaces, which is essential to develop predictive power for the use of microgels in a broad range of applications, including the stabilization of smart emulsions and the versatile patterning of surfaces.

Numerical modelling of non-ionic microgels: an overview

Microgels are complex macromolecules. These colloid-sized polymer networks possess internal degrees of freedom and, depending on the polymer(s) they are made of, can acquire a responsiveness to variations of the environment (temperature, pH, salt concentration, etc.). Besides being valuable for many practical applications, microgels are also extremely important to tackle fundamental physics problems. As a result, these last years have seen a rapid development of protocols for the synthesis of microgels, and more and more research has been devoted to the investigation of their bulk properties. However, from a numerical standpoint the picture is more fragmented, as the inherently multi-scale nature of microgels, whose bulk behaviour crucially depends on the microscopic details, cannot be handled at a single level of coarse-graining. Here we present an overview of the methods and models that have been proposed to describe non-ionic microgels at different length-scales, from the atomistic to the single-particle level. We especially focus on monomer-resolved models, as these have the right level of details to capture the most important properties of microgels, responsiveness and softness. We suggest that these microscopic descriptions, if realistic enough, can be employed as starting points to develop the more coarse-grained representations required to investigate the behaviour of bulk suspensions.

A new look at effective interactions between microgel particles

Thermoresponsive microgels find widespread use as colloidal model systems, because their temperature-dependent size allows facile tuning of their volume fraction in situ. However, an interaction potential unifying their behavior across the entire phase diagram is sorely lacking. Here we investigate microgel suspensions in the fluid regime at different volume fractions and temperatures, and in the presence of another population of small microgels, combining confocal microscopy experiments and numerical simulations. We find that effective interactions between microgels are clearly temperature dependent. In addition, microgel mixtures possess an enhanced stability compared to hard colloid mixtures - a property not predicted by a simple Hertzian model. Based on numerical calculations we propose a multi-Hertzian model, which reproduces the experimental behavior for all studied conditions. Our findings highlight that effective interactions between microgels are much more complex than usually assumed, displaying a crucial dependence on temperature and on the internal core-corona architecture of the particles.

Modelling realistic microgels in an explicit solvent

Thermoresponsive microgels are polymeric colloidal networks that can change their size in response to a temperature variation. This peculiar feature is driven by the nature of the solvent-polymer interactions, which triggers the so-called volume phase transition from a swollen to a collapsed state above a characteristic temperature. Recently, an advanced modelling protocol to assemble realistic, disordered microgels has been shown to reproduce experimental swelling behavior and form factors. In the original framework, the solvent was taken into account in an implicit way, condensing solvent-polymer interactions in an effective attraction between monomers. To go one step further, in this work we perform simulations of realistic microgels in an explicit solvent. We identify a suitable model which fully captures the main features of the implicit model and further provides information on the solvent uptake by the interior of the microgel network and on its role in the collapse kinetics. These results pave the way for addressing problems where solvent effects are dominant, such as the case of microgels at liquid-liquid interfaces.

Force-Induced Unravelling of DNA Origami

The mechanical properties of DNA nanostructures are of widespread interest as applications that exploit their stability under constant or intermittent external forces become increasingly common. We explore the force response of DNA origami in comprehensive detail by combining AFM single molecule force spectroscopy experiments with simulations using oxDNA, a coarse-grained model of DNA at the nucleotide level, to study the unravelling of an iconic origami system: the Rothemund tile. We contrast the force-induced melting of the tile with simulations of an origami 10-helix bundle. Finally, we simulate a recently proposed origami biosensor, whose function takes advantage of origami behavior under tension. We observe characteristic stick-slip unfolding dynamics in our force-extension curves for both the Rothemund tile and the helix bundle and reasonable agreement with experimentally observed rupture forces for these systems. Our results highlight the effect of design on force response: we observe regular, modular unfolding for the Rothemund tile that contrasts with strain-softening of the 10-helix bundle which leads to catastropic failure under monotonically increasing force. Further, unravelling occurs straightforwardly from the scaffold ends inward for the Rothemund tile, while the helix bundle unfolds more nonlinearly. The detailed visualization of the yielding events provided by simulation allows preferred pathways through the complex unfolding free-energy landscape to be mapped, as a key factor in determining relative barrier heights is the extensional release per base pair broken. We shed light on two important questions: how stable DNA nanostructures are under external forces and what design principles can be applied to enhance stability.

On the effect of the thermostat in non-equilibrium molecular dynamics simulations

The numerical investigation of the statics and dynamics of systems in non-equilibrium in general, and under shear flow in particular, has become more and more common. However, not all the numerical methods developed to simulate equilibrium systems can be successfully adapted to out-of-equilibrium cases. This is especially true for thermostats. Indeed, even though thermostats developed to work under equilibrium conditions sometimes display good agreement with rheology experiments, their performance rapidly degrades beyond weak dissipation and small shear rates. Here we focus on gauging the relative performances of three thermostats, Langevin, dissipative particle dynamics, and Bussi-Donadio-Parrinello under varying parameters and external conditions. We compare their effectiveness by looking at different observables and clearly demonstrate that choosing the right thermostat (and its parameters) requires a careful evaluation of, at least, temperature, density and velocity profiles. We also show that small modifications of the Langevin and DPD thermostats greatly enhance their performance in a wide range of parameters.

How to simulate patchy particles⋆

Patchy particles is the name given to a large class of systems of mesoscopic particles characterized by a repulsive core and a discrete number of short-range and highly directional interaction sites. Numerical simulations have contributed significantly to our understanding of the behaviour of patchy particles, but, although simple in principle, advanced simulation techniques are often required to sample the low temperatures and long time-scales associated with their self-assembly behaviour. In this work we review the most popular simulation techniques that have been used to study patchy particles, with a special focus on Monte Carlo methods. We cover many of the tools required to simulate patchy systems, from interaction potentials to biased moves, cluster moves, and free-energy methods. The review is complemented by an educationally oriented Monte Carlo computer code that implements all the techniques described in the text to simulate a well-known tetrahedral patchy particle model.

An Accurate Estimate of the Free Energy and Phase Diagram of All-DNA Bulk Fluids

We present a numerical study in which large-scale bulk simulations of self-assembled DNA constructs have been carried out with a realistic coarse-grained model. The investigation aims at obtaining a precise, albeit numerically demanding, estimate of the free energy for such systems. We then, in turn, use these accurate results to validate a recently proposed theoretical approach that builds on a liquid-state theory, the Wertheim theory, to compute the phase diagram of all-DNA fluids. This hybrid theoretical/numerical approach, based on the lowest-order virial expansion and on a nearest-neighbor DNA model, can provide, in an undemanding way, a parameter-free thermodynamic description of DNA associating fluids that is in semi-quantitative agreement with experiments. We show that the predictions of the scheme are as accurate as those obtained with more sophisticated methods. We also demonstrate the flexibility of the approach by incorporating non-trivial additional contributions that go beyond the nearest-neighbor model to compute the DNA hybridization free energy.

Internal structure and swelling behaviour of in silico microgel particles

Microgels are soft colloids that, by virtue of their polymeric nature, can react to external stimuli such as temperature or pH by changing their size. The resulting swelling/deswelling transition can be exploited in fundamental research as well as for many diverse practical applications, ranging from art restoration to medicine. Such an extraordinary versatility stems from the complex internal structure of the individual microgels, each of which is a crosslinked polymer network. Here we employ a recently-introduced computational method to generate realistic microgel configurations and look at their structural properties, both in real and Fourier space, for several temperatures across the volume phase transition as a function of the crosslinker concentration and of the confining radius employed during the 'in-silico' synthesis. We find that the chain-length distribution of the resulting networks can be analytically predicted by a simple theoretical argument. In addition, we find that our results are well-fitted to the fuzzy-sphere model, which correctly reproduces the density profile of the microgels under study.

In Silico Synthesis of Microgel Particles

Microgels are colloidal-scale particles individually made of cross-linked polymer networks that can swell and deswell in response to external stimuli, such as changes to temperature or pH. Despite a large amount of experimental activities on microgels, a proper theoretical description based on individual particle properties is still missing due to the complexity of the particles. To go one step further, here we propose a novel methodology to assemble realistic microgel particles in silico. We exploit the self-assembly of a binary mixture composed of tetravalent (cross-linkers) and bivalent (monomer beads) patchy particles under spherical confinement in order to produce fully bonded networks. The resulting structure is then used to generate the initial microgel configuration, which is subsequently simulated with a bead-spring model complemented by a temperature-induced hydrophobic attraction. To validate our assembly protocol, we focus on a small microgel test case and show that we can reproduce the experimental swelling curve by appropriately tuning the confining sphere radius, something that would not be possible with less sophisticated assembly methodologies, e.g., in the case of networks generated from an underlying crystal structure. We further investigate the structure (in reciprocal and real space) and the swelling curves of microgels as a function of temperature, finding that our results are well described by the widely used fuzzy sphere model. This is a first step toward a realistic modeling of microgel particles, which will pave the way for a careful assessment of their elastic properties and effective interactions.

Free energy calculations for rings and chains formed by dipolar hard spheres

We employ a method based on Monte Carlo grand-canonical simulations to precisely calculate partition functions of non-interacting chains and rings formed by dipolar hard spheres (DHS) at low temperature. The extended low temperature region offered by such cluster calculations, compared to what had been previously achieved with standard simulations, opens up the possibility of exploring a part of the DHS phase diagram which was inaccessible before. The reported results offer the unique opportunity of verifying well-established theoretical models based on the ideal gas of cluster approximation in order to clarify their range of validity. They also provide the basis for future studies in which cluster-cluster interactions will be included.

Condensation and Demixing in Solutions of DNA Nanostars and Their Mixtures

We present a numerical/theoretical approach to efficiently evaluate the phase diagram of self-assembling DNA nanostars. Combining input information based on a realistic coarse-grained DNA potential with the Wertheim association theory, we derive a parameter-free thermodynamic description of these systems. We apply this method to investigate the phase behavior of single components and mixtures of DNA nanostars with different numbers of sticky arms, elucidating the role of the system functionality and of salt concentration. Specifically, we evaluate the propensity to demix, the gas-liquid phase boundaries and the location of the critical points. The predicted critical parameters compare very well with existing experimental results for the available compositions. The approach developed here is very general, easily extensible to other all-DNA systems, and provides guidance for future experiments.

Small-angle neutron scattering and molecular dynamics structural study of gelling DNA nanostars

DNA oligomers with properly designed sequences self-assemble into well defined constructs. Here, we exploit this methodology to produce bulk quantities of tetravalent DNA nanostars (each one composed of 196 nucleotides) and to explore the structural signatures of their aggregation process. We report small-angle neutron scattering experiments focused on the evaluation of both the form factor and the temperature evolution of the scattered intensity at a nanostar concentration where the system forms a tetravalent equilibrium gel. We also perform molecular dynamics simulations of one isolated tetramer to evaluate the form factor numerically, without resorting to any approximate shape. The numerical form factor is found to be in very good agreement with the experimental one. Simulations predict an essentially temperature-independent form factor, offering the possibility to extract the effective structure factor and its evolution during the equilibrium gelation.

Connectivity, dynamics, and structure in a tetrahedral network liquid

We report a detailed computational study by Brownian dynamics simulations of the structure and dynamics of a liquid of patchy particles which forms an amorphous tetrahedral network upon decreasing the temperature. The highly directional particle interactions allow us to investigate the system connectivity by discriminating the total set of particles into different populations according to a penta-modal distribution of bonds per particle. With this methodology we show how the particle bonding process is not randomly independent but it manifests clear bond correlations at low temperatures. We further explore the dynamics of the system in real space and establish a clear relation between particle mobility and particle connectivity. In particular, we provide evidence of anomalous diffusion at low temperatures and reveal how the dynamics is affected by the short-time hopping motion of the weakly bounded particles. Finally we widely investigate the dynamics and structure of the system in Fourier space and identify two quantitatively similar length scales, one dynamic and the other static, which increase upon cooling the system and reach distances of the order of few particle diameters. We summarize our findings in a qualitative picture where the low temperature regime of the viscoelastic liquid is understood in terms of an evolving network of long time metastable cooperative domains of particles.

Limiting the valence: advancements and new perspectives on patchy colloids, soft functionalized nanoparticles and biomolecules

Artistic representation of limited valance units consisting of a soft core (in blue) and a small number of flexible bonding patches (in orange).

Void-Based Assembly of Colloidal Crystals: Using Structure-Directing Agents to Direct the Assembly of Open Colloidal Crystals

We have used computer simulations to reveal how colloidal crystals may be assembled by engineering the shape and charge of polymeric additives, which act as structure-directing agents (SDAs). Using these agents, only a single desired polymorph may be obtained from a crystallizing mixture of colloids and polymers, which would otherwise result in a defective crystal in the absence of these agents. Building on previous work, which demonstrated this principle in the limiting case of high-density, close-packed crystals, we have now achieved this control over low-density, open crystals, which have broad utility especially in optical applications. These results reveal the general utility of the SDA paradigm in assembling tailored colloidal crystals, which represents a new design method that does not rely on any modification of the colloids themselves.

Bottom-Up Colloidal Crystal Assembly with a Twist

Globally ordered colloidal crystal lattices have broad utility in a wide range of optical and catalytic devices, for example, as photonic band gap materials. However, the self-assembly of stereospecific structures is often confounded by polymorphism. Small free-energy differences often characterize ensembles of different structures, making it difficult to produce a single morphology at will. Current techniques to handle this problem adopt one of two approaches: that of the "top-down" or "bottom-up" methodology, whereby structures are engineered starting from the largest or smallest relevant length scales, respectively. However, recently, a third approach for directing high fidelity assembly of colloidal crystals has been suggested which relies on the introduction of polymer cosolutes into the crystal phase [Mahynski, N.; Panagiotopoulos, A. Z.; Meng, D.; Kumar, S. K. Nat. Commun. 2014, 5, 4472]. By tuning the polymer's morphology to interact uniquely with the void symmetry of a single desired crystal, the entropy loss associated with polymer confinement has been shown to strongly bias the formation of that phase. However, previously, this approach has only been demonstrated in the limiting case of close-packed crystals. Here, we show how this approach may be generalized and extended to complex open crystals, illustrating the utility of this "structure-directing agent" paradigm in engineering the nanoscale structure of ordered colloidal materials. The high degree of transferability of this paradigm's basic principles between relatively simple crystals and more complex ones suggests that this represents a valuable addition to presently known self-assembly techniques.

Direct Simulation of the Self-Assembly of a Small DNA Origami

By using oxDNA, a coarse-grained nucleotide-level model of DNA, we are able to directly simulate the self-assembly of a small 384-base-pair origami from single-stranded scaffold and staple strands in solution. In general, we see attachment of new staple strands occurring in parallel, but with cooperativity evident for the binding of the second domain of a staple if the adjacent junction is already partially formed. For a system with exactly one copy of each staple strand, we observe a complete assembly pathway in an intermediate temperature window; at low temperatures successful assembly is prevented by misbonding while at higher temperature the free-energy barriers to assembly become too large for assembly on our simulation time scales. For high-concentration systems involving a large staple strand excess, we never see complete assembly because there are invariably instances where two copies of the same staple both bind to the scaffold, creating a kinetic trap that prevents the complete binding of either staple. This mutual staple blocking could also lead to aggregates of partially formed origamis in real systems, and helps to rationalize certain successful origami design strategies.

Soft self-assembled nanoparticles with temperature-dependent properties

The fabrication of versatile building blocks that reliably self-assemble into desired ordered and disordered phases is amongst the hottest topics in contemporary materials science. To this end, microscopic units of varying complexity, aimed at assembling the target phases, have been thought, designed, investigated and built. Such a path usually requires laborious fabrication techniques, especially when specific functionalisation of the building blocks is required. Telechelic star polymers, i.e., star polymers made of a number of f di-block copolymers consisting of solvophobic and solvophilic monomers grafted on a central anchoring point, spontaneously self-assemble into soft patchy particles featuring attractive spots (patches) on the surface. Here we show that the tunability of such a system can be widely extended by controlling the physical and chemical parameters of the solution. Indeed, under fixed external conditions the self-assembly behaviour depends only on the number of arms and on the ratio of solvophobic to solvophilic monomers. However, changes in temperature and/or solvent quality make it possible to reliably change the number and size of the attractive patches. This allows the steering of the mesoscopic self-assembly behaviour without modifying the microscopic constituents. Interestingly, we also demonstrate that diverse combinations of the parameters can generate stars with the same number of patches but different radial and angular stiffness. This mechanism could provide a neat way of further fine-tuning the elastic properties of the supramolecular network without changing its topology.

How soft repulsion enhances the depletion mechanism

We investigate binary mixtures of large colloids interacting through soft potentials with small, ideal depletants. We show that softness has a dramatic effect on the resulting colloid-colloid effective potential when the depletant-to-colloid size ratio q is small, with significant consequences on the colloidal phase behaviour. We provide an exact relationship that allows us to obtain the effective pair potential for any type of colloid-depletant interaction in the case of ideal depletants, without having to rely on complicated and expensive full-mixture simulations. We also show that soft repulsion among depletants further enhances the tendency of colloids to aggregate. Our theoretical and numerical results demonstrate that--in the limit of small q--soft mixtures cannot be mapped onto hard systems and hence soft depletion is not a mere extension of the widely used Asakura-Oosawa potential.

Temperature-induced structural transitions in self-assembling magnetic nanocolloids

With the help of a unique combination of density functional theory and computer simulations, we discover two possible scenarios, depending on concentration, for the hierarchical self-assembly of magnetic nanoparticles on cooling.

Gels of DNA nanostars never crystallize

Using state-of-the-art numerical techniques, we show that, upon lowering the temperature, tetravalent DNA nanostars form a thermodynamically stable, fully bonded equilibrium gel. In contrast to atomic and molecular network formers, in which the disordered liquid is always metastable with respect to some crystalline phase, we find that the DNA nanostar gel has a lower free energy than the diamond crystal structure in a wide range of concentrations. This unconventional behavior, here verified for the first time in a realistic model, arises from the large arm flexibility of the DNA nanostars, a property that can be tuned by design. Our results confirm the thermodynamic stability of the recently experimentally realized DNA hydrogels.