Abstract: Using first-principle calculations, we systematically study strain tuning of the electronic properties of two- dimensional graphitic carbon nitride nanosheets with empirical formula CnNm. We found the following: (i) the ferromagnetic ordered state in the metal-free systems (n, m) = (4,3), (10,9), and (14,12) remains stable in the presence of strain of about 6%. However, the system (9,7) loses its ferromagnetic ordering when increasing strain. This is due to the presence of topological defects in the (9,7) system, which eliminates the asymmetry between spin up and spin down of the p(z) orbitals when strain is applied. (ii) By applying uniaxial strain, a band gap opens in systems which are initially gapless. (iii) In semiconducting systems which have an initial gap of about 1 eV, the band gap is closed with applying uniaxial strain.

Abstract: The surface termination of MXenes greatly determines the electrochemical properties and ion kinetics on their surfaces. So far, hydroxyl-, oxygen-, and fluorine-terminated MXenes have been widely studied for energy storage applications. Recently, sulfur-functionalized MXene structures, which possess low diffusion barriers, have been proposed as candidate materials to enhance battery performance. We performed first-principles calculations on the structural, stability, electrochemical, and ion dynamic properties of Li-adsorbed sulfur-functionalized groups 3B, 4B, 5B, and 6B transition-metal (M)-based MXenes (i.e., M2CS2 with M = Sc, Ti, Zr, Hf, V, Nb, Ta, Cr, Mo, and W). We performed phonon calculations, which indicated that all of the above M2CS2 MXenes, except for Sc, are dynamically stable at T = 0 K. The ground-state structure of each M2CS2 monolayer depends on the type of M atom. For instance, while sulfur prefers to sit at the FCC site on Ti2CS2, it occupies the HCP site of Cr-based MXene. We determined the Li adsorption configurations at different concentrations using the cluster expansion method. The highest maximum open-circuit voltages were computed for the group 4B element (i.e., Ti, Zr, and Hf)-based M2CS2, which are larger than 2.1 V, while their average voltages are approximately 1 V. The maximum voltage for the group 6B element (i.e., Cr, Mo, W)-based M2CS2 is less than 1 V, and the average voltage is less than 0.71 V. We found that S functionalization is helpful for capacity improvements over the O-terminated MXenes. In this respect, the computed storage gravimetric capacity may reach up to 417.4 mAh/g for Ti2CS2 and 404.5 mAh/g for V2CS2. Ta-, Cr-, Mo-, and W-based M2CS2 MXenes show very low capacities, which are less than 100 mAh/g. The Li surface diffusion energy barriers for all of the considered MXenes are less than 0.22 eV, which is favorable for high charging and discharging rates. Finally, ab initio molecular dynamic simulations performed at 400 K and bond-length analysis with respect to Li concentration verify that selected promising systems are robust against thermally induced perturbations that may induce structural transformations or distortions and undesirable Li release.

Abstract: Flat bands in the energy spectrum have attracted a lot of attention in recent years because of their unique properties and promising applications. Special arrangement of impurities on monolayer graphene are proposed to generate multiflat bands in the electronic band structure. In addition to the single midgap states in the spectrum of graphene with low hydrogen density, we found closely spaced bands around the Fermi level with increasing impurity density, which are similar to discrete lines in the spectrum of quantum dots, as well as the unusual Landau-level energy spectrum of graphene in the presence of a strong magnetic field. The presence of flat bands crucially depends on whether or not there are odd or even electrons of H(F) atoms bound to graphene. Interestingly, we found that a fully hydrogenated (fluoridated) of a hexagon of graphene sheet with six hydrogen (fluorine) atoms sitting on top and bottom in consecutive order exhibits Dirac cones in the electronic band structure with a 20% smaller Fermi velocity as compared to the pristine graphene. Functionalizing graphene introduces various C-C bond lengths resulting in nonuniform strains. Such a nonuniform strain may induce a giant pseudomagnetic field in the system, resulting in quantum Hall effect.

Abstract: The interplay between a magnetic field and the Coulomb potential from a charged vacancy on the electron states in graphene is investigated within the tight-binding model. The Coulomb potential removes locally Landau level degeneracy, while the vacancy introduces a satellite level next to the normal Landau level. These satellite levels are found throughout the positive-energy region, but in the negative-energy region, they turn into atomic collapse resonances. Crossings between Landau levels with different angular quantum number m are found. Unlike the point impurity system in which an anticrossing occurs between Landau levels of the same m, in this work anticrossing is found between the normal Landau level and the vacancy-induced level. The atomic collapse resonance hybridizes with the Landau levels. The charge at which the lowest Landau level m = -1, N = 1 crosses E = 0 increases with enhancing magnetic field. A Landau level scaling anomaly occurs when the charge is larger than the critical charge beta 0.6 and this critical charge is independent of the magnetic field.

Abstract: Transport measurements are readily used to probe different phases in disordered topological insulators (TIs), where determining topological invariants explicitly is challenging. On that note, universal conductance fluctuations (UCF) theory asserts the conductance G for an ensemble has a Gaussian distribution, and that standard deviation 8G depends solely on the symmetries and dimensions of the system. Using a real-space tight -binding Hamiltonian on a system with Anderson disorder, we explore conductance fluctuations in a thin Bi2Se3 film and demonstrate the agreement of their behavior with UCF hypotheses. We further show that magnetic field applied out-of-plane breaks the time -reversal symmetry and transforms the system's Wigner-Dyson class from root symplectic to unitary, increasing 8G by 2. Finally, we reveal that while Bi2Se3 is a strong TI, weak TI and metallic phases can be stabilized in presence of strain and disorder, and detected by monitoring the conductance fluctuations.

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: In a recent paper [Phys. Rev. B 89, 125418 (2014)], the authors argue that it is possible to map the electronic properties of twisted bilayer graphene to those of bilayer graphene in an in-plane magnetic field. However, their description of the low-energy dynamics of twisted bilayer graphene is restricted to the extended zone scheme and therefore neglects the effects of the superperiodic structure. If the energy spectrum is studied in the supercell Brillouin zone, we find that the comparison with an in-plane magnetic field fails because (i) the energy spectra of the two situations exhibit different symmetries and (ii) the low-energy spectra are very different.

Abstract: Often real samples of graphene consist of islands of both monolayer and bilayer graphene. Bound states in such hybrid quantum dots are investigated for (i) a circular single-layer graphene quantum dot surrounded by an infinite bilayer graphene sheet and (ii) a circular bilayer graphene quantum dot surrounded by an infinite single-layer graphene. Using the continuum model and applying zigzag boundary conditions at the single-layer-bilayer graphene interface, we obtain analytical results for the energy levels and the corresponding wave spinors. Their dependence on perpendicular magnetic and electric fields are studied for both types of quantum dots. The energy levels exhibit characteristics of interface states, and we find anticrossings and closing of the energy gap in the presence of a bias potential.

Abstract: The structural and mechanical properties of N-doped graphene (NG) are investigated using reactive force field (ReaxFF) potentials in large-scale molecular dynamics simulations. We found that ripples, which are induced by the dopants, change the roughness of NG, which depends on the number of dopants and their local arrangement. For any doping ratio N/C, the NG becomes ferroelectric with a net dipole moment. The formation energy increases nonlinearly with N/C ratio, while the Young's modulus, tensile strength, and intrinsic strain decrease with the number of dopants. Our results for the structural deformation and the thermoelectricity of the NG sheet are in good agreement with recent experiments and ab initio calculations.

Abstract: Finding novel atomically thin heterostructures and understanding their characteristic properties are critical for developing better nanoscale optoelectronic devices. In this study, we investigate the electronic and optical properties of a GaS-Ca(OH)(2) heterostructure using first-principle calculations. The band gap of the GaS-Ca(OH)(2) heterostructure is significantly reduced when compared to those of the isolated constituent layers. Our calculations showthat the GaS-Ca(OH)(2) heterostructure is a type-II heterojunction which can be used to separate photoinduced charge carriers where electrons are localized in GaS and holes in the Ca(OH)(2) layer. This leads to spatially indirect excitons which are important for solar energy and optoelectronic applications due to their long lifetime. By solving the Bethe-Salpeter equation on top of a single shot GW calculation (G(0)W(0)), the dielectric function and optical oscillator strength of the constituent monolayers and the heterostructure are obtained. The oscillator strength of the optical transition for the GaS monolayer is an order of magnitude larger than the Ca(OH)(2) monolayer. We also found that the calculated optical spectra of different stacking types of the heterostructure show dissimilarities, although their electronic structures are rather similar. This prediction can be used to determine the stacking type of ultrathin heterostructures.

Abstract: The dynamics of driven pancake vortices in layered superconductors is studied using molecular-dynamics simulations. We found that, with increasing driving force, for strong interlayer coupling, the preexisted vortex lines either directly depin or first transform to two-dimensional (2D) pinned states before they are depinned, depending on the pinning strength. In a narrow region of pinning strengths, we found an interesting repinning process, which results in a negative differential resistance. For weak interlayer coupling, individually pinned pancake vortices first form disordered 2D flow and then transform to ordered three-dimensional (3D) flow with increasing driving force. However, for extremely strong pinning, the random pinning-induced thermal-like Langevin forces melt 3D vortex lines, which results in a persistent 2D flow in the fast-sliding regime. In the intermediate regime, the peak effect is found: With increasing driving force, the moving pancake vortices first crystallize to moving 3D vortex lines, and then these 3D vortex lines are melted, leading to the appearance of a reentrant 2D flow state. Our results are summarized in a dynamical phase diagram.

Abstract: Trions and biexcitons in anisotropic two-dimensional materials are investigated within an effective mass theory. Explicit results are obtained for phosphorene and arsenene, materials that share features such as a direct quasiparticle gap and anisotropic conduction and valence bands. Trions are predicted to have remarkably high binding energies and an elongated electron-hole structure with a preference for alignment along the armchair direction, where the effective masses are lower. We find that biexciton binding energies are also notably large, especially for monolayer phosphorene, where they are found to be twice as large as those for typical monolayer transition metal dichalcogenides.

Abstract: We propose a type of edges arising due to the anisotropy inherent in the puckered structure of a honeycomb system such as in phosphorene. Skewed-zigzag and skewed-armchair nanoribbons are semiconducting and metallic, respectively, in contrast to their normal edge counterparts. Their band structures are tunable, and a metal-insulator transition is induced by an electric field. We predict a field-effect transistor based on the edge states in skewed-armchair nanoribbons, where the edge state is gapped by applying arbitrary small electric field E-z. A topological argument is presented, revealing the condition for the emergence of such edge states.

Abstract: We investigate the superconducting states that arise as a consequence of mesoscopic confinement and a multicomponent order parameter in the Ginzburg-Landau model for p-wave superconductivity. Conventional vortices, but also half-quantum vortices and skyrmions, are found as the applied magnetic field and the anisotropy parameters of the Fermi surface are varied. The solutions are well differentiated by a topological charge that for skyrmions is given by the Hopf invariant and for vortices by the circulation of the superconducting velocity. We revealed several unique states combining vortices and skyrmions, their possible reconfiguration with varied magnetic field, as well as temporal and field-induced transitions between vortical and skyrmionic states.

Abstract: Using the continuum model, we investigate the confined states and the corresponding wave functions of ABC-stacked trilayer graphene (TLG) quantum dots (QDs). First, a general infinite-mass boundary condition is derived and applied to calculate the electron and hole energy levels of a circular QD in both the absence and presence of a perpendicular magnetic field. Our analytical results for the energy spectra agree with those obtained by using the tight-binding model, where a TLG QD is surrounded by a staggered potential. Our findings show that (i) the energy spectrum exhibits intervalley symmetry E-K(e)(m) = -E-K'(h)(m) for the electron (e) and hole (h) states, where m is the angular momentum quantum number, (ii) the zero-energy Landau level (LL) is formed by the magnetic states with m <= 0 for both Dirac valleys, that is different from monolayer and bilayer graphene QD with infinite-mass potential in which only one of the cones contributes, and (iii) groups of three quantum Hall edge states in the tight-binding magnetic spectrum approach the zero LL, which results from the layer symmetry in TLG QDs.