Abstract: A multiscale modeling and simulation approach, including first-principles calculations, ab initio molecular dynamics simulations, and a tight binding approach, is employed to study band flattening of the electronic band structure of oxidized monolayer graphene. The width offlat bands can be tuned by strain, the external electric field, and the density of functional groups and their distribution. A transition to a conducting state is found for monolayer graphene with impurities when it is subjected to an electric field of similar to 1.0 V/angstrom. Several parallel impurity-induced flat bands appear in the low-energy spectrum of monolayer graphene when the number of epoxy groups is changed. The width of the flat band decreases with an increase in tensile strain but is independent of the electric field strength. Here an alternative and easy route for obtaining band flattening in thermodynamically stable functionalized monolayer graphene is introduced. Our work discloses a new avenue for research on band flattening in monolayer graphene.

Abstract: The self rotation of a graphene flake over graphite is controlled by the size, initial misalignment and temperature. Using both ab initio calculations and molecular dynamics simulations, we investigate annealing effects on the self rotation of a graphene flake on a graphene substrate. The energy barriers for rotation and drift of a graphene flake over graphene is found to be smaller than 25 meV/atom which is comparable to thermal energy. We found that small flakes (of about similar to 4 nm) are more sensitive to temperature and initial misorientation angles than larger one (beyond 10 nm). The initial stacking configuration of the flake is found to be important for its dynamics and time evolution of misalignment. Large flakes, which are initially in the AA-or AB-stacking state with small misorientation angle, rotate and end up in the AB-stacking configuration. However small flakes can they stay in an incommensurate state specially when the initial misorientation angle is larger than 2 degrees. Our results are in agreement with recent experiments.

Abstract: Intrinsic ripples in freestanding graphene have been exceedingly difficult to study. Individual ripple geometry was recently imaged using scanning tunnelling microscopy, but these measurements are limited to static configurations. Thermally-activated flexural phonon modes should generate dynamic changes in curvature. Here we show how to track the vertical movement of a one-square-angstrom region of freestanding graphene using scanning tunnelling microscopy, thereby allowing measurement of the out-of-plane time trajectory and fluctuations over long time periods. We also present a model from elasticity theory to explain the very-low-frequency oscillations. Unexpectedly, we sometimes detect a sudden colossal jump, which we interpret as due to mirror buckling. This innovative technique provides a much needed atomic-scale probe for the time-dependent behaviours of intrinsic ripples. The discovery of this novel progenitor represents a fundamental advance in the use of scanning tunnelling microscopy, which together with the application of a thermal load provides a low-frequency nano-resonator.

Abstract: We show that a carbon nanotube decorated with different types of charged metallic nanoparticles exhibits unusual two-dimensional vibrations when actuated by applied electric field. Such vibrations and diverse possible trajectories are not only fundamentally important but also have minimum two characteristic frequencies that can be directly linked back to the properties of the constituents in the considered nanoresonator. Namely, those frequencies and the maximal deflection during vibrations are very distinctively dependent on the geometry of the nanotube, the shape, element, mass and charge of the nanoparticle, and are vastly tunable by the applied electric field, revealing the unique sensing ability of devices made of molecular filaments and metallic nanoparticles.

Abstract: Van der Waals (vdW) interaction between two-dimensional crystals (2D) can trap substances in high pressurized (of order 1 GPa) on nanobubbles. Increasing the adhesion between the 2D crystals further enhances the pressure and can lead to a phase transition of the trapped material. We found that the shape of the nanobubble can depend critically on the properties of the trapped substance. In the absence of any residual strain in the top 2D crystal, flat nanobubbles can be formed by trapped long hydrocarbons (that is, hexadecane). For large nanobubbles with radius 130 nm, our atomic force microscopy measurements show nanobubbles filled with hydrocarbons (water) have a cylindrical symmetry (asymmetric) shape which is in good agreement with our molecular dynamics simulations. This study provides insights into the effects of the specific material and the vdW pressure on the microscopic details of graphene bubbles.

Abstract: Using electrochemical methods a profound enhancement of the capacitance of electric double layer capacitor electrodes was reported when water molecules are strongly confined into the two-dimensional slits of titanium carbide MXene nanosheets [A. Sugahara et al., Nat. Commun., 2019, 10, 850]. We study the effects of hydration on the dielectric properties of nanoconfined water and supercapacitance properties of the cation intercalated MXene. A model for the electric double layer capacitor is constructed where water molecules are strongly confined in two-dimensional slits of MXene. We report an abnormal dielectric constant and polarization of nano-confined water between MXene layers. We found that by decreasing the ionic radius of the intercalated cations and in a critical hydration shell radius the capacitance of the system increases significantly (similar or equal to 200 F g(-1)) which can be interpreted as a negative permittivity. This study builds a bridge between the fundamental understanding of the dielectric properties of nanoconfined water and the capability of using MXene films for supercapacitor technology, and in doing so provides a solid theoretical support for recent experiments.

Abstract: Molecular dynamics simulations are used to investigate the effect of an AFM tip when indenting graphene nano bubbles filled by a noble gas (i.e. He, Ne and Ar) up to the breaking point. The failure points resemble those of viral shells as described by the Foppl-von Karman (FvK) dimensionless number defined in the context of elasticity theory of thin shells. At room temperature, He gas inside the bubbles is found to be in the liquid state while Ne and Ar atoms are in the solid state although the pressure inside the nano bubble is below the melting pressure of the bulk. The trapped gases are under higher hydrostatic pressure at low temperatures than at room temperature.

Abstract: When confined between graphene layers, water behaves differently from the bulk and exhibits unusual properties such as fast water flow and ordering into a crystal. The hydrogen-bonded network is affected by the limited space and by the characteristics of the confining walls. The presence of an extraordinary number of hydronium and hydroxide ions in narrow channels has the following effects: (i) they affect water permeation through the channel, (ii) they may interact with functional groups on the graphene oxide surface and on the edges, and (iii) they change the thermochemistry of water, which are fundamentally important to understand, especially when confined water is subjected to an external electric field. Here we study the physical properties of water when confined between two graphene sheets and containing hydronium and hydroxide. We found that: (i) there is a disruption in the solvation structure of the ions, which is also affected by the layered structure of confined water, (ii) hydronium and hydroxide occupy specific regions inside the nanochannel, with a prevalence of hydronium (hydroxide) ions at the edges (interior), and (iii) ions recombine more slowly in confined systems than in bulk water, with the recombination process depending on the channel height and commensurability between the size of the molecules and the nanochannel height – a decay of 20% (40%) in the number of ions in 8 ps is observed for a channel height of h = 7 angstrom (bulk water). Our work reveals distinctive properties of water confined in a nanocapillary in the presence of additional hydronium and hydroxide ions.

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: 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: 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: 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: 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: 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