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: 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: Hyperbolic polaritons in ultra thin materials such as few layers of van derWaals heterostructures provide a unique control over light-matter interaction at the nanoscale and with various applications in flat optics. Natural hyperbolic surface plasmons have been observed on thin films of WTe2 in the light wavelength range of 16-23 mu m (similar or equal to 13-18 THz) [Nat. Commun. 11, 1158 (2020)]. Using time-dependent density functional theory, it is found that carbon doped monolayer phosphorene (beta-allotrope of carbon phosphide monolayer) exhibits natural hyperbolic plasmons at frequencies above similar or equal to 5 THz which is not observed in its parent materials, i.e. monolayer of black phosphorous and graphene. Furthermore, we found that by electrostatic doping the plasmonic frequency range can be extended to the mid-infrared. (C) 2021 Elsevier Ltd. All rights reserved.

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

Abstract: Using atomistic simulations we investigate the morphological properties of graphene deposited on top of a nanostructured substrate. Sinusoidally corrugated surfaces, steps, elongated trenches, one-dimensional and cubic barriers, spherical bubbles, Gaussian bumps, and Gaussian depressions are considered as support structures for graphene. The graphene-substrate interaction is governed by van der Waals forces and the profile of the graphene layer is determined by minimizing the energy using molecular dynamics simulations. Based on the obtained optimum configurations, we found that (i) for graphene placed over sinusoidally corrugated substrates with corrugation wavelengths longer than 2 nm, the graphene sheet follows the substrate pattern while for supported graphene it is always suspended across the peaks of the substrate, (ii) the conformation of graphene to the substrate topography is enhanced when increasing the energy parameter in the van der Waals model, (iii) the adhesion of graphene into the trenches depends on the width of the trench and on the graphene's orientation, i. e., in contrast to a small-width (3 nm) nanoribbon with armchair edges, the one with zigzag edges follows the substrate profile, (iv) atomic-scale graphene follows a Gaussian bump substrate but not the substrate with a Gaussian depression, and (v) the adhesion energy due to van der Waals interaction varies in the range [0.1-0.4] J/m(2).

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: 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: Large-scale atomistic simulations using the reactive empirical bond order force field approach is implemented to investigate thermal and mechanical properties of single-layer (SL) and multilayer (ML) molybdenum disulfide (MoS2). The amplitude of the intrinsic ripples of SL MoS2 are found to be smaller than those exhibited by graphene (GE). Furthermore, because of the van der Waals interaction between layers, the out-of-plane thermal fluctuations of ML MoS2 decreases rapidly with increasing number of layers. This trend is confirmed by the buckling transition due to uniaxial stress which occurs for a significantly larger applied tension as compared to graphene. For SL MoS2, the melting temperature is estimated to be 3700 K which occurs through dimerization followed by the formation of small molecules consisting of two to five atoms. When different types of vacancies are inserted in the SL MoS2 it results in a decrease of both the melting temperature as well as the stiffness.

Abstract: The available synthesized silicene-like structures have been only realized on metallic substrates which are very different from the standalone buckled silicene, e. g., the Dirac cone of silicene is destroyed due to lattice distortion and the interaction with the substrate. Using graphene bilayer as a scaffold, a route is proposed to synthesize silicene with electronic properties decoupled from the substrate. The buckled hexagonal arrangement of silicene between the graphene layers is found to be very similar to the theoretically predicted standalone buckled silicene which is only very weakly van der Waals coupled to the graphene layers with a graphite-like interlayer distance of 3.42 angstrom and without any lattice distortion. We found that these stacked layers are stable well above room temperature. (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.

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.