Abstract: The buckling of defected rectangular graphene nanoribbons when subjected to axial stress with supported boundary conditions is investigated using atomistic simulations. The buckling strain and mechanical stiffness of monolayer graphene decrease with the percentage of randomly distributed vacancies. The elasticity to plasticity transition in the stress-strain curve, at low percentage of vacancies, are found to be almost equal to the buckling strain thresholds and they decrease with increasing percentage of vacancies.

Abstract: The buckling of graphene nano-ribbons containing a grain boundary is studied using atomistic simulations where free and supported boundary conditions are invoked. We consider the buckling transition of two kinds of grain boundaries with special symmetry. When graphene contains a large angle grain boundary with theta = 21.8 degrees, the buckling strains are larger than those of perfect graphene when the ribbons with free (supported) boundary condition are subjected to compressive tension parallel (perpendicular) to the grain boundary. This is opposite for the results of theta = 32.2 degrees. The shape of the deformations of the buckled graphene nanoribbons depends on the boundary conditions, the presence of the particular used grain boundaries, and the direction of applied in-plane compressive tension. (C) 2012 American Institute of Physics. [http://dx.doi.org/10.1063/1.3692573]

Abstract: The electronic properties of a triaxially strained hexagonal graphene flake with either armchair or zigzag edges are investigated using molecular dynamics simulations and tight-binding calculations. We found that (i) the pseudomagnetic field in strained graphene flakes is not uniform neither in the center nor at the edge of zigzag terminated flakes, (ii) the pseudomagnetic field is almost zero in the center of armchair terminated flakes but increases dramatically near the edges, (iii) the pseudomagnetic field increases linearly with strain, for strains lower than 15% but increases nonlinearly beyond it, (iv) the local density of states in the center of the zigzag hexagon exhibits pseudo-Landau levels with broken sublattice symmetry in the zeroth pseudo-Landau level, and in addition there is a shift in the Dirac cone due to strain induced scalar potentials, and (v) there is size effect in pseudomagnetic field. This study provides a realistic model of the electronic properties of inhomogeneously strained graphene where the relaxation of the atomic positions is correctly included together with strain induced modifications of the hopping terms up to next-nearest neighbors.

Abstract: The spatial dependence of the van der Waals (vdW) energy between graphene and hexagonal boron-nitride (h-BN) is investigated using atomistic simulations. The van der Waals energy between graphene and h-BN shows a hexagonal superlattice structure identical to the observed Moire pattern in the local density of states, which depends on the lattice mismatch and misorientation angle between graphene and h-BN. Our results provide atomistic features of the weak van der Waals interaction between graphene and BN which are in agreement with experiment and provide an analytical expression for the size of the spatial variation of the weak van der Waals interaction. We also found that the A-B-lattice symmetry of graphene is broken along the armchair direction. (C) 2014 AIP Publishing LLC.

Abstract: Ripples in pristine freestanding graphene naturally orient themselves in an array that is alternately curved-up and curved-down; maintaining an average height of zero. Using scanning tunneling microscopy (STM) to apply a local force, the graphene sheet will reversibly rise and fall in height until the height reaches 60%-70% of its maximum at which point a sudden, permanent jump occurs. We successfully model the ripples as a spin-half Ising magnetic system, where the height of the graphene plays the role of the spin. The permanent jump in height, controlled by the tunneling current, is found to be equivalent to an antiferromagnetic-to-ferromagnetic phase transition. The thermal load underneath the STM tip alters the local tension and is identified as the responsible mechanism for the phase transition. Four universal critical exponents are measured from our STM data, and the model provides insight into the statistical role of graphene's unusual negative thermal expansion coefficient.

Abstract: Using atomistic simulations, we determine the roughness and the thermal properties of a suspended graphane sheet. As compared to graphene, we found that (i) hydrogenated graphene has a larger thermal contraction, (ii) the roughness exponent at room temperature is smaller, i.e., ≃ 1.0 versus ≃ 1.2 for graphene, (iii) the wavelengths of the induced ripples in graphane cover a wide range corresponding to length scales in the range 30125 Å at room temperature, and (iv) the heat capacity of graphane is estimated to be 29.32±0.23 J/mol K, which is 14.8% larger than that for graphene, i.e., 24.98±0.14 J/mol K. Above 1500 K, we found that graphane buckles when its edges are supported in the x-y plane.

Abstract: Density-functional tight-binding and classical molecular dynamics simulations are used to investigate the structural deformations and melting of planar carbon nanoclusters C-N with N = 2-55. The minimum-energy configurations for different clusters are used as starting configurations for the study of the temperature effects on the bond breaking and rotation in carbon lines (N < 6), carbon rings (5 < N < 19), and graphene nanoflakes. The larger the rings (graphene nanoflakes) the higher the transition temperature (melting point) with ring-to-line (perfect-to-defective) transition structures. The melting point was obtained by using the bond energy, the Lindemann criteria, and the specific heat. We found that hydrogen-passivated graphene nanoflakes (CNHM) have a larger melting temperature with a much smaller dependence on size. The edges in the graphene nanoflakes exhibit several different metastable configurations (isomers) during heating before melting occurs. DOI: 10.1103/PhysRevB.87.134103

Abstract: Epitaxial graphene is grown on a non-polar n(+) 6H-SiC m-plane substrate and studied using atomic scale scanning tunneling microscopy. Multilayer graphene is found throughout the surface and exhibits rotational disorder. Moire patterns of different spatial periodicities are found, and we found that as the wavelength increases, so does the amplitude of the modulations. This relationship reveals information about the interplay between the energy required to bend graphene and the interaction energy, i.e. van der Waals energy, with the graphene layer below. Our experiments are supported by theoretical calculations which predict that the membrane topographical amplitude scales with the Moire pattern wavelength, L as L-1 + alpha L-2. (C) 2014 Elsevier Ltd. All rights reserved.

Abstract: Nanoindentation of bilayer graphene is studied using molecular-dynamics simulations. We compared our simulation results with those from elasticity theory as based on the nonlinear Föppl-Hencky equations with rigid boundary condition. The force-deflection values of bilayer graphene are compared to those of monolayer graphene. Youngs modulus of bilayer graphene is estimated to be 0.8 TPa which is close to the value for graphite. Moreover, an almost flat bilayer membrane at low temperature under central load has a 14% smaller Youngs modulus as compared to the one at room temperature.

Abstract: Large-scale atomistic simulations using the reactive empirical bond order force field approach is implemented to investigate thermal and mechanical properties of single-layer (SL) and multilayer (ML) molybdenum disulfide (MoS2). The amplitude of the intrinsic ripples of SL MoS2 are found to be smaller than those exhibited by graphene (GE). Furthermore, because of the van der Waals interaction between layers, the out-of-plane thermal fluctuations of ML MoS2 decreases rapidly with increasing number of layers. This trend is confirmed by the buckling transition due to uniaxial stress which occurs for a significantly larger applied tension as compared to graphene. For SL MoS2, the melting temperature is estimated to be 3700 K which occurs through dimerization followed by the formation of small molecules consisting of two to five atoms. When different types of vacancies are inserted in the SL MoS2 it results in a decrease of both the melting temperature as well as the stiffness.

Abstract: Freestanding silicene is predicted to display comparable electronic properties as graphene. However, the yet synthesized silicenelike structures have been only realized on different substrates which turned out to exhibit versatile crystallographic structures that are very different from the theoretically predicted buckled phase of freestanding silicene. This calls for a different approach where silicene is stabilized using very weakly interacting surfaces. We propose here a route by using graphene bilayer as a scaffold. The confinement between the flat graphene layers results in a planar clustering of Si atoms with small buckling, which is energetically unfavorable in vacuum. Buckled hexagonal arrangement of Si atoms similar to freestanding silicene is observed for large clusters, which, in contrast to Si atoms on metallic surfaces, is only very weakly van der Waals coupled to the graphene layers. These clusters are found to be stable well above room temperature. Our findings, which are supported by density-functional tight-binding calculations, show that intercalating bilayer graphene with Si is a favorable route to realize silicene.

Abstract: The strain-induced pseudomagnetic field in supported graphene deposited on top of a nanostructured substrate is investigated by using atomistic simulations. A step, an elongated trench, a one-dimensional barrier, a spherical bubble, a Gaussian bump, and a Gaussian depression are considered as support structures for graphene. From the obtained optimum configurations we found very strong induced pseudomagnetic fields which can reach up to similar to 1000 T due to the strain-induced deformations in the supported graphene. Different magnetic confinements with controllable geometries are found by tuning the pattern of the substrate. The resulting induced magnetic fields for graphene on top of a step, barrier, and trench are calculated. In contrast to the step and trench the middle part of graphene on top of a barrier has zero pseudomagnetic field. This study provides a theoretical background for designing magnetic structures in graphene by nanostructuring substrates. We found that altering the radial symmetry of the deformation changes the sixfold symmetry of the induced pseudomagnetic field.

Abstract: Thermal fluctuations of single layer hydrogenated graphene (graphane) are investigated using large scale atomistic simulations. By analyzing the mean square value of the height fluctuations < h(2)> and the height-height correlation function H(q) for different system sizes and temperatures, we show that hydrogenated graphene is an unrippled system in contrast to graphene. The height fluctuations are bounded, which is confirmed by a H(q) tending to a constant in the long wavelength limit instead of showing the characteristic scaling law q(4-eta)(eta similar or equal to 0.85) predicted by membrane theory. This unexpected behavior persists up to temperatures of at least 900 K and is a consequence of the fact that in graphane the thermal energy can be accommodated by in-plane bending modes, i.e., modes involving C-C-C bond angles in the buckled carbon layer, instead of leading to significant out-of-plane fluctuations that occur in graphene.

Abstract: Using the valence force field model of Perebeinos and Tersoff (2009 Phys. Rev. B 79 241409(R)), different energy modes of suspended graphene subjected to tensile or compressive strain are studied. By carrying out Monte Carlo simulations it is found that: (i) only for small strains (vertical bar epsilon vertical bar (sic) 0.02) is the total energy symmetrical in the strain, while it behaves completely differently beyond this threshold; (ii) the important energy contributions in stretching experiments are stretching, angle bending, an out-of-plane term, and a term that provides repulsion against pi-pi misalignment; (iii) in compressing experiments the two latter terms increase rapidly, and beyond the buckling transition stretching and bending energies are found to be constant; (iv) from stretching-compressing simulations we calculated the Young's modulus at room temperature 350 +/- 3.15 N m(-1), which is in good agreement with experimental results (340 +/- 50 N m(-1)) and with ab initio results (322-353) N m(-1); (v) molar heat capacity is estimated to be 24.64 J mol(-1) K-1 which is comparable with the Dulong-Petit value, i. e. 24.94 J mol(-1) K-1, and is almost independent of the strain; (vi) nonlinear scaling properties are obtained from height-height correlations at finite temperature; (vii) the used valence force field model results in a temperature independent bending modulus for graphene, and (viii) the Gruneisen parameter is estimated to be 0.64.

Abstract: Using first-principles calculations and reactive force field molecular dynamics simulations, we study the structural properties and dynamics of a fluorine (F) atom, either adsorbed on the surface of single layer graphene (F/GE) or between the layers of AB stacked bilayer graphene (F@ bilayer graphene). It is found that the diffusion of the F atom is very different in those cases, and that the mobility of the F atom increases by about an order of magnitude when inserted between two graphene layers. The obtained diffusion constant for F/GE is twice larger than that experimentally found for gold adatom and theoretically found for C-60 molecule on graphene. Our study provides important physical insights into the dynamics of fluorine atoms between and on graphene layers and explains the mechanism behind the separation of graphite layers due to intercalation of F atoms.