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: We investigate the stability of circular monolayer graphene subjected to a radial load using non-equilibrium molecular dynamics simulations. When monolayer graphene is radially stressed, after some small circular strain (~0.4%) it buckles and bends into a new bowl-like shape. Young's modulus is calculated from the linear relation between stress and strain before the buckling threshold, which is in agreement with experimental results. The prediction of elasticity theory for the buckling threshold of a radially stressed plate is presented and its results are compared to the one of our atomistic simulation. The Jarzynski equality is used to estimate the difference between the free energy of the non-compressed states and the buckled states. From a calculation of the free energy we obtain the optimum radius for which the system feels the minimum boundary stress.

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: Recent experiments showed that nonuniform strain can be produced by depositing graphene over pillars. We employed atomistic calculations to study the nonuniform strain and the induced pseudomagnetic field in graphene on top of nanopillars. By decreasing the distance between the nanopillars a complex distribution for the pseudomagnetic field can be generated. Furthermore, we performed tight-binding calculations of the local density of states (LDOS) by using the relaxed graphene configuration obtained from atomistic calculations. We find that the quasiparticle LDOS are strongly modified near the pillars, both at low energies showing sublattice polarization and at high energies showing shifts of the van Hove singularity. Our study shows that changing the specific pattern of the nanopillars allows us to create a desired shape of the pseudomagnetic field profile while the LDOS maps provide an input for experimental verification by scanning tunneling microscopy.

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: 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: The adsorption of ammonia molecule onto a graphene hexagonal flake, aluminum (Al) and boron (B) doped graphene flakes (graphene quantum dots, GQDs) are investigated using density functional theory. We found that NH3 molecule is absorbed to the hollow site through the physisorption mechanism without altering the electronic properties of GQD. However, the adsorption energy of NH3 molecule onto the Al- and B-doped GQDs increases with respect GQD resulting chemisorption. The adsorption of NH3 onto the Al-doped and B-doped GQDs makes graphene locally buckled, i.e., B-doped and Al-doped GQDs are not planar. The adsorption mechanism onto a GQD is different than that of graphene. This study reveals important features of the edge passivation and doping effects of the adsorption mechanism of external molecules onto the graphene quantum dots. (C) 2013 AIP Publishing LLC.

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: The instability of a free-standing one-sided hydrogenated/fluorinated graphene nanoribbon, i.e., graphone/fluorographene, is studied using ab initio, semiempirical, and large-scale molecular dynamics simulations. Free-standing semi-infinite armchairlike hydrogenated/fluorinated graphene (AC-GH/AC-GF) and boatlike hydrogenated/fluorinated graphene (B-GH/B-GF) (nanoribbons which are periodic along the zigzag direction) are unstable and spontaneously transform into spiral structures. We find that rolled, spiral B-GH and B-GF are energetically more favorable than spiral AC-GH and AC-GF which is opposite to the double-sided flat hydrogenated/fluorinated graphene, i.e., graphane/fluorographene. We found that the packed, spiral structures exhibit an unexpected localized highest occupied molecular orbital and lowest occupied molecular orbital at the edges with increasing energy gap during rolling. These rolled hydrocarbon structures are stable beyond room temperature up to at least T = 1000 K within our simulation time of 1 ns. DOI: 10.1103/PhysRevB.87.075448

Abstract: Using atomistic simulations we investigate the thermodynamical properties of a single atomic layer of hexagonal boron nitride (h-BN). The thermal induced ripples, heat capacity, and thermal lattice expansion of large scale h-BN sheets are determined and compared to those found for graphene (GE) for temperatures up to 1000 K. By analyzing the mean-square height fluctuations < h(2)> and the height-height correlation function H(q) we found that the h-BN sheet is a less stiff material as compared to graphene. The bending rigidity of h-BN (i) is about 16% smaller than the one of GE at room temperature (300 K), and (ii) increases with temperature as in GE. The difference in stiffness between h-BN and GE results in unequal responses to external uniaxial and shear stress and different buckling transitions. In contrast to a GE sheet, the buckling transition of a h-BN sheet depends strongly on the direction of the applied compression. The molar heat capacity, thermal-expansion coefficient, and Gruneisen parameter are estimated to be 25.2 J mol(-1) K-1, 7.2 x 10(-6) K-1, and 0.89, respectively.

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: Twisted graphene layers produce a moire pattern (MP) structure with a predetermined wavelength for a given twist angle. However, predicting the membrane corrugation amplitude for any angle other than pure AB-stacked or AA-stacked graphene is impossible using first-principles density functional theory (DFT) due to the large supercell. Here, within elasticity theory, we define the MP structure as the minimum-energy configuration, thereby leaving the height amplitude as the only unknown parameter. The latter is determined from DFT calculations for AB-and AA-stacked bilayer graphene in order to eliminate all fitting parameters. Excellent agreement with scanning tunneling microscopy results across multiple substrates is reported as a function of twist angle.

Abstract: Freestanding graphene membranes were successfully functionalized with platinum nanoparticles (Pt NPs). High-resolution transmission electron microscopy revealed a homogeneous distribution of single-crystal Pt NPs that tend to exhibit a preferred orientation. Unexpectedly, the NPs were also found to be partially exposed to the vacuum with the top Pt surface raised above the graphene substrate, as deduced from atomic-scale scanning tunneling microscopy images and detailed molecular dynamics simulations. Local strain accumulation during the growth process is thought to be the origin of the NP self-organization. These findings are expected to shape future approaches in developing Pt NP catalysts for fuel cells as well as NP-functionalized graphene-based high-performance electronics.

Abstract: Knowledge of and control over the curvature of ripples in freestanding graphene are desirable for fabricating and designing flexible electronic devices, and recent progress in these pursuits has been achieved using several advanced techniques such as scanning tunnelling microscopy. The electrostatic forces induced through a bias voltage (or gate voltage) were used to manipulate the interaction of freestanding graphene with a tip (substrate). Such forces can cause large movements and sudden changes in curvature through mirror buckling. Here we explore an alternative mechanism, thermal load, to control the curvature of graphene. We demonstrate thermal mirror buckling of graphene by scanning tunnelling microscopy and large-scale molecular dynamic simulations. The negative thermal expansion coefficient of graphene is an essential ingredient in explaining the observed effects. This new control mechanism represents a fundamental advance in understanding the influence of temperature gradients on the dynamics of freestanding graphene and future applications with electro-thermal-mechanical nanodevices.