Tailoring the anisotropic effect of Janus In2 XY (X/Y = S, Se, Te) monolayers toward realizing multifunctional optoelectronic device applications

Anisotropic effect of two-dimensional materials remains one of the most attractive properties, introducing an additional degree of freedom for tuning physical performances. In this work, we investigate the anisotropies of electronic, transport, piezoelectric, and optoelectronic properties of Janus In2 XY (X/Y = S, Se, Te) monolayers (J-In2 XY MLs) by performing first-principles calculations. We find that such J-In2 XY MLs possess moderate bandgap (2.07–2.29 eV), high carrier mobility (∼103 cm−2V−1s−1), visible light absorption (∼105 cm−1), large out-of-plane piezoelectric response (d31=0.35 pmV−1) and ultra-soft mechanical nature (C11=32.84 Nm−1) . We construct a kind of J-In2 XY-based phototransistor to investigate the optoelectronic properties under linearly polarized light. We find that the low recombination probability of photogenerated carriers ensured by anisotropic effect enhances the photocatalytic potential of J-In2 XY MLs. And the pivotal role induced by anisotropy in photocurrent can cause a prominent on/off ratio (∼100), considerable responsivity (0.038 AW−1, 0.036 AW−1) and external quantum efficiency (10.6%, 11.5%). Our study provides an avenue for the design of future anisotropic J-In2 XY-based multifunctional optoelectronic device.

Lately, an emerging class of van der Waals layered 2D semiconductors that belong to the III-VI chalcogenides combined by the formula MX (M = Ga, In; X = S, Se, Te), has drawn growing attention due to its extraordinary carrier mobility, excellent mechanical flexibility, superior light-matter interaction, and abundant polytypes [15][16][17][18]. Typically, the first member of monochalcogenides that is 2D InSe possesses tunable direct bandgap and can be synthesized by pulsed laser deposition [19], mechanical exfoliation [20], and chemical vapor deposition [21]. 2D InSe FET exhibits high intrinsic electron mobility of 10 3 cm 2 (V s) −1 and superior I on /I off ratio of 10 8 at room temperature [22]. Feng et al pronounced 2D InSe photodetectors showed a broadband photoresponse from ultraviolet (UV)-visible to near-infrared with the photoresponsivity of 10 4 AW −1 and photodetectivity of 10 13 Jones [23]. In addition, 2D GaTe [24][25][26][27][28] and InTe [29,30] with indirect bandgap (1.34 and 1.29 eV) have been proposed as one of the promising candidates for the next-generation electronic and optoelectronic devices. Moreover, for the group III-VI family, strikingly different from the aforementioned honeycomb-type hexagonal lattice structures, more complicated and lower symmetry monoclinic crystal with C2 m −1 space group can be obtained by choosing different growth patterns [25,26], which exhibit strong anisotropic nonlinear optical properties and can be applied to phototransistors [24]. In particular, the anisotropic effects can add one more level of complexity and realize extraordinary nanodevice applications that are impossible for isotropic structures.
Breaking structural symmetry is an effective strategy to tune the electronic structure of 2D materials. For example, a J-MoSSe ML can be fabricated by fully replacing the top-layer S of ML MoS 2 with Se [31]. This kind of Janus structures bring an additional degree of freedom and exhibit more fascinating characteristics than traditional transition metal dichalcogenides (TMDs), such as large vertical piezoelectric responses [32][33][34][35][36]. Interestingly, for the J-MXY (M = Ga, In; X/Y = S, Se, Te) MLs, some progresses show that their bandgaps are in the range of 0.89-2.03 eV, and possess in-plane piezoelectric coefficients up to 8.47 pm V −1 , four times more than their binary counterpart MLs, and so on [37,38]. However, few reports have uncovered that a single system based on J-MXY (M = Ga, In; X/Y = S, Se, Te) MLs can realize excellent multifunctional optoelectronic devices in terms of anisotropic effect.
In this work, we investigate J-In 2 XY (X/Y = S, Se, Te) MLs to conduct their potential applications in terms of first-principles calculations. Firstly, we examine the dynamical, thermal, mechanical and energy stabilities of three configurations. Subsequently, we explore the electronic, mechanical, and optical anisotropic effects of J-In 2 XY MLs. We find that J-In 2 XY MLs possess excellent out-of-plane piezoelectric responses and flexibility, while the high and anisotropic carrier mobility and light absorption can promote the separation and transfer of photogenerated electrons and holes, ensuring high efficiency of photocatalytic reactions. Lastly, we explore direction-dependent optoelectronic properties of J-In 2 XY p-n junctions, which show a large photocurrent on/off ratio, significant photo responsivities and external quantum efficiencies along the y direction of p-n junctions. Our study reveals the promising applications of J-In 2 XY MLs as a multifunctional building block in optoelectronic devices.

Methodology
All first-principles calculations in this presentation are carried out with the Vienna ab initio simulation package (VASP) [39,40] based on density functional theory (DFT) with the plane-wave basis projector augmented wave method [41]. The generalized gradient approximation (GGA) with Perdew-Burke-Ernzerhof (PBE) formalism is used for describing the exchange-correlation potential [42]. The plane-wave cutoff energy ofh 2 |K + G| 2 /2m = 500 eV is applied to take into account the ion-electron interactions. The Γ-centered 24 × 4 × 1 k-point grids within the Monkhorst-Pack scheme are executed to sample the Brillouin zone in the self-consistent field calculations [43]. The geometric optimization can be terminated until the Hellmann-Feynman forces on each atom become smaller than 0.01 eV Å −1 and the energy difference between two sequential steps is less than 10 −5 eV. The vacuum space of ∼20 Å is introduced along the c-axis to avoid artificial interactions between adjacent layers. Since PBE/GGA formalism systematically underestimates the electronic bandgap, the Heyd-Scuseria-Ernzerhorf (HSE06) hybrid functional is selected to provide a more accurate value for bandgap and band alignments. The dipole correction [44] along the c-axis is employed to handle inherent dipole moments as well as work functions. A 1 × 7 × 1 supercell is selected to acquire the phonon dispersion relations by the small displacement method implemented in the PHONOPY package [45]. The thermal stability of J-In 2 XY MLs is examined by performing ab initio molecular dynamics (AIMD) simulations during 144 atoms at 300 K with a total simulation time of 6 ps and time steps of 3 fs. The charge transfers between individual atoms are performed by using the Bader charge analysis technique [46]. The VASP within the density functional perturbation theory (DFPT) and small displacement methodology are carried out to obtain piezoelectric stress coefficients and linear-elastic properties, respectively. The optical response is especially important relative to photocatalysis and photoelectric applications. The absorption coefficient α(ω) is calculated as the following formula [47]: where c is the velocity of light, ω denotes the angular frequency of light, ε 1 (ω) and ε 2 (ω) refer to the real and imaginary parts of the dielectric function. The ε 1 (ω) is calculated by using the Kramers-Kronig relation [47], and the ε 2 (ω) of the frequency-dependent is determined based on the relation: where the subscripts c and v represent the conduction and valence band states, respectively, ω k stand for the weight of the k-point; and u ck denotes the cell periodic part of the orbitals at the k-point. Here, we mainly evaluate the maximum absorption coefficient and the overall optical properties. In order to save calculation cost, the optical properties under excitonic effect are not considered. The quantum transport calculations of the optoelectronic FET are carried out using the QUANTUM Atomistix ToolKit [48] based on DFT and non-equilibrium Green's function (NEGF) [49]. The electron-photon interaction is considered in the Hamiltonian asĤ =Ĥ 0 + (e/m 0 )A ·P with H 0 being the Hamiltonian without the electron-light interaction, and the e represents the electron charge, m 0 denotes the free electron mass,P and A are the momentum operator and the electromagnetic vector potential, respectively. The density mesh cut-off of 80 Hartree and 1 ×1 ×112 k-points mesh is adopted to the self-consistent calculation. The photocurrent in electrode (α) with spin (σ) is calculated using first-order perturbation theory within the 1st Born approximation, as follows, where e andh refer to the electron charge and reduced Planck's constant, and the effective transmission coefficients T α (E) are obtained as [50]: where G >/< refer to the greater/lesser Green's function including electron-photon interactions, the Γ α and f α stand for the linewidth and Fermi distribution function of the electrode (left or right), respectively.

Results and discussion
The structure of J-In 2 XY MLs is designed by replacing the surface layer of InTe MLs with S and Se atoms. In contrast to its hexagonal lattice phase, the layered monoclinic phase with weak van der Waals interactions is more stable and possesses C2 m −1 space group with 2 m −1 point group symmetry. As illustrated in figure 1(a), the J-In 2 SeTe ML can be viewed as In atomic layer sandwiched between Te and Se atoms, and the rectangular unit cell consists of 12 indium, six selenium, and six tellurium atoms. The optimized lattice parameters of J-In 2 SeTe (a = 25.11 Å; b = 4.23 Å), J-In 2 STe (a = 24.59 Å; b = 4.15 Å) and J-In 2 SSe (a = 24.22 Å; b = 3.99 Å) MLs reveal in-plane anisotropy derived from the low symmetry introduced by the pseudo-one-dimensional structure. As listed in table 1, the bond lengths (d In-In ) of J-In 2 SeTe, J-In 2 STe, and J-In 2 SSe MLs are 2.83, 2.82 and 2.82 Å, the distance between upper and lower surface layers (h) are 6.02, 5.95 and 5.79 Å, respectively, which are very close to the previous results of the isomorphic system [18]. It should be noted that the difference of electronegativity between X and Y atoms leads to interfacial charge transfer and induces a large internal electric field (E int ). The Bader charge analysis confirms that Se and Te atoms obtain electrons and charge transfer is more favorable than Se atoms, as evident in figure 1(a).
To verify the structural stability of J-In 2 XY MLs, the cohesive energy (E coh ) per atom is calculated to quantitatively compare the energy stability as follows: where E In2XY is the total energy of J-In 2 XY ML unit cell, and E In , E X and E Y stand for the single-atom energies of In, X and Y under vacuum, a, b and c represent the atom number of In, X and Y in the unit cell, respectively. From table 1, it is clearly seen that E coh of J-In 2 STe ML (−3.41 eV atom −1 ) is between those of J-In 2 SeTe (−3.27 eV atom −1 ) and J-In 2 SSe (−3.60 eV atom −1 ). To evaluate the dynamical stability, the phonon spectra of J-In 2 SeTe and J-In 2 SY (Y = Se, Te) MLs are calculated as depicted in figures 1(b) and S1 (supplemental materials). One can see that the phonon dispersions of J-In 2 XY MLs are free from any imaginary frequencies through the whole Brillouin zone. Note that all phonon branches are both nondegenerate owing to the in-plane anisotropy that each J-In 2 XY MLs possesses 69 optical phonon branches and three acoustic phonon branches. The thermal stability of both J-In 2 XY MLs is predicted by   carrying out AIMD simulations at 300 K for 6 ps, with a time step of 3 fs. Neither bond breaking nor structural reconstruction is observed in the resultant snapshot, and the total energy fluctuated around a fixed value without rising, as given in figures 1(c) and S1, implying that J-In 2 XY MLs are thermal stability. Albeit the phonon spectra and AIMD simulations provide essential conditions, the mechanical stability of J-In 2 XY MLs is also systematically examined. Such as the calculated elastic stiffness tensors of J-In 2 SeTe ML is C 11 = 32.84, C 12 = 8.74, C 22 = 37.75 and C 66 = 9.15 N m −1 , which satisfies the Born criteria [51] (C 11 C 22 > C 2 12 and C 66 > 0), as listed the table 1. Figure 1(d) shows the atom-dependent projected band structure and projected density of states (PDOSs) of J-In 2 SeTe ML. Clearly, J-In 2 SeTe ML has a direct bandgap of 1.38 eV at the Γ point of Brillouin zone. The main contribution to the conduction-band minimum (CBM) comes from In atoms while that of In, Se and Te atoms to the valence-band maximum (VBM) is comparable. Besides, HSE06 hybrid functional is considered to provide a more accurate bandgap. The band structures of J-In 2 XY MLs at the level of GGA/PBE and HSE06 are illustrated in figures 2(a)-(c). For the HSE06 algorithm, apart from J-In 2 SSe ML, which is predicted to be an indirect bandgap of 2.29 eV, J-In 2 SeTe and J-In 2 STe ML are obviously characterized by a direct bandgap semiconductor with a bandgap of 2.07 and 2.15 eV, respectively. Interestingly, the difference between the VBM of direct and indirect characteristics of J-In 2 XY MLs is only ∼100 meV, hence the bandgap type (G T ) can be easily tuned by strain engineering [52]. Considering that the characteristics of band structures of J-In 2 XY MLs are almost consistent based on the GGA/PBE and HSE06 algorithms, the effect of uniaxial strain along the x-and y-directions from −4% to 4% in the scheme of GGA/PBE on the band structures are investigated, and the indirect-direct band transition is given in figures 2(d) and (e). It is found that J-In 2 SSe ML is always an indirect bandgap semiconductor under x or y uniaxial strain, regardless of tensile or compressive strain. However, for J-In 2 SeTe ML, the bandgap transition from direct-to-indirect occurs at a tensile strain along the x direction, while that is just the opposite for J-In 2 STe ML. Unlike the uniaxial strain along the x direction, J-In 2 STe ML has an indirect bandgap under tensile strain along the y direction, and the bandgap transition from direct-to-indirect occurs from −1% to −2% for J-In 2 SeTe ML. In this case, strain engineering can be used to tune their electronic properties.
Piezoelectricity in layered material with interlayer van der Waals interactions has received considerable interest as the energy conversion bridges between mechanical and electrical energies. However, the piezoelectric response of most 2D layered materials is confined to exist only at the basal plane, which limits the operation of piezoelectric devices [32]. Here J-In 2 XY MLs are purposed to exhibit an additional piezoelectric degree of freedom arising from the out-of-plane symmetry breaking in which an electric dipole moment is generated. Theoretically, the piezoelectric response can be characterized by the relaxed-ion piezoelectric stress tensor e ijk and piezoelectric strain tensor d ijk , which is combined the contributions from ionics and electronics. Based on the Voigt notation, the e ijk and d ijk can be connected by elastic stiffness tensor C ijk , i.e. e ik = d ij C jk .
The piezoelectric anisotropy of in-plane and out-of-plane is also taken into account, thus there are four non-zero elements of e ij for the rectangle supercell J-In 2 XY MLs, i.e. e 11 , e 12 , e 31 and e 32 , where i represents the direction being polarized, j denotes the direction in which the strain is applied, and numbers 1, 2 and 3 stand for x, y, and out-of-plane directions, respectively. Therefore, we have the relations:  Moreover, we consider the net dipole moments of J-In 2 XY MLs along the out-of-plane direction to compare the out-of-plane polarization responses. The net dipole moments along the out-of-plane direction, namely z-direction, are calculated to be 0.18, 0.29 and 0.11 eÅ for J-In 2 SeTe, J-In 2 STe and J-In 2 SSe MLs, respectively. Interestingly, the vacuum energy level of the surface layer can be regulated by the net dipole moment of the Janus structures, which leads to the difference of electrostatic potential energy. In figure 4(a), the discrepancies between two surface layers are 0.34, 0.54 and 0.23 eV, respectively, which allow the oxidation or reduction potential of water splitting to shift downward and better serve photocatalysis, a new reaction mechanism of efficient photocatalysis was first proposed by Li et al [33]. To make it clear, the traditional and the new mechanism illustrations for photocatalytic water splitting are depicted in the left and right sections of figure 4(b). For the former, the standard reduction potential for hydrogen evolution reaction (HER) is −4.44 eV, while that of the oxidation potential for oxygen evolution reaction (OER) is −5.67 eV. Therefore, the bandgap of photocatalytic materials is not less than 1.23 eV, which is necessary to ensure HER and OER. However, for the latter, the internal electric field induced by the out-of-plane polarization bends the vacuum level of layer Y downward by ∆Φ , while the oxidation/reduction potential is still −5.67/−4.44 relative to the vacuum level of layer X/Y, thus the limitation of the bandgap for photocatalytic materials is redefined as E G = 1.23 − ∆Φ eV.
Although numerous semiconductors satisfy the bandgap criterion, their band alignments are not viable for the two half-reactions standard potential of water splitting (2 H + (aq) + 2e − → H 2 (g), H 2 O → 2 H + + 1/2O 2 ). To test the mentioned requirements for oxidation/reduction potential and further clarify the mechanism of photocatalytic water splitting, the band-edge positions of J-In 2 SeTe ML calculated by HSE06 are aligned concerning the vacuum level, as illustrated in figure 4(c). Clearly, the work functions (Φ = E vac − E F ) of Se and Te layers induced by polarization are 5.36 and 5.02 eV, respectively, resulting in the corresponding HER potential (−4.44 eV) relative to the CBM (−3.96 eV) moved down by 0.34 eV. Thus, the J-In 2 SeTe ML is predicted to exhibit the excellent photocatalytic capacity for hydrogen production by water splitting owing to the CBM has greatly exceeded the standard reduction potential (E H + /H2 ). More interesting, the E int directed by Te to Se can accelerate efficient spatial separation and transfer of photogenerated carriers and suppress their recombination, which facilitates for HER and OER. Also, we evaluate the band-edge position of relative vacuum levels of J-In 2 STe and J-In 2 SSe MLs as well as the corresponding binary parent structure InTe ML, as displayed in figure 4(d). Evidently, the band-edge positions of J-In 2 XY MLs straddle the standard redox potentials of water in strongly acidic (pH = 0) conditions. Note that qualified photocatalytic material should still be suitable over a wide pH range, and the standard potential of both E H + /H2 and E H2O/O2 sensitively varies with PH according to the Nernst equations [54]: E red H + /H2 = −4.44 eV + PH × 0.059 eV and E ox H2O/O2 = −5.67 eV + PH × 0.059 eV. As evident in figure 4(d), the CBM and VBM of J-In 2 XY MLs completely straddle both the potential of HER and OER in strongly acidic (pH = 0) and neutral (pH = 7) environments. However, the corresponding binary InTe ML is incapable of covering the potential of OER under pH = 0, indicating that J-In 2 XY MLs exhibit more superior photocatalytic ability for overall water splitting than the corresponding binary parent materials under extreme conditions. A promising photocatalytic material should possess high carrier mobility in addition to the two criteria investigated above, which facilitate the rapid transfer of electrons and holes into HER and OER before recombination. Therefore, we address the carrier mobility of J-In 2 XY MLs along x-and y-directions using the deformation potential theory [55], µ (x/y) = , where m d = m * ix m * iy refers to the average effective mass (i = e for electrons, i = h for holes, and x and y are the transport directions), The effective Table 2. Effective mass mx * (my * ), deformation potential E1x (E1y), 2D elastic modulus C2Dx (C2Dy), and carrier mobility µx (µy) of J-In2XY MLs along Γ-X (Γ-Y) directions for electron (e) and hole (h).
/S 0 stands for the corresponding elastic constant of the 2D system, (E − E 0 ) and S 0 are the total energy shift and the area of initial unit cell, k B andh denote the Boltzmann constant and Planck constant, room temperature T = 300 K is used, E 1(x/y) = ∆E/ ε is the deformation potential, and ∆E is the band edge drift of CBM for electrons and VBM for holes under uniaxial x or y strain. The evaluated values are summarized in table 2 and the corresponding fitting details are given in the supplemental materials ( figure S4). Clearly, all the MLs exhibit outstanding carrier mobility, where the electron mobility of J-In 2 SSe ML along the x-direction amounts to ∼ 4.53 × 10 3 cm −2 V −1 s −1 and the hole in the y-direction is ∼ 769.99 cm −2 V −1 s −1 , is one order of magnitude higher than the theoretically evaluated for MoS 2 ML (∼ 0.2 × 10 3 cm −2 V −1 s −1 and ∼ 72.16 cm −2 V −1 s −1 for electron and hole mobilities, respectively) [56]. Moreover, it is found that the electron and hole mobilities of J-In 2 XY MLs are significantly discrepancy, which may enhance electron-hole separation efficiency and photocatalytic activity [57].
What's more, the carrier mobilities of J-In 2 XY MLs show strong in-plane directional anisotropy. For instance, the electron mobility of J-In 2 SeTe ML along the x-direction is ∼ 1.46 × 10 3 cm −2 V −1 s −1 , much higher than that (∼ 0.43 × 10 3 cm −2 V −1 s −1 ) in the y-direction, and the anisotropic ratio [max(µ x , µ y ) min(µ x , µ y )] of electron mobility is about 3.36, which is comparable to the cases of SiAs ML (4.03) [58] and 2D ReS 2 sheet (3.1) [59]. To more intuitively reflect the anisotropy of electronic structure and carrier mobility, the three-dimensional valleys for VBM and CBM of J-In 2 SeTe ML are depicted in figure 5(a), and the corresponding 2D contour view is given in figures 5(b) and (c). Clearly, the valleys of VBM and CBM both show different dispersion profiles in the K x (Γ − X) and K y (Γ − Y), respectively, which can induce a high degree of anisotropy of electronic properties. Furthermore, we fit the quadratic function curve of the energy with respect to the K x and K y in the near Γ point (see figures 5(d) and (e)). The trend of the curve is consistent with the effective mass approximation, and by unit conversion, the carrier effective mass (m * i(x/y) ) and quadratic coefficient (S VBM/CBM ) can be defined as m * i(x/y) ≈ 3.8m 0 /S VMB/CBM , the details of calculation are shown in the supplemental materials. For the case of J-In 2 SeTe ML, the presence of electron effective mass in the x-direction is 0.30m 0 , and 0.25m 0 for the y-direction, while the hole effective mass along the x-direction (0.26m 0 ) is much smaller than that along the y-direction (5.84m 0 ). In that sense, the high carrier mobility with strong in-plane anisotropy can effectively inhibit the recombination of free electrons and holes [60,61], which can improve the overall efficiency of photocatalytic water splitting.
Additionally, the absorption capacity of solar light is an important indicator to confirm the efficiency of photocatalysis, and also plays a crucial role in the field of photovoltaic devices. In our case, the optical responses of J-In 2 XY MLs are examined by computing their frequency-dependent dielectric functions and absorption coefficients [α(ω)]. The imaginary part of the dielectric function (ε 2 ) and [α(ω)] versus the energy of the incident light polarization along x(E ∥x ) and y(E ∥y ) directions are depicted in figures S5 and 6, respectively. Notably, the significant optical anisotropy can be directly observed in the visible to UV region, and the zoomed-in image around 0.5-1 eV shows that there is still anisotropy in the near-infrared region because of their low lattice symmetry, in which the absorption intensity along the x(E ∥x ) direction is higher than that in the y(E ∥y ) direction in the near-infrared to visible light region for all J-In 2 XY MLs. Moreover, the absorption spectrum of J-In 2 XY MLs overlaps massively under the incident AM1.5 solar flux in the visible region, and the first conspicuous absorption peak also occurs in this area, it is evident that J-In 2 XY MLs are sensitive to the absorption response in the visible range. The broad optical absorption range with high absorption coefficients (10 5 cm −1 ) are comparable to the perovskites that are known to be highly efficient for solar cells [62]. The large absorption coefficients combining the significant optical anisotropy in the visible region can support its superior properties in terms of water-splitting reactions and optoelectronic applications.  To further explore the potential applications of J-In 2 XY MLs in optoelectronic devices based on their excellent optical properties, we construct an abrupt-change J-In 2 XY ML p-n junction phototransistor to investigate the photocurrent along the two main directions. Noticeably, the top view of the schematic diagram is illustrated in figure 7(a), and the side view along the y-direction is given in figure 7(b). The p-n junctions along the transport directions (x and y) with two-probed devices consisting of two semi-infinite electrodes sandwiched a central scattering region. The lengths of the central scattering region in the x-and y-directions are both about 7.5 nm. The source (S) and drain (D) electrodes are doped with p-and n-type concentrations of 10 19 e cm −3 , which similar doping level is commonly applied in 2D-based FET simulations [63,64] and experimentally possible [65,66]. An external bias voltage V ds = 1.0 V with both relative permittivity and relative permeability set to 1 is applied to drive the photo-generated carriers, and the photon-generated current under the linearly polarized light radiation is obtained by using AM1.5 standard solar spectrum with that of photon flux F ph = 1 s −1 Å −2 in the 300 K. The local density of states (LDOS) of J-In 2 XY p-n junctions along the y-direction are illustrated in figure 7(c), A typical p-n profile can be clearly seen along the device, and there is a 1.0 V potential difference between S and D electrodes. The width of the black area refers to the bandgap of the electrode material, which is very close to the results mentioned above. Furthermore, the flat band near the S/D electrode represents the converged screening built-in potential that provides a drift field for the photoexcited carrier. The direction-dependent photogenerated current density (I CD ) is calculated for testing the polarization sensitivity based on the electron-phonon coupling within the DFT +NEGF formalism. As shown in figures 7(d) and (e), the peaks of I CD for the p-n junctions are in the visible and near UV regions, the maxima are at the level of 10 −9 µA mm −1 in the x-direction and 10 −7 µA mm −1 in the y-direction. It manifests that the p-n junctions are prominent photoelectric anisotropy and result in a high on/off ratio, the measured ratio is ∼100 according to the formula [max(I CDx , I CDy )/ min(I CDx , I CDy )]. Furthermore, the I CD of J-In 2 SeTe ML along y-direction is 14.41 ×10 −7 µA mm −1 , which is higher than those of SbBi ML optoelectronic FET [67]. Note the frontier molecular orbitals of the electrode and the center scattering region do not produce resonance peaks at 4 and 5 eV, the photo-excited carriers are localized, resulting in very small photocurrents. To characterize the energy conversion efficiency, the photo responsivity (R ph ) is a critical indicator that is defined as: R ph = J ph P in A , where J ph and A refer to the photocurrent and effective illumination area, P in = eE ph F ph represents the incident photon power density. The obtained R ph of J-In 2 SeTe, J-In 2 STe and J-In 2 SSe MLs in the y-direction are 0.038, 0.036 and 0.015 AW −1 , respectively. Which is close to that SbBi (0.089 AW −1 ) [67], black phosphorus (0.0648 AW −1 ) [68], and higher than those of MoS 2 (0.016 AW −1 ) [66]. Moreover, the external quantum efficiency (EQE) is EQE = R ph hc eλ , where h and c are the Planck's constant and the speed of light, e and λ denote the electron charge and the wavelength of incident light, respectively. The calculated peaks of EQEs in J-In 2 SeTe, J-In 2 STe and J-In 2 SSe MLs are 10.6%, 11.5% and 4.7%. Therefore, the results of high R ph and EQE combined anisotropic photoelectric response further confirm that J-In 2 XY MLs can be suitable for the potential applications in multifunctional optoelectronic devices.

Conclusions
In summary, we investigate J-In 2 XY MLs with low symmetry as a promising multifunctional candidate for overall photocatalytic water splitting and flexible optoelectronic nanodevices. The anisotropy of electronic, piezoelectric, transport, and photoelectric properties are systematically explored by using first-principles calculations with the framework of DFT + NEGF. The favorable stability of J-In 2 XY MLs are confirmed by phonon dispersions, cohesive energies and AIMD simulations. Three types of stable structures are both flexible semiconductors with moderate bandgaps. The positions of valence and conduction-band edges are adequate for OER and HER over a wide PH range (from 0 to 7), owing to internal electric field with the direction of the Y layer pointing to the X layer caused by the out-of-plane polarization. The J-In 2 XY MLs with large out-of-plane polarization could be conducive to facilitating the efficient separation and transfer of photogenerated carriers. Moreover, the out-of-plane piezoelectric coefficient of J-In 2 SeTe ML (d 31 = 0.35 pm V −1 ) is larger by one order of magnitude compared to the J-TMDs. More importantly, highly anisotropic carrier mobility (more than 10 3 cm −2 V −1 s −1 along the x-direction, and 10 2 cm −2 V −1 s −1 in the y-direction) and visible light absorption (approximately 10 5 cm −1 ) herald notable performance of the overall photocatalysis. Moreover, the direction-dependent p-n junction based on J-In 2 XY MLs exhibits considerable photocurrent on/off ratios (∼100). Our results of anisotropic piezoelectric, transport and photoelectric together with suitable band edges for the novel J-In 2 XY MLs suggest that they exhibit attractive prospects for the design of multifunctional optoelectronic devices.

Data availability statement
All data that support the findings of this study are included within the article (and any supplementary files).