Manipulation of graphene's dynamic ripples by local harmonic out-of-plane excitation

Corrugated graphene layers for sea water desalination using capacitive deionization

The effect of the electric field and surface morphology of corrugated graphene (GE) layers on their capacitive deionization process is studied using molecular dynamics simulations. Deionization performances are evaluated in terms of water flow rate and ion adsorption and explained by analysing the water density distribution, radial distribution function and distribution of the ions inside the GE layers. The simulation results reveal that corrugation of GE layers reduces the water flow rate but largely enhances ion adsorption in comparison to the flat GE layers. Such enhancement is mainly due to the adsorption of ions on the GE layers due to the anchoring effect in the regions with wide interlayer distances. Moreover, it reveals that the entrance configuration of the GE layers also has a significant effect on the performance of deionization. Overall, the results from this study will be helpful in designing effective electrode configurations for capacitive deionization.

Three-dimensional surface topography of graphene by divergent beam electron diffraction

There are only a handful of scanning techniques that can provide surface topography at nanometre resolution. At the same time, there are no methods that are capable of non-invasive imaging of the three-dimensional surface topography of a thin free-standing crystalline material. Here we propose a new technique—the divergent beam electron diffraction (DBED) and show that it can directly image the inhomogeneity in the atomic positions in a crystal. Such inhomogeneities are directly transformed into the intensity contrast in the first-order diffraction spots of DBED patterns and the intensity contrast linearly depends on the wavelength of the employed probing electrons. Three-dimensional displacement of atoms as small as 1 angstrom can be detected when imaged with low-energy electrons (50–250 eV). The main advantage of DBED is that it allows visualization of the three-dimensional surface topography and strain distribution at the nanometre scale in non-scanning mode, from a single shot diffraction experiment.

Graphene, and other 2D materials, do not exist as strictly planar sheets but instead have topographic ripples on the sub-nanometre scale. Here, Latychevskaia et al. present a method to non-invasively image ripples in 2D materials with a single-shot, wide-area, electron diffraction measurement.

Rippling ultrafast dynamics of suspended 2D monolayers, graphene

Here, using ultrafast electron crystallography (UEC), we report the observation of rippling dynamics in suspended monolayer graphene, the prototypical and most-studied 2D material. The high scattering cross-section for electron/matter interaction, the atomic-scale spatial resolution, and the ultrafast temporal resolution of UEC represent the key elements that make this technique a unique tool for the dynamic investigation of 2D materials, and nanostructures in general. We find that, at early time after the ultrafast optical excitation, graphene undergoes a lattice expansion on a time scale of 5 ps, which is due to the excitation of short-wavelength in-plane acoustic phonon modes that stretch the graphene plane. On a longer time scale, a slower thermal contraction with a time constant of 50 ps is observed and associated with the excitation of out-of-plane phonon modes, which drive the lattice toward thermal equilibrium with the well-known negative thermal expansion coefficient of graphene. From our results and first-principles lattice dynamics and out-of-equilibrium relaxation calculations, we quantitatively elucidate the deformation dynamics of the graphene unit cell.

Green's function modeling of response of two-dimensional materials to point probes for scanning probe microscopy

A Green's function (GF) method is developed for interpreting scanning probe microscopy (SPM) measurements on new two-dimensional (2D) materials. GFs for the Laplace/Poisson equations are calculated by using a virtual source method for two separate cases of a finite material containing a rectangular defect and a hexagonal defect. The prescribed boundary values are reproduced almost exactly by the calculated GFs. It is suggested that the GF is not just a mathematical artefact but a basic physical characteristic of material systems, which can be measured directly by SPM for 2D solids. This should make SPM an even more powerful technique for characterization of 2D materials.

Nucleobase-functionalized graphene nanoribbons for accurate high-speed DNA sequencing

We propose a water-immersed nucleobase-functionalized suspended graphene nanoribbon as an intrinsically selective device for nucleotide detection. The proposed sensing method combines Watson-Crick selective base pairing with graphene's capacity for converting anisotropic lattice strain to changes in an electrical current at the nanoscale. Using detailed atomistic molecular dynamics (MD) simulations, we study sensor operation at ambient conditions. We combine simulated data with theoretical arguments to estimate the levels of measurable electrical signal variation in response to strains and determine that the proposed sensing mechanism shows significant promise for realistic DNA sensing devices without the need for advanced data processing, or highly restrictive operational conditions.

Anisotropic growth of buckling-driven wrinkles in graphene monolayer

We theoretically and numerically investigate the growth of buckling-driven wrinkles in graphene monolayers. It is found that the growth of buckling-driven wrinkles in a graphene monolayer is remarkably chirality- and size-dependent. In small sizes, the flexural response of a graphene sheet cannot be accurately described by the classical Euler regime, and the non-continuum effect leads to zigzag-along-preferred buckling. With the increase of size, the width/length ratio α of the compressed region plays an important role in the growth of buckling-driven wrinkles. When α < 0.5, the oblique buckling happens in armchair-along compression; when 0.5 < α < 1.0, the effect of edge warp leads to zigzag-along-preferred buckling. When 1.0 < α < 3.0, the potential energy density difference due to chiral bending stiffness leads to armchair-along-preferred buckling. When α > 3.0, the non-continuum effect and chiral bending stiffness can both be neglected, and the buckling in a graphene monolayer is isotropic. The chirality-along-preferred transition of compressed buckling in a graphene monolayer leads to an improved fundamental understanding of the dynamics mechanism of graphene-based nanodevices, especially for the nanodevices with high frequency response.

Shear displacement controlled periodic wrinkles in hexagonal boron nitride sheet

The initiation and geometry pattern of wrinkles in a single-layer hexagonal boron nitride sheet induced by in-plane shear displacement are studied. The periodic wrinkles in the central region are parallel to each other with an angle of approximate 50° to the fix edges, and the longitudinal shape of the wrinkle matches the sinusoidal mode shape well. The wrinkle wavelength decreases with an increase in shear loading, while the amplitude is found to initially increase and then become stable. The dependence of the wrinkle geometry on chirality and shear direction is further elucidated. This study theoretically provides a powerful way to produce uniform wrinkles in two-dimensional membranes and to tune their properties in devices.

Effects of thermal rippling on the frictional properties of free-standing graphene

A unique case of dynamically corrugated surfaces is presented for the first time, attributed to the effect of significant thermally excited flexural waves in atomically thin layers on kinetic friction.

Molecular dynamics simulation of thermal ripples in graphene with bond-order-informed harmonic constraints

We describe the results of atomistic molecular dynamics simulations of thermal rippling in graphene with the use of a generic harmonic constraint model. The distance and angular constraint constants are calculated directly from the second-generation bond-order interatomic potential that describes carbon binding in graphene. We quantify the thermal rippling process in detail by calculating the overall rippling averages, the normal-normal correlation distributions and the height distributions. In addition, we consider the effect of a dihedral angular constraint, as well as the effect of sample size on the simulated rippling averages. The dynamic corrugation morphologies of simulated graphene samples obtained with the harmonic constraint model at various temperatures are, overall, consistent with those obtained with the bond-order potential and are in qualitative accord with previously reported findings. Given the wide availability of the harmonic constraint model in various molecular mechanics implementations, along with its high computational efficiency, our results indicate a possible use for the presented model in multicomponent dynamic simulations, including atomically thin layers.