Abstract: We investigate the effect of an external magnetic field on the carrier states that are localized at a potential kink and a kink-antikink in bilayer graphene. These chiral states are localized at the interface between two potential regions with opposite signs.

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: In a recent paper [Phys. Rev. B 89, 155403 (2014)], the authors investigated the spectrum of a Coulomb impurity in graphene in the presence of magnetic and electric fields using the coupled series expansion approach. In the first part of their paper, they investigated how Coulomb impurity states collapse in the presence of a perpendicular magnetic field. We argue that the obtained spectrum does not give information about the atomic collapse and that their interpretation of the spectrum regarding atomic collapse is not correct. We also argue that the obtained results are only valid up to the dimensionless charge vertical bar alpha vertical bar = 0.5 and, to obtain correct results for alpha > 0.5, a proper regularization of the Coulomb interaction is required. Here we present the correct numerical results for the spectrum for arbitrary values of alpha.

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: 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 electronic and structural properties of excitons and trions in monolayer transition metal dichalcogenides are investigated using both a multiband and a single- band model. In the multiband model we construct the excitonic Hamiltonian in the product base of the single-particle states at the conduction and valence band edges. We decouple the corresponding energy eigenvalue equation and solve the resulting differential equation self-consistently, using the finite element method (FEM), to determine the energy eigenvalues and the wave functions. As a comparison, we also consider the simple single-band model which is often used in numerical studies. We solve the energy eigenvalue equation using the FEM as well as with the stochastic variational method (SVM) in which a variational wave function is expanded in a basis of a large number of correlated Gaussians. We find good agreement between the results of both methods, as well as with other theoretical works for excitons, and we also compare with available experimental data. For trions the agreement between both methods is not as good due to our neglect of angular correlations when using the FEM. Finally, when comparing the two models, we see that the presence of the valence bands in the mutiband model leads to differences with the single- band model when (interband) interactions are strong.

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: 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: We theoretically investigate the ground-state properties of a quantum dot defined on the surface of a strong three-dimensional time-reversal invariant topological insulator. Confinement is realized by ferromagnetic barriers and Coulomb interaction is treated numerically for up to seven electrons in the dot. Experimentally relevant intermediate interaction strengths are considered. The topological origin of the dot has several consequences: (i) spin polarization increases and the ground state exhibits quantum phase transitions at specific angular momenta as a function of interaction strength, (ii) the onset of Wigner correlations takes place mainly in one spin channel, and (iii) the ground state is characterized by a robust persistent current that changes sign as a function of the distance from the center of the dot.

Abstract: Using the tight-binding approach, we investigate confined states in two different hybrid monolayer-bilayer systems: (i) a hexagonal monolayer area surrounded by bilayer graphene in the presence of a perpendicularly applied electric field and (ii) a hexagonal bilayer graphene dot surrounded by monolayer graphene. The dependence of the energy levels on dot size and external magnetic field is calculated. We find that the energy spectrum for quantum dots with zigzag edges consists of states inside the gap which range from dot-localized states, edge states, to mixed states coexisting together, whereas for dots with armchair edges, only dot-localized states are observed.

Abstract: Using the tight-binding approach, we study the influence of a perpendicular magnetic field on the energy levels of hexagonal, triangular, and circular bilayer graphene (BLG) quantum dots (QDs) with zigzag and armchair edges. We obtain the energy levels for AB (Bernal)-stacked BLG QDs in both the absence and the presence of a perpendicular electric field (i.e., biased BLG QDs). We find different regions in the spectrum of biased QDs with respect to the crossing point between the lowest-electron and -hole Landau levels of a biased BLG sheet. Those different regions correspond to electron states that are localized at the center, edge, or corner of the BLG QD. Quantum Hall corner states are found to be absent in circular BLG QDs. The spatial symmetry of the carrier density distribution is related to the symmetry of the confinement potential, the position of zigzag edges, and the presence or absence of interlayer inversion symmetry.

Abstract: We present a systematic study of the energy spectra of graphene quantum rings having different geometries and edge types in the presence of a perpendicular magnetic field. Results are obtained within the tight-binding (TB) and Dirac models and we discuss which features of the former can be recovered by using the approximations imposed by the latter. Energy levels of graphene quantum rings obtained by diagonalizing the TB Hamiltonian are demonstrated to be strongly dependent on the rings geometry and the microscopical structure of the edges. This makes it difficult to recover those spectra by the existing theories that are based on the continuum (Dirac) model. Nevertheless, our results show that both approaches (i.e., TB and Dirac model) may provide similar results, but only for very specific combinations of ring geometry and edge types. The results obtained by a simplified model describing an infinitely thin circular Dirac ring show good agreement with those obtained for hexagonal and rhombus armchair graphene rings within the TB model. Moreover, we show that the energy levels of a circular quantum ring with an infinite mass boundary condition obtained within the Dirac model agree with those for a ring defined by a ring-shaped staggered potential obtained within the TB model.

Abstract: We investigate the magnetic interface states of graphene quantum dots that contain p-n junctions. Within a tight-binding approach, we consider rectangular quantum dots in the presence of a perpendicular magnetic field containing p-n as well as p-n-p and n-p-n junctions. The results show the interplay between the edge states associated with the zigzag terminations of the sample and the snake states that arise at the p-n junction due to the overlap between electron and hole states at the potential interface. Remarkable localized states are found at the crossing of the p-n junction with the zigzag edge having a dumb-bell-shaped electron distribution. The results are presented as a function of the junction parameters and the applied magnetic flux. DOI: 10.1103/PhysRevB.87.035426

Abstract: Magnetotransport in bilayer graphene in the presence of a weak and periodic potential is investigated in the presence of a perpendicular magnetic field B. The modulation broadens the Landau levels into bands and for weak magnetic fields leads to the well-known Weiss oscillations in their bandwidth and their transport coefficients at very low B and to the Shubnikov-de Haas oscillations at larger B. The amplitude of the Weiss oscillations is severely reduced if the periodic potentials applied to the two layers oscillate out of phase. We also contrast some results with those corresponding to single-layer graphene. Relative to them the flat-band condition and the oscillation amplitude differ substantially, due to the interlayer coupling, and agree only when this coupling is extremely weak. We further show that the Hall conductivity exhibits the well-known steps at half-integer and integer multiples of 4e(2)/h in single-layer and bilayer graphene, respectively, even for very weak magnetic fields. The results are pertinent to weak and periodic corrugations when the potential modulation dominates the strain-induced magnetic modulation.

Abstract: An analytical approach, using the Dirac-Weyl equation, is implemented to obtain the energy spectrum and optical absorption of a circular graphene quantum dot in the presence of an external magnetic field. Results are obtained for the infinite-massand zigzag boundary conditions. We found that the energy spectrum of a dot with the zigzag boundary condition exhibits a zero-energy band regardless of the value of the magnetic field, while for the infinite-mass boundary condition, the zero-energy states appear only for high magnetic fields. The analytical results are compared to those obtained from the tight-binding model: (i) we show the validity range of the continuum model and (ii) we find that the continuum model with the infinite-mass boundary condition describes rather well its tight-binding analog, which can be partially attributed to the blurring of the mixed edges by the staggered potential.