Abstract: Recently proposed accurate correlation energies are used to determine the phase diagram of strongly coupled electron-hole graphene bilayers. The control parameters of the phase diagram are the charge carrier density and the insulating barrier thickness separating the bilayers. In addition to the electron-hole superfluid phase we find two new inhomogeneous ground states, a one dimensional charge density wave phase and a coupled electron-hole Wigner crystal. The elementary crystal structure of bilayer graphene plays no role in generating these new quantum phases, which are completely determined by the electrons and holes interacting through the Coulomb interaction. The experimental parameters for the new phases lie within attainable ranges and therefore coupled electron-hole bilayer graphene presents itself as an experimental system where novel emergent many-body phases can be realized.

Abstract: Starting from a Hamiltonian with anharmonic coupling between in-plane acoustic displacements and outof-plane (flexural) modes, we derived coupled equations of motion for in-plane displacements correlations and flexural mode density fluctuations. Linear response theory and time-dependent thermal Green's functions techniques are applied in order to obtain different response functions. As external perturbations we allow for stresses and thermal heat sources. The displacement correlations are described by a Dyson equation where the flexural density distribution enters as an additional perturbation. The flexural density distribution satisfies a kinetic equation where the in-plane lattice displacements act as a perturbation. In the hydrodynamic limit this system of coupled equations is at the basis of a unified description of elastic and thermal phenomena, such as isothermal versus adiabatic sound motion and thermal conductivity versus second sound. The general theory is formulated in view of application to graphene, two-dimensional h-BN, and 2H-transition metal dichalcogenides and oxides.

Abstract: We propose a tight-binding (TB) model, that includes spin-orbit coupling (SOC), to describe the electronic properties of methyl-substituted germanane (GeCH3). This model gives an electronic spectrum in agreement with first principle results close to the Fermi level. Using the Z(2) formalism, we show that a topological phase transition from a normal insulator (NI) to a quantum spin Hall (QSH) phase occurs at 11.6% biaxial tensile strain. The sensitivity of the electronic properties of this system on strain, in particular its transition to the topological insulating phase, makes it very attractive for applications in strain sensors and other microelectronic applications.

Abstract: The electron energy spectrum in monolayer-bilayer-monolayer and in bilayer-monolayer-bilayer graphene structures is investigated and the effects of a perpendicular magnetic field and electric bias are studied. Different types of monolayer-bilayer interfaces are considered as zigzag (ZZ) or armchair (AC) junctions which modify considerably the bulk Landau levels (LLs) when the spectra are plotted as a function of the center coordinate of the cyclotron orbit. Far away from the two interfaces, one obtains the well-known LLs for extended monolayer or bilayer graphene. The LL structure changes significantly at the two interfaces or junctions where the valley degeneracy is lifted for both types of junctions, especially when the distance between them is approximately equal to the magnetic length. Varying the nonuniform bias and the width of this junction-to-junction region in either structure strongly influence the resulting spectra. Significant differences exist between ZZ and AC junctions in both structures. The densities of states (DOSs) for unbiased structures are symmetric in energy whereas those for biased structures are asymmetric. An external bias creates interface LLs in the gaps between the LLs of the unbiased system in which the DOS can be quite small. Such a pattern of LLs can be probed by scanning tunneling microscopy.

Abstract: The energy spectrum and eigenstates of single-layer black phosphorus nanoribbons in the presence of a perpendicular magnetic field and an in-plane transverse electric field are investigated by means of a tight-binding method, and the effect of different types of edges is examined analytically. A description based on a continuum model is proposed using an expansion of the tight-binding model in the long-wavelength limit. Thewave functions corresponding to the flatband part of the spectrum are obtained analytically and are shown to agree well with the numerical results from the tight-binding method for both narrow (10 nm) and wide (100 nm) nanoribbons. Analytical expressions for the critical magnetic field at which Landau levels are formed and the ranges of wave numbers in the dispersionless flatband segments in the energy spectra are derived. We examine the evolution of the Landau levels when an in-plane lateral electric field is applied, and we determine analytically how the edge states shift withmagnetic field. For wider nanoribbons, the conductance is shown to have a characteristic staircase shape in combined magnetic and electric fields. Some of the stairs in zigzag and skewed armchair nanoribbons originate from edge states that are found in the band gap.

Abstract: <script type='text/javascript'>document.write(unpmarked('Two-dimensional materials are known to harbour properties very different from those of their bulk counterparts. Recent years have seen the rise of atomically thin superconductors, with a caveat that superconductivity is strongly depleted unless enhanced by specific substrates, intercalants or adatoms. Surprisingly, the role in superconductivity of electronic states originating from simple free surfaces of two-dimensional materials has remained elusive to date. Here, based on first-principles calculations, anisotropic Eliashberg theory, and angle-resolved photoemission spectroscopy (ARPES), we show that surface states in few-monolayer MgB2 make a major contribution to the superconducting gap spectrum and density of states, clearly distinct from the widely known, bulk-like sigma-and pi-gaps. As a proof of principle, we predict and measure the gap opening on the magnesium-based surface band up to a critical temperature as high as similar to 30 K for merely six monolayers thick MgB2. These findings establish free surfaces as an unavoidable ingredient in understanding and further tailoring of superconductivity in atomically thin materials.'));

Abstract: <script type='text/javascript'>document.write(unpmarked('The motion of Abrikosov vortices in type-II superconductors results in a finite resistance in the presence of an applied electric current. Elimination or reduction of the resistance via immobilization of vortices is the \u0022holy grail\u0022 of superconductivity research. Common wisdom dictates that an increase in the magnetic field escalates the loss of energy since the number of vortices increases. Here we show that this is no longer true if the magnetic field and the current are applied parallel to each other. Our experimental studies on the resistive behavior of a superconducting Mo0.79Ge0.21 nanostrip reveal the emergence of a dissipative state with increasing magnetic field, followed by a pronounced resistance drop, signifying a reentrance to the superconducting state. Large-scale simulations of the 3D time-dependent Ginzburg-Landau model indicate that the intermediate resistive state is due to an unwinding of twisted vortices. When the magnetic field increases, this instability is suppressed due to a better accommodation of the vortex lattice to the pinning configuration. Our findings show that magnetic field and geometrical confinement can suppress the dissipation induced by vortex motion and thus radically improve the performance of superconducting materials.'));

Abstract: The self rotation of a graphene flake over graphite is controlled by the size, initial misalignment and temperature. Using both ab initio calculations and molecular dynamics simulations, we investigate annealing effects on the self rotation of a graphene flake on a graphene substrate. The energy barriers for rotation and drift of a graphene flake over graphene is found to be smaller than 25 meV/atom which is comparable to thermal energy. We found that small flakes (of about similar to 4 nm) are more sensitive to temperature and initial misorientation angles than larger one (beyond 10 nm). The initial stacking configuration of the flake is found to be important for its dynamics and time evolution of misalignment. Large flakes, which are initially in the AA-or AB-stacking state with small misorientation angle, rotate and end up in the AB-stacking configuration. However small flakes can they stay in an incommensurate state specially when the initial misorientation angle is larger than 2 degrees. Our results are in agreement with recent experiments.

Abstract: Two-dimensional crystals with strong interactions between layers has attracted increasing attention in recent years in a variety of fields. In particular, the growth of a single layer of oxide materials (e.g., MgO, CaO, NiO, and MnO) over metallic substrates were found to display different physical properties than their bulk. In this study, we report on the physical properties of a single layer of metallic oxide materials and compare their properties with their bulk and other two-dimensional (2D) crystals. We found that the planar structure of metallic monoxides are unstable whereas the buckled structures are thermodynamically stable. Also, the 2D-MnO and NiO exhibit different magnetic (ferromagnetic) and optical properties than their bulk, whereas band-gap energy and linear stiffness are found to be decreasing from NiO to MgO. Our findings provide insight into oxide thin-film technology applications.

Abstract: The interaction of monolayer graphene with specific substrates may break its sublattice symmetry and results in unidirectional chiral states with opposite group velocities in the different Dirac cones (Zarenia et al 2012 Phys. Rev. B 86 085451). Taking advantage of this feature, we propose a valley filter based on a transversal mass kink for low energy electrons in graphene, which is obtained by assuming a defect region in the substrate that provides a change in the sign of the substrate-induced mass and thus creates a non-biased channel, perpendicular to the kink, for electron motion. By solving the time-dependent Schrodinger equation for the tight-binding Hamiltonian, we investigate the time evolution of a Gaussian wave packet propagating through such a system and obtain the transport properties of this graphene-based substrate-induced quantum point contact. Our results demonstrate that efficient valley filtering can be obtained, provided: (i) the electron energy is sufficiently low, i.e. with electrons belonging mostly to the lowest sub-band of the channel, and (ii) the channel length (width) is sufficiently long (narrow). Moreover, even though the transmission probabilities for each valley are significantly affected by impurities and defects in the channel region, the valley polarization in this system is shown to be robust against their presence.

Abstract: <script type='text/javascript'>document.write(unpmarked('We show that delaminations in bilayer graphene (BLG) with electrostatically induced interlayer symmetry can provide one with ballistic channels for electrons with energies inside the electrostatically induced BLG gap. These channels are formed by a combination of valley-polarized evanescent states propagating along the delamination edges (which persist in the presence of a strong magnetic field) and standing waves bouncing between them inside the delaminated region (in a strong magnetic field, these transform into Landau levels in the monolayers). For inverted stackings in BLGs on the left and right of the delamination (AB-2ML-BA or BA-2ML-AB, where 2ML indicates two decoupled monolayers of graphene), the lowest-energy ballistic channels are gapless, have linear dispersion, and appear to be weakly topologically protected. When BLG stackings on both sides of the delamination are the same (AB-2ML-AB or BA-2ML-BA), the lowest-energy ballistic channels are gapped, with a gap epsilon(g) scaling as epsilon(g) alpha W-1 with delamination width and epsilon(g) alpha delta(-1) with the on-layer energy difference in the delaminated part of the structure. Depending on the width, delaminations may also support several \u0022higher-energy\u0022 waveguide modes. Our results are based on both the analytical study of the wave matching of Dirac states and tight-binding model calculations, and we analyze in detail the dependence of the delamination spectrum on the electrostatic conditions in the structure, such as the vertical displacement field.'));

Abstract: <script type='text/javascript'>document.write(unpmarked('Quantum electrodynamics predicts that heavy atoms (Z \u003E Z(c) approximate to 170) will undergo the process of atomic collapse where electrons sink into the positron continuum and a new family of so-called collapsing states emerges. The relativistic electrons in graphene exhibit the same physics but at a much lower critical charge (Z(c) approximate to 1) which has made it possible to confirm this phenomenon experimentally. However, there exist conflicting predictions on the effect of a magnetic field on atomic collapse. These theoretical predictions are based on the continuum Dirac-Weyl equation, which does not have an exact analytical solution for the interplay of a supercritical Coulomb potential and the magnetic field. Approximative solutions have been proposed, but because the two effects compete on similar energy scales, the theoretical treatment varies depending on the regime which is being considered. These limitations are overcome here by starting from a tight-binding approach and computing exact numerical results. By avoiding special limit cases, we found a smooth evolution between the different regimes. We predict that the atomic collapse effect persists even after the magnetic field is activated and that the critical charge remains unchanged. We show that the atomic collapse regime is characterized: (1) by a series of Landau level anticrossings and (2) by the absence of root B scaling of the Landau levels with regard to magnetic field strength.'));

Abstract: Using the tight-binding approach, we investigate the energy spectrum of square, triangular, and hexagonal MoS2 quantum dots (QDs) in the presence of a perpendicular magnetic field. Novel edge states emerge in MoS2 QDs, which are distributed over the whole edge which we call ring states. The ring states are robust in the presence of spin-orbit coupling (SOC). The corresponding energy levels of the ring states oscillate as a function of the perpendicular magnetic field which are related to Aharonov-Bohm oscillations. Oscillations in the magnetic field dependence of the energy levels and the peaks in the magneto-optical spectrum emerge (disappear) as the ring states are formed (collapsed). The period and the amplitude of the oscillation decrease with the size of the MoS2 QDs.

Abstract: MXenes are the newest class of two-dimensional (2D) materials, and they offer great potential in a wide range of applications including electronic devices, sensors, and thermoelectric and energy storage materials. In this work, we combined the outstanding electrical conductivity, that is essential for battery applications, of graphene with MXene monolayers (M2CX2 where M = Sc, Ti, V and X = OH, O) to explore its potential in Li battery applications. Through first principles calculations, we determined the stable stacking configurations of M2CX2/graphene bilayer heterostructures and their Li atom intercalation by calculating the Li binding energy, diffusion barrier and voltage. We found that: (1) for the ground state stacking, the interlayer binding is strong, yet the interlayer friction is small; (2) Li binds more strongly to the O-terminated monolayer, bilayer and heterostructure MXene systems when compared with the OHterminated MXenes due to the H+ induced repulsion to the Li atoms. The binding energy of Li decreases as the Li concentration increases due to enhanced repulsive interaction between the positively charged Li ions; (3) Ti2CO2/graphene and V2CO2/graphene heterostructures exhibit large Li atom binding energies making them the most promising candidates for battery applications. When fully loaded with Li atoms, the binding energy is -1.43 eV per Li atom and -1.78 eV per Li atom for Ti2CO2/graphene and V2CO2/graphene, respectively. These two heterostructures exhibit a nice compromise between storage capacity and kinetics. For example, the diffusion barrier of Li in Ti2CO2/graphene is around 0.3 eV which is comparable to that of graphite. Additionally, the calculated average voltages are 1.49 V and 1.93 V for Ti2CO2/graphene and V2CO2/graphene structures, respectively; (4) a small change in the in-plane lattice parameters (<1%), interatomic bond lengths and interlayer distances (<0.5 angstrom) proves the stability of the heterostructures against Li intercalation, and the impending phase separation into constituent layers and capacity fading during charge-discharge cycles in real battery applications; (5) as compared to bare M2CX2 bilayers, M2CX2/graphene heterostructures have lower molecular mass, offering high storage capacity; (6) the presence of graphene ensures good electrical conductivity that is essential for battery applications. Given these advantages, Ti2CO2/graphene and V2CO2/graphene heterostructures are predicted to be promising for lithium-ion battery applications.

Abstract: The dependence of the excitonic photoluminescence (PL) spectrum of monolayer transition metal dichalcogenides (TMDs) on the tilt angle of an applied magnetic field is studied. Starting from a four-band Hamiltonian we construct a theory which quantitatively reproduces the available experimental PL spectra for perpendicular and in-plane magnetic fields. In the presence of a tilted magnetic field, we demonstrate that the dark exciton PL peaks brighten due to the in-plane component of the magnetic field and split for light with different circular polarizations as a consequence of the perpendicular component of the magnetic field. This splitting is more than twice as large as the splitting of the bright exciton peaks in tungsten-based TMDs. We propose an experimental setup that will allow for accessing the predicted splitting of the dark exciton peaks in the PL spectrum.