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: 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: Integration of graphene-based two-dimensional materials is essential for nanoelectronics applications. Using density-functional theory, we systematically investigate the electronic properties of vertically stacked graphene-graphitic carbon nitrides (GE/GCN). We also studied the covalently lateral stitched graphene-graphitic carbon nitrides (GE-GCN heterojunctions). The effects of perpendicular electric field on the electronic properties of six different heterostructures, i.e., (i) one layer of GE on top of a layer of CnNm with (n, m) = (3,1), (3,4), and (4,3) and (ii) three heterostructures CnNm/Cn'Nm', where (n, m) not equal (n', m') are elucidated. The most important calculated features are (i) the systems GE/C3N4, C3N/C3N4, GE-C3N, GE-C4N3, and C3N-C3N4 exhibit semiconducting characteristics having small band gaps of Delta(0)=20, 250, 100, 100, 80 meV, respectively while (ii) the systems GE/C4N3, C3N/C4N3, and C3N-C4N3 show ferromagnetic-metallic properties. In particular, we found that, in semiconducting heterostructures, the band gap increases nontrivially with increasing the absolute value of the applied perpendicular electric field. This work is useful for designing heterojunctions and heterostructures made of graphene and other two-dimensional materials such as those proposed in recent experiments [X. Liu and M. C. Hersam Sci. Adv. 5, 6444 (2019)].

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