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.

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.

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.

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.

Cross-scale biophysics modelling of F-actin cytoskeleton in cell

This thesis is a comparative study of the modelling of mechanical behaviours of F-actin cytoskeleton which is an important structural component in living cells. A new granular model was developed for F-actin cytoskeleton based on the concept of multiscale modelling. This framework overcomes difficulties encountered in physical modelling of cytoskeleton in conventional continuum mechanics modelling, and the computational challenges in all-atom molecular dynamics simulation. The thermostat algorithm was further modified to better predict the thermodynamic properties of F-actin cytoskeleton in modelling. This multiscale modelling framework was applied in explaining the physical mechanisms of cytoskeleton responses to external mechanical loads.

Graphene-covered perovskites: An effective strategy to enhance light absorption and resist moisture degradation

The long-term stability of methylammonium lead triiodide (MAPbI3) perovskite in moist environments is a paramount challenge to realise the commercialization of perovskite solar cells. In an attempt to address this concern, we have carried out systematic first-principles studies on the MAPbI3 perovskite with a hydrophobic graphene layer interfaced as a water barrier. We find there is a charge transfer at the graphene/MAPbI3 interface and electrons can be excited from graphene into the perovskite surface, leading to well separated electron–hole pairs, i.e. reduced recombination. By studying the optical properties, we find the hybrid graphene/MAPbI3 nanocomposite displays enhanced light absorption compared with the pristine MAPbI3. Furthermore, from an ab initio molecular dynamics simulation, the graphene/MAPbI3 nanocomposite is confirmed to be able to resist the reaction with water molecules, highlighting a great advantage of this nanocomposite in promoting long-term photovoltaic performance.

Diffusion of flexible, charged, nanoscopic molecules in solution: Size and pH dependence for PAMAM dendrimer

In order to understand self-diffusion (D) of a charged, flexible, and porous nanoscopic molecule in water, we carry out very long, fully atomistic molecular dynamics simulation of PAMAM dendrimer up to eight generations in explicit salt water under varying pH. We find that while the radius of gyration (R-g) varies as N-1/3, the self-diffusion constant (D) scales, surprisingly, as N-alpha, with alpha=0.39 at high pH and 0.5 at neutral pH, indicating a dramatic breakdown of Stokes-Einstein relation for diffusion of charged nanoscopic molecules. The variation in D as a function of radius of gyration demonstrates the importance of treating water and ions explicitly in the diffusion process of a flexible nanoscopic molecule. In agreement with recent experiments, the self-diffusion constant increases with pH, revealing the importance of dielectric friction in the diffusion process. The shape of a dendrimer is found to fluctuate on a nanosecond time scale. We argue that this flexibility (and also the porosity) of the dendrimer may play an important role in determining the mean square displacement of the dendrimer and the breakdown of the Stokes-Einstein relation between diffusion constant and the radius.

Correlation between conductivity or diffusivity and activation energy in amorphous solids

There exist many investigations of ionic transport in a variety of glasses. These studies exhibit strong correlation between ionic conductivity and activation energy: Typically, it is found that higher conductivity is associated with lower activation energies and vice versa. Although there are explanations for this at a phenomenological level, there is no consistent physical picture to explain the correlation between conductivity and activation energy. We have carried out molecular dynamics simulation as a function of the size of the impurity atom or diffusant (both neutral and charged) in a host amorphous matrix. We find that there is a maximum in self-diffusivity as a function of the size of the impurity atom suggesting that there is an appropriate size for which the diffusivity is maximum. The activation energy is found to be the lowest for this size of the impurity. A similar maximum has been previously found in other condensed phases, such as confined fluids and dense liquids, and has its origin in the levitation effect. The implications of this result for understanding ionic conductivity in glasses are discussed. Our results suggest that there is a relation between microscopic structure of the amorphous solid, diffusivity or conductivity, and activation energy. The nature of this relationship is discussed in terms of the levitation parameter showing that diffusivity is maximum when the size of the neck or doorway radius is comparable with the size of the diffusant. Our computational results here are in excellent agreement with independent experimental results of Nascimento et al. [Braz. J. Phys. 35, 626 (2005)] that structural features of the glass are important in determining the ionic conductivity.

Graphene-like two-dimensional ionic boron with double dirac cones at ambient condition

Recently, partially ionic boron (γ-B28) has been predicted and observed in pure boron, in bulk phase and controlled by pressure [Nature, 457 (2009) 863]. By using ab initio evolutionary structure search, we report the prediction of ionic boron at a reduced dimension and ambient pressure, namely, the two-dimensional (2D) ionic boron. This 2D boron structure consists of graphene-like plane and B2 atom pairs, with the P6/mmm space group and 6 atoms in the unit cell, and has lower energy than the previously reported α-sheet structure and its analogues. Its dynamical and thermal stability are confirmed by the phonon-spectrum and ab initio molecular dynamics simulation. In addition, this phase exhibits double Dirac cones with massless Dirac fermions due to the significant charge transfer between the graphene-like plane and B2 pair that enhances the energetic stability of the P6/mmm boron. A Fermi velocity (vf) as high as 2.3 x 106 m/s, which is even higher than that of graphene (0.82 x 106 m/s), is predicted for the P6/mmm boron. The present work is the first report of the 2D ionic boron at atmospheric pressure. The unique electronic structure renders the 2D ionic boron a promising 2D material for applications in nanoelectronics.

The 5-Methyl Group in Thymine Dynamically Influences the Structure of A-Tracts in DNA at the Local and Global Level

DNA sequences containing a stretch of several A:T basepairs without a 5'-TA-3' step are known as A-tracts and have been the subject of extensive investigation because of their unique structural features such as a narrow minor groove and their crucial role in several biological processes. One of the aspects under investigation has been the influence of the 5-methyl group of thymine on the properties of A-tracts. Detailed molecular dynamics simulation studies of the sequences d(CGCAAAUUUGCG) and d(CGCAAATTTGCG) indicate that the presence of the 5-methyl group in thymine increases the frequency of a narrow minor groove conformation, which could facilitate its specific recognition by proteins, and reduce its susceptibility to cleavage by DNase I. The bias toward a wider minor groove in the absence of the thymine 5-methyl group is a static structural feature. Our results also indicate that the presence of the thymine 5-methyl group is necessary for calibrating the backbone conformation and the basepair and dinucleotide step geometry of the core A-tract as well as the flanking CA/TG and the neighboring GC/GC steps, as observed in free and protein-bound DNA. As a consequence, it also fine-tunes the curvature of the longer DNA fragment in which the A-tract is embedded.

Structural Polymorphism in d(T).d(A)*d(T) Triple Helices

DNA triple helices containing two thymine strands and one adenine strand have been studied, using model building followed by energy minimisation, for different orientations of the third strand resulting from variation in the hydrogen bonding between the Watson-Crick duplex and the third strand and the glycosidic torsion angle in the third strand. Our results show that the structure with a parallel orientation of the third strand, in which the third strand base forms Hoogsteen hydrogen bonds with the adenine base in the Watson-Crick duplex, is energetically the most favourable. An antiparallel orientation of the third strand is also possible, in which the third strand base hydrogen bonds to both the bases in the Watson-Crick duplex. This structure is energetically comparable to the parallel structure. For the parallel triplex a 200ps molecular dynamics simulation starting from two different starting structures indicates that at 300K significant structural heterogeneity exists in this tripler structure. The results are compared with existing structural data on this class of triplexes derived from theoretical and NMR techniques.

Study of translational and rotational mobility and orientational preference of ethane in one-dimensional channels

Structural and dynamical properties of ethane in one-dimensional channels of AlPO4-5 and carbon nanotube have been investigated at dilute concentration with the help of molecular dynamics simulation. Density distributions and orientational structure of ethane have been analyzed. Repulsive interactions seem to play an important role when ethane is located in the narrow part of the AlPO4-5 channel. In AlPO4-5, parallel orientation is predominant over perpendicular orientation except when ethane is located in the broader part of the channel. Unlike in the case of single-file diffusion, our results in carbon nanotube show that at dilute concentrations the mean squared displacement, mu(2)(t) approximate to t(alpha), alpha = 1.8. The autocorrelation function for the z-component of angular velocity of ethane in space-fixed frame of reference shows a pronounced negative correlation. This is attributed to the restriction in the movement of ethane along the x- and y- directions. It is seen that the ratio of reorientational correlation times does not follow the Debye model for confined ethane but it is closer to the predictions of the Debye model for bulk ethane.

Allosteric Communication in Cysteinyl tRNA Synthetase A NETWORK OF DIRECT AND INDIRECT READOUT

Protein structure networks are constructed for the identification of long-range signaling pathways in cysteinyl tRNA synthetase (CysRS). Molecular dynamics simulation trajectory of CysRS-ligand complexes were used to determine conformational ensembles in order to gain insight into the allosteric signaling paths. Communication paths between the anticodon binding region and the aminoacylation region have been identified. Extensive interaction between the helix bundle domain and the anticodon binding domain, resulting in structural rigidity in the presence of tRNA, has been detected. Based on the predicted model, six residues along the communication paths have been examined by mutations (single and double) and shown to mediate a coordinated coupling between anticodon recognition and activation of amino acid at the active site. This study on CysRS clearly shows that specific key residues, which are involved in communication between distal sites in allosteric proteins but may be elusive in direct structure analysis, can be identified from dynamics of protein structure networks.