Abstract
Ferroelectric tunnel junctions (FTJs) have great potential for application in high-density non-volatile memories. Recently, α-In2Se3 was found to exhibit robust in-plane and out-of-plane ferroelectric polarizations at a monolayer thickness, which is ideal to serve as a ferroelectric component in miniaturized electronic devices. In this work, we design two-dimensional van der Waals heterostructures composed of an α-In2Se3 ferroelectric and a hexagonal IV–VI semiconductor and propose an in-plane FTJ based on these heterostructures. Our first-principles calculations show that the electronic band structure of the designed heterostructures can be switched between insulating and metallic states by ferroelectric polarization. We demonstrate that the in-plane FTJ exhibits two distinct transport regimes, tunneling and metallic, for OFF and ON states, respectively, leading to a giant tunneling electroresistance effect with the OFF/ON resistance ratio exceeding 1 × 104. Our results provide a promising approach for the high-density ferroelectric memory based on the 2D ferroelectric/semiconductor heterostructures.
Similar content being viewed by others
Introduction
Ferroelectric materials are useful for information storage applications due to their spontaneous polarization which can be switched by an applied electric field. Driven by the demands of device miniaturization and non-volatility, a ferroelectric tunnel junction (FTJ) has been proposed and considered as one of the most promising memory types, as compared to the current commercial dynamic random-access memories1,2,3,4,5. A typical FTJ is composed of two different electrodes separated by an ultrathin ferroelectric barrier layer6,7,8,9. Due to the two electrodes and/or interfaces in a FTJs being asymmetric, the reversal of ferroelectric polarization leads to the variation of the height and/or width of the electron tunneling potential barrier, resulting in two distinct resistance states7,9,10,11,12. This phenomenon is known as the tunneling electroresistance (TER) effect12,13,14. So far, most FTJs utilized oxide perovskite ferroelectrics as barrier layers. Unfortunately, ferroelectricity of perovskite ferroelectric films often deteriorate at the nanoscale due to the well-known size effect15, making the miniaturization of the memory devices challenging.
In recent years, ferroelectricity was found in a few two-dimensional (2D) van der Waals (vdW) materials, such as CuInP2S616,17, WTe218,19, and In2Se320,21,22. 2D vdW ferroelectrics are appealing due to their uniform atomic thickness, absence of dangling bonds, and the ability to be integrated with other vdW materials which allows useful functionalities. For example, in combination with conducting 2D materials, these 2D vdW ferroelectric materials can be used as ultrathin barriers in FTJs. Here we focus on α-In2Se3 which exhibits robust ferroelectricity even at a monolayer thickness and possesses intrinsically coupled out-of-plane (OOP) and in-plane (IP) ferroelectric polarizations20,23. Its sizable polarization and a moderate ferroelectric switching energy barrier make α-In2Se3 attractive for memory devices, especially FTJs24.
Since the discovery of ferroelectricity in 2D α-In2Se3, plenty of experimental and theoretical works have been devoted to exploring its application in non-volatile memories25,26,27,28,29,30,31,32,33,34. Wan et al. demonstrated a ferroelectric field-effect transistor (FeFET) using ultrathin α-In2Se3 as a non-volatile ferroelectric gate and graphene as a conducting channel, which were separated by an insulating hexagonal boron nitride layer35. Si et al. reported an asymmetric metal/α-In2Se3/Si crossbar ferroelectric semiconductor junction with an ON/OFF ratio >104 at room temperature36. Also, several theoretic works have been published focusing on the band alignment of 2D α-In2Se3 heterostructures controlled by the OOP polarization of α-In2Se3, which leads to about 1.1 eV potential difference between its two surfaces37,38,39. Ferroelectric control of the band alignment is important for the design of FTJs and other types of memory devices. Wang et al. summarized different categories of the band alignment between ferroelectric α-In2Se3 and 2D semiconductors or metals38.
A stronger polarization control of the band alignment in ferroelectric α-In2Se3 heterostructures and a simpler cell structure are still in focus for the memory application of 2D ferroelectrics. The recently proposed IP FTJs are especially interesting for the high-density information storage40. Using this concept, Kang et al.41,42,43 explored IP 2D FTJs based on graphene/α-In2Se3 vdW heterostructures, while Ding et al.44 proposed an antiferroelectric-like behavior of bilayer In2X3 (X = S, Se, Te) to enhance the FTJ performance. Relevant to the above studies, Shen et al.45 designed a 2D FTJ by p- and n-type doping at the two terminals of IP ferroelectric SnSe.
In this paper, we design a symmetric FTJs where both the ferroelectric barrier and the electrodes consist of layered vdW heterostructures. In these FTJs, the barrier layer represents an α-In2Se3 ferroelectric monolayer deposited on semiconductor monolayers, SnTe or PbSe, of hexagonal group IV–VI46,47, while the electrodes are composed of a heavily electron-doped In2Se3/SnTe bilayer, (In0.5Sn0.5)2Se3/(Sn0.5Sb0.5)Te. Using first-principles density functional theory (DFT) calculations, we demonstrate a robust ferroelectric polarization-controlled switching behavior between metallic and insulating states in the α-In2Se3/SnTe (and α-In2Se3/PbSe) ferroelectric heterostructures. By utilizing the coupled IP and OOP ferroelectricity of α-In2Se3 and the ferroelectric-controlled band alignment of the designed vdW heterostructures, we propose a two-terminal IP 2D FTJ with a giant TER effect (OFF/ON resistance ratio >1.2 × 104) and low ON state resistance. We infer the physical mechanism for the control of electron tunneling barriers and resulting transport properties for the ON and OFF states to be different and more effective as compared to those in the previous studies of IP FTJs. These results enrich the understanding of the ferroelectric band control in α-In2Se3 heterostructures and indicate a great potential for the 2D or ultrathin ferroelectrics for application in electronic devices.
Results and discussion
2D α-In2Se3/semiconductor vdW heterostructures
We consider heterostructures of 2D ferroelectric α-In2Se3 interfaced with 2D semiconductors SnTe and PbSe. The three-fold rotational symmetry of all these 2D materials results in the same IP triangular cell shapes, as shown in Fig. 1a–c. The optimized IP lattice constants of monolayers α-In2Se3, SnTe and PbSe (4.110, 4.172, and 4.096 Å, respectively) agree well with the existing results20,46,47,48,49, leading to tolerable lattice mismatches of 1.4% (In2Se3/SnTe) and 0.05% (In2Se3/PbSe). Band gaps of monolayers α-In2Se3, SnTe, and PbSe are calculated to be 0.78, 1.87, and 1.83 eV, respectively, which are in agreement with the existing results20,46,47,48,49 (see Supplementary Fig. 1 for their band structures).
Both Sn and Te terminations of SnTe at the interface with In2Se3 are considered in the heterostructures. The two terminations of SnTe in conjunction with two different OOP polarization directions of α-In2Se3 make four types of In2Se3/SnTe heterostructures, as shown in Fig. 1d–g. The atomic configurations of the heterostructures are determined by translating the SnTe layer along the interface plane relative to the In2Se3 layer and fully relaxing the IP lattice parameters and atomic positions. Based on the structure optimization, it is found that each heterostructure has the global energy minimum with the top layer atoms of SnTe stacked directly under the top layer In atoms of In2Se3. The IP lattice parameters for the four heterostructures shown in Fig. 1d–g are 4.113, 4.164, 4.080, and 4.086 Å, respectively. The distances between the top atomic plane of SnTe and the bottom Se atomic plane of In2Se3 for the four heterostructures are 2.822, 2.085, 2.976, and 2.815 Å, respectively.
The total energies of the four heterostructures (Fig. 1d–g) are calculated to be −1.971, −2.127, −1.947, and −1.978 eV per unit cell, respectively (the total energy of separate In2Se3 and SnTe monolayers is set to be 0). It is seen that, the P↑ state (where polarization of In2Se3 pointing away from SnTe) has 156 meV lower energy than the P↓ state (where polarization of In2Se3 pointing into SnTe) for the Sn–Se interface heterostructure. This energy difference is much larger than the energy difference of 31 meV between the two polarization states for the Te–Se interface heterostructure. The larger energy difference between the two polarization states in the Sn–Se interface heterostructure relative to that in the Te–Se interface heterostructure originates from the interfacial Coulomb interaction. Due to the ionic bonding between the Sn cation and Se anion at the Sn–Se interface which reduces the electrostatic energy, this type of heterostructure is sensitive to the polarization direction of In2Se3. When the polarization points away from the interface (P↑ state), a negative polarization charge of the P↑ state further reduces the electrostatic energy and leads to a smaller interlayer distance and an interfacial charge transfer (discussed below). On the contrary, due to the Coulomb repulsion between the Te and Se anions, the interfacial interaction at the Te–Se interface is much weaker. As a result, ferroelectric polarization reversal induced changes in the energy, interlayer distance, and interfacial charge transfer are insignificant. The obtained polarization states are stable and switchable as follows from the polarization switching study discussed in Supplementary Information (Supplementary Note 2 and Supplementary Fig. 2). Based on the same computational methods, In2Se3/PbSe heterostructures exhibit similar energetics and structural properties (Supplementary Note 3 and Supplementary Fig. 3) to those of In2Se3/SnTe heterostructures.
Next, we investigate the band alignment in the In2Se3/SnTe heterostructures. Figure 2 displays the calculated band structures of the four heterostructures shown in Fig. 1d–g. For the P↓ state of the Sn–Se interface heterostructure, the valence band maximum (VBM) and the conduction band minimum (CBM) are mainly contributed by SnTe and In2Se3, respectively, resulting in an indirect band gap of 0.84 eV (Fig. 2a). If the polarization is reversed upward (P↑), the heterostructure becomes metallic (Fig. 2b) with the VBM (mainly contributed by SnTe) being 0.39 eV higher in energy than the CBM (mainly contributed by In2Se3). This transition to the metallic state is caused by the negative ferroelectric polarization charge at the In2Se3/SnTe interface lifting the electrostatic potential energy of SnTe up and resulting in an electron charge transferred from the SnTe valence band to the In2Se3 conduction band. This fact is evident from the atomic layer-resolved density of states (LDOS) and the charge transfer at the interface as illustrated in Supplementary Figs. 4 and 5. The predicted transition between insulating and metallic states for the Sn–Se interface heterostructure driven by polarization reversal from P↓ to P↑ is consistent with the type II to type III band alignment transition as defined in ref. 38. The ferroelectric band control effect predicted in our study appears to be much more pronounced than those found in the similar 2D α-In2Se3/semiconductor heterostructures37,38,50.
The Te–Se interface heterostructure behaves differently. In this case, the negative ionic charge of the Te atom at the interface bends the SnTe bands down thus pushing the In2Se3 up into the band gap of SnTe. As a result, when polarization of In2Se3 is reversed, the change in the polarization charge at the interface is not sufficient to make the heterostructure metallic. We find that the P↓ state exhibits an indirect band gap of 1.08 eV (Fig. 2c) with both VBM and CBM derived mainly from the In2Se3 layer, whereas the P↑ state behaves an indirect band gap of 0.28 eV (Fig. 2d), where the VBM and CBM are contributed mainly by the SnTe and In2Se3, respectively. This transition is consistent with the type I to type II band alignment transition as defined in ref. 38.
The qualitative picture of the band alignment at the In2Se3/SnTe interface, being dependent on ferroelectric polarization of In2Se3 and the interface termination, is described in Supplementary Note 6. The In2Se3/PbSe heterostructures show similar electronic properties (Supplementary Note 7) to those of the In2Se3/SnTe heterostructures.
In-plane ferroelectric tunnel junctions
The significant polarization-induced band tuning effect and the resulting transition between insulating and metallic states, makes the α-In2Se3 ferroelectric heterostructures promising for application in memory and logic devices, such as FeFET and FTJs. In this work, we design an IP FTJ based on the Sn–Se terminated interface In2Se3/SnTe heterostructure. The transport direction of the FTJ is set along the IP polarization direction of α-In2Se3 (see Fig. 1a, b), which allows the simultaneous switching of OOP and IP polarizations through the application of an IP external electric field between the two electrodes. Figure 3 shows the atomic structures of the IP FTJs for the P↓ and P↑ states, respectively. The tunneling barrier in the FTJs consists of integer multiples (N) of the orthorhombic cells (denoted by the green rectangles in Figs. 1a and 3a) stacked along the transport direction. N is set to be ranging from 3 to 6. The lattice constant along the transport direction of each orthorhombic cell is 0.71 nm. Heavily electron-doped In2Se3/SnTe heterostructures, namely (In0.5Sn0.5)2Se3/(Sn0.5Sb0.5)Te, are used as electrodes in the FTJs. Their electronic properties are obtained using the VCA51 by mixing 50% Sn on the In site and 50% Sb on the Sn site (see Supplementary Note 8 for the band structures of the electrodes).
Next, we explore the electron transport and the TER effect in In2Se3/SnTe IP FTJs. Supplementary Fig. 9 shows the calculated transmission (T) as a function of electron energy (E) for the FTJs with N = 5. It should be pointed out that the electronic structures of the electrodes for the P↓ and P↑ states are not the same due to the electrodes being constructed using the VCA from different atomic structure corresponding to P↓ and P↑ (Supplementary Fig. 8). To eliminate the effect of electrodes and obtain the genuine TER effect of the FTJ, the transmission of the pure electrodes is calculated and shown in Supplementary Fig. 10. Then, the resistance of the FTJ is calculated according the formula: R = d(h/e2)[1/TFTJ − 1/Tel]52, where TFTJ and Tel are the calculated transmissions of the FTJ (Supplementary Fig. 9) and electrodes (Supplementary Fig. 10), respectively, and d is the lattice constant perpendicular to the transport direction (0.41 nm). The calculated resistance (R) as a function of electron energy (E) of the FTJ with N = 5 for the P↓ and P↑ states is plotted in Fig. 4a. It is seen that, for the P↓ state, the resistance increases significantly when the electron energy increases from EF −0.6 eV to EF −0.3 eV, and stays high until the energy increases up to EF (EF stands for the Fermi energy). Then the transmission gradually drops from EF to EF +0.2 eV again as the energy keeps increasing. On the contrary, for the P↑ state, the FTJ maintains a low resistance state under the considered electron energy range. The TER ratio (OFF/ON resistance ratio R(P↓)/R(P↑)) of the FTJ at E = EF is estimated to be 2.85 × 103. For E = EF − 0.3 eV, the OFF/ON ratio is found to be as large as 1.28 × 104.
Then, the barrier width dependence of the TER effect is studied by calculating the resistance of FTJs at E = EF with the barrier layer width varying from 3 to 6 orthorhombic unit cells. It is seen from Fig. 4b that the resistance of the ON state (P↑ state) stays low and is enhanced slowly when the barrier width increases. This is because, in the P↑ state, the heterostructure is metallic, allowing the FTJ to behave like a conductor. On the contrary, for the OFF state (P↓ state), the calculated resistance increases exponentially with increasing the barrier width, exhibiting a typical tunneling behavior. As a result, the OFF/ON resistance ratios are calculated to be 1.51 × 103, 2.76 × 103, 2.85 × 103, and 1.21 × 104 for the barrier layer width varying from 3 to 6 orthorhombic cells (Fig. 4b). We conclude, therefore, that the ferroelectric polarization induced TER effect in the designed In2Se3/SnTe IP FTJ can be enhanced significantly by increasing the barrier width.
In order to understand the origin of this transport behavior, the electronic structure of the FTJ is analyzed. Figure 5 shows the calculated layer-resolved density of states (LDOS) along the width direction. Each panel in Fig. 5 displays the DOS projected on half an orthorhombic cell in the electrode and barrier regions for P↓ (Fig. 5a) and P↑ (Fig. 5b) states, respectively. As seen from Fig. 5a, the barrier layer exhibits a band gap from about −0.6 eV to 0.2 eV relative to EF. This feature of the band structure well illustrates why the resistance of the P↓ state increases and decreases at the energy of EF −0.6 eV and EF +0.2 eV, respectively, as shown in Fig. 4a (blue line). On the contrary, for the P↑ state, as shown in Fig. 5b, there is no band gap across the barrier layer, leading to good conductivity under the entire considered energy range (see Fig. 4a, red line). The band structure-controlled insulating and conducting properties of the barrier layer for the P↓ and P↑ states also explain the barrier width dependence of the transport behavior shown in Fig. 4b.
To further clarify the transport mechanism of the designed IP FTJ, the partial charge densities in the real space are analyzed. Figure 6 visualizes the partial charge densities in the range of energies from EF −0.1 eV to EF for P↓ and P↑ states, respectively. It is seen from Fig. 6a that, for the P↓ state, the partial charge density at EF only exists in the electrodes and a few atomic layers of the barrier near the interfaces, agreeing well with the LDOS of the FTJ at the Fermi energy in Fig. 5a. On the contrary, for the P↑ state, as seen from Fig. 6b, there is abundant charge distributed within the barrier region, corresponding to the LDOS in Fig. 5b. From the spatial distribution, the partial charge density is mainly concentrated on the SnTe layer and the In2Se3/SnTe interface in the barrier region, which indicates that the electrons will transport along this channel as in a conductor under the external bias.
We would like to emphasize that the mechanism of TER considered in this work is different from those discussed previously, where asymmetry of an FTJ, i.e., different electrodes and/or different interfaces, was the necessary factor to obtain a large TER value. Here, a metal-insulator transition in the whole barrier layer is driven by the OOP polarization of In2Se3 and the FTJ asymmetry is not required. Contrary to the previous studies, the ferroelectric polarization that controls TER is pointing perpendicular to the transport direction, and the type of electrodes is not essential. Such design of the FTJ may be beneficial for applications because ON/OFF switching can be achieved in two ways: (1) by applying a bias voltage between the two electrodes to reverse the OOP polarization via reversal of the IP polarization as they are intrinsically coupled in ferroelectric In2Se3, or (2) by applying a bias voltage to the top gate electrode to directly switch the OOP polarization as in a three-terminal FeFET.
In conclusion, we have designed vdW heterostructures with α-In2Se3 ferroelectric and hexagonal group IV–VI 2D semiconductor layers (SnTe and PbSe). Based on systematic first-principles studies, we demonstrated that these heterostructures exhibit a pronounced ferroelectric polarization control of the band alignment, which is promising for electronic device applications. We predicted that the In2Se3/SnTe heterostructure with Sn–Se interface termination can be switched from the insulating state with a large band gap to the metallic state, when the OOP ferroelectric polarization of α-In2Se3 is reversed from pointing into to pointing away from 2D semiconductor layer. Driven by the polarization-controlled band alignment and the coupled IP and OOP ferroelectric polarizations in the designed In2Se3/SnTe heterostructures, we proposed a functional two-terminal IP 2D FTJ. We predicted that the designed FTJ has two distinct electron transport mechanisms−tunneling and metallic for the OFF and ON states, respectively, resulting in a giant TER effect with the OFF/ON resistance ratio exceeding 1 × 104. The predicted TER effect is significantly amplified by further increasing the barrier width of the FTJ. Our results are promising for the design of 2D ferroelectric heterostructures and FTJ-based non-volatile memories. We hope therefore that they will stimulate experimental efforts to explore such kind of FTJs in practice and demonstrate the predicted properties.
Methods
Geometry optimization and electric structure calculations
DFT calculations using the projector augmented-wave as implemented in VASP code53 are preformed to optimize the monolayer 2D materials and heterostructures and study their electronic structures. The exchange-correlation effects are described by the Perdew-Burke-Ernzerhof generalized gradient approximation54. The plane-wave cutoff energy is set to be 500 eV. The 12 × 12 × 1 and 21 × 21 × 1 Monkhorst-Pack k-point meshes for the irreducible Brillouin zone sampling are used in the geometric optimization and subsequent self-consistent calculation of the heterostructures. Vacuum layers thicker than 15 Å are set along the normal direction to avoid the interaction between the top and bottom surfaces. The vdW interactions are treated using the DFT-D3 corrections55. The IP lattice parameters and atomic positions are fully relaxed until the Hellmann-Feynman force on each atom is less than 0.005 eV/Å.
Electron transport calculations
The IP FTJs are constructed along the IP polarization direction based on the relaxed In2Se3/SnTe heterostructures. The barrier region consists of several unit cells of In2Se3/SnTe, and the electrodes are treated using the virtual crystal approximation51 where the In site of In2Se3 is occupied by In0.5Sn0.5 and the Sn site of SnTe is occupied by Sn0.5Sb0.5 thus simulating the n-type doping. The electron transmissions across the FTJ are calculated within the general scattering formalism56 as implemented in Quantum ESPRESSO57. The designed IP FTJs and the (In0.5Sn0.5)2Se3/(Sn0.5Sb0.5)Te heterostructures are used as the scattering regions and semi-infinite leads, respectively. A uniform 50 × 10 k||-mesh is used to sample the two-dimension Brillouin zone in transmission calculation.
Data availability
The data that support the findings of the work is in the manuscript’s main text and Supplementary Information. Additional data are available from the corresponding author upon reasonable request.
Code availability
The central codes used in this paper are VASP and Quantum ESPRESSO.
References
Esaki, A. L., Laibowitz, R. & Stiles, P. Polar switch. IBM Tech. Discl. Bull. 13, 114 (1971).
Tsymbal, E. Y. & Kohlstedt, H. Tunneling across a ferroelectric. Science 313, 181–183 (2006).
Kohlstedt, H., Pertsev, N. A., Contreras, J. R. & Waser, R. Theoretical current-voltage characteristics of ferroelectric tunnel junctions. Phys. Rev. B 72, 125341 (2005).
Kim, D. J. et al. Ferroelectric tunnel memristor. Nano Lett. 12, 5697–5702 (2012).
Garcia, V. & Bibes, M. Ferroelectric tunnel junctions for information storage and processing. Nat. Commun. 5, 4289 (2014).
Xi, Z. et al. Giant tunnelling electroresistance in metal/ferroelectric/semiconductor tunnel junctions by engineering the Schottky barrier. Nat. Commun. 8, 15217 (2017).
Wen, Z., Li, C., Wu, D., Li, A. & Ming, N. Ferroelectric-field-effect-enhanced electroresistance in metal/ferroelectric/semiconductor tunnel junctions. Nat. Mater. 12, 617–621 (2013).
Pantel, D. et al. Tunnel electroresistance in junctions with ultrathin ferroelectric Pb(Zr0.2Ti0.8)O3 barriers. Appl. Phys. Lett. 100, 232902 (2012).
Velev, J. P., Burton, J. D., Zhuravlev, M. Y. & Tsymbal, E. Y. Predictive modelling of ferroelectric tunnel junctions. NPJ Comput. Mater. 2, 16009 (2016).
Liu, X., Burton, J. D. & Tsymbal, E. Y. Enhanced tunneling electroresistance in ferroelectric tunnel junctions due to the reversible metallization of the barrier. Phys. Rev. Lett. 116, 197602 (2016).
Wu, Q. et al. Giant tunneling electroresistance induced by ferroelectrically switchable two-dimensional electron gas at nonpolar BaTiO3/SrTiO3 interface. Phys. Rev. B 94, 155420 (2016).
Yang, Q. et al. Ferroelectric tunnel junctions enhanced by a polar oxide barrier layer. Nano Lett. 19, 7385–7393 (2019).
Zhuravlev, M. Y., Sabirianov, R. F., Jaswal, S. S. & Tsymbal, E. Y. Giant electroresistance in ferroelectric tunnel junctions. Phys. Rev. Lett. 94, 246802 (2005).
Zhuravlev, M. Y., Maekawa, S. & Tsymbal, E. Y. Effect of spin-dependent screening on tunneling electroresistance and tunneling magnetoresistance in multiferroic tunnel junctions. Phys. Rev. B 81, 104419 (2010).
Junquera, J. & Ghosez, P. Critical thickness for ferroelectricity in perovskite ultrathin films. Nature 422, 506–509 (2003).
Liu, F. et al. Room-temperature ferroelectricity in CuInP2S6 ultrathin flakes. Nat. Commun. 7, 12357 (2016).
Brehm, J. A. et al. Tunable quadruple-well ferroelectric van der Waals crystals. Nat. Mater. 19, 43–48 (2020).
Sharma, P. et al. A room-temperature ferroelectric semimetal. Sci. Adv. 5, eaax5080 (2019).
Fei, Z. et al. Ferroelectric switching of a two-dimensional metal. Nature 560, 336–339 (2018).
Ding, W. et al. Prediction of intrinsic two-dimensional ferroelectrics in In2Se3 and other III2-VI3 van der Waals materials. Nat. Commun. 8, 14956 (2017).
Xiao, J. et al. Intrinsic two-dimensional ferroelectricity with dipole locking. Phys. Rev. Lett. 120, 227601 (2018).
Zheng, C. et al. Room temperature in-plane ferroelectricity in van der Waals In2Se3. Sci. Adv. 4, eaar7720 (2018).
Cui, C. et al. Intercorrelated in-plane and out-of-plane ferroelectricity in ultrathin two-dimensional layered semiconductor In2Se3. Nano Lett. 18, 1253–1258 (2018).
Su, Y. et al. Van der Waals multiferroic tunnel junctions. Nano Lett. 21, 175–181 (2021).
Xue, F. et al. Gate-tunable and multidirection-switchable memristive phenomena in a van der Waals ferroelectric. Adv. Mater. 31, 1901300 (2019).
Xie, S. et al. Ferroelectric semiconductor junctions based on graphene/In2Se3/graphene van der Waals heterostructures. 2D Mater. 8, 045020 (2021).
Wang, S. et al. Two-dimensional ferroelectric channel transistors integrating ultra-fast memory and neural computing. Nat. Commun. 12, 53 (2021).
Su, H., Hu, T., Wu, F. & Kan, E. Controllable vdW contacts between the ferroelectric In2Se3 monolayer and two-dimensional metals. J. Phys. Chem. C. 125, 10738–10746 (2021).
Mukherjee, S. & Koren, E. Indium selenide (In2Se3)—an emerging van-der-Waals material for photodetection and non-volatile memory applications. Isr. J. Chem. 62, e202100112 (2022).
Hu, H. & Ouyang, G. First-principles calculations of interface engineering for 2D α-In2Se3-based van der Waals multiferroic heterojunctions. Appl. Surf. Sci. 545, 149024 (2021).
Cheng, H.-X., Zhou, J., Wang, C., Ji, W. & Zhang, Y.-N. Nonvolatile electric field control of magnetism in bilayer CrI3 on monolayer In2Se3. Phys. Rev. B 104, 064443 (2021).
Wan, S. et al. Room-temperature ferroelectricity and a switchable diode effect in two-dimensional α-In2Se3 thin layers. Nanoscale 10, 14885–14892 (2018).
Poh, S. M. et al. Molecular-beam epitaxy of two-dimensional In2Se3 and its giant electroresistance switching in ferroresistive memory junction. Nano Lett. 18, 6340–6346 (2018).
Li, Y. et al. Orthogonal electric control of the out-of-plane field-effect in 2D ferroelectric α-In2Se3. Adv. Electron. Mater. 6, 2000061 (2020).
Wan, S. et al. Nonvolatile ferroelectric memory effect in ultrathin α-In2Se3. Adv. Funct. Mater. 29, 1808606 (2019).
Si, M. et al. Asymmetric metal/α-In2Se3/Si crossbar ferroelectric semiconductor junction. ACS Nano 15, 5689–5695 (2021).
Zhou, B. et al. Ferroelectric and dipole control of band alignment in the two dimensional InTe/In2Se3 heterostructure. J. Phys. Condens. Matter 32, 055703 (2019).
Wang, Z. & Zhu, W. Tunable band alignments in 2D ferroelectric α-In2Se3 based van der Waals heterostructures. ACS Appl. Electron. Mater. 3, 5114–5123 (2021).
Ayadi, T. et al. An ab initio study of the ferroelectric In2Se3/graphene heterostructure. Phys. E Low. Dimens. Syst. Nanostruct. 114, 113582 (2019).
Shen, H., Liu, J., Chang, K. & Fu, L. In-plane ferroelectric tunnel junction. Phys. Rev. Appl. 11, 024048 (2019).
Kang, L. et al. Domain-wall induced giant tunneling electroresistance effect in two-dimensional graphene/In2Se3 ferroelectric tunnel junctions. Phys. E Low. Dimens. Syst. Nanostruct. 133, 114783 (2021).
Kang, L. et al. Giant tunnel electroresistance in ferroelectric tunnel junctions with metal contacts to two-dimensional ferroelectric materials. Phys. Rev. B 103, 125414 (2021).
Kang, L. et al. Giant tunneling electroresistance in two-dimensional ferroelectric tunnel junctions with out-of-plane ferroelectric polarization. Phys. Rev. B 101, 014105 (2020).
Ding, J., Shao, D.-F., Li, M., Wen, L.-W. & Tsymbal, E. Y. Two-dimensional antiferroelectric tunnel junction. Phys. Rev. Lett. 126, 057601 (2021).
Shen, X.-W., Fang, Y.-W., Tian, B.-B. & Duan, C.-G. Two-dimensional ferroelectric tunnel junction: the case of monolayer In:SnSe/SnSe/Sb:SnSe homostructure. ACS Appl. Electron. Mater. 1, 1133–1140 (2019).
Kamal, C., Chakrabarti, A. & Ezawa, M. Direct band gaps in group IV-VI monolayer materials: binary counterparts of phosphorene. Phys. Rev. B 93, 125428 (2016).
Tian, X.-Q. et al. Hexagonal layered group IV-VI semiconductors and derivatives: fresh blood of the 2D family. Nanoscale 12, 13450–13459 (2020).
Fatahi, N. et al. 2D Hexagonal SnTe monolayer: a quasi direct band gap semiconductor with strain sensitive electronic and optical properties. Eur. Phys. J. B 93, 32 (2020).
Allaoui, I., Benyoussef, A. & El Kenz, A. Two-dimensional SnTe/Sb van der Waals heterostructure for photovoltaic application. Solid State Sci. 121, 106736 (2021).
Ayadi, T. et al. An ab initio study of the electronic properties of the ferroelectric heterostructure In2Se3/Bi2Se3. Appl. Surf. Sci. 538, 148066 (2021).
Fang, Z., Solovyev, I. V. & Terakura, K. Phase diagram of tetragonal manganites. Phys. Rev. Lett. 84, 3169–3172 (2000).
Xia, K. et al. Interface resistance of disordered magnetic multilayers. Phys. Rev. B 63, 064407 (2001).
Kresse, G. & Joubert, D. From ultrasoft pseudopotentials to the projector augmented-wave method. Phys. Rev. B 59, 1758–1775 (1999).
Perdew, J. P., Burke, K. & Ernzerhof, M. Generalized gradient approximation made simple. Phys. Rev. Lett. 77, 3865–3868 (1996).
Grimme, S., Antony, J., Ehrlich, S. & Krieg, H. A consistent and accurate ab initio parametrization of density functional dispersion correction (DFT-D) for the 94 elements H-Pu. J. Chem. Phys. 132, 154104 (2010).
Smogunov, A., Dal Corso, A. & Tosatti, E. Ballistic conductance of magnetic Co and Ni nanowires with ultrasoft pseudopotentials. Phys. Rev. B 70, 045417 (2004).
Giannozzi, P. et al. QUANTUM ESPRESSO: a modular and open-source software project for quantum simulations of materials. J. Phys. Condens. Matter 21, 395502 (2009).
Acknowledgements
The authors thank Dr. Lingling Tao for helpful discussions. This work was supported by the National Natural Science Foundation of China (Grant Nos. 12072307 and 12175191), the Outstanding Youth Science Foundation of Hunan Province, China (Grant No. 2021JJ20041) and the Research Foundation of Education Bureau of Hunan Province, China (Grant No. 21A0114). The atomic structures were plotted using VESTA software.
Author information
Authors and Affiliations
Contributions
Q.Y. and J.J. designed the idea of the research. Z.L. performed the first-principles calculations and wrote the manuscript. All the authors contributed to the general discussion, review and editing. Q.Y. supervised the research.
Corresponding authors
Ethics declarations
Competing interests
The authors declare no competing interests.
Additional information
Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary information
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Liu, Z., Hou, P., Sun, L. et al. In-plane ferroelectric tunnel junctions based on 2D α-In2Se3/semiconductor heterostructures. npj Comput Mater 9, 6 (2023). https://doi.org/10.1038/s41524-022-00953-x
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41524-022-00953-x
This article is cited by
-
Tuning ferroelectric phase transition temperature by enantiomer fraction
Nature Communications (2024)