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 influence of triaxial in-plane strain on the electronic properties of a hexagonal boron-nitride sheet is investigated using density functional theory. Different from graphene, the triaxial strain localizes the molecular orbitals of the boron-nitride flake in its center depending on the direction of the applied strain. The proposed technique for localizing the molecular orbitals that are close to the Fermi level in the center of boron nitride flakes can be used to actualize engineered nanosensors, for instance, to selectively detect gas molecules. We show that the central part of the strained flake adsorbs polar molecules more strongly as compared with an unstrained sheet.

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: 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: 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: 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: 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: 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: Understanding the appearance of local electrical connectivity in liquid filled layered graphene oxide subjected to an external electric field is important to design electrically controlled smart permeable devices and also to gain insight into the physics behind electrical effects on confined water permeation. Motivated by recent experiments [K. G. Zhou et al. Nature (London) 559, 236 (2018)], we introduce a new model with random percolating paths for electrical connectivity in micron thick water filled layered graphene oxide, which mimics parallel resistors connected across the top and bottom electrodes. We find that a strong nonuniform radial electric field of the order similar to 10-50 mV/nm can be induced between layers depending on the current flow through the formed conducting paths. The maxima of the induced fields are not necessarily close to the electrodes and may be localized in the middle region of the layered material. The emergence of electrical connectivity and the associated electrical effects have a strong influence on the surrounding fluid in terms of ionization and wetting which subsequently determines the permeation properties.

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: 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: The report on square ice sandwiched between two graphene layers by Algara-Siller et al. [Nature (London) 519, 443 (2015)] has generated a large interest in this system. By applying high lateral pressure on nanoconfined water, we found that monolayer ice is transformed to bilayer ice when the two graphene layers are separated by H = 6,7 angstrom. It was also found that three layers of a denser phase of ice with smaller lattice constant are formed if we start from bilayer ice and apply a lateral pressure of about 0.7 GPa with H = 8,9 angstrom. The lattice constant (2.5-2.6 angstrom) in both transitions is found to be smaller than those typical for the known phases of ice and water, i.e., 2.8 angstrom. We validate these results using ab initio calculations and find good agreement between ab initio O-O distance and those obtained from classical molecular dynamics simulations. The reversibility of the mentioned transitions is confirmed by decompressing the systems.

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: The optoelectronic properties of confined water form one of the most active research areas in the past few years. Here we present the multiscale methodology to discern the out-of-plane electronic and dipolar dielectric constants (epsilon(el)(perpendicular to) and epsilon(diP)(perpendicular to)) of strongly confined water. We reveal that epsilon(perpendicular to el) and epsilon(diP)(perpendicular to) become comparable for water confined in angstrom-scale channels (with a height of less than 15 angstrom) within graphene (GE) and hexagonal boron nitride (hBN) bilayers. Channel height (h) associated with a minimum in both epsilon(e)(l)(perpendicular to) and epsilon(dip)(perpendicular to) is linked to the formation of the ordered structure of ice for h approximate to (7 -7.5) angstrom. The recently measured total dielectric constant epsilon(T)(perpendicular to) of nanoconfined water [L. Fumagalli et al., Science 360, 1339 (2018)] is corroborated by our results. Furthermore, we evaluate the contribution from the encapsulating membranes to the dielectric properties, as a function of the interlayer spacing, i.e., the height of the confining channel for water. Finally, we conduct analysis of the optical properties of both confined water and GE membranes, and show that the electron energy loss function of confined water strongly differs from that of bulk water.