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 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: The influence of decoration with impurities and the composition dependent band gap in 2D materials has been the subject of debate for a long time. Here, by using Density Functional Theory (DFT) calculations, we systematically disclose physical properties of chlorinated phosphorene having the stoichiometry of PmCln. By analyzing the adsorption energy, charge density, migration energy barrier, structural, vibrational, and electronic properties of chlorinated phosphorene, we found that (I) the Cl-P bonds are strong with binding energy Eb =-1.61 eV, decreases with increasing n. (II) Cl atoms on phosphorene have anionic feature, (III) the migration path of Cl on phosphorene is anisotropic with an energy barrier of 0.38 eV, (IV) the phonon band dispersion reveal that chlorinated phosphorenes are stable when r <= 0.25 where r = m/n, (V) chlorinated phosphorenes is found to be a photonic crystal in the frequency range of 280 cm-1 to 325 cm-1, (VI) electronic band structure of chlorinated phosphorenes exhibits quasi-flat bands emerging around the Fermi level with widths in the range of 22 meV to 580 meV, and (VII) Cl adsorption causes a semiconducting to metallic/semi-metallic transition which makes it suitable for application as an electroactive material. To elucidate this application, we investigated the change in binding energy (Eb), specific capacity, and open-circuit voltage as a function of the density of adsorbed Cl. The theoretical storage capacity of the chlorinated phosphorene is found to be 168.19 mA h g-1with a large average voltage (similar to 2.08 V) which is ideal number as a cathode in chloride-ion batteries.

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: The influence of triaxial in-plane strain on the electronic properties of a hexagonal boron-nitride sheet is investigated using density functional theory. Different from graphene, the triaxial strain localizes the molecular orbitals of the boron-nitride flake in its center depending on the direction of the applied strain. The proposed technique for localizing the molecular orbitals that are close to the Fermi level in the center of boron nitride flakes can be used to actualize engineered nanosensors, for instance, to selectively detect gas molecules. We show that the central part of the strained flake adsorbs polar molecules more strongly as compared with an unstrained sheet.

Abstract: Dielectric properties of nano-confined water are important in several areas of science, i.e., it is relevant in the dielectric double layer that exists in practically all heterogeneous fluid-based systems. Molecular dynamics simulations are used to predict the in-plane dielectric properties of confined water in planar channels of width ranging from sub-nanometer to bulk. Because of suppressed rotational degrees of freedom near the confining walls, the dipole of the water molecules tends to be aligned parallel to the walls, which results in a strongly enhanced in-plane dielectric constant (epsilon (parallel to)) reaching values of about 120 for channels with height 8 angstrom < h < 10 angstrom. With the increase in the width of the channel, we predict that epsilon (parallel to) decreases nonlinearly and reaches the bulk value for h > 70 angstrom. A stratified continuum model is proposed that reproduces the h > 10 angstrom dependence of epsilon (parallel to). For sub-nanometer height channels, abnormal behavior of epsilon (parallel to) is found with two orders of magnitude reduction of epsilon (parallel to) around h similar to 7.5 angstrom, which is attributed to the formation of a particular ice phase that exhibits long-time (similar to mu s) stable ferroelectricity. This is of particular importance for the understanding of the influence of confined water on the functioning of biological systems.