Supramolecular Packing Drives Morphological Transitions of Charged Surfactant Micelles.

The shape and size of self-assembled structures upon local organization of their molecular building blocks are hard to predict in the presence of long-range interactions. Combining small-angle X-ray/neutron scattering data, theoretical modelling, and computer simulations, sodium dodecyl sulfate (SDS), over a broad range of concentrations and ionic strengths, was investigated. Computer simulations indicate that micellar shape changes are associated with different binding of the counterions. By employing a toy model based on point charges on a surface, and comparing it to experiments and simulations, it is demonstrated that the observed morphological changes are caused by symmetry breaking of the irreducible building blocks, with the formation of transient surfactant dimers mediated by the counterions that promote the stabilization of cylindrical instead of spherical micelles. The present model is of general applicability and can be extended to all systems controlled by the presence of mobile charges.

Hybrid particle-field molecular dynamics under constant pressure.

Hybrid particle-field methods are computationally efficient approaches for modeling soft matter systems. So far, applications of these methodologies have been limited to constant volume conditions. Here, we reformulate particle-field interactions to represent systems coupled to constant external pressure. First, we show that the commonly used particle-field energy functional can be modified to model and parameterize the isotropic contributions to the pressure tensor without interfering with the microscopic forces on the particles. Second, we employ a square gradient particle-field interaction term to model non-isotropic contributions to the pressure tensor, such as in surface tension phenomena. This formulation is implemented within the hybrid particle-field molecular dynamics approach and is tested on a series of model systems. Simulations of a homogeneous water box demonstrate that it is possible to parameterize the equation of state to reproduce any target density for a given external pressure. Moreover, the same parameterization is transferable to systems of similar coarse-grained mapping resolution. Finally, we evaluate the feasibility of the proposed approach on coarse-grained models of phospholipids, finding that the term between water and the lipid hydrocarbon tails is alone sufficient to reproduce the experimental area per lipid in constant-pressure simulations and to produce a qualitatively correct lateral pressure profile.

Mesoscale Electrostatics Driving Particle Dynamics in Nonhomogeneous Dielectrics.

We introduce a density functional-based formalism to compute the electrostatic energy and forces for a mesoscopic system in the condensed phase, described with molecular resolution. The dielectric permittivity is variable in space, and it is dependent on the density fields of the individual particles present in the system. The electrostatic potential is obtained from standard numerical solutions of the generalized Poisson equation. The presence of a particle-dependent varying dielectrics produces the appearance of mesoscopic polarization forces, which are dependent on the local fluctuations of the permittivity, as well as of the electrostatic field. The proposed implementation is numerically robust, with an error on the Coulomb forces that can be systematically controlled by the mesh of spatial grid used for solving the generalized Poisson equation. We show that the method presented here is able to reproduce the concentration-dependent partitioning of an ideal salt in water/oil mixtures, in particular, reproducing the ∝ 1/ϵ dependency of the partition coefficient for the free ions predicted by Born theory. Moreover, this approach reproduces the correct electrostatic features of both dipolar and charged lipid bilayers, with positive membrane and dipole potentials. The sum of both Coulomb and polarization interactions inside the membrane yields a globally repulsive potential of mean force for the ions, independently on their charge. The computational efficiency of the method makes it particularly suitable for the description of large-scale polyelectrolyte soft-matter systems.

Hybrid Particle-Field Molecular Dynamics Simulations of Charged Amphiphiles in an Aqueous Environment.

We develop and test specific coarse-grained models for charged amphiphilic systems such as palmitoyloleoylphosphatidylglycerol (POPG) lipid bilayer and sodium dodecyl sulfate (SDS) surfactant in an aqueous environment, to verify the ability of the hybrid particle-field method to provide a realistic description of polyelectrolytes. According to the hybrid approach, the intramolecular interactions are treated by a standard molecular Hamiltonian, and the nonelectrostatic intermolecular forces are described by density fields. Electrostatics is introduced as an additional external field obtained by a modified particle-mesh Ewald procedure, as recently proposed [Zhu et al. Phys. Chem. Chem. Phys. 2016 , 18 , 9799 ]. Our results show that, upon proper calibration of key parameters, electrostatic forces can be correctly reproduced. Molecular dynamics simulations indicate that the methodology is robust with respect to the choice of the relative dielectric constant, yielding the same correct qualitative behavior for a broad range of values. In particular, our methodology reproduces well the organization of the POPG bilayer, as well as the SDS concentration-dependent change in the morphology of the micelles from spherical to microtubular aggregates. The inclusion of explicit electrostatics with good accuracy and low computational cost paves the way for a significant extension of the hybrid particle-field method to biological systems, where the polyelectrolyte component plays a fundamental role for both structural and dynamical molecular properties.

Self-Assembly of α-Tocopherol Transfer Protein Nanoparticles: A Patchy Protein Model.

We describe the mechanism of self-aggregation of α-tocopherol transfer protein into a spherical nanocage employing Monte Carlo simulations. The protein is modeled by a patchy coarse-grained representation, where the protein-protein interfaces, determined in the past by X-ray diffraction, are represented by simplified two-body interaction potentials. Our results show that the oligomerization kinetics proceeds in two steps, with the formation of metastable trimeric units and the subsequent assembly into the spherical aggregates. Data are in agreement with experimental observations regarding the prevalence of different aggregation states at specific ambient conditions. Finally, our results indicate a route for the experimental stabilization of the trimer, crucial for the understanding of the physiological role of such aggregates in vitamin E body trafficking.

Cholesteric and screw-like nematic phases in systems of helical particles.

Recent numerical simulations of hard helical particle systems unveiled the existence of a novel chiral nematic phase, termed screw-like, characterised by the helical organization of the particle C2 symmetry axes round the nematic director with periodicity equal to the particle pitch. This phase forms at high density and can follow a less dense uniform nematic phase, with relative occurrence of the two phases depending on the helix morphology. Since these numerical simulations were conducted under three-dimensional periodic boundary conditions, two questions could remain open. First, the real nature of the lower density nematic phase, expected to be cholesteric. Second, the influence that the latter, once allowed to form, may have on the existence and stability of the screw-like nematic phase. To address these questions, we have performed Monte Carlo and molecular dynamics numerical simulations of helical particle systems confined between two parallel repulsive walls. We have found that the removal of the periodicity constraint along one direction allows a relatively-long-pitch cholesteric phase to form, in lieu of the uniform nematic phase, with helical axis perpendicular to the walls while the existence and stability of the screw-like nematic phase are not appreciably affected by this change of boundary conditions.

Chiral self-assembly of helical particles.

The shape of the building blocks plays a crucial role in directing self-assembly towards desired architectures. Out of the many different shapes, the helix has a unique position. Helical structures are ubiquitous in nature and a helical shape is exhibited by the most important biopolymers like polynucleotides, polypeptides and polysaccharides as well as by cellular organelles like flagella. Helical particles can self-assemble into chiral superstructures, which may have a variety of applications, e.g. as photonic (meta)materials. However, a clear and definite understanding of these structures has not been entirely achieved yet. We have recently undertaken an extensive investigation on the phase behaviour of hard helical particles, using numerical simulations and classical density functional theory. Here we present a detailed study of the phase diagram of hard helices as a function of their morphology. This includes a variety of liquid-crystal phases, with different degrees of orientational and positional ordering. We show how, by tuning the helix parameters, it is possible to control the organization of the system. Starting from slender helices, whose phase behaviour is similar to that of rodlike particles, an increase in curliness leads to the onset of azimuthal correlations between the particles and the formation of phases specific to helices. These phases feature a new kind of screw order, of which there is experimental evidence in colloidal suspensions of helical flagella.

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.

Left or right cholesterics? A matter of helix handedness and curliness.

Using an Onsager-like theory, we have investigated the relationship between the morphology of hard helical particles and the features (pitch and handedness) of the cholesteric phase that they form. We show that right-handed helices can assemble into right- (R) and left-handed (L) cholesterics, depending on their curliness, and that the cholesteric pitch is a non-monotonic function of the intrinsic pitch of particles. The theory leads to the definition of a hierarchy of pseudoscalars, which quantify the difference in the average excluded volume between pair configurations of helices having (R) and (L)-skewed axes. The predictions of the Onsager-like theory are supported by Monte Carlo simulations of the isotropic phase of hard helices, showing how the cholesteric organization, which develops on scales longer than hundreds of molecular sizes, is encoded in the short-range chiral correlations between the helical axes.

Communication: from rods to helices: evidence of a screw-like nematic phase.

Evidence of a special chiral nematic phase is provided using numerical simulation and Onsager theory for systems of hard helical particles. This phase appears at the high density end of the nematic phase, when helices are well aligned, and is characterized by the C2 symmetry axes of the helices spiraling around the nematic director with periodicity equal to the particle pitch. This coupling between translational and rotational degrees of freedom allows a more efficient packing and hence an increase of translational entropy. Suitable order parameters and correlation functions are introduced to identify this screw-like phase, whose main features are then studied as a function of radius and pitch of the helical particles. Our study highlights the physical mechanism underlying a similar ordering observed in colloidal helical flagella [E. Barry, Z. Hensel, Z. Dogic, M. Shribak, and R. Oldenbourg, Phys. Rev. Lett. 96, 018305 (2006)] and raises the question of whether it could be observed in other helical particle systems, such as DNA, at sufficiently high densities.

The isotropic-to-nematic phase transition in hard helices: theory and simulation.

We investigate the isotropic-to-nematic phase transition in systems of hard helical particles, using Onsager theory and Monte Carlo computer simulations. Motivation of this work resides in the ubiquity of the helical shape motif in many natural and synthetic polymers, as well as in the well known importance that the details of size and shape have in determining the phase behaviour and properties of (soft) condensed matter systems. We discuss the differences with the corresponding spherocylinder phase diagram and find that the helix parameters affect the phase behaviour and the existence of the nematic phase. We find that for high helicity Onsager theory significantly departs from numerical simulations even when a modified form of the Parsons-Lee rescaling is included to account for the non-convexity of particles.