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: 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: Recently, an extremely low lattice thermal conductivity value has been reported for the alkali-based telluride material BaIn2Te4. The value is comparable with low-thermal conductivity metal chalcogenides, and the glass limit is highly intriguing. Therefore, to shed light on this issue, we performed first-principles phonon thermal transport calculations. We predicted highly anisotropic lattice thermal conductivity along different directions via the solution of the linearized phonon Boltzmann transport equation. More importantly, we determined several different factors as the main sources of the predicted ultralow lattice thermal conductivity of this crystal, such as the strong interactions between low-frequency optical phonons and acoustic phonons, small phonon group velocities, and lattice anharmonicity indicated by large negative mode Gruneisen parameters. Along with thermal transport calculations, we also investigated the electronic transport properties by accurately calculating the scattering mechanisms, namely the acoustic deformation potential, ionized impurity, and polar optical scatterings. The inclusion of spin-orbit coupling (SOC) for electronic structure is found to strongly affect the p-type Seebeck coefficients. Finally, we calculated the thermoelectric properties accurately, and the optimal ZT value of p-type doping, which originated from high Seebeck coefficients, was predicted to exceed unity after 700 K and have a direction averaged value of 1.63 (1.76 in the y-direction) at 1000 K around 2 x 1020 cm-3 hole concentration. For n-type doping, a ZT around 3.2 x 1019 cm-3 concentration was predicted to be a direction-averaged value of 1.40 (1.76 in the z-direction) at 1000 K, mostly originating from its high electron mobility. With the experimental evidence of high thermal stability, we showed that the BaIn2Te4 compound has the potential to be a promising mid- to high-temperature thermoelectric material for both p-type and n-type systems with appropriate doping.

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: 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: 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: 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: The effect of the spatially varying spin current on a skyrmion is numerically investigated. It is shown that an inhomogeneous current density induces an elongation of the skyrmion. This elongation can be controlled using current pulses of different strength and duration. Long current pulses lead to a splitting that forms two replicas of the initial skyrmion while for short pulses the elongated skyrmion relaxes back to its initial circular state through rotation in the MHz-GHz frequency range. The frequency is dependent on the strength of the damping coefficient.

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: Recent experimental measurements of light absorption in few-layer black phosphorus (BP) revealed a series of high and sharp peaks, interspersed by pairs of lower and broader features. Here, we propose a theoretical model for these excitonic states in few-layer BP within a continuum approach for the in-plane degrees of freedom and a tight-binding approximation that accounts for interlayer couplings. This yields excitonic transitions between different combinations of the subbands created by the coupled BP layers, which leads to a series of high and low oscillator strength excitonic states, consistent with the experimentally observed bright and dark exciton peaks, respectively. The main characteristics of such subband exciton states, as well as the possibility to control their energies and oscillator strengths via applied electric and magnetic fields, are discussed, towards a full understanding of the excitonic spectrum of few-layer BP and its tunability.

Abstract: Thermal expansion behavior of two-dimensional (2D) nitrides and graphene were studied by ab initio molecular dynamics (MD) simulations as well as quasiharmonic approximation (QHA). Anharmonicity of the acoustic phonon modes are related to the unusual negative thermal expansion (NTE) behavior of the nitrides. Our results also hint that direct ab initio MD simulations are a more elaborate method to investigate thermal expansion behavior of 2D materials than the QHA. Nevertheless, giant NTE coefficients are found for h-GaN and h-AlN within the covered temperature range 100-600 K regardless of the chosen computational method. This unusual NTE of 2D nitrides is reasoned with the out-of-plane oscillations related to the rippling behavior of the monolayers.

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: 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: 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 analyze the phenomenon of back hopping in spin-torque induced switching of the magnetization in perpendicularly magnetized tunnel junctions. The analysis is based on single-shot time-resolved conductance measurements of the pulse-induced back hopping. Studying several material variants reveals that the back hopping is a feature of the nominally fixed system of the tunnel junction. The back hopping is found to proceed by two sequential switching events that lead to a final state P' of conductance close to-but distinct from-that of the conventional parallel state. The P' state does not exist at remanence. It generally relaxes to the conventional antiparallel state if the current is removed. The P' state involves a switching of the sole spin-polarizing part of the fixed layers. The analysis of literature indicates that back hopping occurs only when the spin-polarizing layer is too weakly coupled to the rest of the fixed system, which justifies a posteriori the mitigation strategies of back hopping that were implemented empirically in spin-transfer-torque magnetic random access memories.