Observation of the most H2-dense filled ice under high pressure.

Hydrogen hydrates are among the basic constituents of our solar system's outer planets, some of their moons, as well Neptune-like exo-planets. The details of their high-pressure phases and their thermodynamic conditions of formation and stability are fundamental information for establishing the presence of hydrogen hydrates in the interior of those celestial bodies, for example, against the presence of the pure components (water ice and molecular hydrogen). Here, we report a synthesis path and experimental observation, by X-ray diffraction and Raman spectroscopy measurements, of the most H[Formula: see text]-dense phase of hydrogen hydrate so far reported, namely the compound 3 (or C[Formula: see text]). The detailed characterisation of this hydrogen-filled ice, based on the crystal structure of cubic ice I (ice I[Formula: see text]), is performed by comparing the experimental observations with first-principles calculations based on density functional theory and the stochastic self-consistent harmonic approximation. We observe that the extreme (up to 90 GPa and likely beyond) pressure stability of this hydrate phase is due to the close-packed geometry of the hydrogen molecules caged in the ice I[Formula: see text] skeleton.

Stacking Interactions and Flexibility of Human Telomeric Multimers.

G-quadruplexes (G4s) are helical four-stranded structures forming from guanine-rich nucleic acid sequences, which are thought to play a role in cancer development and malignant transformation. Most current studies focus on G4 monomers, yet under suitable and biologically relevant conditions, G4s undergo multimerization. Here, we investigate the stacking interactions and structural features of telomeric G4 multimers by means of a novel low-resolution structural approach that combines small-angle X-ray scattering (SAXS) with extremely coarse-grained (ECG) simulations. The degree of multimerization and the strength of the stacking interaction are quantitatively determined in G4 self-assembled multimers. We show that self-assembly induces a significant polydispersity of the G4 multimers with an exponential distribution of contour lengths, consistent with a step-growth polymerization. On increasing DNA concentration, the strength of the stacking interaction between G4 monomers increases, as well as the average number of units in the aggregates. We utilized the same approach to explore the conformational flexibility of a model single-stranded long telomeric sequence. Our findings indicate that its G4 units frequently adopt a beads-on-a-string configuration. We also observe that the interaction between G4 units can be significantly affected by complexation with benchmark ligands. The proposed methodology, which identifies the determinants that govern the formation and structural flexibility of G4 multimers, may be an affordable tool aiding in the selection and design of drugs that target G4s under physiological conditions.

On the role of polydispersity on the phase diagram of colloidal rods.

The rich and complex phase diagram typical of anisotropic biological or synthetic nanoparticles, has brought a great deal of interest over the equilibrium phase behaviour of non-spherical colloids. Amongst the class of anisotropic nanoparticles, hard spherocylindrical colloids have been, over the years, extensively studied because of their optical properties, for their rich phase diagrams, and their important industrial applications, as model particles for biological systems (viruses), or for example as potential drug carriers having the ability of surviving the attacks of the immune systems. As real anisotropic nanoparticles are often polydisperse in size and/or in shape, unveiling the effect of such a perturbation over their equilibrium phase diagram is of paramount importance. This work focuses on the effects of polydispersity over the full equilibrium phase diagram of hard spherocylindrical colloids (HSCs). Previous studies showed that a polydispersity in L alters the equilibrium phase diagram of HSCs. With this work we determine, both theoretically as well as computationally, the effects due to a generic polydispersity, namely in D, in L and, in both ones, on the equilibrium phase diagram and introduce a viable theoretical generalisation of the Onsager theory that allows us to get some insight into the observed phase behaviour.

Shape Matters in Magnetic-Field-Assisted Assembly of Prolate Colloids.

An anisotropic colloidal shape in combination with an externally tunable interaction potential results in a plethora of self-assembled structures with potential applications toward the fabrication of smart materials. Here we present our investigation on the influence of an external magnetic field on the self-assembly of hematite-silica core-shell prolate colloids for two aspect ratios ρ = 2.9 and 3.69. Our study shows a rather counterintuitive but interesting phenomenon, where prolate colloids self-assemble into oblate liquid crystalline (LC) phases. With increasing concentration, particles with smaller ρ reveal a sequence of LC phases involving para-nematic, nematic, smectic, and oriented glass phases. The occurrence of a smectic phase for colloidal ellipsoids has been neither predicted nor reported before. Quantitative shape analysis of the particles together with extensive computer simulations indicate that in addition to ρ, a subtle deviation from the ideal ellipsoidal shape dictates the formation of this unusual sequence of field-induced structures. Particles with ρ = 2.9 exhibit a hybrid shape containing features from both spherocylinders and ellipsoids, which make their self-assembly behavior richer than that observed for either of the "pure" shapes. The shape of the particles with higher ρ matches closely with the ideal ellipsoids, as a result their phase behavior follows the one expected for a "pure" ellipsoidal shape. Using anisotropic building blocks and external fields, our study demonstrates the ramifications of the subtle changes in the particle shape on the field-directed self-assembled structures with externally tunable properties.

Self-Assembly of All-DNA Rods with Controlled Patchiness.

Double-stranded DNA (dsDNA) fragments exhibit noncovalent attractive interactions between their tips. It is still unclear how DNA liquid crystal self-assembly is affected by such blunt-end attractions. It is demonstrated that stiff dsDNA fragments with moderate aspect ratio can specifically self-assemble in concentrated aqueous solutions into different types of smectic mesophases on the basis of selectively screening of blunt-end DNA stacking interactions. To this end, this type of attractions are engineered at the molecular level by constructing DNA duplexes where the attractions between one or both ends are screened by short hairpin caps. All-DNA bilayer and monolayer smectic-A type of phases, as well as a columnar phase, can be stabilized by controlling attractions strength. The results imply that the so far elusive smectic-A in DNA rod-like liquid crystals is a thermodynamically stable phase. The existence of the bilayer smectic phase is confirmed by Monte-Carlo simulations of hard cylinders decorated with one attractive terminal site. This work demonstrates that DNA blunt-ends behave as well-defined monovalent attractive patches whose strength and position can be potentially precisely tuned, highlighting unique opportunities concerning the stabilization of nonconventional DNA-based lyotropic liquid crystal phases assembled by all-DNA patchy particles with arbitrary geometry and composition.

Elastic constants of biological filamentous colloids: estimation and implications on nematic and cholesteric tactoid morphologies.

Biological liquid crystals, originating from the self-assembly of biological filamentous colloids, such as cellulose and amyloid fibrils, show a complex lyotropic behaviour that is extremely difficult to predict and characterize. Here we analyse the liquid crystalline phases of amyloid fibrils, and sulfated and carboxylated cellulose nanocrystals and measure their Frank-Oseen elastic constants K1, K2 and K3 by four different approaches. The first two approaches are based on the benchmark of the predictions of: (i) a scaling form and (ii) a variational form of the Frank-Oseen energy functional with the experimental critical volumes at order-order liquid crystalline transitions of the tactoids. The third and the fourth methods imply: (iii) the direct scaling equations of elastic constants and (iv) a molecular theory predicting the elastic constants from the experimentally accessible contour length distributions of the filamentous colloids. These three biological systems exhibit diverse liquid crystalline behaviour, governed by the distinct elastic constants characterizing each colloid. Differences and similarities among the three systems are highlighted and interpreted based on the present analysis, providing a general framework to study dispersed liquid crystalline systems.

Free energy of conformational isomers: The case of gapped DNA duplexes.

Liquid-crystalline phases in all-DNA systems have been extensively studied in the past and although nematic, cholesteric and columnar mesophases have been observed, the smectic phase remained elusive. Recently, it has been found evidence of a smectic-A ordering in an all-DNA system, where the constituent particles, which are gapped DNA duplexes, resemble chain-sticks. It has been argued that in the smectic-A phase these DNA chain-sticks should be folded as a means to suppress aggregate polydispersity and excluded volume. Nevertheless, if initial crystalline configurations are prepared in silico with gapped DNA duplexes either fully unfolded or fully folded by carrying out computer simulations one can end up with two different phases having at the same concentration and temperature the majority of gapped DNA duplexes either folded or unfolded. This result suggests that these two phases have a small free energy difference, since no transition is observed from one to the other within the simulation time span. In the present manuscript, we assess which of these two phases is thermodynamically stable through a suitable protocol based on thermodynamic integration. Our method is rather general and it can be used to discriminate stable states from metastable ones of comparable free energy.

Amyloid Fibrils Length Controls Shape and Structure of Nematic and Cholesteric Tactoids.

Amyloid fibrils offer the possibility of controlling their contour length, aspect ratio, and length distribution, without affecting other structural parameters. Here we show that a fine control in the contour length distribution of ß-lactoglobulin amyloid fibrils, achieved by mechanical shear stresses of different levels, translates into the organization of tactoids of different shapes and morphologies. While longer fibrils lead to highly elongated nematic tactoids in an isotropic continuous matrix, only sufficiently shortened amyloid fibrils lead to cholesteric droplets. The progressive decrease in amyloid fibrils length leads to a linear decrease of the anchoring strength and homogeneous tactoid → bipolar tactoid → cholesteric droplet transitions. Upon fibrils length increase, we first find experimentally and predict theoretically a decrease of the cholesteric pitch, before full disappearance of the cholesteric phase. The latter is understood to arise from the decrease of the energy barrier separating cholesteric and nematic phases over thermal energy for progressively longer, semiflexible fibrils.

Nematic liquid crystals of bifunctional patchy spheres.

Anisotropic interactions can bring about the formation, through self-assembly, of semi-flexible chains, which in turn can give rise to nematic phases for suitable temperatures and concentrations. A minimalist model constituted of hard cylinders decorated with attractive sites has been already extensively studied numerically. Simulation data shows that a theoretical approach recently proposed is able to properly capture the physical properties of these self-assembly-driven liquid crystals. Here, we investigated a simpler model constituted of bifunctional Kern-Frenkel hard spheres which does not possess steric anisotropy but which can undergo a istropic-nematic transition as a result of their self-assembly into semi-flexible chains. For this model we compare an accurate numerical estimate of isotropic-nematic phase boundaries with theoretical predictions. The theoretical treatment, originally proposed for cylinder-like particles, has been greatly simplified and its predictions are in good agreement with numerical results. Finally, we also assess a crucial, and not obvious, hypothesis used in the theory, i.e. the ability of the Onsager trial function to properly model particle orientation in the presence of aggregation, that has not been properly checked yet.

Exploiting limited valence patchy particles to understand autocatalytic kinetics.

Autocatalysis, i.e., the speeding up of a reaction through the very same molecule which is produced, is common in chemistry, biophysics, and material science. Rate-equation-based approaches are often used to model the time dependence of products, but the key physical mechanisms behind the reaction cannot be properly recognized. Here, we develop a patchy particle model inspired by a bicomponent reactive mixture and endowed with adjustable autocatalytic ability. Such a coarse-grained model captures all general features of an autocatalytic aggregation process that takes place under controlled and realistic conditions, including crowded environments. Simulation reveals that a full understanding of the kinetics involves an unexpected effect that eludes the chemistry of the reaction, and which is crucially related to the presence of an activation barrier. The resulting analytical description can be exported to real systems, as confirmed by experimental data on epoxy-amine polymerizations, solving a long-standing issue in their mechanistic description.

Speeding up Monte Carlo simulation of patchy hard cylinders.

The hard cylinder model decorated with attractive patches proved to be very useful recently in studying physical properties of several colloidal systems. Phase diagram, elastic constants and cholesteric properties obtained from computer simulations based on a simple hard cylinder model have been all successfully and quantitatively compared to experimental results. Key to these simulations is an efficient algorithm to check the overlap between hard cylinders. Here, we propose two algorithms to check the hard cylinder overlap and we assess their efficiency through a comparison with an existing method available in the literature and with the well-established algorithm for simulating hard spherocylinders. In addition, we discuss a couple of optimizations for performing computer simulations of patchy anisotropic particles and we estimate the speed-up which they can provide in the case of patchy hard cylinders.

Smectic phase in suspensions of gapped DNA duplexes.

Smectic ordering in aqueous solutions of monodisperse stiff double-stranded DNA fragments is known not to occur, despite the fact that these systems exhibit both chiral nematic and columnar mesophases. Here, we show, unambiguously, that a smectic-A type of phase is formed by increasing the DNA's flexibility through the introduction of an unpaired single-stranded DNA spacer in the middle of each duplex. This is unusual for a lyotropic system, where flexibility typically destabilizes the smectic phase. We also report on simulations suggesting that the gapped duplexes (resembling chain-sticks) attain a folded conformation in the smectic layers, and argue that this layer structure, which we designate as smectic-fA phase, is thermodynamically stabilized by both entropic and energetic contributions to the system's free energy. Our results demonstrate that DNA as a building block offers an exquisitely tunable means to engineer a potentially rich assortment of lyotropic liquid crystals.

Anomalous dynamics of intruders in a crowded environment of mobile obstacles.

Many natural and industrial processes rely on constrained transport, such as proteins moving through cells, particles confined in nanocomposite materials or gels, individuals in highly dense collectives and vehicular traffic conditions. These are examples of motion through crowded environments, in which the host matrix may retain some glass-like dynamics. Here we investigate constrained transport in a colloidal model system, in which dilute small spheres move in a slowly rearranging, glassy matrix of large spheres. Using confocal differential dynamic microscopy and simulations, here we discover a critical size asymmetry, at which anomalous collective transport of the small particles appears, manifested as a logarithmic decay of the density autocorrelation functions. We demonstrate that the matrix mobility is central for the observed anomalous behaviour. These results, crucially depending on size-induced dynamic asymmetry, are of relevance for a wide range of phenomena ranging from glassy systems to cell biology.

Simulation and Theory of Antibody Binding to Crowded Antigen-Covered Surfaces.

In this paper we introduce a fully flexible coarse-grained model of immunoglobulin G (IgG) antibodies parametrized directly on cryo-EM data and simulate the binding dynamics of many IgGs to antigens adsorbed on a surface at increasing densities. Moreover, we work out a theoretical model that allows to explain all the features observed in the simulations. Our combined computational and theoretical framework is in excellent agreement with surface-plasmon resonance data and allows us to establish a number of important results. (i) Internal flexibility is key to maximize bivalent binding, flexible IgGs being able to explore the surface with their second arm in search for an available hapten. This is made clear by the strongly reduced ability to bind with both arms displayed by artificial IgGs designed to rigidly keep a prescribed shape. (ii) The large size of IgGs is instrumental to keep neighboring molecules at a certain distance (surface repulsion), which essentially makes antigens within reach of the second Fab always unoccupied on average. (iii) One needs to account independently for the thermodynamic and geometric factors that regulate the binding equilibrium. The key geometrical parameters, besides excluded-volume repulsion, describe the screening of free haptens by neighboring bound antibodies. We prove that the thermodynamic parameters govern the low-antigen-concentration regime, while the surface screening and repulsion only affect the binding at high hapten densities. Importantly, we prove that screening effects are concealed in relative measures, such as the fraction of bivalently bound antibodies. Overall, our model provides a valuable, accurate theoretical paradigm beyond existing frameworks to interpret experimental profiles of antibodies binding to multi-valent surfaces of different sorts in many contexts.

Hierarchical Propagation of Chirality through Reversible Polymerization: The Cholesteric Phase of DNA Oligomers.

Unveiling the subtle rules that control the buildup of macroscopic chirality starting from chiral molecular elements is a challenge for theory and computations. In this context, a remarkable phenomenon is the formation of helically twisted nematic (cholesteric) phases, with pitch in the micrometer range, driven by self-assembly of relatively small chiral species into supramolecular semiflexible polymers. We have developed a theoretical framework to connect the cholesteric organization to the shape and chirality of the constituents, described with molecular detail, in this kind of system. The theory has been tested against new accurate measurements for solutions of short DNA duplexes. We show that the cholesteric organization is determined by steric repulsion between duplexes, and we identify distinctive features of linear self-assembly in the temperature and concentration dependence of the pitch.

Unusual dynamics of concentration fluctuations in solutions of weakly attractive globular proteins.

The globular protein γB-crystallin exhibits a complex phase behavior, where liquid-liquid phase separation characterized by a critical volume fraction ϕc = 0.154 and a critical temperature Tc = 291.8 K coexists with dynamical arrest on all length scales at volume fractions around ϕ ≈ 0.3-0.35, and an arrest line that extends well into the unstable region below the spinodal. However, although the static properties such as the osmotic compressibility and the static correlation length are in quantitative agreement with predictions for binary liquid mixtures, this is not the case for the dynamics of concentration fluctuations described by the dynamic structure factor S(q,t). Using a combination of dynamic light scattering and neutron spin echo measurements, we demonstrate that the competition between critical slowing down and dynamical arrest results in a much more complex wave vector dependence of S(q,t) than previously anticipated.

Nematic phase characterisation of the self-assembling sphere-cylinders based on the theoretically calculated RDFs.

We put forward a theoretical framework to calculate pair distribution functions in the nematic liquid crystals formed by sphere-cylinders that self-assemble in linear chain structures. For a nematically ordered system, one can distinguish between the spatial correlations in the plane perpendicular to the crystalline axis, and in the direction parallel to the latter. Following this separation, we show that the RDFs in the parallel case can be described using a superposition of a chain model and Onsager distribution, whereas the RDFs in the perpendicular case turn out to be that of the soft disks. Based on this concept, we show how the spatial correlations in the system are influenced by the nematic order parameter. We conclude that even if the nematic ordering is high in the system, the imperfection of the crystal is strongly reflected in the pair distributions.

Modelling the rheology of anisotropic particles adsorbed on a two-dimensional fluid interface.

We present a general approach based on nonequilibrium thermodynamics for bridging the gap between a well-defined microscopic model and the macroscopic rheology of particle-stabilised interfaces. Our approach is illustrated by starting with a microscopic model of hard ellipsoids confined to a planar surface, which is intended to simply represent a particle-stabilised fluid-fluid interface. More complex microscopic models can be readily handled using the methods outlined in this paper. From the aforementioned microscopic starting point, we obtain the macroscopic, constitutive equations using a combination of systematic coarse-graining, computer experiments and Hamiltonian dynamics. Exemplary numerical solutions of the constitutive equations are given for a variety of experimentally relevant flow situations to explore the rheological behaviour of our model. In particular, we calculate the shear and dilatational moduli of the interface over a wide range of surface coverages, ranging from the dilute isotropic regime, to the concentrated nematic regime.

Self-assembly of mesogenic bent-core DNA nanoduplexes.

Short cylinder-like DNA duplexes, comprising 6 to 20 base pairs, self-assemble into semi-flexible chains, due to coaxial stacking interactions between their blunt ends. The mutual alignment of these chains gives rise to macroscopically orientationally ordered liquid crystal phases. Interestingly, experiments show that the isotropic-nematic phase boundary is sequence-dependent. We perform all-atom simulations of several sequences to gain insights into the structural properties of the duplex and correlate the resulting geometric properties with the observed location of the isotropic-nematic phase boundary. We identify in the duplex bending the key parameter for explaining the sequence dependence, suggesting that DNA duplexes can be assimilated to bent-core mesogens. We also develop a coarse-grained model for the different DNA duplexes to evaluate in detail how bending affects the persistence length and excluded volume of the aggregates. This information is fed into a recently developed formalism to predict the isotropic-nematic phase boundary for bent-core mesogens. The theoretical results agree with the experimental observations.

Self-assembly of hard helices: a rich and unconventional polymorphism.

Hard helices can be regarded as a paradigmatic elementary model for a number of natural and synthetic soft matter systems, all featuring the helix as their basic structural unit, from natural polynucleotides and polypeptides to synthetic helical polymers, and from bacterial flagella to colloidal helices. Here we present an extensive investigation of the phase diagram of hard helices using a variety of methods. Isobaric Monte Carlo numerical simulations are used to trace the phase diagram; on going from the low-density isotropic to the high-density compact phases a rich polymorphism is observed, exhibiting a special chiral screw-like nematic phase and a number of chiral and/or polar smectic phases. We present full characterization of the latter, showing that they have unconventional features, ascribable to the helical shape of the constituent particles. Equal area construction is used to locate the isotropic-to-nematic phase transition, and the results are compared with those stemming from an Onsager-like theory. Density functional theory is also used to study the nematic-to-screw-nematic phase transition; within the simplifying assumption of perfectly parallel helices, we compare different levels of approximation, that is second- and third-virial expansions and a Parsons-Lee correction.