Optimal face-to-face coupling for fast self-folding kirigami.

Kirigami-inspired designs can enable self-folding three-dimensional materials from flat, two-dimensional sheets. Hierarchical designs of connected levels increase the diversity of possible target structures, yet they can lead to longer folding times in the presence of fluctuations. Here, we study the effect of rotational coupling between levels on the self-folding of two-level kirigami designs driven by thermal noise in a fluid. Naturally present due to hydrodynamic resistance, we find that this coupling parameter can significantly impact a structure's self-folding pathway, thus enabling us to assess the quality of a kirigami design and the possibility for its optimization in terms of its folding rate and yield.

General Entropy Approach Toward Ultratough Sustainable Plastics.

Entropy is a universal concept across the physics of mixtures. While the role of entropy in other multicomponent materials has been appreciated, its effects in polymers and plastics have not. In this work, it is demonstrated that the seemingly small mixing entropy contributes to the miscibility and performance of polymer alloys. Experimental and modeling studies on over 30 polymer pairs reveal a strong correlation between entropy, morphology, and mechanical properties, while elucidating the mechanism behind: in polymer blends with weak interactions, entropy leads to homogeneously dispersed nanosized domains stabilized by highly entangled chains. This unique microstructure promotes uniform plastic deformation at the interface, thus improving the toughness of conventional brittle polymers by 1-2 orders of magnitude without sacrificing other properties, analogous to high-entropy metallic alloys. The proposed strategy also applies to ternary polymer systems and copolymers, offering a new pathway toward the development of sustainable polymers.

The mechanical and thermal properties of MoS2-WSe2 lateral heterostructures.

The recently fabricated monolayer MoS2-WSe2 lateral heterostructures are promising for many interesting applications, such as p-n diodes, photodetectors, transistors, sensors, light-emitting diodes and thermoelectric and flexible nanodevices. In this work, we study the mechanical and thermal properties of MoS2-WSe2 lateral heterostructures by using molecular dynamics (MD) simulations based on the recently parameterized Stillinger-Weber (SW) potential. It is found that the fracture strength and fracture strain of MoS2-WSe2 lateral heterostructures are dictated by the mechanical properties of MoS2, and are very sensitive to temperature. However, when a crack is introduced in the MoS2-WSe2 heterostructures, failure may occur either in MoS2 or WSe2, depending on the crack length and location. Interestingly, the fracture strengths obtained from our MD simulations are in agreement with those obtained from the Griffith theory. Our MD simulations further reveal that, in addition to the low thermal conductivities of MoS2 and WSe2, the MoS2-WSe2 heterojunctions exhibit a very low interfacial thermal conductance, which is about one order of magnitude lower than that of graphene-hBN heterojunctions.

Thermal conductivity of a h-BCN monolayer.

A hexagonal graphene-like boron-carbon-nitrogen (h-BCN) monolayer, a new two-dimensional (2D) material, has been synthesized recently. Herein we investigate for the first time the thermal conductivity of this novel 2D material. Using molecular dynamics simulations based on the optimized Tersoff potential, we found that the h-BCN monolayers are isotropic in the basal plane with close thermal conductivity magnitudes. Though h-BCN has the same hexagonal lattice as graphene and hexagonal boron nitride (h-BN), it exhibits a much lower thermal conductivity than the latter two materials. In addition, the thermal conductivity of h-BCN monolayers is found to be size-dependent but less temperature-dependent. Modulation of the thermal conductivity of h-BCN monolayers can also be realized by strain engineering. Compressive strain leads to a monotonic decrease in the thermal conductivity while the tensile strain induces an up-then-down trend in the thermal conductivity. Surprisingly, the small tensile strain can facilitate the heat transport of the h-BCN monolayers.

Thermal stability and thermal conductivity of phosphorene in phosphorene/graphene van der Waals heterostructures.

Phosphorene, a new two-dimensional (2D) semiconducting material, has attracted tremendous attention recently. However, its structural instability under ambient conditions poses a great challenge to its practical applications. A possible solution for this problem is to encapsulate phosphorene with more stable 2D materials, such as graphene, forming van der Waals heterostructures. In this study, using molecular dynamics simulations, we show that the thermal stability of phosphorene in phosphorene/graphene heterostructures can be enhanced significantly. By sandwiching phosphorene between two graphene sheets, its thermally stable temperature is increased by 150 K. We further study the thermal transport properties of phosphorene and find surprisingly that the in-plane thermal conductivity of phosphorene in phosphorene/graphene heterostructures is much higher than that of the free-standing one, with a net increase of 20-60%. This surprising increase in thermal conductivity arises from the increase in phonon group velocity and the extremely strong phonon coupling between phosphorene and the graphene substrate. Our findings have an important meaning for the practical applications of phosphorene in nanodevices.

Graphene membranes with nanoslits for seawater desalination via forward osmosis.

Stacked graphene (GE) membranes with cascading nanoslits can be synthesized economically compared to monolayer nanoporous GE membranes, and have potential for molecular separation. This study focuses on investigating the seawater desalination performance of these stacked GE layers as forward osmosis (FO) membranes by using molecular dynamics simulations. The FO performance is evaluated in terms of water flux and salt rejection and is explained by analysing the water density distribution and radial distribution function. The water flow displays an Arrhenius type relation with temperature and the activation energy for the stacked GE membrane is estimated to be 8.02 kJ mol-1, a value much lower than that of commercially available FO membranes. The study reveals that the membrane characteristics including the pore width, offset, interlayer separation distance and number of layers have significant effects on the desalination performance. Unlike monolayer nanoporous GE membranes, at an optimum layer separation distance, the stacked GE membranes with large pore widths and completely misaligned pore configuration can retain complete ion rejection and maintain a high water flux. Findings from the present study are helpful in developing GE-based membranes for seawater desalination via FO.

Atomic vacancies significantly degrade the mechanical properties of phosphorene.

Due to low formation energies, it is very easy to create atomic defects in phosphorene during its fabrication process. How these atomic defects affect its mechanical behavior, however, remain unknown. Here, we report on a systematic study of the effect of atomic vacancies on the mechanical properties and failure behavior of phosphorene using molecular dynamics simulations. It is found that atomic vacancies induce local stress concentration and cause early bond-breaking, leading to a significant degradation of the mechanical properties of the material. More specifically, a 2% concentration of randomly distributed mono-vacancies is able to reduce the fracture strength by ∼40%. An increase in temperature from 10 to 400 K can further deteriorate the fracture strength by ∼60%. The fracture strength of defective phosphorene is also found to be affected by defect distribution. When the defects are patterned in a line, the reduction in fracture strength greatly depends on the tilt angle and the loading direction. Furthermore, we find that di-vacancies cause an even larger reduction in fracture strength than mono-vacancies when the loading is in an armchair direction. These findings provide important guidelines for the structural design of phosphorene in future applications.

In-plane and cross-plane thermal conductivities of molybdenum disulfide.

We investigate the in-plane and cross-plane thermal conductivities of molybdenum disulfide (MoS2) using non-equilibrium molecular dynamics simulations. We find that the in-plane thermal conductivity of monolayer MoS2 is about 19.76 W mK(-1). Interestingly, the in-plane thermal conductivity of multilayer MoS2 is insensitive to the number of layers, which is in strong contrast to the in-plane thermal conductivity of graphene where the interlayer interaction strongly affects the in-plane thermal conductivity. This layer number insensitivity is attributable to the finite energy gap in the phonon spectrum of MoS2, which makes the phonon-phonon scattering channel almost unchanged with increasing layer number. For the cross-plane thermal transport, we find that the cross-plane thermal conductivity of multilayer MoS2 can be effectively tuned by applying cross-plane strain. More specifically, a 10% cross-plane compressive strain can enhance the thermal conductivity by a factor of 10, while a 5% cross-plane tensile strain can reduce the thermal conductivity by 90%. Our findings are important for thermal management in MoS2 based nanodevices and for thermoelectric applications of MoS2.

Failure modes and mechanisms of layered h-BN under local energy injection.

Layered h-BN may serve as an important dielectric and thermal management material in the next-generation nanoelectronics, in which its interactions with electron beam play an important role in device performance and reliability. Previous studies report variations in the failure strength and mode. In this study, using molecular dynamics simulations, we study the effect of local heat injection due to the electron beam and h-BN interaction on the failure start time and failure mode. It is found that at the same heat injection rate, the failure start time decreases with the increase in the layer number. With the introduction of point defects in the heating zone, the failure always starts from the defect site, and the start time can be significantly shortened. For monolayer h-BN, failure always starts within the layer, and once failure starts, its propagation is through melting or vaporization of the h-BN atoms, and no swelling occurs. For multiple layers, once failure starts within the h-BN film, swelling occurs first. With continued heating, the large pressure induced by melting and vaporization can cause the burst of the layers above, leading to the formation of a pit. In the presence of multiple defects within the heating zone, these defects can interact, causing a further reduction in the failure start time. We also reveal the relation of beam power with layer-by-layer failure mode and swelling/pit formation mode. The present work not only reproduces many interesting experimental observations, but also reveal several interesting mechanisms responsible for the failure processes and modes. It is expected that the findings revealed here may provide useful references for the design and engineering of h-BN for device applications.

Mechanical properties of methyl functionalized graphene: a molecular dynamics study.

Molecular dynamics simulations have been performed to study the mechanical properties of methyl (CH(3)) functionalized graphene. It is found that the mechanical properties of functionalized graphene greatly depend on the location, distribution and coverage of CH(3) radicals on graphene. Surface functionalization exhibits a much stronger influence on the mechanical properties than edge functionalization. For patterned functionalization on graphene surfaces, the radicals arranged in lines perpendicular to the tensile direction lead to larger strength deterioration than those parallel to the tensile direction. For random functionalization, the elastic modulus of graphene decreases gradually with increasing CH(3) coverage, while both the strength and fracture strain show a sharp drop at low coverage. When CH(3) coverage reaches saturation, the elastic modulus, strength and fracture strain of graphene drop by as much as 18%, 43% and 47%, respectively.

Exploring the structure-property relationship of three-dimensional hexagonal boron nitride aerogels with gyroid surfaces.

Three-dimensional hexagonal boron nitride aerogels (hBNAGs) are novel porous materials with many promising applications such as energy storage, thermal insulation and sensing. However, the structure-property relationships of hBNAGs in complicated thermo-mechanical coupled environments are still not clear. In this study, we employed a binary phase-field crystal (PFC) model to construct the atomic structures of hBNAGs, upon which the mechanical and thermal behaviors of hBNAGs were systematically investigated using large-scale atomistic simulations. It is found that the hBNAG geometry and topological defects strongly affect the mechanical and thermal properties. For example, the Young's modulus and tensile strength follow the scaling laws of mass density with a power factor of about 1.4 and 1.2, respectively, indicating that the stretching and bending combine toward tensile deformation. In addition, cracks nucleate around the octagon defects, indicating that the tensile strength is also influenced by the topological defects. Under compression, complicated crumpled deformations and ridges in the entire region are observed and the compression strength follows the scaling law of mass density with a power factor above 2.0, which means that a large portion of the hBNAGs do not contribute to the compression load bearing. We find that hBNAGs have a very low thermal conductivity of about two orders of magnitude lower than that of a hBN sheet. Also, the thermal conductivity of hBNAGs increases with increasing mass density, which also follows a scaling law. The power of the scaling law is about 0.5, indicating that the thermal conductivity has a strong nonlinear dependence on the mass density. Our work provides a deep understanding of the structure-property relationships of hBNAGs, which is useful for the engineering applications of hBNAGs.

Molecular-dynamics studies of competitive replacement in peptide-nanotube assembly for control of drug release.

We report molecular-dynamics simulation of carbon-nanotube-based drug delivery and release systems. We show that a peptide encapsulated inside or attached to the outer surface of a carbon nanotube can be released by another nanotube through a competitive replacement process. Energy analysis reveals that the van der Waals interaction plays the key role in this process, and the potential well between two nanotubes drives the competitive replacement. We further show that competitive replacement is a basic principle which may be generally explored for drug release. For example, one type of peptide can be used to replace/release another type of peptide, depending on the difference in their affinity for the nanotube. The effects of the peptide sequence and the nanotube size on the drug release process are also studied in this paper.

Large diffusion anisotropy and orientation sorting of phosphorene nanoflakes under a temperature gradient.

We perform molecular dynamics simulations to investigate the motion of phosphorene nanoflakes on a large graphene substrate under a thermal gradient. It is found that the atomic interaction between the graphene substrate and the phosphorene nanoflake generates distinct rates of motion for phosphorene nanoflakes with different orientations. Remarkably, for square phosphorene nanoflakes, the motion of zigzag-oriented nanoflakes is 2-fold faster than those of armchair-oriented and randomly-oriented nanoflakes. This large diffusion anisotropy suggests that sorting of phosphorene nanoflakes into specific orientations can be realized by a temperature gradient. The findings here provide interesting insights into strong molecular diffusion anisotropy and offer a novel route for manipulating two-dimensional materials.

Some Aspects of Thermal Transport across the Interface between Graphene and Epoxy in Nanocomposites.

Owing to the superior thermal properties of graphene, graphene-reinforced polymer nanocomposites hold great potential as the thermal interface materials (TIMs) dissipating heat for electronic packages. However, this application is greatly hindered by the high thermal resistance at the interface between graphene and polymer. In this paper, some important aspects of the improvement of the thermal transport across the interface between graphene and epoxy in graphene-epoxy nanocomposites, including the effectiveness of covalent and noncovalent functionalization, isotope doping, and acetylenic linkage in graphene are systematically investigated using molecular dynamics (MD) simulations. The simulation results show that the covalent and noncovalent functionalization techniques could considerably reduce the graphene-epoxy interfacial thermal resistance in the nanocomposites. Among different covalent functional groups, butyl is more effective than carboxyl and hydroxyl in reducing the interfacial thermal resistance. Different noncovalent functional molecules, including 1-pyrenebutyl, 1-pyrenebutyric acid, and 1-pyrenebutylamine, yield a similar amount of reductions. Moreover, it is found that the graphene-epoxy interfacial thermal resistance is insensitive to the carbon isotope doping in graphene, while it can be reduced moderately by replacing the sp(2) bonds in graphene with acetylenic linkages.

Thermal conductivities of single- and multi-layer phosphorene: a molecular dynamics study.

As a new two-dimensional (2D) material, phosphorene has drawn growing attention owing to its novel electronic properties, such as layer-dependent direct bandgaps and high carrier mobility. Herein we investigate the in-plane and cross-plane thermal conductivities of single- and multi-layer phosphorene, focusing on geometrical (sample size, orientation and layer number) and strain (compression and tension) effects. A strong anisotropy is found in the in-plane thermal conductivity with its value along the zigzag direction being much higher than that along the armchair direction. Interestingly, the in-plane thermal conductivity of multi-layer phosphorene is insensitive to the layer number, which is in strong contrast to that of graphene where the interlayer interactions strongly influence the thermal transport. Surprisingly, tensile strain leads to an anomalous increase in the in-plane thermal conductivity of phosphorene, in particular in the armchair direction. Both the in-plane and cross-plane thermal conductivities can be modulated by external strain; however, the strain modulation along the cross-plane direction is more effective and thus more tunable than that along the in-plane direction. Our findings here are of great importance for the thermal management in phosphorene-based nanoelectronic devices and for thermoelectric applications of phosphorene.