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: 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: Controlled transport of water molecules through membranes and capillaries is important in areas as diverse as water purification and healthcare technologies(1-7). Previous attempts to control water permeation through membranes (mainly polymeric ones) have concentrated on modulating the structure of the membrane and the physicochemical properties of its surface by varying the pH, temperature or ionic strength(3,8). Electrical control over water transport is an attractive alternative; however, theory and simulations(9-14) have often yielded conflicting results, from freezing of water molecules to melting of ice(14-16) under an applied electric field. Here we report electrically controlled water permeation through micrometre-thick graphene oxide membranes(17-21). Such membranes have previously been shown to exhibit ultrafast permeation of water(17,22) and molecular sieving properties(18,21), with the potential for industrial-scale production. To achieve electrical control over water permeation, we create conductive filaments in the graphene oxide membranes via controllable electrical breakdown. The electric field that concentrates around these current-carrying filaments ionizes water molecules inside graphene capillaries within the graphene oxide membranes, which impedes water transport. We thus demonstrate precise control of water permeation, from ultrafast permeation to complete blocking. Our work opens up an avenue for developing smart membrane technologies for artificial biological systems, tissue engineering and filtration.

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: 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: 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: 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: 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: 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: 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: Molecular dynamics simulations are used to study the effects of different cations on the permeation of charged polymers through flat capillaries with heights below 2 nm. Interestingly, we found that, despite being monovalent, Li+ , Na+ , and K+ cations have different effects on polymer permeation, which consequently affects their transmission speed throughout those capillaries. We attribute this phenomenon to the interplay of the cations' hydration free energies and the hydrodynamic drag in front of the polymer when it enters the capillary. Different alkali cations exhibit different surface versus bulk preferences in small clusters of water under the influence of an external electric field. This paper presents a tool to control the speed of charged polymers in confined spaces using cations.

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