1 Introduction
Ferroelectric (FE) materials with the emergence of spontaneous electric polarization, have a wide range of technological applications, such as non-volatile memories, field effect transistors, sensors, and photovoltaics [
1−
5]. However, miniaturizing FE devices suffers from size limitation of the FE materials due to internal depolarizing electrostatic fields [
6]. To overcome this issue, significant progress has been made in two-dimensional (2D) ferroelectrics. The ferroelectric polarization in 2D materials of atomic thickness can be strong enough to subsist, offering a natural solution for nanoscale ferroelectric applications [
7,
8].
The lateral dimensionality of 2D ferroelectrics is formally infinite. In reality, however, the ferroelectricity may be weakened by edge reconstruction and boundary [
9]. Alternatively, one-dimensional (1D) materials with “no-boundary” feature might be more appropriate for ferroelectrics. So far, several 1D FE nanostructures have been prepared in the laboratory, such as tetrabromo-p-benzoquinone chain [
10,
11], water molecules in quasi-1D beryl crystal [
12], BaTiO
3 and PbTiO
3 nanotubes [
13], which, however, intrinsically contain nonnegligible structural defects. Therefore, it is still desirable to explore 1D ferroelectrics with structurally perfect “boundary” and high stability.
Atomically precise nanoclusters offer new opportunities for exploring excellent properties and can serve as building blocks with tunable electronic structures for assembling materials and devices [
14]. Among them, “superatomic clusters” (namely, “superatoms”) possess outstanding stability and are promising for self-assembling [
15,
16]. Indeed, some superatomic cluster assembled nanomaterials have already been synthesized. For instance, Re
6Se
8 clusters can be assembled into a 2D hierarchical semiconductor via strong interaction by substitutionally labile Cl atoms, whose transport properties are tunable from semiconducting to metallic and to superconducting upon
n-type doping through Cl dissociation [
17−
19]. C
60 together with Co
6Se
8/Cr
6Te
8/Ni
6Te
8 clusters were patterned into three-dimensional crystalline array with CdI
2 or NaCl type [
20]. On the theoretical side, Guo
et al. [
21] proposed a strategy to construct stable 2D lattices using superatomic tetrahedral clusters (i.e., B
4, Al
4, Ga
4, In
4, Cr
4, Mo
4, and Te
4) covalently linked by oxygen atoms. Du
et al. [
22] examined 27 dimers of endohedral cage clusters and predicted that M@X
16 (M = Ti, Zr, or Hf; X = Si, Ge, or Sn) are potential zero-dimensional vdW Lego blocks. Xing
et al. [
23,
24] presented a joint experiment-theory investigation of a rationally designed layered superatomic crystal of Au
6Te
12Se
8 cubes stacked by noncovalent intercube quasibonds, finding a sequential-emerged anisotropic triple-cube charge density wave and polarized metallic states below 120 K. Zhao
et al. [
25] designed a series of fullerene-based 1D chains, namely U
2C@C
80-M (M = Cr, Mn, Mo, and Ru) 1D chains, with both ferroelectric and ferromagnetic properties.
ReNCl
4 cluster with nearly ideal
C4v symmetry has been synthesized in a facile synthesis from perrhenate, sodium azide and HCl gas [
26−
28]. Inspired by the aforementioned cluster assemblies with extraordinary properties, herein we design 1D nanowires of ReNX
4 (X = F, Cl, Br, and I) by cluster self-assembling and extend them to 1D nanowires of MF
5 (M = V, Nb, and Ta) by isoelectronic substitution. The cluster-assembling rule takes into account both the closed electronic shells for ReNX
4 and MF
5 clusters as well as the binding strength between metal and N/F atoms. The molecular orbitals and superatomic characters of these clusters with 40 valence electrons were identified for the first time. Owing to second-order Jahn−Teller effect, these assembled nanowires intrinsically hold superior ferroelectricity with Curie temperature above 300 K, thereby paving a new route for the design of intrinsic 1D ferroelectrics from self-assembling.
2 Computational methods
Density functional theory (DFT) calculations were performed using the Vienna
ab initio simulation package (VASP) [
29], with planewave basis set (energy cutoff of 500 eV), projector augmented wave potentials [
30], and Perdew‒Burke‒Ernzerhof (PBE) exchange–correlation functional [
31]. Based on equilibrium structures, a hybrid HSE06 functional [
32] was further used to compute the electronic band structures. The convergence criteria for total energy and force were set to 10
−7 eV and 0.01 eV·Å
−1, respectively. Uniform
k-point meshes with spacing of ~0.015 Å
−1 were adopted to sample the Brillouin zones. Phonon dispersion analysis was performed using the Phonopy code [
33] based on the density functional perturbation theory. The frontier molecular orbitals and natural population of individual ReNX
4 and MF
5 clusters were analyzed using the PBE functional accompanied with 6-311+G(d) [
34] and SDD [
35] basis sets, implemented in Gaussian16 suite [
36].
3 Results and discussion
We firstly investigate the atomic structure and electronic properties of an individual ReNCl
4 cluster. As shown in Fig.1(a), four equivalent Cl
– atoms and one N
3– atom are bonded with one Re
7+ atom to form a stable N–Re–Cl
4 unit with
C4v symmetry, with N–Re and Re–Cl bond length of 1.67 and 2.30 Å, respectively, and Cl–Re–N angle of 81.24°. According to the valence electronic configurations of Re (5
d56
s2), N (2
s22
p3) and Cl (3
s22
p5) atoms, each ReNCl
4 cluster carry totally 40 valence electrons and achieves a closed electronic shell of 1S
21P
61D
102S
21F
142P
6, which is evidently demonstrated in Fig.1(b). For the first time, the superatomic characteristics [
37] for valence molecular orbitals of ReNCl
4 cluster are identified, suggesting its great potential as a building block to assemble nanomaterial.
According to the geometric feature of ReNCl
4 cluster, Re atom is prone to bond with another N atom in adjacent ReNCl
4 cluster to form distorted corner-sharing ReN
2Cl
4 octahedra without tilt. In this manner, ReNCl
4 units can be self-assembled to construct a ReNCl
4 wire along the
z axial direction with
P4 symmetry and 1D lattice parameter of 4.08 Å. Here, each N
3– ion is bonded in a distorted linear geometry to two Re
7+ atoms and each Cl
− ion is bonded to one Re
7+ atom with bond length of 2.31 Å, resulting in invariable valence states of Re, N and Cl atoms as confirmed via the analysis of on-site Bader charge [
38] (Re: ‒2.1
e; N: 0.9
e; Cl: 0.3
e). There are one shorter (1.69 Å) and one longer (2.30 Å) Re–N bond with bond overlap population [
39] of 0.93 and 0.11, corresponding to nature of strong and weak covalent bond respectively.
To verify the energetic stability of ReNCl4 nanowire, we calculate its formation energy ΔH defined as ΔH = Ewire – Ecluster, where Ewire and Ecluster are the energy for one unit of 1D nanowire and an individual cluster, respectively. ReNCl4 nanowire shows negative ΔH of −0.22 eV per formula (Tab.1), indicating that its formation is exothermic. The kinetic stability of ReNCl4 nanowire is also confirmed by the phonon dispersions in Fig.2(b).
In light of the stability of self-assembled ReNCl4 nanowire, Cl atoms can be substituted by other halogens (X = F, Br and I) and a series of hypothetic 1D ReNX4 compounds might possess analogous structure as ReNCl4 nanowire. More importantly, all these ReNX4 clusters exhibit superatomic characters. Additionally, isoelectronic substitution is another effective strategy to design isostructural materials. Note that MX5 (M = V, Nb and Ta; X = Cl, F, Br and I) hold same number of valence electrons with that of ReNCl4 (40e), resulting in similar geometric structure for both types of clusters, that is, MX5 cluster also displays X–M–X4 arrangement like ReNCl4 cluster. Therefore, 1D distorted MX5 nanowires can be also constructed by self-assembling of MX5 clusters (see Fig.2). The energetic and dynamical stabilities of these predicted nanowires are systematically assessed by their formation energies and phonon dispersions [Fig. S1 of the Electronic Supplementary Materials (ESM)]. Eventually, we screen out ReNF4, ReNCl4, ReNBr4, ReNI4, VF5, NbF5 and TaF5 nanowires with exothermic formation energies ranging from ‒0.01 to ‒0.43 eV per formula, suggesting experimental feasibility for the synthesis of these 1D structures via cluster self-assembling. These results also support that self-assembling superatomic clusters is beneficial for thermodynamic stability. Specifically, the lattice parameters of ReNX4 nanowires in the range of 3.88−4.10 Å are smaller than 4.11−4.63 Å for MF5 nanowires (Tab.1). In addition, the lattice constant of ReNX4 (MF5) nanowire decreases with increasing atomic number of X (M) element. In all situations, metal ions show an off-center displacement along z axial direction by 0.20−0.82 Å. The electronic band structures of ReNX4 and MF5 nanowires are depicted in Fig. S2 of the ESM. Expect ReNI4, all assembled nanowires are semiconductors with moderate or wide band gaps, which allows switching of the polarization a ferroelectric by applying an external field.
Due to non-centrosymmetry of ReNX4 and MF5 nanowires, there is spontaneous polarization along z axial direction. In order to explore the ferroelectric feature of these 1D nanowires, the centrosymmetric phase (P4/m symmetry) without displacement between metal and halide atoms along z axis is firstly identified as the paraelectric (PE) configuration (Fig.2). As a representative, Fig.2(c) displays the phonon dispersion of 1D PE ReNCl4. Here, three imaginary branches including acoustic and optical phonon modes are throughout the entire Brillouin zone, indicating the geometric instability of PE configuration and a tendency of phase transition from ferroelectric state to a paraelectric state below Curie temperature (TC). In addition, PE phase has higher energy than that of FE phase by 0.09‒0.85 eV per formula. Taking ReNCl4 and VF5 nanowires as examples, the scanned free-energy contour is plotted with different values of d2‒d1 and d3 in Fig.3 (d1 and d2 represent the bond length of metal-linked atoms, and d3 represents displacements of metal and halogen atoms along z chain direction). FE and ‒FE phases are equivalent states with lowest energy, separated by the nonpolar PE state (d2 ‒ d1 = d3 = 0) located at the diagonal direction. These indicate the higher stability of FE and ‒FE phases compared with PE phase. The energy maps of other nanowires are presented in Fig. S3 of the ESM, showing similar phenomena.
The central symmetry of 1D FE ReNX
4 and MF
5 nanowires is eliminated due to the cooperative displacements between metal and halogen atoms, resulting in an evident spontaneous polarization
PFE along the chain direction. According to Berry phase through the modern theory of polarization [
40], the calculated
PFE of 1D FE ReNX
4 and MF
5 nanowires are (25.64‒52.03)×10
−10 C/m and (24.64‒55.50)×10
−10 C/m, respectively, which are comparable and even larger than the values of many FE materials, such as transition-metal dichalcogenides [
41], group-IV monochalcogenides [
42,
43], bulk BaTiO
3 [
44,
45] and lead-zirconate-titanate (PZT) [
46,
47].
We further explore the origin of ferroelectricity in these nanowires. 1D ReNX
4 and MF
5 are semiconductors with band gaps in range of 1.43‒1.25 eV and 4.94−6.60 eV, respectively, and display similar features of electronic structure (Fig. S2 of the ESM). Note that ReNF
4 and VF
5 possess direct band gaps. Taking 1D FE ReNCl
4 nanowire as representative, the projected density of states for each element are displayed in Fig.4. The hybridization between Re‒
dxz/
dyz, N‒
pz and Cl‒
pz near the Fermi level is originated from Re–Cl and Re–N bonding. Here, the cooperation of empty
d orbitals for Re
5+ and
p orbitals for N
3− along the polar direction regulates the existence of FE states in 1D ReNCl
4, resulting in the second-order Jahn‒Teller effect aroused by
d0-configuration [
48,
49], which is similar with the well-known FE perovskite oxides [
44,
50]. Analogously, the ferroelectricity of VF
5 nanowire also results from
d0 principle generated by V
5+ empty
d orbitals. Furthermore, more remarkable
pdπ interaction of Re−
dxz/yz/N–
px/y than
pdσ interaction signifies the strong ferroelectric effect, which is also appeared in BaTiO
3 and PbTiO
3 [
44]. Fig.4(c, d) and Fig. S4 of the ESM present the –COOP (Crystal Orbital Overlap Population) integrals for Re−N and M−F bonds, where a positive –COOP denotes favorable bonding interaction and a negative –COOP indicates unfavorable antibonding interaction. The results reveal that the occupied states are bonding states, while the unoccupied states belong to antibonding ones, further confirming that the FE states is stabilized by Re−N and M−F bonding.
Within the Landau−Ginzburg theory [
51], the free energy and polarization should satisfy
. The first two terms are related to the energy contribution to polarization in each unit cell, which can be described by anharmonic double-well shape of the potential profile.
α and
β parameters are obtained by fitting the double-well potentials in Fig.5 and Fig. S5 of the ESM. The third term reflects the nearest-neighbor dipole−dipole coupling within the nanowire, and coefficient γ is fitted using the mean-field theory within nearest-neighbor approximation, which can be plotted as a parabola shape in second-order approximation shown in Fig.5(b) and Fig. S6 of the ESM. The energy barrier
EB for the transition from ferroelectric to paraelectric phases of ReNX
4 and MF
5 nanowires can be determined as 0.74‒0.85 eV and 0.09‒0.50 eV, respectively. It is noteworthy that
EB is one order of magnitude larger than that of other frequently used ferroelectrics such as BaTiO
3 and PbTiO
3 [
44], suggesting high Curie temperature for these ferroelectric nanowires.
Ab initio molecular dynamic (AIMD) simulations [
52] with 1×1×8 supercells are also performed to evaluate the Curie temperature of these nanowires. Fig.5(c) and Fig. S7 of the ESM show
PFE as a function of temperature. Taking ReNCl
4 and VF
5 nanowires as representatives, the static
PFE values of 25.64 C/m and 55.50 C/m are maintained at finite temperature up to about 300 and 100 K, respectively. Then the Curie temperature for transition of FE to PE phases can be fitted as the temperature-dependent
P(
T) =
a/(1+
), where
a and
k are constants. Importantly,
TC value is usually related to the energy barrier
EB and the dipole‒dipole coupling between the adjacent unit cells, that is,
TC increases with increasing
EB. Here, ReNX
4 nanowires with larger
EB results in higher
TC value (224−321 K) compared with MF
5 nanowires (72−95 K), which also can be obtained from the pyroelectric response defined as
p = ‒(d
P/d
T). Moreover, the piezoelectric response
e11 is dependent on a linear fitting of 1D polarization changes with respect to the uniaxial strains
ε11 along the
z direction. A linear relation between polarization change and the strain within ±1% for ReNX
4 and MF
5 nanowires has been demonstrated in Fig.5(d) and Fig. S8 of the ESM. The obtained
e11 is as high as 61×10
−10 C/m, which is an order of magnitude larger than that of MoS
2 monolayer (3.64×10
−10 C/m) [
53]. The active pyroelectric and piezoelectric responses in these nanowires suggest their broad applications in sensors and energy conversion devices. Particularly, the giant pyroelectric coefficients of ReNX
4 and MF
5 nanowires in range of (0.02‒0.13)×10
−10 C/(K·m
2) (Tab.1) are larger than those of most ferroelectric crystals, for instance, 2.0, 8.3, and 1.8×10
−4 C/(K·m
2) for BaTiO
3, LiNbO
3 and LiTaO
3, respectively [
54].
4 Conclusion
In summary, we have designed a series of ReNX4 (X = F, Cl, Br and I) and MF5 (M = V, Nb and Ta) nanowires by self-assembling superatomic clusters, based on the criteria of superatomic orbital characters within the shell model and orbital interaction between adjacent clusters. Self-assembling superatomic clusters can strengthen the stability of these 1D structures and achieve ferroelectric effect. These cluster-self-assembled 1D materials are all dynamically stable and possess a moderate or wide band gap of 1.43‒6.60 eV (expect for metallic ReNI4 nanowire). Notably, ReNX4 and MF5 nanowires exhibit sizable spontaneous polarization and low switching barrier, and ReNX4 nanowires hold high TC above room temperature that is essential for ferroelectric devices. The extraordinary ferroelectric properties in these assembled materials are derived from an electronic mechanism, i.e., the second-order Jahn‒Teller effect with d0-configuration transition-metal ions in the materials. These theoretical results open unprecedented opportunities for the discovery of new ferroelectric materials and point to a feasible route for the exploration of exotic 1D physics from cluster assemblies.