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: 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: 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 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: 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 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 the Comment, the authors state that the separation of the relative and center of mass variables in our work is not correct. Here we point out that there is a typographical error, i.e., qi instead of -e, in two of our equations which, when corrected, makes the Comment redundant. Within the ansatzes mentioned in our paper all our results are correct, in contrast to the claims of the Comment.