• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    Progress on two-dimensional ferrovalley materials

    2024-01-25 07:13:30PingLi李平BangLiu劉邦ShuaiChen陳帥WeiXiZhang張蔚曦andZhiXinGuo郭志新
    Chinese Physics B 2024年1期
    關(guān)鍵詞:李平劉邦

    Ping Li(李平), Bang Liu(劉邦), Shuai Chen(陳帥), Wei-Xi Zhang(張蔚曦), and Zhi-Xin Guo(郭志新),?

    1State Key Laboratory for Mechanical Behavior of Materials,Center for Spintronics and Quantum System,School of Materials Science and Engineering,Xi’an Jiaotong University,Xi’an 710049,China

    2State Key Laboratory for Surface Physics and Department of Physics,Fudan University,Shanghai 200433,China

    3Jinan Aviation Repair Factory,Jinan 250000,China

    4Department of Physics and Electronic Engineering,Tongren University,Tongren 554300,China

    Keywords: ferrovalley,valley polarization,two-dimensional materials,multi-field tunable

    1.Introduction

    For a long time,finding new ways to manipulate electronics to develop new types of electronic components and the next generation of electronics has been a goal that people have been working on.The manipulation of electrons is often inseparable from the intrinsic freedom that it can have.For physics workers,the conceptual discovery and regulation of the unexploited degrees of freedom that electrons possess in solids is a frontier topic.As we all know,electrons have two intrinsic degrees of freedom,spin,and charge.Based on charge and spin degrees of freedom,people have exploited the extensive application of electronic technology and gradually mature spintronics.Do the electrons in the solid still have other degrees of freedom besides spin and charge? The answer is yes.In many crystal structures, such as the metal bismuth, the semiconductor silicon, germanium, gallium arsenide, and the insulating material diamond, there is a new degree of freedom available to the electron called the valley.[1–3]

    In the reciprocal space of Bloch electrons,the valley generally refers to the local energy extremum point in the valence band or conduction band.Although these valleys may degenerate in energy, they are not necessarily equivalent because these valleys tend to be entirely separated in momentum space.For example,the six degenerate but not equivalent valleys (Dirac cones) in graphene are separated in the corner of the first Brillouin region labeledKorK'.The transition of electrons between these discrete valleys needs the involvement of the assistance of impurities or phonons to satisfy the conservation of momentum, which significantly reduces the scattering probability of electronic states between valleys, so that the valley can be seen as the intrinsic discrete degrees of freedom of holes or electrons.The non-equilibrium distribution of electrons in these discrete valleys can be used to realize a new pattern of encoding information.For example,the two logical states of 0 and 1 can correspond to the different occupying states of electrons in these degenerate valleys.The intrinsic freedom of valley can be developed into a potential carrier for information storage and processing in future electronics,related field is called valleytronics.

    In 2004, two-dimensional graphene was experimentally prepared by mechanical stripping.[4]Immediately afterward,the team of ethnic Chinese scientist Qian Niu proposed that the two Dirac cones opened a gap by breaking the spatial inversion symmetry of graphene.It will appear novel quantum transport behavior and valley circularly polarized light dichroism selective adsorption.[5,6]The non-equivalent two valleys(KandK')have the opposite signs of Berry curvatures.Especially, the transportation coefficient and light selection caused by Berry curvature is the opposite,it provides the theoretical foundation for the light field modulation valley polarization.The theoretical scheme has not been experimentally verified,breaking the lattice symmetry of graphene is a huge challenge.Therefore, graphene is not an ideal valley material.Researchers have turned their attention to the 2H phase transition metal chalcogenides (TMDs) represented by MoS2.[7,8]These materials have hexagonal lattices similar to graphene and spatial inversion symmetry breaking, and the valley Hall effect and the optical selection rule of valley resolution can be observed experimentally.[9,10]

    In the 2H-TMDs system, although the unequal valley is separated in momentum space,its energy still remains degenerate and valley polarization cannot be realized.Similar to paraelectric and paramagnetic materials, such valley energy degenerate materials are called paravalley materials.The realization of the degeneracy of the valley is very beneficial to the development of valleytronics.Hence,researchers have gradually developed a series of methods to achieve valley polarization,such as magnetic atom doping,[11–14]magnetic proximity effect,[15–19]external magnetic fields,[20–24]and optical Stark effects.[9,10,25,26]It is worth noting that the valley polarization is achieved by the above means.It is not an intrinsic property of the material,which is not conducive to regulation of the valley degree of freedom.In 2016,Duanet al.predicted the spontaneous valley polarization monolayer 2H-VSe2by the firstprinciples calculations.[27]They named the material ferrovalley material.Until now,a series of ferrovalley materials have been predicted by theory, such as LaBr2,[28]Nb3I8,[29]FeX2(X=Cl,Br),[30,31]VSi2X4(X=N,P,As),[32–35]XY(X=K,Rb, Cs;Y= N, P, As, Sb, Bi),[36]Cr2Se3,[37]MnPTe3,[38]C6N3H3Au,[39]GdI2,[40]BiFeO3,[41]CeI2,[42]and so on.Beyond that, the application of valley freedom cannot be separated from the modulation of the external field.The magnetic coupling of two-dimensional magnetic ferrovalley can effectively tune the reversal of valley polarization,but it is not conducive to the miniaturization and aggregation of valleytronics devices.Therefore, the modulation valley of electrical,stacking, strain, and interface is a new direction to develop valleytronics devices.

    2.Valley freedom and valley polarization

    Two-dimensional materials with hexagonal lattice structures are essential carriers for the investigation of valleytronics.In the section, the basic concepts of valleytronics, valley freedom and valley polarization are introduced through the spatial inversion of broken graphene and 2H-TMDs materials.

    2.1.Valley freedom

    In a two-dimensional system, the concept of valley was first proposed in graphene.Graphene has honeycomb structure that contains two sets of hexagonal sublattices,conventionally known as the A and B sublattices, as shown in Fig.1(a).[43]The two sets of sublattices are connected by spatial inversion symmetry.The corresponding reciprocal lattice is also hexagonal, the orientation is equivalent to the positive lattice being rotated by 90°.The first Brillouin zone is hexagonal, as illustrated in Fig.1(b).In the hexagonal Brillouin zone, two distinct corners (K±=±4π/3ax) are connected by the time reversal symmetry.Other corners of the Brillouin zone are associated withKby translation of lattice vectors.The most striking feature of graphene comes from low energy excitation atKandK'points,where valence and conduction bands touch each other linearly, and these contact points are named Dirac point.

    The semi-classical dynamical equation of Bloch electron under the external electric and magnetic fields can be expressed as[44]

    where?n(k)is used to describe the Berry curvature of Bloch electron.For the system with spatial inversion symmetry,the Berry curvature meets?n(k)=?n(?k);for the system with time reversal symmetry, the Berry curvature meets?n(k)=??n(?k);for the systems protected by both spatial and time reversal symmetry,the Berry curvature is zero everywhere.[45]The Berry curvature is equivalent to the pseudo-magnetic field.When the non-zero Berry curvature exists in the system, the carrier motion will be affected,resulting in the valley Hall effect.Niuet al.found that theK1(K+,K)andK2(K?,K')valleys have opposite non-zero Berry curvature and orbital magnetic moments in graphene with the broken spatial inversion symmetry.[5]As shown in Fig.1(c),the energy extreme point will appear in the six corners of the first Brillouin zone,namely theK1andK2valleys.TheK1andK2valleys are correlated by time reversal symmetry.As shown in Fig.1(d), the carriers ofK1andK2valleys have the same transverse velocity and opposite direction under the applied electric field.They turn around on both sides,respectively.That is the valley Hall effect.However,graphene is protected by both spatial inversion symmetry and time reversal symmetry.Therefore,it isn’t easy to distinguish different valleys and manipulate valley degrees of freedom experimentally,which is greatly limited in practical applications.

    2.2.Valley polarization

    Due to the difficultly distinguishing the different valleys and manipulating valley degrees of freedom in graphene,2HTMDs materials are regarded as more ideal platform for valleytronics investigation.As shown in Figs.2(a)and 2(b),unlike graphene, which only consists of a single planar carbon atom.Monolayer TMDs unit cells are triangular prism composed of two different atoms of transition metal and sulphur,which have naturally broken spatial inversion symmetry.[7]But similar to graphene,monolayer TMDs also have a hexagonal honeycomb structure.Therefore,the valley of degenerate but not equivalent form at theKandK'points due to broken spatial inversion symmetry.

    Fig.2.(a) Coordination environment of Mo (blue ball) in the structure.The golden balls depicted the S atom.[7](b)Top view of the MoS2 monolayer lattice.[7] (c)Valley distribution diagram at the K points.[8] (d)The top valence band is shown blue,and the bottom conduction band is shown pink.The Brillouin zone represents the center hexagon,and the degree of circular polarization is shown by the center hexagon η(k).[7]

    In the 2H TMDs system, although the unequal valley is separated in momentum space,its energy still remains degenerate and valley polarization cannot be realized.Recently,the band splitting is caused by the spin–orbit coupling(SOC)effect of d electron of the heavy metal atom, as shown in Fig.2(c).[8]Caoet al.calculated that the split of the valence band maximum was 0.15 eV.[7,46,47]Moreover,the spin splitting is reversed atKandK'points due to the time reversal symmetry protection.The valley separated in both momentum and energy space is known as valley polarization.It makes it possible to achieve dynamic valley polarization by selective circular dichroism light.Valley polarization can be characterized by the polarization component of photoluminescence(PL)intensity

    As shown in Fig.2(d),the optical selectivity is absolute atK±points,η(K±)=±1.More surprisingly,the selectivity is close to perfect across the entire valley region, with only rapidly changing signs at the valley boundary.It means that the entireK+valley almost completely rejectsσ+photons,whereas the entireK?valley almost completely repelsσ?photons.It indicates that the adoptedσ+(σ?) will produce satisfactory valley polarization.The experimental results show that the valley polarization can not reach 100 % due to the participation of phonons.[9]In addition,the dynamic valley polarization of optical transition usually disappears within a few picoseconds, which greatly limits the practical application of optical induced valley polarization in valleytronics devices.[48]

    3.Intrinsic ferrovalley materials

    Although the degradation of the valley can be achieved by various external means in TMDs materials,the valley polarization cannot be maintained when the external field is removed.Therefore, it is imperative to search for materials with spontaneous valley polarization.In this section, we introduce the investigation progress of ferrovalley materials from forming valley polarization of the orbit.

    3.1.The p-orbital ferrovalley materials

    Firstly,we introduce the p-orbit ferromagnetic ferrovalley material.We predict that theXY(X=K, Rb, Cs;Y=N, P,As,Sb,Bi)monolayer has intrinsic valley polarized quantum anomalous Hall effect.[36]As shown in Figs.3(a) and 3(b),it exhibits a honeycomb lattice with the point groupD3hand space groupPˉ6m2(No.187).The primitive unit cell is composed of anXand aYatom, and its geometry is completely flat as that of graphene.TheYatom has unpaired electrons in its p orbital, leading to appear ferromagnetic.The intrinsic ferromagnetism breaks the time reversal symmetry of the system.As shown in Figs.3(c) and 3(d), the combined interaction of ferromagnetism and SOC causes the monolayerXYto break the energetic degeneracy at theKandK'valleys,inducing spontaneous valley polarization.

    Fig.3.(a)Top and side views of XY monolayer lattice structure.The X and Y atoms are represented by the blue and orange spheres,respectively.(b) The Brillouin zone of the honeycomb lattice.The band structure of CsSb including the SOC for the magnetic moment of Sb along the (c) z and(d)–z axis(out-of-plane),respectively.The Berry curvatures of CsSb(e)over the two-dimensional Brillouin zone and(f)along the high symmetry line.(g) Calculated the AHC σxy as a function of Fermi energy for the CsSb monolayer.The two light blue shades indicate the valley split between K and K' valleys.(h)Schematic diagram of tunable valley polarized quantum anomalous Hall effect in hole doped CsSb monolayer at the K and K' valleys,respectively.The holes are represented by the+symbol.Upward and downward arrows denote the spin up and spin down carriers,respectively.[36]

    To understand the mechanism of the ferrovalley effect inXY, we take CsSb as an example.We try to build an effective Hamiltonian model based on the density functional theory (DFT) calculations.The conduction band minimum(CBM) and valence band maximum (VBM) mainly consist of Sb pxand pyorbitals including the SOC effect.Considering that the little group ofKandK'points belongs toD3hin the out-of-plane magnetic configuration,we employ|ψτv〉=andas the orbital basis for the CBM and VBM,whereτ=±1 shows the valley index referring to theK/K'.Since the CBM and VBM belong to the same spin channel,we take the SOC effect as the perturbation term,which is

    In the out-of-plane magnetic configuration,θ=φ=0, then theterm can be simplified as

    The energy levels of the valleys for the VBM and CBM can be expressed asandrespectively.Then, the valley polarization in the valence and conduction bands can be represented as

    As shown in Figs.3(e) and 3(f), the Berry curvatures of theKandK'points have opposite signs,indicating the typical valley polarization feature.By integrating the Berry curvature over the Brillouin zone, one can further calculate the anomalous Hall conductivity(AHC).As shown in Fig.3(g),a valley polarized Hall conductivity clearly exists in CsSb.Particularly, when the Fermi level lies between the CBM or VBM of theKandK'valleys, the valley polarized Hall conductanceσxycan be obtained.This result confirms the presence of valley polarized quantum anomalous Hall effect(VQAHE)in the CsSb monolayer.Besides, in the hole doping condition, when the magnetic moment direction of CsSb is in +zdirection, the spin-up holes from theK'valley will be generated and accumulated on one boundary of the sample under an in-plane electrical field (upper plane of Fig.3(h)).On the other hand,when the magnetic moment is in the?zdirection,the spin up holes from theKvalley will be generated and accumulated on the opposite boundary of the sample under an in-plane electrical field(lower plane of Fig.3(h)).This characteristic indicates that monolayer CsSb is an ideal candidate for the high-performance valleytronic devices.

    In addition to the p-orbital magnetically induced valley polarization, Jiaet al.found that valley polarization can also be realized by the ferroelectric polarization of p-orbital(monolayer PbO)induced the broken spatial inversion symmetry.[49]As shown in Figs.4(a)and 4(b),the monolayer PbO is orthogonal lattice structure,different from the hexagonal lattice,the unequal valley of monolayer PbO is located on theΓ →XandΓ →Yhigh symmetric paths, respectively.When the system is the para-electric (PE) phase, the valley is energy degenerated.While the system switches to the ferroelectric(FE) phase, the intrinsic FE field breaks the spatial inversion symmetry induced valley degeneracy to disappear.As shown in Figs.4(c)–4(e), the VBM mostly comes from px(aroundVX)and py(aroundVY)orbitals of O, while the CBM mainly contributed px(aroundVX)and py(aroundVY)orbitals of Pb.For the two FE phases, spontaneous polarization actuates the structure into a ferrovalley state: the two pairs of valleys are no longer the same, but move away from the Fermi level in different ways.For theVXandVYvalley of PbO, both spontaneous valley polarization and linear polarization have light dependent optical selectivity.It not only indicates the potential of oxide structures for ultrathin electronic devices,but also provides the platform for future design and realization of twodimensional multiferroics.

    Fig.4.(a) Schematic diagram of monolayer PbO in the non-polar PE phase and FE phase with polarization along the x direction(x-FE)and y direction(y-FE).(b)Brillouin zone of orthogonal lattice.Two valleys are signed VX and VY.(c)–(e)The band structures of monolayer PbO in the x-FE,PE,and y-FE phases,respectively.[49]

    At present, ferrovalley materials of the p-orbital are seldom studied.Usually, the valley splitting of the p-orbital is very significant.Coincidentally, the room-temperature ferromagnetic p-orbital CaCl monolayer has been experimentally synthesized on reduced graphene oxide membranes,[50]which confirms the possibility of synthesizing ferrovalley materials of the p-orbital in experiments.

    3.2.The d-orbital ferrovalley materials

    The concept of ferrovalley was first proposed by Duanet al.in the d-orbital VSe2system.[27]Its structure is the same as that of 2H-MoS2,except for the presence of unpaired electrons in the 3d orbital of the V atom, which makes the monolayer 2H-VSe2ferromagnetic.Ferromagnetism and SOC effect together make monolayer 2H-VSe2decouple the energetically degenerated valleys atK+(K)andK?(K')points,and induce spontaneous valley polarization.The Hamiltonian of monolayer 2H-VSe2can be written as

    whereHsoc(k) andHex(k) represent the SOC effect, and magnetic exchange interaction, respectively.The combination of the SOC effect and ferromagnetic exchange interaction achieves spontaneous valley polarization without external field.The existence of spontaneous valley polarization makes the two unequal valleys of ferrovalley material have different optical band gaps.Optical calculations show that the system has valley-dependent rotation selectivity.As shown in Fig.5(a), when the magnetic moment is along thezaxis,the energetic degeneracy valley breaks.K+valleys can only be excited by left-circularly polarized light, whileK?valley corresponds to right-circularly polarized light,and the absorption peak of right-circularly polarized light shows a little blue shift.As shown in Fig.5(b), while the magnetic moment is along the?zaxis, the optical band gap and the magnitude of valley polarization will be reversed accordingly.Figure 5(c)is a diagram of the interband transition related to the optical band gap.Moreover, the two sets of valleys have opposite and unequal Berry curvatures for monolayer 2H-VSe2, as shown in Fig.5(d).This also makes ferrovalley materials have a measurable Hall voltage, which is similar to the anomalous Hall effect of ferromagnetic materials,called the anomalous valley Hall effect.As shown in Fig.5(e),based on the anomalous valley Hall effect,the magnetically writing and electrically reading memory devices are coming up.The binary information is stored by the valley polarization of the ferrovalley material,which can be modulated by the magnetic moment of the applied magnetic field.

    Moreover,we find that the Janus structure VSiXN4(X=C, Si, Ge, Sn, Pb) can realize multi-valley dependent topological phase transitions under the built-in electric field and strain.[51]As shown in Fig.6(a), with the atomic number ofXincreases, it shows that an interesting topological phase transition from valley semiconductor (VSC) to valley halfsemimetal (VHSM), to valley polarized quantum anomalous Hall insulator(VQAHI),to the VHSM,to the VSC,and to valley metal (VM).Besides, it can be observed that the valence band valley ofK'is obviously lower thanKfor the VSiCN4,VSi2N4, and VSiGeN4, while the conduction band valley ofK'is significantly higher thanKfor VSiSnN4.Through the analysed orbital-resolved band structure, the VBM bands of VSiCN4, VSi2N4, and VSiGeN4are dominated by V dxy/dx2?y2orbitals, while the CBM bands are mainly contributed by dz2orbital of V atom.On the contrary, the VBM bands are mainly dz2orbital, while the CBM are primarily dxy/dx2?y2orbitals for VSiSnN4and VSiPbN4monolayer.More importantly, it can be observed that the atomic number of theXatom induced band inversion between the dz2and dxy/dx2?y2orbitals.Hence,the band structures with the varying strain andXatom are summarized in Figs.6(b)and 6(c).In a word,the valley-related multiple topological phase transitions attribute to the change of the sign of the Berry curvatures atKandK'points.The band gap and topological property with the various strains for the VSiXN4monolayers are summarized in Fig.6(d).

    Fig.5.The band structure with the SOC effect,and the magnetic moments along the(a)z and(b)?z axis,respectively.(c)Schematic diagram of interband transitions dependent on optical band gap.The imaginary parts of the complex dielectric function ε2 for monolayer VSe2. σ+ and σ?represent left and right handed radiation, respectively.(d) Contour maps of Berry curvature in the momentum space for bands are mainly occupied by dxy and dx2?y2 orbits.(e)Diagram of data storage employing hole-doped ferrovalley materials based on anomalous valley Hall effect.[27]

    Through a detailed analysis of the energy band variation with the atomic number of theXelement, it is observed that the topological properties are intimately correlated to the gap at theKandK'points.Figure 6(e)shows the gap at theKandK'points with the increase ofXatomic number,it found that the gap gradually decreases to zero in VSiGeN4.The fitting curve displays that the band gap would reach zero whenXis Ge,which indicates that this could be a phase transition point.It is well known that the built-in electric field was induced when atoms on the upper and lower surfaces are not equal.We define the built-in electric field asEin=(Φ2?Φ1)/Δh,whereΦ1andΦ2represent the electrostatic potential at the bottom and top of VSiXN4, respectively.The Δhrepresents the structural height of VSiXN4.More surprisingly, the variation trend of the built-in electric field is entirely the same as that of the band gap ofKandK', as shown in Fig.6(f).It indicates that the intrinsic of topological phase transition is the strain and built-in electric field induces band inversion between the dz2and dxy/dx2?y2orbitals atKandK'valleys.Our findings not only uncover the mechanism of multiple topological phase transitions,but also provide a new way for the multifield manipulating the spin,valley,and topological physics.In recent years,the investigation of d-orbital ferrovalley materials has developed rapidly,and a series of novel ferrovalley materials have been discovered,includingMBr2(M=Ru,Os),[52]VS2,[53]VSSe,[54–56],VClBr,[57]RuCl2,[58]h-MNX(M=Ti,Zr,Hf;X=Cl,Br)[59]and so on.

    Fig.6.(a)Band structure of the VSiXN4 monolayers including the SOC effect.(b) and (c) Schematic diagrams of valley-dependent topological phase transitions, (b) with the atomic number of X element increasing,and (c) as a function of strain, respectively.(d) Evolution diagram of band structures and Berry curvatures of VSiXN4 with different strains.The valence band,conduction band,and Berry curvatures of VSiXN4 are shown in the green solid line,red solid line,and dotted lines with different strains,respectively.(e)The band gap of VSiXN4 at K and K'points.The dotted line is the fitting K and K'point gap variation trend.(f)The built-in electric field(Ein)as functions of the strain for VSiXN4 monolayers.[50]

    The d-orbital ferrovalley materials have been extensively explored.These systems exhibit extraordinary properties,such as the valley spin Hall effect and valley polarized quantum anomalous Hall effect.The toughest problem is that the room-temperature magnetism cannot be realized in the experiment.Therefore,our unremitting efforts are still needed.

    3.3.The f-orbital ferrovalley materials

    Recently, Shenget al.found the spontaneous valley polarization with f electrons in monolayer CeI2.[42]The monolayer CeI2exhibits a two-dimensional hexagonal lattice with thePˉ6m2 space group.The spatial inversion symmetry is broken in monolayer CeI2.As shown in Fig.7(a),the density of state uncovers that the magnetism originates from the Ce-5d,Ce-4f and I-5p orbital electrons.Figures 7(b) and 7(c) show the band structure of monolayer CeI2including the SOC effect.The addition of the SOC effect eliminates the energy degeneracy of theK+andK?valleys in the top valence band,leading to an intriguing valley polarization.Moreover, it is found that the valley energy difference between theK+andK?valley depends on the direction of magnetization.When the magnetic moment direction is along the+zaxis,the valley polarization of the spin up channel can reach up to 208 meV.Once the magnetization direction is reversed to the?zdirection, the energy of theK?valley is higher than that of theK+valley, exhibiting a valley polarization of?208 meV in the spin down channel.In other words,the valley polarization and spin can be simultaneously flipped by reversing the magnetic moment of Ce atoms, which provides an efficient way to customize the valley characteristics of monolayer CeI2by modulating its direction of magnetization.

    As shown in Fig.7(d),it exhibits the calculated Berry curvature of monolayer CeI2as a contour map of the entire twodimensional Brillouin zone.TheK?andK+points have opposite signs and unequal magnitudes,meaning a robust valleycontrasting Berry curvature.When such p-type CeI2is magnetized along the+zdirection,in the presence of an in-plane electric field,the positive Berry curvature drives the spin down holes from theK+valley to accumulate on the right side of the sample(see Fig.7(e)).Once the magnetization is reversed,the spin-up holes from theK?valley are transferred to the left side of the sample because of the action of negative Berry curvature.In addition to the anomalous valley Hall effect,it is also interesting to note from Fig.7(e)that the charge and spin Hall effects exist in monolayer CeI2,which is much easier to detect in experiments.Such excellent coexistence of charge, valley,and spin Hall current facilitates the integration of electronics,spintronics,and valleytronics.In addition to CeI2,such LaX2(X=Br,I),[28,60]GdXY(X,Y=Cl,Br,I;X/=Y),[61]with felectron were reported.Due to the strong SOC of the f-orbital,rare earth ferrovalley materials of the f-orbital provide another possibility for the experimental studies.

    Fig.7.(a) The orbital-resolved density of state for monolayer CeI2.The SOC band structure with the(b)positive and(c)negative magnetic moment.(d)Berry curvature over the two-dimensional Brillouin zone.(e)Schematic of the anomalous valley Hall effect for the hole doped CeI2 at the K+ and K?valleys under an in-plane electric field.[42]

    4.Multi-field tunable of two-dimensional ferrovalley materials

    The two-dimensional ferrovalley materials maintain spontaneous valley polarization at the monolayer.The discovery of intrinsic ferrovalley materials has greatly promoted the development of valleytronics.The regulation of valley polarization of ferrovalley materials is the basis of information storage in valleytronic devices.The interlayer is bonded by a weak van der Waals interaction,which facilitates to tuning valley polarization by the multi-field(electrical,interlayer stacking, strain, and interface).It can be realized in valleytronics and magnetic memory device applications.

    4.1.Electrical modulation

    The electrically controlled magnet has been widely investigated.[62–65]Recently, the researchers found that the electricity can also modulate valley polarization well.[66–68]Here,we mainly introduce two aspects: ferroelectric substrate and electric field.

    Fig.8.(a)Band structures of the AA'-stacked VSi2N4 bilayers with the interlayer antiferromagnetic coupling.(b) Band structures of the AA'-stacked VSi2N4 bilayers with the interlayer ferromagnetic coupling.The red solid circles and blue hollow circles represent the atomic projected weights of electronic states on the V atoms from the upper and lower monolayer, respectively.(c) and (d) The band and valley splittings as functions of the electric field for the interlayer antiferromagnetic and antiferromagnetic coupling, respectively.The valley splitting is that of the valence bands, while the spin splitting corresponds to the valence band at the K+ valleys.(e) and (f) Schematic depictions of anomalous Hall currents in VSi2N4 bilayers with hole doping.(e)The interlayer antiferromagnetic coupling and(f)the ferromagnetic coupling.[66]

    Firstly, Lianget al.investigated the electric field tune valley degeneracy splitting in VSi2N4bilayer.[66]They found that the valley splitting occurs in the interlayer antiferromagnetism,there are two degenerate valleys of the interlayer ferromagnetism atK±points,as shown in Figs.8(a)and 8(b).For the interlayer antiferromagnetic coupling,the valley splittings do not change with the electric field,because the orbital components of the conduction/valence bands of the two valleys originate from the same monolayer, and the valley splitting of the monolayer is insensitive to the external electric field.Hence, the external electric field results in the spin splittings and has no effect on valley splittings.In addition, the sign of the spin splittings is decided by the field direction.When the field direction changes,the induced spin splittings are accordingly flipped (see Fig.8(c)).In contrast, as illustrated in Fig.8(d), the size of the valley splitting lies on the electric field strength for the interlayer ferromagnetic coupling.While the positive electric field is in the range of 0–0.35 V/nm, the valley splitting is linearly increased, up to 64.1 meV.While the electric field continues enhancing, the valley splitting is unchanged.This trend change is due to the critical field of 0.35 V/nm that leads to a band exchange between the first and the second valence bands at theK?valley.Beyond the critical electric field, the first valence band of two valleys is contributed by the same layer.As a result, the valley-layer coupling is destroyed and the valley splitting is no longer changed with the electric field.In addition, due to the field-induced valley splittings,there is a transition from the direct band gap semiconductor to the indirect band gap semiconductor.Under the action of the vertical electric field, the degeneracy splittings are highly adjustable for both the magnitude and the sign.These intriguing characteristics provide a practical way to design energy-efficient devices based on coupling between multiple electronic degrees of freedom.

    Fig.9.(a) Crystal structures of MnPSe3/Sc2CO2 h-I↓, h-II↓and h-III↓.(b)–(d) Band structure with considering SOC of h-I↓, h-II↓and h-III↓,respectively.(e)–(g)Band structure with considering SOC of h-I↑, h-II↑and h-III↑,respectively.[68]

    Then, Duet al.found that the anomalous valley Hall effect in monolayer MnPSe3can be readily reversed on or off by switching the ferroelectric polarization of monolayer Sc2CO2.[68]As illustrated in Fig.9(a), three typical stacking modes between monolayer MnPSe3and Sc2CO2, where the polarization of Sc2CO2points away from the interface.The band structure including the SOC effect is illustrated in Figs.9(b)–9(d).For the h-I↓and h-II↓configurations,the feature of valley spin splitting of MnPSe3remains unchanged at theKandK'valleys, but the valley spin splitting values atKandK'valleys are no longer degeneracy.The energetic degeneracy of theKandK'valleys in MnPSe3is elevated,resulting in the spontaneous valley polarization.The coexistence of valley polarization and valley spin splitting ensures that the anomalous valley Hall effect is observed in the antiferromagnetic MnPSe3monolayer.In the h-III↓configuration, the nonuniform potentials introduced to Mn1and Mn2atoms by the proximity effect would be rather weak, leading to the small valley spin splitting.When the ferroelectric is switched,all three configurations occur in a semiconductor-tometal transition.As shown in Figs.9(e)–9(g), this transition is also accompanied with the disappearance of the anomalous valley Hall effect and valley physics.Therefore, the anomalous valley Hall effect is ferroelectrically controllable, which is conducive to the development of controllable valleytronic devices.

    4.2.Interlayer stacking mode modulation

    The modulation layer stacking arrangement can realize the control of material properties,such as electronic,magnetic,ferroelectric,and superconductivity.[69–73]Valley polarization also tuned by interlayer stacking mode.[74,75]For the monolayer VSi2P4, it exhibits ferromagnetic with the spontaneous valley polarization.Zhanget al.stacked two monolayers together in the AA mode to construct a bilayer structure.[76]In the bilayer structure,antiferromagnetism usually dominates the interlayer exchange interaction due to competition between the different interlayer orbital hybrids.Combined with the protection ofMzsymmetry, the spontaneous valley polarization is prohibited.For the resulting bilayer structure, a large Berry curvatures with opposite signs would be acquired near theKandK'valleys,respectively.Since theKandK'valleys come from the spin-down channel of the lower layer and the spin-up channel of the upper layer,respectively.

    Under an interlayer translation operation ofor?t3‖), the AA mode could be transformed into an AB (BA)mode whenMzsymmetry is broken.In AB and BA modes,the absence ofMzsymmetry would cause to an out-of-plane electrical polarization.The presence of electric polarization breaks the energetic degeneracy between theKandK'valleys.Therefore,the Berry curvature can be manipulated on a given layer to facilitate the realization of the layer-polarized anomalous Hall effect (LP-AHE).Besides, AB and BA modes are energetically degenerate with opposite electric polarizations,which can be viewed as two FE states.Under the action of interlayer sliding,the two FE states can be reversibly switched,resulting in an interesting sliding ferroelectric state.After the ferroelectric switching,the sign of valley polarization and electric polarization can be reversed at the same time, which has the FE control potential of layer-locked physics.

    Fig.10.Schematic representations of AA',AB'1,AB'2 and AC' modes of a bilayer structure and their low-energy band dispersions.Band structures for(b)AB'1 and(c)AB'2 configuration of bilayer MnBi2Te4 with SOC,respectively.(d)Berry curvatures along high-symmetry points for AB'1 and AB'2 configurations of bilayer MnBi2Te4.(e) The AHC for AB'1 and AB'2 configurations of bilayer MnBi2Te4.(f) Schematic diagrams of the LP-AHE for bilayer MnBi2Te4.[75]

    The low-energy band dispersions near theKandK'valleys for AB and BA modes are acquired from the k·p model.Evidently, the spontaneous valley polarization is realized in AB and BA modes.For the AB(BA)mode, the VBM of theK(K')valley comes from the spin-up(spin-down)channel of the upper (lower) layer, while the CBM of theK'(K) valley comes from the spin-down(spin-up)channel of the lower(upper) layer.When the Fermi level shifts between theKandK'valleys in an in-plane external electric field,the spin-down holes from theKvalley can transversely move to the right edge of the upper layer,resulting in the LP-AHE in AB-stacked bilayer VSi2P4.Unlike the AB mode,the holes in theK'valley for the BA pattern would acquire an inversed anomalous velocity of–υa.While the Fermi level moves between theKandK'valleys,the spin-up holes from theK'valley would be accumulated at the left edge of the lower layer under an in-plane electric field, forming the LP-AHE.This case is also applied to the valleys of the conduction band.

    Moreover, Penget al.also found the LP-AHE in bilayer MnBi2Te4.[72]As shown in Fig.10(a),two configurations can be regarded as non-polar intermediate states, corresponding to two ferroelectric transition paths between AB'1 and AB'2.Figures 10(b)and 10(c)exhibit the band structures,the bands from the bottom and top layers are separated for the AB'1 and AB'2 configurations.In addition, the calculated Berry curvature and the AHC are illustrated in Figs.10(d) and 10(e), it can confirm the LP-AHE in bilayer MnBi2Te4.By shifting the Fermi level, the electrons would be obtained at the right edge of the top layer under an in-plane electric field,as shown in Fig.10(f).The new phenomena and mechanism provide an important new way to realize the LP-AHE and explore its application in electronics.

    4.3.Strain modulation

    Some physical parameters of two-dimensional ferrovalley, such as valley polarization, orbit arrangement and so on, are very sensitive to lattice deformation caused by strain.It is well known that the ferrovalley property mainly originated from the broken symmetry, strain will aggravate symmetry breaking.Huanet al.found that the strain can induce half-valley metals and topological phase transitions inMBr2(M=Ru, Os) monolayer.[52]As shown in Figs.11(a)–11(e),the band structure of RuBr2can be effectively tuned between ferrovalley insulator and half-valley metal many times with the compressive strain.More interestingly,when 1.6%<η <2.8%, the material turns into a Chern insulating state from a ferrovalley insulator.The nature of valley related topological phase transition is the change of band gap at theK+andK?points under the compressive strain.It is summarized in Fig.11(f).The phase transition from a ferrovalley insulator to a quantum anomalous Hall insulator and then to a ferrovalley insulator again can be obviously found.

    In particular,the RuBr2monolayer becomes a half-valley metal with the compressive strainη= 1.6% or 2.8%.As shown in Figs.11(g) and 11(h), a clear edge state connects the bulk valence and conduction bands for the 2% strain.When the magnetization is switched,the calculated edge state changes from theK?toK+valley,as an opposite chirality and spin direction.It means that the chiral-spin–valley was locked for the acquired edge state.

    Their rich valley and topological phase transition properties provide broad prospects for the application of electronic device.As shown in Fig.11(i),without the strain,the system can exhibit anomalous charge/spin/valley Hall effect.While the small compressive strain is employed, 100% spin- and valley-polarization devices can be established based on the half-valley metal states acquired.On this basis,a novel Chern insulator with a ferrovalley edge state is obtained.When larger compressive strain is used,the material again transforms into a half-valley metal but conducted by another valley.It shows that strain engineering can be used as an effective method to modulate the valley and topological properties of the new ferrovalley materials.

    Fig.11.(a)–(e) Band structures of the monolayer RuBr2 including the SOC effect under the compressive strain.The magenta curves represent the Berry curvatures.(f)The phase diagram of the monolayer RuBr2 as a function of compressive strain.(g)and(h)The calculated edge states for a semifinite RuBr2 monolayer at the 2%compressive strain.The magnetic moment along z and ?z,respectively.(i)Schematic representations of the electronic device prototypes made of the monolayer RuBr2 modulated by the compressive strain.[52]

    4.4.Magnetic interface modulation

    The abundant physical phenomena are derived at the interface of two materials with different properties.The interface neighbor effect is an effective way to modulate the valley properties of ferrovalley materials.Maet al.constructed the 2H-VS2/Cr2C heterostructures, as shown in Figs.12(a)–12(d).[77]They investigated the effects of four different magnetic configurations on valley polarization.In the four configurations, all the heterojunctions are semiconducting.In Figs.12(a) and 12(b), the VBM band is the spin-down channel.As shown in Figs.12(c)and 12(d),while the magnetization is switched,the VBM band changes to the spin-up channel.Moreover,while the direction of magnetization is out-ofplane, as illustrated in Figs.12(a) and 12(c), the valley split has the opposite sign but the same value because of the opposite direction of magnetization.More interestingly, while the direction of magnetization turns into in-plane, as shown in Figs.12(b) and 12(d), the valley polarization ofKandK'points disappears.In brief,the switch direction of magnetization causes the valley polarization and spin in reverse.

    Fig.12.(a)–(d) Four different magnetization configurations of 2HVS2/Cr2C and its energy structure.The arrows represent the direction of magnetic moments,and the red dotted line box exhibits the atomic structure of heterjunctions.The circles present the e2 (dxy+dx2?y2)orbital of VS2,the color scale represents the spin projection.[77]

    5.Summary and outlook

    We introduce the investigation progress of twodimensional ferrovalley materials.Firstly, we introduce the background of research related to the valley.Then, the basic concepts of valleytronics,valley freedom and valley polarization are introduced through the spatial inversion of broken graphene and 2H-TMDs materials.Following,we mainly discuss the ferrovalley materials formed by different orbits.Finally, the multi-field tunable of two-dimensional ferrovalley materials is concerned.

    Even though two-dimensional ferrovalley materials have been predicted early by the DFT calculation,the experimental progress has been limited, partly related to the lack of powerful sample characterization and preparation technique.Recently,Guanet al.first reported the successful preparation of intrinsic bulk ferrovalley material in the experiment.[78]For two-dimensional ferrovalley,so far,there has been no experimental progress.

    The key to electronic device design is to construct a binary logic switching state.Similar to spin,valley is a new type of freedom,different valley polarization states can be used to encode and store information.At present,spintronics devices are mainly based on ferromagnetic materials.The investigation of electronic devices based on spin injection semiconductor materials is still in its infancy.Ferrovalley has strong locking between valley and spin, it provides a new platform for the study of spintronics.Moreover,a large number of valleybased spintronic devices have been explored.For example,using valley to transmit information has high fidelity, direction and room temperature operability.[79]More importantly,the neural network calculation using valley transistor has high accuracy.[80]The spin based on the valley is more stable,and the spin control and detection can be realized with the help of valley,which expands the application prospect of the device.

    Acknowledgments

    Project supported by the National Natural Science Foundation of China (Grant Nos.12074301 and 12004295),China’s Postdoctoral Science Foundation funded project(Grant No.2022M722547), the Open Project of State Key Laboratory of Surface Physics (Grant No.KF202209), and the Natural Science Foundation of Guizhou Provincial Education Department(Grant No.ZK[2021]034).

    猜你喜歡
    李平劉邦
    李平,漫游在文壇的夢境
    金橋(2020年11期)2020-12-14 07:52:52
    四季的美
    第三個小板凳
    THE OSCILLATION OF THE POISSON SEMIGROUP ASSOCIATED TO PARABOLIC HERMITE OPERATOR?
    漢初劉邦集團“入都關(guān)中”事跡查驗
    聯(lián)墨雙馨
    對聯(lián)(2017年2期)2017-05-22 06:12:41
    劉邦記嫂之恨
    劉邦斬丁公
    李平藝術(shù)作品欣賞
    劉邦不識字嗎?
    中文字幕av电影在线播放| 亚洲国产欧美日韩在线播放| 日日摸夜夜添夜夜添小说| av免费在线观看网站| 亚洲精品一区蜜桃| av天堂在线播放| av天堂在线播放| 亚洲 国产 在线| 纯流量卡能插随身wifi吗| 男人爽女人下面视频在线观看| 两性午夜刺激爽爽歪歪视频在线观看 | 久久久国产欧美日韩av| 啦啦啦 在线观看视频| 中亚洲国语对白在线视频| 久久久久国产一级毛片高清牌| 亚洲情色 制服丝袜| 国产精品久久久久成人av| 99久久人妻综合| 十八禁网站网址无遮挡| 人人妻人人澡人人看| av国产精品久久久久影院| 动漫黄色视频在线观看| 热99久久久久精品小说推荐| 国产精品久久久久久人妻精品电影 | 捣出白浆h1v1| 精品亚洲乱码少妇综合久久| 久久国产精品影院| 超碰成人久久| 桃花免费在线播放| 一区在线观看完整版| 777米奇影视久久| 亚洲专区字幕在线| 久久狼人影院| 电影成人av| 成人三级做爰电影| 欧美成人午夜精品| 十八禁网站免费在线| 一区二区三区精品91| 成人影院久久| 飞空精品影院首页| 日韩三级视频一区二区三区| 久久久国产成人免费| 高潮久久久久久久久久久不卡| 亚洲五月色婷婷综合| 午夜免费观看性视频| 欧美精品一区二区大全| 搡老乐熟女国产| 亚洲精品成人av观看孕妇| 考比视频在线观看| 国产色视频综合| 国产精品一区二区在线不卡| 每晚都被弄得嗷嗷叫到高潮| 亚洲视频免费观看视频| 国产精品影院久久| 欧美日韩亚洲高清精品| 老汉色∧v一级毛片| 极品少妇高潮喷水抽搐| 国产主播在线观看一区二区| 中国美女看黄片| 国产成人啪精品午夜网站| 丝袜脚勾引网站| avwww免费| 欧美精品一区二区大全| 亚洲情色 制服丝袜| 丰满饥渴人妻一区二区三| 久久av网站| 成人手机av| 久久综合国产亚洲精品| 欧美日韩亚洲高清精品| 香蕉丝袜av| 十分钟在线观看高清视频www| 9191精品国产免费久久| 少妇精品久久久久久久| 老鸭窝网址在线观看| √禁漫天堂资源中文www| 97精品久久久久久久久久精品| 国产精品久久久人人做人人爽| 男女边摸边吃奶| 正在播放国产对白刺激| 黄网站色视频无遮挡免费观看| 亚洲精品一二三| 久久久久久久久久久久大奶| 成人影院久久| 黑人巨大精品欧美一区二区蜜桃| 国产精品自产拍在线观看55亚洲 | 日韩中文字幕欧美一区二区| 丝袜在线中文字幕| videos熟女内射| 免费观看人在逋| 精品欧美一区二区三区在线| 一边摸一边做爽爽视频免费| 美女高潮到喷水免费观看| 精品久久久久久电影网| 啦啦啦视频在线资源免费观看| 久久人人97超碰香蕉20202| 一个人免费在线观看的高清视频 | 飞空精品影院首页| 精品免费久久久久久久清纯 | 国产97色在线日韩免费| 水蜜桃什么品种好| 大香蕉久久网| 桃花免费在线播放| 人妻人人澡人人爽人人| 国产在线观看jvid| a级毛片在线看网站| 国产人伦9x9x在线观看| 欧美+亚洲+日韩+国产| 久久国产精品大桥未久av| 一本—道久久a久久精品蜜桃钙片| 久久精品国产综合久久久| 国产亚洲精品第一综合不卡| 啪啪无遮挡十八禁网站| 窝窝影院91人妻| 99国产极品粉嫩在线观看| 免费观看人在逋| 国产男女内射视频| 亚洲精品在线美女| 中文字幕人妻丝袜制服| 伦理电影免费视频| 亚洲 国产 在线| 国内毛片毛片毛片毛片毛片| 久久久精品94久久精品| 亚洲三区欧美一区| 亚洲av日韩精品久久久久久密| 菩萨蛮人人尽说江南好唐韦庄| 国产精品久久久久成人av| 久久精品aⅴ一区二区三区四区| 日韩大码丰满熟妇| 9191精品国产免费久久| 1024香蕉在线观看| 日韩大片免费观看网站| 韩国高清视频一区二区三区| 成人亚洲精品一区在线观看| 2018国产大陆天天弄谢| 妹子高潮喷水视频| 肉色欧美久久久久久久蜜桃| 人人澡人人妻人| 久久久久久久国产电影| cao死你这个sao货| 青青草视频在线视频观看| 欧美+亚洲+日韩+国产| 国产免费福利视频在线观看| 黑人猛操日本美女一级片| 亚洲精品中文字幕在线视频| 国产精品av久久久久免费| 久久久国产一区二区| 建设人人有责人人尽责人人享有的| 啦啦啦 在线观看视频| 国产区一区二久久| 国产精品免费视频内射| 亚洲精品国产精品久久久不卡| 十分钟在线观看高清视频www| 欧美精品亚洲一区二区| 日韩,欧美,国产一区二区三区| 女人精品久久久久毛片| 亚洲精品国产区一区二| 国产福利在线免费观看视频| 久久 成人 亚洲| 无遮挡黄片免费观看| 男女边摸边吃奶| 啦啦啦视频在线资源免费观看| 欧美av亚洲av综合av国产av| 母亲3免费完整高清在线观看| 天堂8中文在线网| 精品人妻一区二区三区麻豆| 免费人妻精品一区二区三区视频| 无限看片的www在线观看| 另类亚洲欧美激情| 日本撒尿小便嘘嘘汇集6| 一边摸一边抽搐一进一出视频| 国产又色又爽无遮挡免| 精品久久久久久电影网| 亚洲精品在线美女| 欧美在线一区亚洲| 免费一级毛片在线播放高清视频 | 麻豆av在线久日| 日日夜夜操网爽| 巨乳人妻的诱惑在线观看| a在线观看视频网站| 国产成人精品久久二区二区免费| kizo精华| 中文字幕色久视频| 欧美 日韩 精品 国产| 精品免费久久久久久久清纯 | 啪啪无遮挡十八禁网站| 一本色道久久久久久精品综合| 欧美大码av| 97精品久久久久久久久久精品| 国产淫语在线视频| 狠狠狠狠99中文字幕| 老司机午夜十八禁免费视频| 午夜免费鲁丝| 激情视频va一区二区三区| 精品国产一区二区三区四区第35| 国产欧美亚洲国产| 如日韩欧美国产精品一区二区三区| 一本久久精品| 男女床上黄色一级片免费看| 成人三级做爰电影| 亚洲少妇的诱惑av| 麻豆国产av国片精品| 99久久99久久久精品蜜桃| 80岁老熟妇乱子伦牲交| 成人国语在线视频| 国产欧美日韩一区二区三区在线| 国产成人影院久久av| 99九九在线精品视频| 两个人免费观看高清视频| 亚洲人成电影免费在线| av免费在线观看网站| 国产亚洲欧美在线一区二区| 中文字幕最新亚洲高清| 欧美在线黄色| 日韩制服骚丝袜av| 国产一级毛片在线| 电影成人av| 亚洲国产欧美网| 国产一卡二卡三卡精品| 老司机福利观看| 亚洲精品中文字幕在线视频| 视频在线观看一区二区三区| 欧美精品啪啪一区二区三区 | 黑人猛操日本美女一级片| 伊人久久大香线蕉亚洲五| av超薄肉色丝袜交足视频| av有码第一页| 国产精品一区二区在线不卡| 日韩精品免费视频一区二区三区| 三级毛片av免费| 久久国产精品影院| 中文字幕另类日韩欧美亚洲嫩草| 91精品伊人久久大香线蕉| 80岁老熟妇乱子伦牲交| 大码成人一级视频| 18禁观看日本| 久久久久国内视频| 曰老女人黄片| 国产一区二区三区在线臀色熟女 | 国产成人精品在线电影| 脱女人内裤的视频| 国产成人系列免费观看| 久久精品久久久久久噜噜老黄| 欧美日韩精品网址| 亚洲七黄色美女视频| 脱女人内裤的视频| 深夜精品福利| 一级片'在线观看视频| 蜜桃国产av成人99| 亚洲国产欧美一区二区综合| 老司机深夜福利视频在线观看 | 又黄又粗又硬又大视频| 国产欧美亚洲国产| 国产精品一区二区免费欧美 | 91老司机精品| 国产一卡二卡三卡精品| 欧美国产精品一级二级三级| 国产在线观看jvid| 最新在线观看一区二区三区| 三上悠亚av全集在线观看| 亚洲国产欧美一区二区综合| 欧美性长视频在线观看| 亚洲九九香蕉| 99九九在线精品视频| 国产精品二区激情视频| 亚洲精品国产色婷婷电影| 精品人妻一区二区三区麻豆| 欧美成狂野欧美在线观看| 老司机影院成人| 女性生殖器流出的白浆| 国产99久久九九免费精品| 久久毛片免费看一区二区三区| 亚洲av男天堂| 精品亚洲乱码少妇综合久久| 午夜老司机福利片| 在线亚洲精品国产二区图片欧美| 国产一区二区三区av在线| 9热在线视频观看99| 99精国产麻豆久久婷婷| 精品高清国产在线一区| 五月开心婷婷网| 久久ye,这里只有精品| 午夜影院在线不卡| 久久久久视频综合| 国产野战对白在线观看| 欧美日韩一级在线毛片| 黄色 视频免费看| 亚洲精品自拍成人| 少妇 在线观看| 欧美日韩亚洲综合一区二区三区_| 午夜激情久久久久久久| 国产区一区二久久| 性色av乱码一区二区三区2| 中文字幕最新亚洲高清| 久久久国产欧美日韩av| 如日韩欧美国产精品一区二区三区| 精品国产国语对白av| 考比视频在线观看| 亚洲免费av在线视频| 激情视频va一区二区三区| 中文字幕人妻丝袜制服| 天堂8中文在线网| 亚洲欧美色中文字幕在线| 久久性视频一级片| 久久精品亚洲熟妇少妇任你| 日韩欧美一区视频在线观看| 别揉我奶头~嗯~啊~动态视频 | 99精国产麻豆久久婷婷| 国产精品欧美亚洲77777| 麻豆乱淫一区二区| 久久久水蜜桃国产精品网| 一本色道久久久久久精品综合| 99国产精品一区二区三区| 黄片小视频在线播放| 免费在线观看黄色视频的| 少妇 在线观看| 国产成人精品在线电影| 成年女人毛片免费观看观看9 | 精品第一国产精品| 久久av网站| 欧美精品av麻豆av| 国产91精品成人一区二区三区 | 麻豆乱淫一区二区| 桃红色精品国产亚洲av| 亚洲精品第二区| 69av精品久久久久久 | 精品视频人人做人人爽| 日韩中文字幕视频在线看片| 日本欧美视频一区| 亚洲国产av影院在线观看| 久久国产精品影院| 久久女婷五月综合色啪小说| www.熟女人妻精品国产| 久久精品国产亚洲av香蕉五月 | 操美女的视频在线观看| 伊人久久大香线蕉亚洲五| 99香蕉大伊视频| 国产亚洲一区二区精品| 欧美亚洲 丝袜 人妻 在线| 国产精品99久久99久久久不卡| 亚洲欧美精品自产自拍| 国产精品久久久久成人av| 美女脱内裤让男人舔精品视频| 热re99久久国产66热| 国产亚洲欧美在线一区二区| 亚洲av片天天在线观看| 国内毛片毛片毛片毛片毛片| 免费在线观看视频国产中文字幕亚洲 | 极品人妻少妇av视频| 最新在线观看一区二区三区| bbb黄色大片| 久久青草综合色| 久久亚洲国产成人精品v| a级毛片在线看网站| 亚洲 欧美一区二区三区| 国产精品麻豆人妻色哟哟久久| 啪啪无遮挡十八禁网站| 亚洲情色 制服丝袜| 天堂俺去俺来也www色官网| 久久久久精品人妻al黑| 亚洲免费av在线视频| 免费日韩欧美在线观看| 国产精品国产av在线观看| 中文字幕人妻熟女乱码| av电影中文网址| 亚洲男人天堂网一区| 亚洲国产av新网站| 欧美国产精品va在线观看不卡| 新久久久久国产一级毛片| 91av网站免费观看| 欧美精品av麻豆av| 国产亚洲一区二区精品| 国产一区二区激情短视频 | 少妇 在线观看| 国产精品av久久久久免费| 久久久久久久久免费视频了| 久久久久精品人妻al黑| 中文字幕高清在线视频| 日韩欧美免费精品| 亚洲精品成人av观看孕妇| 男女免费视频国产| 一级a爱视频在线免费观看| 中亚洲国语对白在线视频| 又紧又爽又黄一区二区| 久久九九热精品免费| 亚洲精品在线美女| 女警被强在线播放| 日本五十路高清| 99久久精品国产亚洲精品| 婷婷色av中文字幕| 日韩人妻精品一区2区三区| xxxhd国产人妻xxx| 午夜免费成人在线视频| 中文字幕另类日韩欧美亚洲嫩草| 啦啦啦视频在线资源免费观看| av在线app专区| 久久性视频一级片| 在线十欧美十亚洲十日本专区| 国产三级黄色录像| 国产一区二区 视频在线| 高清视频免费观看一区二区| 日日摸夜夜添夜夜添小说| 午夜福利免费观看在线| 亚洲 欧美一区二区三区| 亚洲成人国产一区在线观看| 成人国产av品久久久| 久9热在线精品视频| 制服人妻中文乱码| 无遮挡黄片免费观看| 国产熟女午夜一区二区三区| 黑人巨大精品欧美一区二区mp4| 久久精品成人免费网站| 妹子高潮喷水视频| 精品一区二区三区四区五区乱码| 在线观看一区二区三区激情| 丰满人妻熟妇乱又伦精品不卡| 国产精品久久久人人做人人爽| 久久影院123| 亚洲久久久国产精品| 窝窝影院91人妻| 老司机影院毛片| 欧美人与性动交α欧美软件| 午夜福利在线观看吧| 黄网站色视频无遮挡免费观看| 成人黄色视频免费在线看| 伊人久久大香线蕉亚洲五| 最新在线观看一区二区三区| 人妻 亚洲 视频| 日本猛色少妇xxxxx猛交久久| 狠狠精品人妻久久久久久综合| 国产成人影院久久av| 超碰成人久久| 少妇裸体淫交视频免费看高清 | 人人妻人人澡人人爽人人夜夜| 一区二区三区乱码不卡18| 最近最新中文字幕大全免费视频| 国产亚洲av高清不卡| 亚洲天堂av无毛| 亚洲全国av大片| 在线av久久热| 他把我摸到了高潮在线观看 | 亚洲 国产 在线| 18禁裸乳无遮挡动漫免费视频| 日韩 欧美 亚洲 中文字幕| 免费av中文字幕在线| 99re6热这里在线精品视频| 91精品国产国语对白视频| 最近中文字幕2019免费版| 热99国产精品久久久久久7| 欧美精品人与动牲交sv欧美| 老汉色∧v一级毛片| 极品人妻少妇av视频| 欧美av亚洲av综合av国产av| 两个人看的免费小视频| 欧美激情极品国产一区二区三区| 国产成+人综合+亚洲专区| 美女大奶头黄色视频| 动漫黄色视频在线观看| 性少妇av在线| av欧美777| 黑人欧美特级aaaaaa片| 久久女婷五月综合色啪小说| 下体分泌物呈黄色| 久久久国产成人免费| 新久久久久国产一级毛片| 午夜久久久在线观看| 在线精品无人区一区二区三| 女人爽到高潮嗷嗷叫在线视频| 啪啪无遮挡十八禁网站| 12—13女人毛片做爰片一| 久久久久久亚洲精品国产蜜桃av| 热99久久久久精品小说推荐| 少妇猛男粗大的猛烈进出视频| 日本a在线网址| 一二三四在线观看免费中文在| 免费观看a级毛片全部| 91成人精品电影| av一本久久久久| 大香蕉久久网| 色老头精品视频在线观看| 国产精品成人在线| videos熟女内射| 久久久精品94久久精品| 五月开心婷婷网| 国产老妇伦熟女老妇高清| 成人手机av| 精品人妻一区二区三区麻豆| 国产片内射在线| 亚洲中文av在线| 老司机影院成人| 成年人黄色毛片网站| 涩涩av久久男人的天堂| 国产片内射在线| 性高湖久久久久久久久免费观看| a级毛片黄视频| 亚洲天堂av无毛| 少妇精品久久久久久久| 亚洲情色 制服丝袜| 成在线人永久免费视频| 黑人巨大精品欧美一区二区mp4| 国产国语露脸激情在线看| 亚洲国产看品久久| 国产av国产精品国产| 51午夜福利影视在线观看| 国产精品国产三级国产专区5o| 国产国语露脸激情在线看| 精品视频人人做人人爽| 日日爽夜夜爽网站| 欧美在线一区亚洲| 午夜精品国产一区二区电影| 免费久久久久久久精品成人欧美视频| 国产欧美亚洲国产| 王馨瑶露胸无遮挡在线观看| 国产男女超爽视频在线观看| 国产日韩一区二区三区精品不卡| 99久久综合免费| 两个人免费观看高清视频| 国产人伦9x9x在线观看| 黄片小视频在线播放| 青青草视频在线视频观看| 女人久久www免费人成看片| 王馨瑶露胸无遮挡在线观看| 国产亚洲精品久久久久5区| 欧美性长视频在线观看| 啪啪无遮挡十八禁网站| 97精品久久久久久久久久精品| 亚洲五月婷婷丁香| 婷婷成人精品国产| 亚洲一卡2卡3卡4卡5卡精品中文| 午夜精品国产一区二区电影| 亚洲av男天堂| 天天影视国产精品| 高清在线国产一区| 深夜精品福利| 日日爽夜夜爽网站| 大香蕉久久网| 国产免费一区二区三区四区乱码| 国产av精品麻豆| 亚洲精品成人av观看孕妇| 日本撒尿小便嘘嘘汇集6| 久久久欧美国产精品| 黑人操中国人逼视频| 亚洲精华国产精华精| 91字幕亚洲| 欧美中文综合在线视频| 午夜福利免费观看在线| 精品少妇一区二区三区视频日本电影| 免费不卡黄色视频| 高清视频免费观看一区二区| 热99国产精品久久久久久7| 久久久国产成人免费| 精品视频人人做人人爽| 亚洲中文日韩欧美视频| 亚洲精品久久久久久婷婷小说| videos熟女内射| 淫妇啪啪啪对白视频 | 国产精品 欧美亚洲| 91成年电影在线观看| 少妇 在线观看| 桃花免费在线播放| www.999成人在线观看| 亚洲va日本ⅴa欧美va伊人久久 | 国产精品一二三区在线看| 一本色道久久久久久精品综合| 十八禁网站免费在线| 下体分泌物呈黄色| 极品人妻少妇av视频| 美女高潮喷水抽搐中文字幕| 免费在线观看视频国产中文字幕亚洲 | 18禁观看日本| av一本久久久久| 欧美日韩亚洲综合一区二区三区_| 999久久久精品免费观看国产| 桃花免费在线播放| 国产欧美日韩精品亚洲av| 一边摸一边抽搐一进一出视频| 三级毛片av免费| 亚洲国产精品999| 久久久久网色| 极品人妻少妇av视频| 99精品欧美一区二区三区四区| 五月天丁香电影| 精品国产一区二区久久| 1024香蕉在线观看| 免费日韩欧美在线观看| 热99re8久久精品国产| 亚洲第一青青草原| 久久久久网色| 精品久久久久久电影网| 美女高潮喷水抽搐中文字幕| 老熟妇乱子伦视频在线观看 | 久久久久久人人人人人| 男女国产视频网站| 性色av乱码一区二区三区2| 999久久久国产精品视频| av天堂久久9| 亚洲国产精品一区二区三区在线| 嫩草影视91久久| 久久久久久久久免费视频了| 精品福利永久在线观看| 涩涩av久久男人的天堂| av国产精品久久久久影院| 亚洲一区二区三区欧美精品| 丁香六月天网| 中文字幕av电影在线播放| 免费高清在线观看视频在线观看| 日韩视频在线欧美| 热re99久久国产66热| 制服诱惑二区| 日韩欧美一区二区三区在线观看 | 少妇的丰满在线观看| 性高湖久久久久久久久免费观看| 亚洲三区欧美一区| 国产男女内射视频| 精品国产乱码久久久久久男人| 亚洲精品国产av蜜桃|