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.

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.

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.

How Molecular Chiralities of Bis(mandelato)borate Anions Affect Their Binding Structures With Alkali Metal Ions and Microstructural Properties in Tetraalkylphosphonium Ionic Liquids.

Spiroborate anion-based inorganic electrolytes and ionic liquids (ILs) have fascinating electrochemical and tribological properties and have received widespread attention in industrial applications. The molecular chiralities of spiroborate anions have a significant effect on the microstructures and macroscopic functionalities of these ionic materials in application and thus deserve fundamental consideration. In the current work, we performed quantum chemistry calculations to address the binding strength and coordination structures of chiral bis(mandelato)borate ([BMB]) anions with representative alkali metal ions, as well as the electronic properties of alkali metal ion-[BMB] ion pair complexes. The optimized [BMB] conformers are categorized into V-shaped, bent, and twisted structures with varied electrostatic potential contours and conformational energies and distinct alkali metal ion-[BMB] binding structures. Alkali metal ions have additional associations with phenyl groups in V-shaped [BMB] conformers owing to preferential cation-π interactions. Furthermore, the effects of the molecular chiralities of [BMB] anions on the thermodynamics and microstructural properties of tetraalkylphosphonium [BMB] ILs were studied by performing extensive atomistic interactions. Oxygen atoms in [BMB] anions have competitive hydrogen bonding interactions with hydrogen atoms in cations depending on the molecular chiralities and steric hindrance effects of [BMB] anions. However, the molecular chiralities of [BMB] anions have a negligible effect on the liquid densities of tetraalkylphosphonium [BMB] ILs and the spatial distributions of boron atoms in anions around phosphorous atoms in cations. Enlarging tetraalkylphosphonium cation sizes leads to enhanced cation-anion intermolecular hydrogen bonding and Coulombic interactions due to enhanced segregation of polar groups in apolar networks in heterogeneous IL matrices, as verified by scattering structural functions.

Feature vector clustering molecular pairs in computer simulations.

A clustering framework is introduced to analyze the microscopic structural organization of molecular pairs in liquids and solutions. A molecular pair is represented by a representative vector (RV). To obtain RV, intermolecular atom distances in the pair are extracted from simulation trajectory as components of the key feature vector (KFV). A specific scheme is then suggested to transform KFV to RV by removing the influence of permutational molecular symmetry on the KFV as the predicted clusters should be independent of possible permutations of identical atoms in the pair. After RVs of pairs are obtained, a clustering analysis technique is finally used to classify all the RVs of molecular pairs into the clusters. The framework is applied to analyze trajectory from molecular dynamics simulations of an ionic liquid (trihexyltetradecylphosphonium bis(oxalato)borate ([P6,6,6,14 ][BOB])). The molecular pairs are successfully categorized into physically meaningful clusters, and their effectiveness is evaluated by computing the product moment correlation coefficient (PMCC). (Willett, Winterman, and Bawden, J. Chem. Inf. Comput. Sci. 1986, 26, 109-118; Downs, Willett, and Fisanick, J. Chem. Inf. Comput. Sci. 1994, 34, 1094-1102) It is observed that representative configurations of two clusters are related to two energy local minimum structures optimized by density functional theory (DFT) calculation, respectively. Several widely used clustering analysis techniques of both nonhierarchical (k-means) and hierarchical clustering algorithms are also evaluated and compared with each other. The proposed KFV technique efficiently reveals local molecular pair structures in the simulated complex liquid. It is a method, which is highly useful for liquids and solutions in particular with strong intermolecular interactions. © 2019 Wiley Periodicals, Inc.

Enhancement of surface nonwettability by grafting loops.

We present a computer simulation study on the nonwettability of a flat surface tethered with deformable looped polymer chains. Two kinds of loops are studied: monodispersed loops (loops with the same length) and polydispersed loops (loops with different lengths). Both kinds of loops include two arrangements: with regularly tethered sites and with randomly tethered sites. Regularly grafted loops form typical grooves on the surface, while randomly grafted loops form a more rugged surface. For monodispersed loops, we analyze the factors that influence the nonwettability when varying the rigidity of the loops. The loops are divided into two categories based on their rigidity according to our previous analysis procedure (Phys. Chem. Chem. Phys., 2016, 18, 18767-18775): rigid loops and flexible loops. It is found that the loop can partially form a re-entrant-like structure, which is helpful to increase the nonwettability of the surface. The surfaces with grafted loops have increased nonwettability, especially those grafted with flexible chains. However, the contact angle on the loop structure cannot further increase for the rigid chains due to a large top layer density (Phys. Chem. Chem. Phys., 2016, 18, 18767-18775). For polydispersed loops, the contact angle is highly related to the rigidity of the long loops that contact the droplet. Different from monodispersed loops, the mechanism of the nonwettability of polydispersed loops is attributed to the supporting ability (rigidity) of long loops.

Understanding the wettability of a hairy surface: effect of hair rigidity and topology.

We present a computer simulation study on the wettability of a hairy surface with different topological structures such as single hairs, hair bundles and network structure. Superficially, for end-tethered rigid hairs or flexible hairs, the nonwettability of the substrate should be analyzed in completely different ways. For rigid hairs, the contact angle is dominantly dependent on the top layer density of hairs. A larger top layer density leads to a larger interaction between droplet and surface, as well as a lower contact angle. For flexible hairs, the nonwettability is determined by the typical properties of hairs right below the droplet, e.g., the chemistry of the surface, the topography and strength of the hair bundle/network or nonwetted area below the projection of the droplet. Nevertheless, it is also possible to generalize these aspects into a uniform procedure, which implies an intrinsic consistent mechanism of the dewetting behavior for droplets on such hairy surfaces. Counterintuitively, we also suggest that the surface which can strongly resist the transition to the Wenzel state does not necessarily lead to a large contact angle, especially in a system where the droplet is treated as liquid bulk. This study helps to build up guidelines for the design of nonwetting surface materials.

Tuning surface wettability by designing hairy structures.

We present a molecular dynamics simulation study on the controlling factors that influence the wettability of a hairy surface. By adopting the hairs with appropriate grafting density, hair length, and hair rigidity, the hairy surface shows good performance on droplet repellency. When the droplet sits on the hairy surface, the flexible hairs can spontaneously bundle with the appropriate amount of neighboring hairs to enhance the surface hydrophobicity, thus providing a new possibility to control the surface wettability. The hairy surface with tunable grafting density and hair rigidity, bridges the gap between surfaces with soft polymer brushes and surfaces that are completely hard but porous.