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.

Abstract: We investigate the stability of circular monolayer graphene subjected to a radial load using non-equilibrium molecular dynamics simulations. When monolayer graphene is radially stressed, after some small circular strain (~0.4%) it buckles and bends into a new bowl-like shape. Young's modulus is calculated from the linear relation between stress and strain before the buckling threshold, which is in agreement with experimental results. The prediction of elasticity theory for the buckling threshold of a radially stressed plate is presented and its results are compared to the one of our atomistic simulation. The Jarzynski equality is used to estimate the difference between the free energy of the non-compressed states and the buckled states. From a calculation of the free energy we obtain the optimum radius for which the system feels the minimum boundary stress.

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: 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: An external electric field changes the physical properties of polar liquids due to the reorientation of their permanent dipoles. Using molecular dynamics simulations, we predict that an in-plane electric field applied parallel to the channel polarizes water molecules which are confined between two graphene layers, resulting in distinct ferroelectricity and electrical hysteresis. We found that electric fields alter the in-plane order of the hydrogen bonds: Reversing the electric field does not restore the system to the nonpolar initial state, instead a residual dipole moment remains in the system. The square-rhombic structure of 2D ice is transformed into two rhombic-rhombic structures. Our study provides insights into the ferroelectric state of water when confined in nanochannels and shows how this can be tuned by an electric field.

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: 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: 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: The dielectric properties of confined water is of fundamental interest and is still controversial. For water confined in channels with height smaller than h = 8 angstrom, we found a commensurability effect and an extraordinary decrease in the out-of-plane dielectric constant down to the limit of the dielectric constant of optical water. Spatial resolved polarization density data obtained from molecular dynamics simulations are found to be antisymmetric across the channel and are used as input in a mean-field model for the dielectric constant as a function of the height of the channel for h > 15 angstrom. Our results are in excellent agreement with a recent experiment [L. Fumagalli et al., Science 360, 1339 (2018)].

Abstract: The buckling of defected rectangular graphene nanoribbons when subjected to axial stress with supported boundary conditions is investigated using atomistic simulations. The buckling strain and mechanical stiffness of monolayer graphene decrease with the percentage of randomly distributed vacancies. The elasticity to plasticity transition in the stress-strain curve, at low percentage of vacancies, are found to be almost equal to the buckling strain thresholds and they decrease with increasing percentage of vacancies.

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