Symmetry-protected solitons and bulk-boundary correspondence in generalized Jackiw–Rebbi models

We investigate the roles of symmetry and bulk-boundary correspondence in characterizing topological edge states in generalized Jackiw–Rebbi (JR) models. We show that time-reversal (T), charge-conjugation (C), parity (P), and discrete internal field rotation (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$Z_n$$\end{document}Zn) symmetries protect and characterize the various types of edge states such as chiral and nonchiral solitons via bulk-boundary correspondence in the presence of the multiple vacua. As two representative models, we consider the JR model composed of a single fermion field having a complex mass and the generalized JR model with two massless but interacting fermion fields. The JR model shows nonchiral solitons with the \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$Z_2$$\end{document}Z2 rotation symmetry, whereas it shows chiral solitons with the broken \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$Z_2$$\end{document}Z2 rotation symmetry. In the generalized JR model, only nonchiral solitons can emerge with only \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$Z_2$$\end{document}Z2 rotation symmetry, whereas both chiral and nonchiral solitons can exist with enhanced \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$Z_4$$\end{document}Z4 rotation symmetry. Moreover, we find that the nonchiral solitons have C, P symmetries while the chiral solitons do not, which can be explained by the symmetry-invariant lines connecting degenerate vacua. Finally, we find the symmetry correspondence between multiply-degenerate global vacua and solitons such that T, C, P symmetries of a soliton inherit from global minima that are connected by the soliton, which provides a novel tool for the characterization of topological solitons.


Introduction
The bulk-boundary correspondence and symmetry play pivotal roles in understanding topological materials.The bulk-boundary correspondence is that those of the bulk govern the topological properties of the edge modes, and it has been confirmed in many topological materials 1,2 .Symmetry classifies and protects the topological properties of such materials.For instance, topological insulators and superconductors are classified by 10-fold Altland-Zirnbauer classes [3][4][5] , based on time-reversal (T ), chargeconjugation (C), and chiral symmetries, and topological crystalline insulators are protected by space-group symmetries [6][7][8][9] .These studies strongly imply that analyzing the bulk and symmetry is important to understand the physical properties of topological edge states.
Chiral edge states are ubiquitous in nature and often protected by topology and symmetry 1,2,39 .Such chiral topological excitations are expected to show promise for information technology due to their robustness against external perturbations [40][41][42] , and had been treated as the hallmark of 2D and 3D topological insulators 1,2 .Recently, topological solitons having chirality or chiral solitons were experimentally realized in a 1D electronic system with Z 4 topological symmetry 35,43 .Chiral switching using such solitons was demonstrated, which implies that chiral and nonchiral solitons are expected to be used as multi-digit information carriers 44 .However, the universal conditions for the emergence of chiral and nonchiral solitons in 1D are not disclosed yet.

arXiv:2108.00240v1 [cond-mat.mes-hall] 31 Jul 2021
The purpose of this Letter is to discover the role of symmetry and bulk-boundary correspondence in characterizing various topological solitons in 1D systems using JR models.In this work, we show that the cooperation of time-reversal, chargeconjugation, parity, and discrete field rotation (Z n ) symmetries protects and identifies the various types of soliton states via bulk-boundary correspondence.We consider two representative models: The JR model, composed of a single fermion field with a complex mass, possesses phase rotation symmetry.The generalized JR model with two massless but interacting fermion fields has a discrete internal field rotation symmetry.The JR model shows charge-conjugate and parity-invariant solitons (equivalently, nonchiral solitons) in the presence of Z 2 field rotation symmetry, otherwise it shows chiral solitons with broken C and P symmetries.The generalized JR model allows more abundant kinds of solitons-nonchiral (NC), right chiral (RC), and left chiral (LC) solitons [Fig.2(d)]-depending on the controllable physical parameters.In the generalized JR model, the Z 2 rotation symmetry protects NC solitons imposing the equivalence between an NC soliton and its anti-NC soliton.The enhanced Z 4 rotation symmetry promotes the symmetry of global minima from Z 2 to Z 4 , which endows the chirality to solitons.Therefore, we reveal that the internal symmetry enriches the variety of topological solitons in 1D systems.Further T , C, P analysis discloses that a pair of chiral solitons forms a particle-antiparticle pair, while a nonchiral soliton is its antiparticle.Finally, we discover the symmetry correspondence between multiply-degenerate global vacua and solitons such that T , C, P symmetries of a soliton inherit from global minima that are connected by the soliton.

Results and Discussion
The low-energy effective theory of the SSH (m z = 0) and the Rice-Mele (RM, m z = 0) models 15,16 is described by the Lagrangian density of the JR model 12 with a fermion mass m z [Fig.1(a,b)]: where ψ(x) and φ (x) are two-component spinor and bose fields, and γ 0 = σ y , γ 1 = −iσ z , γ 5 = σ x .The mass terms can be represented as φ (x) + iγ 5 m z = m(x)e iγ 5 θ (x) , where m(x) = φ 2 (x) + m 2 z and tan θ (x) = m z φ (x) .For a homogeneous bose field, φ (x) = φ 0 , this model has gapped spectra of E = ± k 2 x + φ 2 0 + m 2 z and two degenerate vacua.In the complex plane, the two degenerate vacua are located at θ and π − θ with tan θ = m z φ 0 , which are denoted as A and B (A and B ) for the SSH (RM) model [Fig.1(c)].If the bose field φ (x) becomes a solitonic background field, φ (x) = φ s (x), that interpolates the two degenerate vacua, there appears an isolated soliton mode in the gap.Depending on m z , the emergent soliton modes are different.The soliton modes in the SSH model are zero-energy modes, while the RM type soliton modes are not [Fig.1(a,b)].
The JR model in Eq. ( 1) has the following global rotation symmetry: Because we are interested in the rotation symmetry that exchanges the two global minima, we set θ = π − 2θ .The rotation symmetry in Eq. (2) becomes Z 2 symmetry (equivalently, π symmetry) only when θ = 0 (or m z = 0): We also analyze T , C, and P symmetries.Because the JR model is spinless, T symmetry exists.On the other hand, the Lagrangian density in Eq. ( 1) has C and P symmetries only when m z = 0: Therefore, Z 2 rotation (or π rotation) symmetry protects the C and P symmetries, and vice versa.Moreover, when the bose field is a soliton field, φ (x) = φ s (x), the soliton system-composed of ψ s (x) and φ s (x)-is self-charge conjugate and parityinvariant [10][11][12] .Thus, the JR model shows charge-conjugate and parity-invariant solitons (or, equivalently, nonchiral solitons) in the presence of the enhanced Z 2 rotation symmetry, otherwise it shows chiral solitons with broken C and P symmetries.
The solitons in the JR model can also be distinguished by topology.With T , C, P and Z 2 symmetries, the two minima in the JR model are distinguished by their quantized Berry phases; hence, the topological phase transition between the two minima occurs through the SSH type soliton.With broken C, P, and Z 2 symmetries, in the RM model, the topology of the two minima can be continuously connected due to the fermion mass; thus, the soliton does not involve such phase transitions.Thus, an enhanced rotation symmetry can endow richer topological structures to solitons [Fig.1(d)].
Schematics of atom configurations and band structures for (a) SSH and (b) RM models.φ is a dimerization displacement that corresponds to a bose field φ (x) in JR model.Red and blue circles represent isolated soliton and antisoliton modes in the gap, respectively.(c) The degenerate vacua for SSH and RM models in the order parameter space of (φ , m z ): A and B (A and B ) for the SSH (RM) model.(d) Atland-Zirnbauer class 3 and π rotation symmetry with respect to the fermion mass m z .R π (/ R π ) indicates that the π rotation symmetry is preserved (broken).Now, we consider the generalized JR model composed of two massless fermion fields and two interfield couplings (t 1 , t 2 ).The Lagrangian density is given by where Ψ j and Φ j are a two-component spinor field and a real scalar bose field localized in the j-th wire ( j = 1, 2), respectively.
L B corresponds to the potential of the bose fields, and K is a positive constant.Figure 2 shows a schematic, the total energy profile, and solitons for several t 2 /t 1 values.When t 1 = t 2 , there are two global and two local energy minima, and the atomic configuration in a global minimum has a staggered dimerization ordering.As a realistic example, this case can describe the low-energy effective theory of the quasi-1D systems of polyacetylene chains [45][46][47] .At the critical point t 1 = t 2 , the two local minima transform to global minima resulting in Z 4 symmetric four degenerate global minima.For example, this case can describe the low-energy effective theory of the coupled double Peierls chain model and explains the charge density wave system of indium atomic wires 35,48,49 .The generalized JR model allows various kinds of solitons depending on t 1 and t 2 .When t 1 = t 2 , only nonchiral solitons can exist, whereas chiral solitons can additionally exist at the critical point (t 1 = t 2 ), as shown in Fig. 2(d).The nature of chiral and nonchiral solitons will be explained shortly.We find that the enhanced Z 4 rotation symmetry supports the emergence of chiral and nonchiral solitons.Using Clifford matrices, the Lagrangian density in Eq. ( 6) can be rewritten as follows: where where e iΓ 14 θ = cos θ sin θ − sin θ cos θ and θ = nπ 2 .Similar to the JR model, the Z 2 rotation exchanges two global minima that are located oppositely with respect to the origin in order parameter space [Fig.2(b)].Hence, the Z 2 rotation symmetry implies that an NC soliton and its anti-NC soliton systems are identical, which indicates that the Lagrangian for an NC soliton system satisfies C and P symmetries as discussed later.While the action is invariant under the Z 2 rotation regardless of t 1 and t 2 , the Z 4 rotation symmetry is realized only at the critical point t 1 = t 2 [Sec.S2.6. in Supplementary Information].The Z 4 rotation rotates the energy profile counterclockwise by π/2 and hence the Z 4 symmetry guarantees four degenerate global minima, which indicates that the Z 4 rotation symmetry is necessary for the simultaneous emergence of the chiral and nonchiral solitons.In addition, the Z 4 rotation symmetry supports the equivalence relations among the same type of solitons because it changes a soliton to the soliton rotated counterclockwise by the same angle in order parameter space [Fig.2(d)].
We analyze T , C, P symmetries in the generalized JR model [see details in Sec.S2.5 of Supplementary Information].T symmetry is always present as the JR model is spinless.On the other hand, C and P symmetries exist only along with the  9) has the following symmetries: For the AC line, and for the BD line, In the viewpoint of energetics, only one symmetry-invariant line is realized when t 1 = t 2 , whereas both symmetry-invariant lines are realized at the critical point t 1 = t 2 .This indicates that the Z 4 rotation symmetry gives the coexistence of the two symmetry-invariant lines.Now we unveil the symmetries of soliton systems by explicitly analyzing the Lagrangian density for soliton systems [see details in Sec.S2.5 of Supplementary Information].An NC soliton system has T,C, P symmetries because the bose fields are localized in the symmetry-invariant line (AC or BD line).This indicates that an AC soliton is its antiparticle having no chirality.On the other hand, in RC and LC soliton systems, C and P symmetries are broken because the bose fields deviate from the symmetry-invariant line, which naturally endows the chirality to solitons.However, the charge conjugated and parity transformed fields of an RC soliton system can be described by the fields of an LC soliton system, and vice versa: where i and j in the subscript indicate two global minima that the corresponding soliton field connects, and the unitary matrices U are given by This indicates that the RC and LC solitons are charge-conjugate and parity partner, i.e., they satisfy particle-antiparticle duality.Similarly, in the JR model, particle-antiparticle duality between soliton and antisoliton can be proved [Sec.S1.4. in Supplementary Information].Similar to the bulk-boundary correspondence in topological insulators, we discover symmetry correspondences between multiple global minima and various types of solitons.T , C, P symmetries of a soliton are determined by the symmetries of two global minima connected by the soliton.Thus, the symmetry properties of solitons can be obtained without definitive proof for the soliton system.Figure 3 shows the diagrams of T , C, P symmetries for the JR and generalized JR models.In the SSH soliton system, two global minima are located in a symmetry-invariant line, and hence both two global minima and solitons connecting them have T,C, P symmetries.In the RM soliton system, two global minima do not have C and P symmetries, hence the solitons do not have C and P symmetries.In the generalized JR model with the Z 4 rotation symmetry, each global minima have T,C and P symmetries.However, the forms of C and P transformations are different according to the symmetry-invariant AC or BD lines, as discussed before.When a soliton interpolates two global minima in the AC or BD line, the soliton has T,C, P symmetries, which results in a self-charge-conjugated and parity-invariant NC soliton.When a soliton connects two global minima out of the same symmetry-invariant line, the soliton does not have C and P symmetries, which endows chirality to the RC and LC solitons as shown in Fig. 3(b).
Such C duality between chiral soliton systems and C symmetry in a non-chiral soliton system appear in physical observables such as energy spectra and fermion numbers of isolated soliton modes.In the generalized JR model, two isolated soliton modes for each soliton system appear, and they are the results of fractionalization or hybridization of zero modes [Fig.4].For RC and LC soliton systems, one isolated energy mode is the primary mode, and the other is the induced mode [Fig.4(a,b)].When a soliton field exists in one fermion field resulting in the primary mode, an effective soliton field is induced in the other fermion field due to the interfield couplings, which leads to the emergence of the induced mode.Thus, a zero mode is fractionalized into primary and induced modes.For an NC soliton system, the two isolated energy modes are distinguishable as bonding and antibonding modes regardless of t 1 and t 2 when t 1 t 2 = 0 [Fig.4(c)].They are the result of hybridization between two zero modes.
We calculate the energy spectra and fermion numbers of such isolated soliton modes using the effective one-particle Hamiltonian and Goldstone-Wilczek method 50 .See detailed calculations in Sec.S2. of Supplemental Information.Figure 4(d-f) shows the energy spectra with respect to the strength of the interwire coupling at the critical point (t 1 = t 2 ).The bonding and antibonding modes in an NC soliton system are symmetrically located with respect to E = 0 because of C symmetry [Fig.4(f)].Because the nature of an NC soliton is protected by the Z 2 rotation symmetry and symmetry-invariant lines, the symmetric spectra of an NC soliton are maintained even if t 1 = t 2 [see Fig. S1 in Supplementary Information].However, primary and induced modes in each RC and LC soliton system are not symmetrically located because C symmetry is broken.Instead, the 5/16   The spectra are normalized by the energy gap for each t 1 .The fermion numbers for soliton modes with respect to t 1 for (g) RC, (h) LC, and (i) NC systems.P, I, AB, and B represent the primary, induced, antibonding, and bonding modes, respectively.Analytical results (AR; lines) and numerical results (NR; squares) are calculated at the critical point isolated energy modes in the RC soliton system are located symmetrically with the isolated energy modes in the LC soliton system [Fig.4(d,e)].This indicates that the RC and LC solitons form a charge-conjugate pair and satisfy C duality. Figure 4(g-i) shows the fermion number of each isolated mode with respect to the strength of the interfield coupling.We now consider the Fermi-level at E = 0, i.e., the bonding mode in an NC soliton system and soliton modes in an RC soliton system are filled, while the other modes are empty.Similar to the energy spectra of RC and LC solitons, the fermion numbers of RC and LC solitons are also opposite [Fig.4(g,h)].The fermion number for the bonding and antibonding modes in an NC soliton system are also oppositely positioned, as shown in Fig. 4(i).This behavior is maintained even for t 1 = t 2 case [see Sec.

Conclusion
In summary, we have studied the roles of symmetry and bulk-boundary correspondence in the JR and generalized JR models.We showed that the cooperation between T , C, P, and internal discrete field rotation symmetries characterizes the physical properties of topological solitons in 1D systems.Furthermore, we found bulk-boundary correspondence between the multiply degenerate global minima and solitons such that the symmetry of a soliton inherits the symmetries of two global minima that the soliton interpolates.Our work can be generally applied to 1D systems such as conducting polymers 51 , charge density wave systems 52 , engineered atomic chain systems 21 , photonic crystal systems 53 , ultracold atomic systems 54 , and higher dimensional systems.

S1.2. Energy spectrum
The energy spectrum can be obtained by the one-particle Hamiltonian.For the homogeneous bose field in a minimum φ(x) = φ 0 , the gapped energy spectrum of a JR model is given by E = ± k 2 + φ 2 0 + m 2 z .On the other hand, when the bose field is a soliton field [φ(x) = φ s (x)] that varies spatially and connects the two degenerate minima [φ s (±∞) = ±φ 0 ], there appears the isolated energy mode E s = m z in the gap.For an anti-soliton field [φ(x) = φ s(x) = −φ s (x)], an isolated energy mode appears at E s = −m z .For the SSH model, all isolated energy mode are zero modes.

S1.3. T , C, P symmetries
For the SSH model, the Lagrangian density in Eq. ( 1) has T , C, P symmetries: (ii) Charge conjugation When the bose field is a soliton field, an additional parity transformation for the bose field Pφ s (x)P −1 = φ s (−x) should be accompanied.Therefore, the soliton system composed of ψ s (x) and φ s (x) is self-charge conjugate [2].
For the RM model, C and P symmetries are broken due to the fermion mass m z , while T symmetry is preserved.

S1.4. Duality of JR solitons
We prove that a soliton system composed of ψ s (x) and φ s (x) and the antisoliton system composed of ψ s(x) and φ s(x) in the JR model form a charge conjugation and parity pair.That is, the charge conjugated field (ψ C s ) for a soliton system can be described with the help of the field (ψ s) in the antisoliton system and vice versa: This can be proved by the equation of motion.
Similarly, the parity transformed field for a soliton system can be described by the help of the field in an antisoliton system, and vice versa: The Lagrangian density of the generalized JR model is where Ψ j and Φ j are two-component spinor and real scalar bose fields (j = 1, 2).In this model, we minimally coupled two JR models by adding the inter-field coupling between Dirac fermion fields.Therefore, the potential for the Bose fields ) has a O(2) rotational symmetry.However, the total energy potential has a discrete symmetry (Z 2 or Z 4 ) depending on the interfield couplings, as discussed in the main text.For simplicity, we do not consider the possible interactions between bosons because we focus on the fermion fields and treat the bose fields as static background fields.
Using Clifford matrices, the Lagrangian density in Eq. ( 12) can be rewritten as follows: where Here, Γ a and Γ ab (a, b = 1, • • • 5) are gamma matrices.They are chosen as Here, we calculate the energy spectrum in a global minimum.From Eq. ( 16), we obtain 4 × 4 one-particle Hamiltonian.When the bose fields are in A or C minimum in Fig. 4(c), for instance, the one-particle Hamiltonian can be written as where +Φ 0 (−Φ 0 ) is for A (C) minimum.The energy spectrum is given by which gives the following energy spectrum.
From this energy spectrum, we quickly check that the bands are degenerated when t 1 = t 2 = 0, while the bands are split when t 1 = 0 or t 2 = 0. Similarly, for the B and D minima, the same energy spectrum can be obtained except for the exchange between t 1 and t 2 .
The energy spectrum can also be calculated by using 2×2 effective one-particle Hamiltonian for each JR fields.From Eq. ( 12), we get the following equations of motion.
By inserting Eq. ( 21) into Eq.( 12), the Lagrangian density of fermions can be described as two effective Lagrangian densites For a constant or slowly-varying soliton bose field, the gradient of Φ j is zero or very small.Then, the effective Lagrangian density in Eq. ( 23) can be approximated as From each effective Lagrangian density in Eq. ( 24), the following effective one-particle Hamiltonian H j for the positive frequency is obtained.

S2.2.2. Soliton modes
The isolated energy spectrum of the fermion field in the background soliton bose fields can be easily obtained with the effective Hamiltonians.
As an example, we consider the RC soliton system of which bose fields interpolate from A to B minima as shown in Fig. 2(d), Φ 1 (±∞) = Φ s (±∞) = ∓Φ 0 and Φ 2 (x) = Φ 0 .Since the RC soliton can exist only at the critical point t 1 = t 2 , the effective one-particle Hamiltonian in Eq. ( 25) becomes As discussed in the main text, there are two isolated soliton modes.For an RC soliton system, the primary mode appears mainly in (Ψ 1 , Φ 1 ) while the induced mode does in (Ψ 2 , Φ 2 ).Thus, H 1 and H 2 can be used to calculate the energy spectra of primary and induced modes, respectively.We first calculate the spectrum of the primary mode.When |x| 1, from Eq. ( 26), the energy eigenvalue equation for the primary mode can be approximated as because the amplitude of the isolated soliton mode vanishes as it moves away from the soliton center.By solving Eq.( 27), the energy spectrum and wavefunction of the primary mode are given by where ∆ s (x) = Similarly, the energy spectrum of the induced mode can be obtained as 1 ) using H 2 , where one can also use a proper unitary transformation of ψ I → ψ I = 1 √ 2 (1 − iσ x )ψ I for convenience.In the same way, the energy spectra of soliton modes in LC soliton systems can be obtained.
For NC soliton systems, the energy spectra of the hybridized bonding and antibonding soliton modes can also be calculated regardless of t 1 and t 1 in a similar way.For a more straightforward calculation, one can use the hybridized bonding and antibonding fields.For instance, can be used in the NC soliton system of which bose fields interpolate from A to C minima.The calculated spectra are plotted in Fig. 3 for t 1 = t 2 case and in Fig. S1 for t 1 = t 2 case.

S2.2.3. Lorentz violation terms
Here, we discuss the Lorentz violating terms in the effective Lagrangian density and their effects on the spectra of the soliton modes.When t 1 = t 2 , from Eq. ( 24), the effective Lagrangian density for the jth field is given by This effective Lagrangian has two types of Lorentz violating terms: A chiral mass term which is proportional to 2t 2 1 Φ 3−j γ 5 and two space-time deformation terms which are proportional to 2t 2 1 γ 1 and 2t 2 1 γ 0 , respectively.Even though both can give corrections to the physical quantities, the space-time deformation terms give higher-order corrections and hence can be ignored in the leading order.For example, consider the energy spectrum of the primary mode in an RC soliton system.If we calculate the energy spectrum ignoring the space-time deformation terms, we have E P = − .Energy spectra for bonding and antibonding modes in an NC soliton system with respect to t2/t1.The upper and lower lines represent the energy spectra for antibonding and bonding modes, respectively.The spectra are normalized by the energy gap.Here, t1 = 0.1.When t2/t1 = 0, bonding and antibonding become two degenerate zero modes even in the presence of the interfield coupling, which is consistent with the ladder-SSH model [5].

S2.3. Fermion number
The zero mode in the SSH model has the quantized 1 2 fermion number [1, 3].Similarly, the isolated soliton modes in the generalized JR model have fractional fermion numbers.We calculate the fermion number of each isolated soliton mode using the Goldstone-Wilczek method [6] and effective Lagrangian densities.
As an example, we calculate the fractional fermion number of the isolated modes in the RC soliton system of which bose fields interpolate from A to B minima, i.e, Φ 1 (±∞) = Φ s (±∞) = ∓Φ 0 and Φ 2 (x) = Φ 0 .First, we consider the fermion number of the primary mode, which mainly exists in Ψ 1 .From Eq. ( 29), the effective Lagrangian density for Ψ 1 can be approximated as Here, up to leading order, we ignore space-time deformations because they give the higher order correction similar to the energy spectrum as we discussed in S2.2.3.

2
. Then, the induced current j µ E± (x) and the fermion number Q(E ± ) for each bands E ± are given by Q(E ± ) = dx j µ E± (x), where Hence, the fermion number of the primary mode is given by Here, 1 2 is a factor for the total fermion number conservation between zero interfield coupling case (t 1 = 0) and nontrivial interfield coupling case (t 1 = 0).When t 1 = 0, the primary mode in the RC soliton system is reduced to the zero mode in the SSH model, which has a fractional fermion number − 1 2 .Hence, using the properties of tangent function, the fermion number in Eq. ( 34) is modified to and the form of Clifford matrices Γ i are given in Sec.S2.1.in Supplemental Information.Depending on interfield couplings, this Lagrangian density can have Z 2 rotation symmetry (n = 2; π rotation symmetry) or Z 4 rotation symmetry (n = 1; π/2 rotation symmetry):

Figure 2 .
Figure 2. (a) Schematic for the generalized JR model composed of two JR wires and two interwire couplings (t 1 and t 2 ).In each ith JR wire, the ith fermion and bose fields, Ψ i (t, x) and Φ i (t, x), reside.(b) Total energy profiles, (c) the corresponding lattice models, and (d) diagrams of possible solitons in the order parameter space (Φ 1 , Φ 2 ) for three t 2 /t 1 values.In (d), nonchiral (NC), right chiral (RC), and left chiral (LC) solitons are represented by grey, red and blue arrows, respectively.A, B, C, and D indicate global and local minima.

Figure 3 .
Figure 3. Diagrams for C and P symmetries of global minima and solitons for the (a) JR and (b) generalized JR models.C and P indicate the symmetries, while / C and / P do broken symmetries.The dashed lines represent C and P symmetry-invariant lines.In (a), the φ axis and the symmetry-invariant line are overlapped.The gray arrow represents solitons in the SSH model, and red and blue arrows do in the RM model.A and B (A and B ) are two degenerate vacua for the SSH (RM) model.In (b), the prime in superscript represents that the corresponding symmetry transformation has a different form comparing with the unprimed one.The red, blue, and gray arrows represent right-chiral, left-chiral, and non-chiral solitons, respectively.

Figure 4 .
Figure 4. Schematics for the soliton modes of (a) RC, (b) LC, and (c) NC soliton systems.Energy spectra for soliton modes for (d) RC, (e) LC, and (f) NC systems with respect to t 1 .The spectra are normalized by the energy gap for each t 1 .The fermion numbers for soliton modes with respect to t 1 for (g) RC, (h) LC, and (i) NC systems.P, I, AB, and B represent the primary, induced, antibonding, and bonding modes, respectively.Analytical results (AR; lines) and numerical results (NR; squares) are calculated at the critical point t 1 = t 2 .
S2.3. in Supplementary Information].Similarly, for the JR model, the spectra of isolated soliton and antisoliton modes show the C duality and C symmetry [Fig.1(a,b)].Therefore, chiral soliton shows C duality, while nonchiral solitons do C symmetry.

2t 2 1 Φ0t 1 EFig. S1
Fig.S1.Energy spectra for bonding and antibonding modes in an NC soliton system with respect to t2/t1.The upper and lower lines represent the energy spectra for antibonding and bonding modes, respectively.The spectra are normalized by the energy gap.Here, t1 = 0.1.When t2/t1 = 0, bonding and antibonding become two degenerate zero modes even in the presence of the interfield coupling, which is consistent with the ladder-SSH model[5].