Abstract: By taking account of the electric-field-induced charge screening, a self-consistent calculation within the framework of the tight-binding approach is employed to obtain the electronic band structure of gated multilayer phosphorene and the charge densities on the different phosphorene layers. We find charge density and screening anomalies in single-gated multilayer phosphorene and electron-hole bilayers in dual-gated multilayer phosphorene. Due to the unique puckered lattice structure, both intralayer and interlayer charge screenings are important in gated multilayer phosphorene. We find that the electric-field tuning of the band structure of multilayer phosphorene is distinctively different in the presence and absence of charge screening. For instance, it is shown that the unscreened band gap of multilayer phosphorene decreases dramatically with increasing electric-field strength. However, in the presence of charge screening, the magnitude of this band-gap decrease is significantly reduced and the reduction depends strongly on the number of phosphorene layers. Our theoretical results of the band-gap tuning are compared with recent experiments and good agreement is found.

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: Excitonic superfluidity in double phosphorene monolayers is investigated using the BCS mean-field equations. Highly anisotropic superfluidity is predicted where we found that the maximum superfluid gap is in the Bose-Einstein condensate (BEC) regime along the armchair direction and in the BCS-BEC crossover regime along the zigzag direction. We estimate the highest Kosterlitz-Thouless transition temperature with maximum value up to similar to 90 K with onset carrier densities as high as 4 x 10(12) cm(-2). This transition temperature is significantly larger than what is found in double electron-hole few-layers graphene. Our results can guide experimental research toward the realization of anisotropic condensate states in electron-hole phosphorene monolayers.

Abstract: We use a variational technique to study the many-body phase diagram of electrons and holes in n-doped and p-doped monolayer transition metal dichalcogenides (TMDs). We find a total of four different phases. (i) A fully spin polarized and valley polarized ferromagnetic state. (ii) A state with no global spin polarization but with spin polarization in each valley separately, i.e., spin-valley locking. (iii) A state with spin polarization in one of the valleys and little to no spin polarization in the other valley. (iv) A paramagnetic state with no valley polarization. These phases are separated by first-order phase transitions and are determined by the particle density and the dielectric constant of the substrate. We find that in the presence of a perpendicular magnetic field the four different phases persist. In the case of n-doped MoS2, a fifth phase, which is completely valley polarized but not spin polarized, appears for magnetic fields larger than 7 T and for magnetic fields larger than 23 T completely replaces the second phase.

Abstract: Starting from the single-particle Dirac Hamiltonian for charge carriers in monolayer transition metal dichalcogenides (TMDs), we construct a four-band Hamiltonian describing interlayer excitons consisting of an electron in one TMD layer and a hole in the other TMD layer. An expression for the electron-hole interaction potential is derived, taking into account the effect of the dielectric environment above, below, and between the two TMD layers as well as polarization effects in the transition metal layer and in the chalcogen layers of the TMD layers. We calculate the interlayer exciton binding energy and average in-plane interparticle distance for different TMD heterostructures. The effect of different dielectric environments on the exciton binding energy is investigated and a remarkable dependence on the dielectric constant of the barrier between the two layers is found, resulting from competing effects as a function of the in-plane and out-of-plane dielectric constants of the barrier. The polarization effects in the chalcogen layers, which in general reduce the exciton binding energy, can lead to an increase in binding energy in the presence of strong substrate effects by screening the substrate. The excitonic absorbance spectrum is calculated and we show that the interlayer exciton peak depends linearly on a perpendicular electric field, which agrees with recent experimental results.

Abstract: Using first-principles calculations, we study the dependence of the electronic and vibrational properties of multilayered PbI2 crystals on the number of layers and focus on the electronic-band structure and the Raman spectrum. Electronic-band structure calculations reveal that the direct or indirect semiconducting behavior of PbI2 is strongly influenced by the number of layers. We find that at 3L thickness there is a direct-to-indirect band gap transition (from bulk-to-monolayer). It is shown that in the Raman spectrum two prominent peaks, A(1g) and E-g, exhibit phonon hardening with an increasing number of layers due to the interlayer van der Waals interaction. Moreover, the Raman activity of the A(1g) mode significantly increases with an increasing number of layers due to the enhanced out-of-plane dielectric constant in the few-layer case. We further characterize rigid-layer vibrations of low-frequency interlayer shear (C) and breathing (LB) modes in few-layer PbI2. A reduced monatomic (linear) chain model (LCM) provides a fairly accurate picture of the number of layers dependence of the low-frequency modes and it is shown also to be a powerful tool to study the interlayer coupling strength in layered PbI2.

Abstract: In the quest for miniaturizing magnetic devices, the effects of boundaries and surfaces become increasingly important. Here we show how the recently predicted boundary-induced Dzyaloshinskii-Moriya interaction (DMI) affects the magnetization of ferromagnetic films with a C-infinity v symmetry and a perpendicular magnetic anisotropy. For an otherwise uniformly magnetized film, we find a surface twist when the magnetization in the bulk is canted by an in-plane external field. This twist at the surfaces caused by the boundary-induced DMI differs from the common canting caused by internal DMI observed at the edges of a chiral magnet. Furthermore, we find that the surface twist due to the boundary-induced DMI strongly affects the width of the domain wall at the surfaces. We also find that the skyrmion radius increases in the depth of the film, with the average size of the skyrmion increasing with boundary-induced DMI. This increase suggests that the boundary-induced DMI contributes to the stability of the skyrmion.

Abstract: The charge density is among the most fundamental solid state properties determining bonding, electrical characteristics, and adsorption or catalysis at surfaces. While atomic-scale charge densities have as yet been retrieved by solid state theory, we demonstrate both charge density and electric field mapping across a mono-/bilayer boundary in 2D MoS2 by momentum-resolved scanning transmission electron microscopy. Based on consistency of the four-dimensional experimental data, statistical parameter estimation and dynamical electron scattering simulations using strain-relaxed supercells, we are able to identify an AA-type bilayer stacking and charge depletion at the Mo-terminated layer edge.

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: We demonstrate that combining standing wave (SW) excitation with resonant inelastic x-ray scattering (RIXS) can lead to depth resolution and interface sensitivity for studying orbital and magnetic excitations in correlated oxide heterostructures. SW-RIXS has been applied to multilayer heterostructures consisting of a superconductor La1.85Sr0.15CuO4 (LSCO) and a half-metallic ferromagnet La0.67Sr0.33MnO3 (LSMO). Easily observable SW effects on the RIXS excitations were found in these LSCO/LSMO multilayers. In addition, we observe different depth distribution of the RIXS excitations. The magnetic excitations are found to arise from the LSCO/LSMO interfaces, and there is also a suggestion that one of the dd excitations comes from the interfaces. SW-RIXS measurements of correlated-oxide and other multilayer heterostructures should provide unique layer-resolved insights concerning their orbital and magnetic excitations, as well as a challenge for RIXS theory to specifically deal with interface effects.

Abstract: We discuss positron-annihilation lifetimes for a set of illustrative bulk materials within the framework of the weighted-density approximation (WDA). The WDA can correctly describe electron-positron correlations in strongly inhomogeneous systems, such as surfaces, where the applicability of (semi-)local approximations is limited. We analyze the WDA in detail and show that the electrons which cannot screen external charges efficiently, such as the core electrons, cannot be treated accurately via the pair correlation of the homogeneous electron gas. We discuss how this problem can be addressed by reducing the screening in the homogeneous electron gas by adding terms depending on the gradient of the electron density. Further improvements are obtained when core electrons are treated within the LDA and the valence electron using the WDA. Finally, we discuss a semiempirical WDA-based approach in which a sum rule is imposed to reproduce the experimental lifetimes.

Abstract: The Dzyaloshinskii-Moriya interaction (DMI) is a chiral interaction that favors formation of domain walls. Recent experiments and ab initio calculations show that there are multiple ways to modify the strength of the interfacially induced DMI in thin ferromagnetic films with perpendicular magnetic anisotropy. In this paper we reveal theoretically the effects of spatially varied DMI on the magnetic state in thin films. In such heterochiral 2D structures we report several emergent phenomena, ranging from the equilibrium spin canting at the interface between regions with different DMI, over particularly strong confinement of domain walls and skyrmions within high-DMI tracks, to advanced applications such as domain tailoring nearly at will, design of magnonic waveguides, and much improved skyrmion racetrack memory.

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 investigate the effects related to vortex-core deformations when vortices approach each other. As a result of these vortex-core deformations, the vortex-vortex interaction effectively acquires an attractive component leading to a variety of vortex patterns typical for systems with nonmonotonic repulsive-attractive interaction, such as stripes and labyrinths. The core deformations are anisotropic and can induce frustration in the vortex-vortex interaction. In turn, this frustration has an impact on the resulting vortex patterns, which are analyzed in the presence of additional random pinning, as a function of the pinning strength. This analysis can be applicable to vortices in multiband superconductors or to vortices in Bose-Einstein condensates.

Abstract: We investigate scattering of the topological surface state of a three-dimensional time-reversal invariant topological insulator when graphene is deposited on the topological-insulator surface. Specifically, we consider the (111) surface of a Bi2Se3-like topological insulator. We present a low-energy model for the graphene-topological insulator heterostructure and we calculate the transmission probability at zigzag and armchair edges of the deposited graphene, and the conductance through graphene nanoribbon barriers, and show that its features can be understood from antiresonances in the transmission probability.