Revealing inverted chirality of hidden domain wall states in multiband systems without topological transition

Chirality, a fundamental concept from biological molecules to advanced materials, is prevalent in nature. Yet, its intricate behavior in specific topological systems remains poorly understood. Here, we investigate the emergence of hidden chiral domain wall states using a double-chain Su-Schrieffer-Heeger model with interchain coupling specifically designed to break chiral symmetry. Our phase diagram reveals single-gap and double-gap phases based on electronic structure, where transitions occur without topological phase changes. In the single-gap phase, we reproduce chiral domain wall states, akin to chiral solitons in the double-chain model, where chirality is encoded in the spectrum and topological charge pumping. In the double-gap phase, we identify hidden chiral domain wall states exhibiting opposite chirality to the domain wall states in the single-gap phase, where the opposite chirality is confirmed through spectrum inversion and charge pumping as the corresponding domain wall slowly moves. By engineering gap structures, we demonstrate control over hidden chiral domain states. Our findings open avenues to investigate novel topological systems with broken chiral symmetry and potential applications in diverse systems.


INTRODUCTION
Chirality and topology are concepts of great importance that lead to novel physical properties and potential applications in various fields.Examples include topological surface states in topological insulators 1,2 , Majorana fermions in topological superconductors 3,4 , chiral stacking orders in charge density waves [5][6][7] , and topological lasers in photonic systems 8,9 .As prototypical systems, the Su-Schrieffer-Heeger (SSH) 10 and Rice-Mele 11 models exhibit exotic topological properties such as highly robust Jackiw-Rebbi domain wall zero-energy states 12 , charge fractionalization 13 , and spin-charge separation 14 .As a coupled SSH model, the double-chain (DC) model with broken chiral symmetry shows chiral solitons having topological chiral degrees of freedom and Z 4 topological algebraic operation 15 , where the chirality manifests as a spectrum of the chiral soliton and topological charge pumping observed during the adiabatic process as the chiral soliton slowly moving.
Such SSH and Rice-Mele models have been experimentally realized in various physical systems-polyacetylene 10,16 , cold atomic systems 17,18 , artificial electronic lattices 19,20 , photonic systems 21,22 , and acoustic systems 23,24 .Their quantized Berry phases 25 are shown to be consistent with the bulk-boundary correspondence 26,27 .While many coupled SSH chain systems with nontrivial topology have been extensively studied [28][29][30][31] , most possess chiral symmetry, precluding the emergence of chirality as chirality necessitates symmetry breaking.In contrast, the DC model with broken chiral symmetry has been demonstrated in limited physical systems such as selfassembled indium nanowires and artificial atomic chains, exhibiting distinct chiral domain wall states 20,32,33 .Despite being in the same topological class with preserved time-reversal and broken chiral symmetries 15, [34][35][36] , a comprehensive understanding of such multiband systems remains elusive.Therefore, this work endeavors to present a unified framework elucidating the chirality, topology, and bulk-boundary correspondence underlying the emergence of chiral domain wall states in the coupled SSH chain systems with broken chiral symmetry.
Utilizing a representative DC model with interchain coupling whose chiral symmetry is broken, we unveil the emergence of hidden chiral domain wall states possessing inverted chirality even without necessitating any topological phase transition.Our investigation yields a phase diagram revealing single-and double-gap phases depending on the dimerization of each SSH chain and the strength of the interchain coupling.Within the single-gap phase, the chiral domain wall states manifest two localized states akin to chiral solitons observed in the DC model 32 .In the double-gap phase, we observe the emergence of hidden chiral domain wall states characterized by opposite chirality compared to the preexisting domain wall states in the single-gap phase.
We physically verify this opposite chirality through spectrum inversion of the domain wall state and counter-directional charge pumping observed during the adiabatic process as the domain wall state slowly moves.Using the extended two-dimensional effective Hamiltonian corresponding to the adiabatic process and the Berry curvature distribution, we topologically confirm the chirality of hidden chiral domain wall states.Furthermore, by engineering the gap structure via tuning of the interchain coupling, we successfully control the emergence of the hidden chiral domain state.
Our results not only provide insights into the fundamental physics of multiband SSH systems with broken chiral symmetry but also have important implications for the design of novel devices based on chiral domain wall states.

RESULTS AND DISCUSSION
Double-chain model.First, we introduce the DC model consisting of two SSH chains with interchain coupling (Fig. 1a).The interchain coupling acts as a tuning parameter that controls the electronic structure of this model while it was treated as a small perturbation in the previous works 15,32 .Combining two SSH Hamiltonians, we get the Hamiltonian of the DC model: where spin is abbreviated.The superscript (i = 1, 2) represents the upper and lower chains.
n ) denotes a creation (annihilation) operator for the nth site of the ith chain.t indicates the horizontal nearest hopping integral for the ith chain, where t (> 0) and ∆ (i) represent a hopping amplitude in the absence of the dimerization and the energyvalued dimerization displacement of the i-th chain, respectively.α denotes the interchain coupling strength between the lower and upper SSH chains.Since the A and B dimerized states are two degenerate groundstates for each SSH chain, the DC model naturally leads to four degenerate groundstates 32 , which are denoted as AA, AB, BA, and BB states (Fig. 1a).For instance, the AA groundstate is characterized by ∆ (1) = ∆ (2) = δ > 0, while the BB groundstate exhibits FIG. symmetry 15,32 , while the SSH model belongs to the BDI class due to the preserved time-reversal and chiral symmetries.Such chiral symmetry breaking of the DC model provides the realization of the chirality of the nontrivial domain wall states 15,32 .
Figure 1b-e shows the calculated phase map as a function of α/t and δ/t and three representative electronic band structures.Depending on the number of gaps, there are two large distinct regions: single-and double-gap phases (Fig. 1b). Figure 1c shows a single gap between the second and third bands from the bottom, while Fig. 1d,e has an additional gap between the third and fourth bands.Furthermore, the region of the double-gap phase is divided into two subregions depending on whether the additional gap is direct or indirect: the upper gap between the red solid and black dashed curves is indirect (Fig. 1d), while the upper gap becomes direct above the black dashed curve (Fig. 1e).
In the SSH model, a one-dimensional one-band metallic chain at half filling undergoes the Peierls dimerization, which results in a two-band topological insulator 1,37,38 .Similarly, the DC model becomes a four-band insulator after the dimerization 15,32 , which leads to the gap opening between the second and third bands near the Fermi level.On the other hand, the gap-opening mechanism between the third and fourth bands is different due to the strong interchain coupling.As the interchain coupling increases, the energy eigenvalue at k x = 0 of the third band decreases while the energy eigenvalue at k x = π/a 0 (a 0 is the unit cell size) of the fourth band increases (Fig. 1c-e).Such behavior eventually generates another gap between the third and fourth bands when the interchain coupling is larger than the phase boundary (red line in Fig. 1b).Analytically, the phase boundary between single-and double-gap phases is given by . Moreover, the system undergoes an indirect-direct gap transition with increasing interchain coupling.The indirect-direct gap transition boundary (dashed line in Fig. 1b) reads as Additionally, the inversion symmetry protects the gap-closing between the first and second bands at the Brillouin zone boundary regardless of interchain coupling, the details of which are provided in Supplementary Note 1.
Geometric configurations and quantum spectra of chiral domain walls.We now discuss geometric configurations and quantum energy spectra for all possible domain wall states connecting different groundstates.When two of four groundstates are connected, we find only three distinct types of geometric configurations for nontrivial domain wall states (inset in Fig. 2a) due to the equivalence between the same geometric configuration of domain walls 15,32 .To distinguish such nontrivial geometries, we introduce the chirality and denote the AA-BA and AA-AB type configurations as right-chiral (RC) and left-chiral (LC) domain walls and the AA-BB type configuration as an achiral (AC) domain wall, following the notation of the previous works 15, 32 .For such three geometric configurations of nontrivial domain walls, we obtain the energy spectra and local density of states (LDOS) using tight-binding methods for single-and double-gap phases (Fig. 2).
Even though the same geometric domain wall configurations are employed for single-and double-gap phases, the electronic features of the localized domain wall states are quite different from each other.For a single-gap phase, only two in-gap states (denoted as '2' and '3' in Fig. 2a,c) exist as localized domain wall states for all three types of domain walls.Two in-gap states for the AA-BA (AA-AB) geometric configuration are located below (above) the midgap, while two ingap states for the AA-BB geometric configuration are located symmetrically with respect to the midgap.The nontrivial positioning of the localized electronic states of the domain wall results from chiral symmetry breaking 32,39 .This chiral symmetry breaking confers chirality upon the domain wall states, with chirality defined by the spectrum and topological charge pumping.These findings shed light on the interplay between chiral symmetry breaking and localized domain wall states, further enriching our understanding of the electronic behavior within such systems.
For a double-gap phase (Fig. 2b,d), an otherwise hidden in-gap state (denoted as '4') emerges in the upper gap in addition to two in-gap states in the lower gap (denoted as '2' and '3').The two in-gap states of each domain wall in the lower gap appear similar to those of the single-gap phase.
Surprisingly, the in-gap state in the upper gap is located oppositely to those in the lower gap for the right-and left-chiral domain walls.To clearly indicate such spectrum inversion of domain wall states between upper and lower gaps with respect to each midgap, we adopt the term 'chirality inversion'.The chiral inversion also occurs in the achiral AA-BB domain wall even though the ingap states for the lower and upper gaps seem not to be inverted due to the symmetrical positioning of in-gap states with respect to each midgap.Note that the topological meaning of the chirality inversion is the counter-directional charge pumping observed during the adiabatic process as the domain wall state slowly moves, which will be discussed in the next subsection.
The chirality inversion of the in-gap states of domain walls between the upper and lower gaps is more clearly seen if we plot energy spectra as a function of interchain coupling.Figure 3 shows the evolution of spectra of domain wall states with increasing interchain coupling via the singleto double-gap phase transition.Regardless of the interchain coupling, two nontrivial domain wall states always exist in the lower gap.On the other hand, additional domain wall states emerge when the upper gap opens.Even though additional domain wall states are hidden in the region of α < α 1 , the otherwise hidden domain wall states eventually emerge as the upper gap opens in the region of α 1 < α < α 2 .In the region of α ≥ α 2 where the upper gap is direct, the emerging domain wall states maintain their relative energy positions within the gap.Therefore, each domain wall state's chirality (or the spectrum feature) is preserved during the evolution.
Before proceeding further, we briefly discuss the state 1 appearing in Figs.Indirect gaps exist between dashed lines (α 1 < α < α 2 ) while direct gaps exist when α ≥ α 2 .Here, δ/t = 0.3.As the hidden domain wall state approaches either conduction band minima or valence band maxima, pinpointing its exact position becomes challenging.
between the first and second bands at the Brillouin zone boundary, as discussed in the previous subsection (see also Fig. 1).Therefore, the hidden domain wall states labeled 1 cannot manifest as in-gap states.Consequently, we will focus on the other states from now on.However, it is worth noting that the hidden domain wall states labeled 1 can indeed emerge as in-gap states by introducing a symmetry-breaking mechanism that opens up the gap in the system, as demonstrated in Supplementary Figure 1.
Bulk-boundary correspondence.We now investigate the correspondence between the electronic states and the topological properties of the domain walls using Berry phase and Berry curvature via bulk-boundary correspondence 1,2 .In the context of one-dimensional systems, the electronic spectra of a domain wall state are related to the Berry phase difference between two groundstates that the domain wall interpolates 11,40 .
Table I shows the calculated Berry phases for the four groundstates up to the second and third bands from the lowest one.The well-separated electronic bands depicted in Fig. 1c-e  electronic properties and their topological aspects.All Berry phases are quantized as integer multiples of π/2, due to the Z 4 symmetry of the system 15,32 (The mathematical details are provided in Supplementary Notes 2 and 3).The Berry phases up to the second band decrease as 0, −π, −2π, and −3π for AA, BA, BB, and AB groundstates.On the other hand, the Berry phases to the third band increase as 0, π/2, π, and 3π/2 for AA, BA, BB, and AB groundstates.In contrast to the SSH model, some Berry phases are larger than 2π.It is noteworthy that the Berry phase is defined within a range of modulo 4π instead of the conventional 2π reflecting the evolution of the Wannier center and the Z 4 symmetry of the system 15,32 .Such adjustment also accommodates the charge pumping phenomena occurring during the relevant adiabatic process, as discussed below.
To clarify such Berry phases, as shown in Fig. 4, we plot the continuous change of the Berry phase (or the evolution of the Wannier charge center 15,41 ) up to the second and third bands for three types of cyclic adiabatic processes.The three types of cyclic adiabatic processes are generated as the same types of chiral domain walls move from right to left along the one-dimensional chain.
Above all, we will discuss the Berry phase under RC and LC adiabatic processes.During the RC adiabatic process (Fig. 4a), the total Berry phase up to the second band evolves from 0 to 2π while the one up to the third band changes from 0 to −4π.Conversely, the LC adiabatic process shows the opposite change of the Berry phase (Fig. 4b).As a result, two intriguing observations emerge.First, within a given adiabatic process type, the tendencies for the Berry phase up to the second band and the Berry phase up to the third band are diametrically opposite.Second, if we hold the chemical potential fixed, thereby keeping the band occupation constant, the variation tendencies of the Berry phases between the RC and LC adiabatic processes are also in opposition.
Such findings provide compelling topological evidence for the inversion of chirality between the upper and lower gaps, as well as between RC and LC chiral domain wall states.From a physical perspective, this opposite variation in the Berry phase is tantamount to the counter-directional charge pumping observed during the adiabatic process, which will be elaborated upon in the subsequent paragraph.
From the viewpoint of topology, the chirality of the in-gap state of a domain wall is determined by the direction of topological charge pumping under the adiabatic evolution from one groundstate to the other groundstate when the two groundstates are interpolated by the domain wall 15,32,34,42,43 .
If the direction of the topological charge pumping is negative (positive), the electronic states are located below (above) the midgap.Because the local information of such topological charge pumping is encoded in the Berry curvature during the adiabatic process, we can also identify the chirality of the in-gap state of a domain wall using Berry curvature distribution.Notice that the direction of charge pumping (or the moving direction of the Wannier charge center) and the sign of Berry curvature are opposite.See more details in Supplementary Note 2.
First, let us consider the in-gap states in the lower gap.For a right-chiral (left-chiral) domain wall connecting from AA to BA (from AA to AB) groundstates, the sign of Berry curvature distribution up to the second band from the lowest one, as shown in the first (second) panel of Fig. 5b, is positive (negative), and hence, the in-gap states in the lower gap are located below (above) the midgap, as shown in Fig. 3a,b.Thus, the corresponding charge pumping under the RC and LC adiabatic processes are also opposite, as shown in Supplementary Figure 2. Next, let us consider the in-gap states in the upper gap.In this case, the directions of topological pumping are opposite compared to the lower gap case (Supplementary Figure 2), which is consistent with the chirality inversion of chiral domain walls between upper and lower gaps.
Thus, for a right-chiral (left-chiral) domain wall connecting from AA to BA (from AA to AB), the sign of the Berry curvature distribution up to the third band in the first (second) panels of Fig. 5a is negative (positive), and hence the in-gap states in the upper gap are located above (below) the midgap as shown in Fig. 3a,b.
Finally, for the in-gap states of achiral domain walls in both upper and lower gaps, there is no charge pumping because of the zero total Berry curvature during the adiabatic process (AA→BB→AA) as shown in the third panels of Fig. 5a,b, which is consistent with the change of the Berry phase in Fig. 4c.This gives the symmetrically located electronic states, as indicated by the purple lines in Fig. 3c.
Note that the same (opposite) Berry curvature distribution is repeated during the chiral (achiral) adiabatic process due to the system's Z 4 symmetry 15,32 as shown in Fig. 5, which leads to the quantized Berry phases of the groundstates.Furthermore, we find that the quantized Berry phases for the four groundstates in Table I are independent of the interchain coupling and dimerization strength.This strongly implies that both the hidden chiral domain wall state and the in-gap chiral domain wall state possess consistent topological properties independent of the interchain coupling and dimerization.Usually, a topological phase transition occurs when the gap between two energy bands closes and reopens with the change of topological invariants.Surprisingly, our study does not observe such a topological transition even after the upper gap opens, as the topological invariants remain unchanged.Note that the structure of Dirac points remains constant regardless of the interchain coupling in the absence of dimerization, as shown in Fig. 6.This observation highlights the robustness of the topological properties in the double-chain model, independent of the interchain coupling.

CONCLUSION
In summary, we have studied the emergence of the hidden topological domain wall states via gap engineering without requiring any topological phase transition using a representative doublechain SSH model, specifically designed to break the chiral symmetry.By varying the dimerization and interchain coupling, we constructed the phase diagram composed of single-and double-gap phases.For a small interchain coupling, we found the chiral domain wall states with two localized in-gap states in the single gap.However, for a larger interchain coupling, hidden domain wall states emerge, featuring only a single localized state in an additional gap.Intriguingly, the chirality of these emergent domain wall states in the second gap was found to be opposite compared to that of the original domain wall states in the first gap.We validated this chirality inversion through spectrum inversion of the domain wall state and the observation of opposite charge pumping during the adiabatic process observed as the domain wall state moved slowly.The topological confirmation of the chirality of hidden chiral domain wall states was further supported through the analysis of the Berry curvature distribution.
This kind of chirality emergence is plentiful in nature and has many applications such as chirality-dependent light-matter devices 44 , chiral quantum optics 45 , and chiral magnetic domain memory devices 46,47 .Therefore, we expect our theoretical approach can be applied to diverse topological systems such as In/Si(111) 32,48 , Cl vacancies on Cl/Cu(100) 20,33 , and photonic lattices as well as topological laser systems 8,9 .For instance, our model system can be used as a multidigit topological information carrier 39 by engineering the gap structure and Fermi level.We also expect a new type of multi-frequency topological laser, where the topological single-mode lasing frequency can be selectively controlled.

METHODS
The band structures and phase diagrams in Figs. 1 and 6 were studied using the Bloch Hamiltonian of the DC model, which is given by H(∆ (1) , ∆ (2) , k + e ikx a 0 2 + t where t (i) ± = t ± ∆ (i) with energy-valued dimerization ∆ (i) for the i-th chain.To obtain the spectra and LDOS for the RC, LC, and AC chiral domain walls states in Figs. 2   and 3, we used the tight-binding method for the finite system having 4n + 3, 4n + 1, and 4n + 2 atoms with n = 200, respectively.Such boundary conditions ensure the absence of localized edge states at both ends.The dimerization patterns for the domain wall states were simulated using the position-dependent dimerizations and hyperbolic tangent functions: ∆ (i) (x) = ±δ tanh (x/ξ), with ξ being the characteristic width of the domain wall, where ξ = 1.5a 0 in Figs. 2 and 3.
For the Berry phase and Berry curvature in Figs. 4 and 5, we took into account a cyclic adiabatic process of a 1D Hamiltonian, denoted as H(k x , τ ), where τ is time for the adiabatic evolution.
By extending the 1D lattice system into a 2D lattice system, we replace the time evolution with momentum k y in an extra dimension.Then, we constructed the 2D Hamiltonian H 2D (k x , k y ) such that H 2D (k x , k y = 0) = H(k x , τ = 0) and H 2D (k x , k y = 2π) = H(k x , τ = T ), where T represents the period of the corresponding cyclic adiabatic process.Using this 2D Hamiltonian, we have calculated the Berry phase and Berry curvature.Further comprehensive information can be found in Supplementary Notes 2-4.
Correspondence and requests for materials should be addressed to Tae-Hwan Kim or Sangmo Cheon.

FIG. 2 .
FIG. 2. Spectra and local density of state of localized domain wall states.a (b) Electronic spectra and c (d) local density of state (LDOS) of the AA-BA, AA-AB, and AA-BB domain wall states in the singlegap (double-gap) phase.Red, blue, and purple dots denote the AA-BA, AA-AB, and AA-BB domain wall states, respectively, while gray dots indicate the bulk spectra.In the single-gap phase, states 2 and 3 represent the in-gap domain-wall states, while states 1 and 4 do the hidden ones.The state 1 in a and c is hard to distinguish, as discussed in the main text.In the double-gap phase, states 2-4 represent the in-gap domain-wall states, while state 1 does the hidden one.The LDOS plots show that all domain wall states are confined near the centers of the domain walls.The parameter sets (α/t, δ/t) are given by (0.44, 0.3) and (0.7, 0.3) for a,c and b,d, respectively.The inset above a shows the geometric configurations for AA-BA, AA-AB, and AA-BB domain wall states.The insets in b show the close-up of the states 2 and 3.In c,d, the color bar indicates the intensity of the LDOS.

2 and 3 .
The states labeled 1 are also possible hidden domain wall states between the first and second bands.The LDOS maps in Fig.2c,dclearly show the localized feature except for the AA-BA configuration in Fig.2c.In the AA-BA configuration, the domain wall state's spectrum lies too close to the top of the second band for given parameters, making it hard to discern between such a hidden domain wall state and the bulk state.Moreover, the inversion symmetry protects the gap-closing

3 FIG. 3 .
FIG.3.Evolution of energy spectra of domain wall states with respect to the interchain coupling.Energy spectra for a AA-BA, b AA-AB, and c AA-BB domain walls.Solid and dotted lines represent the in-gap domain wall states and hidden domain wall states, respectively.Two domain wall states exist in the lower gap regardless of α/t.When the upper gap opens, a hidden domain wall state appears for all cases.

AAFIG. 4 .
FIG. 4. Evolution of Berry phase under cyclic adiabatic processes.Evolution of Berry phase for a right-chiral, b left-chiral, and c achiral domain walls under the corresponding cyclic adiabatic processes.These processes are generated as the corresponding chiral domain walls move from right to left along the one-dimensional chain.Black (green) color denotes the total Berry phase up to the third (second) bands from the lowest band.Here, (α/t, δ/t) = (0.7, 0.7).

FIG. 5 .
FIG. 5. Berry curvature distributions.Berry curvature distributions up to a third and b second bands under a cyclic adiabatic evolution using extended 2D Hamiltonians, normalized by the maximum absolute magnitude.The color bar indicates the intensity of the normalized Berry curvature.The cyclic adiabatic process is represented in terms of momentum k y through the dimensional extension 15,32 , enabling the calculation of Berry curvature within the two-dimensional Brillouin zone using extended 2D Hamiltonians.The black arrows indicate singular points.Here, (α/t, δ/t) = (0.7, 0.7).

FIG. 6 .
FIG. 6. Band structures without horizontal dimerization.a The same phase diagram in Fig. 1b.b,c Band structures for open circles in a: b for the orange circle and c for the blue circle.In the absence of dimerization, the formation of the four Dirac points is the same for single-and double-gap phases as shown in b and c.The parameters (α/t, δ/t) are (0.2, 0.0) for b and (0.8, 0.0) for c.

TABLE I .
Berry phases up to the 2nd and 3rd bands for groundstates.The Berry phase is defined within a range of modulo 4π instead of the conventional 2π in accordance with the Z 4 property and the