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: Using first-principle calculations, we systematically study strain tuning of the electronic properties of two- dimensional graphitic carbon nitride nanosheets with empirical formula CnNm. We found the following: (i) the ferromagnetic ordered state in the metal-free systems (n, m) = (4,3), (10,9), and (14,12) remains stable in the presence of strain of about 6%. However, the system (9,7) loses its ferromagnetic ordering when increasing strain. This is due to the presence of topological defects in the (9,7) system, which eliminates the asymmetry between spin up and spin down of the p(z) orbitals when strain is applied. (ii) By applying uniaxial strain, a band gap opens in systems which are initially gapless. (iii) In semiconducting systems which have an initial gap of about 1 eV, the band gap is closed with applying uniaxial strain.

Abstract: We report the formation of nanobubbles on graphene with a radius of the order of 1 nm, using ultralow energy implantation of noble gas ions (He, Ne, Ar) into graphene grown on a Pt(111) surface. We show that the universal scaling of the aspect ratio, which has previously been established for larger bubbles, breaks down when the bubble radius approaches 1 nm, resulting in much larger aspect ratios. Moreover, we observe that the bubble stability and aspect ratio depend on the substrate onto which the graphene is grown (bubbles are stable for Pt but not for Cu) and trapped element. We interpret these dependencies in terms of the atomic compressibility of the noble gas as well as of the adhesion energies between graphene, the substrate, and trapped atoms.

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: 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 melting of fluorographene is very unusual and depends strongly on the degree of fluorination. For temperatures below 1000 K, fully fluorinated graphene (FFG) is thermomechanically more stable than graphene but at T-m approximate to 2800 K FFG transits to random coils which is almost 2 times lower than the melting temperature of graphene, i.e., 5300 K. For fluorinated graphene up to 30% ripples causes detachment of individual F-atoms around 2000 K, while for 40%-60% fluorination large defects are formed beyond 1500 K and beyond 60% of fluorination F-atoms remain bonded to graphene until melting. The results agree with recent experiments on the dependence of the reversibility of the fluorination process on the percentage of fluorination.

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.

Abstract: Moire patterns in the pseudo-magnetic field and in the strain profile of graphene (GE) when put on top of a hexagonal lattice substrate are predicted from elasticity theory. The van der Waals interaction between GE and the substrate induces out-of-plane deformations in graphene which results in a strain field, and consequently in a pseudo-magnetic field. When the misorientation angle is about 0.5 degrees, a three-fold symmetric strain field is realized that results in a pseudo-magnetic field very similar to the one proposed by F. Guinea, M. I. Katsnelson, and A. K. Geim [Nature Phys. 6, 30 (2010)]. Our results show that the periodicity and length of the pseudo-magnetic field can be tuned in GE by changing the misorientation angle and substrate adhesion parameters and a considerable energy gap (23 meV) can be obtained due to out-of-plane deformation of graphene which is in the range of recent experimental measurements (20-30 meV). (C) 2014 AIP Publishing LLC.

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: 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 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: 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: Water inside a nanocapillary exhibits unconventional structural and dynamical behavior due to its ordered structure. The confining walls, density, and lateral pressures control profoundly the microscopic structure of trapped water. Here we study the electrostriction of confined water subjected to pressures of the order of GPa for two different setups: (i) a graphene nanochannel containing a constant number of water molecules independent of the height of the channel, (ii) an open nanochannel where water molecules can be exchanged with those in a reservoir. For the former case, a square-rhombic structure of confined water is formed when the height of the channel is d = 6.5 angstrom having a density of rho = 1.42 g cm(-3). By increasing the height of the channel, a transition from a flat to a buckled state occurs, whereas the density rapidly decreases and reaches the bulk density for d congruent to 8.5 angstrom. When a perpendicular electric field is applied, the water structure and the lateral pressure change. For strong electric fields (similar to 1 V/angstrom), the square-rhombic structure is destroyed. For an open setup, a solid phase of confined water consisting of an imperfect square-rhombic structure is formed. By applying a perpendicular field, the density and phase of confined water change. However, the density and pressure inside the channel decrease as compared to the first setup. Our study is closely related to recent experiments on confined water, and it reveals the sensitivity of the microscopic structure of confined water to the size of the channel, the external electric field, and the experimental setup.

Abstract: Piezo- and flexoelectricity are manifestations of electromechanical coupling in solids with potential applications in nanoscale materials. Naumov etal. [Phys. Rev. Lett. 102, 217601 (2009)] have shown by first principles calculations that a monolayer BN sheet becomes macroscopically polarized in-plane when in a corrugated state. Here, we investigate the interplay of layer corrugation and in-plane polarization by atomistic lattice dynamics. We treat the coupling between static flexural modes and in-plane atomic ion displacements as an anharmonic effect, similar to the membrane effect that is at the origin of negative thermal expansion in layered crystals. We have derived analytical expressions for the corrugation-induced static in-plane strains and the optical displacements with the resulting polarization response functions. Beyond h-BN, the theory applies to transition metal dichalcogenides and dioxides. Numerical calculations show that the effects are considerably stronger for 2D h-BN than for 2H-MoS2.

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.