Graphene ripples generated by grain boundaries in highly ordered pyrolytic graphite

Atomic scale periodic ripples that extend for several nanometers on the surface of adjacent graphitic grains have been observed for the first time on highly ordered pyrolitic graphite by UHV-STM. The ripples emanate from a grain boundary, and are explained in terms of a mechanical deformation due to the elastic strain accumulated along the GB, which is relieved out-of-plane in the topmost graphene layer. We present a molecular dynamics model that accounts for the formation of similar ripples as result of the lattice mismatch induced by two different grain orientations.

Thermal conductivity of Si nanowires with faulted stacking layers

Faulted stacking layers are ubiquitously observed during the crystal growth of semiconducting nanowires (NWs). In this paper, we employ the reverse non-equilibrium molecular dynamics simulation to elucidate the effect of various faulted stacking layers on the thermal conductivity (TC) of silicon (Si) NWs. We find that the stacking faults can greatly reduce the TC of the Si NW. Among the different stacking faults that are parallel to the NW's axis, the 9R polytype structure, the intrinsic and extrinsic stacking faults (iSFs and eSFs) exert more pronounced effects in the reduction of TC than the twin boundary (TB). However, for the perpendicularly aligned faulted stacking layers, the eSFs and 9R polytype structures are observed to induce a larger reduction to the TC of the NW than the TB and iSFs. For all considered NWs, the TC does not show a strong relation with the increasing number of faulted stacking layers. Our studies suggest the possibility of tuning the thermal properties of Si NWs by altering the crystal structure via the different faulted stacking layers.

Tailoring the resonance of bilayer graphene sheets by interlayer sp3 bonds

Graphene-based resonators are envisioned to build the ultimate limit of two-dimensional nanoelectromechanical system due to their ultrasensitive detection of mass, force, pressure and charge. However, such application has been greatly impeded by their extremely low quality factor. In the present work, we explore, using the large-scale molecular dynamics simulation, the possibility of tailoring the resonance properties of a bilayer graphene sheet (GS) with interlayer sp3 bonds. For the bilayer GS resonator with interlayer sp3 bonds, we discovered that the sp3 bonds can either degrade or enhance the resonance properties of the resonator depending on their density and location. It is found that the distribution of sp3 bonds only along the edges of either pristine or hydrogenated bilayer GS, leads to a greatly enhanced quality factor. A quality factor of ~1.18×105 is observed for a 3.07×15.31 nm2 bilayer GS resonator with sp3 bonds, which is more than 30 times larger comparing with that of a pristine bilayer GS. The present study demonstrates that the resonance properties of a bilayer GS resonator can be tuned by introducing sp3 bonds. This finding provides a useful guideline for the synthesis of the bilayer GS for its application as a resonator component.

Atomistic simulation of interfacial behaviour in graphene-polymer nanocomposites

Graphene–polymer nanocomposites have promising properties as new structural and functional materials. The remarkable mechanical property enhancement in these nanocomposites is generally attributed to exceptional mechanical property of graphene and possible load transfer between graphene and polymer matrix. However, the underlying strengthening and toughening mechanisms have not been well understood. In this work, the interfacial behavior of graphene-polyethylene (PE) was investigated using molecular dynamics (MD) method. The interfacial shear force (ISF) and interfacial shear stress (ISS) between graphene and PE matrix were evaluated, taking into account graphene size, the number of graphene layers and the structural defects in graphene. MD results show that the ISS at graphene-PE interface mainly distributes at each end of the graphene nanofiller within the range of 1 nm, and much larger than that at carbon nanotube (CNT)-PE interface. Moreover, it was found that the ISS at graphene-PE interface is sensitive to the layer number.

Effect of defects on fracture strength of graphene sheets

With a hexagonal monolayer network of carbon atoms, graphene has demonstrated exceptional electrical 22 and mechanical properties. In this work, the fracture of graphene sheets with Stone–Wales type defects and vacancies were investigated using molecular dynamics simulations at different temperatures. The initiation of defects via bond rotation was also investigated. The results indicate that the defects and vacancies can cause significant strength loss in graphene. The fracture strength of graphene is also affected by temperature and loading directions. The simulation results were compared with the prediction from the quantized fracture mechanics.

A stochastic thermostat algorithm for coarse-grained thermomechanical modeling of large-scale soft matters : theory and application to microfilaments

As all-atom molecular dynamics method is limited by its enormous computational cost, various coarse-grained strategies have been developed to extend the length scale of soft matters in the modeling of mechanical behaviors. However, the classical thermostat algorithm in highly coarse-grained molecular dynamics method would underestimate the thermodynamic behaviors of soft matters (e.g. microfilaments in cells), which can weaken the ability of materials to overcome local energy traps in granular modeling. Based on all-atom molecular dynamics modeling of microfilament fragments (G-actin clusters), a new stochastic thermostat algorithm is developed to retain the representation of thermodynamic properties of microfilaments at extra coarse-grained level. The accuracy of this stochastic thermostat algorithm is validated by all-atom MD simulation. This new stochastic thermostat algorithm provides an efficient way to investigate the thermomechanical properties of large-scale soft matters.

Surface functionalization on the thermal conductivity of graphene–polymer nanocomposites

Exploring thermal transport in graphene-polymer nanocomposite is significant to its applications with better thermal properties. Interfacial thermal conductance between graphene and polymer matrix plays a critical role in the improvement of thermal conductivity of graphene-polymer nanocomposite. Unfortunately, it is still challenging to understand the interfacial thermal transport between graphene nanofiller and polymer matrix at small material length scale. To this end, using non-equilibrium molecular dynamics simulations, we investigate the interfacial thermal conductance of graphene-polyethylene (PE) nanocomposite. The influence of functionalization with hydrocarbon chains on the interfacial thermal conductance of graphene-polymer nanocomposites was studied, taking into account of the effects of model size and thermal conductivity of graphene. An analytical model is also used to calculate the thermal conductivity of nanocomposite. The results are considered to contribute to development of new graphene-polymer nanocomposites with tailored thermal properties.

Atomistic simulation of surface functionalization on the interfacial properties of graphene-polymer nanocomposites

Graphene has been increasingly used as nano sized fillers to create a broad range of nanocomposites with exceptional properties. The interfaces between fillers and matrix play a critical role in dictating the overall performance of a composite. However, the load transfer mechanism along graphene-polymer interface has not been well understood. In this study, we conducted molecular dynamics simulations to investigate the influence of surface functionalization and layer length on the interfacial load transfer in graphene polymer nanocomposites. The simulation results show that oxygen-functionalized graphene leads to larger interfacial shear force than hydrogen-functionalized and pristine ones during pull-out process. The increase of oxygen coverage and layer length enhances interfacial shear force. Further increase of oxygen coverage to about 7% leads to a saturated interfacial shear force. A model was also established to demonstrate that the mechanism of interfacial load transfer consists of two contributing parts, including the formation of new surface and relative sliding along the interface. These results are believed to be useful in development of new graphene-based nanocomposites with better interfacial properties.

Mechanics of microfilaments networks : a cross-scales study

The mechanical properties of microfilament networks are systematically summarized at different special scales in this paper. We have presented the mechanical models of single microfilaments and microfilament networks at microscale. By adopting a coarse-grained simulation strategy, the mechanical stability of microfilaments related cellular structures are analysed. Structural analysis is conducted to microfilament networks to understand the stress relaxation under compression. The nanoscale molecular mechanisms of the microfilaments deformation is also summarized from the viewpoint of molecular dynamics simulation. This paper provides the fundaments of multiscale modelling framework for the mechanical behaviours simulation of hierarchical microfilament networks.

Tuneable resonance properties of graphene by nitrogen-dopant

Doping as one of the popular methods to manipulate the properties of nanomaterials has received extensive application in deriving different types of graphene derivates, while the understanding of the resonance properties of dopant graphene is still lacking in literature. Based on the large-scale molecular dynamics simulation, reactive empirical bond order potential, as well as the tersoff potential, the resonance properties of N-doped graphene were studied. The studied samples were established according to previous experiments with the N atom’s percentage ranging from 0.43%-2.98%, including three types of N dopant locations, i.e., graphitic N, pyrrolic N and pyridinic N. It is found that different percentages of N-dopant exert different influence to the resonance properties of the graphene, while the amount of N-dopant is not the only factor that determines its impact. For all the considered cases, a relative large percentage of N-dopant (2.98% graphitic N-dopant) is observed to introduce significant influence to the profile of the external energy, and thus lead to an extremely low Q-factor comparing with that of the pristine graphene. The most striking finding is that, the natural frequency of the defective graphene with N-dopant appears uniformly larger than that of the pristine defective graphene. While for the perfect graphene, the N-dopant shows less influence to its natural frequency. This study will enrich the current understanding of the influence of dopants on graphene, which will eventually shed lights on the design of different molecules-doped graphene sheet.

Hydrodynamic theory of liquid slippage on a solid substrate near a moving contact line

In this Letter a hydrodynamic theory of liquid slippage on a solid substrate near a moving contact line is proposed. A family of spatially varying slip lengths in the Navier slip law recovers the results of past formulations for slip in continuum theories and molecular dynamics simulations and is consistent with well-established experimental observations of complete wetting. This formulation gives a general approach for continuum hydrodynamic theories. New fluid flow behaviors are also predicted yet to be seen in experiment. © 2013 American Physical Society.

Structure-mediated thermal transport of monolayer graphene allotropes nanoribbons

We report the study of the thermal transport management of monolayer graphene allotrope nanoribbons (size ∼20 × 4 nm2) by the modulation of their structures via molecular dynamics simulations. The thermal conductivity of graphyne (GY)-like geometries is observed to decrease monotonously with increasing number of acetylenic linkages between adjacent hexagons. Strikingly, by incorporating those GY or GY-like structures, the thermal performance of graphene can be effectively engineered. The resulting hetero-junctions possess a sharp local temperature jump at the interface, and show a much lower effective thermal conductivity due to the enhanced phonon–phonon scattering. More importantly, by controlling the percentage, type and distribution pattern of the GY or GY-like structures, the hetero-junctions are found to exhibit tunable thermal transport properties (including the effective thermal conductivity, interfacial thermal resistance and rectification). This study provides a heuristic guideline to manipulate the thermal properties of 2D carbon networks, ideal for application in thermoelectric devices with strongly suppressed thermal conductivity.

Engineering the mechanical properties of graphene nanotube hybrid structures through structural modulation

The excellent multi-functional properties of carbon nanotube (CNT) and graphene have enabled them as appealing building blocks to construct 3D carbon-based nanomaterials or nanostructures. The recently reported graphene nanotube hybrid structure (GNHS) is one of the representatives of such nanostructures. This work investigated the relationships between the mechanical properties of the GNHS and its structure basing on large-scale molecular dynamics simulations. It is found that increasing the length of the constituent CNTs, the GNHS will have a higher Young’s modulus and yield strength. Whereas, no strong correlation is found between the number of graphene layers and Young’s modulus and yield strength, though more graphene layers intends to lead to a higher yield strain. In the meanwhile, the presences of multi-wall CNTs are found to greatly strengthen the hybrid structure. Generally, the hybrid structures exhibit a brittle behavior and the failure initiates from the connecting regions between CNT and graphene. More interestingly, affluent formations of monoatomic chains and rings are found at the fracture region. This study provides an in-depth understanding of the mechanical performance of the GNHSs while varying their structures, which will shed lights on the design and also the applications of the carbon-based nanostructures.

Numerical investigation of graphene and graphene-polymer nanocomposites

This thesis is a comprehensive and deep investigation on graphene and graphene-polymer nanocomposites. It explores the strong structure-property relationships in both graphene and graphene-based polymeric nanocomposites. A number of significant conclusions, including failure mechanism in graphene, interfacial load transfer and thermal transport mechanisms in graphene-polymer nanocomposites, have been drawn through both atomistic simulations and theoretical analysis. These results can provide direct guidelines for development of new graphene-based materials and devices.

Hydrolysis of genotoxic methyl-substituted oxiranes : Experimental kinetic and semiempirical studies

The kinetics of acid-catalyzed hydrolysis of seven methylated aliphatic epoxides - R1R2C(O)CR3R4 (A: R1=R2=R3=R4=H; B: R1=R2=R3=H, R4=Me; C: R1=R2=H, R3=R4=Me; D: R1=R3=H, R2=R4=Me(trans); E: R1=R3=H, R2=R4=Me(cis); F: R1=R3=R4=Me, R2=H; G: R1=R2=R3=R4=Me) - has been studied at 36 ± 1.5°C. Compounds with two methyl groups at the same carbon atom of the oxirane ring exhibit highest rate constants (k(eff) in reciprocal molar concentration per second: 11.0 ± 1.3 for C, 10.7 ± 2.1 for F, and 8.7 ± 0.7 for G as opposed to 0.124 ± 0.003 for B, 0.305 ± 0.003 for D, and 0.635 ± 0.036 for E). Ethylene oxide (A) displays the lowest rate of hydrolysis (0.027 M-1 s-1). The results are consistent with literature data available for compounds A, B, and C. To model the reactivities we have employed quantum chemical calculations (MNDO, AM1, PM3, and MINDO/3) of the main reaction species. There is a correlation of the logarithm k(eff) with the total energy of epoxide ring opening. The best correlation coefficients (r) were obtained using the AM1 and MNDO methods (0.966 and 0.957, respectively). However, unlike MNDO, AM1 predicts approximately zero energy barriers for the oxirane ring opening of compounds B, C, E and G, which is not consistent with published kinetic data. Thus, the MNDO method provides a preferential means of modeling the acidic hydrolysis of the series of methylated oxiranes. The general ranking of mutagenicity in vitro, A > B > C, is in line with the concept that this sequence also gradually leaves the expoxide reactivity optimal for genotoxicity toward reactivities leading to higher biological detoxifications.