1 Introduction
Ferroelectric materials possess stable and switchable spontaneous polarization attributes, rendering them suitable for deployment in diverse electronic devices, including sensors, memories, and piezoelectric components [
1]. Nevertheless, the advancement of microelectronics technology necessitates the fabrication of low-power devices with higher density. The conventional perovskite ferroelectrics [
2‒
4] encounter limitations due to depolarization fields [
5,
6] and incongruity with the substrate lattice [
7]. Hence, the advancement of two-dimensional ferroelectric layered materials presents a viable approach to address this unparalleled challenge. Primarily, the atomic-level thickness facilitates the utilization of two-dimensional ferroelectric materials in the quest for enhanced integration and reduced power consumption. Additionally, a pristine van der Waals interface fosters the amalgamation of crystals and diverse substrates, emerging as a significant advantage for integration in semiconductor circuits [
8]. Furthermore, the presence of a robust intralayer force and a relatively weaker interlayer force interaction suggests the potential for reduced interference from depolarization effects and the preservation of ferroelectricity at the atomic layer level. Additionally, the existence of distinctive dangling bonds, surface states, and the ability to withstand significant strain in two-dimensional (2D) materials offer significant advantages for the integration of two-dimensional ferroelectric materials with semiconductor technologies [
9].
The theoretical anticipation of two-dimensional ferroelectric materials can be traced back to 2013 [
10], effectively predicting potential candidates with inherent ferroelectricity, such as 1T-MoS
2 [
11,
12], In
2Se
3 [
13,
14], CuInP
2S
6 [
15,
16], and MX (M = Ge, Sn, and X = S, Se, Te) [
17‒
19]. The prediction of two-dimensional ferroelectricity is advantageous for the exploration of two-dimensional ferroelectric materials and their potential application in nano electronic devices. We have successfully synthesized a two-dimensional ferroelectric material, CuInP
2S
6, and observed switchable polarization in a thin film with a thickness of 4 nm at room temperature [
20,
21]. However, the attribution of this phenomenon can be ascribed to the active Cu element, which poses challenges in achieving stable synthesis [
15]. Furthermore,
α-In
2Se
3 has been successfully fabricated as a ferroelectric controllable heterojunction, leading to a selective improvement in optical response [
22,
23]. Despite the experimental discovery of certain 2D ferroelectric materials mentioned earlier, the achievement of stable ferroelectricity in 2D materials at room temperature remains a formidable obstacle.
As mentioned above, the experimental realization of two-dimensional room-temperature ferroelectric materials remains limited. Here, we report a novel 2D room-temperature ferroelectric material: SnP2S6 atomic layers, which can be prepared by the chemical vaper transport (CVT) growth of bulk SnP2S6 crystals and subsequent exfoliation bulk crystals into nanosheets. RT stable ferroelectric order can be achieved in 7 nm SnP2S6 atomic layers. The ferroelectric field effect transistor with SnP2S6 as the top gate insulator and WTe0.6Se1.4 as the channel material exhibits an obvious clockwise hysteresis loop, further identifying the ferroelectricity of SnP2S6. Furthermore, multilayer Graphene/SnP2S6/multilayer Graphene van der Waals vertical heterostructure phototransistor was successfully fabricated, possessing improved optoelectronic performances with responsivity (R) of 2.9 A/W and detectivity (D) of 1.4 × 1012 Jones. The experimental results have proven that SnP2S6 is an excellent 2D ferroelectric material, which is expected to provide new choices for low-power and highly integrated precision devices.
2 Experimental section
Sample preparation. Bulk SnP2S6 crystals were synthesized by the chemical vapor transport technique. At first, high pure Sn, P, and S powders (≥99.9%) with a molar ratio of 1:2:6 were mixed together (totally 0.5 g), and iodine was used as a transport agent. Then, they were sealed into a vacuum quartz tube (about 10−3 Pa) through high-temperature torch. Subsequently, the quartz tube was placed horizontally in a dual-temperature zone tubular furnace. The temperature of the source zone and growth zone were heated to 750 °C and 700 °C within 24 hours and kept at this temperature for 7 days so that the reaction could proceed completely. After heating, the quartz tube was naturally cooled to room temperature with the furnace. Finally, the plate-like SnP2S6 bulk crystals were successfully obtained after completing this growth process and 2D SnP2S6 flakes with various thicknesses were cleaved from the bulk crystals by mechanical exfoliation method.
XRD characterization. The powder X-ray diffraction patterns were collected using an X-ray diffractometer (D/MAX-2200, Rigaku) with monochromatized Cu-Karadiation at room temperature in the 2θ range of 10°‒85° with a scan speed of 10°/min.
EDS characterization. Semi-quantitative microprobe analyses on the single crystals were performed with the aid of a field emission scanning electron microscope (Quanta 400 F, FEI) equipped with an energy dispersive X-ray spectroscope (EDS, Oxford INCA), which was applied to obtain semi-quantitative analysis results, like atomic ratio and elements dispersion mapping.
Optical characterization. SHG measurements were performed in reflection geometry with 100 fs pulses at 786 nm and a repetition rate of 76 MHz, which were focused to a spot size of ~1 μm by a 40 × 0.6 NA objective lens (Olympus). The PL spectra were recorded by a confocal microscope spectrometer (Alpha300 Raman, WITec). The optical absorption spectra were collected by a microscopic spectrophotometer (MSV-5200, Jasco).
Device fabrication. The devices were fabricated by laser direct imaging (μpg501, Heisenberg), and the corresponding Ti/Au (10 nm/100 nm) electrodes were deposited by electron beam evaporation (DE400, wavetest). (Opto)electronic measurements were performed under ambient conditions using a Keithley 4200 semiconductor parameter analyzer in a probe station. 405, 532, and 635 nm monochromatic lasers with tunable power were applied to investigate the photoresponse of SnP2S6.
AFM and PFM. Atomic force microscopy (AFM, Asylum Research Cypher S) in a tapping mode was used to characterize the morphology of the heterostructures and devices. The thickness of the SnP2S6 was identified by AFM. Out-of-plane PFM measurements were carried out on the Cypher S AFM in the DART mode. Off-field PFM hysteresis loops were measured by recording the piezoresponse amplitude and phase signals after the individual DC pulse was turned off. P−E hysteresis loop measurements were carried out using a Radiant ferroelectric tester with an applied voltage of a triangular waveform at 1 kHz. Static transport properties were measured with an Agilent B1500A Semiconductor Device Parameter Analyzer in a vacuum chamber of 10−2 torr. The dynamic and writing speed tests were performed with the top-gate electrode connected to a RIGOL DG1032Z signal generator and source and drain connected to the Agilent B1500A.
TEM sample preparation and image simulations. The conventional wet-transfer method with the help of PMMA was used to transfer the exfoliated SnP2S6 flakes onto a Cu grid with a carbon film for preparation as a TEM sample. Z-contrast STEM imaging was carried out with a modified JEOL 2100F equipped with a delta probe corrector, which was used to correct for aberration up to the fifth order, resulting in a probe size of 1.4 Å. The imaging was conducted at an acceleration voltage of 60 kV. The convergent angle for illumination was ~35 mrad, with a collection detector angle ranging from 62 to 200 mrad.
3 Results and discussion
SnP
2S
6 crystallizes in a 2D atomic structure with
C3 point group and space group
R3 (No. 146) [Fig.1(a, b)] in which
a =
b = 5.999 Å,
c = 19.424 Å,
α =
β = 90°,
γ = 120°,
V = 605.38 Å
3, which was first reported in 1995 [
24]. SnP
2S
6 exhibits a unique pore free structure as it lacks half of the metal ions compared to the parent Sn
2P
2S
6 structure [
25,
26]. Due to its inherent 2D nanoporous structure, recent experimental studies have shown that SnP
2S
6 exhibits strong nonlinear optical response related to its inversion symmetry broken structure [
25,
27]. It is worth noting that SnP
2S
6 monolayer consists of close-packed sulfur atomic skeleton, in which one third of vacancies with octahedral form are occupied by Sn (IV) ions in the form of [SnS
6]
8–, one third of vacancies are occupied by P–P pairs in the form of [P
2S
6]
4– and the remaining ones are unoccupied. SnP
2S
6 bulk crystals was synthesized by the chemical vaper transport (CVT) method using iodine as a transport agent [Fig.1(c)], see the Method Section for more details), which have a larger size than 5 mm and exhibit single-crystal morphologies, implying successful growth of high-quality single crystals [Fig.1(d)]. The X-ray diffraction spectrum (XRD) of as-prepared SnP
2S
6 bulk crystals is in good agreement with standard PDF card (79818-icsd), suggesting high-purity phase of the as-prepared samples. Additionally, it is also noted that the XRD pattern of single-crystal SnP
2S
6 give strong peaks of (
), (
), (
), and (
), arising from the strong grain orientation along
c axis [Fig. S1 of the Electronic Supplementary Materials (ESM)] [
25]. To mitigate the impact of grain orientation on X-ray diffraction (XRD) analysis, the as-prepared SnP
2S
6 bulk crystals were pulverized into a fine powder. The resulting XRD pattern of the powder exhibits a satisfactory correspondence with the standard PDF card (79818-icsd) [Fig.1(e)], wherein the intensity of four prominent grain-orientation peaks is diminished, while the emergence of peaks from alternative crystal planes is observed.
The layered structural behavior of SnP
2S
6 endows as-grown bulk crystals with sheet-like morphology with the
c axis perpendicular to the surface plane of the sheets [Fig.2(a)]. The energy dispersive X-ray spectrum (EDS) results show that the atomic ratio of Sn:P:S in as-grown bulk crystals is 0.97:2.03:5.99 [Fig.2(a)], which is basically consistent with the stoichiometric ratio of 1:2:6, implying high qualities of as-grown bulk crystals. The atomic structure was further characterized and proven to be consistent with the atoms in the
a‒
b plane by scanning transmission electron microscopy (STEM) image [Fig.1(b) and Fig. S2 of the ESM]. SnP
2S
6 atomic layers with thickness gradient were prepared on SiO
2/Si (100) substrate via mechanical exfoliation method [
28], in which the size of the monolayer with uniform thickness can reach 13 μm [Fig.2(c)]. Micro-regions UV–Vis–NIR absorption spectrum is used to measure the bandgaps of SnP
2S
6 atomic layers, in which the absorption of SnP
2S
6 decreases with decreasing thicknesses [Fig.2(d)] as well as the bandgaps can be changed from 1.96 eV (bulk sample) to 2.58 eV (7.6 nm nanosheets) [Fig.2(e)]. In order to evaluate the interlayer interaction of SnP
2S
6, photoluminescence spectra (PL) of SnP
2S
6 atomic layers at different temperatures were measured [Fig.2(f)‒(h)], in which the PL peaks of the trilayer, bilayer and monolayer are located between 715 nm and 735 nm, basically consistent with the previous reports [
29]. Additionally, under the conditions of 78 K and 150 K, PL peaks can be observed at 725 nm and 735 nm, respectively, indicating that the band gap of thin-layer SnP
2S
6 decreases with increasing temperature. It is also worth noting that as the thickness decreases, the intensity of the PL peak also decreases, which is consistent with the absorption peak mentioned earlier. Second nonlinear optical characterizations of SnP
2S
6 were performed to identify the non-centrosymmetric structure of SnP
2S
6. Under the excitation of 800 nm laser, a second harmonic generation (SHG) signal with 400 nm was stably generated, further confirming the breaking inversion symmetry of as-prepared SnP
2S
6 [Fig.2(i)]. Interestingly, the resulting polarization of SHG intensities displays the six-fold rotational symmetry [Fig.2(j, k)], identifying that as-grown samples belong to the
C3 point group with broken inversion symmetry [
30,
31]. In addition, the Raman spectra of bulk SnP
2S
6 were also measured, in which there were three obvious Raman peaks identified as
P1 (140 cm
−1),
P2 (166 cm
−1), and
P3 (266 cm
−1), consistent with previous reports [
25]. It is found that, for few-layer SnP
2S
6, only weakened P
3 peak can be observed and no P
1 and P
2 peaks can be found (Fig. S3 of the ESM), which is also observed in previous reports [
25]. Such weak Raman intensities can be attributed to weak optical absorption of few-layer SnP
2S
6.
In order to explore the ferroelectricity of SnP
2S
6 thin films, the atomic force microscope (AFM) and piezoresponse force microscope (PFM) tests were carried out on the exfoliated SnP
2S
6 few layers (see the Method Section for details). A 15-nm-thick SnP
2S
6 nanosheet [Fig.3(a)] was fabricated via mechanical exfoliation method to study its ferroelectricity, which possesses good uniformity in thickness [Fig.3(d)], making it an ideal sample to eliminate the interference of surface fluctuations on piezoelectric response force at the nanoscale. There are dendritic and partially black dot patterns in both PFM amplitude [Fig.3(b)] and PFM phase [Fig.3(c)], where PFM amplitude reflects the modulus of local piezoelectric response and phase reflects the polarization direction of each region. The values of amplitude [Fig.3(e)] and phase [Fig.3(f)] are extracted from the same position of the film, indicating their good consistency. The PFM signals obtained from the clean surface indicates the existence of spontaneous polarized ferroelectric domains on the film, which is the direct evidence of the ferroelectricity of SnP
2S
6 at room temperature. Examples of more ferroelectric domains in SnP
2S
6 with different thicknesses are given in Fig. S4 of the ESM. Although the ferroelectricity of layered SnP
2S
6 has been proven, the switchable polarization of ferroelectrics also needs to be proven because some effects, such as depolarization effect [
32], interface pinning effect [
33], and “dead layer” effect [
34] on thin films, will lead to the absence of the ferroelectricity. The desired pattern was obtained in the PFM test by applying voltage between the conductive tip and the substrate. For a flat sheet [Fig.3(g)] with a thickness of 7 nm (Fig. S5 of the ESM), a DC voltage of +2 V, −2 V and +2 V is applied from the outside to the inside. The clearly layered rectangular amplitude pattern shows that SnP
2S
6 nanosheet has controllable spontaneous polarization property [Fig.3(h) and Fig. S6 of the ESM]. In addition, the clear ferroelectric hysteresis loops were also observed on the phase diagram and amplitude diagram [Fig.3(i)], further demonstrating the switchable spontaneous polarization of SnP
2S
6 as a ferroelectric material.
Despite SnP
2S
6 as an ultra-thin ferroelectric material has been demonstrated by PFM tests, only one way to identify the RM ferroelectricity is inadequate. Therefore, a ferroelectric field-effect transistor (Fe-FET), with SnP
2S
6 atomic layers as the top ferrielectric gate, thin h-BN layer as the top gate insulator, and p-type semiconductor WTe
0.6Se
1.4 as a channel material at the bottom [Fig.4(a)], was fabricated to further confirm the RM ferroelectricity of SnP
2S
6, in which the introduction of h-BN effectively improves the interface contact performance and prevents gate current breakdown [
35,
36]. In addition, graphene is used to connect the channel with the metal electrodes, which, indeed, also act as the electrodes. Fig.4(b) shows an optical image of as-fabricated ferroelectric transistor, from which it can be clearly seen that the gate completely covers the channel material.
The output curve of the Fe-FET was measured at room temperature with a top-gate
Vg varying from ‒0.6 V to 0.6 V at a step of 0.2 V [Fig.4(c)], which are almost straight, indicating that the device has good ohmic contact. When the top-gate voltages as a floating gate are changed from 0.5 V to ‒0.5 V, then ‒0.5 V to 0.5 V, the Fe-FET exhibit a clear clockwise hysteresis loop in transfer characteristics, also demonstrating ferroelectric properties of SnP
2S
6. Meanwhile, it is worth noting that a small change of
Vg leads to the change of drain current, revealing that, in SnP
2S
6-based Fe-FET, the ferroelectric state switching is easy and can work under low power consumption. Moreover, the transfer curve of Fe-FET shows typical p-type semiconductor characteristics of few-layer WTe
0.6Se
1.4, in which the ferroelectric resistance switching can realizes more than two orders of magnitude changes [Fig.4(d)]. The switching of the state comes from the switching of ferroelectric polarization, which changes the carrier of the channel material [
37]. When the Fe-FET is applied in a negative top-gate voltage, the valence band edge (
EV) will bend above the Fermi level (
EF) of the p-type semiconductor, leading to the generation of holes and the formation of low channel resistance states (LRS) in the semiconductor layer. On the contrary, the positive
Vg will lead to the generation of moving electrons due to the occurrence of polarization transformation, further cause the generation of high channel resistance states (HRS). Therefore, owning to the existence of residual polarization, the channel resistances will have a clockwise hysteresis loop under a floating gate. The characterizations of Fe-FET identify that SnP
2S
6 is a novel 2D ferroelectric material again, which can be used in nonvolatile memory device.
To investigate the optoelectronic properties of SnP
2S
6 thin film, a parallel photodetector based on 75 nm SnP
2S
6 nanosheet was prepared [Fig.5(a)], which was excited by a 405 nm laser (with a spot diameter of 1 cm, much larger than the size of the device). The output curve of the parallel phototransistor transistor [Fig.5(b)] exhibits a quasi-Ohmic contact, and the photocurrent is related to the power density of the laser. Under 405 nm laser irradiation, the photocurrent tends to saturate with increasing light power density. In order to further evaluate the performances of the device, responsivity (
R) and detectivity (
D) are introduced. The responsivity is defined by
R =
Iph/(
PS), where
Iph is photocurrent (
Iph =
Ilight ‒
Idark),
P is light power density,
S is active channel area,
Ilight and
Idark are the drain currents under light and dark. Detectivity is a powerful parameter for comparing the capability to detect weak light signal of photodetectors with various geometries. It is defined as
D =
RS1/2/(2
qIdark)
1/2, where
R is the responsivity,
S is the effective area of the device, and
Idark is the dark current. The parallel device shows a photoresponsivity (
R) of 3.6 × 10
−2 A/W and a detectivity (
D) of 4.2 × 10
10 Jones at a light power density of 41.3 μW·cm
−2 and a bias voltage of 2 V. Additionally, it was also observed that the
R and
D decrease with increasing power densities of the lasers [Fig.5(c)], attributed to severe frequent carrier recombination and shorter carrier lifetime at stronger illumination. The value of
α, extracted from
Iph ∝
Pα, is introduced to judge the mechanism of photocurrent generation, where
Iph is drain current under light and
P is light power density [
38]. The factor
α of 0.59 is achieved in our parallel photodetector, which is smaller than 1, indicating that the recombination kinetics of photocarriers involves both trap states and interactions between photogenerated carriers [
39], and it may have a certain impact on the reaction speed of the device. Generally, the rise time and the decay time are used to evaluate the response speed of devices, fitted by the following equation:
Irise =
I0 –
A exp(–(
t –
t1)/
τ1) and
Idecay =
I0 –
B exp(–(
t –
t2)/
τ2), where
τ is the time constant,
t is the time when laser is switched on or off, and
A and
B are the scaling constants. Therefore, the fitted characteristic photoresponse time coefficients,
τ1 and
τ2, are 100 ms and 200 ms, respectively. It is worth noting that the rise time is faster than the decay time because the photocarriers can be trapped in the trap states to further result in an increase in the lifetime [
40‒
43]. In addition, the optoelectronic performances of our parallel device is also assessed under illumination of 532-nm and 635-nm lasers, in which the photoresponsivities and the rise times are 2.7 × 10
−3 A/W and 260 ms for the former [Fig. S7(b) of the ESM], 1.4 × 10
−3 A/W and 2 s for the latter [Fig. S7(c) of the ESM], respectively. The poor optoelectronic performances under 532-nm and 635-nm lasers should be attributed to the poor optical absorption.
Based on our previous reports that fully vertical van der Waals heterostructure photodetector can efficiently improve the optoelectronic performances, fully vertical van der Waals heterostructure photodetector of CuBiP
2Se
6 [
44], a multilayer graphene/SnP
2S
6/multilayer graphene van der Waals vertical heterostructure phototransistor, was fabricated with 15 nm channel width [Fig.6(a)]. Photoresponse is observed at a low light power density of 1.9 μW·cm
−2 and a stronger photocurrent was exhibited [Fig.6(b)]. It is worth noting that the great improved
R of 2.9 A/W and
D of 1.4 × 10
12 Jones are achieved because the short channel effect of our vertical device efficiently reduces the scattering of photocarriers in the channel [
45,
46]. The
α of 0.41 indicates that vertical devices also do not change the influence of their intrinsic trap states on photogenerated charge carriers [Fig.6(d)]. The rise time of vertical device is significantly improved to 20 ms due to the introduction of heterojunctions and nanoscale channel length. However, the improved decay time of 40 ms is still longer than the rise time because of its intrinsic recombination kinetics. Furthermore, the optoelectronic performances of our vertical device are also performed under illumination of 532-nm and 635-nm lasers, in which the photoresponsivities and the rise times are 0.02 A/W and 200 ms for the former [Fig. S8(b) of the ESM], 0.0015 A/W and 300 ms for the latter [Fig. S7(c) of the ESM], respectively. As in the parallel device, the poor optoelectronic performances of the vertical device under 532-nm and 635-nm lasers should be also attributed to the poor optical absorption.
4 Conclusion
In this work, we developed a novel two-dimensional ferroelectric material, van der Waals SnP2S6, with room-temperature ferroelectricity in ~7 nm thick flake. A Fe-FET with ferroelectric SnP2S6 was designed and fabricated successfully, which exhibits a clear clockwise hysteresis loop in transfer characteristics, and the ferroelectric behavior has the potential to be applied to non-volatile storage device. Meanwhile, we also designed and fabricated a multilayer graphene/SnP2S6/multilayer graphene van der Waals vertical heterostructure phototransistor with improved R of 2.9 A/W and D of 1.4 × 1012 Jones. 2D ferroelectric SnP2S6, with its ability to integrate with mature silicon-based platforms, enriches the family of ferroelectric materials and the functionality of two-dimensional materials, making it an excellent candidate for the next generation of nano-electronic devices.