Abstract: Two-dimensional graphene-based van der Waals heterostructures have received considerable interest because of their intriguing characteristics compared with the constituent single-layer two-dimensional materials. Here, we investigate the interfacial characteristics, Schottky contact, and optical performance of graphene/Ga2SSe van der Waals (vdW) heterostructure using first-principles calculations. The effects of stacking patterns, electric gating, and interlayer coupling on the interfacial properties of graphene/Ga2SSe heterostructures are also examined. Our results demonstrate that the Dirac cone of graphene is well preserved at the F point in all stacking patterns due to the weak vdW interactions, which keep the heterostructures feasible such that they can be obtained in further experiments. Moreover, depending on the stacking patterns, a small band gap of about 13-17 meV opens in graphene and has a high carrier mobility, indicating that the graphene/Ga2SSe heterostructures are potential candidates for future high-speed nanoelectronic applications. In the ground state, the graphene/Ga2SSe heterostructures form an n-type Schottky contact. The transformation from an n-type to a p-type Schottky contact or to an Ohmic contact can be forced by electric gating or by varying the interlayer coupling. Our findings could provide physical guidance for designing controllable Schottky nanodevices with high electronic and optical performances.

Abstract: Atomic scale defects in semiconductors enable their technological applications and realization of different quantum states. Using scanning tunneling microscopy and spectroscopy complemented by ab initio calculations we determine the nature of defects in the anisotropic van der Waals layered semiconductor ReS2. We demonstrate the in-plane anisotropy of the lattice by directly visualizing chains of rhenium atoms forming diamond-shaped clusters. Using scanning tunneling spectroscopy we measure the semiconducting gap in the density of states. We reveal the presence of lattice defects and by comparison of their topographic and spectroscopic signatures with ab initio calculations we determine their origin as oxygen atoms absorbed at lattice point defect sites. These results provide an atomic-scale view into the semiconducting transition metal dichalcogenides, paving the way toward understanding and engineering their properties.

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: We calculate the Curie temperature of layered ferromagnets, chromium tri-iodide (CrI3), chromium tri-bromide (CrBr3), chromium germanium tri-telluride (CrGeTe3), and the Ned temperature of a layered antiferromagnet iron di-chloride (FeCl2), using first-principles density functional theory calculations and Monte Carlo simulations. We develop a computational method to model the magnetic interactions in layered magnetic materials and calculate their critical temperature. We provide a unified method to obtain the magnetic exchange parameters (J) for an effective Heisenberg Hamiltonian from first principles, taking into account both the magnetic ansiotropy as well as the out-of-plane interactions. We obtain the magnetic phase change behavior, in particular the critical temperature, from the susceptibility and the specific-heat, calculated using the three-dimensional Monte Carlo (METROPOLIS) algorithm. The calculated Curie temperatures for ferromagnetic materials (CrI3, CrBr3, and CrGeTe3), match well with experimental values. We show that the interlayer interaction in bulk CrI3 with R (3) over bar stacking is significantly stronger than the C2/m stacking, in line with experimental observations. We show that the strong interlayer interaction in R (3) over bar CrI3 results in a competition between the in-plane and the out-of-plane magnetic easy axes. Finally, we calculate the Ned temperature of FeCl2 to be 47 +/- 8 K and show that the magnetic phase transition in FeCl2 occurs in two steps with a high-temperature intralayer ferromagnetic phase transition and a low-temperature interlayer antiferromagnetic phase transition.

Abstract: Atomic collapse in graphene nanoribbons behaves in a fundamentally different way as compared to monolayer graphene due to the presence of multiple energy bands and the effect of edges. For armchair nanoribbons we find that bound states gradually transform into atomic collapse states with increasing impurity charge. This is very different in zigzag nanoribbons where multiple quasi-one-dimensional bound states are found that originates from the zero-energy zigzag edge states. They are a consequence of the flat band and the electron distribution of these bound states exhibits two peaks. The lowest-energy edge state transforms from a bound state into an atomic collapse resonance and shows a distinct relocalization from the edge to the impurity position with increasing impurity charge.

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: The dependence of the Hall and bend resistances on a perpendicular magnetic field and on vacancy defects in a four-terminal phosphorene single layer Hall bar is investigated. A tight-binding model in combination with the Landauer-Buttiker formalism is used to calculate the energy spectrum, the lead-to-lead transmissions, and the Hall and bend resistances of the system. It is shown that the terminals with zigzag edge orientation are responsible for the absence of quantized plateaus in the Hall resistance and peaks in the longitudinal resistance. A negative bend resistance in the ballistic regime is found due to the presence of high- and low-energy transport modes in the armchair and zigzag terminals, respectively. The system density of states, with single vacancy defects, shows that the presence of in-gap states is proportional to the number of vacancies. Quantized plateaus in the Hall resistance are only formed in a sufficiently clean system. The effects of different kinds of vacancies where the plateaus are destroyed and a diffusive regime appears in the bend resistance are investigated.

Abstract: The discovery of the hydrodynamic electron liquid (HEL) in graphene [D. Bandurin et al., Science 351, 1055 (2016) and J. Crossno et al., Science 351, 1058 (2016)] has marked the birth of the solid-state HEL which can be probed near room temperature in a table-top setup. Here we examine the terahertz (THz) magneto-optical (MO) properties of a graphene HEL. Considering the case where the magnetic length l(B) = root h/eB is comparable to the mean-free path l(ee) for electron-electron interaction in graphene, the MO conductivities are obtained by taking a momentum balance equation approach on the basis of the Boltzmann equation. We find that when l(B) similar to l(ee), the viscous effect in a HEL can weaken significantly the THz MO effects such as cyclotron resonance and Faraday rotation. The upper hybrid and cyclotron resonance magnetoplasmon modes omega(+/-) are also obtained through the RPA dielectric function. The magnetoplasmons of graphene HEL at large wave-vector regime are affected by the viscous effect, and results in red-shifts of the magnetoplasmon frequencies. We predict that the viscosity in graphene HEL can affect strongly the magneto-optical and magnetoplasmonic properties, which can be verified experimentally.

Abstract: Using the continuum model, we investigate the electronic properties of two types of bilayer graphene (BLG) quantum ring (QR) geometries: (i) An isolated BLG QR and (ii) a monolayer graphene (MLG) with a QR put on top of an infinite graphene sheet (hybrid BLG QR). Solving the Dirac-Weyl equation in the presence of a perpendicular magnetic field and applying the infinite mass boundary condition at the ring boundaries, we obtain analytical results for the energy levels and corresponding wave spinors for both structures. In the case of isolated BLG QR, we observe a sizable and magnetically tunable band gap which agrees with the tight-binding transport simulations. Our analytical results also show the intervalley symmetry EeK (m) = ???EK??? h (m) between the electron (e) and the hole (h) states (m is the angular momentum quantum number) for the energy spectrum of the isolated BLG QR. The presence of interface boundary in a hybrid BLG QR modifies drastically the energy levels as compared with that of an isolated BLG QR. Its energy levels are tunable from MLG dot to isolated BLG QR and to MLG Landau energy levels as the magnetic field is varied. Our predictions can be verified experimentally using different techniques such as by magnetotransport measurements.

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: A Raman scattering study covering the entire substitution range of the FeSe1-xSx solid solution is presented. Data were taken as a function of sulfur concentration x for 0 <= x <= 1, of temperature and of scattering symmetry. All types of excitations including phonons, spins, and charges are analyzed in detail. It is observed that the energy and width of the iron-related B-1g phonon mode vary continuously across the entire range of sulfur substitution. The A(1g) chalcogenide mode disappears above x = 0.23 and reappears at a much higher energy for x = 0.69. In a similar way the spectral features appearing at finite doping in A(1g) symmetry vary discontinuously. The magnetic excitation centered at approximately 500 cm(-1) disappears above x = 0.23 where the A(1g) lattice excitations exhibit a discontinuous change in energy. The low-energy mode associated with fluctuations displays maximal intensity at the nematostructural transition and thus tracks the phase boundary.

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 quantum anomalous Hall effect was first observed experimentally by doping the Bi2Se3 materials family with chromium, where 5% doping induces an exchange field of around 0.1 eV. In ultrathin films, a topological phase transition from a normal insulator to a Chern insulator can be induced with an exchange field proportional to the hybridization gap. Subsequent transitions to states with higher Chern numbers require an exchange field larger than the (bulk) band gap, but are prohibited in practice by the detrimental effects of higher doping levels. Here, we show that threshold doping for these phase transitions in thin films is controllable by strain. As a consequence, higher Chern states can be reached with experimentally feasible doping, sufficiently dilute for the topological insulator to remain structurally stable. Such a facilitated realization of higher Chern insulators opens prospects for multichannel quantum computing, higher-capacity circuit interconnects, and energy-efficient electronic devices at elevated temperatures.

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: 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.