• <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ù)作品欣賞
    劉邦不識字嗎?
    日本爱情动作片www.在线观看| 久久午夜综合久久蜜桃| 亚洲国产看品久久| 又大又黄又爽视频免费| 黄色怎么调成土黄色| 成人手机av| 欧美激情极品国产一区二区三区| 国产精品久久久久久精品古装| 97精品久久久久久久久久精品| 国产 精品1| 日韩中文字幕欧美一区二区 | 久久久欧美国产精品| 久久婷婷青草| 中文字幕制服av| 建设人人有责人人尽责人人享有的| 一本久久精品| 丰满少妇做爰视频| 黄色毛片三级朝国网站| 婷婷色综合www| 免费看av在线观看网站| 国产精品免费视频内射| 大片免费播放器 马上看| 亚洲欧洲精品一区二区精品久久久 | videosex国产| 欧美变态另类bdsm刘玥| 一个人免费看片子| 另类精品久久| 校园人妻丝袜中文字幕| 亚洲综合色惰| 狠狠婷婷综合久久久久久88av| 精品一品国产午夜福利视频| 亚洲精品第二区| 国产免费视频播放在线视频| 最近最新中文字幕免费大全7| 久久精品国产亚洲av高清一级| 久久狼人影院| 三级国产精品片| 搡老乐熟女国产| 男人操女人黄网站| 中文字幕人妻丝袜制服| 丝袜喷水一区| 又黄又粗又硬又大视频| 久久午夜综合久久蜜桃| 国产 一区精品| 久久精品人人爽人人爽视色| av在线app专区| 亚洲国产精品成人久久小说| 制服人妻中文乱码| 18禁裸乳无遮挡动漫免费视频| 久久99精品国语久久久| 成年美女黄网站色视频大全免费| 热99久久久久精品小说推荐| 久久久久久久国产电影| 一级毛片我不卡| 免费av中文字幕在线| 香蕉国产在线看| 午夜福利在线免费观看网站| 精品一区二区三卡| 天天躁夜夜躁狠狠躁躁| 亚洲欧美日韩另类电影网站| 免费在线观看视频国产中文字幕亚洲 | 高清黄色对白视频在线免费看| 久久精品人人爽人人爽视色| 人妻 亚洲 视频| 国产亚洲精品第一综合不卡| 久久精品国产a三级三级三级| 99热全是精品| 国产老妇伦熟女老妇高清| 久久午夜福利片| 免费观看性生交大片5| 在线观看国产h片| 777米奇影视久久| 国产精品亚洲av一区麻豆 | 夫妻性生交免费视频一级片| 免费大片黄手机在线观看| 亚洲国产看品久久| 欧美av亚洲av综合av国产av | 韩国av在线不卡| 婷婷色综合www| 性色avwww在线观看| 国产日韩一区二区三区精品不卡| 丰满乱子伦码专区| 中文精品一卡2卡3卡4更新| 久久精品亚洲av国产电影网| 国产女主播在线喷水免费视频网站| 男女午夜视频在线观看| 亚洲国产色片| 丰满饥渴人妻一区二区三| 少妇 在线观看| 狂野欧美激情性bbbbbb| 久久国产亚洲av麻豆专区| 国产人伦9x9x在线观看 | 999精品在线视频| 麻豆乱淫一区二区| 大片免费播放器 马上看| 亚洲欧美一区二区三区久久| 亚洲精品av麻豆狂野| 久久人人爽av亚洲精品天堂| 91在线精品国自产拍蜜月| av国产久精品久网站免费入址| www日本在线高清视频| 一级黄片播放器| 国产乱来视频区| 久久久久网色| 人人澡人人妻人| 高清av免费在线| 日韩三级伦理在线观看| 青春草亚洲视频在线观看| 一级毛片我不卡| 精品少妇内射三级| 一区二区av电影网| 亚洲国产欧美在线一区| 在线观看免费高清a一片| 99热全是精品| 少妇的逼水好多| 狠狠婷婷综合久久久久久88av| 午夜av观看不卡| 国产男人的电影天堂91| 高清不卡的av网站| 99热全是精品| 久久精品aⅴ一区二区三区四区 | 久久免费观看电影| 免费av中文字幕在线| 国产精品熟女久久久久浪| 国产亚洲精品第一综合不卡| 免费观看av网站的网址| 一本色道久久久久久精品综合| 亚洲色图综合在线观看| 香蕉精品网在线| 国产在线免费精品| 国产黄频视频在线观看| 国产精品一区二区在线不卡| 青春草国产在线视频| 久久婷婷青草| 国产精品一区二区在线观看99| 18在线观看网站| 久久久久久伊人网av| 久久韩国三级中文字幕| 亚洲精品第二区| 新久久久久国产一级毛片| www.av在线官网国产| 90打野战视频偷拍视频| 亚洲成av片中文字幕在线观看 | 人人澡人人妻人| 伦理电影免费视频| 久久精品亚洲av国产电影网| 欧美人与性动交α欧美精品济南到 | 母亲3免费完整高清在线观看 | 中文字幕人妻熟女乱码| 亚洲欧洲国产日韩| 久久久国产一区二区| 久久韩国三级中文字幕| www日本在线高清视频| 热re99久久精品国产66热6| 久久久久久人人人人人| 少妇的逼水好多| 伦理电影大哥的女人| 中文字幕另类日韩欧美亚洲嫩草| 制服丝袜香蕉在线| 麻豆av在线久日| 色吧在线观看| 中国国产av一级| 老司机亚洲免费影院| 亚洲欧美中文字幕日韩二区| 亚洲一区中文字幕在线| 久久国内精品自在自线图片| 婷婷色综合www| 中文乱码字字幕精品一区二区三区| 久久久久网色| 欧美日韩精品网址| 9热在线视频观看99| 色婷婷久久久亚洲欧美| 日本欧美国产在线视频| 亚洲av免费高清在线观看| 不卡av一区二区三区| 中文乱码字字幕精品一区二区三区| 欧美精品亚洲一区二区| 少妇熟女欧美另类| 日韩一区二区三区影片| 国产成人欧美| 国产一级毛片在线| 国产精品蜜桃在线观看| 国产成人欧美| 男人舔女人的私密视频| 国产日韩欧美亚洲二区| 波多野结衣av一区二区av| 久久狼人影院| 成人毛片a级毛片在线播放| av在线app专区| 美女国产高潮福利片在线看| 久久午夜综合久久蜜桃| 波多野结衣一区麻豆| 成人免费观看视频高清| 精品人妻一区二区三区麻豆| 天天躁日日躁夜夜躁夜夜| 18禁裸乳无遮挡动漫免费视频| 纵有疾风起免费观看全集完整版| 一区福利在线观看| 少妇人妻精品综合一区二区| 日本欧美国产在线视频| 啦啦啦在线免费观看视频4| 美女国产视频在线观看| 成人毛片60女人毛片免费| 亚洲成色77777| 高清在线视频一区二区三区| 亚洲三区欧美一区| 国产成人欧美| 欧美精品高潮呻吟av久久| 久久久欧美国产精品| 人妻一区二区av| 大陆偷拍与自拍| 亚洲精品久久久久久婷婷小说| 桃花免费在线播放| 我的亚洲天堂| 交换朋友夫妻互换小说| 人人妻人人澡人人看| 啦啦啦在线观看免费高清www| 国产亚洲午夜精品一区二区久久| 日韩中文字幕视频在线看片| av天堂久久9| 日本av手机在线免费观看| 校园人妻丝袜中文字幕| 欧美日韩国产mv在线观看视频| 亚洲av免费高清在线观看| 另类精品久久| 日本欧美视频一区| 人人妻人人爽人人添夜夜欢视频| 色网站视频免费| 丝袜人妻中文字幕| 亚洲欧美成人综合另类久久久| 女人精品久久久久毛片| 国产成人精品婷婷| 综合色丁香网| 国产精品熟女久久久久浪| 啦啦啦啦在线视频资源| 国产xxxxx性猛交| 日本黄色日本黄色录像| 免费女性裸体啪啪无遮挡网站| 一级片免费观看大全| 免费观看a级毛片全部| 国产成人精品婷婷| 日本av手机在线免费观看| 不卡av一区二区三区| 晚上一个人看的免费电影| 男女无遮挡免费网站观看| 免费在线观看完整版高清| 精品酒店卫生间| 少妇的逼水好多| 久久精品国产a三级三级三级| 国产精品久久久久久精品电影小说| 国产人伦9x9x在线观看 | 亚洲欧美成人综合另类久久久| 天天操日日干夜夜撸| 精品人妻熟女毛片av久久网站| 国产精品香港三级国产av潘金莲 | 国产女主播在线喷水免费视频网站| 日产精品乱码卡一卡2卡三| 可以免费在线观看a视频的电影网站 | 岛国毛片在线播放| 一级,二级,三级黄色视频| 久久这里只有精品19| 亚洲av男天堂| 久久久久人妻精品一区果冻| 9191精品国产免费久久| 欧美日韩av久久| 午夜影院在线不卡| 大陆偷拍与自拍| 国产男女内射视频| 日韩一区二区三区影片| 亚洲男人天堂网一区| 亚洲精品视频女| 欧美精品一区二区免费开放| 各种免费的搞黄视频| 色视频在线一区二区三区| 男人爽女人下面视频在线观看| 高清在线视频一区二区三区| 另类亚洲欧美激情| 精品国产一区二区三区久久久樱花| 又大又黄又爽视频免费| 国产亚洲最大av| 中文字幕精品免费在线观看视频| 欧美成人午夜免费资源| 男女无遮挡免费网站观看| 人人妻人人澡人人爽人人夜夜| 国产成人精品久久久久久| 一区在线观看完整版| 精品一区在线观看国产| av女优亚洲男人天堂| 欧美日韩亚洲高清精品| 久久精品国产亚洲av涩爱| 2022亚洲国产成人精品| 制服人妻中文乱码| 99热全是精品| 不卡视频在线观看欧美| 菩萨蛮人人尽说江南好唐韦庄| 中文字幕精品免费在线观看视频| 少妇被粗大的猛进出69影院| 午夜av观看不卡| 免费观看性生交大片5| 少妇人妻精品综合一区二区| 精品国产超薄肉色丝袜足j| 成年av动漫网址| 国产免费又黄又爽又色| videosex国产| 久久精品亚洲av国产电影网| 日韩av免费高清视频| 少妇猛男粗大的猛烈进出视频| 少妇 在线观看| 男人添女人高潮全过程视频| 国产成人精品在线电影| 成人国语在线视频| 王馨瑶露胸无遮挡在线观看| 最近2019中文字幕mv第一页| 亚洲欧美色中文字幕在线| videosex国产| 日本av免费视频播放| 97在线视频观看| videossex国产| 纯流量卡能插随身wifi吗| 亚洲欧美成人精品一区二区| 伊人亚洲综合成人网| 一本—道久久a久久精品蜜桃钙片| 免费观看无遮挡的男女| 日韩中文字幕欧美一区二区 | 午夜福利在线观看免费完整高清在| 国产极品天堂在线| 国产男女内射视频| 国产精品久久久久成人av| 久久久a久久爽久久v久久| 国产伦理片在线播放av一区| 精品人妻偷拍中文字幕| 亚洲av男天堂| 热re99久久国产66热| 建设人人有责人人尽责人人享有的| 日韩视频在线欧美| 99国产综合亚洲精品| 久久精品国产自在天天线| 男女国产视频网站| 最近手机中文字幕大全| 女人被躁到高潮嗷嗷叫费观| 秋霞伦理黄片| 久久精品国产亚洲av高清一级| 大片电影免费在线观看免费| 国产精品99久久99久久久不卡 | 亚洲精品在线美女| h视频一区二区三区| 大话2 男鬼变身卡| 99久久中文字幕三级久久日本| 亚洲精品av麻豆狂野| 三级国产精品片| 秋霞伦理黄片| 交换朋友夫妻互换小说| 18禁裸乳无遮挡动漫免费视频| 最近手机中文字幕大全| 在线观看人妻少妇| 午夜福利乱码中文字幕| 欧美激情极品国产一区二区三区| 免费大片黄手机在线观看| 美女脱内裤让男人舔精品视频| 狠狠婷婷综合久久久久久88av| 亚洲av.av天堂| 亚洲婷婷狠狠爱综合网| 亚洲av综合色区一区| 精品少妇黑人巨大在线播放| 如日韩欧美国产精品一区二区三区| 少妇猛男粗大的猛烈进出视频| 99久久人妻综合| 捣出白浆h1v1| 国产日韩欧美亚洲二区| 最近中文字幕2019免费版| 成人国产麻豆网| 一二三四中文在线观看免费高清| 妹子高潮喷水视频| 欧美激情 高清一区二区三区| 欧美日韩成人在线一区二区| 丰满乱子伦码专区| 九色亚洲精品在线播放| 一区福利在线观看| 欧美精品人与动牲交sv欧美| 久久午夜综合久久蜜桃| 中文字幕亚洲精品专区| 国产成人精品福利久久| 少妇被粗大猛烈的视频| 在线观看免费高清a一片| 欧美亚洲日本最大视频资源| 黄色一级大片看看| 国产成人免费观看mmmm| 下体分泌物呈黄色| 深夜精品福利| 男人操女人黄网站| 国精品久久久久久国模美| 老熟女久久久| 久久人妻熟女aⅴ| 少妇的丰满在线观看| 亚洲精品日本国产第一区| 久热这里只有精品99| 免费观看av网站的网址| 国产精品亚洲av一区麻豆 | 国产精品三级大全| 亚洲av综合色区一区| 亚洲成色77777| 日韩大片免费观看网站| 国产精品99久久99久久久不卡 | 天天影视国产精品| 永久免费av网站大全| 在线观看国产h片| 18禁动态无遮挡网站| 免费观看在线日韩| 中国三级夫妇交换| 亚洲精品第二区| 黄片小视频在线播放| 日韩精品有码人妻一区| 日日爽夜夜爽网站| 丝袜喷水一区| 久久久久久久久久久久大奶| 国产亚洲午夜精品一区二区久久| 亚洲精品日本国产第一区| 久久久精品免费免费高清| 亚洲激情五月婷婷啪啪| 久久综合国产亚洲精品| 国产精品蜜桃在线观看| 看十八女毛片水多多多| 精品午夜福利在线看| 97精品久久久久久久久久精品| 亚洲中文av在线| xxx大片免费视频| 国产精品蜜桃在线观看| 边亲边吃奶的免费视频| 日日爽夜夜爽网站| 精品国产一区二区三区四区第35| 久久99精品国语久久久| 欧美 亚洲 国产 日韩一| 国产不卡av网站在线观看| 最新中文字幕久久久久| 成年人午夜在线观看视频| 久久久久国产网址| 国产精品一区二区在线不卡| 亚洲伊人色综图| 国产黄色免费在线视频| 亚洲色图综合在线观看| 新久久久久国产一级毛片| 汤姆久久久久久久影院中文字幕| 色播在线永久视频| 国产av一区二区精品久久| 少妇人妻久久综合中文| 国产精品蜜桃在线观看| 一级,二级,三级黄色视频| 国产精品不卡视频一区二区| 久久国产精品男人的天堂亚洲| 亚洲精品自拍成人| 亚洲成人手机| 国产淫语在线视频| 成人手机av| 婷婷色综合www| 考比视频在线观看| 欧美日韩一级在线毛片| 老司机影院成人| 女性生殖器流出的白浆| 午夜久久久在线观看| 男人舔女人的私密视频| 边亲边吃奶的免费视频| 日韩成人av中文字幕在线观看| 新久久久久国产一级毛片| 免费高清在线观看视频在线观看| 日韩免费高清中文字幕av| 亚洲国产精品999| 少妇的丰满在线观看| 国产无遮挡羞羞视频在线观看| 国产免费福利视频在线观看| 精品少妇久久久久久888优播| 美女午夜性视频免费| 一区在线观看完整版| 欧美 亚洲 国产 日韩一| 午夜免费男女啪啪视频观看| 免费女性裸体啪啪无遮挡网站| 国产亚洲欧美精品永久| 成年av动漫网址| 一级片免费观看大全| 亚洲国产欧美网| 亚洲av在线观看美女高潮| 人人妻人人澡人人爽人人夜夜| 2018国产大陆天天弄谢| 一本久久精品| 日韩一区二区视频免费看| 国产国语露脸激情在线看| 老司机亚洲免费影院| 日韩,欧美,国产一区二区三区| 精品人妻一区二区三区麻豆| 永久网站在线| 卡戴珊不雅视频在线播放| 午夜福利视频在线观看免费| 香蕉国产在线看| 中文乱码字字幕精品一区二区三区| 午夜福利,免费看| 精品久久蜜臀av无| 婷婷成人精品国产| 国产福利在线免费观看视频| 久久精品久久精品一区二区三区| 97在线人人人人妻| 亚洲av国产av综合av卡| 亚洲国产欧美网| 久久毛片免费看一区二区三区| 日本vs欧美在线观看视频| 亚洲三区欧美一区| 日本-黄色视频高清免费观看| 黄片播放在线免费| 一级爰片在线观看| 日产精品乱码卡一卡2卡三| 性色avwww在线观看| 黄片小视频在线播放| 日韩大片免费观看网站| 国产成人欧美| 香蕉国产在线看| 如何舔出高潮| 久久99蜜桃精品久久| 久久久久国产精品人妻一区二区| 五月开心婷婷网| 黄色怎么调成土黄色| 另类亚洲欧美激情| 国产亚洲午夜精品一区二区久久| 国产一区有黄有色的免费视频| 一区在线观看完整版| 99热全是精品| 国产成人精品久久二区二区91 | 一级,二级,三级黄色视频| 久久午夜福利片| 亚洲第一青青草原| 中国国产av一级| 亚洲三级黄色毛片| 七月丁香在线播放| 亚洲国产欧美日韩在线播放| 亚洲精品aⅴ在线观看| 天天躁日日躁夜夜躁夜夜| 捣出白浆h1v1| 亚洲av电影在线观看一区二区三区| 天天影视国产精品| 精品久久久久久电影网| 超碰97精品在线观看| 在线观看免费日韩欧美大片| av女优亚洲男人天堂| 美女脱内裤让男人舔精品视频| 黄片播放在线免费| 久久狼人影院| 午夜激情久久久久久久| 一个人免费看片子| 日韩大片免费观看网站| 97在线人人人人妻| 最新中文字幕久久久久| 啦啦啦啦在线视频资源| 欧美日韩一级在线毛片| 天天躁夜夜躁狠狠躁躁| 日本色播在线视频| 亚洲av综合色区一区| 日韩一本色道免费dvd| 国语对白做爰xxxⅹ性视频网站| av有码第一页| 在线观看美女被高潮喷水网站| 在线精品无人区一区二区三| 啦啦啦视频在线资源免费观看| 欧美日本中文国产一区发布| 国产高清国产精品国产三级| 少妇被粗大的猛进出69影院| www.精华液| 一区二区日韩欧美中文字幕| 日韩制服骚丝袜av| 18禁动态无遮挡网站| 欧美亚洲 丝袜 人妻 在线| 一级a爱视频在线免费观看| 永久免费av网站大全| 免费日韩欧美在线观看| 熟女电影av网| 免费日韩欧美在线观看| a级毛片在线看网站| 人妻系列 视频| 这个男人来自地球电影免费观看 | 高清不卡的av网站| 老汉色av国产亚洲站长工具| 我的亚洲天堂| 伦精品一区二区三区| 久久人人爽av亚洲精品天堂| 亚洲人成网站在线观看播放| 自拍欧美九色日韩亚洲蝌蚪91| av在线老鸭窝| 一级毛片 在线播放| videos熟女内射| 日日撸夜夜添| 波多野结衣av一区二区av| 久久久久久伊人网av| 极品少妇高潮喷水抽搐| av免费在线看不卡| 欧美97在线视频| 国产高清国产精品国产三级| 午夜福利,免费看| 欧美日韩视频高清一区二区三区二| 亚洲久久久国产精品| 亚洲 欧美一区二区三区| 亚洲中文av在线| 黄色毛片三级朝国网站| 女性生殖器流出的白浆| 免费观看无遮挡的男女| 精品酒店卫生间| 精品久久久精品久久久| 满18在线观看网站| 一二三四在线观看免费中文在| 亚洲av国产av综合av卡| 精品国产超薄肉色丝袜足j| 欧美日韩精品成人综合77777| 婷婷成人精品国产| 菩萨蛮人人尽说江南好唐韦庄| 国产成人精品一,二区| 制服人妻中文乱码| 国产精品一区二区在线观看99| 热99国产精品久久久久久7| 1024视频免费在线观看| 久久久久久人人人人人|