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 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: 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: 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: Large-scale atomistic simulations using the reactive force field approach are implemented to investigate the thermomechanical properties of fluorinated graphene (FG). A set of parameters for the reactive force field potential optimized to reproduce key quantum mechanical properties of relevant carbon-fluorine cluster systems are presented. Molecular dynamics simulations are used to investigate the thermal rippling behavior of FG and its mechanical properties and compare them with graphene, graphane and a sheet of boron nitride. The mean square value of the height fluctuations < h(2)> and the height-height correlation function H(q) for different system sizes and temperatures show that FG is an unrippled system in contrast to the thermal rippling behavior of graphene. The effective Young's modulus of a flake of fluorinated graphene is obtained to be 273 N/m and 250 N/m for a flake of FG under uniaxial strain along armchair and zigzag directions, respectively. DOI: 10.1103/PhysRevB.87.104114

Abstract: The electronic properties of boron nitride nanoribbons (BNNRs) doped with a line of carbon atoms are investigated using density functional calculations. By replacing a line of alternating B and N atoms with carbons, three different configurations are possible depending on the type of the atoms which bond to the carbons. We found very different electronic properties for these configurations: (i) the NCB arrangement is strongly polarized with a large dipole moment having an unexpected direction, (ii) the BCB and NCN arrangements are nonpolar with zero dipole moment, (iii) the doping by a carbon line reduces the band gap regardless of the local arrangement of the borons and the nitrogens around the carbon line, and (iv) the polarization and energy gap of the carbon-doped BNNRs can be tuned by an electric field applied parallel to the carbon line. Similar effects were found when either an armchair or zigzag line of carbon was introduced.

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: 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: Atomistic simulations are used to study the bending of rectangular graphene nanoribbons subjected to axial stress both for free boundary and supported boundary conditions. The shapes of the deformations of the buckled graphene nanoribbons, for small values of the stress, are sine waves where the number of nodal lines depend on the longitudinal size of the system and the applied boundary condition. The buckling strain for the supported boundary condition is found to be independent of the longitudinal size and estimated to be 0.86%. From a calculation of the free energy at finite temperature we find that the equilibrium projected two-dimensional area of the graphene nanoribbon is less than the area of a flat sheet. At the optimum length the boundary strain for the supported boundary condition is 0.48%.

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 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: 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: The ability of different materials to display self-limiting growth has recently attracted an enormous amount of attention because of the importance of nanoscale materials in applications for catalysis, energy conversion, (opto)-electronics, and so forth. Here, we show that the electrochemical deposition of palladium (Pd) between graphene oxide (GO) sheets result in the self-limiting growth of 5-nm-thick Pd nanosheets. The self-limiting growth is found to be a consequence of the strong interaction of Pd with the confining GO sheets, which results in the bulk growth of Pd being energetically unfavorable for larger thicknesses. Furthermore, we have successfully carried out liquid exfoliation of the resulting Pd-GO laminates to isolate Pd nanosheets and have demonstrated their high efficiency in continuous flow catalysis and electrocatalysis.

Abstract: Atoms start behaving as waves rather than classical particles if confined in spaces commensurate with their de Broglie wavelength. At room temperature this length is only about one angstrom even for the lightest atom, hydrogen. This restricts quantum-confinement phenomena for atomic species to the realm of very low temperatures(1-5). Here, we show that van der Waals gaps between atomic planes of layered crystals provide angstrom-size channels that make quantum confinement of protons apparent even at room temperature. Our transport measurements show that thermal protons experience a notably higher barrier than deuterons when entering van der Waals gaps in hexagonal boron nitride and molybdenum disulfide. This is attributed to the difference in the de Broglie wavelengths of the isotopes. Once inside the crystals, transport of both isotopes can be described by classical diffusion, albeit with unexpectedly fast rates comparable to that of protons in water. The demonstrated angstrom-size channels can be exploited for further studies of atomistic quantum confinement and, if the technology can be scaled up, for sieving hydrogen isotopes.

Abstract: The recent advantages of the fabrication of artificial nanochannels enabled new research on the molecular transport, permeance, and selectivity of various gases and molecules. However, the physisorption/chemisorption of the unwanted molecules (usually hydrocarbons) inside nanochannels results in the alteration of the functionality of the nanochannels. We investigated contamination due to hydrocarbon molecules, nanochannels made of graphene, hexagonal boron nitride, BC2N, and molybdenum disulfide using molecular dynamics simulations. We found that for a certain size of nanochannel (i.e., h = 0.7 nm), as a result of the anomalous hydrophilic nature of nanochannels made of graphene, the hydrocarbons are fully adsorbed in the nanochannel, giving rise to full uptake. An increasing temperature plays an important role in unclogging, while pressure does not have a significant role. The results of our pioneering work contribute to a better understanding and highlight the important factors in alleviating the contamination and unclogging of nanochannels, which are in good agreement with the results of recent experiments.