发表论文

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.

Spider silk is one of the most robust natural materials, which has extremely high strength in combination with great toughness and good elasticity. Inspired by spider silk but beyond it, a healable and recyclable supramolecular elastomer, possessing superhigh true stress at break (1.21 GPa) and ultrahigh toughness (390.2 MJ m−3), which are, respectively, comparable to and ≈2.4 times higher than those of typical spider silk, is developed. The elastomer has the highest tensile strength (ultimate engineering stress, 75.6 MPa) ever recorded for polymeric elastomers, rendering it the strongest and toughest healable elastomer thus far. The hyper-robust elastomer exhibits superb crack tolerance with unprecedentedly high fracture energy (215.2 kJ m−2) that even exceeds that of metals and alloys, and superhigh elastic restorability allowing dimensional recovery from elongation over 12 times. These extraordinary mechanical performances mainly originate from the meticulously engineered hydrogen-bonding segments, consisting of multiple acylsemicarbazide and urethane moieties linked with flexible alicyclic hexatomic spacers. Such hydrogen-bonding segments, incorporated between extensible polymer chains, aggregate to form geometrically confined hydrogen-bond arrays resembling those in spider silk. The hydrogen-bond arrays act as firm but reversible crosslinks and sacrificial bonds for enormous energy dissipation, conferring exceptional mechanical robustness, healability, and recyclability on the elastomer.

In this study, we investigate the motion of polymer segments in polymer/nanoparticle composites by varying nanoparticle (NP) volume fractions. By studying the probability distribution of segment displacement, segment trajectory, and the square displacement of segment, we find the intermittent motion of segments, accompanied with the coexistence of slow and fast segments in polymer nanocomposites (PNCs). The displacement distribution of segments exhibits an exponential tail, rather than a Gaussian form. The intermittent dynamics of chain segments is comprised of a long-range jump motion and a short-range localized motion, which is mediated by the weakly attractive interaction between NP and chain segment and the strong confinement induced by NPs. Meanwhile, the intermittent motion of chain segments can be described by the adsorption-desorption transition at low particle loading and confinement effect at high particle loading. These findings may provide important information for understanding the anomalous motion of polymer chains in the presence of NPs.

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.

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.

We systematically investigated the roles of tail length on the self-assembly of shape amphiphiles composed of a hydrophobic polymer chain (tail) and a hydrophilic nanoparticle in selective solvent using Brownian dynamics simulations. The shape amphiphiles exhibited a variety of self-assembled aggregate morphologies which can be tuned by changing tail length (n) in combination with amphiphile concentration (φ) and system temperature (T*). Specifically, at high φ with T*=1.4, the morphology varied following the sequence “spheres → cylinders → vesicles” upon increasing n, agreeing well with experimental observations. At low φ with T*=1.4 or at high φ with T*=1.2, the morphology sequence becomes “spheres or spheres and cylinders mixture → cylinders → vesicles → spheres” upon increasing n, which has not been found experimentally. Two morphological phase diagrams depending on n and φ were constructed for T*=1.4 and 1.2, respectively. The rich phase behaviors on varying tail length could provide the feasible routes to fabricate target aggregate morphologies in various applications, especially for the vesicles with tunable thickness of membranes that are crucial in drug and gene delivery.

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 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.

Inspired by the distinct functions of desert beetles with efficient droplet nucleation and lotus leaves with excellent droplet removal, an integrated method is presented for the design of a superhydrophobic surface decorated with hydrophilic groups that can efficiently nucleate and remove water droplets. We constructed a cellulose-based superhydrophobic surface containing numerous olefin terminal groups by solvent exchange and spray coating. This surface is different from most of the reported biomimicking water harvesting surfaces that rely on complicated lithography and micropatterning techniques requiring special instruments. The obtained superhydrophobic surface was further modified using various thiol compounds via a thiol–ene reaction to manipulate the water harvesting property. The modified surfaces containing hydrophobic groups (e.g., 1-octadecanethiol and 1H,1H,2H,2H-perfluorodecanethiol) or a strong hydrophilic group (e.g., 3-mercaptopropionic acid and 6-mercapto-1-hexanol) exhibited insufficient fog collecting abilities due to poor water droplet nucleation or strong water adhesion. By contrast, the modified surface decorated with moderately hydrophilic amino groups combines the advantages of biological surfaces with distinct wetting features (such as fog-harvesting beetles and water-repellent lotus leaves), resulting in accelerated water nucleation and less compromise of the water removal efficiency. Molecular dynamic simulations revealed that the efficient droplet nucleation is attributed to the hydrophilic amino groups whereas the rapid droplet removal is due to the maintained superhydrophobicity of the amino group-modified surface. This strategy of decorating a superhydrophobic surface with moderately hydrophilic functional groups provides insight into the manipulation of droplet nucleation and removal for water collection efficiency.

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.

Synthesizing polymers with tailor-made molecular weight distribution (MWD) is an essential step toward better control and design of functional polymer materials. We propose a novel one-pot reaction strategy that can facilitate the inverse design of the shape, breadth, and skew of the MWD in a controlled polymerization. By a multistep initiator addition scheme that involves a sequence of addition operations with determined amounts of initiators and addition times, the polymers with target MWD can be possibly synthesized. This strategy is in principle suitable to reproduce most target MWDs, even with a multimodal profile or with large breadth and/or skew. As compared to previous relevant methods by blending polymers with determined molar masses and ratios, this one-pot reaction strategy avoids most of the tedious intermediate steps for controlling precision in blending and would be more convenient and timesaving. Our study supplies the inspiration of better control of synthesizing polymers with designated MWD in controlled polymerization.