Entropy-favorable adsorption of polymer-grafted nanoparticles at fluid-fluid interfaces.

The adsorption of polymer-grafted nanoparticles at interfaces is a problem of fundamental interest in physics and soft materials. This adsorption behavior is governed by the interplay between interaction potentials and entropic effects. Here, we use molecular dynamics simulations and umbrella sampling methods to study the adsorption behavior of a Janus-like homopolymer-grafted nanoparticle at fluid-fluid interfaces. By calculating the potential of the mean force as the particle moves from fluid A to the interface, the adsorption energy Ea can be obtained. When two homopolymer chains with types A and B are grafted to the opposite poles of the particle, Ea shows a scaling behavior with respect to chain length N: Ea ∝ N0.598. This is determined by the interactions between polymers and fluids. The enthalpy dominates, and the entropy effects mainly come from the rotational entropy loss of the polymer-grafted nanoparticle at interfaces, which disfavors the stabilization of particles at interfaces. When the grafted polymer number m is large, the adsorption energy exhibits a linear dependence on m. While the enthalpy dominates the behavior, the entropy becomes significant at a larger chain length of N = 15, where the configurational entropy of the polymer chains dominates the entropy of the system. The globule-coil transition occurs when polymers move from poor solvents to good solvents, increasing the configurational entropy and favoring the stabilization of particles at interfaces. Our study provides novel insights into the stabilization mechanism of polymer-grafted nanoparticles at interfaces and reveals the stabilization mechanism favored by the configurational entropy of grafted polymer chains.

Relative kinetic stability of defect patterns in two-dimensional nematic liquid crystals with rectangular confinement.

Guiding and dynamically modulating topological defects are critical challenges in defect engineering of liquid crystals. Here, we employ molecular dynamics simulations to investigate the transition dynamics and relative kinetic stability of defect patterns in two-dimensional nematic Gay-Berne liquid crystals confined within rectangular geometries. We observe the formation of various defect patterns including long-axis, diagonal, X-shaped, composite, and bend configurations under different confinement conditions. The competition between boundary effects and the uniformity of nematic orientation induces the continuous realignment of liquid crystal molecules, facilitating the spatially continuous transformation of defect patterns over time. This transition involves changes in both defect types and their locations, typically initiating from defect regions. Furthermore, we demonstrate that the relative stability of these defect patterns can be effectively controlled by adjusting confinement parameters and external field conditions. Our findings provide fundamental insights into the transition kinetics of defect patterns in confined nematic liquid crystals, thereby enhancing our ability to manipulate topological defects for advanced applications.

Deposition patterns formed by the evaporation of linear diblock copolymer solution nanodroplets on solid surfaces.

The evaporation-induced deposition pattern of the linear diblock copolymer solution has attracted attention in recent years. Given its critical applications, we study deposition patterns of the linear diblock copolymer solution nanodroplet on a solid surface (the wall) by molecular dynamics simulations. This study focuses on the influence of the nonbonded interaction strength, including the interaction between the wall and polymer blocks (ɛAW and ɛBW), the interaction between the solvent and the wall (ɛSW), and the interaction between polymer blocks (ɛAB). Conditions leading to diverse deposition patterns are explored, including the coffee-ring and the volcano-like structures. The formation of the coffee-ring structure is attributed to receding interfaces, the heterogeneity inside the droplet, and the self-assembly of polymer chains. This study contributes to the establishment of guidelines for designing deposition patterns of the linear diblock copolymer solution nanodroplet, which facilitates practical applications such as inkjet printing.

Structural heterogeneity in tetra-armed gels revealed by computer simulation: Evidence from a graph theory assisted characterization.

Designing homogeneous networks is considered one typical strategy for solving the problem of strength and toughness conflict of polymer network materials. Experimentalists have proposed the hypothesis of obtaining a structurally homogeneous hydrogel by crosslinking tetra-armed polymers, whose homogeneity was claimed to be verified by scattering characterization and other methods. Nevertheless, it is highly desirable to further evaluate this issue from other perspectives. In this study, a coarse-grained molecular dynamics simulation coupled with a stochastic reaction model is applied to reveal the topological structure of a polymer network synthesized by tetra-armed monomers as precursors. Two different scenarios, distinguished by whether internal cross-linking is allowed, are considered. We introduce the Dijkstra algorithm from graph theory to precisely characterize the network structure. The microscopic features of the network structure, e.g., loop size, dispersity, and size distribution, are obtained via the Dijkstra algorithm. By comparing the two reaction scenarios, Scenario II exhibits an overall more idealized structure. Our results demonstrate the feasibility of the Dijkstra algorithm for precisely characterizing the polymer network structure. We expect this work will provide a new insight for the evaluation and description of gel networks and further help to reveal the dynamic process of network formation.

PolyNC: a natural and chemical language model for the prediction of unified polymer properties.

Language models exhibit a profound aptitude for addressing multimodal and multidomain challenges, a competency that eludes the majority of off-the-shelf machine learning models. Consequently, language models hold great potential for comprehending the intricate interplay between material compositions and diverse properties, thereby accelerating material design, particularly in the realm of polymers. While past limitations in polymer data hindered the use of data-intensive language models, the growing availability of standardized polymer data and effective data augmentation techniques now opens doors to previously uncharted territories. Here, we present a revolutionary model to enable rapid and precise prediction of Polymer properties via the power of Natural language and Chemical language (PolyNC). To showcase the efficacy of PolyNC, we have meticulously curated a labeled prompt-structure-property corpus encompassing 22 970 polymer data points on a series of essential polymer properties. Through the use of natural language prompts, PolyNC gains a comprehensive understanding of polymer properties, while employing chemical language (SMILES) to describe polymer structures. In a unified text-to-text manner, PolyNC consistently demonstrates exceptional performance on both regression tasks (such as property prediction) and the classification task (polymer classification). Simultaneous and interactive multitask learning enables PolyNC to holistically grasp the structure-property relationships of polymers. Through a combination of experiments and characterizations, the generalization ability of PolyNC has been demonstrated, with attention analysis further indicating that PolyNC effectively learns structural information about polymers from multimodal inputs. This work provides compelling evidence of the potential for deploying end-to-end language models in polymer research, representing a significant advancement in the AI community's dedicated pursuit of advancing polymer science.

Exploring mechanical properties and failure mechanisms of aramid and PBO crystals through molecular dynamics simulations.

Molecular dynamics simulations were used to analyze the mechanical properties and failure processes of poly(p-phenylene-terephthalamide) (PPTA), poly(p-phenylene-benzimidazole-terephthalamide) (PBIA), PBIA-PPTA (formed by 1:1 copolymerization of PPTA and PBIA), and poly(p-phenylene-benzobisoxazole) (PBO) crystals at different strain rates and temperatures. The failure stress and strain were found to be linear with the temperature and logarithmic strain rate. Moreover, based on the kinetic theory of fracture and the comprehensive simulation results, we formulated a model that describes the failure stress of the aforementioned crystals under varying strain rates and temperatures. Through the analysis of the failure process, we found that in PPTA, PBIA, and PBIA-PPTA crystals, the bond failure probability is correlated with the strain rate and temperature. The examination of bond lengths and angles unveiled that bonds with larger initial aligning angles are more susceptible to failure during the strain process. Intriguingly, the stretching process induced a conformational change in the PBO molecular chain, leading to a deviation from the linear relation in its stress-strain curve.

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.

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.

Automatic Multiscale Method of Building up a Cross-linked Polymer Reaction System: Bridging SMILES to the Multiscale Molecular Dynamics Simulation.

An automatic method is introduced to generate the initial configuration and input file from SMILES for multiscale molecular dynamics (MD) simulation of cross-linked polymer reaction systems. Inputs are a modified version of SMILES of all the components and conditions of coarse-grained (CG) and all-atom (AA) simulations. The overall process comprises the following steps: (1) Modified SMILES inputs of all the components are converted to 3-dimensional coordinates of molecular structures. (2) Molecular structures are mapped to the coarse-grained scale, followed by a CG reaction simulation. (3) CG beads are backmapped to the atomic scale after the CG reaction. (4) An AA productive run is finally performed to analyze volume shrinkage, glass transition, and atomic detail of network structure. The method is applied to two common epoxy resin reactions, that is, the cross-linking process of DGEVA (diglycidyl ether of vanillyl alcohol) and DHAVA (dihydroxyaminopropane of vanillyl alcohol) and that of DGEBA (diglycidyl ether of bisphenol A) and DETA (diethylenetriamine). These components form network structures after the CG cross-linking reaction and are then backmapped to calculate properties in the atomic scale. The result demonstrates that the method can accurately predict volume shrinkage, glass transition, and all-atom structure of cross-linked polymers. The method bridges from SMILES to MD simulation trajectories in an automatic way, which shortens the time of building up cross-linked polymer reaction model and suitable for high-throughput computations.

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.

The advantages of nanoparticle surfactants over Janus nanoparticles on structuring liquids.

The nanoparticle (NP) surfactants generated in situ by binding NPs and polymers can assemble into an elastic NP monolayer at the interface of two immiscible liquids, structuring the liquids. Janus NPs can be more strongly bound to the interface than the NP surfactants, but they are unable to structure liquids into complex shapes due to the difficulty of assembling the jamming arrays. By molecular dynamics simulations, we give an insight into the better performance of NP surfactants than Janus NPs on dynamically structuring liquids. The high energy binding of Janus NPs to the interface will drive the Janus NPs to assemble into micelles in binary liquids. The micelles are stabilized in one liquid by encapsulating a little of the other liquid, hindering interfacial adsorption when the interface is marginally extended upon liquid deformation. In contrast, the in situ formed NP surfactants can rapidly fill the enlarged interfacial area to arrest the consecutive shape changes of the liquids. Moreover, NP surfactants can be designed with an appropriate coverage ratio (≤50%) of NP surface bearing host-guest sites to avoid dissolution and impart a desirable mechanical elasticity to their assembly.

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.

Multiple 2D crystal structures in bilayered lamellae from the direct self-assembly of 3D systems of soft Janus particles.

Numerous crystals and Frank-Kasper phases in two-dimensional (2D) systems of soft particles have been presented by theoretical investigations. How to realize 2D crystals or Frank-Kasper phases via the direct self-assembly of three-dimensional (3D) systems remains an important issue. Here, through numerical simulations, we report the surprising finding of multiple 2D crystal structures in bilayered lamellae from the direct self-assembly of 3D systems of soft Janus particles. With varying the patch size and particle density, soft Janus particles, which exhibit very similar self-assembly behavior to giant amphiphiles, spontaneously form ordered bilayered lamellae. Within each layer of the bilayered lamellae, we find abundant highly-ordered 2D crystals including the Frank-Kasper σ phase and open kagome lattice. The kinetic mechanisms of the formation of these 2D crystals within the layers are revealed, and include a classical one-step nucleation mechanism and a two-step nucleation mechanism. Our findings suggest a simple route towards 2D crystals via the direct self-assembly of 3D systems of amphiphilic Janus building blocks.

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.

Transition kinetics of defect patterns in confined two-dimensional smectic liquid crystals.

Topological defects in liquid crystals under confined geometries have attracted extensive research interests. Here, we perform molecular dynamics simulations to investigate the formation and transition of defect patterns in two-dimensional smectic Gay-Berne liquid crystals with a simple rectangular confinement boundary. Two typical types of defect patterns, bridge and diagonal defect patterns, are observed, which can be transformable continuously between each other over time. The transition usually starts from the line or point defect regions, and the competition between neighboring and opposite boundary effects induces the continuous realignments of the smectic layers to connect the neighboring or opposite walls. The relative stability of these two defect patterns can be controlled by changing the confinement conditions. These results deepen our understanding of transition kinetics of defect patterns in confined liquid crystals.

Intercluster Exchange-Stabilized Novel Complex Colloidal χc Phase.

Designing complex cluster crystals with a specific function using simple colloidal building blocks remains a challenge in materials science. Herein, we propose a conceptually new design strategy for constructing complex cluster crystals via hierarchical self-assembly of simple soft Janus colloids. A novel and previously unreported colloidal cluster-χ (χc) phase, which resembles the essential structural features of α-manganese but at a larger length scale, is obtained through molecular dynamics simulations. The formation of the χc phase undergoes a remarkable two-step self-assembly process, that is, the self-assembly of clusters with specific size dispersity from Janus colloids, followed by the highly ordered organization of these clusters. More importantly, the dynamic exchange of particles between these clusters plays a critical role in stabilizing the χc phase. Such a conceptual design framework based on intercluster exchange has the potential to effectively construct novel complex cluster crystals by hierarchical self-assembly of colloidal building blocks.

Solvent-Evaporation Induced and Mechanistic Entropy-Enthalpy-Balance Controlled Polymer Patch Formation on Nanoparticle Surfaces.

The formation of polymer-patch nanoparticles (PNPs) involves a condensation process of grafted chains on a nanoparticle (NP) surface, which is conventionally achieved via a fine-tuning of the solvent quality. However, such a critical solvent condition differs dramatically between polymers, and the formation mechanism of different patchy structures remains under debate. In this study, we demonstrate by a combined simulation and experimental study that such a surface-patterning process can be easily achieved via a simple solvent evaporation process, which creates a natural nonsolvent condition and is, in principle, adaptable for all polymers. More importantly, we find that patchy structures are controlled by a delicate balance between enthalpic interaction and the entropy penalty of grafted chains. A small variation of cohesive energy density can lead to a dramatic change in patch structure. This work offers a robust yet easy approach for the fabrication of PNPs and provides new insights into polymer segregation on spherical surfaces.

Softness-Enhanced Self-Assembly of Pyrochlore- and Perovskite-like Colloidal Photonic Crystals from Triblock Janus Particles.

It remains extremely challenging to build three-dimensional photonic crystals with complete photonic bandgaps by simple and experimentally realizable colloidal building blocks. Here, we demonstrate that particle softness can enhance both the self-assembly of pyrochlore- and perovskite-like lattice structures from simple deformable triblock Janus colloids and their photonic bandgap performances. Dynamics simulation results show that the region of stability of pyrochlore lattices can be greatly expanded by appropriately increasing softness, and the perovskite lattices are unexpectedly obtained at enough high softness. Photonic calculations show that the direct pyrochlore lattices formed from overlapping soft triblock Janus particles exhibit even larger photonic bandgaps than the ideal nonoverlapping pyrochlore lattice, and proper overlap arising from softness can also dramatically improve the photonic properties of the inverse pyrochlore and perovskite lattices. Our study offers a new and feasible self-assembly path toward three-dimensional photonic crystals with large and robust photonic bandgaps.

Mechanism of periodic field driven self-assembly process.

Dissipative self-assembly, a ubiquitous type of self-assembly in biological systems, has attracted a lot of attention in recent years. Inspired by nature, dissipative self-assembly driven by periodic external fields is often adopted to obtain controlled out-of-equilibrium steady structures and materials in experiments. Although the phenomena in dissipative self-assembly have been discovered in the past few decades, fundamental methods to describe dynamical self-assembly processes and responsiveness are still lacking. Here, we develop a theoretical framework based on the equations of motion and Floquet theory to reveal the dynamic behavior changing with frequency in the periodic external field driven self-assembly. Using the dissipative particle dynamics simulation method, we then construct a block copolymer model that can self-assemble in dilute solution to confirm the conclusions from the theory. Our theoretical framework facilitates the understanding of dynamic behavior in a periodically driven process and provides the theoretical guidance for designing the dissipative conditions.