Efficient skyrmion transport mediated by a voltage controlled magnetic anisotropy gradient.

Despite the inefficiencies associated with current-induced spin torques, they remain the predominant mode of skyrmion propulsion. In this work, we demonstrate numerically that skyrmions can be transported much more efficiently with a voltage-controlled magnetic anisotropy (VCMA) gradient. An analytical model was developed to understand the underlying skyrmion dynamics on a track under the VCMA conditions. Our calculations reveal that the repulsive skyrmion-edge interaction not only prevents the skyrmion from annihilating but also generates most of the skyrmion propulsion. A multiplexed array of gate electrodes can be used to create discrete anisotropy gradients over a long distance, leading to the formation of a series of translatable skyrmion potential wells. Due to the strong confining potentials, skyrmions are transported at a 70% higher packing density. Finally, we demonstrated that this form of skyrmion propulsion can also be implemented on almost any 2D geometry, providing improved versatility over current-induced methods.

Topological state transport in topological insulators under the influence of hexagonal warping and exchange coupling to in-plane magnetizations.

A hexagonal warping term has been proposed recently to explain the experimentally observed 2D equal energy contours of the surface states of the topological insulator Bi2Te3. Differing from the Dirac fermion Hamiltonian, the hexagonal warping term leads to the opening up of a band gap by an in-plane magnetization. We study the transmission between two Bi2Te3 segments subjected to different in-plane magnetizations and potentials. The opening up of a bandgap, and the accompanying displacement and distortion of the constant energy surfaces from their usual circular shapes by the in-plane magnetizations, modify the transverse momentum overlap between the two Bi2Te3 segments, and strongly modulate the transmission profile. The strong dependence of the TI surface state transport of Bi2Te3 on the magnetization orientation of an adjacent ferromagnetic layer may potentially be utilized in, e.g., a memory readout application.

Giant Faraday and Kerr rotation with strained graphene.

Polarized electromagnetic waves passing through (reflected from) a dielectric medium parallel to a magnetic field undergo Faraday (Kerr) rotation of their polarization. Recently, Faraday rotation angles as much as 0.1 rad were observed for terahertz waves propagating through graphene over a SiC substrate. We show that the same effect is observable with the magnetic field replaced by an in-plane strain field which induces a pseudomagnetic field in graphene. With two such sheets a rotation of π/4 can be achieved, which is the required rotation for an optical diode. Similarly a Kerr rotation of 1/4 rad is predicted from a single reflection from a strained graphene sheet.

Charge fractionalization in biased bilayer graphene.

Fractional charge may arise when fermionic zero modes exist in a topological background field. In biased bilayer graphene (BBLG), the bias plays the role of the nontrivial background field. When semi-infinite BBLG with a zigzag edge is used, the dynamics induces an odd number of zero-energy modes, which, together with the conjugation symmetry between positive- and negative-energy states, are the requisite conditions for fractionalization. Exploiting the trigonal interaction to isolate a given zero-energy mode on the zigzag edge, we consider extended and localized modes (the latter being obtained from a localized wavepacket generated by prior irradiation of the sample with an electromagnetic vortex). The valley degeneracy is lifted by a layer asymmetry, while an edge-induced spin polarization breaks the spin degeneracy. We describe scenarios for the detection of charge-1/2 edge states.

The effect of magnetic field and disorders on the electronic transport in graphene nanoribbons.

We developed a unified mesoscopic transport model for graphene nanoribbons, which combines the nonequilibrium Green's function (NEGF) formalism with the real-space π-orbital model. Based on this model, we probe the spatial distribution of electrons under a magnetic field, in order to obtain insights into the various signature Hall effects in disordered armchair graphene nanoribbons (AGNR). In the presence of a uniform perpendicular magnetic field (B[Symbol: see text]-field), a perfect AGNR shows three distinct spatial current profiles at equilibrium, depending on its width. Under nonequilibrium conditions (i.e. in the presence of an applied bias), the net electron flow is restricted to the edges and occurs in opposite directions depending on whether the Fermi level lies within the valence or conduction band. For electrons at an energy level below the conduction window, the B[Symbol: see text]-field gives rise to local electron flux circulation, although the global flux is zero. Our study also reveals the suppression of electron backscattering as a result of the edge transport which is induced by the B[Symbol: see text]-field. This phenomenon can potentially mitigate the undesired effects of disorder, such as bulk and edge vacancies, on the transport properties of AGNR. Lastly, we show that the effect of [Formula: see text]-field on electronic transport is less significant in the multimode compared to the single-mode electron transport.

High and tunable spin current induced by magnetic-electric fields in a single-mode spintronic device.

We proposed that a viable form of spin current transistor is one to be made from a single-mode device which passes electrons through a series of magnetic-electric barriers built into the device. The barriers assume a wavy spatial profile across the conduction path due to the inevitable broadening of the magnetic fields. Field broadening results in a linearly increasing vector potential across the conduction channel, which increases spin polarization. We have identified that the important factors for generating high spin polarization and conductance modulation are the low source-drain bias, the broadened magnetic fields, and the high number of FM gates within a fixed channel length.

Spin polarization of tunneling current in barriers with spin-orbit coupling.

We present a general method for evaluating the maximum transmitted spin polarization and optimal spin axis for an arbitrary spin-orbit coupling (SOC) barrier system, in which the spins lie in the azimuthal plane and finite spin polarization is achieved by wavevector filtering of electrons. Besides momentum filtering, another prerequisite for finite spin polarization is asymmetric occupation or transmission probabilities of the eigenstates of the SOC Hamiltonian. This is achieved most efficiently by resonant tunneling through multiple SOC barriers. We apply our analysis to common SOC mechanisms in semiconductors: pure bulk Dresselhaus SOC, heterostructures with mixed Dresselhaus and Rashba SOC and strain-induced SOC. In particular, we find that the interplay between Dresselhaus and Rashba SOC effects can yield several advantageous features for spin filter and spin injector functions, such as increased robustness to wavevector spread of electrons.

Analytical solution to transport in Brownian ratchets via the gambler's ruin model.

We present an analogy between the classic gambler's ruin problem and the thermally activated dynamics in periodic Brownian ratchets. By considering each periodic unit of the ratchet as a site chain, we calculated the transition probabilities and mean first passage time for transitions between energy minima of adjacent units. We consider the specific case of Brownian ratchets driven by Markov dichotomous noise. The explicit solution for the current is derived for any arbitrary temperature, and is verified numerically by Langevin simulations. The conditions for current reversal in the ratchet are obtained and discussed.

Advanced modeling techniques for micromagnetic systems.

We present a review of micromagnetic and magnetotransport modeling methods which go beyond the standard model. We first give a brief overview of the standard micromagnetic model, which for (i) the steady-state (equilibrium) solution is based on the minimization of the free energy functional, and for (ii) the dynamical solution, relies on the numerical solution of the Landau-Lifshitz-Gilbert (LLG) equation. We present three complements to the standard model, i.e., (i) magnetotransport calculations based on ohmic conduction in the presence of the anisotropic magnetoresistance (AMR) effect, (ii) magnetotransport calculations based on spin-dependent tunneling in the presence of single charge tunneling (Coulomb blockade) effect, and (iii) stochastic micromagnetics, which incorporates the effects of thermal fluctuations via a white-noise thermal field in the LLG equation. All three complements are of practical importance: (i) magnetotransport model either in the ohmic or tunneling transport regimes, enables the conversion of the micromagnetic results to the measurable quantity of magnetoresistance ratio, while (ii) stochastic modeling is essential as the dimensions of the micromagnetic system reduces to the deep submicron regime and approaches the superparamagnetic limit.

Mapping the Monte Carlo scheme to Langevin dynamics: a Fokker-Planck approach.

We propose a general method of using the Fokker-Planck equation (FPE) to link the Monte Carlo (MC) and the Langevin micromagnetic schemes. We derive the drift and diffusion FPE terms corresponding to the MC method and show that it is analytically equivalent to the stochastic Landau-Lifshitz-Gilbert (LLG) equation of Langevin-based micromagnetics. Subsequent results such as the time-quantification factor for the Metropolis MC method can be rigorously derived from this mapping equivalence. The validity of the mapping is shown by the close numerical convergence between the MC method and the LLG equation for the case of a single magnetic particle as well as interacting arrays of particles. We also find that our Metropolis MC method is accurate for a large range of damping factors alpha, unlike previous time-quantified MC methods which break down at low alpha, where precessional motion dominates.