1. Introduction
Ferroelectric materials stand out for their ability to retain an electric polarization even without an external electric field. These materials can exist in two stable polarization states, which are switchable using electric fields [
1]. Within a ferroelectric sample, regions with different polarization orientations, called domains, coexist and are separated by domain walls (DWs)—boundaries between regions with different polarization orientations (
Figure 1). These DWs are essentially ultra-thin interfaces, just a few nanometers thick [
2,
3]. Interestingly, a domain wall can accumulate a net polarization charge, transforming it into a charged domain wall (CDW). This transformation unlocks a fascinating property: significantly enhanced conductivity. This phenomenon, akin to a localized metal–insulator transition [
4,
5], has been observed in various ferroelectric semiconductors like BiFeO
3, BaTiO
3, LiNbO
3, and PbZr
xTi
1−xO
3, in both thin films and crystals [
6,
7,
8,
9,
10]. The enhanced conductivity of CDWs is attributed to the uncompensated polarization charges (positive or negative) residing within them. These polarization charges generate a local electric field that attracts free carriers, boosting the overall conductivity of DWs [
11]. Notably, head-to-head (HH) CDWs exhibit a remarkable ability to amass a large number of carriers, leading to a dramatic increase in conductivity—sometimes by several orders of magnitude. For instance, BaTiO
3 DWs show a staggering 10
9-fold enhancement [
7,
12]. In LiNbO
3, this translates to a conductance 10 orders of magnitude higher than the bulk material, with charge carriers confined to nanoscale sheets [
13]. Another intriguing aspect of DWs is their controllability. By applying an electric field, one can create, erase, or manipulate the position of CDWs within the material [
14,
15]. This precise control, coupled with the dramatic conductivity enhancement at these nanoscale regions, makes charged domain walls highly attractive for developing future electronic devices [
16]. These devices could include junctions, memory elements, diodes, logic gates, and memristors [
17,
18,
19,
20].
The precise mechanisms governing DW conduction remain a topic of active research, with material-specific variations influencing the dominant factors [
21]. Researchers have broadly categorized these mechanisms into intrinsic and extrinsic effects. Intrinsic mechanisms propose that structural changes within the DW create a polarization discontinuity. This discontinuity leads to disruptions in the electrostatic potential, manifesting as band bending or band gap reduction within the DW region [
2,
22,
23]. The type of polarization discontinuity (HH or tail-to-tail (TT)) determines whether electron or hole conduction is favored through band modulation [
24,
25,
26]. The intrinsic mechanism, relying solely on band modulation within the DW, offers potential advantages such as ultrafast response times and inherent resistance to degradation from defects, compared to the extrinsic mechanism. In contrast, extrinsic DW conduction relies on point defects. These charged defects can segregate towards or away from DWs, attracting or repelling charge-compensating carriers (electrons or holes) to maintain electrical neutrality [
9,
27,
28,
29]. However, real materials often exhibit a complex interplay between both intrinsic and extrinsic mechanisms, making it challenging to isolate their respective contributions. Furthermore, the evaluation of DW-specific responses is hindered by the lack of comprehensive data on fundamental electronic properties like carrier concentrations (electron
n and hole
p) and mobilities (
µ) [
30].
Understanding the transport properties of ferroelectric DWs is critical for their application in next-generation devices. These properties, particularly the concentration and mobility of majority charge carriers, remain elusive due to limitations in current experimental techniques [
13,
30]. Quantifying carrier concentration and mobility within DWs remains a significant challenge [
13,
18,
31,
32,
33,
34]. Accurately measuring even the most common parameter, intrinsic DW conductivity, is also hindered by several factors [
35]. In ferroelectric thin film devices, the interface between the film and the electrode forms a Schottky barrier. This barrier complicates the separation of true DW conduction from the contribution of the interface to the overall measured current. This effect becomes more pronounced in thicker films. The current flowing through DWs can vary significantly between devices fabricated with the same ferroelectric material. This inconsistency is attributed to the influence of randomly distributed defects within the film. These defects lead to poor reproducibility and an unreliable ON/OFF current ratio. The defect concentration on the thin film surface not only alters the state of the interfacial barrier but also significantly affects how defects inject and trap charge under an applied electric field. This, in turn, strongly influences the measured DW current. Therefore, theoretical research that sheds light on the concentration and mobility of the carriers within DWs is highly desirable. Such research can complement experimental efforts and pave the way for the development of next-generation ferroelectric devices.
This study presents a groundbreaking approach that leverages first-principles calculations in conjunction with Boltzmann transport theory. This combined method allows us to predict the mobility and concentration of charge carriers within both HH and TT CDWs. Importantly, the simulations applied here capture the intrinsic properties of the materials, eliminating the influence of confounding factors like defects or dopants. By calculating both carrier concentration and mobility, one can comprehensively determine the conductivity of the CDWs. This study represents the first systematic first-principles investigation of conductance in ferroelectric CDWs, offering valuable insights into the fundamental behavior of charge carriers within these materials.
2. Definitions and Methods
First-principles calculation was performed within the framework of density functional theory (DFT) [
36,
37] based on the Vienna ab initio simulation package [
38,
39,
40]. Exchange and correlation effects were treated with the local density approximation (LDA) [
41] to simulate the eight-unit periodic structure containing both HH and TT 180° CDWs in PbTiO
3. Furthermore, the generalized gradient approximations (GGA) [
42] with Perdew–Burke–Ernzerhof (PBE) parameterization was employed as the exchange-correlation potential to calculate the polarization in the ferroelectric domain. Specifically, the valence electron configurations of O (2s
2 2p
4), Ti (3s
2 3p
6 3d
2 4s
2), and Pb (5s
2 5p
6 5d
10 6s
2 6p
2) were considered. A plane-wave cutoff energy of 630 eV was used. As shown in
Figure 1a, a periodically repeated tetragonal (1 × 1 ×
l) supercell with the composition (PbO-TiO
2)
l (
l = 8) was employed to model both DW structures. Initially, an ideal structure was constructed by stacking m unit cells of PbTiO
3 along the [001] direction in a tetragonal ferroelectric configuration. Both sides possessed polarization vectors converging towards the center. Subsequently, a stable structure (
Figure 2a) was obtained by relaxing the domain walls under constrained shape/volume conditions. Periodic boundary conditions were applied to the end faces of the supercell.
Boltzmann transport theory, with the relaxation time approximation [
43], was employed to calculate the electronic transport properties. In semiconductors, conductivity (
σ) is expressed as:
where
e is the elementary charge,
μn (
μp) denotes the mobilities of the electrons (holes), and
n (
p) signifies their concentrations. Carrier concentration was determined using the Fermi–Dirac distribution for the occupation of Kohn–Sham electronic eigenstates at room temperature. The values for
n and
p were calculated from the density of states using the following formulas [
44]:
where
is the unit cell volume,
f0 is the Fermi–Dirac distribution,
u is the chemical potential,
is the electron energy, and
is the density of states (DOS). The Drude model, within the framework of effective mass theory, describes carrier mobility as
, where
τ is the relaxation time during which electrons lose momentum due to scattering from defects, impurities, or vibrations (phonons) [
45], and
is the carrier’s effective mass. Accounting for anisotropy, mobility in the
β direction (
) can be expressed as:
where
is the average relaxation time for the
β direction [
46]. This relationship allows us to calculate mobility upon obtaining the effective mass and average relaxation time. The effective mass of carriers is calculated based on the band structure of the five-atom unit cell of PbTiO
3, incorporating structural distortions and strains extracted from the HH or TTDW regions. Under an external electric field (treated classically), the effective mass is defined by:
, where
α and
β represent reciprocal lattice vector components,
is the dispersion relation for the
ith band, and
is the reduced Planck’s constant. The required derivatives are evaluated numerically using a finite difference method (five-point stencil) [
47]. This method offers high accuracy of O(
h4) and avoids limitations associated with the single parabolic band approximation [
46]. Determining the relaxation time is a more complex process and will be outlined in detail in a separate section.
The probability at which electrons scatter due to imperfections in the material determines their relaxation time. A higher scattering probability translates to a shorter relaxation time and lower mobility. Fermi’s golden rule was employed to calculate the probability of an electron transitioning from an initial state (
) to a final state (
) after a scattering event [
46]:
Here,
i and
j represent band indexes,
,
are electron wavevectors,
ε is the electron energy,
δ is the Dirac delta function, and
is the scattering matrix element. The scattering matrix element (
) describes scattering from state (
) to state (
) by the deviation potential arising from atomic displacement associated with phonons or the perturbation potential caused by defects or impurities. This equation applies to perfectly elastic scattering events [
43] where electrons maintain their energy. For relaxation time in a specific field direction
β (
τβ), the equation is adjusted by a weighting factor that accounts for the scattering angle. This factor is represented by
, where
(
) is the group velocity [
48]. Calculating the scattering matrix element is a key challenge. Let us focus on the dominant scattering mechanism: electron–phonon interactions, specifically with acoustic phonons. In this case, the scattering matrix element depends on the deformation potential (DP) of the material. DP theory assumes a linear relationship between the change in the lattice potential and the volume change caused by thermal motions [
49]. When considering only acoustic phonons, the matrix element for scattering an electron (or hole) from state
to
is expressed as:
, where
. At high temperatures where lattice waves are fully excited, the amplitude of the wave is given by:
, where
is the velocity of the acoustic wave,
m is the total mass of the lattice per unit volume,
kB is the Boltzmann constant, and
T is the temperature. Finally, the average scattering probability is obtained by averaging the squared matrix element [
46,
49]:
, where
is the elastic constant for longitudinal strain in the propagation direction of the longitudinal acoustic (LA) wave (
β) and
ρ is the mass density.
is the DP constant of the
i-th band. The DP constant,
, and the elastic constant,
, are crucial parameters for calculating the average scattering probability.
can be calculated using a total-energy versus a strain approach. Here, the crystal volume is isotropically deformed by ±1% to determine its response to strain.
represents the shift in the band edge (valence band maximum or conduction band minimum) per unit strain [
46,
50,
51,
52].
3. Results and Discussion
Relaxing the eight-unit cell domain wall structure resulted in an unexpected outcome: the most stable configuration lacked the desired HH and TT CDWs and reverted to a centrosymmetric unpolarized structure [
53]. However, by partially constraining the structure, a meta-stable state that contained the targeted CDWs is achieved, as illustrated in
Figure 2a. The meta-stable structure exhibited lattice parameters of a = 3.822 Å, b = 3.822 Å, and c = 33.778 Å for the entire eight-unit supercell. In terms of a single unit cell, these dimensions translate to a = 3.822 Å, b = 3.822 Å, and c = 4.222 Å. Notably, the typical lattice parameters for unconstrained bulk PbTiO
3 are a = 3.88 Å, b = 3.88 Å, and c = 4.29 Å. This comparison reveals that the simulations used required a slight compressive strain (approximately 1.5%) to be applied in all directions to achieve a stable configuration. This suggests that the system is effectively under high hydrostatic pressure. The numbering scheme in
Figure 2a corresponds to the stacking order of atomic layers within the crystal structure. For instance, ‘1’ represents the PbO layer on the left-hand side, and ‘2’ represents the adjacent TiO
2 layer.
Figure 2b presents the calculated layer-by-layer effective polarization profile within the eight-unit PbTiO
3 supercell containing alternating 180° HH and TT CDWs. The Born effective charge method is used to estimate the polarization for each layer. The profile exhibits a remarkably flat polarization plateau within the domain interiors (about 136.5 µC cm
−2). This value exceeds those previously reported (45–92 µC cm
−2) [
53,
54] but remains lower than the giant polarization observed in super-tetragonal PbTiO
3 thin films (236.3 µC cm
−2) [
55]. Although the chosen PBE functional might slightly overestimate polarization [
56], the dominant factor is the compressive strain, consistent with findings that high pressure enhances PbTiO
3’s ferroelectricity through electronic mechanisms, not traditional ionic ones [
57,
58]. The authors propose that pressure alters the balance between long-range Coulomb interactions and short-range electronic effects. Beyond a critical point, PbTiO
3 becomes more ferroelectric to minimize the overlap between the Ti 3d and O 2s orbitals, naturally increasing due to pressure-induced contraction [
57]. At the DWs, the polarization deviates from the constant value due to localized, strong, and non-homogeneous atomic distortions. As shown in
Figure 2b, the HHDW exhibits zero net polarization at the central TiO
2 layer, unlike a previous study where the PbO layer displayed this behavior [
53]. This discrepancy can arise from the differing boundary conditions, with the prior study including a vacuum on one end. For TTDWs, the zero net polarization layer is the edge PbO layer. These non-uniform polarization regions span roughly two-unit cells, allowing estimation of a surprisingly small DW width of only two-unit cells—significantly smaller than previous results [
53]. Specifically, as depicted in
Figure 2, the HHDW encompasses layers 6–10, while the TTDW involves layers 1–3 and 15–16. This predicted thickness aligns well with the prior theoretical prediction of 10
−7 cm [
59] but remains smaller than the six-unit cell DW width reported by Sifuna et al. [
53] The underlying physical mechanism for stabilizing this intrinsic domain wall and achieving such a small width is an intriguing question.
Figure 2b presents the calculated layer-by-layer effective polarization profile within the eight-unit PbTiO
3 supercell containing alternating 180° HH and TT CDWs. The Born effective charge method is used to estimate the polarization for each layer. The profile exhibits a remarkably flat polarization plateau within the domain interiors (about 136.5 µC cm
−2). This value exceeds those previously reported (45–92 µC cm
−2) [
53,
54] but remains lower than the giant polarization observed in super-tetragonal PbTiO
3 thin films (236.3 µC cm
−2) [
55]. Although the chosen PBE functional might slightly overestimate polarization [
56], the dominant factor is the compressive strain, consistent with findings that high pressure enhances PbTiO
3’s ferroelectricity through electronic mechanisms, not traditional ionic ones [
57,
58]. The authors propose that pressure alters the balance between long-range Coulomb interactions and short-range electronic effects. Beyond a critical point, PbTiO
3 becomes more ferroelectric to minimize the overlap between the Ti 3d and O 2s orbitals, naturally increasing due to pressure-induced contraction [
57]. At the DWs, the polarization deviates from the constant value due to localized, strong, and non-homogeneous atomic distortions. As shown in
Figure 2b, the HHDW exhibits zero net polarization at the central TiO
2 layer, unlike a previous study where the PbO layer displayed this behavior [
53]. This discrepancy can arise from the differing boundary conditions, with the prior study including a vacuum on one end. For TTDWs, the zero net polarization layer is the edge PbO layer. These non-uniform polarization regions span roughly two-unit cells, allowing estimation of a surprisingly small DW width of only two-unit cells—significantly smaller than previous results [
53]. Specifically, as depicted in
Figure 2, the HHDW encompasses layers 6–10, while the TTDW involves layers 1–3 and 15–16. This predicted thickness aligns well with the prior theoretical prediction of 10
−7 cm [
59] but remains smaller than the six-unit cell DW width reported by Sifuna et al. [
53] The underlying physical mechanism for stabilizing this intrinsic domain wall and achieving such a small width is an intriguing question.
Figure 2b presents the calculated layer-by-layer effective polarization profile within the eight-unit PbTiO
3 supercell containing alternating 180° HH and TT CDWs. The Born effective charge method is used to estimate the polarization for each layer. The profile exhibits a remarkably flat polarization plateau within the domain interiors (about 136.5 µC cm
−2). This value exceeds those previously reported (45–92 µC cm
−2) [
53,
54] but remains lower than the giant polarization observed in super-tetragonal PbTiO
3 thin films (236.3 µC cm
−2) [
55]. Although the chosen PBE functional might slightly overestimate polarization [
56], the dominant factor is the compressive strain, consistent with findings that high pressure enhances PbTiO
3’s ferroelectricity through electronic mechanisms, not traditional ionic ones [
57,
58]. The authors propose that pressure alters the balance between long-range Coulomb interactions and short-range electronic effects. Beyond a critical point, PbTiO
3 becomes more ferroelectric to minimize the overlap between the Ti 3d and O 2s orbitals, naturally increasing due to pressure-induced contraction [
57]. At the DWs, the polarization deviates from the constant value due to localized, strong, and non-homogeneous atomic distortions. As shown in
Figure 2b, the HHDW exhibits zero net polarization at the central TiO
2 layer, unlike a previous study where the PbO layer displayed this behavior [
53]. This discrepancy can arise from the differing boundary conditions, with the prior study including a vacuum on one end. For TTDWs, the zero net polarization layer is the edge PbO layer. These non-uniform polarization regions span roughly two-unit cells, allowing estimation of a surprisingly small DW width of only two-unit cells—significantly smaller than previous results [
53]. Specifically, as depicted in
Figure 2, the HHDW encompasses layers 6–10, while the TTDW involves layers 1–3 and 15–16. This predicted thickness aligns well with the prior theoretical prediction of 10
−7 cm [
59] but remains smaller than the six-unit cell DW width reported by Sifuna et al. [
53] The underlying physical mechanism for stabilizing this intrinsic domain wall and achieving such a small width is an intriguing question.
Opposing polarization vectors within these 180° DWs generate bound charges, creating a localized region of high electrostatic energy and inherent instability. However,
Figure 2 demonstrates a stable domain structure, suggesting an intrinsic screening mechanism for the bound charges. The current study calculations, absent of defects or dopants, reveal that this screening mechanism in PbTiO
3 DWs involves electronic reconstruction leading to the formation of a two-dimensional electron/hole gas [
53]. The accumulation of electrons or holes at the DWs effectively neutralizes the bound charges, resulting in the formation of thin, quasi-two-dimensional conductive layers. Subsequent analysis of carrier concentration at the domain walls further supports this conclusion. A critical question persists: what mechanism governs the generation of free carriers that perfectly screen the domain wall’s polarization charges?
Figure 3 depicts the layer-resolved partial density of states (PDOS) for the atomic orbitals of Ti and O in the TiO
2 layers and Pb and O in the PbO layers of the eight-unit PbTiO
3 supercell containing HHDWs and TTDWs. The energy scale is shown in electron volts (eVs), with the Fermi level (
EF) set to 0 eV for reference. Consistent with the labeling, the bottom curve corresponds to the first layer (PbO), and the top curve corresponds to the 16th layer (TiO
2). Notably, all layers exhibit a pseudo energy gap, despite variations in the Fermi level position.
Figure 4 and
Figure 5 depict the electronic structure of the HH and TT CDWs in the PbTiO
3 supercell, respectively. They present both the total DOS and the PDOS for each atomic species within the DW region. Below the
EF, the bands in both figures primarily from O 2p orbitals. Above the
EF, these bands arise from the hybridization of O 2p, Ti 3d, and Pb 6p orbitals. The position of the
EF is critical for understanding the electronic properties. In the HH CDWs (
Figure 4), the
EF lies above the bottom of the conduction band, leading to a partially filled conduction band and a degenerate electronic state. This state is further supported by the non-vanishing PDOS at the
EF, indicating the presence of free electrons. These electrons may originate from the TT CDWs, as evidenced by the crossing of the
EF with the PDOS at the valence band top in
Figure 5. This suggests the formation of holes (missing electrons) at the TT CDWs. Conversely, in the TT CDWs (
Figure 5), the
EF falls within the band gap, specifically below the top of the valence band. This positioning also contributes to a degenerate state with holes populating the top of the valence band. These holes seem to originate from the transfer of electrons to the HH CDWs, as observed in
Figure 4. The local metallization is confined to the CDW region. Moving away from the interface, the PDOS of both the conduction and valence bands converge rapidly towards their bulk counterparts and nearly vanish at the Fermi level, indicating a return to semiconductor behavior. Furthermore, the observed rigid shift in the PDOS of each layer suggests the presence of an internal remnant depolarization field. This field acts to electrostatically confine the free charges within the CDW region [
60]. As shown in
Figure 6, the bending of the energy bands can be influenced by this field, playing a crucial role in determining the atomic orbitals involved in the electronic reconstruction [
53].
Building upon the electronic structure analysis (
Figure 3,
Figure 4 and
Figure 5), one can now investigate the distribution and origin of charge carriers within CDWs. Semiconductor theory dictates that the concentration of electrons depends on the DOS near the conduction band minimum, while the hole concentration relies on the DOS near the valence band maximum.
Figure 4 and
Figure 5 offer deeper insights by decomposing the orbital contributions to these electron and hole charges. Analysis reveals that, in the HH CDW (
Figure 4), the electron gas primarily originates from Ti 3d orbitals, with a minor contribution from Pb 6p orbitals. Conversely, the TT CDW (
Figure 5) exhibits a hole gas mainly composed of O 2p orbitals from both the PbO and the TiO
2 layers. These findings support the concept of electronic reconstruction at the CDWs to screen the polarization charges [
16,
30]. The discontinuity in polarization at the DWs acts as the driving force for this reconstruction. Charge transfer between HH and TT CDWs further enhances screening, stabilizing the overall two-CDW configuration. Both CDWs exhibit localization effects due to the remnant electric potential energy within the domains. As depicted in
Figure 6, these potentials (∨ for HH and ∧ for TT) confine the respective carriers (electrons for HH and holes for TT) within potential wells. Within these wells, a stable, degenerate electron (or hole) gas forms, effectively screening the polarization divergences in the corresponding CDW [
61]. This scenario involves a significant amount of free charge populating the band edges.
To quantify the charge transfer between HH and TT CDWs for screening purposes, one can calculate the concentration of free carriers (electrons and holes) within each layer of PbO and TiO
2 in the periodic DW structure (
Figure 7a). Within the framework of band theory, the distribution of free carriers in ferroelectric semiconductors follows the temperature-dependent Fermi–Dirac distribution. As temperature increases, the number of free carriers increases, leading to a rise in intrinsic conductivity. Here, the Fermi–Dirac distribution was employed to account for the occupancy of Kohn–Sham electronic eigenstates, assuming a room temperature of 300 K for both the HH and TT configurations. These calculations indicate that an electron (or hole) concentration of approximately 10
21~10
22 cm
−3 is necessary for effective screening of the HH (TT) walls in the PbTiO
3 supercell. This carrier concentration aligns with previous reports by Vul et al. (10
21 cm
−3) and Sifuna et al. (10
21–10
22 cm
−3) [
53,
59]. To confirm the observed non-local charge neutrality of the DW system, one can integrate the average free charge densities profiles along the (001) direction. This integration provides the total free carrier concentration per unit volume (
n and
p). The HH CDW exhibits an additional electron concentration populating the bottom of the conduction band, quantified as
= 1.16 × 10
22 electrons cm
−3 (
Figure 7a). These electrons seem to originate from the TT CDW, where the integration of the free hole charge profile yields a value of
= 8.69 × 10
21 holes cm
−3. This observation aligns with the global charge neutrality condition, which dictates that
=
, where
represents the volume of the CDWs. In simpler terms, the additional holes at the TT CDWs balance the extra electrons populating the conduction band at the HH CDWs, maintaining overall electrical neutrality [
53].
The electronic transport properties of the CDWs are calculated using Boltzmann transport theory, effective mass theory, and the relaxation time approximation. The results reveal a mobility of 5.09 × 10
−2 cm
2 V
−1 s
−1 for electrons in the HH CDW (
Table 1). This value aligns with the free carrier mobility of 1.49 × 10
−2 cm
2 V
−1 s
−1 reported for the CDW region in BiFeO
3 [
35]. In contrast, the calculated hole mobility for the TT CDW is significantly higher, at 3.09 cm
2 V
−1 s
−1.
Figure 7b presents the calculated conductivities for each layer of the PbTiO
3 supercell containing both HH and TT CDWs. The average electron conductivity at the HH CDW is 9.48 × 10
3 S m
−1, while the average hole conductivity at the TT CDW is 5.05 × 10
5 S m
−1. Notably, the calculated conductivity of the DWs, ranging from 10
3 to 10
5 S m
−1, is significantly higher than the conductivity of the typical semiconductor silicon (4.35 × 10
−4 S m
−1). This finding confirms good enough CDW conduction, although it remains lower than the conductivity of a typical metal like gold (4.11× 10
7 S m
−1). As shown in
Figure 7b, the hole conductance from the TT CDWs is ten times greater than the electron conductance from the HH CDWs. This difference primarily arises from the larger hole mobility observed at the TT CDWs. Similar observations of significantly higher TT CDW conductance compared to HH CDW conductance have been reported for ErMnO
3 CDWs [
24]. A comprehensive analysis reveals that domain wall conductance originates from the accumulation of charge carriers, with carrier concentration playing a crucial role. However, the difference in conductance between HH CDWs and TT CDWs is primarily driven by the disparity in carrier mobility.
Traditionally, CDWs display metallic-like conductivity, but this property is often hampered by inherent instability. This instability arises from the high electrostatic energy associated with bound charges induced by polarization within the CDW [
62]. To achieve stability and a robust domain wall structure, a sufficiently high density of free carriers is necessary to screen the bound charges. However, attaining this density in practical materials often proves challenging, leading to a reliance on defects for DW stabilization [
7]. Stable, naturally occurring ferroelectric CDWs have only been observed in specific improper ferroelectrics like YMnO
3 [
63] and ErMnO
3 [
24]. However, recent studies suggest that CDWs can form in defect-free PbTiO
3 ferroelectrics with sufficiently thick domains [
53,
61]. The current research indicates the formation of stable charged HHDWs and TTDWs in defect-free, nanoscale PbTiO
3 under high pressure. These intrinsic CDWs achieve stability via self-generated free carriers, eliminating the need for external defects [
59,
61,
64]. The system considered here exhibits remarkable stability despite the absence of external defects. Two key factors contribute to this: a lowered carrier generation barrier and the coexistence of opposite-sign carriers. While the overall band gap shape remains largely intact (
Figure 3), a non-zero, but extremely low, DOS persists at the band gap bottom. The near-zero DOS suggests an extremely low concentration of free carriers, resulting in weak metallicity. The ferroelectricity is evidenced by the large polarization calculated in this paper and supported by existing literature [
57,
58]. Notably, ferroelectricity and metallicity can coexist, as observed in ferroelectric metals, with the potential for exciting applications [
65,
66,
67]. Unlike traditional metals with only electrons, PbTiO
3 under high pressure seems possessing both electrons and holes due to the covalent nature of Pb–O and Ti–O bonds (
Figure 4 and
Figure 5). When an electron breaks free from a covalent bond, it leaves behind a hole. Conversely, consider a hypothetical scenario where a high-pressure environment transforms the covalent bond into a metallic bond, as described by the Drude model with its sole reliance on electron carriers. In such a metallic scenario, an accumulation of electrons at the HHDW due to the shielding of its polarization. However, the TTDW would retain its negative polarization without a corresponding positive counterpart (hole). This inherent imbalance in charge distribution signifies an unstable system, contradicting the calculated stability observed in this study. The presence of holes is crucial for maintaining electrostatic balance within the system [
68]. High pressure creates this non-zero DOS in the band gap, indicating weak metallicity, and eliminates the forbidden band, reducing the energy barrier for the generation of electron-hole pairs. This facilitates the creation of free carriers that can screen the bound charges. However, the extremely low DOS limits the overall carrier concentration, resulting in quite weak metallicity. This aligns with the carrier concentration, which matches values reported for relaxed single-DW systems [
53], indicating that relatively high pressure does not significantly increase carrier numbers. This suggests a critical state between semiconducting and metallic behavior, leading to a semiconductor-like system with quite a small number of electrons and holes. The ultra-low DOS in the pseudo-band gap justifies the use of semiconductor theory for analysis. In summary, the stability of the high-pressure defect-free system arises from a combination of preserved band gap character, weak metallicity, and the presence of electron–hole pairs. As experimentally observed, smaller band gaps lead to more stable CDWs [
21,
69]. In this case, the weak metallicity of system used here accelerates the stabilization of small-scale domain walls.
For traditional semiconductors, such as silicon, germanium, and gallium arsenide, the mobility of electrons is greater than the mobility of holes [
68]. However, in both HHDWs and TTDWs, the carrier mobility is larger for the holes than for the electrons, as shown in
Table 1. Let us attribute this behavior to the anisotropic nature of the effective mass for electrons and holes within the domain wall structures. The effective masses of electrons and holes at both DWs are referenced in
Table 2. At the HHDW, the electron effective mass perpendicular to the domain wall plane (out-plane effective mass,
) is significantly larger (four orders of magnitude) compared to the in-plane effective mass (
and
). This substantial increase suggests that electrons encounter significant difficulty moving perpendicular to the domain wall. The underlying reason is that the volume of the unit cell at HHDW shrinks along the out-plane direction due to HH polarization. In the HHDW region, the lattice constant c = 4.042 Å, while in the TTDW region, the lattice constant c = 4.626 Å. It is understandable that the size of the unit cell at the HHDW is significantly smaller than that at the TTDW. This volume reduction leads to a higher probability of electron scattering, consequently increasing the effective mass. Considering all three dimensions, the large
contributes to a higher density-of-states effective mass
for electrons at the HHDW, as shown in
Table 1. A larger
translates to a lower average scattering time (〈
τ〉), ultimately reducing electron mobility. Conversely, holes in the HHDW exhibit a lower effective mass in the out-plane direction compared to the in-plane direction. This allows for quite easier hole movement perpendicular to the domain wall. Additionally, the lower density-of-states effective mass for holes (about 2
me) further enhances their mobility. In contrast to the HHDW region, the TTDW region experiences a volume expansion along the direction perpendicular to the domain wall due to the reversed polarization direction. This creates a more favorable pathway for electron movement in the out-plane direction, where they have a lower effective mass. However, the in-plane effective mass of electrons within the TTDW remains relatively high, resulting in a large density-of-states effective mass (33.79
me) and, consequently, lower mobility. The lattice distortion at the TTDW also introduces anisotropy in the hole’s effective mass, with different values in all three directions. The highest value (
) indicates a greater impediment for hole motion perpendicular to the domain wall compared to other directions. However, the overall density-of-states effective mass for holes remains relatively low, and the effective masses in all directions are not significantly high. This combination leads to a higher mobility for holes compared to electrons in both DWs.
The existence of “strongly correlated” CDWs, like those observed in PbTiO
3, hinges on near-perfect compensation of polarization charges. This compensation leads to the formation of a degenerate, quasi-two-dimensional electron gas exhibiting metallic-like carrier concentrations [
7]. The potential applications for CDWs are vast, but significant hurdles remain in translating this potential into practical devices. One promising approach involves “frustrated poling” techniques to create stable CDWs during device fabrication [
7,
70]. However, careful selection of electrode materials is critical. Here, the contact potential difference between the ferroelectric and metal electrodes plays a crucial role, and it is directly influenced by the work function difference of the materials involved [
59]. Additionally, the inherent coupling between ferroelectric and ferroelastic domain states can be harnessed for the design and creation of stable CDWs [
7]. This study offers valuable insights into achieving stable CDWs. The studied system, under hydrostatic pressure, undergoes a pressure-induced metallization transition. At the critical state of this transition, electron–hole pairs are readily generated, effectively shielding the domain wall charge and promoting its stability. This suggests that hydrostatic pressure can be a powerful tool for generating stable charged intrinsic domain walls in PbTiO
3-like materials. This paper presents a novel approach to generating stable CDWs without relying on defects.