Abstract: The permeation of water between neighboring graphene oxide (GO) flakes, i.e., 2D nanochannels, are investigated using a simple model for the GO membrane. We simulate the hydrophilic behavior of nanocapillaries and study the effect of surface charge on the dynamical properties of water flow and the influence of Na+ and Cl- ions on water permeation. Our approach is based on extensive equilibrium molecular dynamics simulations to obtain a better understanding of water permeation through charged nanochannels in the presence of ions. We found significant change in the slippage dynamics of confined water such as a profound increase in viscosity/slip length with increasing charges over the surface. The slip length decreases one order of magnitude (i.e., 1/30) with increasing density of surface charge, while it increases by a factor of 2 with ion concentration. We found that commensurability induced by nanoconfinement plays an important role on the intrinsic dynamical properties of water.

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: We show that a carbon nanotube decorated with different types of charged metallic nanoparticles exhibits unusual two-dimensional vibrations when actuated by applied electric field. Such vibrations and diverse possible trajectories are not only fundamentally important but also have minimum two characteristic frequencies that can be directly linked back to the properties of the constituents in the considered nanoresonator. Namely, those frequencies and the maximal deflection during vibrations are very distinctively dependent on the geometry of the nanotube, the shape, element, mass and charge of the nanoparticle, and are vastly tunable by the applied electric field, revealing the unique sensing ability of devices made of molecular filaments and metallic nanoparticles.

Abstract: Van der Waals (vdW) interaction between two-dimensional crystals (2D) can trap substances in high pressurized (of order 1 GPa) on nanobubbles. Increasing the adhesion between the 2D crystals further enhances the pressure and can lead to a phase transition of the trapped material. We found that the shape of the nanobubble can depend critically on the properties of the trapped substance. In the absence of any residual strain in the top 2D crystal, flat nanobubbles can be formed by trapped long hydrocarbons (that is, hexadecane). For large nanobubbles with radius 130 nm, our atomic force microscopy measurements show nanobubbles filled with hydrocarbons (water) have a cylindrical symmetry (asymmetric) shape which is in good agreement with our molecular dynamics simulations. This study provides insights into the effects of the specific material and the vdW pressure on the microscopic details of graphene bubbles.

Abstract: The report on square ice sandwiched between two graphene layers by Algara-Siller et al. [Nature (London) 519, 443 (2015)] has generated a large interest in this system. By applying high lateral pressure on nanoconfined water, we found that monolayer ice is transformed to bilayer ice when the two graphene layers are separated by H = 6,7 angstrom. It was also found that three layers of a denser phase of ice with smaller lattice constant are formed if we start from bilayer ice and apply a lateral pressure of about 0.7 GPa with H = 8,9 angstrom. The lattice constant (2.5-2.6 angstrom) in both transitions is found to be smaller than those typical for the known phases of ice and water, i.e., 2.8 angstrom. We validate these results using ab initio calculations and find good agreement between ab initio O-O distance and those obtained from classical molecular dynamics simulations. The reversibility of the mentioned transitions is confirmed by decompressing the systems.

Abstract: We report subnanometer, high-bandwidth measurements of the out-of-plane (vertical) motion of atoms in freestanding graphene using scanning tunneling microscopy. By tracking the vertical position over a long time period, a 1000-fold increase in the ability to measure space-time dynamics of atomically thin membranes is achieved over the current state-of-the-art imaging technologies. We observe that the vertical motion of a graphene membrane exhibits rare long-scale excursions characterized by both anomalous mean-squared displacements and Cauchy-Lorentz power law jump distributions.

Abstract: The rate of water flow through hydrophobic nanocapillaries is greatly enhanced as compared to that expected from macroscopic hydrodynamics. This phenomenon is usually described in terms of a relatively large slip length, which is in turn defined by such microscopic properties as the friction between water and capillary surfaces and the viscosity of water. We show that the viscosity of water and, therefore, its flow rate are profoundly affected by the layered structure of confined water if the capillary size becomes less than 2 nm. To this end, we study the structure and dynamics of water confined between two parallel graphene layers using equilibrium molecular dynamics simulations. We find that the shear viscosity is not only greatly enhanced for subnanometer capillaries, but also exhibits large oscillations that originate from commensurability between the capillary size and the size of water molecules. Such oscillating behavior of viscosity and, consequently, the slip length should be taken into account in designing and studying graphene-based and similar membranes for desalination and filtration.

Abstract: Water confined between two graphene layers with a separation of a few A forms a layered two-dimensional ice structure. Using large scale molecular dynamics simulations with the adoptable ReaxFF interatomic potential we found that flat monolayer ice with a rhombic-square structure nucleates between the graphene layers which is nonpolar and nonferroelectric. We provide different energetic considerations and H-bonding results that explain the interlayer and intralayer properties of two-dimensional ice. The controversial AA stacking found experimentally [Algara-Siller et al., Nature (London) 519, 443 (2015)] is consistent with our minimum-energy crystal structure of bilayer ice. Furthermore, we predict that an odd number of layers of ice has the same lattice structure as monolayer ice, while an even number of ice layers exhibits the square ice AA stacking of bilayer ice.

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: 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: 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: 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: The electronic properties of graphene nano-flakes (GNFs) with different edge passivation are investigated by using density functional theory. Passivation with F and H atoms is considered: C-Nc X-Nx (X = F or H). We studied GNFs with 10 < N-c < 56 and limit ourselves to the lowest energy configurations. We found that: (i) the energy difference Delta between the highest occupied molecular orbital and the lowest unoccupied molecular orbital decreases with N-c, (ii) topological defects (pentagon and heptagon) break the symmetry of the GNFs and enhance the electric polarization, (iii) the mutual interaction of bilayer GNFs can be understood by dipole-dipole interaction which were found sensitive to the relative orientation of the GNFs, (iv) the permanent dipoles depend on the edge terminated atom, while the energy gap is independent of it, and (v) the presence of heptagon and pentagon defects in the GNFs results in the largest difference between the energy of the spin-up and spin-down electrons which is larger for the H-passivated GNFs as compared to F-passivated GNFs. Our study shows clearly the effect of geometry, size, termination, and bilayer on the electronic properties of small GNFs. This study reveals important features of graphene nano-flakes which can be detected using Raman spectroscopy. (C) 2014 AIP Publishing LLC.

Abstract: The spatial dependence of the van der Waals (vdW) energy between graphene and hexagonal boron-nitride (h-BN) is investigated using atomistic simulations. The van der Waals energy between graphene and h-BN shows a hexagonal superlattice structure identical to the observed Moire pattern in the local density of states, which depends on the lattice mismatch and misorientation angle between graphene and h-BN. Our results provide atomistic features of the weak van der Waals interaction between graphene and BN which are in agreement with experiment and provide an analytical expression for the size of the spatial variation of the weak van der Waals interaction. We also found that the A-B-lattice symmetry of graphene is broken along the armchair direction. (C) 2014 AIP Publishing LLC.