Abstract: Electrically gating of trilayer graphene (TLG) opens a band gap offering the possibility to electrically engineer TLG quantum dots. We study the energy levels of such quantum dots and investigate their dependence on a perpendicular magnetic field B and different types of stacking of the graphene layers. The dots are modeled as circular and confined by a truncated parabolic potential which can be realized by nanostructured gates or position-dependent doping. The energy spectra exhibit the intervalley symmetry E-K(e) (m) = -E (h)(K') (m) for the electron (e) and hole (h) states, where m is the angular momentum quantum number and K and K' label the two valleys. The electron and hole spectra for B = 0 are twofold degenerate due to the intervalley symmetry E-K (m) = E-K' [-(m + 1)]. For both ABC [alpha = 1.5 (1.2) for large (small) R] and ABA (alpha = 1) stackings, the lowest-energy levels show approximately a R-alpha dependence on the dot radius R in contrast with the 1/R-3 one for ABC-stacked dots with infinite-mass boundary. As functions of the field B, the oscillator strengths for dipole-allowed transitions differ drastically for the two types of stackings.

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: Using the effective continuum model, the electron energy spectrum of gated bilayer graphene with a step-like region of decoupled graphene layers at the edge of the sample is studied. Different types of coupled-decoupled interfaces are considered, i.e., zigzag (ZZ) and armchair junctions, which result in significant different propagating states. Two non-valley-polarized conducting edge states are observed for ZZ type, which are mainly located around the ZZ-ended graphene layers. Additionally, we investigated both BA-BA and BA-AB domain walls in the gated bilayer graphene within the continuum approximation. Unlike the BA-BA domain wall, which exhibits gapped insulating behaviour, the domain walls surrounded by different stackings of bilayer regions feature valley-polarized edge states. Our findings are consistent with other theoretical calculations, such as from the tight-binding model and first-principles calculations, and agree with experimental observations. Published by AIP Publishing.

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 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: Using the tight-binding approach coupled with mean-field Hubbard model, we theoretically study the effect of mechanical deformations on the magnetic properties of bilayer graphene (BLG) quantum dots (QDs). Results are obtained for AA-and AB(Bernal)-stacked BLG QDs, considering different geometries (hexagonal, triangular and square shapes) and edge types (armchair and zigzag edges). In the absence of strain, our results show that (i) the magnetization is affected by taking different dot sizes only for hexagonal BLG QDs with zigzag edges, exhibiting different critical Hubbard interactions, and (ii) the magnetization does not depend on the interlayer hopping energies, except for the geometries with zigzag edges and AA stacking. In the presence of in-plane and uniaxial strain, for all geometries we obtain two different magnetization regimes depending on the applied strain amplitude. The appearance of such different regimes is due to the breaking of layer and sublattice symmetries in BLG QDs.

Abstract: The influence of a perpendicular magnetic field on the binding energy and structural properties of excitons, trions, and biexcitons in monolayers of semiconducting transition metal dichalcogenides (TMDs) is investigated. The stochastic variational method (SVM) with a correlated Gaussian basis is used to calculate the different properties of these few-particle systems. In addition, we present a simplified variational approach which supports the SVM results for excitons as a function of magnetic field. The exciton diamagnetic shift is compared with recent experimental results, and we extend this concept to trions and biexcitons. The effect of a local potential fluctuation, which we model by a circular potential well, on the binding energy of trions and biexcitons is investigated and found to significantly increase the binding of those excitonic complexes.

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: We propose two coupled electron-hole sheets of few-layer graphene as a new nanostructure to observe superfluidity at enhanced densities and enhanced transition temperatures. For ABC stacked few-layer graphene we show that the strongly correlated electron-hole pairing regime is readily accessible experimentally using current technologies. We find for double trilayer and quadlayer graphene sheets spatially separated by a nano-thick hexagonal boron-nitride insulating barrier, that the transition temperature for electron-hole superfluidity can approach temperatures of 40 K.

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.