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.

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.

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.

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.

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.

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.

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.

Mechanisms of Defect Correction by Reversible Chemistries in Covalent Organic Frameworks.

Reversible chemistries have been extensively explored to construct highly crystalline covalent organic frameworks (COFs) via defect correction. However, the mechanisms of defect correction that can explain the formation of products as single crystals, polycrystal/crystallites, or amorphous solids remain unknown. Herein, we employed molecular dynamics simulations combined with a polymerization model to investigate the growth kinetics of two-dimensional COFs. By virtue of the Arrhenius two-state model describing reversible reactions, we figured out the conditions in terms of active energy and binding energy for different products. Specifically, the ultraslow growth of COFs under high reversibility of reactions corresponding to low binding energies resulted in a single crystal by inhibiting the emergence of nuclei as well as correcting defects through continually dropping small defective fragments off at crystal boundaries. High bonding energies responsible for the high nucleation rate and rapid growth that incorporated defects in crystals and caused the division of crystals through defect correcting processes led to small crystallites or polycrystals. The insights into the mechanisms help us to understand and further control the growth kinetics by exploiting reversible conditions to synthesize COFs of higher quality.

A controlling parameter of topological defects in two-dimensional covalent organic frameworks.

Synthesis of covalent organic frameworks with long-range molecular ordering is an outstanding challenge due to the fact that defects against predesigned topological symmetries are prone to form and break crystallization. The physical origins and controlling parameters of topological defects remain scarcely understood. By virtue of molecular dynamics simulations, we found that pentagons for combination [C4 + C4] and [C4 + C2] and heptagons for [C3 + C3] and [C3 + C2] were initial defects for growth dynamics with both uncontrolled and suppressed nucleation, further inducing more complex defects. The defects can be significantly reduced by achieving the growth with monomers added to a single nucleus, agreeing well with previous simulations and experiments. To understand the nature of defects, we proposed a parameter φ to describe the range of biased rotational angle between two monomers, within which chemical reactions are allowed. The parameter φ shows a monotonic relationship with defect population, which is demonstrated to be highly computable by using density functional theory calculations. When φ < 20, we can even observe defect-free growth for the four combinations, irrespective of growth dynamics. The results are essential for screening and designing condensation reactions for the synthesis of single crystals of high quality.

Brownian Diffusion of Individual Janus Nanoparticles at Water/Oil Interfaces.

Janus nanoparticles could exhibit a higher interfacial activity and adsorb stronger to fluid interfaces than homogeneous nanoparticles of similar sizes. However, little is known about the interfacial diffusion of Janus nanoparticles and how it compares to that of homogeneous ones. Here, we employed fluorescence correlation spectroscopy to study the lateral diffusion of ligand-grafted Janus nanoparticles adsorbed at water/oil interfaces. We found that the diffusion was significantly slower than that of homogeneous nanoparticles. We carried out dissipative particle dynamic simulations to study the mechanism of interfacial slowdown. Good agreement between experimental and simulation results has been obtained only provided that the flexibility of ligands grafted on the nanoparticle surface was taken into account. The polymeric ligands were deformed and oriented at an interface so that the effective radius of Janus nanoparticles is larger than the nominal one obtained by measuring the diffusion in bulk solution. These findings highlight further the critical importance of the ligands grafted on Janus nanoparticles for applications involving nanoparticle adsorption at an interface, such as oil recovery or two-dimensional self-assembly.