Quantized topological pumping of solitons in nonlinear photonics and ultracold atomic mixtures

Exploring the interplay between topological band structures and tunable nonlinearities has become possible with the development of synthetic lattice systems. In this emerging field of nonlinear topological physics, an experiment revealed the quantized motion of solitons in Thouless pumps and suggested that this phenomenon was dictated by the Chern number of the band from which solitons emanate. Here, we elucidate the origin of this nonlinear topological effect, by showing that the motion of solitons is established by the quantized displacement of the underlying Wannier functions. Our general theoretical approach, which fully clarifies the central role of the Chern number in solitonic pumps, provides a framework for describing the topological transport of nonlinear excitations in a broad class of physical systems. Exploiting this interdisciplinarity, we introduce an interaction-induced topological pump for ultracold atomic mixtures, where solitons of impurity atoms experience a quantized drift resulting from genuine interaction processes with their environment.

Quantized responses have been a central theme throughout the realm of topological physics, which was initiated with the discovery of the quantum Hall effects in two-dimensional electron gases 1,2 . A wide variety of topological band structures have been revealed over the last decades, leading to the identification of various forms of quantized responses, from quantized Faraday and Kerr rotations in threedimensional topological insulators 3 to quantized circular dichroism 4 and topological Bloch oscillations 5,6 in two-dimensional ultracold atomic gases. An emblematic and minimal instance of quantized topological transport concerns the adiabatic motion of a quantum particle moving in a slowly-varying periodic potential, an effect known as the Thouless pump 7 . In this setting, the center-of-mass motion is quantized according to the Chern number of the underlying band structure, as defined over a hybrid momentum-time space 8 . The realization of synthetic lattice systems has allowed for the experimental implementation of Thouless pumps and for the observation of the related quantized motion, in both photonics [9][10][11][12][13][14] and ultracold gases setups 15,16 .
Interestingly, synthetic topological systems 17,18 can operate beyond the linear regime of the Schrödinger equation, hence opening the door to nonlinear topological physics 19 . In this emerging framework, a central topic concerns the possible interplay between nonlinear excitations, known as solitons, and the underlying topological band structure [20][21][22][23][24][25][26][27][28][29][30][31] . Interestingly, exact correspondences between topological indices and nonlinear modes have been identified in mechanical systems 32 and for the Korteweg-de-Vries equation of fluid dynamics 33 , hence allowing for a formal topological classification of nonlinear excitations in certain special cases. In the context of nonlinear topological photonics, a recent experimental study reported on the quantized motion of solitons in a lattice system undergoing a Thouless pump sequence 34 . Despite the presence of considerable nonlinearity, these observations suggest that the quantization of the solitons displacement is dictated by the Chern numbers of the underlying band structure. This nonlinear topological effect is particularly intriguing: On the one hand, the soliton is associated with a non-uniform occupation of a Bloch band, hence violating the standard condition for quantized motion in conventional (linear) Thouless pumps 7,8 . On the other hand, the nonlinearity prevents the dispersive evolution of a localized wave packet through the formation of a soliton. As already pointed out in ref. 34, these facts highlight the crucial role of nonlinearities in this emerging context of soliton Thouless pumps. In this work, we elucidate and explore the quantized transport of solitons in nonlinear topological Thouless pumps. Inspired by the experiment of ref. 34, we address this topic by considering a general class of one-dimensional systems described by the discrete nonlinear Schrödinger equation (DNLS) 35,36 . Following ref. 37, we represent the solitons in terms of maximally localized Wannier functions of the corresponding Bloch band from which the soliton bifurcates. In the so obtained Wannier representation, the adiabatic motion of the soliton can be deduced from an ordinary (scalar) DNLS; from this, we show that the quantized motion of the soliton is directly related to the quantized displacement of Wannier functions upon pumping, which is known to be set by the Chern number of the band 15,38,39 . This general approach allows us to mathematically demonstrate the topological nature of nonlinear Thouless pumps, by relating the quantized motion of solitons to the Chern number of the underlying Bloch band. More generally, these developments introduce a theoretical framework by which a broad class of nonlinear topological phenomena can be formulated in terms of topological band indices. We then broaden the scope by applying this nonlinear topological framework to the realm of quantum gases. By considering an instructive mapping to a Bose-Bose atomic mixture on a lattice [40][41][42] , which encompasses the aforementioned DNLS as its semiclassical limit, we identify a scenario by which a topological pump emerges from inter-particle interaction processes: a soliton of impurity atoms is dragged by the driven majority atoms, hence leading to interaction-induced topological transport.

Topological pumps of solitons: General theory
The theoretical description of soliton pumping relies on the topological character of Wannier functions, namely, the displacement of Wannier centers per pump cycle by the Chern number of the associated Bloch band, see Fig. 1a. An intuitive reason for soliton pumping then is solitons are dragged by Wannier functions, resulting in their topologically quantized motion, see Fig. 1b. Interestingly, nonlinearities can induce this drag effect in topologically trivial systems. To demonstrate this effect, we consider a 1D atomic Bose-Bose mixture in a species-selective optical lattice, where the majority (minority) atoms experience a topological (trivial) lattice. We show that a soliton of minority atoms undergoes quantized displacement by activating a Thouless pump sequence for the majority atoms; see Fig. 1c. In the following, we elaborate on the theoretical description of soliton motion in nonlinear Thouless pumps and the interaction-induced topological pumps for atomic mixtures.
Our theoretical framework concerns a generic class of lattice models governed by the DNLS, where the field ϕ i,α is defined at the lattice site α of the ith unit cell; H(t) is a time-dependent Hamiltonian matrix, which includes a Thouless pump sequence 7,38 ; and g > 0 is the (onsite) nonlinearity strength. Equation (1) preserves the norm of the field, which we set to ∑ α,i |ϕ i,α | 2 = 1, without loss of generality. An illustrative model, used below to validate the general theory, is provided by the two-band Rice-Mele model 15 : a 1D chain with alternating couplings J 1,2 (t) and staggered potential ±Δ(t) (Methods). Considering the nonlinear Rice-Mele model, Eq. (1) takes the more explicit form Here, the Thouless pump cycle corresponds to a loop in the parameter space spanned by (J 2 − J 1 ) and Δ, which encircles the origin (J 1 = J 2 , Δ = 0); see Methods. When g = 0, the Bloch bands defined in momentum-time space are associated with a Chern number C = ±1. This topological invariant is known to determine the quantized displacement for a filled band upon each cycle of the pump 38 .
Our analysis starts by studying the adiabatic evolution associated to the general Eq. (1), which is characterized by the period of the pump T (exceeding all other time scales). To simplify notations, we use the multi-index i = (i, α) and write H ij H αβ ij ðtÞ. Introducing the adiabatic time s = t/T, Eq. (1) takes the form iε∂ s ϕ i = ∑ j H ij (s) ϕ j − g|ϕ i | 2 ϕ i , where ε = 1/T. The solutions to the adiabatic DNLS can be well approximated by stationary states of the form ϕ i / e Àiθ s φ i , where θ s is a time-dependent phase factor and φ i is an instantaneous solution to the stationary nonlinear Schrödinger equation (see Methods and refs. 43,44) where the instantaneous eigenvalue μ s explicitly depends on the adiabatic time s. Equation (3) admits (bright) solitons as stationary state solutions, which are stable localized structures in the bulk. For sufficiently weak nonlinearity, solitons predominantly occupy the band from which they bifurcate 45 , while increasing nonlinearity leads to band mixing. In real space, solitons are immobile without external forcing, and are degenerate modulo a lattice translation set by the translational  In a nonlinear setting, the motion of a soliton follows the quantized Wannier drift. c In a Bose-Bose atomic mixture, quantized pumping can be induced by interactions: a soliton of impurity atoms is dragged by the driven majority atoms, leading to interaction-induced topological transport. symmetry of the system. By adiabatically changing the Hamiltonian H ij (s), a single soliton undergoes smooth deformation, and after one period, it is mapped to the manifold of initial solutions, implying translation by an integer multiple of the unit cell. The observations of ref. 34 suggest that solitons bifurcating from a single Bloch band undergo a quantized displacement dictated by the Chern number of the band 7 over each pump cycle. Demonstrating this intriguing relation between the transport of nonlinear excitations and topological band indices is at the core of the present work.
To elucidate the topological nature of nonlinear pumps, we follow ref. 37 and represent the solitons of Eq. (3) in the basis of maximally localized Wannier states, where ω ðnÞ l = 1=N P NÀ1 k = 0 expði ð2π=NÞ k lÞ ϵ ðnÞ k is the Fourier transform of the nth Bloch band ϵ ðnÞ k associated with H ij (s); N is the number of unit cells; n = ðn, n 1 , n 2 , n 3 Þ, l = ðl, l 1 , l 2 , l 3 Þ; and W ðnÞ ðlÞ are the following Wannier overlaps The Wannier states of a Bloch band are not unique, as they depend on the gauge choice for the Bloch functions 46 . Nevertheless, a unique set of maximally localized Wannier functions is provided by the eigenstates of the position operator's projection onto the associated band. Since such Wannier functions are exponentially localized, the contribution to the Wannier overlaps in Eq. (6) from Wannier functions corresponding to different unit cells are negligible. The Wannier overlaps can thus be simplified as W ðnÞ l = W ðnÞ δ ll 1 δ l 1 l 2 δ l 2 l 3 , where W ðnÞ = P j w ðnÞ* j ðlÞ w ðlÞ; we point out that W ðnÞ does not depend on the index l, because of translational invariance.
Moreover, in the regime of weak nonlinearity, the initial state soliton occupies a single band 34,37,47 , which allows us to neglect interband contributions to Eq. (5). We note that this simplification holds throughout the evolution of the pump, during which the soliton adiabatically follows the same band.
Under those realistic assumptions, the Wannier representation of the DNLS reduces to the form where W ðnÞ = P j |w ðnÞ j ðlÞ| 4 . Equation (7) has the form of a scalar DNLS on a simple lattice with one degree of freedom per unit cell labeled by Wannier indices l, with hopping terms involving nearest and beyondnearest neighbors. The properties of such scalar DNLS are well established 36,48,49 : Equation (7) admits inter-site solitons, with maxima on two adjacent sites, and on-site solitons, with their maximum on a single site. The inter-site solitons are known to be unstable against small perturbations, we thus restrict ourselves to the stable on-site solitons. Crucially, on-site solitons are always peaked around a single site (l) throughout the pumping cycle, as there is a finite energy (Peierls-Naborro) barrier for delocalization 36,49 . Interestingly, the Peierls-Naborro barrier plays a role analogous to the "gap condition" of conventional topological physics, by forbidding transitions to other stable states during the adiabatic time evolution (Methods). This observation suggests that solitons are dragged by Wannier states upon pumping, hence exhibiting a quantized displacement in real space established by the Chern number 15,38,39 ; see Fig. 1a, b.
To firmly prove the topological nature of the nonlinear Thouless pump, we evaluate the solitons center-of-mass displacement after one period s = 1 (Methods) Δhφ ðnÞ , X φ ðnÞ i = Δhw ðnÞ ð0Þ, Xw ðnÞ ð0Þi where X is the position operator of the lattice; h f , gi P i f * i g i is the inner product of fields on the lattice; and Δ(⋅) ≡ (⋅) s=1 − (⋅) s=0 . The first term in Eq. (8) reflects the displacement of Wannier functions upon one pump cycle, which is known to correspond to the Chern number of the band 15,38,39 ; the additional terms displayed on the second line are possible corrections due to the finite overlap of different Wannier states. Importantly, we find that these small interference effects are periodic in time (Methods), such that these correction terms in Eq. (8) do not contribute to the solitons center-of-mass displacement over a pump cycle. Altogether, this completes the proof: the displacement of solitons is indeed quantized according to the Chern number of the band from which they emanate.

Numerical validation
We now demonstrate the validity of our assumptions by solving the nonlinear Rice-Mele model [Eq. (2)]. In Fig. 2a, b, we compare the on-site soliton solution of the simplified Eq. (7), which emerges from the lowest band, with the Wannier representation of the exact soliton obtained by We depict the motion of the exact soliton in Fig. 3, as obtained by solving Eq. (3) over two pump cycles s ∈ [0, 2], and we compare this trajectory with the drift of its underlying Wannier function, i.e., the Wannier state that contributes the most to the expansion (4). In order to obtain a contiguous path for the Wannier center, we relabeled the Wannier functions whenever the Wannier centers met discontinuities; this smoothing corresponds to a singular gauge transformation of the corresponding Bloch states, and has no physical implication. Figure 3 indicates that the trajectories of the soliton and Wannier center differ at intermediate times (s ≠ integer), which we attribute to the aforementioned interference effects involving different Wannier states (Methods); however, in agreement with our theoretical predictions, this deviation remains small and time-periodic over the whole pump cycle, and does not introduce any (integer) correction to the quantized center-of-mass displacement.

An interaction-induced topological pump for ultracold atomic mixtures
The theoretical framework presented in this work is based on the general DNLS in Eq. (1), and hence, it applies to a broad range of nonlinear lattice systems. In particular, this equation corresponds to the Gross-Pitaevskii equation describing a weakly-interacting Bose gas evolving on a moving lattice potential. In this section, we propose to go beyond the paradigm of nonlinear pumps for single-component bosonic systems, by introducing a mapping to an imbalanced Bose-Bose atomic mixture, which encompasses the DNLS in Eq. (1) as its semiclassical limit (within the Thomas-Fermi approximation). As we explain below, this original approach reveals an interaction-induced topological pump, where solitons of impurity atoms undergo a quantized drift resulting from genuine interaction processes with their environment.
We start from a microscopic theory for an imbalanced Bose-Bose atomic mixture on a 1D lattice 42 , as described by the second-quantized whereφ i andσ i are bosonic field operators on the lattice; note that we use the same conventions for indices i = (i, α) as before. Specifically, the first line describes single-body processes (i.e., nearestneighbor hopping and onsite potentials) and intra-species contact interaction processes for the majority atoms, which are described by the field operatorφ i ; the second line describes single-body processes and intra-species contact interactions for impurity atoms, represented by the field operatorσ i ; and the third line describes inter-species interaction processes. We assume that the intra-species interaction strengths are both repulsive, (U σσ , U ϕϕ > 0), whereas the inter-species interaction strength is attractive (U ϕσ < 0). In the semi-classical limit, where quantum fluctuations are suppressed for both species, this Bose-Bose mixture setting is well described by two coupled nonlinear Schrödinger equations (Methods and ref. 42) where ϕ i and σ i denote classical fields satisfying the constraints , with N ϕ and N σ the particle number of majority and impurity species, respectively; the interaction parameters are defined as g αβ = U αβ (N ϕ + N σ ), with α, β = (ϕ, σ); μ ϕ,σ denote the chemical potentials; and ω 0 is the eigenvalue of the nonlinear Eqs. (10).
Considering the case of heavy impurities, we neglect their kineticenergy contributions (H ðσÞ ij ) to Eq. (10), the so-called Thomas-Fermi approximation. In this regime, one can relate the impurity mean-field profile to the majority profile as and Eq. (10) simplifies to the DNLS (Methods) where g = Àg ϕϕ + g 2 ϕσ =g σσ . Interestingly, Eq. (12) is formally equivalent to the DNLS in Eq. (3): the majority atoms described by the field ϕ i can form a soliton and undergo a quantized motion upon driving a Thouless pump sequence in the corresponding lattice Hamiltonian, i.e., H ðϕÞ ij ðsÞ. Importantly, according to Eq. (11), the impurity atoms also form a soliton and undergo a quantized motion: the impurities exhibit topological pumping from genuine interaction processes with the majority atoms. In particular, this interaction-induced topological pumping occurs even when the lattice felt by the impurities H ðσÞ ij is associated with a trivial (non-topological) band structure. This intriguing phenomenon, which could be implemented in ultracold atomic mixtures in optical lattices [40][41][42] , is reminiscent of topological polarons [50][51][52][53][54][55] , in the sense that impurities inherit the topological properties of their environment through genuine interaction processes. We first analyze this interaction-induced topological effect by considering the Thomas-Fermi approximation. It appears from Eq. (12) that u MF acts as an effective potential for the majority atoms; a soliton then emerges as the bound state of the impurity field. In the context of highly-imbalanced mixtures with strong impurity-majority coupling, i.e., in the strong-coupling Bose polaron regime, it is customary to assume a variational ansatz describing the profile of the impurity and majority fields 56 ; the majority field is then found as the bound state of the impurity potential u MF using the first relation in Eq. (12). Here, the variational problem for obtaining u MF reduces to one for ϕ, because of the constraint u MF = g|ϕ| 2 . As before, we express ϕ in the Wannier basis, and the variational problem is then solved simultaneously for both u MF and ϕ using the ansatz a l = η sechðξ ðl À l 0 ÞÞ for the Wannier coefficients of ϕ. The bound state of the resulting impurity potential u MF = g|ϕ| 2 then corresponds to the soliton (Methods). Figure 4 a, b show the adiabatic evolution of the amplitude η and width ξ of the variational solution a l = η sechðξ ðl À l 0 ÞÞ used for the Wannier coefficients of ϕ. We compare these results with the amplitude and width extracted from the bound-state solution associated with the impurity potential u MF = g|ϕ| 2 , as well as to those extracted from the exact soliton of Eq. (3) expressed in Wannier representation. We also show the dependence of these parameters on the nonlinearity g in Fig. 4c, for both the exact soliton and the variational solution. These results validate our variational approach, as well as the boundstate picture of our soliton.
The minimum-energy solutions obtained from the variational ansatz are realized for integer values of the Wannier index l 0 , and thus correspond to stable on-site solitons. Moreover, this Wannier index l 0 remains constant over a pump cycle. Hence, this again suggests that the real-space motion of the soliton should follow the quantized Wannier drift, as established by the Chern number. This is verified in Fig. 4d, where the center-of-mass displacement of the calculated bound state is shown to be quantized over a pump cycle (compare with Fig. 3).
In order to demonstrate the validity of our results, in particular, the robustness of the interaction-induced topological pump away from the Thomas-Fermi limit, we solve Eq. (10) numerically for a massbalanced mixture, thus including the effects of the impurities' kinetic energy. We again use the Rice-Mele model, but consider two different pump sequences for the majority and impurity species: the majority feels the same (topological) pump sequence as in Fig. 3, while we apply a trivial sequence for the impurity species. We obtain the steady state solution of Eq. (10) over two pump cycles, where the majority particles predominantly occupy the lowest Bloch band. The corresponding Center-of-mass (CM) trajectories of both species in a Bose-Bose mixture. Here, the interaction strengths are set as g ϕϕ ≃ 0.226 J 0 , g ϕσ ≃ − 11.32 J 0 and g σσ ≃ 2.26 J 0 . The majority atoms undergo the pumping sequence as in Fig. 3, while impurities feel a trivial sequence (Methods). Impurity atoms undergo quantized transport through interactions with their environment.
trajectories of the CM of both species are depicted in Fig. 5, where the impurity CM is shown to be dragged by the majority particles. While the exact form of the CM trajectories depend on the details of the model and pumping sequence, the CM displacement after one pump cycle is dictated by the Chern number of the topological band occupied by the majority species (C = −1 in this case). Although the impurity atoms experience a topologically trivial lattice, they are shown to undergo topological pumping through genuine interaction effects with their environment.

Implementation in ultracold atoms
The interaction-induced topological pump introduced above could be experimentally implemented in ultracold atomic gases involving two bosonic species. In fact, the parameters values incorporated in our numerical simulations of Eq. (10), and displayed in Fig. 5, are compatible with an experimental realization based on bosonic 7 Li − 7 Li mixtures, with two different hyperfine states of 7 Li as "majority" and "impurity" atoms; we note that the formation of solitons in Lithium gases was previously investigated, both theoretically and experimentally 57,58 . Following ref. 59, the scattering lengths between atoms in state (F = 1, m F = 1)-"impurity" atoms-and (F = 1, m F = 0)-" majority" atoms-can be set to a ϕϕ ≃ 0.154 a 0 , a ϕσ ≃ −7.57 a 0 , a σσ ≃ 1.514 a 0 , at a magnetic field B ≃ 575 G, where a 0 is the Bohr radius (a 0 = 0.0529 nm); we note that these scattering lengths are highly tunable thanks to a broad Feshbach resonance. As further discussed below, this configuration is compatible with the interaction parameters (g ϕϕ , g σσ , g ϕσ ) used in our numerics. The lattice structure and pump sequence can be designed within a time-dependent optical lattice. For instance, following ref. 16, the atoms can be loaded in a potential landscape comprised of two superimposed optical lattices, with a long-wavelength lattice (λ l = 1064 nm) and a shorter lattice (λ s = λ l /2), with different amplitudes (V l = 3.0 E R and V s = 1.0 E R , with E R = h 2 =ð2m λ 2 l Þ the recoil energy of the long lattice). Such an optical lattice potential takes the form V ðx, ϕÞ = À V l cos 2 ð2πx=λ l À ϕÞ À V s cos 2 ð2πx=λ s Þ, and it implements the Rice-Mele lattice considered in our numerics: the Thouless pump sequence is simply realized by sweeping the phase ϕ from 0 to 2π. The relevant parameters of the Rice-Mele model can be extracted from a tight-binding analysis of the optical lattice potential 16 , and the resulting pump sequence is described by the following elliptic path in parameter space: ððJ 1 À J 2 Þ=aÞ 2 + ðΔ=bÞ 2 = 1, with a ≃ 0.19 E R and b ≃ 0.475 E R . In our numerics, we choose a closely related pumping sequence with a = 0.15 E R and b = 0.5 E R ; this choice does not affect our final conclusions, since topological pumping is robust against smooth deformations of the pumping sequence. Finally, to reveal the interactioninduced topological transport for impurities, we propose to implement a trivial pump sequence for that species only [see Fig. 5]; this could be realized by designing a state-dependent optical lattice 60 , for instance, using the Floquet-engineering scheme of ref. 61. The particle numbers of the two species can be set to N ϕ ≃ 1500 62 and N σ /N ϕ ≃ 1/30. With this choice, we obtain the interaction parameters according to the relation g αβ =E R = ðN σ + N ϕ Þ ffiffiffiffiffiffiffiffi ffi 8=π p k l a αβ ðV s =E R Þ 3=440 , where α, β = (ϕ, σ) and k l = 2π/λ l . Setting the pump parameter J 0 = 0.5E R , the numerical values for the interaction parameters are obtained as g ϕϕ ≃ 0.226 J 0 , g ϕσ ≃ −11.32 J 0 and g σσ ≃ 2.26 J 0 , which are the values used in our numerical simulations [ Fig. 5].

Discussion
In this work, we outlined a general theoretical framework that connects Bloch band's topology to nonlinear excitations, hence elucidating the topological transport of solitons in the context of nonlinear Thouless pumps. Solitons are stable states of nonlinear lattice systems described by the paradigmatic discrete nonlinear Schrödinger equation (DNLS), which is central in describing nonlinear phenomena in a wide range of physical settings, from nonlinear optics and photonics, to ultracold quantum matter, fluid dynamics and plasma physics. In this sense, characterizing the influence of Bloch band's topology on the behavior of the stable states of DNLS is of prime importance. This program is particularly challenging due to the lack of generic theoretical approaches connecting notions of topological physics to nonlinear systems and vice versa. Furthermore, introducing nonlinearities in more sophisticated topological systems, such as higher-dimensional settings, or lattices exhibiting higher-order topology and symmetry-protected features, could lead to exotic phenomena exhibited by the nonlinear modes of the system; see ref. 63 and references therein. By providing a scheme that naturally connects topological indices of band structures to nonlinear excitations, our work opens the door to the exploration of novel nonlinear topological phenomena.
We also illustrated the universality of our approach, by introducing a topological pump for Bose-Bose atomic mixtures, where one species (impurity atoms) experience a quantized drift through genuine interaction processes with the other species (the surrounding majority atoms). Importantly, the impurity atoms inherit the topological properties of their environment through inter-species interactions. We note that such interaction-induced topology has been previously studied in the context of topological polarons, namely, in mixtures with strong population imbalance, where individual topological excitations can bind to mobile impurities 50,51,53,55 . The present scheme extends those concepts to more complex majority-impurity states, such as coupled coherent states within a superfluid phase. We also point out that the proposed scheme can be implemented using available coldatom technologies, and the quantized transport of impurities can be measured in-situ, using state-selective imaging techniques 64 . Besides, the Chern number characterizing the interaction-induced topological pump could also be directly extracted by interferometry 51 .
During the preparation of this manuscript, the authors became aware of a related work by M. Jürgensen and M. C. Rechtsman 47 , and also ref. 65.

Adiabatic theorem for NLS
The adiabatic theorem for NLS (both continuous and discrete forms), follows closely the formulation of its linear counterpart 43,44 . For a system with a time-dependent Hamiltonian H(t), which varies on a time scale T much larger than all the time scales in the problem, the timedependent NLS takes the following form (see main text) where s = t/T is the adiabatic time and ε = 1/T the rate of change. The stationary state solutions of Eq. (13) are of the form where φ s is the instantaneous solution of the stationary NLS, and θ s = 1=ε R s 0 ds 0 μ s 0 À γ s À Á is a global phase factor consisting of a dynamical contribution and a Berry phase, and it can be ignored. The correction term δ φ s accounts for non-adiabatic variations, and for ε → 0, it behaves as ||δ φ||~ε, hence vanishes in the adiabatic limit ε → 0. The relevant dynamical information is therefore encoded in the instantaneous solutions of Eq. (15).

The Rice-Mele model and pump sequence
Throughout this work, we illustrate the general concepts and results using the Rice-Mele model, with periodic boundary conditions. This simple two-band model, which is reviewed in some detail below, is known to exhibit a topological (Thouless) pump sequence.
The Rice-Mele model is a 1D tight-binding model with alternating nearest-neighbor tunneling matrix elements (J 1 , J 2 , J 1 , J 2 ,…), and a staggered on-site potential. We denote the two sites within each unit cell by α = A, B and the unit cells by i, 0 ≤ i ≤ N − 1, where N is the number of unit cells. The hopping matrix element between sites A and B within each unit cell (resp. between adjacent unit cells) is written as J 1 = − J(1 + δ) (resp. J 2 = −J(1 − δ)) and the magnitude of the staggered potential on site A (resp. B) equals Δ (resp. −Δ). The Hamiltonian of the Rice-Mele model thus reads The simulations shown in the main text were performed on a lattice with N = 100 unit cells, and using the following pump sequence with J 0 = 0.5 and δ 0 = 0.6, corresponding to a topological pump with Chern number C = −1. The nonlinear Rice-Mele model, which is used in our simulations, is obtained by adding an on-site nonlinearity to this lattice model; see Eq. (2).
In order to demonstrate the interaction-induced topological pumping in the Bose-Bose mixture setting, we assume that the two species experience the same Rice-Mele lattice described above, but with different pump sequences: the majority atoms experience the topological pumping sequence in Eq. (18), while the impurity atoms experience a trivial sequence with Jδ = Δ = 0. The resulting center-ofmass displacement of both species are depicted in Fig. 5 of the main text.

Derivation of the scalar DNLS
We outline the derivation of the simplified scalar DNLS from the original lattice DNLS, The Wannier functions are related to the Bloch waves of the Hamiltonian by the following relations: where ψ ðnÞ j ðkÞ = e ið2π=NÞkðjÞ u ðnÞ j ðkÞ is the Bloch wave of band n with momentum k and u ðnÞ j ðkÞ is the corresponding Bloch function, which is periodic over the unit cells and does not depend on j. To represent the Hamiltonian part in Wannier basis, we evaluate the matrix elements of the Hamiltonian over the Wannier states where ω ðnÞ l = 1=N P NÀ1 k = 0 e ið2π=NÞkðlÞ ϵ ðnÞ k is the Fourier transform of the Bloch band ϵ ðnÞ k ; see main text. Next, we express the nonlinearity in terms of Wannier functions, hw ðnÞ ðlÞ, |ϕ| 2 ϕi = X n 1 ,n 2 ,n 3 X l 1 ,l 2 ,l 3 Taking the inner product of Eq. (19) with w ðnÞ l and using Eqs. (21) and (22), we obtain the following DNLS where we used Eq. (24) in the last equality. The first term in the last equality of Eq. (25) reduces to hw ðnÞ ð0Þ,Xw ðnÞ ð0Þi s since we normalized the soliton intensity to unity, N ϕ = P l |a ðnÞ l | 2 = 1. The second term in the last expression is the mean value of the position of the Wannier functions indices, which is constant since the on-site solution is always peaked around a Wannier label and remains symmetric around it. Its contribution to the displacement over a pump cycle thus vanishes. The third term contains products of the form P l a ðnÞ* l + δl a ðnÞ l hw ðnÞ ðδlÞ, Xw ðnÞ ð0Þi s and its treatment requires more care. The coefficient P l a ðnÞ* l + δl a ðnÞ l is time-periodic, since a ðnÞ l is, by assumption, the solution of the scalar DNLS in Eq. (7), in the main text. To investigate the behavior of hw ðnÞ ðδlÞ, Xw ðnÞ ð0Þi s , we note that after a pump cycle, the Wannier functions are displaced by the Chern number, w ðnÞ ðlÞ| s = 1 = w ðnÞ ðl + C n Þ| s = 0 , with C n the Chern number of where we used Eq. (24) in the last step. This proves that the quantity 〈w (n) (δl), Xw (n) (0)〉| s , in the last equality of Eq. (25), is a time-periodic quantity.
Altogether, the third term in Eq. (25) is also time-periodic, and the soliton's center-of-mass displacement over a pump cycle is given by Δhφ ðnÞ , X φ ðnÞ i = Δ hw ðnÞ ð0Þ, Xw ðnÞ ð0Þi: ð27Þ This result directly relates the soliton's displacement to the displacement of Wannier functions upon one pump cycle, as dictated by the Chern number of the band 15,38,39 . This proves the quantized pumping of the soliton according to the Chern number.

Derivation of the Bose-Bose mixture equations
In order to derive the equations governing the coherent state profiles of the two species in the mixture, we start from the microscopic Hamiltonian in Eq. (9). The coherent-state action of the system takes the following form (ℏ = 1), with the Lagrangian To proceed, we seek stationary state solutions for the coherent state fields of the form ϕ ðssÞ i ðtÞ = e Àiω 0 t ϕ i and σ ðssÞ i ðtÞ = e Àiω 0 t σ i , which minimize L½ ϕ, ϕ; σ, σ. Such solutions are the saddle-point solutions of the quantum mechanical action, giving the mean-field stable states of the system. The Lagrangian then takes the time-independent form To minimize the Lagrangian, the corresponding Euler-Lagrange equations are derived from δL=δ ϕ i = 0 and δL=δ σ i = 0 , which leads to the two coupled equations in Eq. (10) in the main text.
Variational ansatz for the state of Bose-Bose mixture in the Thomas-Fermi limit The variational treatment of Eqs. (11) and (12) accounts to minimizing the following energy functional for the field ϕ From the knowledge obtained from the soliton solutions of the DNLS in the main text, we assume that ϕ i belongs to a single band and expand it in terms of the Wannier functions of the corresponding band, ϕ i = P l a ðnÞ l w ðnÞ ðlÞ. We then use a sech variational ansatz for the coefficient amplitudes, a ðnÞ l = η= sech ξðl À l 0 Þ À Á . The variational energy functional takes the following form: For the simulations presented in the main text [ Fig. 4], we assume that N ϕ = 1 ; see refs. 36, 48 for more details on variational ansätze for DNLS. From the solution of Eqs. (35) and (36) we then obtain the boson field, ϕ i , which is then used to obtain the effective attractive potential u MF i = g|ϕ i | 2 ; see Eq. (12).

Potential barrier preventing soliton delocalization
The nonlinear term in Eq. (1) [resp. in Eq. (7)] leads to the formation of localized soliton solutions, which do not satisfy the lattice (resp. Wannier lattice) translational symmetry. While the stable-state solitons are not translationally invariant, they can be mapped to one another through lattice translations. It is known that an effective potential barrier exists for continuous deformations of each stable-state soliton to a neighboring one. This potential barrier is reminiscent of the Peierls-Naborro barrier (PNB) known in the theory of dislocation dynamics in crystals 66 . Under adiabatic evolution, a soliton in Wannier representation will always remain peaked on a single Wannier index since the potential barrier rules out the existence of solutions that interpolate continuously between two on-site solitons. The strength of this potential barrier can be estimated in terms of the model parameters using the variational ansatz a ðnÞ l = η=sech ξðl À l 0 Þ À Á for the Wannier soliton in band n, and the corresponding energy functional in Eq. (35), The estimated Δ Barrier depends on the model parameters via the Fourier transform of the dispersion relation at l = 0, ω ðnÞ 0 = 1=N P NÀ1 k = 0 ϵ ðnÞ k , and the interaction parameter g. We verified that the expression in Eq. (37) is in agreement with a result found in ref. 36 for DNLS equations with nearest-neighbor hopping.

Data availability
The data that support the findings of this study are available from the corresponding author upon reasonable request.