• <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ù)作品欣賞
    劉邦不識字嗎?
    欧美3d第一页| 久久久久久久久久黄片| 黄色丝袜av网址大全| 久久久久九九精品影院| 国产午夜精品久久久久久一区二区三区 | 成人一区二区视频在线观看| 在线免费观看不下载黄p国产 | 久久欧美精品欧美久久欧美| 国产欧美日韩精品一区二区| 久久久国产精品麻豆| 身体一侧抽搐| 国内精品一区二区在线观看| 中文字幕熟女人妻在线| 黄色成人免费大全| 18美女黄网站色大片免费观看| 国产高潮美女av| 久久人妻av系列| 黄片大片在线免费观看| 91字幕亚洲| 每晚都被弄得嗷嗷叫到高潮| 少妇高潮的动态图| 三级男女做爰猛烈吃奶摸视频| 成人一区二区视频在线观看| 成年免费大片在线观看| 精华霜和精华液先用哪个| 国内少妇人妻偷人精品xxx网站| 校园春色视频在线观看| 嫩草影院入口| 国产精品一及| 欧美日本亚洲视频在线播放| 日本黄大片高清| 色综合婷婷激情| 一级毛片高清免费大全| 狂野欧美白嫩少妇大欣赏| 国产精品久久久久久人妻精品电影| 中文字幕人妻熟人妻熟丝袜美 | 午夜免费激情av| 久久久久久久久久黄片| 欧美精品啪啪一区二区三区| 成人av在线播放网站| 国产亚洲精品久久久com| 色老头精品视频在线观看| 夜夜爽天天搞| 男人和女人高潮做爰伦理| 免费人成视频x8x8入口观看| 丰满人妻一区二区三区视频av | 亚洲成人精品中文字幕电影| 亚洲精品日韩av片在线观看 | 久久这里只有精品中国| 婷婷精品国产亚洲av| 亚洲人成网站高清观看| 久久久久性生活片| 性欧美人与动物交配| 1024手机看黄色片| 久久精品国产自在天天线| 国产69精品久久久久777片| 搡老熟女国产l中国老女人| 午夜视频国产福利| netflix在线观看网站| 国产精品野战在线观看| 小蜜桃在线观看免费完整版高清| 国产亚洲欧美98| 很黄的视频免费| 国产淫片久久久久久久久 | 成人av在线播放网站| 久久久久久大精品| 男女做爰动态图高潮gif福利片| av天堂中文字幕网| 亚洲午夜理论影院| 18禁黄网站禁片午夜丰满| 一个人免费在线观看电影| 91在线观看av| 欧美丝袜亚洲另类 | 欧美精品啪啪一区二区三区| 久久久久久久亚洲中文字幕 | 亚洲精品日韩av片在线观看 | 欧美中文日本在线观看视频| 欧美黑人欧美精品刺激| 欧美+日韩+精品| a级毛片a级免费在线| 两性午夜刺激爽爽歪歪视频在线观看| 国产精华一区二区三区| 国产高清视频在线播放一区| 夜夜躁狠狠躁天天躁| 成年女人永久免费观看视频| a在线观看视频网站| 国产视频一区二区在线看| 欧美大码av| or卡值多少钱| 啦啦啦免费观看视频1| 国产精品野战在线观看| 国产精品久久电影中文字幕| 日韩高清综合在线| 国产精品免费一区二区三区在线| 精品一区二区三区视频在线观看免费| 国产乱人伦免费视频| 久久精品国产综合久久久| 又爽又黄无遮挡网站| 亚洲成人久久爱视频| 免费观看的影片在线观看| 琪琪午夜伦伦电影理论片6080| 久久精品人妻少妇| 免费av不卡在线播放| 一级毛片高清免费大全| 一区二区三区免费毛片| 亚洲av免费在线观看| 亚洲国产精品999在线| 欧美日韩黄片免| 熟女电影av网| 夜夜看夜夜爽夜夜摸| 日本成人三级电影网站| 亚洲专区国产一区二区| 黄色成人免费大全| 激情在线观看视频在线高清| 一级黄片播放器| 亚洲片人在线观看| 国产精品,欧美在线| 欧美日韩黄片免| 精品人妻偷拍中文字幕| 国产精品香港三级国产av潘金莲| 日本 av在线| 天堂动漫精品| 香蕉丝袜av| 91久久精品国产一区二区成人 | 久久6这里有精品| 在线观看午夜福利视频| 十八禁人妻一区二区| 国产麻豆成人av免费视频| 在线观看日韩欧美| av片东京热男人的天堂| 国产精品久久电影中文字幕| 国产探花在线观看一区二区| 99国产精品一区二区三区| 又紧又爽又黄一区二区| 在线观看av片永久免费下载| 丰满人妻一区二区三区视频av | 99热这里只有是精品50| 两个人的视频大全免费| 欧美一区二区亚洲| 国产中年淑女户外野战色| 欧美日韩国产亚洲二区| 99国产极品粉嫩在线观看| 男女床上黄色一级片免费看| 一级黄片播放器| 小蜜桃在线观看免费完整版高清| 一进一出抽搐gif免费好疼| 国产亚洲欧美在线一区二区| 欧美色欧美亚洲另类二区| 欧美另类亚洲清纯唯美| 一进一出抽搐gif免费好疼| 九色成人免费人妻av| 亚洲五月天丁香| xxx96com| 又粗又爽又猛毛片免费看| 亚洲av美国av| 欧美成人免费av一区二区三区| 久久精品亚洲精品国产色婷小说| 日本 av在线| 日韩亚洲欧美综合| 亚洲成av人片在线播放无| 男女午夜视频在线观看| 亚洲五月天丁香| 精品久久久久久久久久久久久| 亚洲专区国产一区二区| 国产黄片美女视频| 亚洲内射少妇av| 男女视频在线观看网站免费| 尤物成人国产欧美一区二区三区| 久久久久久人人人人人| 日韩欧美免费精品| www国产在线视频色| 国产老妇女一区| 高清日韩中文字幕在线| 精品欧美国产一区二区三| 亚洲午夜理论影院| 黄色女人牲交| 久久性视频一级片| 99国产精品一区二区三区| 日韩 欧美 亚洲 中文字幕| 在线观看午夜福利视频| 成人18禁在线播放| 欧美不卡视频在线免费观看| 亚洲精品日韩av片在线观看 | 小蜜桃在线观看免费完整版高清| 亚洲va日本ⅴa欧美va伊人久久| 美女免费视频网站| 亚洲欧美日韩高清专用| 又黄又粗又硬又大视频| 国产一区二区在线观看日韩 | 精品久久久久久久久久久久久| 在线观看免费视频日本深夜| 成人特级av手机在线观看| 中文字幕精品亚洲无线码一区| 可以在线观看毛片的网站| 亚洲成人中文字幕在线播放| 国产精品久久久久久人妻精品电影| 久久草成人影院| 欧美绝顶高潮抽搐喷水| 午夜影院日韩av| 亚洲欧美日韩卡通动漫| 亚洲在线观看片| 国内精品美女久久久久久| 国产欧美日韩精品一区二区| 欧美精品啪啪一区二区三区| 日韩欧美三级三区| 国产中年淑女户外野战色| 免费av不卡在线播放| 亚洲片人在线观看| 母亲3免费完整高清在线观看| 免费观看精品视频网站| 久久人人精品亚洲av| 亚洲狠狠婷婷综合久久图片| 国产v大片淫在线免费观看| 婷婷丁香在线五月| 在线a可以看的网站| 亚洲18禁久久av| 中文字幕精品亚洲无线码一区| 欧美乱色亚洲激情| 亚洲av日韩精品久久久久久密| aaaaa片日本免费| 中文资源天堂在线| 亚洲精品美女久久久久99蜜臀| 午夜激情福利司机影院| 国产精品98久久久久久宅男小说| 免费看美女性在线毛片视频| 黑人欧美特级aaaaaa片| 国产精品香港三级国产av潘金莲| 18禁黄网站禁片午夜丰满| h日本视频在线播放| 香蕉av资源在线| 97碰自拍视频| 偷拍熟女少妇极品色| 国产精品久久久久久精品电影| 亚洲五月天丁香| 久久国产精品人妻蜜桃| 黄色女人牲交| 黄片大片在线免费观看| 国产欧美日韩精品一区二区| 国产97色在线日韩免费| 久久精品91蜜桃| 男女之事视频高清在线观看| 免费看日本二区| 日韩精品中文字幕看吧| 国产av不卡久久| 欧美色欧美亚洲另类二区| 成人特级黄色片久久久久久久| 亚洲熟妇熟女久久| 99久国产av精品| 国产乱人视频| 亚洲欧美日韩卡通动漫| av专区在线播放| www.色视频.com| 在线观看日韩欧美| 午夜福利在线观看免费完整高清在 | 亚洲人成网站高清观看| 欧美又色又爽又黄视频| 18美女黄网站色大片免费观看| 久久久精品大字幕| 少妇丰满av| 毛片女人毛片| 久久久久免费精品人妻一区二区| 国产亚洲欧美在线一区二区| 精品99又大又爽又粗少妇毛片 | 少妇的丰满在线观看| 三级国产精品欧美在线观看| 嫩草影院精品99| 亚洲不卡免费看| 亚洲一区二区三区色噜噜| 99国产综合亚洲精品| 最新在线观看一区二区三区| 国产精品98久久久久久宅男小说| 哪里可以看免费的av片| 精品日产1卡2卡| 亚洲天堂国产精品一区在线| 午夜免费成人在线视频| 性色avwww在线观看| 一级毛片高清免费大全| 国产高清激情床上av| 久久久色成人| 欧美精品啪啪一区二区三区| 亚洲国产欧美人成| 91麻豆av在线| 美女大奶头视频| 18禁裸乳无遮挡免费网站照片| a级一级毛片免费在线观看| 亚洲国产精品999在线| 午夜精品在线福利| 天天一区二区日本电影三级| 少妇的逼好多水| 可以在线观看的亚洲视频| 在线天堂最新版资源| av在线天堂中文字幕| 国产真实乱freesex| 天天一区二区日本电影三级| 久久香蕉国产精品| 亚洲国产精品久久男人天堂| 亚洲精品亚洲一区二区| 男女视频在线观看网站免费| 午夜免费男女啪啪视频观看 | 亚洲美女黄片视频| 99在线视频只有这里精品首页| 欧美精品啪啪一区二区三区| 九色成人免费人妻av| 欧美日韩亚洲国产一区二区在线观看| 老汉色∧v一级毛片| 一a级毛片在线观看| 成年人黄色毛片网站| 香蕉久久夜色| 久久中文看片网| 美女高潮的动态| 中文亚洲av片在线观看爽| 亚洲av电影在线进入| 精品久久久久久久久久久久久| 97碰自拍视频| 18美女黄网站色大片免费观看| 午夜日韩欧美国产| 99热只有精品国产| 亚洲久久久久久中文字幕| 精品人妻偷拍中文字幕| 色精品久久人妻99蜜桃| 午夜免费男女啪啪视频观看 | 亚洲av一区综合| 亚洲第一电影网av| 午夜福利在线在线| 黄片大片在线免费观看| 伊人久久精品亚洲午夜| 亚洲av美国av| 亚洲av电影不卡..在线观看| 每晚都被弄得嗷嗷叫到高潮| 90打野战视频偷拍视频| 国产成人av激情在线播放| 91av网一区二区| xxxwww97欧美| 亚洲成人精品中文字幕电影| 免费看日本二区| 国产午夜精品久久久久久一区二区三区 | 国产综合懂色| 岛国视频午夜一区免费看| 一本精品99久久精品77| 亚洲av美国av| 中文字幕久久专区| 国产三级黄色录像| 在线国产一区二区在线| 叶爱在线成人免费视频播放| 欧美日韩国产亚洲二区| 精品免费久久久久久久清纯| 色老头精品视频在线观看| 51午夜福利影视在线观看| 黑人欧美特级aaaaaa片| 国产精品av视频在线免费观看| 99热6这里只有精品| 亚洲成人免费电影在线观看| 欧美在线一区亚洲| 在线播放无遮挡| 宅男免费午夜| 99久久精品国产亚洲精品| 全区人妻精品视频| 日本 欧美在线| 啦啦啦韩国在线观看视频| avwww免费| 老司机午夜十八禁免费视频| 国产精品1区2区在线观看.| 午夜福利免费观看在线| 精品久久久久久久久久免费视频| 精华霜和精华液先用哪个| 午夜激情福利司机影院| av天堂中文字幕网| 亚洲av第一区精品v没综合| 国内精品久久久久精免费| 亚洲片人在线观看| 欧美3d第一页| 日韩免费av在线播放| 欧美xxxx黑人xx丫x性爽| 一本综合久久免费| 亚洲av成人av| 观看免费一级毛片| 男女做爰动态图高潮gif福利片| 国产日本99.免费观看| 国产99白浆流出| 最近最新免费中文字幕在线| 一个人看视频在线观看www免费 | 深夜精品福利| 九色国产91popny在线| 亚洲,欧美精品.| 欧美日韩亚洲国产一区二区在线观看| 国产免费一级a男人的天堂| 国产精品 国内视频| 又黄又粗又硬又大视频| 听说在线观看完整版免费高清| 精品一区二区三区av网在线观看| 日韩欧美一区二区三区在线观看| 日本撒尿小便嘘嘘汇集6| 亚洲精品一卡2卡三卡4卡5卡| 亚洲国产精品成人综合色| 国产午夜精品论理片| 有码 亚洲区| 国产 一区 欧美 日韩| 亚洲 欧美 日韩 在线 免费| 亚洲五月天丁香| 免费看光身美女| 哪里可以看免费的av片| 亚洲国产欧美网| 欧美性感艳星| 手机成人av网站| 久久亚洲真实| 少妇高潮的动态图| 丝袜美腿在线中文| 欧美绝顶高潮抽搐喷水| 色综合站精品国产| 99精品在免费线老司机午夜| 亚洲激情在线av| 18禁裸乳无遮挡免费网站照片| 欧美日韩瑟瑟在线播放| 免费看美女性在线毛片视频| 天堂av国产一区二区熟女人妻| 亚洲av美国av| 中文字幕人妻熟人妻熟丝袜美 | 精品熟女少妇八av免费久了| 亚洲中文字幕日韩| 久久香蕉国产精品| 特大巨黑吊av在线直播| 日韩成人在线观看一区二区三区| 中文字幕精品亚洲无线码一区| 三级国产精品欧美在线观看| tocl精华| 亚洲中文字幕一区二区三区有码在线看| 久久精品综合一区二区三区| 久久久久久久亚洲中文字幕 | 美女大奶头视频| 国产精品美女特级片免费视频播放器| 亚洲一区二区三区色噜噜| 日韩欧美精品免费久久 | 亚洲第一欧美日韩一区二区三区| 国产真实伦视频高清在线观看 | 亚洲国产精品合色在线| 天美传媒精品一区二区| 亚洲中文日韩欧美视频| 欧美日韩中文字幕国产精品一区二区三区| 亚洲av成人精品一区久久| 精品国产超薄肉色丝袜足j| 久久久精品欧美日韩精品| 国产精品久久电影中文字幕| 欧美日韩综合久久久久久 | 亚洲国产精品999在线| 欧美色欧美亚洲另类二区| 99久久久亚洲精品蜜臀av| 俺也久久电影网| av福利片在线观看| 18美女黄网站色大片免费观看| 欧美国产日韩亚洲一区| 中文亚洲av片在线观看爽| 丝袜美腿在线中文| 少妇高潮的动态图| 亚洲av一区综合| 日本一本二区三区精品| 好男人电影高清在线观看| 成年人黄色毛片网站| 97碰自拍视频| 美女cb高潮喷水在线观看| 国产高潮美女av| 桃色一区二区三区在线观看| 国产主播在线观看一区二区| 一级毛片女人18水好多| 免费观看人在逋| 黄色成人免费大全| 在线观看免费午夜福利视频| 91av网一区二区| 天堂影院成人在线观看| 成人特级av手机在线观看| 日本熟妇午夜| 禁无遮挡网站| 欧美又色又爽又黄视频| 日韩欧美在线二视频| 免费高清视频大片| 国产免费男女视频| 黄片大片在线免费观看| 午夜免费成人在线视频| 日本三级黄在线观看| www.999成人在线观看| 韩国av一区二区三区四区| av视频在线观看入口| 女同久久另类99精品国产91| 国产美女午夜福利| 少妇丰满av| 国产成人a区在线观看| 久久久精品欧美日韩精品| 白带黄色成豆腐渣| 成人av在线播放网站| 亚洲国产精品999在线| 中文在线观看免费www的网站| 九色成人免费人妻av| 老司机深夜福利视频在线观看| 欧美成人a在线观看| 国产高潮美女av| 久久草成人影院| 美女黄网站色视频| 在线观看舔阴道视频| 最近最新免费中文字幕在线| 精品人妻偷拍中文字幕| 精品久久久久久成人av| 亚洲熟妇熟女久久| 757午夜福利合集在线观看| 国产亚洲精品久久久com| 精品人妻1区二区| 亚洲不卡免费看| 国产一区二区在线av高清观看| 老司机深夜福利视频在线观看| 俺也久久电影网| 久久久久亚洲av毛片大全| 精品99又大又爽又粗少妇毛片 | 搡女人真爽免费视频火全软件 | 亚洲av免费在线观看| 亚洲国产日韩欧美精品在线观看 | 国产乱人视频| 法律面前人人平等表现在哪些方面| АⅤ资源中文在线天堂| www日本在线高清视频| 久久久色成人| 五月伊人婷婷丁香| 国产探花在线观看一区二区| 欧美成人性av电影在线观看| 老汉色∧v一级毛片| 亚洲精品456在线播放app | 国产免费av片在线观看野外av| 国产欧美日韩精品一区二区| 日韩av在线大香蕉| 精品电影一区二区在线| 亚洲不卡免费看| 丰满乱子伦码专区| 欧美最新免费一区二区三区 | 成人亚洲精品av一区二区| 老汉色av国产亚洲站长工具| 岛国在线观看网站| 国产亚洲av嫩草精品影院| 色在线成人网| 亚洲在线观看片| 窝窝影院91人妻| 日日干狠狠操夜夜爽| 免费在线观看影片大全网站| 久久草成人影院| www.999成人在线观看| 在线十欧美十亚洲十日本专区| 日韩亚洲欧美综合| 欧美bdsm另类| 午夜视频国产福利| 成人无遮挡网站| 99国产极品粉嫩在线观看| 一本精品99久久精品77| 在线观看美女被高潮喷水网站 | 岛国在线观看网站| 久99久视频精品免费| 精品国产三级普通话版| 18+在线观看网站| a在线观看视频网站| 村上凉子中文字幕在线| 最近最新中文字幕大全电影3| 香蕉丝袜av| 色哟哟哟哟哟哟| 亚洲av电影不卡..在线观看| 久久久久久久午夜电影| 免费无遮挡裸体视频| 可以在线观看的亚洲视频| 九九久久精品国产亚洲av麻豆| 免费av观看视频| 一本久久中文字幕| 日韩欧美精品v在线| 久久久精品大字幕| 日韩av在线大香蕉| 日韩欧美国产一区二区入口| 制服人妻中文乱码| 免费看美女性在线毛片视频| 麻豆久久精品国产亚洲av| 精品久久久久久成人av| a级毛片a级免费在线| 中文亚洲av片在线观看爽| 少妇的丰满在线观看| 亚洲内射少妇av| 波野结衣二区三区在线 | 国内少妇人妻偷人精品xxx网站| 久久伊人香网站| 俄罗斯特黄特色一大片| 国产极品精品免费视频能看的| 女人被狂操c到高潮| 一二三四社区在线视频社区8| 老熟妇乱子伦视频在线观看| 真人一进一出gif抽搐免费| 午夜免费成人在线视频| 精品99又大又爽又粗少妇毛片 | 在线免费观看的www视频| 国产精品乱码一区二三区的特点| 欧美日韩一级在线毛片| 亚洲激情在线av| 99国产精品一区二区三区| 国产美女午夜福利| 色吧在线观看| 日韩精品中文字幕看吧| 老司机午夜十八禁免费视频| 国产高清视频在线观看网站| 欧美激情久久久久久爽电影| 中文字幕精品亚洲无线码一区| 欧美黄色片欧美黄色片| 搡老熟女国产l中国老女人| 好男人在线观看高清免费视频| 亚洲美女视频黄频| 亚洲午夜理论影院| 1024手机看黄色片| 国产精品久久久人人做人人爽| 特级一级黄色大片| 久久午夜亚洲精品久久| 婷婷六月久久综合丁香| 熟妇人妻久久中文字幕3abv| 久久久国产成人免费| 香蕉av资源在线| 又黄又爽又免费观看的视频| 99久久精品国产亚洲精品|