Non-Equilibrium Dynamics of Vesicles and Micelles by Self-Assembly of Block Copolymers with Double Thermoresponsivity

A Versatile Approach to Stabilize Liquid-Liquid Interfaces using Surfactant Self-Assembly

Stabilizing liquid-liquid interfaces, whether between miscible or immiscible liquids, is crucial for a wide range of applications, including energy storage, microreactors, and biomimetic structures. In this study, a versatile approach for stabilizing the water-oil interface is presented using the morphological transitions that occur during the self-assembly of anionic, cationic, and nonionic surfactants mixed with fatty acid oils. The morphological transitions underlying this approach are characterized and extensively studied through small-angle X-ray scattering (SAXS), rheometry, and microscopy techniques. Dissipative particle dynamics (DPD) as a simulation tool is adopted to investigate these morphological transitions both in the equilibrium ternary system as well as in the dynamic condition of the water-oil interface. Such a versatile strategy holds promise for enhancing applications such as liquid-in-liquid 3D printing. Moreover, it has the potential to revolutionize a wide range of fields where stabilizing liquid-liquid interfaces not only offers unprecedented opportunities for fine-tuning nanostructural morphologies but also imparts interesting practical features to the resulting liquid shapes. These features include perfusion capabilities, self-healing, and porosity, which could have significant implications for various industries.

Self-assembly of sequence-regulated amphiphilic copolymers with alternating rod and coil pendants

We conducted a computational study on the self-assembly behavior of sequence-controlled amphiphilic copolymers with alternating rod and coil pendants. Complex self-assembled morphologies, such as onion-like vesicles with two layers, can be generated by introducing rod pendants. The amphiphilic alternating copolymers self-assemble into onion-like vesicles through a fusion process of tiny micelles and a bending operation of disk-like micelles with double layers. A stimuli-responsive simulation shows that the cylindrical vesicles can transform into onion-like vesicles by a rod-to-coil conformation transition of rigid pendants. Inspired by this finding, we conducted a drug-loading simulation by adding two reactive drugs at different stages and found that the onion-like vesicles can almost completely isolate two drugs. This work provides theoretical guidance on the self-assembly of amphiphilic alternating copolymers with rod and coil pendants for future experimental design.

Engulfing Behavior of Nanoparticles into Thermoresponsive Microgels: A Mesoscopic Simulation Study

The engulfing of nanoparticles into microgels provides a versatile platform to design nano- and microstructured materials with various shape anisotropies and multifunctional properties. Manipulating the spontaneous engulfment process remains elusive. Herein, we report a mesoscopic simulation study on the engulfing behavior of nanoparticles into thermoresponsive microgels. The effects of the multiple parameters, including binding strength, temperature, and nanoparticle size, are examined systematically. Our simulation results disclose three engulfing states at different temperatures, namely full-engulfing, half-engulfing, and surface contact. The engulfing depth is determined by the complementary balance of interfacial elastocapillarity. Specifically, the van der Waals interaction of hybrid microgel-nanoparticle offers the capillary force while the internally networked structure of microgel reinforces the elasticity repulsion. Our study, validated by relevant experimental results, provides a mechanistic understanding of the interfacial elastocapillarity for nanoparticle-microgels.

Rod-coil block copolymer aggregates via polymerization-induced self-assembly

Polymerization-induced self-assembly (PISA), incorporating the polymerization with in situ self-assembly, can achieve nano-objects efficiently. However, the cooperative polymerization and self-assembly lead to unclear polymerization kinetics and aggregation behavior, especially for the systems forming rigid chains. Here, we used dissipative particle dynamics simulations with a probability-based reaction model to explore the PISA behavior of rod-coil block copolymer systems. The impact of the length of macromolecular initiators, the targeted length of rigid chains, and the reaction probability on the PISA behavior, including polymerization kinetics and self-assembly, were examined. The difference between PISA and traditional self-assembly was revealed. A comparison with experimental observations shows that the simulation can capture the essential feature of the PISA. The present work provides a comprehensive understanding of rod-coil PISA systems and may provide meaningful information for future experimental research.

Self-assembly of amphiphilic alternating copolymers with stimuli-responsive rigid pendant groups

Amphiphilic alternating copolymers (AACs) possess unique self-assembly behaviours owing to their unique regular architecture.

Implicit-solvent coarse-grained modeling for polymer solutions via Mori-Zwanzig formalism

We present a bottom-up coarse-graining (CG) method to establish implicit-solvent CG modeling for polymers in solution, which conserves the dynamic properties of the reference microscopic system. In particular, tens to hundreds of bonded polymer atoms (or Lennard-Jones beads) are coarse-grained as one CG particle, and the solvent degrees of freedom are eliminated. The dynamics of the CG system is governed by the generalized Langevin equation (GLE) derived via the Mori-Zwanzig formalism, by which the CG variables can be directly and rigorously linked to the microscopic dynamics generated by molecular dynamics (MD) simulations. The solvent-mediated dynamics of polymers is modeled by the non-Markovian stochastic dynamics in GLE, where the memory kernel can be computed from the MD trajectories. To circumvent the difficulty in direct evaluation of the memory term and generation of colored noise, we exploit the equivalence between the non-Markovian dynamics and Markovian dynamics in an extended space. To this end, the CG system is supplemented with auxiliary variables that are coupled linearly to the momentum and among themselves, subject to uncorrelated Gaussian white noise. A high-order time-integration scheme is used to solve the extended dynamics to further accelerate the CG simulations. To assess, validate, and demonstrate the established implicit-solvent CG modeling, we have applied it to study four different types of polymers in solution. The dynamic properties of polymers characterized by the velocity autocorrelation function, diffusion coefficient, and mean square displacement as functions of time are evaluated in both CG and MD simulations. Results show that the extended dynamics with auxiliary variables can construct arbitrarily high-order CG models to reproduce dynamic properties of the reference microscopic system and to characterize long-time dynamics of polymers in solution.

Switch It Inside-Out: "Schizophrenic" Behavior of All Thermoresponsive UCST-LCST Diblock Copolymers

This feature article reviews our recent advancements on the synthesis, phase behavior, and micellar structures of diblock copolymers consisting of oppositely thermoresponsive blocks in aqueous environments. These copolymers combine a nonionic block, which shows lower critical solution temperature (LCST) behavior, with a zwitterionic block that exhibits an upper critical solution temperature (UCST). The transition temperature of the latter class of polymers is strongly controlled by its molar mass and by the salt concentration, in contrast to the rather invariant transition of nonionic polymers with type II LCST behavior such as poly(N-isopropylacrylamide) or poly(N-isopropyl methacrylamide). This allows for implementing the sequence of the UCST and LCST transitions of the polymers at will by adjusting either molecular or, alternatively, physical parameters. Depending on the location of the transition temperatures of both blocks, different switching scenarios are realized from micelles to inverse micelles, namely via the molecularly dissolved state, the aggregated state, or directly. In addition to studies of (semi)dilute aqueous solutions, highly concentrated systems have also been explored, namely water-swollen thin films. Concerning applications, we discuss the possible use of the diblock copolymers as "smart" nanocarriers.

All-In-One "Schizophrenic" Self-Assembly of Orthogonally Tuned Thermoresponsive Diblock Copolymers

Smart, fully orthogonal switching was realized in a highly biocompatible diblock copolymer system with variable trigger-induced aqueous self-assembly. The polymers are composed of nonionic and zwitterionic blocks featuring lower and upper critical solution temperatures (LCSTs and UCSTs). In the system investigated, diblock copolymers from poly( N-isopropyl methacrylamide) (PNIPMAM) and a poly(sulfobetaine methacrylamide), systematic variation of the molar mass of the latter block allowed for shifting the UCST of the latter above the LCST of the PNIPMAM block in a salt-free condition. Thus, successive thermal switching results in "schizophrenic" micellization, in which the roles of the hydrophobic core block and the hydrophilic shell block are interchanged depending on the temperature. Furthermore, by virtue of the strong electrolyte-sensitivity of the zwitterionic polysulfobetaine block, we succeeded to shift its UCST below the LCST of the PNIPMAM block by adding small amounts of an electrolyte, thus inverting the pathway of switching. This superimposed orthogonal switching by electrolyte addition enabled us to control the switching scenarios between the two types of micelles (i) via an insoluble state, if the LCST-type cloud point is below the UCST-type cloud point, which is the case at low salt concentrations or (ii) via a molecularly dissolved state, if the LCST-type cloud point is above the UCST-type cloud point, which is the case at high salt concentrations. Systematic variation of the block lengths allowed for verifying the anticipated behavior and identifying the molecular architecture needed. The versatile and tunable self-assembly offers manifold opportunities, for example, for smart emulsifiers or for sophisticated carrier systems.

Effect of temperature and hydrophilic ratio on the structure of poly(N-vinylcaprolactam)-block-poly(dimethylsiloxane)-block-poly(N-vinylcaprolactam) polymersomes

Nanosized polymeric vesicles (polymersomes) assembled from ABA triblock copolymers of poly(N-vinylcaprolactam)-poly(dimethylsiloxane)-poly(N-vinylcaprolactam) (PVCL-PDMS-PVCL) are a promising platform for biomedical applications, as the temperature-responsiveness of the PVCL blocks enables reversible vesicle shrinkage and permeability of the polymersome shell at elevated temperatures. Herein, we explore the effects of molecular weight, polymer block weight ratios, and temperature on the structure of these polymersomes via electron microscopy, dynamic light scattering, small angle neutron scattering (SANS), and all-atom molecular dynamic methods. We show that the shell structure and overall size of the polymersome can be tuned by varying the hydrophilic (PVCL) weight fraction of the polymer: at room temperature, polymers of smaller hydrophilic ratios form larger vesicles that have thinner shells, whereas polymers with higher PVCL content exhibit interchain aggregation of PVCL blocks within the polymersome shell above 50 °C. Model fitting and model-free analysis of the SANS data reveals that increasing the mass ratio of PVCL to the total copolymer weight from 0.3 to 0.56 reduces the temperature-induced change in vesicle diameter by a factor of 3 while simultaneously increasing the change in shell thickness by a factor of 1.5. Finally, by analysis of the shell structures and overall size of polymersomes with various PVCL weight ratios and those without temperature-dependent polymer components, we bring into focus the mechanism of temperature-triggered drug release reported in a previous study. This work provides new fundamental perspectives on temperature-responsive polymersomes and elucidates important structure-property relationships of their constituent polymers.

Dynamic principle for ensemble control tools

Dynamical equations describing physical systems in contact with a thermal bath are commonly extended by mathematical tools called "thermostats." These tools are designed for sampling ensembles in statistical mechanics. Here we propose a dynamic principle underlying a range of thermostats which is derived using fundamental laws of statistical physics and ensures invariance of the canonical measure. The principle covers both stochastic and deterministic thermostat schemes. Our method has a clear advantage over a range of proposed and widely used thermostat schemes that are based on formal mathematical reasoning. Following the derivation of the proposed principle, we show its generality and illustrate its applications including design of temperature control tools that differ from the Nosé-Hoover-Langevin scheme.

Stimulated Transitions of Directed Nonequilibrium Self-Assemblies

Near-equilibrium stimulus-responsive polymers have been used extensively to introduce morphological variations in dependence of adaptable conditions. Far-less-well studied are triggered transformations at constant conditions. These require the involvement of metastable states, which are either able to approach the equilibrium state after deviation from metastability or can be frozen on returning from nonequilibrium to equilibrium. Such functional nonequilibrium macromolecular systems hold great promise for on-demand transformations, which result in substantial changes in their material properties, as seen for triggered gelations. Herein, a diblock copolymer system consisting of a hydrophilic block and a block that is responsive to both pressure and temperature, is introduced. This species demonstrates various micellar transformations upon leaving equilibrium/nonequilibrium states, which are triggered by a temperature deflection or a temporary application of hydrostatic pressure.

Refined energy-conserving dissipative particle dynamics model with temperature-dependent properties and its application in solidification problem

It has been observed previously that the physical behaviors of Schmidt number (Sc) and Prandtl number (Pr) of an energy-conserving dissipative particle dynamics (eDPD) fluid can be reproduced by the temperature-dependent weight function appearing in the dissipative force term. In this paper, we proposed a simple and systematic method to develop the temperature-dependent weight function in order to better reproduce the physical fluid properties. The method was then used to study a variety of phase-change problems involving solidification. The concept of the "mushy" eDPD particle was introduced in order to better capture the temperature profile in the vicinity of the solid-liquid interface, particularly for the case involving high thermal conductivity ratio. Meanwhile, a way to implement the constant temperature boundary condition at the wall was presented. The numerical solutions of one- and two-dimensional solidification problems were then compared with the analytical solutions and/or experimental results and the agreements were promising.

GPU-accelerated Red Blood Cells Simulations with Transport Dissipative Particle Dynamics

Mesoscopic numerical simulations provide a unique approach for the quantification of the chemical influences on red blood cell functionalities. The transport Dissipative Particles Dynamics (tDPD) method can lead to such effective multiscale simulations due to its ability to simultaneously capture mesoscopic advection, diffusion, and reaction. In this paper, we present a GPU-accelerated red blood cell simulation package based on a tDPD adaptation of our red blood cell model, which can correctly recover the cell membrane viscosity, elasticity, bending stiffness, and cross-membrane chemical transport. The package essentially processes all computational workloads in parallel by GPU, and it incorporates multi-stream scheduling and non-blocking MPI communications to improve inter-node scalability. Our code is validated for accuracy and compared against the CPU counterpart for speed. Strong scaling and weak scaling are also presented to characterizes scalability. We observe a speedup of 10.1 on one GPU over all 16 cores within a single node, and a weak scaling efficiency of 91% across 256 nodes. The program enables quick-turnaround and high-throughput numerical simulations for investigating chemical-driven red blood cell phenomena and disorders.