Transport properties through graphene with sequence of alternative magnetic barriers and wells in the presence of time-periodic scalar potential

We investigate the electronic transport properties of a graphene sheet under the magnetic barriers and wells through the oscillating scalar potential combined with the static scalar potential barrier having two types of uniform and alternative profiles. We compute the total sideband transmission of the system by additional sidebands at energy, in presence of oscillating potential, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$V_1$$\end{document}V1, using the transfer-matrix formalism and the Floquet sidebands series. The oscillating potential, generally, suppresses the Klein tunneling and the confinement of the charge carriers. In the absence of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$V_1$$\end{document}V1, both profiles show the wave vector filtering effect for the carriers by controlling the energy E relative to the potential barrier height, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$V_0$$\end{document}V0. The \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$(N-1)$$\end{document}(N-1)-fold resonance splittings are observed through a region around \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$E=V_0$$\end{document}E=V0 with reduction of the transmission. The transmission vanishes in this region upon increasing the number of magnetic blocks N, strength of the magnetic field B in both configurations. We present an estimate relation for the width of the reduction region expressed in terms of E, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$V_0$$\end{document}V0, B and the angle of incidence of the quasiparticles. We observe, in the second profile, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$(N-1)$$\end{document}(N-1)-fold resonances in the transmission for special values of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$E=V_0$$\end{document}E=V0 with a separation depending on the width of the magnetic blocks. The magnetic field and the width of the magnetic blocks have critical values, where the transmission reduces to zero. All the features observed in the transmission reflect to the conductance. In both configurations, there are some peaks in the conductance corresponding to the resonances of the transmission. The oscillations of the conductance are obtained which was observed in the experimental results. We, also, find the possibility for switching the transport properties of the system by changing the characteristic parameters of the magnetic system.


Model and theoretical method
We consider a monolayer graphene sheet in an external magnetic field subjected to a scalar potential energy. The magnetic field is along the z-direction as B = B(x)e z and has an alternating profile consisting of N magnetic barriers with heights B and widths d B and N magnetic wells with depth B and width d −B so that each barrier is followed by a well, as shown in Fig. 1a. The Landau gauge is chosen by considering A = A(x)e y with B z = ∂ x A(x) , as illustrated in Fig. 1b. We introduce the magnitude of the vector potential, for our model, by where x n = n(d B + d −B ) , running n from 0 to N. Here D = d B − d −B and DB means the net magnetic flux through each block formed by a magnetic barrier and a magnetic well 6 . The total magnetic flux, NDB, is a controlling parameter for the transmission and the conductance. The scalar potential is zero in the non-magnetic regions while it oscillates in time around V 0 (x) with the frequency ω inside the magnetic region, as We consider two types of the static scalar potential profiles, V 0 (x) . The first profile, shown in Fig. 1c, consists of a uniform constant static scalar potential, given by It will be referred to this profile as the uniform static scalar potential and discussed in "Results and discussion" section. In the second configuration, illustrated in Fig. 1d, the static scalar potential changes, alternatively, from zero to V 0 as It will be referred to this profile as the alternative static scalar potential and discussed in "Conclusion" section. The model described above can be realized if we first grow the monolayer graphene on a designed ferromagnetic sample and implant a top gate on it, then we apply the uniform or alternative voltage through the gate 37 . The two band DW Hamiltonian, including the magnetic field and the scalar potential energy is written as x n ≤ x < x n + d B ; (nD + 2d B + x n − x)B, x n + d B ≤ x < x n+1 ; NDB, x ≥ x N , where v F is the Fermi velocity ( ∼ c/300 ), σ = (σ x , σ y ) are the Pauli matrices 1 . Because of the conservation of the y-component of the momentum and the time periodicity of the Hamiltonian in Eq. (5), in accordance with the Floquet theory 25,39,40 , the full wave function can be written as where E F is the Floquet energy and f(t) is a periodic function with period T = 2π/ω . Here ψ 1 (x) and ψ 2 (x) are the spatial parts of the full wave function for each sub-lattice. From now on, we express all the variables in dimensionless form using a typical strength of the magnetic field, B 0 , and the magnetic length scale, l 0 = √ /(eB 0 ) . So, we represent all the lengths, time, ω , k y , A(x) and all the energy scales in the units of l 0 , l 0 /v F , v F /l 0 , l −1 0 , l 0 B 0 and v F /l 0 , respectively. Using these considerations, the DW equation takes the form In the magnetic region, inserting Eqs. (2) and (6) in Eq. (7), the following two coupled differential equations are obtained for ψ 1 (x) and ψ 2 (x) . The variables can be separated by introducing the Floquet sideband energy as where m is an integer called the Floquet sideband index in the magnetic region 26,37 . So, the time dependent part of the wave function is obtained by solving Eq. (9) to get f (t) ∝ exp [−i(E m − E F )t − i(V 1 /ω) sin (ωt)] . The Jacob-Anger identity, exp(iz sin θ) = ∞ n=−∞ J n (z) exp(inθ) , and periodicity of f(t) yield E m = E F + mω . Then, the time dependent part is written as (6) �(x, y, t) = e ik y y e −iE F t/ f (t) ψ 1 (x) ψ 2 (x) , The vector potential profile A(x) corresponding to N magnetic blocks of (a). (c) The uniform static scalar potential barrier profile of the first configuration with width x N and height V 0 that is oscillating sinusoidally with amplitude V 1 and frequency ω . (d) The alternative static scalar potential barriers and wells, both oscillating sinusoidally with amplitude V 1 and frequency ω , corresponding to the second profile. www.nature.com/scientificreports/ and the wave function in the magnetic region is The coupled Eqs. (8) can be rewritten, using Eq. (9), to get the following uncoupled differential equation for the spinor ψ m 1 (x) . Replacing the variable x with q = √ 2/B[k y + A(x)] , Eq. (12) takes the form which has two independent solutions D −a m −1/2 (q) and D −a m −1/2 (−q) , known as the Parabolic Cylinder functions 41 . Here a m = −(E m − V 0 (x)) 2 /(2B) ± 1/2 , where the upper (lower) sign corresponds to the magnetic barrier (well) region. So, the spinors ψ m 1 (x) and ψ m 2 (x) takes the forms in the barriers, while for the wells they are where b m 1 , b m 2 , b m 3 and b m 4 are constant coefficients. In the non-magnetic region, solution of the matrix Eq. (7), for x < 0 , is and for x > x N , where a l 1 , a l 2 , c l 1 and c l 2 are constants and k l = sgn(E l ) E 2 l − [k y + A(x)] 2 . The energy consists of l modes as E l = E 0 + lω , where E 0 is the lowest Floquet energy. For convenience we set E F = E 0 . Now we apply the boundary conditions at the boundaries In the numerical calculation of the infinite summations in Eqs. (11), (16) and (17) the terms are terminated by running l and m from −n to n where n is called the number of Floquet levels. In order to avoid the problem of divergency and get the accurate result, due to the oscillatory nature of the Bessel functions, it is necessary that the argument of the Bessel functions to be small, i.e. n > V 1 /ω . It means that the central time dependent potential region is weakly coupled to the other parts of the system. Using the orthogonality of the set of the oscillating functions exp (ilωt) , dropping the t and y parts, the wave functions can be written as  (16) and (17) are considered as a l 1 = δ l0 (because a single electron income from left side), a l 2 = r l (reflection amplitude), c l 1 = t l (transmission amplitude) and c l 2 = 0. The transfer matrix is written, in the matrix form, as where T ij are square matrices of order 2n + 1 . Then and t l is defined through the Eq. (27). The factor cos θ l / cos ϕ guaranties the flux conservation. Then, the total transmission probability from the set of barriers and wells is Now, knowing the transmission, the zero-temperature conductance can be calculated through the Landauer-Buttiker formalism 42 with G 0 = 2e 2 EL y /(πh) as the conductance unit where L y is the length of the graphene sample along the y-direction. Finally, we summarize the main parameters used in this section in Table 1.

Results and discussion
In this section, the physical results of our numerical analysis on the transmission and conductance, given by Eqs. (29) and (30) are presented. The main interest is to analyse the effects of the incident angle, strength of the magnetic field, widths of the magnetic barriers and wells, energy, static and oscillating scalar potentials on the behavior of the transmission and conductance through the system shown in Fig. 1. The numerical results are generalized to the former and latter profiles in "Uniform static scalar potential" section and "Alternative static scalar potential" section, respectively. The resonance and Klein tunneling effects are studied in "Resonance and Klein" tunneling section and the conductance is investigated in "Conductance" section, for both configurations. The terms of the series in Eq. (29) are kept up to n = 5 . It has been considered B 0 = 0.1 T for the scaling value of the magnetic field. In this case, the length, angular frequency and energy scales are obtained to be l 0 = 81.13 nm, v F /l 0 = 12.325 × 10 12 rad/s and v F /l 0 = 8.113 meV, respectively. All of the coming figures are color online.
Uniform static scalar potential. In this subsection, we concentrate on the graphene-based system with the magnetic profile given by Eq. (1) and the scalar potential profile described by Eqs. (2)   www.nature.com/scientificreports/ V 1 . This limitation can also be viewed for the energy of the quasiparticles as E >| NBD/2 | . So, removing the oscillating scalar potential, turns off the transmission for ϕ < ϕ c or E ≤| NBD/2 | . Then, it is possible to confine the DW quasiparticles by turning off V 1 in the described magnetic structure. A perfect transmission is observed in θ = 30 • , for d B = 1 , however the second lobe is very small. For the rest of this section, we concentrate on the Figure 3a shows the angular dependence of transmission in the, so called, equal-barrier case, i.e. E = V 0 , for one magnetic block with B = 1 and d B = d −B = 1 . The curves are deflected towards negative angles. The transmission lobes become sharper by increasing the energy E. In Fig. 3b, it has been zoomed on the angular dependence of the transmission in the case of E = V 0 = 12 and V 1 = 0 for N = 1, 6 . The increment of the number of magnetic blocks makes the beams shorter and thinner. This induces a strong wave vector filtering and suppression in the transmission for the structure. Figure 4a depicts that for E = 1 and in the absence of the oscillating scalar potential, V 1 = 0 , by applying the static potential V 0 the transmission increases and the Klein tunneling effect turns on. For the special value of V 0 = 12 , a perfect transmission is observed in a wide range of the incident angle. In the absence of V 0 , turning on V 1 , the transmission increases. On the contrary, the transmission decreases when V 1 is turned on, in the presence of V 0 . The Klein tunneling is observed for V 0 = 0 . For E = 5 , some changes appear in the mentioned behaviors (see Fig. 4b). In the absence of the scalar potential ( V 0 = V 1 = 0 ) there is a wide angular range with the perfect transmission. Applying V 0 , V 1 or both of them, the transmission decreases over some parts of this range. The Klein tunneling is observed here, moreover, there are some angles indicating the resonance effects 6 .
The width of the magnetic blocks is one of the important factors in the system which can affect the transmission and the Klein tunneling effects. Figure 5 shows the angular dependence of the transmission for N = 1 , E = 5 ,   Alternative static scalar potential. In this subsection, we consider the graphene-based system with the magnetic profile given by Eq. (1), similar to the case considered in "Results and discussion" section, while the scalar potential profile is described by Eqs. (2) and (4) which is shown in Fig. 1a, b and d. In Fig. 7, the effect of non-zero magnetic flux is shown on the incident angle dependence of the transmission probability for V 1 = 0, 1.98 . The whole parameters are the same as in Fig. 2 in the previous configuration. The transmission is decreased in comparison to the uniform static scalar potential case. There is the same limitation for the incident angle as in the first profile, i.e. ϕ > − sin −1 and E >| NBD/2 | for the energy in absence of V 1 , due to the conservation of k y . Again, turning V 1 on, this confinement is removed. In contrary with the first configuration, the transmission increases by turning V 1 on. From now on, we consider d B = d −B . Corresponding to Fig. 3 in the previous configuration, Fig. 8 is obtained for the angular dependence of transmission at E = V 0 , in this profile. Here, also, increment of the energy makes the lobes sharper in negative angles and applying V 1 makes the lobes wider, which is apparent in Fig. 8a. The lobes become sharper for N > 1 , as it is seen in Fig. 8b. So, a strong wave vector filtering can be achieved by choosing the large values for the energy in the case of multiblock magnetic system. Now we analyze the angular dependence of the transmission with respect to the variation of V 0 and V 1 . The same results shown before in Fig. 4 is, exactly, obtained for V 0 = 0 . Again, there is a small range of the incident angles with low transmission for E = 1 and a wide perfect transmission range for E = 5 and V 1 = 0 . In Fig. 9, the angular dependence of the transmission is depicted varying V 0 (non-zero) and V 1 for E = 5 . The Klein tunneling is observed for all values of V 0 and V 1 . There is a wide angular range with perfect transmission for V 1 = 0 . Turning

Resonance and Klein tunneling.
In this subsection, we study the dependence of the total sideband transmission probability on the parameters of the system as the energy, the scalar potential and strength of the magnetic field. Our aim is investigation of resonances and Klein tunneling effect in the normal incidence for two profiles investigated in "Results and discussion" section and "Conclusion" section. In Fig. 11, the normal transmission is plotted versus the energy for the former (Fig. 11a and c) and latter ( Fig. 11b and d) configurations, respectively. In both profiles, in the absence of V 1 , there is a transmission drop region in the energy interval, given by |E − V 0 (x)| < |k y + B| , inside which the transmission experiences some resonance peaks 34 . This limitation can be expressed in terms of the energy, magnetic field or incident angle by sharing between this condition for different values of vector potential. This is due to the fact that the vector potential, A, takes values between zero and Bd B , for d B = d −B . The eigenstates in the barriers are evanescent and they propagate in the well to form quasibound states, for N ≥ 2 . If the incident energy coincide with this bound state energies in the well, the transmission resonances occurs 43 . So, in this region, the transmission reduces except in resonance peaks. The   www.nature.com/scientificreports/ degenerate eigenlevels in the wells split because of the coupling between the wells via tunneling in the barriers and it leads to the ( N − 1)-fold resonance splittings for N magnetic blocks (see Fig. 11a and b). Far from this region, the eigenstates are propagating states and the DW quasiparticles can transmit perfectly, so the transmission approaches to unity and the Klein tunneling is governed in the magnetic system. There are, also, (N − 1) distinct peaks, associated to the Fabry-Pérot interference, because of interaction between the static and magnetic barriers and wells. The resonance splitting effect exists in the magnetic superlattice (the magnetic system with N ≥ 2 ) versus energy, through an electrostatic barrier with the suppression of Klein tunneling. Applying V 1 Other parameters are the same as in Fig. 3. In the first configuration (uniform static scalar potential) the curve of transmission versus energy exhibits a cusp at E = V 0 for B = 1 and non-zero V 1 . The cusp is removed for V 1 = 0 where the transmission is suppressed (see Fig. 11c). The cusps, induced by the oscillating potential, for more than one magnetic block grow up by the amplitude of the oscillating potential and the number of blocks. They are attributed to the Fabry-Pérot fringes. The transmission around E = V 0 can be turned on (off), for multiple magnetic blocks with N ≥ 4 , just by turning on (off) V 1 . The transmission in E = V 0 decreases by decreasing V 1 and it approaches zero for V 1 = 1.98 and large number of blocks, i.e. N > 40 . By increasing the difference between V 0 and E, the Klein tunneling appears and this occurs in higher differences by increasing V 1 . Increasing B can also turn off the normal transmission (see the inset of Fig. 11a in the absence and the inset of Fig. 11c in the presence of V 1 ). Figure 11b and d show similar effects for the second profile (alternative static scalar potential). Here, unlike the first configuration, the energy E and the static scalar potential V 0 have not the same roles and the transmission shows completely nonsymmetric behavior in two sides of E = V 0 . In the absence of V 1 , the transmission around E = V 0 reduces to zero by increasing the number of blocks ( N ≥ 6 ) or B. There are not any resonance peak for E < V 0 . For V 0 ≥ 5 , a Klein tunneling region is appeared, centered around E = V 0 /2 , which its width increases by increasing V 0 . Here, the barrier and well modes are similar to each other due to the fact |E − V 0 | = E which was, also, seen in the first configuration. The transmission decreases from unity in two areas in both sides of this region, given by |E − V 0 (x)| < |k y + A| . Out of this range, for E > V 0 , the Klein tunneling appears again and also alternative sets of (N − 1)-fold resonances are observed. Turning V 1 on, the transmission around E = V 0 increases and the resonance peaks are changed to the usual peaks. The considerations studied here, can be used in designing switching on or off instruments for the transmission of the charge carriers in the graphene based systems.
The counter plots for the normal transmission for the first configuration are sketched in Fig. 12 versus the energy and the static scalar potential for N = 1, 4 in both cases of absence and presence of the oscillating scalar potential. The energy E and the static scalar potential V 0 have similar roles in this configuration, which is not the case in second profile. In the single block case (Fig. 12a and b), there are ribbon bands in the main diameters, i.e. E = V 0 , which indicates the non-zero minimum for the transmission which its width increases by increasing V 1 . Moving away from the main diameter, the difference between E and V 0 increases and the normal transmission grows towards the unity which leads to the Klein tunneling effect. From the semiclassical point of view 44 , the Dirac fermions subjected to a perpendicular magnetic field are rotating around a circular orbit with the radius of cyclotron radius. Increasing the difference between energy and static scalar potential increases the cyclotron www.nature.com/scientificreports/ radius and the Dirac fermions can easily complete their cyclotron orbits, so the transmission and conductance are increased. In the case of four magnetic blocks ( Fig. 12c and d) the 3-fold resonance lines are appeared close the main diameters. Here, the minimum transmission is zero for multiblock magnetic system, in the absence of V 1 (see Fig. 12c). So, the possibility for switching the total normal transmission to zero is provided by choosing V 0 and E close to each other for N > 3 in the absence of V 1 . The number of Fabry-Pérot fringe patterns (resonances and peaks) increases in the presence of V 1 for the multiple magnetic blocks, symmetrically on both sides of the ribbon band. The line in the center of the main diameter in Fig. 12d corresponds to a cusp for E = V 0 in Fig. 11 which is created and grows by increasing N for B = 1 . The appearance of the second (N − 1)-fold resonance lines are, also, observed in both sides of the main diameter, in Fig. 12d, by applying V 1 = 1.98. In Fig. 13a dependence of the normal transmission to V 0 is investigated in the alternative static scalar potential case. The transmission drop region for V 0 is E − |k y + B| < V 0 < E + |k y + B| which consists of the fluctuations as (N − 1)-fold resonances centered at V 0 = E . In the presence of V 1 the resonance, asymmetrically, alters to the peaks and the fluctuations outside this region are exacerbated. For even number of blocks, the central peak does not include perfect transmission like the others. Increasing N, increases the amplitude of the transmission drop. These behaviors for the transmission are strongly affected by the energy. Depending on the value of the energy, it shows two types of behaviors. In Fig. 13a the energy was E = 5 and, in its inset, the behavior of the transmission is shown for the energies E = 4, 7 with N = 3 and V 1 = 0 . For these energies, two (N − 1)-fold resonances are observed in both sides of V 0 = E and the transmission drops in V 0 = E . In Fig. 13b the behavior of the transmission is studied versus V 0 = E . The transmission, for one magnetic block, increases up to the energies around V 0 = E = 3 and, then, is going to be satisfied in τ = 0.8 . For multiple blocks, the transmission takes oscillatory behavior due to the multireflections in the walls of the magnetic barriers. The period of this oscillations depends on d B = d −B . The (N − 1)-fold resonances are, also, observed as before. So, unlike the first configuration, the transmission is not decreased for V 0 = E but the perfect transmissions are, also, observed in the resonance peaks. Figure 14 shows the total normal transmission of the Dirac fermions versus the magnetic field strength for different values of energy, oscillating scalar potential and number of magnetic blocks in the former and latter profiles. In the former configuration, as is shown in Fig. 14a, the Klein tunneling is observed for weak magnetic fields and it is suppressed in the first critical magnetic field with the magnitude B c1 . The transmission decreases and vanishes around a second critical field B c2 . The radius of the cyclotron orbit decreases with the increment of B and the transmission vanishes 16 . The critical fields B c1 and B c2 depend on E, V 0 , V 1 and N. They both decrease with the increment of the energy, for E < V 0 , due to the spreading of the scalar potential barrier in the whole magnetic region. Applying V 1 , the first (second) critical field decreases (increases) and the difference B c2 − B c1 increases. B c2 decreases by increasing N, as is shown in Fig. 14c. The same analysis is performed for the latter configuration in Fig. 14b and d with the same parameters of the first profile. The first critical magnetic field is almost zero but the second is decreased, in comparison with the former configuration. B c2 increases by the increment of the energy and is independent of V 1 . In both profiles, by adding number of the magnetic blocks, two sets of (N − 1)-fold resonances are appeared; the first between B c1 and B c2 and the second after B c2 . It has been observed that they are transformed, in the presence of V 1 , to the fluctuations with non-perfect transmission. The perfect and zero transmission regions are separated from each other by a drop region which its width decreases by increment of N. For large number of magnetic blocks this region becomes narrower and the zero transmission region is obtained by the condition |E − V (x)| < |k y + B| which yields B > 7 and B > 5 in the first (Fig. 14c) and second (Fig. 14d) configurations, respectively.
The resonance splitting angle spaces are shown in Fig. 15 for different values of energy and number of the magnetic blocks with V 0 = 10 in the former and latter profiles. In the angular interval given by www.nature.com/scientificreports/ Klein scattering induced by the electrostatic barrier. It yields −90 • < ϕ < 53.1 • for E = 5 and −41.8 • < ϕ < 30 • for E = 6 . For the special value of E = V 0 /2 for the energy, the angular profiles of the former and latter profiles coincide (see Fig. 15a), while they are separated for other values of the energy, as shown in Fig. 15b. The (N − 1) -fold resonances are observed out of the angular interval of the propagating state's lobe. It has been observed that, by applying V 1 , some fluctuation peaks are created with non-perfect transmission. The presence and absence of the oscillating scalar potential can be used as a characteristic of the transport properties of the studied magnetic systems. In the whole figures, up to now, the frequency has been set on ω = 2 . Changing ω does not affect the whole discussed physical results. In Fig. 16, the normal transmission has been plotted for ω = 2, 10 in both configurations. The (N − 1)-fold peaks are observed for both frequencies.
To avoid the divergency due to the oscillatory nature of the Bessel functions discussed before, it should be noted that ω > V 1 /n.  www.nature.com/scientificreports/ Conductance. The effects discussed up to now related to the transmission is, also, reflected in the total sideband conductance of the studied magnetic system. In Fig. 17 the total sideband conductance of the system is plotted versus the energy for N = 2 , V 0 = 10 , V 1 = 0 and different value of B for the former (Fig. 17a) and latter ( Fig. 17b) configurations, respectively. In the first profile, the conductance experiences a non-zero minimum in E = V 0 like the normal transmission. This minimum conductance approaches to zero by increasing the number of blocks up to N ∼ 6 in the absence of V 1 while this happens, in the presence of V 1 , for N ∼ 40 . In the second configuration, as is shown in Fig. 17b, the conductance has two minimums around E = 0 and E = V 0 which are attributed to the well and barrier modes and the second one has the same behaviors mentioned for the first configuration. The conductance experiences a resonance peak after E = V 0 . The conductance reaches to a maximum around E = V 0 /2 which decreases in the presence of V 1 and by increasing B. The conductance decreases by increasing B but there is a region around E = V 0 and after E = V 0 in the former and latter configurations, respectively, where the conductance is almost insensitive to the increment of the magnetic field. Applying V 1 , removes the insensitivity of the conductance, as is shown in the insets in Fig. 17. The oscillations of the conductance, particularly shown in Fig. 17a, have already been observed in a similar experimental work 45 . Figure 18 shows the dependence of the conductance to the static scalar potential for the alternative static scalar potential profile. As it is clear in Fig. 18a, the conductance comes to a minimum in V 0 = E except for some special values for the energy with a period depending on d B . For E = 2, 5, 8, 11, · · · , the conductance reaches to a maximum at V 0 = E . Switching V 1 on, the conductance increases around V 0 = E and decreases out of this region. In Fig. 18b the behavior of the conductance is studied for the equal-barrier case, i. e. V 0 = E . For a single block structure, the conductance reaches uniformly to a maximum in low energies and then reduces to zero. For multiple blocks, the conductance has an oscillatory behavior with damping amplitude and a period depending on d B . Increasing N, creates peaks in the maxima which was seen in Fig. 13 for normal transmission. The conductance increases by introducing V 1 for all values of V 0 = E . In the inset of Fig. 18b, the effect of B and www.nature.com/scientificreports/ d B is studied for equal-barrier case. Increasing B does not affect the position of the peaks unless the amplitude of the first maximum is reduced but increasing d B reduces the amplitude and period of the conductance oscillations. In spite of the damping behavior of the transmission and conductance in the equal-barrier case in the first configuration and other perviously studied researches 16 , they have oscillatory behaviors in the second profile of the present work. This increasing behavior can be helpful in controlling and interrupting the transport. Figure 19 illustrates the conductance as a function of the strength of the magnetic field for the first (Fig. 19a) and second (Fig. 19b) configurations, respectively. In accordance with the behavior of normal transmission, studied through Fig. 14, the conductance is suppressed after the second critical magnetic field B c2 for both profiles. Its magnitude is about 1.55 T and 0.69 T by consideration of B 0 = 0.1 T for N = 1 , E = 1 and V 1 = 0 for the former and latter configurations, respectively. It decreases by increment of N, as is shown in the insets of Fig. 19. The critical field B c2 for the latter profile is less than the case for the former configuration with the same parameters. In spite of the fact that the normal transmission is unity for weak magnetic fields ( B < B c1 ) in the first configuration, the conductance decreases by applying V 1 which can be useful in designing of graphene-based nanostructures. The oscillations in the insets are associated with the Fabry-Pérot fringes.

Conclusion
In summary, we have studied the transport properties of a graphene sheet subjected to the magnetic blocks (barriers and wells) with the time-dependent scalar potential barriers. Two configurations of uniform and alternative static scalar potential have been investigated through the magnetic barriers and wells. It has been shown that the non-zero magnetic flux can confine DW quasiparticles in both configurations and the presence of the oscillating scalar potential,V 1 , suppresses this confinement. We concentrated on the states with zero magnetic flux. The wave vector filtering effect has been observed for DW quasiparticles, in the absence of V 1 for the equal-barrier Other parameters are the same as in Fig. 3. www.nature.com/scientificreports/ case (with the energy equal to the static scalar potential, i.e. E = V 0 ) which is enhanced by increasing E and the number of magnetic blocks, N. The transmission is suppressed over increasing the width of the magnetic blocks, in both profiles, and the critical width is associated with the cyclotron orbit diameter for DW fermions in the presence of the magnetic field. The normal transmission and the conductance vanishes, drastically, for both configurations in equal-barrier case by increasing N, B and decreasing V 1 . The resonance effects have been observed specially in the former configuration. The magnetic system in the first configuration becomes full transparent in a wide angular range with the suitable choice for the parameters. The (N − 1)-fold resonance splittings, also, have been shown for normal transmission around the equal barrier energy, which satisfies the non-degenerate eigenlevels and the region arising from the investigated model. So, the transport properties of the magnetic graphene system creates a reduction region by adjusting E and V 0 . Out of this region in the first configuration, the Klein tunneling occurs but in the second one it is observed in a finite region around E = V 0 /2 . The (N − 1)-fold resonances have been observed in the transmission and correspondingly some peaks in the conductance, for special values of E = V 0 , in the second configuration which their distances depends on the width of the magnetic blocks. These (N − 1)-fold resonances are, also, observed by variation of the strength of the magnetic field and the angle of incidence in both configurations. An estimation has been found for the place of this transmission reduction region. It has been shown that it is possible to switch the transport properties of the system by changing the system's characteristic parameters as V 1 , B and N. The phenomena studied here as realistic models are important from experimental point of view and we hope that our results can be useful in designment of graphene based electronic devices.