Spatial inversion symmetry breaking of vortex current in biased-ladder superfluid

We investigate the quench dynamics of interacting bosons on a two-leg ladder in presence of a uniform Abelian gauge field. The model hosts a variety of emergent quantum phases, and we focus on the superfluid biased-ladder phase breaking the $Z_{2}$ symmetry of two legs. We observe an asymmetric spreading of vortex current and particle density, i.e., the current behaves particle-like on the right and wave-like on the left, indicating spontaneous breaking of the spatial inversion symmetry. By decreasing the repulsion strength, it is found the particle-like current is more robust than the wave-like one. The evolution of entanglement entropy manifests logarithmic growth with time suggesting many-body localization matters.

We investigate the quench dynamics of interacting bosons on a two-leg ladder in presence of a uniform Abelian gauge field.The model hosts a variety of emergent quantum phases, and we focus on the superfluid biased-ladder phase breaking the Z2 symmetry of two legs.We observe an asymmetric spreading of vortex current and particle density, i.e., the current behaves particlelike on the right and wave-like on the left, indicating spontaneous breaking of the spatial inversion symmetry.By decreasing the repulsion strength, it is found the particle-like current is more robust than the wave-like one.The evolution of entanglement entropy manifests logarithmic growth with time suggesting many-body localization matters.
As shown in a landmark study based on bosonization in 2001 60 , the Meissner and vortex phases, which are reminiscent of a type-II superconductor, were found in the two-leg flux-ladder model.A variety of emergent quantum phases are hosted by the bosonic flux ladders in presence of on-site interaction 52,[60][61][62][63][64][65][66][67][68][69][70][71][72][73][74][75][76] , including Meiss- * Electronic address: yaoyao2016@scut.edu.cnner phase, vortex-liquid phase, vortex lattice phase, the charge-density wave (CDW) phase, and the biased-ladder phase (BLP).The local current configurations serve as one of the most important quantities to distinguish the phases.The Meissner phase possesses a finite uniform chiral leg current encircling the ladder while the rung current vanishes.Currents in CDW and BLP are very similar with that in the Meissner phase but would break some discrete symmetries 66 .In the vortex phases, the rung currents on the inner rungs develop, forming vortices in the system depressing the chiral current.Except CDW, other phases can be further divided into two phases, superfluid and Mott insulator, which can be characterized by calculating the central charge, as well as the entropy 58,66 .More interestingly, the main characteristics of BLP turns out to be the finite leg-population imbalance, stabilized by the inter-chain interactions 58 .In the thermodynamic limit the ground state would thus be twofold degenerate, and the subsequent Z 2 symmetry in terms of inversion of two legs and the signs of relevant flux is spontaneously broken 66 .Considering these remarkable features lead to exotic dynamical effects, therefore, we focus on the superfluid BLP phase in the present work.
As sketched in Fig. 1, the Hamiltonian of a paradigmatic two-leg flux-ladder model is given by which is residing on a ladder with L rungs.Herein, the local operator a † ℓ,r (a ℓ,r ) creates (annihilates) a boson on the lower (ℓ = 1) or the upper (ℓ = 2) leg of the r-th rung; J and J ⊥ are nearest-neighbor hopping constants along legs and rungs, respectively; φ is the flux per ladder plaquette; the bosons also have on-site repulsion U and inter-chain repulsion V .Throughout this work, we set φ = 0.85π, J = 2, J ⊥ = 3, U = 2, ρ = 0.8, the lattice constant a = 1 and = 1.By this setting, the system has been determined to reside in the superfluid BLP phase 66 , which is essential to produce the symmetry breaking as discussed below, and we always choose the ground state with lower leg population larger than upper.
It is worth noting that, the model possesses a gauge freedom that one can choose different Peierls phase factors as long as the total flux of a single ladder plaquette remains invariant.The chosen gauge described by Hamiltonian (1) is so-called rung gauge in which the hopping matrix elements on the legs are real and on the rungs are complex.
By the Heisenberg equation of motion r ], we define a local current operator on legs as and the current operator on the rung as It is also convenient to define a background chiral current to characterize the average current circulating the ladder along the legs.The ground state results are calculated by using the density matrix renormalization group (DMRG) method [77][78][79][80] .We simulate the flux ladder up to L = 200 rungs and the bond dimension is typically up to 1000.
As the existence of repulsive interactions, in the numerical calculations it is safe to restrict the local basis to small values.We keep at most four bosons per site and have also checked that with six bosons per site which produces consistent results.The subsequent time evolution is then simulated via time evolving block decimation (TEBD) method [81][82][83] .We typically employ bond dimension of 1000 to ensure convergence of local observables and the bond dimension up to 5000 is used in the calculation of entropy.
Let us first describe the main scenario we are investigating as displayed in Fig. 1.The BLP ground state initially hosts a finite uniform anticlockwise chiral leg current encircling the ladder and vanishing rung current.A single rung in the middle is then excited to generate vortex current at time t = 0 by applying a rung current operator without phase factor to the superfluid BLP ground state, i.e., i a † 1,L/2 a 2,L/2 − a † 2,L/2 a 1,L/2 |ψ 0 .The excitation creates accumulation of populations on the upper leg (ℓ = 2) which we call it 'particle' (bigger dots) and three sites with exhausting populations on the lower leg (ℓ = 1) as we call 'holes' (circles).This leads to current change on top of background chiral currents which generates two adjacent clockwise vortices (red squares) with current (black arrow) on the edge.On the right side, the 'particle' moves along the same direction with the upper leg current of clockwise vortex making them constructive and move outwards like a solitary wave.On the left side, however, the concentration of 'hole' and 'particle' oscillate between two neighboring sites on the same rung, which is similar to the recombination process of majority and minority carriers in semiconductors.The corresponding rung current oscillates between up and down legs in the same fashion, giving rise to a vortex on the initial position alternating between clockwise and counterclockwise.In addition, vertices with opposite directions are destructive, so they move forward following with quick decay.This resembles how plane waves spread.In short, vortex currents behave particle-like on the right and wave-like on the left in just one ladder.
Snapshots of current configurations in the time evolution after rung current excitation are shown in Fig. 2. The vortex currents generated by the initial excitation on the upper leg j 2,r and lower leg j 1,r can be regarded as 'particle' and 'hole' current, respectively.The initial single wave packets of both 'particle' and 'hole' then broaden, which as expected by particle-hole symmetry are shown to be almost identical in magnitude but the signs are different.This means 'particles' go left and 'holes' go right, implying a net charge current appears.Hence, in order to manifest the symmetry breaking of spatial inversion, we define a 'charge' current j c r = j 1,r − j 2,r to fingerprint the vertices.It is observed that, the central 'charge' current behaves like a source, emitting vertices outward and moving.Positive and negative vortices alternate and oscillate on the left side, while a big solitary wavepacket visibly moves on the right side: The asymmetry appears.
We relevantly calculate the population configurations on each leg n ℓ,r (not shown) which symmetrically spread out from the middle site.To see the asymmetric influence of chiral current, we define the relevant 'charge' concentration as n c r = n 1,r + n 2,r .As shown in Fig. 3, after the wave packet splits into two counter-propagating parts, the main part of wavepacket on the left is blocked by the chiral current and the minor spreading part decays quickly.On the other hand, the decay of the right wavepacket is much slower.More importantly, regardless of the slight broadening, a distinguishable solitary wavepacket shape is always reserved.This visible asymmetry significantly exhibits that the spreading of 'charge' is inversion symmetry broken in space.Although not shown, we can also parallel define the 'spin' concentration whose spreading will be different from that of 'charge' concentration and a similar effect of spin-charge separation could be observed.
To figure out the origin of this asymmetric spreading, the parameter of inter-chain repulsion is changed.Fig. 4 shows the net current j ν = j c r − j ch of the two counterpropagating wave packets and the maximum point of 'charge' concentration n c r .From Fig. 4(a)-(c), one can clearly see that following the increase of V , the wave on the left is fading while the right particle-like current is fairly stable.That is, with small V , the net current of the left wave is larger than the right one, and with large V , the left current fades and the spreading then prefers the right side.The asymmetric fading facilitates the imbalance of the two sides which mainly leads to the asymmetric spreading.In addition, the speeds on the two sides are roughly equal and linear with time.
As described above, the excitation creates the 'particle' and 'holes' which spread from the central rung creating a counterclockwise vortex (blue in Fig. 1) moving left and a clockwise vortex (red in Fig. 1) moving right.It is essential that the counterclockwise vortex enlarges the chiral current j ch while the clockwise vertex suppresses it.For cases with small V , the spreading of 'particle' and 'hole' will weaken the imbalance of leg-population so that the system is almost close to the Meissner phase.Since the chiral current grows monotonically from BLP to Meissner phase 58,66 , the counterclockwise vortex (blue in Fig. 1) will be mutually enhanced with the chiral current.Following the increase of V , the chiral current is suppressed and the imbalance of leg-population is then stabilized.As a consequence, the counterclockwise vortex is decreased, and when V exceeds a certain value, the 'particles' and 'holes' on the left side are almost blocked and difficult to spread out.In this situation, the concentrations of 'hole' and 'particle', together with the corresponding rung current, start to oscillate on the same rung.On the other hand, considering that the clockwise vortex always decreases the chiral current, V does not matter on it.Therefore, the clockwise vortex moving right overtakes, suggesting the particle-like vortex current is more robust in a counterintuitive manner.
In a free space, the nature of microscopic particle is normally a plane wave due to the thermalization hypothesis [84][85][86] .To comprehend the robustness of particle-like current, we further calculate time evolution of the von Neumann entanglement entropy, defined as S vN = −Tr[ρln(ρ)] with ρ being the reduced density matrix by cutting both two legs at the same points to divide the ladder into two halves, as shown in Fig. 5.As the needed bond dimensions in such simulation grow exponentially with time, a smaller system size L = 80 and a cutoff value 10 −6 in the truncation are used with bond dimension up to 5000 states.One can find that as well fitted by a logarithmic function, the entropy grows logarithmically with time which is a typical signature of many-body localization [87][88][89][90][91] .This then interprets why the particle-like shape can be stable.
In summary, we have studied the time evolution of vor- tex current and concentration in a superfluid BLP phase.An asymmetric spreading of the vortex current is found determined by the direction of background chiral current.We find the right-moving particle-like current is more robust than the left-moving wave-like one, which stems from the many-body localization as indicated by the logarithm of entanglement entropy.The particle-like current therefore has got great potential in applications, such as quantum computations.

FIG. 1 :
FIG. 1: Sketch of the rung current excitation creating asymmetric spreading process.Typical states on the ladder at four instant points are drawn from top to bottom.At t = 0, the rung current excitation operator (green arrow) is applied to the BLP ground state.The local population is denoted by the size of dots and the exhausted population is denoted by empty circles.The current deviated from background chiral current (not shown) is represented by black arrows.The red and blue squares label the clockwise and counterclockwise vertices, respectively, and the color depth indicates the magnitude of the vortex current.

8 FIG. 2 :
FIG.2: Snapshots of the time evolution of local leg current configurations j ℓ,r with V = 3 after the initial rung current excitation.'Charge' current j c r is shown in the rightmost column, whose deviation from mean values quantifies the leg current of the vortices.

FIG. 3 :
FIG.3: Snapshots of time evolution of the charge concentration configurations n c r = n1,r + n2,r with V = 3.

FIG. 4 :
FIG. 4: The net 'charge' current jν and 'charge' concentration n c r in spreading of wave packets with various V .The two counter-propagating wave packets are shown in different colors.(a)-(c) time evolution of jν with V = 0.125, V = 1.125, and V = 3.The maximum points of n c r and jν at t = 30 are shown in (d) and (e).

6 FIG. 5 :
FIG. 5: Time evolution of the von Neumann entanglement entropy for the central bipartition SvN(L/2).The blue line represents the obtained curve by fitting the data via a logarithmic function a ln (t + b) + c with a = 3.1829, b = 3.2637, and c = −1.4585.The system size is L = 80 and V = 3.