Twist-assisted optoelectronic phase control in two-dimensional (2D) Janus heterostructures

Atomically thin two-dimensional (2D) Janus materials and their Van der Waals heterostructures (vdWHs) have emerged as a new class of intriguing semiconductor materials due to their versatile application in electronic and optoelectronic devices. Herein, We have invstigated most probable arrangements of different inhomogeneous heterostructures employing one layer of transition metal dichalcogenide, TMD (MoS2, WS2, MoSe2, and WSe2) piled on the top of Janus TMD (MoSeTe or WSeTe) and investigated their structural, electronic as well as optical properties through first-principles based calculations. After that, we applied twist engineering between the monolayers from 0\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$^{\circ } \rightarrow$$\end{document}∘→ 60\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$^{\circ }$$\end{document}∘ twist angle, which delivers lattice reconstruction and improves the performance of the vdWHs due to interlayer coupling. The result reveals that all the proposed vdWHs are dynamically and thermodynamically stable. Some vdWHs such as MoS2/MoSeTe, WS2/WSeTe, MoS2/WSeTe, MoSe2/MoSeTe, and WS2/MoSeTe exhibit direct bandgap with type-II band alignment at some specific twist angle, which shows potential for future photovoltaic devices. Moreover, the electronic property and carrier mobility can be effectively tuned in the vdWHs compared to the respective monolayers. Furthermore, the visible optical absorption of all the Janus vdWHs at \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\theta$$\end{document}θ = 0\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$^{\circ }$$\end{document}∘ can be significantly enhanced due to the weak inter-layer coupling and redistribution of the charges. Therefore, the interlayer twisting not only provides an opportunity to observe new exciting properties but also gives a novel route to modulate the electronic and optoelectronic properties of the heterostructure for practical applications.

Tremendous research attention is currently being focused on two-dimensional (2D) nanomaterials and their heterostructures due to their rich physical properties and diverse technological applications [1][2][3] .The transition metal dichalcogenides (TMDs) such as MoS 2 , WS 2 , MoSe 2 , and WSe 2 are the most intriguing members of the 2D family due to their excellent electronic, thermal, and optical properties 4 .In addition to having graphene-like structures, these materials have direct band gaps in the visible-near IR range and provide large carrier mobility 5 , flexibility 6 , transparency, and fast lithium diffusion 7 .These properties make them suitable for digital circuits 8 , logic transistors 9 , light-emitting diodes (LEDs) 10 , photovoltaic devices 11 , chemical sensors 12 , and Li-ion battery anodes 13 etc.These atomically thin TMDs are large-gap higher-order topological crystalline insulators, which are protected by C 3 rotational symmetry 14 .One of the critical aspect of TMDs is that they can be deposited onto flexible substrates and survive the stress and strain compliance of flexible supports 15 .These fascinating properties of TMDs originate from the quantum confinement and surface effects that arise during the transition of an indirect bandgap to a direct bandgap semiconductor when bulk materials are scaled down to their monolayers.The TMD monolayers are dynamically and thermodynamically stable compared to other 2D materials.Therefore, they were widely synthesized by various experimental methods, including chemical vapor deposition (CVD) 16 , physical vapor deposition (PVD) 17 , mechanical exploitation (ME) 18 , and hydrothermal synthesis 19 etc.In addition to this, 2D Janus materials (such as 2D Janus TMDs) are also dragging tremendous attention due to their unique structural, physical, and chemical properties 20 , which differ from conventional 2D materials.Recent studies also demonstrate the successful synthesis of 2D Janus MoSSe and WSSe monolayers by modified CVD technique 21 .These materials offer a novel property due to their broken mirror symmetry with distinct atomic types on their upper and lower side.The general formula of artificial 2D Janus TMDs is MXY, where M is a transition metal atom, which is sandwiched between two different chalcogen atoms X, and Y.As a result of the different atomic radii and electronegativities of the X and Y elements, the charge distributions between M-X and M-Y layers are

Methodology
First-principles based density functional theory (DFT) calculations were performed using the Quantum ATK 44 , using the projector augmented wave (PAW) method.The plane-wave cut-off energy is set to 790 eV and the Monkhorst-Pack 45 mesh of k-points with 9 × 9 × 1 and 11 × 11 × 1 points are used to sample the Brillouin zone for monolayers and heterostructures, respectively.A 15 Å vacuum layer along the z-direction is employed to avoid interactions between the neighboring layers.The lattice constant and atomic positions are relaxed until the force on each atom is less than 0.01 eV/Å.The convergence criterion of the total energy is set at 10 −5 eV/atom.The phonon band structure was calculated with the finite displacement method by using the PHONOPY code 46 , where a 3 × 3 × 1 supercell was used.The Kubo-Greenwood formula 47 was used to calculate the susceptibility tensor for estimating the optical spectrum as given in Eq. (1).
Here, p i nm and p j mn belong to ith and jth component of the dipole matrix element between state n and m respectively, V is the volume of the material, δ is broadening and f nk indicates Fermi function evaluated at the band energy E n (k).

Results and discussion
In this present study, we have considered two Janus (MoSeTe and WSeTe) monolayers and four TMD monolayers ( MoS 2 , MoSe 2 , WS 2 and WSe 2 ).We have systematically investigated the stability, electrical, and optical proper- ties of Janus MoSeTe and WSeTe monolayers, as well as vertically stacked heterostructures, made up of Janus monolayers and TMDs.The weak van der Waals forces between the monolayers do not cause substantial atomicscale alterations and usually maintain the original features of the monolayer in the heterostructure.Hence, the combination of two different monolayers into a heterostructure can offer some unique and fascinating properties as compared to its monolayer. (1) Structural and electronic properties of the monolayers.The top and side views of the optimized Janus MoSeTe and WSeTe monolayers are shown in Fig. 1a,b,e,f.The structure of Janus monolayers (XSeTe, X = Mo, and W) is constructed in three planes, where the X atom is sandwiched between the layer composed of Se and Te layers.The bond length between the Mo-Se, Mo-Te, W-Se, and W-Te are different due to the difference in atomic number between the atoms.The TMDs monolayers ( XY 2 , X = Mo, W, and Y = S, Se) also consist of three layers, where the X layer is stacked between two similar Y layers (as shown in Fig. S1 in the supplementary information (SI)).The computed lattice parameters for all the TMDs and Janus monolayers are presented in  S3a in the SI, the projected density of states (PDOS) of MoSeTe shows that d xy (Mo) is the main contribu- tor to the valence band (VB), with small contributions from the p x (Te) states, whereas the conduction band (CB) originates from the p z (Te) and d yz (Mo) states.In the case of WSeTe monolayer, the valence band mainly originates from the d xy (W) and p y (Te) states, but the conduction band is mostly contributed by the d x 2 −y 2 (Te) with a negligible contribution of p y (Te).
Heterostructures and their properties.We have proposed nine different types of vertically stacked heterostructure (TMD/Janus: 8 and Janus/Janus: (1) since the TMDs and Janus monolayers are well-studied material for electronic and optoelectronic applications.In addition, the structural similarity between the TMDs and Janus monolayers allows the formation of better quality interface.All the nine vdWHs are: MoS 2 /MoSeTe, MoS 2 / WSeTe, WS 2 /MoSeTe, WS 2 /WSeTe, MoSe 2 /MoSeTe, MoSe 2 /WSeTe, WSe 2 /MoSeTe, WSe 2 /WSeTe, and MoSeTe/ WSeTe.In addition to this, the relative twist angle between the monolayers also gives rise to lattice reconstruction and enhances the performance of the vdWHs due to the complex interlayer coupling.We considered all the (2) possible twist angles between the monolayers in the range of 0-60 • .The twist angles were chosen based on the number of atoms present in the supercell of the heterostructure and the minimum lattice mismatch between the monolayers.All nine heterostructures fall into two categories, allowing four rotation angles for XSe 2 /XSeTe (X = Mo, W) and MoSeTe/WSeTe and three rotation angles for XS 2 /XSeTe (X = Mo, W) heterostructure, by satisfying the above criteria.Figure S4a in the SI shows the top and side views of the Janus MoSeTe/WSeTe heterostructure at θ = 0 • , whereas Fig. S5a-d in the SI represents additional schematic representations of the heterostructure for different rotation angles.We have also calculated the interlayer distance for each rotation angle by considering the minimum energy criteria.An example of estimation of interlayer distance for Janus WSe 2 /WSeTe vdWH is shown in Fig. S6 in the SI.The estimated number of atoms present in the supercell, lattice mismatch, and interlayer distance for Janus XS 2 /XSeTe, XSe 2 /XSeTe, and MoSeTe/WSeTe vdWHs are presented in the SI of Tables S2-S4 respectively.
Stability.The dynamical stability of 2D vdWHs plays a very vital role in its practical application in nanoelectronic and optoelectronic devices.We have calculated phonon bandstructure along the high symmetric points ( Ŵ-Y-X-Ŵ ) for all the Janus vdWHs at θ = 0 • .The phonon band structures of Janus XS 2 /XSeTe and XSe 2 /XSeTe vdWHs are presented in Fig. 2a-h and Janus MoSeTe/WSeTe vdWH is shown in Fig. S4b in the SI.The absence of imaginary vibration frequencies confirms the dynamical stability of the heterostructure.Interestingly, the soft acoustic mode of Janus MoSeTe and WSeTe monolayer has a frequency of 3 THz and 3.3 THz at the Y-point, respectively.This frequency value becomes softer when the Janus monolayers from a heterostructure with TMDs monolayer.
To address the thermodynamic stabilities of the Janus heterostructure, the interface cohesive energies were calculated using Eq. ( 3).
Here, E CE and E X/Y represent the cohesive energy and total energy of the vdWHs, respectively.The E X and E Y show the energy of the constituent monolayers.The variation of cohesive energy for all the Janus vdWHs at different rotation angles is presented in Fig. 3.The low value of cohesive energies indicates that the Janus and TMD monolayers can form a stable interface.The Janus vdWHs are quite favorable because all TMDs and Janus monolayers have comparable structural geometry.The minimum cohesive energy obtained for Janus MoSeTe/ WSeTe vdWHs is − 137.56 meV/Å −2 at θ = 21.79 • .It is worth mentioning that the cohesive energy is minimum, at which the lattice mismatch is also minimum, allowing each monolayer to make a better and strongly coupled interface.Overall, the cohesive energy varies between 11 and 138 meV/Å −2 for all the Janus vdWHs.This value is also well comparable with other similar type of 2D heterostructures such as XSSe/Mg(OH) 2 (X = Mo, W) 51 , Graphene/MoS 2 52 , and WTe 2 /Graphene 53 .To check the mechanical stability, we computed the elastic constants of Janus vdWHs using the Born-Huang criteria 54  (3) It is observed that Young's modulus is higher along the x-direction than in the y and z directions.There was a maximum value of 137.71 N/m for MoS 2 /WSeTe, which is lower than MoSi 2 As 2 monolayer 55 , but comparable to other 2D materials, such as Ti 2 C (137 N/m) 56 , B 2 S (137 N/m to 143 N/m) 57 , TiC 2 (131.17N/m) 58 and higher than the Phosphorene (92.4 N/m) 59 , silicene (62 N/m) 57 .The ZX and ZY components of Poission's ratio has a higher value along the in-plane directions, where the highest value has been estimated to be 0.2949, which is approximately the same as that of monolayer MoSi 2 As 4 55 .
Electronic properties.The energy band structure of the Janus vdWHs is critically important for determining the performance of the material for different technological applications.The band structures of Janus MoSeTe/ WSeTe vdWH at θ = 0 • and 21.79 • are calculated and presented in Fig. 4a,c respectively.The Janus MoSeTe/ WSeTe vdWH has the same indirect bandgap semiconducting property as monolayers but with a reduced value.The calculated band gaps for this vdWHs are 0.67, 0.49, 0.78, 0.71, and 0.46 eV at θ = 0 • , 16.10 • respectively.To characterize the behavior of electrons at the band edges, the projected density of states (PDOSs) has been calculated and plotted in Fig. 4b,d.The PDOS analysis reveals that at θ = 0 • , the VBM is dominated by the d x 2 −y 2 orbital of W, and the CBM is contributed by the d z 2 orbital Mo atom.In contrast, the reverse nature of PDOS was observed for the rest of the rotation angles, which means the VBM and CBM reside on d-orbitals of Mo and W, respectively.The origin of the CBM and VBM is due to the presence of two different Janus monolayers, implying the MoSeTe/WSeTe vdWHs have type-II band alignment with the staggered gaps, which is also supported by the projected band structure of the WSeTe/MoSeTe heterostructure ( θ = 0 • , as illustrated in Fig. S11a in the SI), in which orange band lines represent the contribution of the WSeTe layer and blue lines indicate the MoSeTe contribution.We can easily see that WSeTe creates valence band maxima, while MoSeTe creates conduction band minima.As shown in Fig. S11b,c in SI, the partial charge density above and below the Fermi level also suggested type-II band alignment.A dense cloud of electrons surround the W atom below the Fermi level, while a cloud of electrons surround the Mo atom above the Fermi level.Such materials are beneficial for photocatalytic water splitting and enable high-efficiency solar energy conversion.The variation of band gap for all the Janus vdWHs at different twisting angles are shown in Fig. 5.The bandgap values of Janus vdWHs lie between 0.2 and 1.5 eV.The energy bandstructure along the high symmetric paths for the two different classes of Janus vdWHs ( MoS 2 /MoSeTe and WSe 2 /WSeTe) are depicted in Figs.S7 and S8 in the SI, respectively.Mostly, the CBM lies in between X and Y-points of the BZ for all twisting angles, but the VBM resides on X-point for θ = 0 • and at Ŵ-point for the rest of the twisting angles.The maximum bandgap obtained is 1.46 eV for Janus WSe 2 /WSeTe vdWH at θ = 21.79 • .Similarly, we observed a minimum band gap of 0.24 eV for Janus MoS 2 /MoSeTe vdWH at θ = 0 • .We also observed that the Janus XS 2 /XSeTe vdWHS have a lower bandgap as compared to Janus XSe 2 /XSeTe vdWHs.This may be due to the increased bond length and atomic number   www.nature.com/scientificreports/ between X-S to X-Se.As shown in Fig. S9 in the SI, the band diagram of MoSeTe/WSeTe heterostructure at θ = 0 • includes vdW and dipole corrections.The band structure shows semiconducting behavior with an indirect bandgap of 0.67 eV, which is approx.equal to the PBE calculation.We have also compared the bandgap of Janus vdW heterostructures at θ = 0 • between the PBE calculation and vdW + dipole corrections, as shown in Fig. S10 in the SI.We have not found any significant difference in the electronic structure for all 9 Janus vdW heterostructures between PBE calculation and PBE + vdW + dipole correction.The PDOS analysis confirms that all the Janus vdWHs at different rotation angle have type-II band alignment.Surprisingly, MoS 2 /MoSeTe ( θ = 19.11 ) have direct bandgap with type-II band alignment, which is highly desirable to enhance the efficiency of photovoltaic solar cells.The bandstructure difference among the vdWHs mainly arises from the change in interlayer coupling induced by the twisting angle.The Janus vdWHs show an electronic phase cross-over from direct to indirect bandgap semiconductor since the hybridization of atomic orbitals changes significantly at the Fermi level due to the interlayer rotation.The hybridization extent determines the value of the band gap.The lower band gap ofthe Janus vdWHs can improve the optical absorption in visible light and near-ultraviolet regions.
The main reason behind the change in bandgap for Janus vdWHs at different rotation angles can be described by calculating interlayer charge transfer.
To understand the charge transfer in the Janus vdWHs, we plot the charge density difference, as shown in Fig. 6a.The planar-averaged charge density difference is calculated using the Eq. ( 4).
Here, ρ H (x, y, z)dxdy , ρ MoSeTe (x, y, z)dxdy , and ρ WSeTe (x, y, z)dxdy are the charge density at the (x,y,z) point for the Janus MoSeTe/WSeTe vdWHs, MoSeTe, and WSeTe monolayers respectively.Based on the Bader charge analysis, the calculated value of the charge transfer for Janus MoSeTe/WSeTe heterostructure at θ = 0 • through the interface region is estimated at �ρ = 0.812 × 10 −3 e (as shown in Fig. 6a).The positive values denote charge accumulation, and the negative values denote charge depletion.The results show that a charge rearrangement occurs, and electrons accumulate at the interface.We also calculated the Barder charge for all the heterostructures at different twist angles and found that the value is positive for all the twist angles.To get a clear picture of charge accumulation at the interface, we have plotted isosurface plot, as shown in Fig. 6b.The cyan and pink color represent the charge depletion and accumulation, respectively.As can be seen, the charge accumulation mostly occurs near the interface region of MoSeTe and WSeTe monolayer.However, the charge redistribution strength gets decreased as we further move away from the interface region.This is due to the weak van der Waals interaction between the Janus monolayers.The accumulated charge forms an interlayer electric field, which ultimately affects the charge transfer between the monolayers.The most significant difference between the interlayer electronic interaction in Janus vdWHs is that there is discernible interlayer charge transfer due to the presence of different X and chalcogen atoms in monolayers.The interlayer charge transfer will induce an electrostatic interlayer interaction, which will give rise to an energy level shift in opposite directions for the two constituent monolayers.Furthermore, the twisting angle, interlayer distance, lattice mismatch, and geometry of the Janus vdWHs determine the amount of charge transfer between the monolayers.These Janus vdWHs have direct and indirect bandgap with type-II band alignment to facilitate the separation of electrons and holes, which is useful for optoelectronic and photovoltaic devices 60 .
The higher carrier mobility is another important indicator for determining the performance of the heterostructure in photovoltaic and photocatalytic applications.We have quantitatively estimated the carrier mobility for the Janus vdWHs in the x-direction by using the acoustic phonon limited method by using Eq. ( 5).
Here, T is the temperature, m * is the effective mass in the transport direction, and is the elastic modulus.E is the total energy after applying strain ( δ = l/l 0 ) and S 0 is the area at equilibrium.In addition E l ( E/δ ) is the diffusion constant, where E is the energy shift of the band edge of CBM or VBM with respect to the vacuum level.The estimated electron mobility of the Janus MoSeTe and WSeTe monolayers are 1.07 × 10 4 and 2.33 × 10 4 cm 2 /Vs and for holes are 3.11 × 10 4 and 7.79 × 10 4 cm 2 /Vs respectively.The hole mobilities of Janus monolayers are higher than their electron mobilities due to the relatively smaller deformation potentials.Earlier studies revealed that the conventional TMD monolayers have low carrier mobility 61,62 .The electron mobilities of the Janus vdWHs can also be enhanced when the Janus monolayer is stacked into vertical heterostructures with TMDs monolayer.High electron mobility ( ∼ 10 4 cm 2 V −1 s −1 ) have been observed for MoS 2 /WSeTe, MoSe 2 /MoSeTe, and WSe 2 /WSeTe vdWHs, as well as high hole mobility ( ∼ 10 4 cm 2 V −1 s −1 ) for MoSeTe/WSeTe, MoS 2 /MoSeTe, MoSe 2 /WSeTe, and WS 2 /MoSeTe vdWHs.The variation of electron and hole mobilities for Janus vdWHs at θ = 0 • is depicted in Fig. 6c.Moreover, high and low electron mobilities reach 4.61 × 10 4 cm 2 /Vs for Janus WSe 2 /MoSeTe and WS 2 /WSeTe vdWH at θ = 0 • respectively.These calculated carrier mobilities are quite high as compared to conventional 2D vdWHs.The high carrier mobility in vdWHs is attributed to the low effective mass (Fig. S12 in SI) and more dispersed band at the Fermi level, which would induce fast migration of the carriers and suppress the recombination of photo generated electron and hole pairs.Optical properties.The optical properties of the vertically stacked Janus vdWHs at different twisting angles are studied by computing the complex dielectric function as given below, (4) �ρ(z) = ρ H (x, y, z)dxdy − ρ MoSeTe (x, y, z)dxdy − ρ WSeTe (x, y, z)dxdy. (5) Here, ε 1 (ω ) and ε 2 (ω ) represents the real and imaginary part of dielectric constant respectively.The ε 2 (ω ) can be calculated by summing all the empty states in the Brillouin zone, while the ε 1 (ω ) of the dielectric func- tion can be obtained by using the Kramers-Kronig relation 63 .We have calculated the optical properties in the energy window of 0-5 eV for all the Janus vdWHs.The dielectric function is calculated along xx/yy (E ⊥ z) and zz (E z) directions.All the Janus vdWHs exhibit strong anisotropic nature due to their crystal symmetry, but we have considered only the maximum contribution, which is along yy-directions for the majority of the case.
The real static part of the dielectric function is the direct measure of refractive index (n).The higher the value of n, the better will be the probability of absorbing light.The refractive index is 1.27, 1.51, 2.84, 2.68, 2.63, and 2.38 for MoSeTe, WSeTe, MoS 2 , MoSe 2 , WS 2 , and WSe 2 monolayers respectively.Table 1 lists the estimated refractive index for nine Janus vdWHs at various rotation angles.The value of n behaves differently for different heterostructures as well as rotation angles.We have also observed that the maximum enhancement in n is at θ = 0 • for all the Janus vdWHs.This enrichment is due to the perfect interlayer coupling between the TMDs and Janus monolayers at that rotation angle.In case of Janus MoSeTe/WSeTe vdWHs, the n increases up to a value of 3.87 at θ = 0 • .The rotation angles give an extra freedom to use the same material for different applications with a low and high value of the refractive index.The low value of refractive index materials has a great potential in many photonic device applications due to their high reflective behavior, whereas the high value of refractive index materials are useful in storing energies such as solar cell.The imaginary part of the dielectric function reflects the transitions from occupied to unoccupied bands.The optical absorption coefficients of the heterostructures are obtained using the Eq. ( 7): Here, ω and α are utilized to express angular frequency and absorption coefficient.
(a) (b) (c) The absorption coefficient for Janus MoSeTe/WSeTe heterostructure is plotted in Fig. 7.The starting point of the absorption edge represents the optical bandgap which means the minimum amount of energy required for the transition of electron from VBM to CBM.For Janus XS 2 /XSeTe vdWHs, the order of absorption coefficient with respect to twist angle is θ = 0 • > 30  for Janus XSe 2 /XSeTe (MoSeTe/WSeTe) vdWHs.Thus, we can conclude that the maximum contribu- tion of absorption coefficient is at θ = 0 • for all the Janus vdWHs.Figure S16 in the SI represents the variation of absorption coefficients at θ = 0 • for a photon energy of 3 eV.The maximum absorption coefficient obtained is 4.34 × 10 5 cm −1 for Janus WSe 2 /MoSeTe VdWH at θ = 0 • , enabling the Janus vdWH to be a promising solar light absorber material in the visible region.This value of absorption coefficient is also quite high as compared to other 2D heterostructures.This blue shifting nature of the abortion spectrum at θ = 0 • is due to the weak interlayer coupling and the formation of a suitable interface between the monolayers.Other twisting angles such as 16.10  , and 43.90 • are active in between 1 and 3 eV and thereafter become completely transparent.The absorption peaks arise due to the electronic transition from the occupied states to the empty states.For example, The transition from d x 2 −y 2 (W) to d z 2 (Mo) in MoSeTe/WSeTe ( θ = 0 • ) vdWH resulted in a prominent peak at a photon energy of 2.6 eV.The value of the absorption coefficient also well matches the refractive index.As we can see from Table 1, with the twist between the monolayers, the refractive index decreases, which results into the low value of the absorption coefficient.The optical properties of the Janus vdWHs completely depend on the hybridization of orbitals at the Fermi level due to the interlayer coupling.Furthermore, the Janus MoS 2 /MoSeTe ( θ = 19.11 ) in the infrared region with direct bandgap, which is an essential requirement for infrared sensors.Therefore, the relative twist between the monolayers provides a new route to manipulate the properties of a heterostructure.

Conclusion
In summery, we have systematically studied the stability, electronic, and optical properties of nine different Janus vdWHs ( MoS 2 /MoSeTe, MoS 2 /WSeTe, WS 2 /MoSeTe, WS 2 /WSeTe, MoSe 2 /MoSeTe, MoSe 2 /WSeTe, WSe 2 /MoSeTe, WSe 2 /WSeTe, and MoSeTe/WSeTe) through first-principles calculations.Our results show that when TMD monolayers are stacked vertically with Janus monolayers, an intrinsic electric field appears because of the lack of mirror symmetry and charge accumulation, which leads to the origin of versatile properties in Janus vdWHs.In addition, the interlayer twist also induces an electronic phase crossover from indirect to direct bandgap at some specific rotation angles.It is interesting to note that the Janus vdWHs have direct bandgap with type-II band alignment, which is an essential parameter for photocatalytic water splitting applications.The interfacial interaction also induces high carrier mobility, efficient charge separation and high value of refractive index.Furthermore, the Janus vdWHs exhibit dramatically reduced direct band gaps and more efficient optical absorption from the UV to visible light range.The enhancement rate of absorption coefficient of the Janus vdWHs at θ = 0 • ) is quite high ( ∼ 10 5 cm −1 ) in the visible region due to the perfect inter- layer coupling between the monolayers.Our findings not only provide a compelling platform for exploring the behavior of the Janus vdWHs but also provides theoretical guidance for the designing of novel 2D nanodevices.

Figure 1 .
Figure 1.(a,b,e,f) Represent the top, and side view, (c,g) indicate the phonon dispersion curve, and (d,h) represents the bandstructure and density of states (DOS) of Janus MoSeTe and WSeTe monolayer respectively.The green and purple color represent the spin-↑ and spin-↓ states, respectively.The black dotted line stands for the Fermi level.

Figure 3 .
Figure 3.The variation of cohesive energy of Janus vdWHs as a function of rotation angles.

Figure 6 .
Figure 6.(a) Plane-integrated charge density difference and (b) three-dimensional (3D) charge density difference with the isosurface value of 4.9 × 10 −5 eÅ −3 of Janus MoSeTe/WSeTe vdWH at θ = 0 • .The dash- dot lines represent the interface.The cyan and pink color represent the charge depletion and accumulation respectively.(c) The variation of mobility of the Janus vdWHs at θ = 0 • .
Table S1 of the SI.The phonon dispersion of Janus MoSeTe and WSeTe monolayer is plotted in Fig.

Table 1 .
The variation of refractive index (n) at different rotation angle ( θ ) of Janus vdWHs.