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: Sodium-ion batteries (SIBs) are potential cost-effective solutions for stationary energy storage applications. Unavailability of suitable anode materials, however, is one of the important barriers to the maturity of SIBs. Here, we report a Na2+xTi4O9/C composite as a promising anode candidate for SIBs with high capacity and cycling stability. This anode is characterized by a capacity of 124 mAh g(-1) (plus 11 mAh g(-1) contributed by carbon black), an average discharge potential of 0.9 V vs Na/Na+, a good rate capability and a high stability (89% capacity retention after 250 cycles at a rate of 1 degrees C). The mechanisms of sodium insertion/deinsertion and of the formation of Na2+xTi4O9/C are investigated with the aid of various ex/in situ characterization techniques. The in situ formed carbon is necessary for the formation of the reduced sodium titanate. This synthesis method may enable the convenient synthesis of other composites of crystalline phases with amorphous carbon.

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: 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: Bilayer hexagonal borophene, which is bound together through pillars, is a novel topological semimetal. Using density functional theory, we investigate its electronic band structure and show that it is a Dirac material which exhibits a nodal line. A tight-binding model was constructed based on the Slater-Koster approach, which accurately models the electronic spectrum. We constructed an effective four-band model Hamiltonian to describe the spectrum near the nodal line. This Hamiltonian can be used as a new platform to study the new properties of nodal line semimetals. We found that the nodal line is created by edge states and is very robust against perturbations and impurities. Breaking symmetries can split the nodal line, but cannot open a gap.

Abstract: The spectrum of a Coulomb impurity in bilayer graphene is investigated as function of the strength of a perpendicular magnetic field for different values of the angular quantum number m and for different values of the gate voltage. We point out fundamental differences between the results from the two-band and four-band model. The supercritical instability and fall-to-center phenomena are investigated in the presence of a magnetic field. We find that in the four-band model the fall-to-center phenomenon occurs as in monolayer graphene, while this is not the case in the two-band model. We find that in a magnetic field the supercritical instability manifests itself as a series of anticrossings in the hole part of the spectrum for states coming from the low-energy band. However, we also find very distinct anticrossings in the electron part of the spectrum that continue into the hole part, which are related to the higher energy band of the four-band model. At these anticrossings, we find a very sharp peak in the probability density close to the impurity, reminiscent for the fall-to-center phenomenon. In this paper, these peculiar and interesting effects are studied for different magnetic field, interlayer coupling, and bias potential strengths.

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: 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: 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 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: Coupled 2D sheets of electrons and holes are predicted to support novel quantum phases. Two experiments of Coulomb drag in electron-hole (e-h) double bilayer graphene (DBLG) have reported an unexplained and puzzling sign reversal of the drag signal. However, we show that this effect is due to the multiband character of DBLG. Our multiband Fermi liquid theory produces excellent agreement and captures the key features of the experimental drag resistance for all temperatures. This demonstrates the importance of multiband effects in DBLG: they have a strong effect not only on superfluidity, but also on the drag.

Abstract: In this study, we predicted two new stable metallic Re-C based monolayer structures with a rectangular (r-ReC2) and a hexagonal (h-Re2C) crystal symmetry using first-principle calculations based on density functional theory. Our results obtained from mechanical and phonon calculations and high-temperature molecular dynamic simulations clearly proved the stability of these two-dimensional (2D) crystals. Interestingly, Re-C monolayers in common transition metal carbide structures (i.e. MXenes) were found to be unstable, contrary to expectations. We found that the stable structures, i.e. r-ReC2 and h-Re2C, display superior mechanical properties over the well-known 2D materials. The Young's modulus for r-ReC2 and h-Re2C are extremely high and were calculated as 351 (1310) and 617 (804) N/m (GPa), respectively. Both materials have larger Young's modulus values than the most of the well-known 2D materials. We showed that the combination of the short strong directional p-d bonds, the high coordination number of atoms in the unit-cell and high valence electron density result in strong mechanical properties. Due to its crystal structure, the r-ReC2 monolayer has anisotropic mechanical properties and the crystallographic direction parallel to the C-2 dimers is stiffer compared to perpendicular direction due to strong covalent bonding within C-2 dimers. h-Re2C was derived from the corresponding bulk structure for which we determined the critical thickness for the dynamically stable bulk-derived monolayer structures. In addition, we also investigated the electronic of these two stable structures. Both exhibit metallic behavior and Re-5d orbitals dominate the states around the Fermi level. Due to their ultra high mechanical stability and stiffness, these novel Re-C monolayers can be exploited in various engineering applications.

Abstract: Electrenes, an atomically thin form of layered electrides, are very recent members of the 2D materials family. In this work, we employed first principle calculations to determine stable, exfoliatable, and application-promising 2D electrene materials among possible M2X compounds, where M is a group II-A metal and X is a nonmetal element (C, N, P, As, and Sb). The promise of stable electrene compounds for battery applications is assessed via their exfoliation energy, adsorption properties, and migration energy barriers toward relevant Li, Na, K, and Ca atoms. Our calculations revealed five new stable electrene candidates in addition to previously known Ca2N and Sr2N. Among these seven dynamically stable electrenes, Ba2As, Ba2P, Ba2Sb, Ca2N, Sr2N, and Sr2P are found to be very promising for either K or Na ion batteries due to their extremely low migration energy barriers (5-16 meV), which roughly demonstrates 105 times higher mobility than graphene and two to four times higher mobility than other promising 2D materials such as MXene (Mo2C).