Ultrasmall Single-Chain Nanoparticles Derived from Amphiphilic Alternating Copolymers.

The collapse or folding of an individual polymer chain into a nanoscale particle gives rise to single-chain nanoparticles (SCNPs), which share a soft nature with biological protein particles. The precise control of their properties, including morphology, internal structure, size, and deformability, are a long-standing and challenging pursuit. Herein, a new strategy based on amphiphilic alternating copolymers for producing SCNPs with ultrasmall size and uniform structure is presented. SCNPs are obtained by folding the designed alternating copolymer in N,N-dimethylformamide (DMF) and fixing it through a photocatalyzed cycloaddition reaction of anthracene units. Molecular dynamics simulation confirms the solvophilic outer corona and solvophobic inner core structure of SCNPs. Furthermore, by adjusting the length of PEG units, precise control over the mean size of SCNPs is achieved within the range of 2.8 to 3.9 nm. These findings highlight a new synthetic strategy that enables enhanced control over morphology and internal structure while achieving ultrasmall and uniform size for SCNPs.

Exploring the interplay of liquid crystal orientation and spherical elastic shell deformation in spatial confinement.

The application of liquid crystal technology typically relies on the precise control of molecular orientation at a surface or interface. This control can be achieved through a combination of morphological and chemical methods. Consequently, variations in constrained boundary flexibility can result in a diverse range of phase behaviors. In this study, we delve into the self-assembly of liquid crystals within elastic spatial confinement by using the Gay-Berne model with the aid of molecular dynamics simulations. Our findings reveal that a spherical elastic shell promotes a more regular and orderly alignment of liquid crystals compared to a hard shell. Moreover, during the cooling process, the hard-shell confined system undergoes an isotropic-smectic phase transition. In contrast, the phase behavior within the spherical elastic shell closely mirrors the isotropic-nematic-smectic phase transition observed in bulk systems. This indicates that the orientational arrangement of liquid crystals and the deformations induced by a flexible interface engage in a competitive interplay during the self-assembly process. Importantly, we found that phase behavior could be manipulated by altering the flexibility of the confined boundaries. This insight offers a fresh perspective for the design of innovative materials, particularly in the realm of liquid crystal/polymer composites.

Small Things Make a Big Difference: Conductive Cross-Linking Sodium Alginate@MXene Binder Enables High-Volumetric-Capacity and High-Mass-Loading Li-S Battery.

Binders are crucial for maintaining the integrity of an electrode, and there is a growing need for integrating multiple desirable properties into the binder for high-energy batteries, yet significant challenges remain. Here, we successfully synthesized a new binder by cross-linking sodium alginate (SA) with MXene materials (Ti3C2Tx). Besides the improved adhesion and mechanical properties, the integrated SA@Ti3C2Tx binder demonstrates much improved electronic conductivity, which enables ruling out the fluffy conductive additive from the electrode component with enhanced volumetric capacity. When SA@Ti3C2Tx is used to fabricate sulfur (S) cathodes, the conductive-additive-free electrode demonstrates extremely high capacity (1422 mAh cm-3/24.5 mAh cm-2) under an S loading of 17.2 mg cm-2 for Li-S batteries. Impressively, the SA@Ti3C2Tx binder shows remarkable feasibility in other battery systems such as Na-S and LiFePO4 batteries. The proposed strategy of constructing a cross-linking conductive binder opens new possibilities for designing high-mass-loading electrodes with high volumetric capacity.

Ionic-Nanophase Hybridization of Nafion by Supramolecular Patching for Enhanced Proton Selectivity in Redox Flow Batteries.

Nafion, as the mostly used proton exchange membrane material in vanadium redox flow batteries (VRFBs), encounters serious vanadium permeation problems due to the large size difference between its anionic nanophase (3-5 nm) and cationic vanadium ions (∼0.6 nm). Bulk hybridization usually suppresses the vanadium permeation at the expense of proton conductivity since conventional additives tend to randomly agglomerate and damage the nanophase continuity from unsuitable sizes and intrinsic incompatibility. Here, we report the ionic-nanophase hybridization strategy of Nafion membranes by using fluorinated block copolymers (FBCs) and polyoxometalates (POMs) as supramolecular patching additives. The cooperative noncovalent interactions among Nafion, interfacial-active FBCs, and POMs can construct a 1 nm-shrunk ionic nanophase with abundant proton transport sites, preserved continuity, and efficient vanadium screeners, which leads to a comprehensive enhancement in proton conductivity, selectivity, and VRFB performance. These results demonstrate the intriguing potential of the supramolecular patching strategy in precisely tuning nanostructured electrolyte membranes for improved performance.

Covalent Organic Framework with Multiple Redox Active Sites for High-Performance Aqueous Calcium Ion Batteries.

Organic materials are promising for cation storage in calcium ion batteries (CIBs). However, the high solubility of organic materials in an electrolyte and low electronic conductivity remain the key challenges for high-performance CIBs. Herein, a nitrogen-rich covalent organic framework with multiple carbonyls (TB-COF) is designed as an aqueous anode to address those obstacles. TB-COF demonstrates a high reversible capacity of 253 mAh g-1 at 1.0 A g-1 and long cycle life (0.01% capacity decay per cycle at 5 A g-1 after 3000 cycles). The redox mechanism of Ca2+/H+ co-intercalated in COF and chelating with C═O and C═N active sites is validated. In addition, a novel C═C active site was identified for Ca2+ ion storage. Both computational and empirical results reveal that per TB-COF repetitive unit, up to nine Ca2+ ions are stored after three staggered intercalation steps, involving three distinct Ca2+ ion storage sites. Finally, the evolution process of radical intermediates further elucidates the C═C reaction mechanism.

The design of highly conductive and stretchable polymer conductors with low-load nanoparticles.

Highly conductive and stretchable polymer conductors fabricated from conductive fillers and stretchable polymers are urgently needed in flexible electronics, implants, soft robotics, etc. However, polymer conductors encounter the conductivity-stretchability dilemma, in which high-load fillers needed for high conductivity always result in the stiffness of materials. Herein, we propose a new design of highly conductive and stretchable polymer conductors with low-load nanoparticles (NPs). The design is achieved by the self-assembly of surface-modified NPs to efficiently form robust conductive pathways. We employ computer simulations to elucidate the self-assembly of the NPs in the polymer matrices under equilibrium and tensile states. The conductive pathways retain 100% percolation probability even though the loading of the NPs is lowered to ∼2% volume. When the tensile strain reaches 400%, the percolation probability of the ∼2% NP system is still greater than 25%. The theoretical prediction suggests a way for advancing flexible conductive materials.

Defect transition of smectic liquid crystals confined in spherical cavities.

The formation and transformation of defects in confined liquid crystals are fascinating fundamental problems in soft matter. Here, we use molecular dynamics (MD) simulations to study ellipsoidal liquid crystals (LCs) confined in a spherical cavity, which significantly affects the orientation and translation of LC molecules near the surface. The liquid-crystal droplet can present the isotropic to smectic-B phase transition through the smectic-A phase, as the number density of the LC molecules increases. We further find the change of LC structure from bipolar to watermelon-striped during the phase transition from smectic-A (SmA) to smectic-B (SmB) phases. Our results reveal the transition from bipolar defects to the inhomogeneous structures with the coexistence of nematic and smectic phases in smectic liquid-crystal droplets. We also study the influence of the sphere size in the range of 10σ0 ≤ Rsphere ≤ 50σ0 on the structural inhomogeneities. It shows a weak dependence on the sphere size. We further focus on how the structures can be affected by the interaction strength εGB-LJ. Interestingly, we find the watermelon-striped structure can be changed into a configuration with four defects at the vertices of a tetrahedron upon increasing the interaction strength. The liquid crystals at a strong interaction strength of εGB-LJ = 10.0ε0 show the two-dimensional nematic phase at the surface. We further present an explanation for the origin of the striped-pattern formation. Our results highlight the potential for using confinement to control these defects and their associated nanostructural heterogeneity.

Helical structures of achiral liquid crystals under cylindrical confinement.

Confined liquid crystals (LCs) exhibit complex and intriguing structures, which are fascinating fundamental problems in soft matter. The helical structure of cylindrical cavities is of great importance in LC studies, particularly for their application in optical devices. In this study, we employ molecular dynamics simulations to explore the behavior of achiral smectic-B LCs confined in narrow cylindrical cavities, where geometric frustration plays an important role. By increasing the cylinder size, LCs exhibit a transition from multi-helical to layered structures. Notably, we observe two stable structures, namely the helical structure and the layered structure, at moderate cylinder size. We also investigate the effects of the arrangement of cylindrical wall particles (hexagonal or square array) and anchoring strength on the LC structure. Our findings reveal that both the hexagonal array and strong anchoring strength promote the formation of helical structures. Our study provides novel insights into the confinement physics of LCs and highlights the potential for achieving helical structures in achiral LCs, which will expand the future applications of LCs.

External field induced defect transformation in circular confined Gay-Berne liquid crystals.

Normally, defects in two-dimensional, circular, confined liquid crystals can be classified into four types based on the position of singularities formed by liquid crystal molecules, i.e., the singularities located inside the circle, at the boundary, outside the circle, and outside the circle at infinity. However, it is considered difficult for small aspect ratio liquid crystals to generate all these four types of defects. In this study, we use molecular dynamics simulation to investigate the defect formed in Gay-Berne, ellipsoidal liquid crystals, with small aspect ratios confined in a circular cavity. As expected, we only find two types of defects (inside the circle and at the boundary) in circular, confined, Gay-Berne ellipsoids under static conditions at various densities, aspect ratios, and interactions between the wall and liquid crystals. However, when introducing an external field to the system, four types of defects can be observed. With increasing the strength of the external field, the singularities in the circular, confined system change from the inside to the boundary and the outside, and the farthest position that the singularities can reach depends on the strength of the external field. We further introduce an alternating, triangular wave, external field to the system to check if we can observe the transformation of different defects within an oscillating period. We find that the position of the singularities greatly depends on the oscillating intensity and oscillating period. By changing the oscillating intensity and oscillating period of the external field, the defect types can be adjusted, and the transformation between different defects can be easily observed. This provides a feasible way to modulate liquid crystal defects and investigate the transformation between different defects.

Non-equilibrium Nanoassemblies Constructed by Confined Coordination on a Polymer Chain.

Biological systems employ non-equilibrium self-assembly to create ordered nanoarchitectures with sophisticated functions. However, it is challenging to construct artificial non-equilibrium nanoassemblies due to lack of control over assembly dynamics and kinetics. Herein, we design a series of linear polymers with different side groups for further coordination-driven self-assembly based on shape-complementarity. Such a design introduces a main-chain confinement which effectively slows down the assembly process of side groups, thus allowing us to monitor the real-time evolution of lychee-like nanostructures. The function related to the non-equilibrium nature is further explored by performing photothermal conversion study. The ability to observe and capture non-equilibrium states in this supramolecular system will enhance our understanding of the thermodynamic and kinetic features as well as functions of living systems.

Interactions on Proteins Arising from the Self-Assembly of a Polyelectrolyte Brush.

Surfaces grafted with polyelectrolyte chains for excellent performance in protein antifouling are highly desired in many applications, such as biomedical implants and devices. In general, the adsorbing/resisting behaviors of proteins can be mainly attributed to the electrostatic interactions that are associated with the charge properties of proteins and polyelectrolytes. By coarse-grained molecular dynamics simulations, we examined the self-assembled structures of polyanion and polyzwitterion brushes as well as the interactions on negatively and positively charged proteins. We found that in addition to charges, the structural polarization induced by self-assembly with a certain charge distribution shows significant influences on protein behavior. The large-scale dipole-dipole interactions between brushes and proteins can dominate the behavior of proteins on the brushes under certain circumstances. To ensure simulation accuracy, we compared two models and found a polar Martini model that explicitly treats electrostatic interactions as long-ranged ones, giving a more reasonable structural description compared with the normal Martini model that truncates electrostatic interactions.

Colloidal cubic diamond photonic crystals through cooperative self-assembly.

Colloidal cubic diamond crystals with low-coordinated and staggered structures could display a wide photonic bandgap at low refractive index contrasts, which makes them extremely valuable for photonic applications. However, self-assembly of cubic diamond crystals using simple colloidal building blocks is still considerably challenging, due to their low packing fraction and mechanical instability. Here we propose a new strategy for constructing colloidal cubic diamond crystals through cooperative self-assembly of surface-anisotropic triblock Janus colloids and isotropic colloidal spheres into superlattices. In self-assembly, cooperativity is achieved by tuning the interaction and particle size ratio of colloidal building blocks. The pyrochlore lattice formed by self-assembly of triblock Janus colloids acts as a soft template to direct the packing of colloidal spheres into cubic diamond lattices. Numerical simulations show that this cooperative self-assembly strategy works well in a large range of particle size ratio of these two species. Moreover, photonic band structure calculations reveal that the resulting cubic diamond lattices exhibit wide and complete photonic bandgaps and the width and frequency of the bandgaps can also be easily adjusted by tuning the particle size ratio. Our work will open up a promising avenue toward photonic bandgap materials by cooperative self-assembly employing surface-anisotropic Janus or patchy colloids as a soft template.

A chiral smectic phase induced by an alternating external field.

Using simple achiral building blocks modulated by an external field to achieve chiral liquid crystal phases remains a challenge. In this study, a chiral helix liquid crystal phase is obtained for a simple Gay-Berne ellipsoid model under an alternating external field by using molecular dynamics simulations. Our results show that the chiral helix liquid crystal phase can be observed in a wide range of external field strengths when the oscillation period is smaller than the rotational characteristic diffusion timescale of ellipsoids. In addition, we find that the pitch and tilt angle of the helix structure can also be adjusted by changing the strength and oscillation period of the applied alternating external field. This may provide a feasible route for the regulation of chiral liquid crystal phases by an alternating external field.

The ENUF method-Ewald summation based on nonuniform fast Fourier transform: Implementation, parallelization, and application.

Computer simulations of model systems are widely used to explore striking phenomena in promising applications spanning from physics, chemistry, biology, to materials science and engineering. The long range electrostatic interactions between charged particles constitute a prominent factor in determining structures and states of model systems. How to efficiently calculate electrostatic interactions in simulation systems subjected to partial or full periodic boundary conditions has been a grand challenging task. In the past decades, a large variety of computational schemes has been proposed, among which the Ewald summation method is the most reliable route to accurately deal with electrostatic interactions between charged particles in simulation systems. In addition, extensive efforts have been done to improve computational efficiencies of the Ewald summation based methods. Representative examples are approaches based on cutoffs, reaction fields, multi-poles, multi-grids, and particle-mesh schemes. We sketched an ENUF method, an abbreviation for the Ewald summation method based on the nonuniform fast Fourier transform technique, and have implemented this method in particle-based simulation packages to calculate electrostatic energies and forces at micro- and mesoscopic levels. Extensive computational studies of conformational properties of polyelectrolytes, dendrimer-membrane complexes, and ionic fluids demonstrated that the ENUF method and its derivatives conserve both energy and momentum to floating point accuracy, and exhibit a computational complexity of O N log N with optimal physical parameters. These ENUF based methods are attractive alternatives in molecular simulations where high accuracy and efficiency of simulation methods are needed to accelerate calculations of electrostatic interactions at extended spatiotemporal scales.

The coarse-grained models of poly(ethylene oxide) and poly(propylene oxide) homopolymers and poloxamers in big multipole water (BMW) and MARTINI frameworks.

Polyethylene oxide (PEO) and poly(propylene oxide) (PPO), especially their tri-block copolymers PEO-PPO-PEO (poloxamers), have a broad range of applications in biotechnology and medical science. Understanding their specific interactions with biomembranes is the key to unveil the unique features of poloxamers either as membrane-healing or membrane pore-forming agents. Based on the coarse-graining convention of the MARTINI force field and the big multipole water (BMW) model, which has a three charged site topology and can reproduce the correct dipole moment of four-water clusters, we generated coarse-grained (CG) models with analytical and numerical potentials for PEO and PPO homopolymers and poloxamers in dilute solution. The effective bonded interaction potentials between CG beads were determined from the probability distributions of bond lengths, angles and dihedrals that are determined from atomistic simulations. The nonbonded interaction parameters were fine-tuned to reproduce the conformational properties of atomistic PEO and PPO homopolymers and poloxamers via extensive CG simulations of PEO and PPO homopolymers and poloxamers in a BMW water environment. The reported CG models provide a promising framework for a comprehensive understanding of the microstructural, conformational, and dynamic properties of poloxamers and their delicate interactions with other species in an explicit water environment.

Effect of the Self-Assembled Structures of Hydrated Polyzwitterionic and Polyanionic Brushes on Their Self-Cleaning Capabilities.

The capability of polyelectrolyte brushes to spontaneously clean oil fouling via water is determined by factors including water wettability and the self-assembled structures of hydrated polyelectrolytes. Although the charged groups of polyelectrolytes provide the original source of water wettability, the self-assembled structures play a significant role in the self-cleaning performances. Here, we employ coarse-grained molecular dynamics simulations to study the general self-cleaning characteristics of two types of surface-grafted polyelectrolyte brushes (i.e., zwitterionic and anionic polyelectrolytes). It has been found that the high grafting density is favorable to fouling reduction for both polyzwitterions and polyanions. To be specific, the hydrated polyzwitterions form an intermolecular cross-linked network via zwitterionic complexes, resulting in better self-cleaning capabilities than the polyanions at lower grafting densities. However, polyanions form bundles with each consisting of several chains via hydrophobic interactions and electrostatic repulsions presenting better self-cleaning performances than the polyzwitterions at higher grafting densities.

Coupling and decoupling between translational and rotational dynamics in supercooled monodisperse soft Janus particles.

We perform dynamics simulations to investigate the translational and rotational glassy dynamics in a glass-forming liquid of monodisperse soft Janus particles. We find that, with decreasing temperature, the mean-square angular displacement shows no clear plateau in the caging region, in contrast with the apparent caging behavior of translational motion. By defining a reorientational mean-square angular displacement, the caging behavior of rotational motion can be recognized. On approaching the glass transition (decreasing temperature), the coupling between translational and rotational relaxation increases, while the coupling between translational and rotational diffusion decreases, whereas the coupling between translational and reorientational diffusion increases. The strong decoupling between translational and rotational diffusion is due to the suppressed translational mobility but promoted rotational mobility of soft Janus particles. We think that the low-T SE and SED decoupling is mainly attributed to hopping motion of soft Janus particles, whereas the high-T SE and SED decoupling is mainly attributed to collective cage motion of soft Janus particles. Our results demonstrate that interaction anisotropy has a critical effect on the translational and rotational dynamics of soft Janus particles.

Heterogeneous dynamics of unentangled chains in polymer nanocomposites.

We present a systematic investigation on the effect of adding nanoparticles on the dynamics of polymer chains by using coarse-grained molecular dynamics simulation. The dynamics is characterized by three aspects: molecular motion, relaxation at different length scales, and dynamical heterogeneity. It is found that the motion of polymer chains slows down and the deviation from Gaussian distribution becomes more pronounced with increasing nanoparticle volume fractions. For polymer nanocomposites with R ≤ Rg, the relaxation at the wave vector q = 7.0 displays multistep decay, consistent with the previous reports in strongly interacting polymer nanocomposites. Moreover, a qualitatively universal law is established that dynamic heterogeneity at whole chain's scale follows a nonmonotonic increase with increasing nanoparticle loadings, where the volume fraction of the maximum dynamic heterogeneity corresponds to the particle loading when the average distance between nanoparticles is equal to the Kuhn length of polymer chains. We show that the decoupling between whole chain's dynamics and segment dynamics is responsible for the nonmonotonic behavior of dynamic heterogeneity of whole chains.

Electrostatic interactions in soft particle systems: mesoscale simulations of ionic liquids.

Computer simulations provide a unique insight into the microscopic details, molecular interactions and dynamic behavior responsible for many distinct physicochemical properties of ionic liquids. Due to the sluggish and heterogeneous dynamics and the long-ranged nanostructured nature of ionic liquids, coarse-grained meso-scale simulations provide an indispensable complement to detailed first-principles calculations and atomistic simulations allowing studies over extended length and time scales with a modest computational cost. Here, we present extensive coarse-grained simulations on a series of ionic liquids of the 1-alkyl-3-methylimidazolium (alkyl = butyl, heptyl-, and decyl-) family with Cl, [BF4], and [PF6] counterions. Liquid densities, microstructures, translational diffusion coefficients, and re-orientational motion of these model ionic liquid systems have been systematically studied over a wide temperature range. The addition of neutral beads in cationic models leads to a transition of liquid morphologies from dispersed apolar beads in a polar framework to that characterized by bi-continuous sponge-like interpenetrating networks in liquid matrices. Translational diffusion coefficients of both cations and anions decrease upon lengthening of the neutral chains in the cationic models and by enlarging molecular sizes of the anionic groups. Similar features are observed in re-orientational motion and time scales of different cationic models within the studied temperature range. The comparison of the liquid properties of the ionic systems with their neutral counterparts indicates that the distinctive microstructures and dynamical quantities of the model ionic liquid systems are intrinsically related to Coulombic interactions. Finally, we compared the computational efficiencies of three linearly scaling O(N log N) Ewald summation methods, the particle-particle particle-mesh method, the particle-mesh Ewald summation method, and the Ewald summation method based on a non-uniform fast Fourier transform technique, to calculate electrostatic interactions. Coarse-grained simulations were performed using the GALAMOST and the GROMACS packages and hardware efficiently utilizing graphics processing units on a set of extended [1-decyl-3-methylimidazolium][BF4] ionic liquid systems of up to 131 072 ion pairs.

General patchy ellipsoidal particle model for the aggregation behaviors of shape- and/or surface-anisotropic building blocks.

We present a general patchy ellipsoidal particle model suitable for conducting dynamics simulations of the aggregation behaviors of various shape- and/or surface-anisotropic colloids, especially patchy ellipsoids with continuously variable shape and tunable patchiness. To achieve higher computational efficiency in dynamics simulations, we employ a multi-GPU acceleration technique based on a domain decomposition algorithm. The validation and performance evaluation of this GPU-assisted model are performed by simulating several typical benchmark systems of non-patchy and patchy ellipsoids. Given the generality and efficiency of our GPU-assisted patchy ellipsoidal particle model, it will provide a highly feasible dynamics simulation framework to investigate the aggregation behaviors of anisotropic soft matter systems comprised of shape- and/or surface-anisotropic building blocks.