Quadruple-well ferroelectricity and moderate switching barrier in defective wurtzite α-Al2S3: a first-principles study

doi:10.21203/rs.3.rs-4680101/v1

Download PDF

Article

Quadruple-well ferroelectricity and moderate switching barrier in defective wurtzite α-Al2S3: a first-principles study

https://doi.org/10.21203/rs.3.rs-4680101/v1

This work is licensed under a CC BY 4.0 License

You are reading this latest preprint version

Wurtzite-type ferroelectrics are highly promising for next-generation microelectronic devices due to their ferroelectric properties and integration with exiting semiconductors. However, their high coercive fields, which are close to breakdown electric fields, need to be lowered. To deal with this issue and secure device reliability, much effort has been devoted to exploring novel wurtzite compounds with lower polarization switching barriers and implementing doping strategies. Here, we report first-principles calculations on polarization switching in cation-vacancy ordered wurtzite α-Al₂S₃, unveiling its uniaxial quadruple-well ferroelectricity and moderate switching barrier, 51 meV/cation, which is much lower than that of conventional wurtzite ferroelectrics. There are three important features relevant to the Al vacancies leading to the uncommon quadruple-well ferroelectricity and the moderate switching barrier: mitigation of cation-cation repulsion, structural flexibility that alleviates an in-plane lattice expansion, and formation of s-like bonding states consisting of Al 3p_z and S 3p_z orbitals. Biaxial compressive strain and Ga doping lower the switching barriers by up to 40%. This study encourages experimental investigation of the ferroelectric properties for defective wurtzite α-Al₂S₃ as a new promising material with unconventional and intriguing ferroelectricity and suggests a potential strategy for reducing switching barriers in wurtzite ferroelectrics: introducing cation vacancies.

Physical sciences/Materials science/Materials for devices/Electronic devices

Physical sciences/Chemistry/Materials chemistry/Electronic materials

Ferroelectric materials, characterized by electrically and/or mechanically switchable spontaneous polarization, are utilized in various devices such as piezoelectric actuators, pyroelectric sensors, and capacitors ^1–3. Recently, wurtzite-type compounds have been attracting much attention. This class of materials was known to be piezoelectric and pyroelectric owing to their polar crystal structures with P6₃mc space group symmetry. The possible ferroelectricity of wurtzite-type compounds is highly promising for cutting-edge computing and data storage devices due to their high spontaneous polarization. However, their ferroelectricity was not observed prior to dielectric breakdown because of their high coercive fields until Fichtner et al. clearly demonstrated ferroelectric switching for wurtzite-type Sc-doped AlN in 2019 ⁴. Their work hinted at doping strategies that facilitate ferroelectric switching and ignited rapid progress in the exploration of wurtzite ferroelectrics. Consequently, ferroelectric switching has been reported for other doped wurtzite compounds such as Al_1-xB_xN, Ga_1-xSc_xN, Zn_xMg_1-xO, and Al_1-xY_xN ^5–8. However, their high coercive fields, close to breakdown electric fields, reduce device reliability ^9,10. To address this issue, considerable effort has been made to explore novel wurtzite compounds with lower switching barriers and substantial breakdown electric fields^10–12 as well as to implement the doping strategy as mentioned above ^4–8.

Wurtzite structures consist of an hcp anion arrangement with cations occupying half of the 4-fold coordinated tetrahedral sites (Fig. 1b and d). The filled tetrahedra are all oriented upwards (or all oriented downwards), sharing a vertex. It is commonly considered that upon ferroelectric switching in typical wurtzite-type compounds such as pristine AlN, all the cations move collectively from the upwards-oriented tetrahedra to neighboring downwards-oriented tetrahedra along the c-axis direction by passing through anion basal triangles, which involves a bottleneck of switching ^13–15. The saddle point of the minimum energy pathway (MEP) is a nonpolar hexagonal-boron nitride (h-BN)-like structure with space group symmetry of P6₃/mmc, in which all the cations locate at 5-fold coordinated trigonal bipyramidal sites ¹⁴. Doping isovalent cations and applying biaxial tensile strain reduce the switching barriers by destabilizing the polar wurtzite structures and/or stabilizing the saddle-point h-BN-like structures to achieve switching prior to breakdown ^16,17. Recently, Lee et al. reported a first-principles study suggesting that certain ternary wurtzite-type compounds containing two types of cations exhibit polarization switching not in a "collective" way via a h-BN-like structure, but in an "individual" or "stepwise" way via multiple local-minimum intermediate structures ¹⁸, highlighting that avoiding the saddle-point h-BN-like structures enables those multinary wurtzite compounds to have low switching barriers.

Here, we propose introducing cation vacancies as a potential strategy to lower the switching barriers through a first-principles study on the polarization switching in a defective wurtzite compound. γ-In₂Se₃ is known as a cation-vacancy ordered wurtzite (Fig. 1a and c) ^19,20. About 30 compounds with γ-In₂Se₃-type structures are registered in the Inorganic Crystal Structure Database (ICSD). In this study, we focus on α-Al₂S₃, which consists of the most abundant elements found in the γ-In₂Se₃-type compounds. Our first-principles calculations in conjunction with solid-state nudged elastic band (SS-NEB) methods²¹ unveiled MEPs for polarization switching in α-Al₂S₃, which manifest itself as uniaxial quadruple-well potential curves with four local-minimum states. The intermediate local-minimum structures encompass 5-fold trigonal bipyramidal coordination. The calculated switching barrier is much lower than that of conventional wurtzite ferroelectrics. Detailed analysis of the evolution of the atomic arrangements and chemical bonding during the polarization switching indicates that the Al vacancies play important roles in yielding the unconventional quadruple-well ferroelectricity and the moderate switching barrier. Our calculations also predict that biaxial compressive strain and Ga doping enable the tunability of switching barriers. The calculated piezoelectric constants suggest that piezoresponse force microscopy (PFM) allows us to distinguish the four local-minimum states with different polarization values experimentally. This work showcases introducing cation vacancies as a potential route to lower the switching barriers in wurtzite ferroelectrics as well as predicts the unconventional ferroelectricity in α-Al₂S₃.

2.1 Crystal Structure of α-Al₂S₃ and Polarization Switching Behavior

The crystal structures of cation-vacancy ordered wurtzite α-Al₂S₃ and conventional wurtzite ZnS are summarized in Fig. 1. α-Al₂S₃ shows a polar and chiral P6₁ space group symmetry. Two-thirds of the tetrahedral sites in wurtzite-type structures are occupied by the Al atoms, while the remaining one-third are vacant ^22,23. In contrast to conventional wurtzite-type compounds, α-Al₂S₃ has two kinds of crystallographically inequivalent cation sites, both of which are located at Wyckoff positions 6a. These sites can be readily distinguished based on the position of the nearest-neighbor Al vacant sites, i.e., whether it locates either in the sulfur vertex or basal-triangle side along the c axis. Figure 1e and f extract a column of the coordination tetrahedra along the c axis from the whole structure of α-Al₂S₃ and ZnS, respectively. Hereafter, the tetrahedral sites with a nearest-neighbor vacant site in the vertex or basal-triangle side are referred to as T1 and T2, respectively. The suffixes "u" and "d" indicate the upwards-oriented or downwards-oriented tetrahedra, respectively.

The calculated lattice constants (a = 6.41 Å, c = 17.80 Å) for α-Al₂S₃ are in good agreement with the experimental ones (a = 6.44 Å, c = 17.90 Å) ²³, substantiating that the GGA-PBEsol functional used in this study well reproduces crystal structures. The calculated electric polarization of α-Al₂S₃, 66 µC/cm², is comparable to or smaller than those of wurtzite ferroelectrics (e.g., 65 and 135 µC/cm² for ZnS and AlN, respectively) ^13,15. An MEP between the two oppositely polarized states with all the tetrahedra oriented upwards or downwards was calculated by using SS-NEB methods to elucidate the switching barrier and behavior (Fig. 2a). Interestingly, α-Al₂S₃ has four local energy minima in the MEP, referred to as + HP, +LP, –LP, and –HP in the order of the magnitude of polarization. Here, the states with high, low, and zero absolute polarization values are denoted as HP, LP, and ZP, respectively, and the prefix "+" and "–" symbols indicate positive and negative polarization values, respectively. Figure 3 illustrates these five structures providing the local energy minima and the saddle point. (An animation of the structural evolution associated with the switching is available in the Supporting Information and would help to understand the switching behavior.) Both + HP and + LP states possess P6₁ space group symmetry. The P6₁ symmetry is preserved during the switching pathway, except in the saddle-point ZP state, which has zero polarization and nonpolar P6₁22 symmetry. The 6₁ symmetry element is preserved during the entire polarization switching process, indicating that the polarization direction remains along the c axis. In contrast to quadruple-well ferroelectrics such as BiFeO₃, which have multiple independent polarization axes ²⁴, the four local-minimum states in the MEP of α-Al₂S₃ are polarized along the c axis, revealing that α-Al₂S₃ is an unusual ferroelectric with a uniaxial quadruple-well potential as observed in CuInP₂S₆ ²⁵.

Figure 2b illustrates the structural evolution of α-Al₂S₃ during the polarization switching as the columns of coordination tetrahedra along the c axis for simplicity. All the six columns included in the unit cell are related to each other by the 6₁ symmetry operation, which is preserved during the polarization switching. The atomic displacement in the column reflects the entire polarization switching process. In the initial high-polarization + HP state, the two inequivalent Al atoms, denoted as Al1 and Al2, occupy the T1u and T2u sites, respectively, enclosed by the upwards-oriented tetrahedra. Beyond a small potential hill, there emerges the + LP state with lower total energy and polarization (+ 30 µC/cm²), where the Al1 atoms occupy the T1u sites whereas the Al2 atoms occupy 5-fold coordinated bipyramidal sites, referred to as B1 sites. In the structural evolution from the + HP to + LP states, the Al2 atoms move from the T2u sites towards the vacant sites and migrate to the bipyramidal B1 sites, whereas the Al1 atoms remain at the T1u sites. Upon transitioning from the + LP to ZP states, the Al2 atoms move from the bipyramidal B1 sites to the T1d sites in the downward-oriented tetrahedra, whereas the Al1 atoms still remain at the T1u sites, resulting in a half-switched state where half of the tetrahedra and the other half are oriented upward and downward with the polarization canceled out.

In the latter half of polarization switching from the ZP state to the –LP state to the final –HP state, the Al1 atoms at the T1u sites migrate to the T2d sites via the bipyramidal B2 sites. It should be noted that the Al1 (Al2) atoms sitting at the T1u (T2u) sites in the initial + HP state occupy the T2d (T1d) sites in the final –HP states, indicating the swapping of the crystallographically inequivalent sites during the polarization switching between the T1 and T2 sites. This highlights the nonconventional ferroelectricity of α-Al₂S₃.

Our first-principles phonon calculation revealed that both the HP and LP states were dynamically stable (Figure S1 and S2). Notably, Al₂S₃ with the LP structure has been synthesized by chemical vapor transport ²³, while solid-state reaction methods yield α-Al₂S₃, which adopts the HP structure. These facts imply that these isosymmetric polymorphs are energetically antagonized at room temperature or above, although our first-principles calculations predict that the LP state is more stable than the HP state at 0 K.

2.2 Comparison of Switching Behavior and Barriers in Wurtzite Compounds

In this study, the switching barriers for quadruple-well potential surfaces are defined as the highest energy barrier between a valley and its neighboring peak towards the switching direction ¹⁸. The switching barrier of α-Al₂S₃ corresponds to the total energy difference between the LP and ZP states. The switching barrier of α-Al₂S₃, 51 meV/cation, is approximately one-tenth that of AlN (523 meV/cation) and one-third that of Al_15/16B_1/16N (150 meV/cation) ²⁶. The moderate switching barrier is anticipated to enable polarization switching prior to electric breakdown. We compare the switching behavior in α-Al₂S₃ and other wurtzite ferroelectrics to understand the moderate switching barrier in α-Al₂S₃ below.

In the binary wurtzite ferroelectrics such as pristine AlN, all the cations displace collectively from the upwards-oriented to the downwards-oriented tetrahedra during polarization switching ^12,14. In ternary wurtzite systems such as Li₂SiO₃, for which two-step polarization switching has been predicted by first-principles calculations, more electronegative atoms move first, followed by the migration of less electronegative atoms, with the first switching barrier being the highest ¹⁸. In contrast to these cases, in α-Al₂S₃, consisting of only one type of cations, the Al atoms situated at the crystallographically different sites exhibit individual motion; the Al2 atoms exhibit displacement preceding that of the Al1 atoms, as shown in Fig. 2b, predominantly because the Al2 atoms can displace towards the vacancy sites apart from the Al1 atoms, thereby mitigating cation-cation electrostatic repulsion.

In the saddle-point structures for MEPs in binary and ternary wurtzite-type compounds, cations locate at the trigonal bipyramidal sites ^10,12,26. The basal triangles of the bipyramids do not offer enough space for accommodating the cations, resulting in the expansion of the anion triangles. The resultant in-plane lattice expansion destabilizes the saddle-point structures. In sharp contrast, α-Al₂S₃ has the lowest total energy when half of the Al atoms locate at the bipyramidal sites (i.e., the LP states), highlighting the switching behavior quite different from that for other wurtzite ferroelectrics. Upon transitioning from the + HP to + LP states, the Al2 atoms move to the center of S triangles, accompanying an expansion of the triangle (See Figure S3). In conventional wurtzite ferroelectrics, such triangle expansion increases the in-plane lattice constants typically by 5% ^14,16, destabilizing the h-BN-like saddle-point structures. In α-Al₂S₃, the in-plane lattice constant for the LP state is only 1% larger than that of the HP state (See Figure S4). Figure 3b and d depict the LP structure, in which the AlS₅ and AlS₄ polyhedra tilt with respect to the c axis, resulting in the buckling of close-packed sulfur basal planes, i.e., non-flat sulfur basal planes. The buckling alleviates the in-plane lattice expansion. In conventional wurtzite-type compounds, four anion tetrahedra are connected to each other via an S atom, "locking" the anion polyhedral network. Meanwhile, in α-Al₂S₃, just two or three tetrahedra are connected to each other due to the Al vacancies, leading to a flexibility of the polyhedral network.

Thus, in contrast to conventional wurtzite-type ferroelectrics, the cation vacancies mitigate the electrostatic repulsion between cations in α-Al₂S₃ (Fig. 2b). Also, the disconnection of polyhedral network induced by the Al vacancies produces structural flexibility alleviating the in-plane lattice expansion. These features are considered as the primary reasons why the switching barrier of the defective wurtzite α-Al₂S₃ is much smaller than those of "filled" wurtzite-type compounds.

2.3 Stabilization Mechanism of AlS₅ Trigonal Bipyramidal Coordination

We have discussed above how the Al vacancies stabilize the intermediate structures including the LP states from the structural aspects. The half of Al atoms occupy the trigonal bipyramidal sites in the most stable intermediate LP structures for α-Al₂S₃, whereas, in binary and ternary wurtzite ferroelectrics, cations occupy the trigonal bipyramidal sites in the unstable saddle-point structures. It remains unclear why the LP states are stable in α-Al₂S₃ from the viewpoint of chemical bonding. To elucidate the underlying stabilization mechanism of the AlS₅ bipyramids in terms of chemical bonding, bond valence (BV) and crystal orbital Hamilton population (COHP) between the Al and S atoms were calculated. Here, BV and COHP describe the strength of chemical bonding based on the bond length and the interactions between atomic orbitals, respectively. The negative and positive values of COHP indicate bonding and antibonding interactions, respectively. We utilize negative COHP (–COHP) integrated with respect to energy (–ICOHP) within the valence bands, which indicates the magnitude of net energy gain due to bonding and anti-bonding interactions. Figure 4a labels the S atoms composing of AlS₄ tetrahedra or AlS₅ bipyramids as follows: the axial S atoms located at the + c and –c sides with respect to the Al atoms as ax1 and ax2, respectively, and the equatorial S atoms composing of basal triangles as Eq. 1, Eq. 2, and Eq. 3. Figure 4b and c present the BVs and –ICOHP, respectively, for the Al1 and Al2 atoms against the S atoms for the ZP, +LP, and + HP states.

Transitioning from the + HP to + LP states, the Al2 atoms migrate from the tetrahedral T2u sites to the bipyramidal B1 sites (Fig. 2b), as described above. In the + HP state, the BV sums and –ICOHP of both Al1 and Al2 atoms are primarily contributed to by four S atoms, ax1, Eq. 1, Eq. 2, and Eq. 3, with minimal contribution from ax2, confirming 4-fold tetragonal coordination (Fig. 4b and c). In contrast, in the + LP state, the BV sum and –ICOHP of Al2 atom are significantly contributed to by the five S atoms, corroborating 5-fold coordination rather than 3-fold coordination.

Let us pay attention to the chemical bonding evolution for the Al2 atoms, as it is crucial to understand the stabilization mechanism of the + LP state. Figure 4b reveals that the BV sum of the Al2 atoms is significantly smaller for the + LP state than for the + HP state, indicating that the ionic bonding between Al2 and S atoms become weaker upon transitioning from the + HP to + LP states. As shown in Fig. 4c, however, the sum of –ICOHP for the Al2 atom is comparable in the + LP and + HP states, revealing that the energy gain due to covalent bonding remains upon the structural evolution. This motivates us to unravel the Al2-S covalent bonding that compensates for the poor ionic bonding in more detail.

Figure 4d and e depict the –COHP for the Al2-S bonds as a function of energy in the + LP and + HP states, respectively. A remarkable difference is seen just below the Fermi level; the axial and equatorial S atoms show bonding and antibonding contributions, respectively, in the energy range from − 1 to 0 eV for the + LP state, yielding a net bonding contribution, whereas negligible –COHP is found in this energy range for the + HP state, indicating no bonding contribution. This remarkable difference is considered as a major factor stabilizing the LP states. Figure 4f and g illustrate the wavefunctions of eigenstates just below the Fermi level for the + LP and + HP states, respectively, both of which mainly consist of S 3p_z orbitals. In the eigenstate of the + LP state, the 3p_z orbitals of the ax2 S atoms elongate towards the Al2 atoms to overlap with the Al2 3p_z orbitals, clearly indicating a s-like bonding interaction between these orbitals. Meanwhile, the S 3p_z orbitals of the ax2 S atoms are localized and nonbonding in the eigenstate for the + HP state. The ax2 S atoms are 2-fold coordinated by Al atoms in the + HP state, rendering one of the 3p orbitals nonbonding. The Al2 displacements to the bipyramidal sites yield an additional Al coordination to the ax2 S atoms, causing the ax2 S 3p_z orbitals to take part in the bonding states with the Al 3p_z states. The formation of the bonding states in the LP states is considered as another factor stabilizing the LP state with respect to the HP state.

2.4 Effects of Epitaxial Strain and Chemical Doping on Switching Barriers

For wurtzite ferroelectrics, lowering the switching barrier is a key priority for wider practical applications. It has been demonstrated that epitaxial biaxial strain and chemical doping lower the switching barriers for ZnO and AlN by theory and experiments ^16,17,27. Here, the effects of biaxial strain and chemical doping on MEPs are investigated for Al₂S₃ to give insights into the experimental control of ferroelectric switching.

The reduction of switching barrier in α-Al₂S₃ requires the stabilization of the ZP states and/or the destabilization of the LP states. Without the constraint of strain, the in-plane lattice constant of the LP state is larger than that of the ZP state, as shown in Figure S4, implying that the MEP can be controlled by harnessing biaxial strain. Figure 5a shows calculated MEPs under biaxial strain. The total energies of the LP states become higher with respect to those of the ZP states under compressive (negative) biaxial strain, which is consistent with the in-plane lattice constant of the LP state larger than that of the ZP state (Figure S4). Figure 5b shows the switching barriers under strain, which correspond to the total energy difference between the LP and ZP states. The switching barriers decrease by 8.5% with 2% compressive strain. The biaxial strain dependence of switching barriers for α-Al₂S₃ is contrary to that for wurtzite ferroelectrics such as ZnO and AlN, where the switching barriers are reduced by tensile biaxial strain ¹⁴. This is because the in-plane lattice constant decreases in α-Al₂S₃ and increases in typical wurtzite ferroelectrics when climbing the potential hills for the polarization switching.

Next, we consider the chemical doping effects on the switching barriers. It has been reported that the switching barrier is successfully reduced for AlN by doping B atoms, which favor planar triangle coordination. It is not likely that B-atom doping helps lower the switching barrier for α-Al₂S₃, because it can stabilize the LP states with trigonal bipyramidal coordination. Here, we focus on doping of Ga atoms. Ga is considered as a prime candidate for the following three reasons; (1) Ga atoms are isovalent to Al. (2) α-Ga₂S₃ adopts γ-In₂Se₃ structure similarly to α-Al₂S₃ ²⁸. (3) The ionic radius of Ga³⁺ is bigger than that of Al³⁺. The cation-anion radius ratio of Al to S (r_Al/r_S = 0.212) is less than 0.225, which is the minimum value for 4-fold tetrahedral coordination according to the Pauling's third rule. Meanwhile, r_Ga/r_S is 0.255, indicating that Ga atoms favor tetrahedral coordination of S atoms ²⁹. It is expected from these facts that Ga-doping stabilize the ZP and HP states and destabilize the LP states, reducing the switching barrier.

Figure 5c shows the MEPs for representative structural models of Al_(12-x)/6Ga_x/6S₃ with x = 0, 2, 4, 6, 8, 10, 12. An increase in x stabilizes the ZP and HP states with respect to the LP states, as expected. In the low-x range (0 $\:\text{≤}$ x $\:\text{≤}$ 6), the pathway from the +LP to ZP to –LP states corresponds to the highest potential hill, whereas, in the high-x range (6 $\:\text{≤}$ x $\:\text{≤}$ 12), the +HP-+LP pathway includes the highest potential hill. Figure 5d shows the box-and-whisker plot of the switching barriers against doping concentration x in Al_(12-x)/6Ga_x/6S₃. The switching barrier decreases with an increase in x, and shows a minimum at x = 6, followed by an increase in the barrier above x = 6. The barrier is about 40% smaller for x = 6 compared to pristine α-Al₂S₃. Thus, our calculations predict that Ga-doping facilitates the polarization switching for α-Al₂S₃. Furthermore, Ga-doping modulates the relative energy relationship of the four different polarization states, which enables the control of stable phases.

2.5 Piezoelectric Constants

α-Al₂S₃ was found to be a quite rare example of quadruple-well ferroelectrics. Its polarization is always oriented along the c axis in the MEP with four local energy minima. The four polar states ($\:\text{±}$HP and $\:\text{±}$LP) can be distinguished by PFM if they have distinct piezoelectric constants. Here, we calculated the piezoelectric constants for the +HP and +LP states. Table 1 summarizes the piezoelectric stress constants (e₃₃), elastic stiffness coefficients (C₃₃), and piezoelectric constants (d₃₃) for the + HP and + LP states. Their piezoelectric constants are comparable to that of pristine AlN (~5 pC/N) ^30,31. The piezoelectric constants of the LP and HP states are distinctly different so that these two states are distinguished for the c-plane cleaved single crystal samples by PFM, as in CuInP₂S₆ ²⁵. The four polar states can be detected by PFM since the positively and negatively polarized states show piezoresponse with opposite signs.

Table 1

Piezoelectric stress constants (e₃₃), elastic stiffness coefficients (C₃₃), and piezoelectric constants (d₃₃) for the + HP and + LP states of α-Al₂S₃.
State	e₃₃ (C/m²)	C₃₃ (GPa)	d₃₃ (pC/N)
+HP	0.654	90	7.3
+LP	0.246	65	3.8

We unveiled an unusual ferroelectricity in a defective wurtzite α-Al₂S₃ using first-principles calculations. α-Al₂S₃ is predicted to be a rare example of quadruple-well ferroelectrics, which have four local energy minima in the MEPs. The intermediate lower-polarization states contain 5-fold coordinated AlS₅ trigonal bipyramids. The polarization switching barrier of α-Al₂S₃ (51 meV/cation, 1.0 meV/Å³) is one order of magnitude smaller than that of a typical wurtzite ferroelectric AlN. The Al vacancies alleviate electrostatic repulsion between Al atoms and bring into structural flexibility mitigating elastic energy penalty during the polarization switching. The bonding interactions between Al and S 3p_z states play a role in stabilizing the AlS₅ bipyramidal coordination. Biaxial compressive strain and Ga-atom doping destabilize the intermediate lower-polarization structures, facilitating the polarization switching. In particular, 50% Ga substitution is predicted to reduce the switching barrier by about 40%. Our calculated piezoelectric constants revealed that PFM enables to distinguish the four polarized states. Overall, this study encourages the experimental investigation of an unconventional ferroelectric Al₂S₃ as a new ferroelectric material promising for computing and data storage devices. Notably, the predicted quadruple-well potential surface as well as its fine tunability with chemical doping enables innovative devices such as a multi-valued high-density memory. This work also provides a new strategy for reducing the switching barrier in wurtzite ferroelectrics: introducing cation defects.

4.1 Density functional theory calculations

First-principles calculations were carried out based on density functional theory (DFT). We used the projector augmented-wave (PAW) method^32,33 and the GGA-PBEsol functional^34–36 as implemented in the Vienna Ab-initio Simulation Package (VASP 5.4.4) ^37–40. A plane-wave cutoff energy of 300 eV was used. The radial cutoffs of PAW data sets for Al, Ga, and S are of 1.4, 1.2, and 1.2 Å, respectively. Al 3s, 3p; Ga 3d, 4s, 4p; and S 3s, 3p states are treated as valence electrons. Γ-centered 3 × 3 × 2 k-point mesh sampling was employed. The lattice constants and internal coordinates were optimized until residual stress and forces converged to 0.01 GPa and 1 meV/Å, respectively. Born effective charge tensors and piezoelectric stress tensors were obtained by using density functional perturbation theory (DFPT) calculations. For DFPT calculations, a plane-wave cutoff energy of 600 eV and Γ-centered 6 × 6 × 4 k-point mesh sampling were used.

4.2 Minimum energy pathways

First, the higher symmetry structure of α-Al₂S₃ was searched using the spglib code ⁴¹. The atomic displacements of the polar structure with respect to the higher symmetric structure were obtained using the STRUCTURE RELATIONS code in BCS ^42–44. By reversing the displacements, the polar structural models with polarization in the opposite direction were created. Structural relaxation was performed for the polar end structures. The intermediate images were generated by linear interpolation between the relaxed end structures. To represent a complicated pathway, we employed 32 intermediate images ¹¹. Solid-state nudged elastic band (SS-NEB) method²¹ was employed to determine the switching pathways and barriers using VASP Transition State Theory (VTST) tools (VTST code-198) developed by Henkelman and Jonsson^45,46 except for the calculations under biaxial epitaxial strain. The polarization was calculated for each image using its structure and Born effective charges.

The effects of Ga doping on the switching barriers were also examined. Ga-doped structures were thoroughly searched using the CLUPAN code ⁴⁷. As mentioned above, the polarization switching in α-Al₂S₃ involves site exchange between the T1 and T2 sites. Depending on doped structural models, the initial structures are not equivalent to the final structures, leading to asymmetric MEPs. In this study, for simplicity, we chose the doped structural models for which the initial and final structures are equivalent to each other, and calculated MEPs by SS-NEB to obtain the switching barriers. We used the algorithm implemented in the pymatgen code⁴⁸ to check the consistence of the structures before and after the switching. Structural relaxation was performed for the obtained end structural models. After preparing initial pathways by the linear interpolation of the end structures, MEPs and switching barriers were calculated using SS-NEB methods. The number of data is 6 for x = 2 and 10; 15 for x = 4 and 8; and 14 for x = 6.

To examine the strain dependence of the switching barriers, MEPs were also calculated under the constraint of fixed in-plane lattice constants using NEB methods implemented in the VASP code. We defined in-plane biaxial stain as s = (a – a₀)/a₀, where a₀ is the in-plane lattice constant of the unstrained α-Al₂S₃ structure. The calculations were carried out within the s range from − 2 to 2%. Out-of-plane lattice constants and internal coordinates were optimized for the polar end structures of the switching pathways with fixed in-plane lattice constants. The initial pathways were created by the linear interpolation of the end structures.

4.3 Piezoelectric constants

In this study, we focused on the diagonal component of piezoelectric (strain) constant, d₃₃. The piezoelectric constants were derived from piezoelectric stress tensors and elastic constants according to the procedure described in Ref. ⁴⁹. We calculated elastic constants using the strain-energy relationship. In the Voigt notation, the elastic energy can be written as ^50,51

$$\:E=\:\frac{1}{2}{\epsilon\:}_{p}{C}_{pq}{\epsilon\:}_{p}$$

where C_pq is the elastic constants, and e_p is the strain. Consider the point group 6 for α-Al₂S₃, the elastic constant matrix C is represented as ⁵²

$\:\left(\begin{array}{cccccc}{C}_{11}&\:{C}_{12}&\:{C}_{13}&\:0&\:0&\:0\\\:{C}_{12}&\:{C}_{11}&\:{C}_{13}&\:0&\:0&\:0\\\:{C}_{13}&\:{C}_{13}&\:{C}_{33}&\:0&\:0&\:0\\\:0&\:0&\:0&\:{C}_{44}&\:0&\:0\\\:0&\:0&\:0&\:0&\:{C}_{44}&\:0\\\:0&\:0&\:0&\:0&\:0&\:\frac{1}{2}\left({C}_{11}-{C}_{12}\right)\end{array}\right)$ (2).

The strain tensor $\:\epsilon\:=\left({\epsilon\:}_{1},\:{\epsilon\:}_{2},\:{\epsilon\:}_{3},\:{\epsilon\:}_{4},\:{\epsilon\:}_{5},\:{\epsilon\:}_{6}\right)$ is defined as

$\:\epsilon\:=\left(\begin{array}{ccc}{\epsilon\:}_{1}&\:\frac{{\epsilon\:}_{6}}{2}&\:\frac{{\epsilon\:}_{5}}{2}\\\:\frac{{\epsilon\:}_{6}}{2}&\:{\epsilon\:}_{2}&\:\frac{{\epsilon\:}_{4}}{2}\\\:\frac{{\epsilon\:}_{5}}{2}&\:\frac{{\epsilon\:}_{4}}{2}&\:{\epsilon\:}_{3}\end{array}\right)$ (3).

By applying monoaxial strain $\:\epsilon\:=\left(0,\:0,\:\delta\:,\:0,\:0,\:0\right)$ to the crystal, we obtained the elastic energy:

$\:E=\:\frac{1}{2}{C}_{33}{\delta\:}^{2}$ (4).

Therefore, the elastic constants $\:{C}_{33}$ was determined from total energies versus strain along the c axis.

4.4 Others

We used the LOBSTER code to perform COHP analysis ^53–57. Bond valence was calculated using the bond valence parameters reported by Brese and O'Keeffe ⁵⁸. Phonon band structures were calculated using the PHONOPY code ^59,60. Space group symmetry was determined using the spglib code ⁴¹. The VESTA code was used to visualize crystal structures ⁶¹.

ASSOCIATED CONTENT

Supporting Information.

The Supporting Information is available free of charge at the website. It includes calculated phonon band structures of α-Al₂S₃ in the +HP and +LP states and changes in the area of the basal triangles involved with the Al1 and Al2 atoms and the lattice constants a and c during the polarization switching in α-Al₂S₃ (PDF). An animation of structural evolution during the polarization switching of α-Al₂S₃ is also available (MP4).

AUTHOR CONTRIBUTIONS

Y.S. conducted the calculations under the supervision of H.A. Y.S. and H.A. drafted the manuscript, and all authors edited it.

ACKNOWLEDGMENT

This research was supported by Japan Society of the Promotion of Science (JSPS) KAKENHI Grants Nos. JP17K19172, JP18H01892, JP19H00883, JP21K19027, JP21H05568, JP21H04619, JP23H02069, and JP23H01869. H.A. appreciates Murata Science Foundation and Collaborative Research Project of Laboratory for Materials and Structures, Institute of Innovative Research, Tokyo Institute of Technology. S.O. gratefully acknowledges the Toyota Riken for financial support through a Rising Fellow Program. The computation was carried out using the computer resource offered under the category of General Projects by Research Institute for Information Technology, Kyushu University.

COMPETING INTERESTS

The authors declare that they have no competing interests.

DATA AVAILABILITY

The datasets used and/or analysed during the current study available from the corresponding author on reasonable request.

Kishi, H., Mizuno, Y. & Chazono, H. Base-metal electrode-multilayer ceramic capacitors: Past, present and future perspectives. Japanese Journal of Applied Physics, Part 1: Regular Papers and Short Notes and Review Papers 42, 1–5 (2003).
Setter, N. et al. Ferroelectric thin films: Review of materials, properties, and applications. J Appl Phys 100, 51606 (2006).
Kohli, M. et al. Pyroelectric thin-film sensor array. Sens Actuators A Phys 60, 147–153 (1997).
Fichtner, S., Wolff, N., Lofink, F., Kienle, L. & Wagner, B. AlScN: A III-V semiconductor based ferroelectric. J Appl Phys 125, 114103 (2019).
Hayden, J. et al. Ferroelectricity in boron-substituted aluminum nitride thin films. Phys Rev Mater 5, 044412 (2021).
Wang, D., Wang, P., Wang, B. & Mi, Z. Fully epitaxial ferroelectric ScGaN grown on GaN by molecular beam epitaxy. Appl Phys Lett 119, 111902 (2021).
Ferri, K. et al. Ferroelectrics everywhere: Ferroelectricity in magnesium substituted zinc oxide thin films. J Appl Phys 130, (2021).
Wang, D. et al. Ferroelectric YAlN grown by molecular beam epitaxy. Appl Phys Lett 123, (2023).
Kim, K. H., Karpov, I., Olsson, R. H. & Jariwala, D. Wurtzite and fluorite ferroelectric materials for electronic memory. Nature Nanotechnology 2023 18:5 18, 422–441 (2023).
Moriwake, H. et al. A computational search for wurtzite-structured ferroelectrics with low coercive voltages. APL Mater 8, 121102 (2020).
Lee, C.-W. et al. Emerging Materials and Design Principles for Wurtzite-type Ferroelectrics. (2023) doi:10.26434/CHEMRXIV-2023-HF60W.
Dai, Y. & Wu, M. Covalent-like bondings and abnormal formation of ferroelectric structures in binary ionic salts. Sci Adv 9, (2023).
Moriwake, H. et al. Ferroelectricity in wurtzite structure simple chalcogenide. Appl Phys Lett 104, 242909 (2014).
Konishi, A. et al. Mechanism of polarization switching in wurtzite-structured zinc oxide thin films. Appl Phys Lett 109, 102903 (2016).
Dreyer, C. E., Janotti, A., Van de Walle, C. G. & Vanderbilt, D. Correct implementation of polarization constants in wurtzite materials and impact on III-nitrides. Phys Rev X 6, 021038 (2016).
Liu, Z., Wang, X., Ma, X., Yang, Y. & Wu, D. Doping effects on the ferroelectric properties of wurtzite nitrides. Appl Phys Lett 122, 122901 (2023).
Noor-A-Alam, M., Olszewski, O. Z. & Nolan, M. Ferroelectricity and Large Piezoelectric Response of AlN/ScN Superlattice. ACS Appl Mater Interfaces 11, 20482–20490 (2019).
Lee, C.-W., Yazawa, K., Zakutayev, A., Brennecka, G. L. & Gorai, P. Switching it Up: New Mechanisms Revealed in Wurtzite-type Ferroelectrics. (2023) doi:10.26434/CHEMRXIV-2023-XDP95.
Likforman, A., Carré, D., Hillel, R. & IUCr. Structure cristalline du séléniure d’indium In2Se3. urn:issn:0567-7408 34, 1–5 (1978).
Pfitzner, A. & Lutz, H. D. Redetermination of the Crystal Structure of γ-In2Se3by Twin Crystal X-Ray Method. J Solid State Chem 124, 305–308 (1996).
Sheppard, D., Xiao, P., Chemelewski, W., Johnson, D. D. & Henkelman, G. A generalized solid-state nudged elastic band method. Journal of Chemical Physics 136, (2012).
Eisenmann, B. Crystal structure of a-dialuminium trisulfide, Al2S3. Zeitschrift fur Kristallographie - New Crystal Structures 198, 307–308 (1992).
Krebs, B., Schiemann, A. & Läge, M. Synthese und Kristallstruktur einer Neuen hexagonalen Modifikation von Al2S3 mit fünffach koordiniertem Aluminium. Z Anorg Allg Chem 619, 983–988 (1993).
Baek, S. H. et al. Ferroelastic switching for nanoscale non-volatile magnetoelectric devices. Nature Materials 2010 9:4 9, 309–314 (2010).
Brehm, J. A. et al. Tunable quadruple-well ferroelectric van der Waals crystals. Nature Materials 2019 19:1 19, 43–48 (2019).
Sebastian Calderon, V. et al. Atomic-scale polarization switching in wurtzite ferroelectrics. Science (1979) 380, 1034–1038 (2023).
Clima, S. et al. Strain and ferroelectricity in wurtzite ScxAl1−xN materials. Appl Phys Lett 119, 172905 (2021).
Tomas, A., Pardo, M. P., Guittard, M., Guymont, M. & Famery, R. Determination des structures des formes α et β de Ga2S3 structural determination of α and β Ga2S3. Mater Res Bull 22, 1549–1554 (1987).
Shannon, R. D. Revised effective ionic radii and systematic studies of interatomic distances in halides and chalcogenides. Acta Crystallographica Section A 32, 751–767 (1976).
Lueng, C. M., Chan, H. L. W., Surya, C. & Choy, C. L. Piezoelectric coefficient of aluminum nitride and gallium nitride. J Appl Phys 88, 5360–5363 (2000).
Akiyama, M. et al. Enhancement of Piezoelectric Response in Scandium Aluminum Nitride Alloy Thin Films Prepared by Dual Reactive Cosputtering. Advanced Materials 21, 593–596 (2009).
Blöchl, P. E. Projector augmented-wave method. Phys Rev B 50, 17953 (1994).
Kresse, G. & Joubert, D. From ultrasoft pseudopotentials to the projector augmented-wave method. Phys Rev B 59, 1758 (1999).
Perdew, J. P. et al. Restoring the density-gradient expansion for exchange in solids and surfaces. Phys Rev Lett 100, 136406 (2008).
Perdew, J. P., Burke, K. & Ernzerhof, M. Generalized Gradient Approximation Made Simple [Phys. Rev. Lett. 77, 3865 (1996)]. Phys Rev Lett 78, 1396 (1997).
Perdew, J. P., Burke, K. & Ernzerhof, M. Generalized Gradient Approximation Made Simple. Phys Rev Lett 77, 3865 (1996).
Kresse, G. & Furthmüller, J. Efficient iterative schemes for ab initio total-energy calculations using a plane-wave basis set. Phys Rev B 54, 11169 (1996).
Kresse, G. & Hafner, J. Ab initio molecular dynamics for open-shell transition metals. Phys Rev B 48, 13115 (1993).
Kresse, G. & Hafner, J. Ab initio molecular dynamics for liquid metals. Phys Rev B 47, 558 (1993).
Kresse, G. & Furthmüller, J. Efficiency of ab-initio total energy calculations for metals and semiconductors using a plane-wave basis set. Comput Mater Sci 6, 15–50 (1996).
Togo, A., Shinohara, K. & Tanaka, I. $\texttt{Spglib}$: a software library for crystal symmetry search. (2018).
Aroyo, M. I. et al. Crystallography online: Bilbao crystallographic server. Bulgarian Chemical Communications 43, 183–197 (2011).
Aroyo, M. I. et al. Bilbao Crystallographic Server: I. Databases and crystallographic computing programs. Zeitschrift fur Kristallographie 221, 15–27 (2006).
Aroyo, M. I., Kirov, A., Capillas, C., Perez-Mato, J. M. & Wondratschek, H. Bilbao Crystallographic Server. II. Representations of crystallographic point groups and space groups. urn:issn:0108-7673 62, 115–128 (2006).
Henkelman, G., Uberuaga, B. P. & Jónsson, H. A climbing image nudged elastic band method for finding saddle points and minimum energy paths. J Chem Phys 113, 9901–9904 (2000).
Henkelman, G. & Jónsson, H. Improved tangent estimate in the nudged elastic band method for finding minimum energy paths and saddle points. J Chem Phys 113, 9978–9985 (2000).
Seko, A., Koyama, Y. & Tanaka, I. Cluster expansion method for multicomponent systems based on optimal selection of structures for density-functional theory calculations. Phys Rev B Condens Matter Mater Phys 80, 165122 (2009).
Ong, S. P. et al. Python Materials Genomics (pymatgen): A robust, open-source python library for materials analysis. Comput Mater Sci 68, 314–319 (2013).
Wu, X., Vanderbilt, D. & Hamann, D. R. Systematic treatment of displacements, strains, and electric fields in density-functional perturbation theory. Phys Rev B 72, 035105 (2005).
Sarasamak, K., Limpijumnong, S. & Lambrecht, W. R. L. Pressure-dependent elastic constants and sound velocities of wurtzite SiC, GaN, InN, ZnO, and CdSe, and their relation to the high-pressure phase transition: A first-principles study. Phys Rev B Condens Matter Mater Phys 82, 035201 (2010).
Wright, A. F. Elastic properties of zinc-blende and wurtzite AlN, GaN, and InN. J Appl Phys 82, 2833–2839 (1997).
Newnham, R. E. Properties of Materials: Anisotropy, Symmetry, Structure. (2004) doi:10.1093/OSO/9780198520757.001.0001.
Dronskowski, R. & Blöchl, P. E. Crystal orbital hamilton populations (COHP). Energy-resolved visualization of chemical bonding in solids based on density-functional calculations. Journal of Physical Chemistry 97, 8617–8624 (1993).
Deringer, V. L., Tchougréeff, A. L. & Dronskowski, R. Crystal orbital Hamilton population (COHP) analysis as projected from plane-wave basis sets. Journal of Physical Chemistry A 115, 5461–5466 (2011).
Maintz, S., Deringer, V. L., Tchougréeff, A. L. & Dronskowski, R. Analytic projection from plane-wave and PAW wavefunctions and application to chemical-bonding analysis in solids. J Comput Chem 34, 2557–2567 (2013).
Maintz, S., Deringer, V. L., Tchougréeff, A. L. & Dronskowski, R. LOBSTER: A tool to extract chemical bonding from plane-wave based DFT. J Comput Chem 37, 1030–1035 (2016).
Nelson, R. et al. LOBSTER: Local orbital projections, atomic charges, and chemical-bonding analysis from projector-augmented-wave-based density-functional theory. J Comput Chem 41, 1931–1940 (2020).
Brese, N. E. & O’Keeffe, M. Bond-valence parameters for solids. Acta Crystallographica Section B 47, 192–197 (1991).
Togo, A. First-principles Phonon Calculations with Phonopy and Phono3py. https://doi.org/10.7566/JPSJ.92.012001 92, (2022).
Togo, A., Chaput, L., Tadano, T. & Tanaka, I. Implementation strategies in phonopy and phono3py. Journal of Physics: Condensed Matter 35, 353001 (2023).
Momma, K. & Izumi, F. VESTA 3 for three-dimensional visualization of crystal, volumetric and morphology data. urn:issn:0021-8898 44, 1272–1276 (2011).

(Not answered)

Al2S3switching.mp4
An animation of structural evolution during the polarization switching of α-Al2S3
Al2S3SIShimomuraetal.pdf
Supplemental Material

Download PDF

Editorial decision: revise
03 Sep, 2024
Review #2 received at journal
01 Sep, 2024
Review #3 received at journal
01 Sep, 2024
Review #1 received at journal
18 Aug, 2024
Reviewer #3 agreed at journal
16 Aug, 2024
Reviewer #2 agreed at journal
01 Aug, 2024
Reviewer #1 agreed at journal
31 Jul, 2024
Reviewers invited by journal
27 Jul, 2024
Submission checks completed at journal
11 Jul, 2024
First submitted to journal
10 Jul, 2024
Unknown event
10 Jul, 2024
Editor assigned by journal
03 Jul, 2024

You are reading this latest preprint version

Quadruple-well ferroelectricity and moderate switching barrier in defective wurtzite α-Al2S3: a first-principles study

Status:

Version 1

Abstract

Figures

1. Introduction

2. Result and Discussion

2.1 Crystal Structure of α-Al₂S₃ and Polarization Switching Behavior

2.2 Comparison of Switching Behavior and Barriers in Wurtzite Compounds

2.3 Stabilization Mechanism of AlS₅ Trigonal Bipyramidal Coordination

2.4 Effects of Epitaxial Strain and Chemical Doping on Switching Barriers

2.5 Piezoelectric Constants

3. Conclusions

4. Method

4.1 Density functional theory calculations

4.2 Minimum energy pathways

4.3 Piezoelectric constants

4.4 Others

Declarations

References

Additional Declarations

Supplementary Files

Status:

Version 1

Quadruple-well ferroelectricity and moderate switching barrier in defective wurtzite α-Al2S3: a first-principles study

Status:

Version 1

Abstract

Figures

1. Introduction

2. Result and Discussion

2.1 Crystal Structure of α-Al2S3 and Polarization Switching Behavior

2.2 Comparison of Switching Behavior and Barriers in Wurtzite Compounds

2.3 Stabilization Mechanism of AlS5 Trigonal Bipyramidal Coordination

2.4 Effects of Epitaxial Strain and Chemical Doping on Switching Barriers

2.5 Piezoelectric Constants

3. Conclusions

4. Method

4.1 Density functional theory calculations

4.2 Minimum energy pathways

4.3 Piezoelectric constants

4.4 Others

Declarations

References

Additional Declarations

Supplementary Files

Status:

Version 1

2.1 Crystal Structure of α-Al₂S₃ and Polarization Switching Behavior

2.3 Stabilization Mechanism of AlS₅ Trigonal Bipyramidal Coordination