Hybrid Lattice Boltzmann/Dynamic Self-Consistent Field Simulations of Microphase Separation and Vesicle Formation in Block Copolymer Systems

Polymersomes in Drug Delivery─From Experiment to Computational Modeling

Polymersomes, composed of amphiphilic block copolymers, are self-assembled vesicles that have gained attention as potential drug delivery systems due to their good biocompatibility, stability, and versatility. Various experimental techniques have been employed to characterize the self-assembly behaviors and properties of polymersomes. However, they have limitations in revealing molecular details and underlying mechanisms. Computational modeling techniques have emerged as powerful tools to complement experimental studies and enabled researchers to examine drug delivery mechanisms at molecular resolution. This review aims to provide a comprehensive overview of the state of the art in the field of polymersome-based drug delivery systems, with an emphasis on insights gained from both experimental and computational studies. Specifically, we focus on polymersome morphologies, self-assembly kinetics, fusion and fission, behaviors in flow, as well as drug encapsulation and release mechanisms. Furthermore, we also identify existing challenges and limitations in this rapidly evolving field and suggest possible directions for future research.

Hybrid Molecules Consisting of Lysine Dendrons with Several Hydrophobic Tails: A SCF Study of Self-Assembling

In this article, we used the numerical self-consistent field method of Scheutjens-Fleer to study the micellization of hybrid molecules consisting of one polylysine dendron with charged end groups and several linear hydrophobic tails attached to its root. The main attention was paid to spherical micelles and the determination of the range of parameters at which they can appear. A relationship has been established between the size and internal structure of the resulting spherical micelles and the length and number of hydrophobic tails, as well as the number of dendron generations. It is shown that the splitting of the same number of hydrophobic monomers from one long tail into several short tails leads to a decrease in the aggregation number and, accordingly, the number of terminal charges in micelles. At the same time, it was shown that the surface area per dendron does not depend on the number of hydrophobic monomers or tails in the hybrid molecule. The relationship between the structure of hybrid molecules and the electrostatic properties of the resulting micelles has also been studied. It is found that the charge distribution in the corona depends on the number of dendron generations G in the hybrid molecule. For a small number of generations (up to G=3), a standard double electric layer is observed. For a larger number of generations (G=4), the charges of dendrons in the corona are divided into two populations: in the first population, the charges are in the spherical layer near the boundary between the micelle core and shell, and in the second population, the charges are near the periphery of the spherical shell. As a result, a part of the counterions is localized in the wide region between them. These results are of potential interest for the use of spherical dendromicelles as nanocontainers for drug delivery.

Using reactive dissipative particle dynamics to understand local shape manipulation of polymer vesicles

Biological cells have long been of interest to researchers due to their capacity to actively control their shape. Accordingly, there is significant interest in generating simplified synthetic protocells that can alter their shape based on an externally or internally generated stimulus. To date, most progress has been made towards controlling the global shape of a protocell, whereas less is known about generating a local shape change. Here, we seek to better understand the possible mechanisms for producing local morphological changes in a popular protocell system, the block copolymer vesicle. Accordingly, we have combined Dissipative Particle Dynamics (DPD) and the Split Reactive Brownian Dynamics algorithm (SRBD) to produce a simulation tool that is capable of modeling the dynamics of self-assembled polymer structures as they undergo chemical reactions. Using this Reactive DPD or RDPD method, we investigate local morphological change driven by either the microinjection of a stimulus or an enzymatically-produced stimulus. We find that sub-vesicle-scale morphological change can be induced by either a solvent stimulus that swells the vesicle membrane, or by a reactant stimulus that alters the chemistry of the block polymer in the membrane corona. Notably, the latter method results in a more persistent local deformation than the former, which we attribute to the slower diffusion of polymer chains relative to the solvent. We quantify this deformation and show that it can be modulated by altering the interaction parameter of the parts of the polymer chain that are affected by the stimulus.

Dynamic Self-Consistent Field Approach for Studying Kinetic Processes in Multiblock Copolymer Melts

The self-consistent field theory is a popular and highly successful theoretical framework for studying equilibrium (co)polymer systems at the mesoscopic level. Dynamic density functionals allow one to use this framework for studying dynamical processes in the diffusive, non-inertial regime. The central quantity in these approaches is the mobility function, which describes the effect of chain connectivity on the nonlocal response of monomers to thermodynamic driving fields. In a recent study, one of us and coworkers have developed a method to systematically construct mobility functions from reference fine-grained simulations. Here we focus on melts of linear chains in the Rouse regime and show how the mobility functions can be calculated semi-analytically for multiblock copolymers with arbitrary sequences without resorting to simulations. In this context, an accurate approximate expression for the single-chain dynamic structure factor is derived. Several limiting regimes are discussed. Then we apply the resulting density functional theory to study ordering processes in a two-length scale block copolymer system after instantaneous quenches into the ordered phase. Different dynamical regimes in the ordering process are identified: at early times, the ordering on short scales dominates; at late times, the ordering on larger scales takes over. For large quench depths, the system does not necessarily relax into the true equilibrium state. Our density functional approach could be used for the computer-assisted design of quenching protocols in order to create novel nonequilibrium materials.

Studying polymer solutions with particle-based models linked to classical density functionals: co-non-solvency

We demonstrate the potential of hybrid particle-based models, where interactions are introduced through functionals of local order parameters, in describing multicomponent polymer solutions. The link to a free-energy-like functional is advantageous for controlling the thermodynamics of the model. We focus on co-non-solvency - the collapse of polymer chains in dilute mixtures with two miscible good solvents, having different affinities towards the polymer. We employ a simple model where polymers and solvents are represented, respectively, by worm-like chains and single particles. Non-bonded interactions are captured by a polynomial which is third order in local densities and can, therefore, describe liquid-vapour coexistence. The parameterisation of the functional benefits from an elementary mean-field approximation to the statistical mechanics of the model. The model provides a framework for Monte Carlo simulations using a particle-to-mesh algorithm. Studies with conventional generic bead-spring and all-atom models have demonstrated that co-non-solvency is caused by preferential binding of the better solvent (termed cosolvent) with polymer. Hence, segmental loops bridged by cosolvent molecules are formed, initiating polymer collapse. The mesoscopic hybrid model differs conceptually from the conventional microscopic descriptions. Yet, it reproduces the same co-non-solvency mechanism supporting its universality. Films of adsorbed ternary solutions, showing co-non-solvency in the dilute regime, are considered at high concentrations. In this case, chains do not collapse. The properties of loops and tails of the adsorbed polymer agree with early theoretical predictions obtained for concentrated binary solutions.

Self-Assembly of Lysine-Based Dendritic Surfactants Modeled by the Self-Consistent Field Approach

Implementing a united atom model, we apply self-consistent field theory to study structure and thermodynamic properties of spherical micelles composed of surfactants that combine an alkyl tail with a charged lysine-based dendritic headgroup. Following experiments, the focus was on dendron surfactants with varying tail length and dendron generations G0, G1, G2. The heads are subject to acetylation modification which reduces the charge and hydrophilicity. We establish a reasonable parameter set which results in semiquantitative agreement with the available experiments. The critical micellization concentration, aggregation number, and micelle size are discussed. The strongly charged dendronic surfactants micelles are stable for generation numbers G0 and G1, for progressively higher ionic strengths. Associates of G2 surfactants are very small and can only be found at extreme surfactant concentration and salt strengths. Micelles of corresponding weaker charged acetylated variants exist up to G2, tolerate significantly lower salt concentrations, but lose the spherical micelle topology for G0 at high ionic strengths.

Well-ordered self-assembled nanostructures of block copolymer films via synergistic integration of chemoepitaxy and zone annealing

The integrated chemical template/zone annealing method has the capability to rapidly fabricate well-aligned and well-oriented nanostructures over a macroscopic area.

Hybrid particle-continuum simulations coupling Brownian dynamics and local dynamic density functional theory

We present a multiscale hybrid particle-field scheme for the simulation of relaxation and diffusion behavior of soft condensed matter systems. It combines particle-based Brownian dynamics and field-based local dynamics in an adaptive sense such that particles can switch their level of resolution on the fly. The switching of resolution is controlled by a tuning function which can be chosen at will according to the geometry of the system. As an application, the hybrid scheme is used to study the kinetics of interfacial broadening of a polymer blend, and is validated by comparing the results to the predictions from pure Brownian dynamics and pure local dynamics calculations.

Structure of electric double layers in capacitive systems and to what extent (classical) density functional theory describes it

Ongoing scientific interest is aimed at the properties and structure of electric double layers (EDLs), which are crucial for capacitive energy storage, water treatment, and energy harvesting technologies like supercapacitors, desalination devices, blue engines, and thermocapacitive heat-to-current converters. A promising tool to describe their physics on a microscopic level is (classical) density functional theory (DFT), which can be applied in order to analyze pair correlations and charge ordering in the primitive model of charged hard spheres. This simple model captures the main properties of ionic liquids and solutions and it predicts many of the phenomena that occur in EDLs. The latter often lead to anomalous response in the differential capacitance of EDLs. This work constructively reviews the powerful theoretical framework of DFT and its recent developments regarding the description of EDLs. It explains to what extent current approaches in DFT describe structural ordering and in-plane transitions in EDLs, which occur when the corresponding electrodes are charged. Further, the review briefly summarizes the history of modeling EDLs, presents applications, and points out limitations and strengths in present theoretical approaches. It concludes that DFT as a sophisticated microscopic theory for ionic systems is expecting a challenging but promising future in both fundamental research and applications in supercapacitive technologies.

Combining cell-based hydrodynamics with hybrid particle-field simulations: efficient and realistic simulation of structuring dynamics

We have extended an existing hybrid MD-SCF simulation technique that employs a coarsening step to enhance the computational efficiency of evaluating non-bonded particle interactions. This technique is conceptually equivalent to the single chain in mean-field (SCMF) method in polymer physics, in the sense that non-bonded interactions are derived from the non-ideal chemical potential in self-consistent field (SCF) theory, after a particle-to-field projection. In contrast to SCMF, however, MD-SCF evolves particle coordinates by the usual Newton's equation of motion. Since collisions are seriously affected by the softening of non-bonded interactions that originates from their evaluation at the coarser continuum level, we have devised a way to reinsert the effect of collisions on the structural evolution. Merging MD-SCF with multi-particle collision dynamics (MPCD), we mimic particle collisions at the level of computational cells and at the same time properly account for the momentum transfer that is important for a realistic system evolution. The resulting hybrid MD-SCF/MPCD method was validated for a particular coarse-grained model of phospholipids in aqueous solution, against reference full-particle simulations and the original MD-SCF model. We additionally implemented and tested an alternative and more isotropic finite difference gradient. Our results show that efficiency is improved by merging MD-SCF with MPCD, as properly accounting for hydrodynamic interactions considerably speeds up the phase separation dynamics, with negligible additional computational costs compared to efficient MD-SCF. This new method enables realistic simulations of large-scale systems that are needed to investigate the applications of self-assembled structures of lipids in nanotechnologies.

Multiscale Molecular Simulations of Polymer-Matrix Nanocomposites: or What Molecular Simulations Have Taught us About the Fascinating Nanoworld

Following the substantial progress in molecular simulations of polymer-matrix nanocomposites, now is the time to reconsider this topic from a critical point of view. A comprehensive survey is reported herein providing an overview of classical molecular simulations, reviewing their major achievements in modeling polymer matrix nanocomposites, and identifying several open challenges. Molecular simulations at multiple length and time scales, working hand-in-hand with sensitive experiments, have enhanced our understanding of how nanofillers alter the structure, dynamics, thermodynamics, rheology and mechanical properties of the surrounding polymer matrices.

Ordering kinetics of lamella-forming block copolymers under the guidance of various external fields studied by dynamic self-consistent field theory

We theoretically engineer a new scheme, which integrates a permanent field for pattern registration and a dynamic external field for defect annihilation, to direct the self-assembly of block copolymers.

Insights into ordered microstructures and ordering mechanisms of ABC star terpolymers by integrating dynamic self-consistent field theory and variable cell shape methods

A theoretical approach coupling dynamic self-consistent field (SCF) theory for inhomogeneous polymeric fluids and variable cell shape (VCS) method for automatically adjusting cell shape and size is developed to investigate ordered microstructures and the ordering mechanisms of block copolymer melts. Using this simulation method, we first re-examined the microphase separation of the simplest AB diblock copolymers, and tested the validity and efficiency of the novel method by comparing the results with those obtained from the dynamic SCF theory. An appropriate relaxation parameter of the VCS method effectively accelerates the system towards a stable morphology without distortions or defects. The dynamic SCF/VCS method is then applied to identify the richness morphologies of ABC star terpolymers and explore the ordering mechanisms of star terpolymer melts quenched from homogenous states. A diverse range of ordered microstructures, including two-dimensional tiling patterns, hierarchical structures and ordinary microstructures, are predicted. Three types of ordering mechanisms, namely, one-step, quick-slow and step-wise procedures, are discovered in the disorder-to-order transition of ABC star terpolymers. The procedures of microphase separation in the ABC star terpolymer melts are remarkably affected by the composition of star terpolymers and the strength of interaction parameters.

Harnessing Anisotropic Nanoposts to Enhance Long-Range Orientation Order of Directed Self-Assembly Nanostructures via Large Cell Simulations

Self-assembly behaviors of cylinder-forming diblock copolymers directed by an array of anisotropic nanoposts with elliptical shape are explored by large cell simulations of self-consistent field theory. The strategy using elliptical nanoposts allows us to generate long-range order cylinders with single orientation by suppressing other selections of cylinder alignment. The anisotropy of nanoposts plays a significant role in improving the tolerance of commensurability conditions between the dimensions of nanopost lattices and the period of cylinders. Moreover, the local defect structures could be regulated through varying the spacing and orientation of elliptical nanoposts. This work may provide useful guidelines for designing the topographical templates and lay the groundwork for fabricating well-ordered nanostructures of block copolymer lithography.

Self-consistent field approach for cross-linked copolymer materials

A generalized self-consistent field approach for polymer networks with a fixed topology is developed. It is shown that the theory reproduces the localization of cross-links, which is characteristic for gels. The theory is then used to study the order-disorder transition in regular networks of end-linked diblock copolymers. Compared to diblock copolymer melts, the transition is shifted towards lower values of the incompatibility parameter χ (the Flory-Huggins parameter). Moreover, the transition becomes strongly first order already at the mean-field level. If stress is applied, the transition is further shifted and finally vanishes in a critical point.