1.
Introduction

Owning the switchable and non-volatile electric polarizations at room temperature, conventional ferroelectrics, such as complex perovskite oxides, have broad applications in electronics, including high energy density capacitor^[1], tunable diode^[2], ferroelectric tunneling junction (FTJ)^[3], and ferroelectric field-effect transistors (FeFET)^[4]. With the further miniaturization of electronic devices beyond the Moore’s law, there is a need to explore ferroelectric thin films with nanometer thickness. However, due to the imperfect charge screening, local chemical environment, defects, and misfit strain at the interface, there is depolarization field arising in ultra-thin ferroelectrics, which leads to the so-called critical size effect^[5–9]. The electric polarizations of conventional ferroelectrics will reduce or even disappear when the thickness reaches to the single atomic layer limit. As a consequence, their application in modern electronics is strongly hindered. In parallel, since the rise of graphene, the family of two-dimensional (2D) materials has been greatly expanded with members owning numerous novel physical properties such as superconductivity^[10], charge density wave^[11], valley polarization^[12], and exotic magnetism^[13]. These van der Waals (vdW) materials have saturated interface chemical environment, weak interaction between the layers and the benefit of easily preparing stable atomically thin slabs, which offer a novel platform to study some emerging phenomena that have never been observed in their bulk form counterpart. Therefore, it offers a feasible way to explore ferroelectricity at atomic scale in 2D materials, especially for the out-of-plane electric polarizations, which is more important than the in-plane electric polarization from the technological perspective.



In the early days of the research on vdW materials, 2D ferroelectricity is lack of attention with rare findings. The reason is that when the thickness is reduced to nanometer scale, the competition between crystal structure stability and the long-range ordering of electric dipoles will become significant. Thus, for a long time, it is hard to get stable crystal phase with switchable electric polarizations in vdW materials under ambient conditions. As ferroelectrics, the prerequisite is the inversion symmetry breaking in structure. For most of the commonly studied 2D materials, taking graphene as a typical example^[14], the centrosymmetric lattice naturally rules out the possibility of being ferroelectric. Even in 2H-MoS₂^{[15, 16]}, where the inversion symmetry is explicitly broken in its monolayer structure, the presence of the mirror symmetry along its middle plane of metal Mo atoms makes it keeping symmetric Mo-S chemical bonds in the out-of-plane direction, resulting void electric dipoles. Given these facts, seeking 2D ferroelectricity in vdW materials is of challenges but of great values.



The 2D nature of vdW materials allows to explore ferroelectricity in both the in-plane and out-of-plane directions. The in-plane vdW ferroelectrics provide a platform to study domain physics at 2D limit, while the out-of-plane ferroelectric polarizations lay down the fundamental for device applications. In recent years, several candidates of 2D ferroelectrics have been theoretically predicated and verified experimentally. The first example is the layered CuInP₂S₆. In 2015, Belianinov et al.^[17] reported the observation of room-temperature out-of-plane ferroelectricity in bulk CuInP₂S₆ with film thickness more than 100 nm. The electric polarization stemmed from the off-center ordering of the Cu and In cations in each unit cell with its polar axis along the vertical direction relative to the 2D plane. Although the collinear two-sublattice structure of CuInP₂S₆ makes it a ferrielectric system, this finding also inspires further exploration on 2D ferroelectrics. Later in 2016, Liu et al.^[18] observed the out-of-plane electric polarization in few-layer CuInP₂S₆ with much thinner thickness of 4 nm under ambient conditions. Utilizing the second harmonic generation (SHG) technique, they clearly showed that the phase transition temperature (T_c) from the ferroelectric phase to the paraelectric phase is 320 K, which is slightly higher than the room temperature. This is crucial for their practical device application. At the same time, Chang et al.^[19] from Tsinghua University experimentally reported the first 2D ferroelectric with single atomic layer thickness. They studied the monolayer SnTe with distorted lattice at cryogenic temperature and verified the in-plane ferroelectricity with evidences such as spontaneous domains and electric polarization induced band bending. These pioneering findings open up the route towards 2D ferroelectricity with vdW materials. However, for layered SnTe, the nature as an in-plane ferroelectric strongly hinders its application in practical device application, while for CuInP₂S₆, it is ferrielectric, resulting relatively weak saturated electric polarizations. From the past researches on ferroelectricity, as a consensus, their practical application relies on the technically more important out-of-plane electric polarization, high T_c and simple lattice structure for material realization. Based on these requirements, very recently Ding et al.^[20] predicted the intrinsic room-temperature 2D ferroelectricity with intercorrelated in-plane and out-of-plane electric polarization in the 2D semiconducting α-In₂Se₃. Stimulated by this unprecedented inter-locking of electric dipoles in α-In₂Se₃, intense experimental efforts were devoted with various techniques^[21–28]. Clear ferroelectric domains, ferroelectric hysteresis loop and piezoelectricity have been observed in ultrathin samples with the thickness down to 2D limit^{[23, 26]}. Besides, several practical device applications based on the few layer α-In₂Se₃ including switchable ferroelectricdiode^[24], FeFETs^[29], and ferroelectric semiconductor FET (FeS-FET)^[30] were demonstrated. These breakthroughs point to the possibility of stable ferroelectric phase in layered vdW materials.



In view of the rapid growing research interest in the 2D ferroelectricity of α-In₂Se₃, we review some recent relevant progresses in this article. We begin with the introduction of the crystal structure of α-In₂Se₃ and the corresponding ferroelectricity. The experimental verifications of the in-plane and out-of-plane electric polarizations in ultra-thin α-In₂Se₃ will be highlighted. Finally, we present two prototypes of ferroelectric devices to demonstrate the potential device application of 2D ferroelectric α-In₂Se₃.

2.
Crystal structure of layered α-In₂Se₃

Layered In₂Se₃ is an available 2D semiconductor with remarkable optical and electrical properties^[31–33]. It has long been studied in the application of photo sensors and phase-change memories^[34–36]. Previous studies have found that the phases of In₂Se₃ are very rich. In addition to the α phase, there exists crystal structures with different lattice constants, which are labeled as β, γ or δ phase. In all these phases, the α-In₂Se₃ is recognized as the most stable layered structure at room temperature. As shown in Fig. 1, its single layer is composed of alternating Se or In atomic layer via covalent bonds, constituting a Se–In–Se–In–Se quintuple layer (QL). These QLs are vertically stacked via the weak van der Waals force. However, the identification of the crystal structure for α phase is highly controversial at early stage^{[31–35, 37–39]}, which is due to the small energy barrier in structure stability between the α and β phase.






onerror="this.onerror=null;this.src='http://www.jos.ac.cn/fileBDTXB/journal/article/jos/2019/6/PIC/19030026-1.jpg'"
class="figure_img" id="Figure1"/>



Download



Larger image


PowerPoint slide






Figure1.
(Color online) Layered structures of In₂Se₃. (a) Three-dimensional crystal structure of layered In₂Se₃, with the In atoms in blue and Se atoms in red, and a quintuple layer (QL) is indicated by the black dashed square. (b) Top view of the system along the vertical direction. Each atomic layer in a QL contains only one elemental species, with the atoms arranged in one of the triangular lattices A, B or C as illustrated. (c–h) Side views of several representative structures of one QL In₂Se₃, among which the (c) to (e) structures are derived from the zincblende, wurtzite and fcc crystals, respectively. In (f), the interlayer spacings between the central Se layer and the two neighboring In layers are displayed. The black arrows in (f) and (g) indicate the directions of the spontaneous electric polarization (P) in the FE-ZB' and FE-WZ' structures, respectively. The FE-ZB’ and FE-WZ’ structures in (f) and (g) are identified as α phase, while the metastable structure fcc and fcc’ in (e) and (h) are identified as β phase. (Courtesy of Ref. [20])




To clarify the crystal structure of In₂Se₃, Ding et al.^[20] has performed detailed numerical calculations with density function theory (DFT). They examined the possible atomic arrangements in layered In₂Se₃, and obtained two lattice structures with lowest total energy in thermodynamics. As shown in Figs. 1(f) and 1(g), the FE-ZB' and the FE-WZ' structures were identified as ground states, which belonged to α phase. From the top view of its crystal structure, there are three positions marked as A, B, and C (represent one corner of the triangular lattices in Fig. 1(b)) that can be occupied by Se or In atoms. Along the out-of-plane direction, the atoms of a given QL in FE-ZB' and FE-WZ' structure are arranged with ABBCA and ABBAC sequence respectively. For the ABCAB structure (Fig. 1(e)), it is identified as β phase, which is metastable at room temperature. With the confirmation of its crystal structure, they predicted the emergence of 2D ferroelectricity in α-In₂Se₃, which was attributed to the subtle atomic configuration of the central Se atoms. Taking the FE-ZB’ structure as an example, in each QL, the Se and In atoms in the lower three Se–In–Se atomic layers form a regular tetrahedral structure, while the in the upper Se–In–Se layers they form the octahedral structure. Thus, the Se atoms in the middle layer are in different chemical environments with respect to its neighboring layers. This breaks the central inversion symmetry thereby produces the electric polarizations along both the out-of-plane and in-plane direction. Moreover, in α-In₂Se₃, the in-plane and out-of-plane electric dipoles are inherently coupled, which not only provides a novel approach of electric polarization manipulation with external electric field but also helps the stabilization of the out-of-plane electric polarizations against the depolarization electric field in ultrathin samples. In contrast, the β phase of In₂Se₃ presented in Fig. 1(e) follows the ABCAB layer sequence with the central inversion symmetry, which forbids the spontaneous formation of net electric polarizations. Meanwhile, due to the structural phase transition, the phonon modes of α and β phases of In₂Se₃ were predicted to have distinct features, such as the softening of A1 phonon mode in the ferroelectric phase. These Raman modes can be utilized as an effective tool to identify the α-In₂Se₃ in following experimental studies, including the results in our group^{[21–22, 24, 25, 29]}.



As for ferroelectrics, one important signature is the switchable electric polarization with external stimulus such as electric field. For α-In₂Se₃, intuitively, the degenerate polarization states with opposite polarity can be achieved by in-plane 180° rotation of a single QL without changing the crystal structure and symmetry. It has been theoretically suggested that the transformation of these two polarization states is kinetically feasible^[20]. The flipping process is proposed as "three-step cooperated polarization inversion process", where the electron polarization reversion is achieved by the metastable fcc' structure (β phase) as the intermediate transition state. The energy barrier for polarization flipping is only 0.066 eV per unit cell, which is comparable with that in conventional ferroelectric material PbTiO₃ (0.07 eV/ unit cell). This feasible polarization flipping process in α-In₂Se₃ further evidence its ferroelectric nature.

3.
Experimental verification of coupled ferroelectricity in 2D α-In₂Se₃

Inspired by the previous theoretical prediction, several groups have aimed to experimentally verify the room temperature ferroelectricity in α-phase In₂Se₃^{[21–26, 28]}. Ferroelectrics naturally have the piezoelectric effect, in which deformation of the crystals can be realized by an external electric field. Therefore, the detection of ferroelectricity in vdW materials can be realized via surface sensitive techniques such as Piezo force microscopy (PFM). In most of recent studies on the 2D ferroelectricity in α-In₂Se₃, PFM was chosen as the versatile tool. Fig. 2 shows the typical PFM study on ultra-thin α-In₂Se₃ with thickness about 20 nm under ambient condition^[24]. The samples were mechanically exfoliated from bulk and post transferred to conductive substrate such as heavily doped Si or Au. Clear out-of-plane spontaneous electric polarizations were visualized through the contrasted PFM phase image as shown in Fig. 2(b). The phase difference between adjacent ferroelectric domains was found to be 180°, indicating that there were antiparallel electric dipoles (Fig. 2(c)) in the out-of-plane direction. In the single poling measurement, the butterfly-like ferroelectric hysteresis is observed in the out-of-plane PFM amplitude (see Fig. 2(d)), and the transition between the two different ferroelectric states can also be driven by an external electric field, with two different critical field strengths during forward and backward sweeping (see Fig. 2(e)). These results show that the electric polarization could be flipped under external electric field. The coercive field (E_c) of α-In₂Se₃ was found to be about 200 kV/cm^[24], which was much lower than that of layered ferroelectric CuInP₂S₆ (700 kV/cm)^[18]. The relatively low E_c benefits the low operating power ferroelectric device based on α-In₂Se₃. Other important features of ferroelectrics include the domain engineering and the retention of electric polarizations. In order to verify these features, artificially creating the ferroelectric domain structure with the box-in-box pattern (Fig. 2(f)) has been tested in α-In₂Se₃. In detail, a square area of ??1 × 1 μm² is controllably written with positive bias via the PFM tip to establish the downward polarization. After that, the negative sample/tip bias is applied in the center 0.5 × 0.5 μm² area to reverse or erase the polarization. As a result of the ferroelectric nature of α-In₂Se₃, an artificial domain structure with the box-in-box pattern is created. It is found that artificial ferroelectric domain could sustain more than 24 h under ambient condition^[24].






onerror="this.onerror=null;this.src='http://www.jos.ac.cn/fileBDTXB/journal/article/jos/2019/6/PIC/19030026-2.jpg'"
class="figure_img" id="Figure2"/>



Download



Larger image


PowerPoint slide






Figure2.
(Color online) Ferroelectricity of α-In₂Se₃ thin layers. (a) The surface topography of α-In₂Se₃ thin layers (~20 nm) on the heavily doped Si substrate. The scale bar is 2 μm. (b) The corresponding PFM phase image in the out-of-plane direction, showing clear ferroelectric domains. (c) The phase profile of different ferroelectric domains as sketched by the red dashed line in (b). A phase contrast of 180° is observed, which indicates the antiparallel directions of out-of-plane electric polarization between the adjacent domains. The arrows indicate the directions of electric polarization. (d) PFM amplitude and (e) PFM phase hysteresis loop measured from α-In₂Se₃ thin layers. (f) PFM amplitude image of domain engineering in α-In₂Se₃ with a film thickness of 12 nm. The scale bar is 370 nm. (Courtesy of Ref. [24])




Another method to study the 2D ferroelectricity in vdW materials is the optical SHG technique, which is a direct probe of the lattice inversion symmetry breaking. With this technique, Xiao et al.^[25] has verified the out-of-plane and in-plane electric dipole inter-locking in α-In₂Se₃, which is consistent with previous PFM study^[22]. The rotation of in-plane electric dipole could be achieved by an out-of-plane electric field, which demonstrates the inherent locking of these two orthogonal polarizations. As shown in Fig. 3, in a 3-nm-thick sample with confirmed electric polarization, they first artificially construct the box-in-box domain in the out-of-plane direction with the PFM tip. In their SHG measurements, the normal incidence set-up is chosen, which is only sensitive to the in-plane electric polarization. By comparing the SHG mappings before and after the out-of-plane bias poling (Figs. 3(c) and 3(d)), clear dark lines (less SHG intensity) correspond to the boundaries of the patterned area by PFM tip bias. This observation can be understood as the destructive optical interference arising from the reversal of the in-plane nonlinear optical polarization and corresponding in-plane lattice asymmetry. These findings clearly demonstrate the interlocking of the out-of-plane and in-plane electric dipoles in 2D α-In₂Se₃. Moreover, the temperature-dependent SHG is also studied, which show that the SHG intensity is robust and is almost independent of temperature in a wide range of temperature range up to 500 K. The transition from the ferroelectric phase to paraelectric phase happens at about 700 K, beyond which a sharp decrease of SHG signal can be found due to the recovery of structural inversion symmetry. These results confirm that the T_c for ferroelectric α-In₂Se₃ is much higher than that of the other 2D ferroelectrics such as SnTe and CuInP₂S₆, laying the fundamental basis for its practical device application.






onerror="this.onerror=null;this.src='http://www.jos.ac.cn/fileBDTXB/journal/article/jos/2019/6/PIC/19030026-3.jpg'"
class="figure_img" id="Figure3"/>



Download



Larger image


PowerPoint slide






Figure3.
(Color online) Electrically switching the out-of-plane ferroelectric polarization and corresponding in-plane atomic configuration through dipole locking. (a) The hysteresis of remnant out-of-plane polarization of a 3-nm-thick In₂Se₃ crystal on conductive SrRuO₃, as a function of perpendicular poling voltage. Black, red, blue curves represent the normalized piezoresponse measured with V_ac = 0.5, 1, and 1.5 V, respectively. The collapse of the hysteresis loop is similar to the behavior of the conventional field-switchable ferroelectrics, but totally different from charging artifact of dielectrics. (b) Polarized domain patterned by electrically biased scanning probe and measured by PFM. The inner box corresponds to positive applied voltage (+6 V) with positive piezoresponse while the outer box to negative voltage (?6 V) with negative piezoresponse. (c) SHG intensity mapping on another trilayer In₂Se₃ sample before PFM reversed poling. The area enclosed by dashed line was then scanning by a negatively biased AFM tip. The color bar is in linear scale with arbitrary unit. (d) SHG mapping after the electrical reversed poling. It shows dark lines at the boundary of the patterned area resulting from destructive interference, which indicates the reversal of in-plane crystal orientation and corresponding nonlinear optical polarization after reversed electrical poling. (Courtesy of Ref. [25])

4.
Ferroelectric device application of atomically thin α-In₂Se₃

The vacuum potential difference between the top and bottom surfaces of a given α-In₂Se₃ QL can be as high as 1.4 eV^[20]. The sign of this potential difference can be controlled by the direction of polarization. Therefore, when the 2D ferroelectric α-In₂Se₃ is artificially combined with other vdW materials to form multi-layer vdW heterojunction, the physical properties of the entire heterojunction will be manipulated by the electric polarization. Tuning the properties of the ferroelectric heterojunction as a whole can be simply achieved by reversing the polarity of the ferroelectrics with an external electric field. With this feature, it is convenient to develop multifunctional ferroelectric device with 2D ferroelectrics.



We firstly introduce the ferroelectric diode application based on 2D α-In₂Se₃. With the out-of-plane ferroelectricity of α-In₂Se₃, a polarity switchable ferroelectric diode with the vertical few-layer graphene/α-In₂Se₃/graphene heterojunction is proposed and experimentally demonstrated^[24]. Figs. 4(a) and 4(b) show the device structure and optical image of the proposed ferroelectric diode, in which the thickness of the ferroelectric layer α-In₂Se₃ is as thin as 5 nm. Graphene and few-layer graphene serve as the bottom and top electrodes, respectively. By comparing work functions of α-In₂Se₃, monolayer graphene (p-type) and few-layer graphene (intrinsic), the Schottky barrier at bottom graphene/α-In₂Se₃ is found to be higher than α-In₂Se₃/top graphene. Therefore, the initial state of the diode is in the forward direction, thus the current flows from α-In₂Se₃ to the FLG. When the device is polarized with +2 V pulse voltage, as shown in Fig. 4(c), it shows an enhanced forward rectifying behavior. However, in the case of –2 V polarized state, the device functions as a backward diode, allowing the current conduction from α-In₂Se₃ to bottom graphene (Fig. 4(d)). This polarity-controlled rectification effect could be understood through the band structure of the device (Figs. 4(f) and 4(g)). In the α-In₂Se₃ ferroelectric diode, the polarization reversal in α-In₂Se₃ changes the direction of the built-in electric field. The Schottky barrier height at the interfaces of the upper and lower sides of α-In₂Se₃ changes significantly, leading to the switch of the conduction direction of the diode. Further studies on the I–V characteristic of the ferroelectric diode under high DC bias are also performed. The I–V curves are measured by sweeping the bias voltage from –2 to +2 V and then backward to –2 V. As shown in Fig. 4(e), under high DC bias with both positive and negative voltages, the device is at its “on” state due to the modification of Schottky barrier. The on/off ratio of the ferroelectric diode was found to be at the order of ~ 10⁵.






onerror="this.onerror=null;this.src='http://www.jos.ac.cn/fileBDTXB/journal/article/jos/2019/6/PIC/19030026-4.jpg'"
class="figure_img" id="Figure4"/>



Download



Larger image


PowerPoint slide






Figure4.
(Color online) Switchable ferroelectric diode based on α-In₂Se₃ thin layers. (a) Schematic and (b) optical image of the device. (c) and (d) I–V curves of the ferroelectric diode with switchable rectifying behavior. (e) I–V curves measured under high DC bias, showing clear hysteresis characteristics. The arrows indicate the voltage sweeping sequence. (f) and (g) Schematic of energy band diagrams of the graphene/In₂Se₃ heterostructure, illustrating the evolution of the Schottky barrier in the polarized state of the ferroelectric. The positive and negative charges on the vertical grey lines stand for the polarization charges on the top and bottom sides of the α-In₂Se₃ thin layer. The screening charges are visualized in the graphene/FLG electrodes. (Courtesy of Ref. [24])




A second typical example of the ferroelectric device application is the nonvolatile memory via the device so called FeFET, which is a kind of the three-terminal device by separating the reading operation in the conducting channel from the writing operation in and the ferroelectric gate. FeFETs are non-destructive with no requirement on ferroelectric film thickness. However, the thickness of ferroelectric used in conventional FeFET devices are usually of several hundred nanometers due to the constraint of Critical size effect. Therefore, relatively large operation gate voltage is always required. In contrast, 2D ferroelectrics will significantly scale down the dimensions of the device and meanwhile the voltage to flip the electric polarization could be lowered down dramatically, resulting a low-power consumption memory device.



Fig. 5 shows the schematic and optical image of a FeFET device based on graphene and ferroelectric α-In₂Se₃^[29]. Atomically thin hexagonal boron nitride (hBN), monolayer or bilayer, was chosen as the buffer and insulating layer. The introduction of hBN layer not only improvs the device interface properties, such as suppression of ion diffusion and leakage current^{[40, 41]}, but also enhances the intrinsic conductivity of graphene by reducing the phonon and Coulomb scatterings. The thinnest thickness of the ferroelectric film in the device is 2.6 nm, which is roughly about the thickness of bilayer α-In₂Se₃. This thickness is much thinner than conventional ferroelectric layer. The resistance of graphene channel is controllable and retentive due to the electrostatic doping, which stems from the electric polarization of the ferroelectric α-In₂Se₃. As shown in Fig. 5(b), when sweeping the ferroelectric gate voltage cyclically from negative to positive and back to negative, the resistance of the graphene channel could be efficiently modulated with large and stable hysteresis. By applying the top-gate voltage, the electric dipoles in α-In₂Se₃ are electrically polarized in the out-of-plane direction, either upward or downward. The dipoles would significantly dope graphene with different types of carriers via the induced screening charges, leading to the change of its Fermi level (E_F) accompanied by the resistance modulation of the FeFET device. When the Fermi level is tuned to the Dirac point, the conducting graphene channel reaches its maximum resistance (R_Max) state. Two R_Max states are found at different top gate voltages in the forward (from negative to positive) and backward (from positive to negative) gate voltage sweeping. The apparent differences of gate voltage at R_Max are due to the requirement of coercive electric field for domain reversion in ferroelectric α-In₂Se₃. However, even under ±5 V gate bias, equivalent to an electrical field strength at the order of 10⁵ V/cm in this device, the ferroelectric polarization in the top α-In₂Se₃ thin layers does not get saturated. Moreover, it is necessary to point out that the observed hysteretic dependence of resistance on V_TG allowed us to explore the reversal processes of the electric dipoles in ferroelectric α-In₂Se₃.






onerror="this.onerror=null;this.src='http://www.jos.ac.cn/fileBDTXB/journal/article/jos/2019/6/PIC/19030026-5.jpg'"
class="figure_img" id="Figure5"/>



Download



Larger image


PowerPoint slide






Figure5.
(Color online) Structure and optical characterization of the 2D FeFET. (a) 3D schematic diagram of the FeFET. The FeFET is fabricated by vertically stacking graphene, hBN, and α-In₂Se₃ thin layers in sequence. The white arrows indicate the direction of electric polarization. The zoomed area shows the crystal structure of ferroelectric α-In₂Se₃. (b) The hysteretic ferroelectric loop in 2D α-In₂Se₃ based FeFET. The resistance follows a butterfly-like dependence on gate voltage. Inset shows the optical image and topography of the FeFET. The substrate is Si wafer with 300 nm fused SiO₂ on top. The graphene, hBN, and ultrathin α-In₂Se₃ are indicated by black, white, and red-dashed frames, respectively. From the AFM topography, the thinnest area of ultrathin α-In₂Se₃ in the device is 2.6 nm. (Courtesy of Ref. [29])

5.
Conclusions and outlooks

The vdW materials are promising for realizing 2D ferroelectricity which is the long-sought goal in conventional ferroelectrics. Recent progresses, especially the studies on ultra-thin ferroelectric α-In₂Se₃, have confirmed their practicability in device applications. They offer a versatile platform to study the ferroelectricity, domain physics, and ferroelectric device at nanometer scale. For α-In₂Se₃, due to its unique lattice structure, the in-plane and out-of-plane polarizations are inherently coupled, which provides a novel opportunity to control the out-of-plane ferroelectricity by in-plane electric field. However, one of the major bottlenecks for developing 2D ferroelectrics lies in the imperfect sample quality and uniformity. Large-scale ferroelectric thin-film synthesis techniques, such as molecular beam epitaxy (MBE) and chemical vapor deposition (CVD), are highly desired for future practical device application.

Acknowledgments

This work was supported by the National Key Research and Development Program of China (Grant Nos. 2017YFA0205004, 2018YFA03066004, and 2016YFA0301700), the National Natural Science Foundation of China (Grant Nos. 11674295 and 11774328), the Fundamental Research Funds for the Central Universities (Grant No. WK2340000082), Anhui Initiative in Quantum Information Technologies (Grant No. AHY170000), the USTC start-up funding and the China Government Youth 1000-Plan Talent Program.