Abstract: Edge functionalization of graphene nanoribbons with nitrogen atoms for various adatom configurations at armchair and zigzag edges are investigated. We provide comprehensive information on the electronic and magnetic properties and investigate the stability of the various systems. Two types of rippling of the nanoribbons, namely edge and bulk rippling depending on the sign of edge stress induced at the edge, are found. They are found to play the decisive role for the stability of the structures. We also propose a type of edge decoration in which every third nitrogen adatom at the zigzag edges is replaced by an oxygen atom. In this way, the electron count is compatible with a full aromatic structure, leading to additional stability and a disappearance of magnetism that is usually associated with zigzag nanoribbons.

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: We investigate the energy spectrum, wave functions, and local density of states of an electrical dipole placed on a sheet of gapped graphene as function of the charge strength Z alpha for different sizes of the dipole and for different regularization parameters. The dipole is modeled as consisting of a positive and negative charge. Bound states are found within the gap region with some energy levels that anticross and others that cross as function of the impurity strength Z alpha. The anticrossings are more pronounced and move to higher charges Z alpha when the length of the dipole decreases. These energy levels turn into atomic collapse states when they enter the positive (or negative) energy continuum. A smooth transition from the single-impurity behavior to the dipole one is observed: The states diving towards the continuum in the single-impurity case are gradually replaced by a series of anticrossings that represent a continuation of the diving states in the single-impurity case. By studying the local density of states at the edge of the dipole we show how the series of anticrossings persist in the positive and negative continuum.

Abstract: Bilayer graphene samples may exhibit regions where the two layers are locally delaminated forming a so-called quanttun blister in the graphene sheet. Electron and hole states can be confined in this graphene quantum blisters (GQB) by applying a global electrostatic bias. We scrutinize the electronic properties of these confined states under the variation of interlayer bias, coupling, and blister's size. The spectra display strong anti-crossings due to the coupling of the confined states on upper and lower layers inside the blister. These spectra are layer localized where the respective confined states reside on either layer or equally distributed. For finite angular momentum, this layer localization can be at the edge of the blister and corresponds to degenerate modes of opposite momenta. Furthermore, the energy levels in GQB exhibit electron-hole symmetry that is sensitive to the electrostatic bias. Finally, we demonstrate that confinement in GQB persists even in the presence of a variation in the interlayer coupling.

Abstract: The possibility of an inhomogeneous charge density wave phase is investigated in a system of two coupled electron and hole monolayers separated by a hexagonal boron nitride insulating layer. The charge-density-wave state is induced through the assumption of negative compressibility of electron/hole gases in a Coulomb drag configuration between the electron and hole sheets. Under equilibrium conditions, we derive analytical expressions for the density oscillation along the zigzag and armchair directions. We find that the density modulation not only depends on the sign of the compressibility but also on the anisotropy of the low-energy bands. Our results are applicable to any two-dimensional system with anisotropic parabolic bands, characterized by different effective masses. For equal effective masses, i.e., isotropic energy bands, our results agree with Hroblak et al. [Phys. Rev. B 96, 075422 (2017)]. Our numerical results are applied to phosphorene.

Abstract: The effect of anisotropy in the energy spectrum on the binding energy and structural properties of excitons, trions, and biexcitons is investigated. To this end we employ the stochastic variational method with a correlated Gaussian basis. We present results for the binding energy of different excitonic complexes in black phosphorus (bP) and TiS3 and compare them with recent results in the literature when available, for which we find good agreement. The binding energies of excitonic complexes in bP are larger than those in TiS3. We calculate the different average interparticle distances in bP and TiS3 and show that excitonic complexes in bP are strongly anisotropic whereas in TiS3 they are almost isotropic, even though the constituent particles have an anisotropic energy spectrum. This is also confirmed by the correlation functions.

Abstract: The effects of strain, induced by a Gaussian bump, on the magnetic field dependent transport properties of a graphene Hall bar are investigated. The numerical simulations are performed using both classical and quantum mechanical transport theory and we found that both approaches exhibit similar characteristic features. The effects of the Gaussian bump are manifested by a decrease of the bend resistance, RB, around zero-magnetic field and the occurrence of side-peaks in RB. These features are explained as a consequence of bump-assisted scattering of electrons towards different terminals of the Hall bar. Using these features we are able to give an estimate of the size of the bump. Additional oscillations in RB are found in the quantum description that are due to the population/depopulation of Landau levels. The bump has a minor influence on the Hall resistance even for very high values of the pseudo-magnetic field. When the bump is placed outside the center of the Hall bar valley polarized electrons can be collected in the leads.

Abstract: The magnetic-field dependence of the energy spectrum, wave function, binding energy, and oscillator strength of exciton states confined in a circular graphene quantum dot (CGQD) is obtained within the configuration interaction method. We predict that (i) excitonic effects are very significant in the CGQD as a consequence of a combination of geometric confinement, magnetic confinement, and reduced screening; (ii) two types of excitons (intravalley and intervalley) are present in the CGQD because of the valley degree of freedom in graphene; (iii) the intravalley and intervalley exciton states display different magnetic-field dependencies due to the different electron-hole symmetries of the single-particle energy spectra; (iv) with increasing magnetic field, the exciton ground state in the CGQD undergoes an intravalley to intervalley transition accompanied by a change of angular momentum; (v) the exciton binding energy does not increase monotonically with the magnetic field due to the competition between geometric and magnetic confinements; and (vi) the optical transitions of the intervalley and intravalley excitons can be tuned by the magnetic field, and valley-dependent excitonic transitions can be realized in a CGQD.

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: 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: Two-dimensional (2D) carbon nitride materials play an important role in energy-harvesting, energy-storage and environmental applications. Recently, a new carbon nitride, 2D polyaniline (C3N) was proposed [PNAS 113 (2016) 7414-7419]. Based on the structure model of this C3N monolayer, we propose two new carbon nitride monolayers, named dumbbell (DB) C4N-I and C4N-II. Using first-principles calculations, we systematically study the structure, stability, and band structure of these two materials. In contrast to other carbon nitride monolayers, the orbital hybridization of the C/N atoms in the DB C4N monolayers is sp(3). Remarkably, the band structures of the two DB C4N monolayers have a Dirac cone at the K point and their Fermi velocities (2.6/2.4 x 10(5) m/s) are comparable to that of graphene. This makes them promising materials for applications in high-speed electronic devices. Using a tight-binding model, we explain the origin of the Dirac cone. (C) 2017 Elsevier Ltd. All rights reserved.

Abstract: The elastic constant C-11 and piezoelectric stress constant e(1),(11) of two-dimensional (2D) dielectric materials comprising h-BN, 2H-MoS2, and other transition-metal dichalcogenides and dioxides are calculated using lattice dynamical theory. The results are compared with corresponding quantities obtained with ab initio calculations. We identify the difference between clamped-ion and relaxed-ion contributions with the dependence on inner strains which are due to the relative displacements of the ions in the unit cell. Lattice dynamics allows us to express the inner-strain contributions in terms of microscopic quantities such as effective ionic charges and optoacoustical couplings, which allows us to clarify differences in the piezoelectric behavior between h-BN and MoS2. Trends in the different microscopic quantities as functions of atomic composition are discussed.

Abstract: We compare classical versus quantum electron transport in recently investigated magnetic focusing devices (Bhandari et al 2016 Nano Lett. 16 1690) exposed to the perturbing potential of a scanning gate microscope (SGM). Using the Landauer-Buttiker formalism for a multi-terminal device, we calculate resistance maps that are obtained as the SGM tip is scanned over the sample. There are three unique regimes in which the scanning tip can operate (focusing, repelling, and mixed regime) which are investigated. Tip interacts mostly with electrons with cyclotron trajectories passing directly underneath it, leaving a trail of modified current density behind it. Other (indirect) trajectories become relevant when the tip is placed near the edges of the sample, and current is scattered between the tip and the edge. We point out that, in contrast to SGM experiments on gapped semiconductors, the STM tip can induce a pn junction in graphene, which improves contrast and resolution in SGM. We also discuss possible explanations for spatial asymmetry of experimentally measured resistance maps, and connect it with specific configurations of the measuring probes.

Abstract: We introduce a new approach for the correlation energy of one- and two-valley two-dimensional electron gas (2DEG) systems. Our approach is based on an interpolation between two limits, a random phase approximation at high densities and a classical approach at low densities which gives excellent agreement with available Quantum Monte Carlo (QMC) calculations. The two-valley 2DEG model is introduced to describe the electron correlations in monolayer transition metal dichalcogenides (TMDs). We study the zero-temperature transition from a Fermi liquid to a quantum Wigner crystal phase in monolayer TMDs. Consistent with QMC, we find that electrons crystallize at r(s) = 31 in one-valley 2DEG. For two valleys, we predict Wigner crystallization at r(s) = 30, implying that valley degeneracy has little effect on the critical r(s), in contrast to an earlier claim.

Abstract: Spin transitions driven by a periodically varying electric potential in dilute fluorinated graphene quantum dots are investigated. Flakes of monolayer graphene as well as electrostatic electron traps induced in bilayer graphene are considered. The stationary states obtained within the tight-binding approach are used as the basis for description of the system dynamics. The dilute fluorination of the top layer lifts the valley degeneracy of the confined states and attenuates the orbital magnetic dipole moments due to current circulation within the flake. The spin-orbit coupling introduced by the surface deformation of the top layer induced by the adatoms allows the spin flips to be driven by the AC electric field. For the bilayer quantum dots the spin flip times is substantially shorter than the spin relaxation. Dynamical effects including many-photon and multilevel transitions are also discussed.