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: 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 study the influence of confinement on the atomic collapse due to a Coulomb impurity placed at the center of a graphene quantum dot of radius R. We apply the zigzag or infinite-mass boundary condition and consider both a point-size and a finite-size impurity. As a function of the impurity strength Za, the energy spectra are discrete. In the case of the zigzag boundary condition, the degenerate (with respect to the angular momentum m) zero-energy levels are pulled down in energy as Z alpha increases, and they remain below epsilon = – Z alpha. Our results show that the energy levels exhibit a 1/R dependence in the subcritical regime [Z alpha < |km + 1/2|, k = 1 (-1) for the K (K') valley]. In the supercritical regime (Z alpha > |km + 1/2|) we find a qualitatively very different behavior where the levels decrease as a function of R in a nonmonotonic manner. While the valley symmetry is preserved in the presence of the impurity, we find that the impurity breaks electron-hole symmetry. We further study the energy spectrum of zigzag quantum dots in gapped graphene. Our results show that as the gap increases, the lowest electron states are pushed into the gap by the impurity.

Abstract: We present a polarization-resolved, high-resolution Raman scattering study of the three consecutive charge density wave (CDW) regimes in 1T-TaS2 single crystals, supported by ab initio calculations. Our analysis of the spectra within the low-temperature commensurate (C-CDW) regime shows P (3) over bar symmetry of the system, thus excluding the previously proposed triclinic stacking of the “star-of-David” structure, and promoting trigonal or hexagonal stacking instead. The spectra of the high-temperature incommensurate (IC-CDW) phase directly project the phonon density of states due to the breaking of the translational invariance, supplemented by sizable electron-phonon coupling. Between 200 and 352 K, our Raman spectra show contributions from both the IC-CDW and the C-CDW phases, indicating their coexistence in the so-called nearly commensurate (NC-CDW) phase. The temperature dependence of the symmetry-resolved Raman conductivity indicates the stepwise reduction of the density of states in the CDW phases, followed by a Mott transition within the C-CDW phase. We determine the size of the Mott gap to be Omega(gap) approximate to 170-190 meV, and track its temperature dependence.

Abstract: Edge functionalization of graphene nanoribbons with nitrogen atoms for various adatom configurations at armchair and zigzag edges are investigated. We provide comprehensive information on the electronic and magnetic properties and investigate the stability of the various systems. Two types of rippling of the nanoribbons, namely edge and bulk rippling depending on the sign of edge stress induced at the edge, are found. They are found to play the decisive role for the stability of the structures. We also propose a type of edge decoration in which every third nitrogen adatom at the zigzag edges is replaced by an oxygen atom. In this way, the electron count is compatible with a full aromatic structure, leading to additional stability and a disappearance of magnetism that is usually associated with zigzag nanoribbons.

Abstract: Using the tight-binding model, we investigate the influence of vacancy disorder on electrical transport in graphene Hall bars in the presence of quantizing magnetic fields. Disorder, induced by a random distribution of monovacancies, breaks the graphene sublattice symmetry and creates states localized on the vacancies. These states are observable in the bend resistance, as well as in the total DOS. Their energy is proportional to the square root of the magnetic field, while their localization length is proportional to the cyclotron radius. At the energies of these localized states, the electron current flows around the monovacancies and, as we show, it can follow unexpected paths depending on the particular arrangement of vacancies. We study how these localized states change with the vacancy concentration, and what are the effects of including the next-nearest-neighbor hopping term. Our results are also compared with the situation when double vacancies are present in the system. Double vacancies also induce localized states, but their energy and magnetic field dependencies are different. Their localization energy scales linearly with the magnetic field, and their localization length appears not to depend on the field strength.

Abstract: The optoelectronic properties of confined water form one of the most active research areas in the past few years. Here we present the multiscale methodology to discern the out-of-plane electronic and dipolar dielectric constants (epsilon(el)(perpendicular to) and epsilon(diP)(perpendicular to)) of strongly confined water. We reveal that epsilon(perpendicular to el) and epsilon(diP)(perpendicular to) become comparable for water confined in angstrom-scale channels (with a height of less than 15 angstrom) within graphene (GE) and hexagonal boron nitride (hBN) bilayers. Channel height (h) associated with a minimum in both epsilon(e)(l)(perpendicular to) and epsilon(dip)(perpendicular to) is linked to the formation of the ordered structure of ice for h approximate to (7 -7.5) angstrom. The recently measured total dielectric constant epsilon(T)(perpendicular to) of nanoconfined water [L. Fumagalli et al., Science 360, 1339 (2018)] is corroborated by our results. Furthermore, we evaluate the contribution from the encapsulating membranes to the dielectric properties, as a function of the interlayer spacing, i.e., the height of the confining channel for water. Finally, we conduct analysis of the optical properties of both confined water and GE membranes, and show that the electron energy loss function of confined water strongly differs from that of bulk water.

Abstract: The effect of anisotropy in the energy spectrum on the binding energy and structural properties of excitons, trions, and biexcitons is investigated. To this end we employ the stochastic variational method with a correlated Gaussian basis. We present results for the binding energy of different excitonic complexes in black phosphorus (bP) and TiS3 and compare them with recent results in the literature when available, for which we find good agreement. The binding energies of excitonic complexes in bP are larger than those in TiS3. We calculate the different average interparticle distances in bP and TiS3 and show that excitonic complexes in bP are strongly anisotropic whereas in TiS3 they are almost isotropic, even though the constituent particles have an anisotropic energy spectrum. This is also confirmed by the correlation functions.

Abstract: The electronic and transport properties of a three-dimensional (3D) Dirac system are investigated theoretically, which is motivated by recent experimental measurements on quantum and transport mobilities in the 3D Dirac semimetal Na3Bi by J. Xiong et al. [Science 350, 413 (2015); Europhys. Lett. 114, 27002 (2016)]. The electron Hamiltonian is taken from a simplified k center dot p approach. From the obtained electronic band structure and the Fermi energy, we explain why the anomalous effect induced by the chiral anomaly and the Berry curvature in the energy band can be observed experimentally in magnetotransport coefficients in both low-and high-density samples. Moreover, the quantum and transport mobilities are calculated on the basis of the momentum-balance equation derived from a semiclassical Boltzmann equation with the electron-impurity interaction. The quantum and transport mobilities obtained from this study agree both qualitatively and quantitatively with those measured experimentally. We also examine the electron mobilities along different crystal directions in Na3Bi and find them largely anisotropic. The theoretical findings from this work can be helpful in gaining an in-depth understanding of the experimental results and of the basic electronic and transport properties of newly developed 3D Dirac systems.

Abstract: Using spin- and angle-resolved photoemission spectroscopy and relativistic many-body calculations, we investigate the evolution of the electronic structure of (Bi1-xInx)(2)Se-3)(2)Se-3 bulk single crystals around the critical point of the trivial to topological insulator quantum-phase transition. By increasing x, we observe how a surface gap opens at the Dirac point of the initially gapless topological surface state of Bi2Se3, leading to the existence of massive fermions. The surface gap monotonically increases for a wide range of x values across the topological and trivial sides of the quantum-phase transition. By means of photon-energy-dependent measurements, we demonstrate that the gapped surface state survives the inversion of the bulk bands which occurs at a critical point near x = 0.055. The surface state exhibits a nonzero in-plane spin polarization which decays exponentially with increasing x, and which persists in both the topological and trivial insulator phases. Our calculations reveal qualitative agreement with the experimental results all across the quantum-phase transition upon the systematic variation of the spin-orbit coupling strength. A non-time-reversal symmetry-breaking mechanism of bulk-mediated scattering processes that increase with decreasing spin-orbit coupling strength is proposed as explanation.

Abstract: Considering a noncentrosymmetric pinning texture composed of a square array of triangular holes, the magnetic flux penetration and expulsion are investigated experimentally and theoretically. A direct visualization of the magnetic landscape obtained using a magneto-optical technique on a Nb film is complemented by a multiscale numerical modeling. This combined approach allows the magnetic flux dynamics to be identified from the single flux quantum limit up to the macroscopic electromagnetic response. Within the theoretical framework provided by time-dependent Ginzburg-Landau simulations, an estimation of the in-plane current anisotropy is obtained and its dependence with the radius of the curvature of hole vertices is addressed. These simulations show that current crowding plays an important role in channeling the flux motion, favoring hole-to-hole flux hopping rather than promoting interstitial flux displacement in between the holes. The resulting anisotropy of the critical current density gives rise to a distinct pattern of discontinuity lines for increasing and decreasing applied magnetic fields, in sharp contrast to the invariable patterns reported for centrosymmetric pinning potentials. This observation is partially accounted for by the rectification effect, as demonstrated by finite-element modeling. At low temperatures, where magnetic field penetration is dominated by thermomagnetic instabilities, highly directional magnetic flux avalanches with a fingerlike shape are observed to propagate along the easy axis of the pinning potential. This morphology is reproduced by numerical simulations. Our findings demonstrate that anisotropic pinning landscapes and, in particular, ratchet potentials produce subtle modifications to the critical state field profile that are reflected in the distribution of discontinuity lines.

Abstract: Using a particle-based model we examine the depinning motion of solitons in skyrmion chains in quasi -onedimensional (1D) and two-dimensional (2D) systems containing embedded 1D interfaces. The solitons take the form of a particle or hole in a commensurate chain of skyrmions. Under an applied drive, just above a critical depinning threshold, the soliton moves with a skyrmion Hall angle of zero. For higher drives, the entire chain depins, and in a 2D system we observe that both the solitons and chain move at zero skyrmion Hall angle and then transition to a finite skyrmion Hall angle as the drive increases. In a 2D system with a 1D interface that is at an angle to the driving direction, there can be a reversal of the sign of the skyrmion Hall angle from positive to negative. Our results suggest that solitons in skyrmion systems could be used as information carriers in racetrack geometries that would avoid the drawbacks of finite skyrmion Hall angles. The soliton states become mobile at significantly lower drives than the depinning transition of the skyrmion chains themselves.

Abstract: We investigate a Josephson junction in an electron-hole superfluid in a double-layer transition metal dichalco-genide heterostructure. The observation of a critical tunneling current is a clear signature of superfluidity. In addition, we find the BCS-BEC crossover physics in the narrow barrier region controls the critical current across the entire system. The corresponding critical velocity, which is measurable in this system, has a maximum when the excitations pass from bosonic to fermionic. Remarkably, this occurs for the density at the boundary of the BEC to BCS-BEC crossover regime determined from the condensate fraction. This provides, in a semiconductor system, an experimental way to determine the position of this boundary.

Abstract: Magnetic skyrmions are topological spin configurations in materials with chiral Dzyaloshinskii-Moriya interaction (DMI), that are potentially useful for storing or processing information. To date, DMI has been found in few bulk materials, but can also be induced in atomically thin magnetic films in contact with surfaces with large spin-orbit interactions. Recent experiments have reported that isolated magnetic skyrmions can be stabilized even near room temperature in few-atom-thick magnetic layers sandwiched between materials that provide asymmetric spin-orbit coupling. Here we present the minimum-energy path analysis of three distinct mechanisms for the skyrmion collapse, based on ab initio input and the performed atomic-spin simulations. We focus on the stability of a skyrmion in three atomic layers of Co, either epitaxial on the Pt(111) surface or within a hybrid multilayer where DMI nontrivially varies per monolayer due to competition between different symmetry breaking from two sides of the Co film. In laterally finite systems, their constrained geometry causes poor thermal stability of the skyrmion toward collapse at the boundary, which we show to be resolved by designing the high-DMI structure within an extended film with lower or no DMI.

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.