Self-Assembly of Ternary Cubic, Hexagonal, and Lamellar Mesophases Using the Lattice-Boltzmann Kinetic Method

pH- and temperature-responsive supramolecular assemblies with highly adjustable viscoelasticity: a multi-stimuli binary system

Stimuli-responsive materials are increasingly needed for the development of smart electronic, mechanical, and biological devices and systems relying on switchable, tunable, and adaptable properties. Herein, we report a novel pH- and temperature-responsive binary supramolecular assembly involving a long-chain hydroxyamino amide (HAA) and an inorganic hydrotrope, boric acid, with highly tunable viscous and viscoelastic properties. The system under investigation demonstrates a high degree of control over its viscosity, with the capacity to achieve over four orders of magnitude of control through the concomitant manipulation of pH and temperature. In addition, the transformation from non-Maxwellian to Maxwellian fluid behavior could also be induced by changing the pH and temperature. Switchable rheological properties were ascribed to the morphological transformation between spherical vesicles, aggregated/fused spherical vesicles, and bicontinuous gyroid structures revealed by cryo-TEM studies. The observed transitions are attributed to the modulation of the head group spacing between HAA molecules under different pH conditions. Specifically, acidic conditions induce electrostatic repulsion between the protonated amino head groups, leading to an increased spacing. Conversely, under basic conditions, the HAA head group spacing is reduced due to the intercalation of tetrahydroxyborate, facilitated by hydrogen bonding.

Mesoscale models of dispersions stabilized by surfactants and colloids

In this paper we discuss and give an outlook on numerical models describing dispersions, stabilized by surfactants and colloidal particles. Examples of these dispersions are foams and emulsions. In particular, we focus on the potential of the diffuse interface models based on a free energy approach, which describe dispersions with the surface-active agent soluble in one of the bulk phases. The free energy approach renders thermodynamic consistent models with realistic sorption isotherms and adsorption kinetics. The free energy approach is attractive because of its ability to describe highly complex dispersions, such as emulsions stabilized by ionic surfactants, or surfactant mixtures and dispersions with surfactant micelles. We have classified existing numerical methods into classes, using either a Eulerian or a Lagrangian representation for fluid and for the surfactant/colloid. A Eulerian representation gives a more coarse-grained, mean field description of the surface-active agent, while a Lagrangian representation can deal with steric effects and larger complexity concerning geometry and (amphiphilic) wetting properties of colloids and surfactants. However, the similarity between the description of wetting properties of both Eulerian and Lagrangian models allows for the development of hybrid Eulerian/Lagrangian models having advantages of both representations.

High-order thermal lattice Boltzmann models derived by means of Gauss quadrature in the spherical coordinate system

We use the spherical coordinate system in the momentum space and an appropriate discretization procedure to derive a hierarchy of lattice Boltzmann (LB) models with variable temperature. The separation of the integrals in the momentum space into angular and radial parts allows us to compute the moments of the equilibrium distribution function by means of Gauss-Legendre and Gauss-Laguerre quadratures, as well as to find the elements of the discrete momentum set for each LB model in the hierarchy. The capability of the high-order models in this hierarchy to capture specific effects in microfluidics is investigated through a computer simulation of Couette flow by using the Shakhov collision term to get the right value of the Prandtl number.

Small-Angle X-ray Scattering Investigations of Biomolecular Confinement, Loading, and Release from Liquid-Crystalline Nanochannel Assemblies

This Perspective explores the recent progress made by means of small-angle scattering methods in structural studies of phase transitions in amphiphilic liquid-crystalline systems with nanochannel architectures and outlines some future directions in the area of hierarchically organized and stimuli-responsive nanochanneled assemblies involving biomolecules. Time-resolved small-angle X-ray scattering investigations using synchrotron radiation enable monitoring of the structural dynamics, the modulation of the nanochannel hydration, as well as the key changes in the soft matter liquid-crystalline organization upon stimuli-induced phase transitions. They permit establishing of the inner nanostructure transformation kinetics and determination of the precise sizes of the hydrophobic membraneous compartments and the aqueous channel diameters in self-assembled network architectures. Time-resolved structural studies accelerate novel biomedical, pharmaceutical, and nanotechnology applications of nanochannel soft materials by providing better control of DNA, peptide and protein nanoconfinement, and release from diverse stimuli-responsive nanocarrier systems.

Phenyl-ring-bearing cationic surfactants: effect of ring location on the micellar structure

A series of isomeric cationic surfactants (S1-S5) bearing a long alkyl chain that carries a 1,4-phenylene unit and a trimethyl ammonium headgroup was synthesized; the location of the phenyl ring within the alkyl tail was varied in an effort to understand its influence on the amphiphilic properties of the surfactants. The cmc's of the surfactants were estimated using ionic conductivity measurements and isothermal calorimetric titrations (ITC); the values obtained by the two methods were found to be in excellent agreement. The ITC measurements provided additional insight into the various thermodynamic parameters associated with the micellization process. Although all five surfactants have exactly the same molecular formula, their micellar properties were seen to vary dramatically depending on the location of the phenyl ring; the cmc was seen to decrease by almost an order of magnitude when the phenyl ring was moved from the tail end (cmc of S1 is 23 mM) to the headgroup region (cmc of S5 is 3 mM). In all cases, the enthalpy of micellization was negative but the entropy of micellization was positive, suggesting that in all of these systems the formation of micelles is both enthalpically and entropically favored. As expected, the decrease in cmc values upon moving the phenyl ring from the tail end to the headgroup region is accompanied by an increase in the thermodynamic driving force (ΔG) for micellization. To understand further the differences in the micellar structure of these surfactants, small-angle neutron scattering (SANS) measurements were carried out; these measurements reveal that the aggregation number of the micelles increases as the cmc decreases. This increase in the aggregation number is also accompanied by an increase in the asphericity of the micellar aggregate and a decrease in the fractional charge. Geometric packing arguments are presented to account for these changes in aggregation behavior as a function of phenyl ring location.

Petascale lattice-Boltzmann studies of amphiphilic cubic liquid crystalline materials in a globally distributed high-performance computing and visualization environment

We present very large-scale rheological studies of self-assembled cubic gyroid liquid crystalline phases in ternary mixtures of oil, water and amphiphilic species performed on petascale supercomputers using the lattice-Boltzmann method. These nanomaterials have found diverse applications in materials science and biotechnology, for example, in photovoltaic devices and protein crystallization. They are increasingly gaining importance as delivery vehicles for active agents in pharmaceuticals, personal care products and food technology. In many of these applications, the self-assembled structures are subject to flows of varying strengths and we endeavour to understand their rheological response with the objective of eventually predicting it under given flow conditions. Computationally, our lattice-Boltzmann simulations of ternary fluids are inherently memory- and data-intensive. Furthermore, our interest in dynamical processes necessitates remote visualization and analysis as well as the associated transfer and storage of terabytes of time-dependent data. These simulations are distributed on a high-performance grid infrastructure using the application hosting environment; we employ a novel parallel in situ visualization approach which is particularly suited for such computations on petascale resources. We present computational and I/O performance benchmarks of our application on three different petascale systems.

Derivation and thermodynamics of a lattice Boltzmann model with soluble amphiphilic surfactant

We derive a lattice Boltzmann model of a ternary fluid mixture, one component of which consists of amphiphilic molecules with a rotational degree of freedom. The model was first introduced elsewhere where it was constructed on the basis of the earlier Shan-Chen type lattice Boltzmann model for binary fluids. We provide a rigorous derivation of the model from an underlying continuum kinetic theory. In particular, we show how the model can be interpreted as a discretization of a Vlasov-Boltzmann type kinetic theory for a fluid composed of two species of oppositely charged monopoles and a composite, dipolar molecule. We also derive a free-energy functional for the model, including the contribution from the amphiphiles, and confirm that the free energy is lowered during phase separation and during the formation of a stable emulsion.

Real science at the petascale

We describe computational science research that uses petascale resources to achieve scientific results at unprecedented scales and resolution. The applications span a wide range of domains, from investigation of fundamental problems in turbulence through computational materials science research to biomedical applications at the forefront of HIV/AIDS research and cerebrovascular haemodynamics. This work was mainly performed on the US TeraGrid 'petascale' resource, Ranger, at Texas Advanced Computing Center, in the first half of 2008 when it was the largest computing system in the world available for open scientific research. We have sought to use this petascale supercomputer optimally across application domains and scales, exploiting the excellent parallel scaling performance found on up to at least 32 768 cores for certain of our codes in the so-called 'capability computing' category as well as high-throughput intermediate-scale jobs for ensemble simulations in the 32-512 core range. Furthermore, this activity provides evidence that conventional parallel programming with MPI should be successful at the petascale in the short to medium term. We also report on the parallel performance of some of our codes on up to 65 636 cores on the IBM Blue Gene/P system at the Argonne Leadership Computing Facility, which has recently been named the fastest supercomputer in the world for open science.

Rheological response and dynamics of the amphiphilic diamond phase from kinetic lattice–Boltzmann simulations

The purpose of the present paper is to report on the first computational study of the dynamical and rheological response of a self-assembled diamond mesophase under Couette flow in a ternary mixture composed of oil, water and an amphiphilic species. The amphiphilic diamond mesophase arises in a wide range of chemical and biological systems, and a knowledge of its rheological response has important implications in materials science and biotechnological applications. The simulations reported here are performed using a kinetic lattice–Boltzmann method. Lyotropic liquid crystals exhibit characteristic rheological responses in experiments that include shear-banding and a non-Newtonian flow curve as well as viscoelasticity under oscillatory shear. Their behaviour under steady and oscillatory shear is correctly reproduced in our simulations. On cessation of shear, as the morphology returns to the diamond phase, the relaxation of the stress response follows a stretched-exponential form for low initial strain rates.

Effect of tail architecture on self-assembly of amphiphiles for polymeric micelles

Brownian dynamics simulations were carried out to explore the self-assembly of amphiphilic copolymers composed of a linear hydrophilic head and a hydrophobic tail with different architectures. In order to investigate the effect of architecture of hydrophobic tail on self-assembling behavior, these architectures of linear, branched, starlike, and dendritic tails were selected for comparison, and the branching parameter of the tail was employed to characterize the tail architectures. The critical micelle concentration (cmc), dynamics of aggregation, aggregate distribution, gyration radius distribution, density profiles of micelle, shape anisotropy, and thermal stability were examined for the four typical types of copolymers. The calculated results reveal that the self-assembly of linear tail copolymer has the lowest cmc, and the consequently formed polymeric micelles have narrow dispersion and greater aggregate size, and the micelle is closer to spherical shape. It was found that the cmc is inversely proportional to the branching parameter. Linear tail aggregates in solution to form polymeric micelles with higher physical stability, compared to other architectures of tail. The size of polymeric micelle increases with the increase of the branching parameter of the tail, and it exhibits an exponential relationship with the branching parameter. In addition, the micelles formed from copolymers with a high branching parameter of the tail were found to have higher thermal stability. This work provides useful information on designing self-assembling systems for preparing polymeric micelles applied to drug delivery.