• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    Quantum dynamics of Gaudin magnets

    2022-10-22 08:14:44WenBinHeStefanoChesiHaiQingLinandXiWenGuan
    Communications in Theoretical Physics 2022年9期

    Wen-Bin He,Stefano Chesi,Hai-Qing Lin and Xi-Wen Guan

    1 Beijing Computational Science Research Center,Beijing 100193,China

    2 State Key Laboratory of Magnetic Resonance and Atomic and Molecular Physics,Wuhan Institute of Physics and Mathematics,APM,Chinese Academy of Sciences,Wuhan 430071,China

    3 The Abdus Salam International Center for Theoretical Physics,Strada Costiera 11,I-34151 Trieste,Italy

    4 Department of Physics,Beijing Normal University,Beijing 100875,China

    5 NSFC-SPTP Peng Huanwu Center for Fundamental Theory,Xi’an 710127,China

    6 Department of Theoretical Physics,Research School of Physics and Engineering,Australian National University,Canberra ACT 0200,Australia

    Abstract Quantum dynamics of many-body systems is a fascinating and significant subject for both theory and experiment.The question of how an isolated many-body system evolves to its steady state after a sudden perturbation or quench still remains challenging.In this paper,using the Bethe ansatz wave function,we study the quantum dynamics of an inhomogeneous Gaudin magnet.We derive explicit analytical expressions for various local dynamic quantities with an arbitrary number of flipped bath spins,such as: the spin distribution function,the spin–spin correlation function,and the Loschmidt echo.We also numerically study the relaxation behavior of these dynamic properties,gaining considerable insight into coherence and entanglement between the central spin and the bath.In particular,we find that the spin–spin correlations relax to their steady value via a nearly logarithmic scaling,where as the Loschmidt echo shows an exponential relaxation to its steady value.Our results advance the understanding of relaxation dynamics and quantum correlations of long-range interacting models of the Gaudin type.

    Keywords: Gaudin magnets,central spin model,Bethe ansatz,spin polarization,spin–spin correlation,Loschmidt echo

    1.Introduction

    Although mature paradigms,like the celebrated Fermi liquid theory [1–3]and Luttinger liquid theory [4,5],enable us to understand well a wide class of many-body systems in the stationary state,their dynamical evolution is much less understood.In fact,quantum dynamics of many-body systems always presents itself with formidable challenges,due to the difficulty of deriving eigenfunctions analytically and the exponentially growing complexity of numerics.Despite these theoretical difficulties,significant progresses have been made in the last decade,which greatly improves our understanding of dynamical properties of many-body systems.For example,the Kibble–Zurek mechanism [6]reveals power-law relations between the typical length scale ξ of defect domains,the relaxation time τ,and the rate of change of the driving parameter λ,namely ξ~λ-νand τ~λ-zν,which has been observed in ion Coulomb crystals[7].Thermalization of isolated many-body systems is also a difficult topic.The eigenstate thermalization hypothesis(ETH) [8,9]asserts that,for a closed systems evolving to the steady state,the diagonal and micro-canonical ensembles are equivalent.In recent experiments with ultracold atoms[10,11],Kinoshita et al showed that Bose gases in one-dimensional traps do not thermalize,even after thousands of collisions [12].The ETH can be violated by integrable models [13],due to the existence of infinitely many conserved charges,or in the presence of many-body localization[14].The study of dynamics is of crucial importance for many-body quantum systems out of equilibrium,as it provides important clues on if and how their equilibrium state can be reached.

    Therefore,while thermalization usually concerns the steady-state properties of many-body systems after a longtime evolution,an interesting question is on how the manybody system evolves at intermediate times.For a system with a short-range interaction,like the Heisenberg spin chain and the Lieb–Lininger model,light-cone dynamics dominates the relaxation of local observables to steady-state values,such as spin polarization and correlation functions[15].However,the relaxation process of local and global quantities to the steady state is less explored in models with long-range interactions.In this paper,we use Gaudin magnets,i.e.the inhomogeneous central spin models,to study the relaxation dynamics of many-body systems.Such magnets with long-range interactions may describe the decoherence of a spin-qubit due to the interaction with surrounding spins,for example an electron interacting with nuclear spins in quantum dots or defect centers[16,17].Quantum dynamics of the Gaudin model has been extensively studied with a variety of methods as,for example,exact diagonalization and quantum many-body expansions [17–22],the Bethe ansatz method for the polarized bath with one flipped bath spin [23],hybrid methods combining the Bethe ansatz and Monte Carlo sampling[24,25],the DMRG method and semiclassical approaches[26–29].

    Recent work described in detail the collapse and revival dynamics of the homogeneous central spin model [30],and the partition of eigenstates between entangled and separable states [31,32].

    While an early description of Gaudin magnets using the Bethe ansatz technique considered one flipped bath spin,i.e.M=1[23],here we study the time evolution of the central spin problem with an arbitrary number of flipped spins M,starting from an initial statewhere |?〉 is the fully polarized state and ai=0,1,2,...N label the positions of spins in the system.The quantum dynamics of the model with arbitrary number of flipped spins has been already studied through Bethe ansatz [24,25],as well as various other theoretical approaches,leading to important insights such as a description of quantum decoherence by the Chebyshev expansion method[33],non-Markovian dynamics[20],persistent spin correlations [34],etc.Here we study the time evolutions of several interesting quantities analytically and numerically: the spin distribution function,the spin–spin correlation function,and the Loschmidt echo.Evaluating these,we find that the time dependence of the spin polarization function reveals a breakdown of thermalization,such that the steady-state of the Gaudin magnets cannot be described by an equilibrium ensemble.Instead,the relaxation to steady-state values of the spin–spin correlation function suggests a logarithmic decay,and we find that at intermediate times the Loschmidt echo relaxes exponentially to its steady-state value,using both the Bethe ansatz solution and matrix product states (MPS) numerical approach.

    2.Model and Bethe ansatz equations

    The Gaudin magnet describes a central spin at position ‘0’coupled to bath spins through long-range interactions.The Hamiltonian reads

    where the bath spins are labelled with 1→N.For convenience,we parametrize the magnetic field as B=-2/g.Furthermore,we write Aj=1/(∈0-∈j) for the inhomogeneous interaction strength determined in experiments.Setting ∈0=0,the ∈j(with j=1,…,N)correspond to the energy levels of a discrete BCS model associated with H,which can be constructed from the conserved quantitiesHj=Bsjz+∑k≠jsj·sk(∈j-∈k) [35–37].Following the algebraic Bethe ansatz method,and considering a subspace where M spins are flipped-down with respect to the reference state |?〉,the eigenstates are given below [35,36]:

    where the M parameters{να}should satisfy the following Bethe ansatz equations

    with α=1,2,…,M.The equation (3) are also called Richardson–Gaudin equations.There areCNM+1sets of solutions to equation (3),which correspond to the number of choices of flipping M spins [38].TheseCNM+1states |ν1,…,νM〉 span the subspace with M spins flipped-down.In appendix D we give details of the numerical method for solving the Bethe Ansatz equation (3).Moreover,their eigenenergy is given by

    We exploited the above eigenfunctions to study the quantum dynamics of the central spin model.The initial state(see introduction)is a simple product state|Φ0〉=|a1,…,aM〉,where {a1,…,aM} is a set of mutually unequal site indexes which can range from‘0’to‘N’.After the interaction is switched on,the wave function evolves as

    where we introduced the orthonormalized wave functions∣φk〉=Nνk∣ν1,k,…,νM,k〉,withNνkthe normalization coeffciient7See the appendix for details.There we present brief derivations for the time evolution of the spin distribution function,the spin–spin correlation function,and the Loschmidt echo for the Gaudin magnet with arbitrary number of flipped spins M.(see appendix B).The index ‘k’ labels all the states corresponding to the roots of equation(3).Since the initial state has M spins which are flipped-down,we only need to consider the Bethe ansatz basis in this particular subspace.A critical quantity in equation (5) is the overlap between the initial state and the eigenfunctions,given as follows:

    where ‘P’means summing over all permutation of indexes{a1,…,aM},see appendices A and B.The overlap can be also written as 〈Φ0∣φk〉=Nνkperm (1/(να,k-∈Pα)),where ‘perm’is the permanent of matrix 1/(να,k-∈Pα),similar to the determinant [35,36].It is worth noting that both permanent and determinant representation can give the exact quantum dynamics of the Gaudin magnets.However,the numerical tasks from the two representations is quite different as determinants can be evaluated in polynomial time.The interested readers can find more details about permanent and determinant representations in[39–41].For either representation,it still remains challenging to get dynamical quantities for the Gaudin magnets with a large system size,and values up to N=48 were reached [24].

    From expressions (5) and (6),the evolution of observables (for example the spin distribution) can be calculated rather straightforwardly.Based on this approach,we will present in the rest of the paper our key analytical and numerical results on the quantum dynamics of Gaudin magnets.

    3.Spin distribution

    By making use of equation (5),the time-dependent polarization of site j is immediately written as

    wherewkk′=Ek′-Eksimplifies toWhile the overlaps are given in equation (6),the remaining difficulty is to obtain the matrix elements of the observable in the eigenfunction basis.As detailed in appendix B,one can insert the complete set of states |a1,…,aM〉,thus writing the matrix elements in terms of overlaps of type 〈Φ0|φk〉.We finally obtain the exact expression:

    where ,here and throughout the article,‘P’ and ‘Q’ indicate summing over all the permutations of indexes {a1,…,aM}and {j1,…,jM},respectively.

    Although it is not immediately obvious from equation (8),sjz(t) recovers the correct initial value when we take the limit t→0.Detailed proof can be found in appendix B.There we also show that,as expected,mz==(N+1)/2-Mis conserved by equation (8).The evolution of spin distribution function for a more general initial state is just an extension of above result equation (8).

    Figure 1.Time evolution of the spin polarization of individual spins,obtained from equation (3) assuming an exponentially decaying coupling constant Aj=A N exp(-j N).We take N=10 and use dimensionless units,setting A=1.In (a) (c) we choose an initial state|ΦA(chǔ)〉=|0,7,8,9,10〉.In(b)(d)the initial state is|ΦB〉=|0,2,4,6,8〉.Figures (a) and (b) present the dynamical polarization of central spins with different magnetic fields for the two initial states.Figures (c) and (d) show the dynamical evolution of bath spins at B=0.5 for the two initial states.

    We have evaluated equation(8)for two representative initial conditions,i.e.the domain-wall stateand the anti-ferromagnetically ordered stateand present in figure 1 the resulting polarization of individual spins.The dynamical evolution of the central spin is shown in panels(a)and(b),and is allowed by the interaction with the bath spins,which leads to change of the orientation of the central spin through flip-flop processes.It is characterized by a quick decay to the steady-state value,with the persistence of an oscillatory behavior due to finite-size effects [24].There is a competition between magnetic field and coupling to the bath spins,evident from the evolution curves at different values of B shown in figures 1(a) and (b).The physical mechanism behind the suppressed evolution at large B is the increase of the energy gap between states with s0z=± 1/2,which inhibits changes of orientation of the central spin.In appendix C,we also present the detailed derivation of the reduced density matrix and von Neumann entropy for the central spin.The general behavior of the central spin polarization in figure 1 is very similar to the oscillatory time dependence observed for the von Neumann entropy,see figure 7.

    For models with short-range interactions,it is wellknown that the time-evolution of non-equilibrium states is characterized by a light-cone dynamics of local observables[42–45].On the other hand,the Gaudin magnets is a longrange model and,as shown in figures 1(c) and (d),the spin polarization of the bath spins evolve simultaneously into their steady-state value.Being infinite-range,the interaction of the Gaudin magnet can take effect without any propagation time.

    Comparing the two initial states,we see that with |ΦA(chǔ)〉the bath spins need a much longer time than |ΦB〉 to evolve into a steady state.In fact,for the initial state |ΦA(chǔ)〉 the two groups of bath spins remain well distinct,and the oscillation in the spin polarization nearly approach periodically their initial values (1/2 or-1/2).On the other hand,the polarizations for the initial state |ΦB〉 show an initial quick decay,leading to oscillations around 0.The evolution for the initial state|ΦA(chǔ)〉displays a strong memory of the initial domain-wall order.In order to further understand such peculiar dynamics,in later sections we will study the time evolution of the spin–spin correlations and the Loschmidt echo.

    Figure 2.(a)Time average of the central spin polarization,computed for different values of B and N=10.The solid lines are for the initial state |ΦA(chǔ)〉 while the dashed lines are for the initial state |ΦB〉.(b) Dependence of -(i.e.the difference in the time-averaged central spin polarization between the two initial conditions) on the bath spin number N.Here we use a fixed value of B=0.5.The red circles are from the exact Bethe ansatz equations,where as the blue diamonds are obtained by the matrix product state (MPS) method(see section 4).

    4.Spin–spin correlations

    The spin–spin correlation function is defined asGmzn(t)=〈ψ(t)∣smzsnz∣ψ(t)〉and,as discussed in appendix B,which can be computed in a way similar to the spin density.The explicit expression reads:

    which characterizes longitudinal correlation between spins.

    In figure 3,we show the spin–spin correlation function at different times,both for the initial state|ΦA(chǔ)〉(left panels)and|ΦB〉 (right panels).On a relatively short timescale (t<50),the correlation function is very similar to t=0 for both initial states and when the evolution time approaches t~102the correlation function has decreased to smaller values.However,it can be noted that in the left panels of figure 3(referring to|ΦA(chǔ)〉),the domain-wall order has not faded away completely even after a long evolution time t~103.This is in contrast to panel (h) on the right side,where the anti-ferromagnetic alignment has nearly disappeared.This confirms our previous observation that the domain-wall configuration is more favourable to retain memory of the initial state(see also figure 5 for the Loschmidt echo).

    From the contour plot of the correlations in figure 3,the correlations between the bath spins Gmznwith m,n >0 tend to zero.Meanwhile the correlations between the central spin and bath spinsGmz0with m >0 oscillate near initial value.After long enough time evolution (t~103),the correlations between the bath spins almost evolve toward to zero for the state |ΦB〉.

    Considering in more detail the right side of figure 3,we see in panel (g) that the correlations Gmznwith m,n ?5 are more robust,and also the values ofGnz0with larger n remain rather similar to the initial values.This behavior is a consequence of the small coupling strength Ajwhen j is large,which implies a longer timescale for those bath spins.However,after a sufficiently long time evolution (t~103) the correlations between all spins are close to zero.The vanishing of correlations in panel (h) indicates that the steady state is like a paramagnetic state.

    To further analyze the evolution ofGmzn(t),we characterize the typical correlation strength for spin m as

    where the self-correlation terms Gmzmare excluded.In figure 4,we show the absolute value of correlation function and the typical correlation strength for the two initial states,i.e.the left two panels(a),(c)are for the initial state|ΦA(chǔ)〉and the right two panels (b),(d) for |ΦB〉.For both ferromagnetic and anti-ferromagnetic alignment,the correlation strength decays from initial maximum value of 1/4 towards smaller values,with Tm=0 corresponding to disordered spin alignment.We can see that there are two stages in the relaxation of correlations.After the initial decay (up to t10,see the first black dash line),there is an intermediate regime in which the correlation strength has a logarithmic relaxation (up to t~102,near the second black dash line),i.e. ΔTm∝-lnt.In this interval the linear decay shown in figure 4 (on a semilogarithmic scale) has its steepest slope.Finally,in the longtime limit the central spin is strongly entangled with bath spins and the value of Tmreaches its steady-state value.Due to finite-size effects,the decay is modified by oscillations,see especially panel (b),but the qualitative form and timescales are quite robust with respect to the value of m and the type of initial state.

    Figure 3.Left panels: the spin–spin correlation functionGmzm at different times,with a bath size N=10,magnetic field B=0.5,and initial state |ΦA(chǔ)〉=|0,7,8,9,10〉.Right panels: spin–spin correlation functionGmzm with the initial state |ΦB〉=|0,2,4,6,8〉.

    Figure 4.Panels (a) and (b) plot the spin–spin correlation function∣G0 zn∣for spins n=1,N.Panels (c) and (d) average the spin–spin correlation function according to equation (10).The bath size is N=10 and the magnetic field B=0.5.Left panels: initial state|ΦA(chǔ)〉=|0,7,8,9,10〉.Right panels: initial state |ΦB〉=|0,2,4,6,8〉.

    5.Loschmidt echo and MPS approach

    The Loschmidt echo,defined as L(t)=|〈Φ0|ψ(t)〉|2,quantifies the memory of the initial state [46],thus can provide us a clear understanding of the difference between the evolution for the two initial states.By using equation (5),we express L(t) as:

    where the overlaps are given in equation (6).Finally,we obtain:

    We show representative examples of the Loschmidt echo in figure 5.As seen in panel (a),for a weak magnetic field(B~0.5) the initial state |ΦA(chǔ)〉 leads to a long-time evolution with large oscillations around a mean value(t)~0.25.Instead,panel (b) shows that L(t)~0 in the long-time limit,thus the system nearly loses its memory of the initial state|ΦB〉.

    In the first part of the time evolution,besides the presence of fast small amplitude oscillations,the Loschmidt echo displays a clear exponential dependence ΔL(t) ∝exp(-γt),up to an intermediate time scale t~102.According to our numerical calculation,the decay coefficients are similar for both initial states,i.e.γA?γB~0.03 with B=0.5.Such exponential scaling behavior of the Loschmidt echo has been widely found in dynamical evolution of quantum many-body system [47–50].

    Figure 5.Time dependence of the Loschmidt echo for two different initial state.Panel (a) with the initial state |ΦA(chǔ)〉,and panel (b) with initial state |ΦB〉.

    In order to confirm the scaling relaxation behavior of Loschmidt echo for a larger system size,we further perform an approximate calculation based on the MPS method[51,52],as implemented in the Itensor library[53].There are two important input parameters of the MPS simulations: the variational tolerance ∈ of the MPS wavefunction,which represents the error of the MPS state in approximating the true wavefunction (we use ∈=10-9in our simulations),and the bond dimension χ,which indicates the maximum dimension of the matrices entering the variational ansatz.In the given implementation of the MPS method,the bond dimension is automatically increased to meet the desired tolerance.In figure 6,we show the Loschmidt echo for a bath size N=50 and the two types of initial states considered in this work,|ΦA(chǔ)〉 (red lines) and |ΦB〉 (cyan lines).Also in this case the Loschmidt echo is characterized by an initial exponential decay,and this behavior persists up to intermediate times(t ≈102).

    Unfortunately,we are not able to access the long-time dynamics at N=50 based on the MPS approach,as the entanglement in the wavefunction grows with time.This reflects itself on a rapid increase of the bond dimension χ,which eventually reaches numerically intractable values.The time dependence of χ in the simulation of figure 6 is shown in the inset.As seen,the bond dimension χ of initial state |ΦB〉grows faster in time than |ΦA(chǔ)〉.Nevertheless,the range of times we can simulate is sufficient to establish conclusively the initial exponential form of decay of Loschmidt echo.In classical systems,the exponential decay of Loschmidt echo usually implies chaotic dynamics[49,54,55].For the central spin model,the exponential decay of Loschmidt echo may indicate the dynamics is chaotic here as well [56].

    Figure 6.Loschmidt echo as function of time t for two different initial states by using an MPS method with N=50.The red line is for the initial states|ΦA(chǔ)〉,and the cyan line is for the initial state|ΦB〉.

    6.Conclusion

    In summary,based on Bethe Anstaz techniques,we have analytically studied the time evolution of several dynamic quantities of Gaudin magnets with arbitrary number M of flipped-down spins,such as spin distribution function,spin–spin correlation function,and Loschmidt echo.Furthermore,we have numerically studied the scaling behavior of the relaxation dynamics,up to relatively large system sizes.Significantly,the correlation function relaxes to its steady value with a logarithmic dependence,and the Loschmidt echo reveals an exponential loss of the memory of the initial state.Our results highlight several interesting features of the dynamical evolution of a quantum many-body system with long range interactions,thus advancing the understanding of decoherence for single qubits coupled to a spin bath at the many-body level.With respect to universal features of quantum many-body systems,much effort is still necessary to better understand the non-equilibrium quantum dynamics[57].In future studies,it would be desirable to extend our methods to investigate scrambling in many-body system[58,59],and the out-of-time-of-correlation function (OTOC)[60]of central-spin models.The details of such dynamical evolution of spin correlations may be accessible experimentally,e.g.through NMR techniques [61].

    Acknowledgments

    We would like to thank R Fazio for useful discussions and F Iemini and D Ferreira their kind help with the MPS simulation package.WBH acknowledges support from NSAF(Grant No.U1930402).XWG is supported by the key NSFC grant No.12134015 and No.11 874 393,and the National Key R&D Program of China No.2017YFA0304500.SC acknowledges support from NSFC (Grants No.11 974 040 and No.12150610464)and the National Key R&D Program of China No.2016YFA0301200.HQL acknowledges financial support from National Science Association Funds U1930402 and NSFC 11 734 002,as well as computational resources from the Beijing Computational Science Research Center.

    Appendix A.Bethe ansatz basis

    The eigenstates of equation (1) can be obtained with Bethe ansatz(BA)approach,as discussed in detail in[36].Here we discuss several properties of these states which will find wide use in the following derivations.Since the Hilbert space of the Gaudin model can be decomposed into the direct sum of subspaces with M flipped-down spins,we introduce two basis sets for such subspaces.One is the natural basisand another is the Bethe ansatz basisThe corresponding completeness relations are simply given by:

    where ∣Nνk∣are the normalization coefficients for the Bethe ansatz basis.The two sets of states |j1,…,jM〉 and |ν1,k,…,νM,k〉are both complete in the subspace with M flipped-down spins.Including all the M,the complete basis of the whole Hilbert space is obtained.

    The explicit form of the Bethe ansatz eigenfunctions reads

    where ,as in the main text,Q are the permutations of indexes{j1,…,jM}.Note also that any of the two ‘j’ indexes are necessarily unequal.The overlap of state |a1,…,aM〉 with a BA eigenstate|ν1,k,…,νM,k〉can be simply obtained from the above expression as:

    where P are the permutations of indexes {a1,…,aM}.The above formula can be also written in terms of the permanent of the matrix[40].

    The scalar product between eigenstates can be obtained by insertion of

    where we used the overlaps equation (A4) between the two basis sets.The above expression gives the normalization coefficient ∣Nνk∣2.

    We can also derive another useful formula:

    which is obtained by inserting the resolution of the identity in terms of the BA eigenstates.

    Appendix B.Evaluation of physical observables

    We now discuss the explicit evaluation of equation(7),where we recall that∣φk〉=Nνk∣ν1,k,…,νM,k〉are the normalized eigenfunctions.Therefore,for an initial state in the form∣Φ0〉=sa-1…sa-M∣?〉=∣a1,…,aM〉,the overlap matrix elements 〈Φ0|φk〉 are immediately found from equation (A4),and are given in equation (6).Instead,to compute the matrix element ofsjz,we use the fact that such operator is diagonal in the natural basis |j1,…,jM〉 and:

    By inserting the completeness relation in terms of the|j1,…,jM〉states,and computing the overlaps with eigenstates|φk〉as in equation (6),it is easy to find that:

    Finally,combining overlaps and matrix element gives:

    In the above formulas,as usual,P are permutations of the indexes {a1,…,aM} (from the initial state) andQ are permutations of the indexes {j1,…,jM}.The term proportional to 1/2 (from the parenthesis in the first line) can be further simplified as follows.First we note that this term contains a factor of type equation (A5),which is∝δkk′.Therefore,we can simply set coswkk′t=1in this particular contribution.After noticing this,we can use equation (A6) to perform the summations over k andk′ and finally recover the 1/2 appearing in equation (8) of the main text.

    Using a similar method,we can derive the spin–spin correlation function,given by

    While the overlaps are found in equation (6),the matrix elements ofsmzsnzcan be obtained in a way completely analogous to equation (B2):

    leading to equation (9) in the main text.

    Finally,we discuss how to recover the initial conditions,and prove the conservation of magnetization.When time t→0,the spin distribution becomes

    which is in a form suitable to equation (A6).Performing the summations over k andk′ we get:

    which recovers the initial spin distribution.For the total magnetization,computingmz=∑j sjz(t) leads to∑j,αδjjα=Min the second term of equation(7).Finally,the summation overk,k′ can be computed as discussed after equation (B3),leading to the expected result mz=(N+1)/2-M.

    Appendix C.Reduced density matrix of the central spin

    The density matrix ρs=|ψ(t)〉〈ψ(t)| is given by

    Figure C1.Time dependence of the von Neumann entropy at different magnetic fields.Upper subplot:initial state|ΦA(chǔ)〉=|0,7,8,9,10〉.Lower subplot: initial state |ΦB〉=|0,2,4,6,8〉.

    where ,by using equation (A3):

    HereR indicates the permutations of {l1,…,lM}.

    Now we consider the reduced density matrix of central spin,defined byρcs=Tr1→N[ρs].With our choice of initial state,there are no off-diagonal elements and the result is in the form:

    The first contribution arises from the terms of equation (C2)where ji=li>0 (i.e.site 0 is not flipped) and A(t) has the following expression:

    The second contribution corresponds to j1=l1=0 (i.e.the central spin is flipped) and D(t) reads:

    With the help of equations (A5) and (A6),we can prove A+D=1.The von Neumann entropy is simply given asS=-AlnA-DlnD,and a representative plot is shown in figure C1.

    Appendix D.Numerical method for solving BA roots

    In our work we solve the BA equations (3) using the numerical method presented in [38].For convenience of the reader,in this appendix we present the numerical method in detail.For general M,the asymptotic solution to equation (3)with g→0 reads:

    We then introduce the quantities

    and transform Eq.(3) to the following form:

    The advantage of using the Λjvariables is in avoiding the singular dependence of ναon system parameters.With να(0)as initial value,we use Newtonian iteration to find the solution of the above equations at fixed g.

    After obtaining Λj,the νjcan be found as follows.We first transform equation (D2) into

    The above equation (D4) gives a linear system in thepivariables:

    which are elementary symmetric polynomials of the BA parameters ναand p0can be defined as 1.Finally,after obtaining the values of piwe consider the Mth-order polynomial:

    whose zeroes give the desired BA parameters να.An example of the dependence of the BA parameters {ν}on g is shown in figure D1.Close to g=0 the BA parameters ν are all real,but complex conjugate pairs are formed by neighboring roots at larger values of|g|.The interested reader can also see[38]for the description of this method.

    Figure D1.Dependence of BA parameters on g.The top (bottom)panel shows the real (imaginary) part.System parameters: N=15,M=8,and A=2.

    成人国语在线视频| 丰满饥渴人妻一区二区三| 亚洲精品粉嫩美女一区| 一二三四社区在线视频社区8| 国产精品 国内视频| 男女无遮挡免费网站观看| 飞空精品影院首页| 亚洲 欧美一区二区三区| 国产精品秋霞免费鲁丝片| 视频在线观看一区二区三区| 亚洲精品美女久久av网站| 欧美午夜高清在线| 97在线人人人人妻| 日韩免费av在线播放| av不卡在线播放| 黄色毛片三级朝国网站| 国产在线精品亚洲第一网站| 三级毛片av免费| 国产一区二区三区综合在线观看| 自线自在国产av| 亚洲精品成人av观看孕妇| 中文字幕精品免费在线观看视频| 啦啦啦视频在线资源免费观看| 久久久精品94久久精品| 十分钟在线观看高清视频www| 久久99一区二区三区| 男女床上黄色一级片免费看| 日韩欧美一区二区三区在线观看 | h视频一区二区三区| 亚洲中文日韩欧美视频| 精品少妇久久久久久888优播| 中文字幕最新亚洲高清| 免费久久久久久久精品成人欧美视频| 国产麻豆69| 久9热在线精品视频| 亚洲精华国产精华精| 亚洲 国产 在线| 丰满饥渴人妻一区二区三| 美国免费a级毛片| 人人妻人人添人人爽欧美一区卜| 国产精品一区二区在线不卡| 美女扒开内裤让男人捅视频| 欧美成人免费av一区二区三区 | 亚洲精品美女久久av网站| 日韩免费av在线播放| 免费在线观看视频国产中文字幕亚洲| 欧美日韩亚洲国产一区二区在线观看 | 亚洲人成电影免费在线| 欧美日韩成人在线一区二区| 交换朋友夫妻互换小说| 51午夜福利影视在线观看| 肉色欧美久久久久久久蜜桃| 久久人妻熟女aⅴ| 亚洲久久久国产精品| 国产在线视频一区二区| 极品少妇高潮喷水抽搐| 一区二区三区激情视频| 女同久久另类99精品国产91| 在线观看66精品国产| 成年版毛片免费区| 曰老女人黄片| 国产精品98久久久久久宅男小说| 色婷婷久久久亚洲欧美| 99久久99久久久精品蜜桃| 亚洲色图综合在线观看| 黄色a级毛片大全视频| 久久久精品94久久精品| av天堂久久9| 精品少妇久久久久久888优播| av国产精品久久久久影院| 中文字幕av电影在线播放| 久久99热这里只频精品6学生| 国产成人av激情在线播放| 精品国内亚洲2022精品成人 | 国产黄色免费在线视频| 丝袜美腿诱惑在线| 国产精品二区激情视频| 久久久水蜜桃国产精品网| 美女高潮喷水抽搐中文字幕| 亚洲久久久国产精品| 国产一区二区激情短视频| 在线天堂中文资源库| av有码第一页| 极品教师在线免费播放| 亚洲七黄色美女视频| 一级,二级,三级黄色视频| 久久天躁狠狠躁夜夜2o2o| 1024视频免费在线观看| 亚洲欧美一区二区三区黑人| 亚洲精品中文字幕一二三四区 | 两个人看的免费小视频| 制服人妻中文乱码| 亚洲欧洲精品一区二区精品久久久| 最新美女视频免费是黄的| 91大片在线观看| 最新美女视频免费是黄的| 色视频在线一区二区三区| av又黄又爽大尺度在线免费看| 亚洲综合色网址| 亚洲精品国产色婷婷电影| 久久久久久久久免费视频了| 久久人人爽av亚洲精品天堂| 淫妇啪啪啪对白视频| 日本av手机在线免费观看| 丁香六月欧美| 欧美日韩国产mv在线观看视频| 黄色毛片三级朝国网站| 免费日韩欧美在线观看| 97人妻天天添夜夜摸| 久热爱精品视频在线9| 一边摸一边做爽爽视频免费| 少妇粗大呻吟视频| 精品亚洲成a人片在线观看| 天天躁狠狠躁夜夜躁狠狠躁| 80岁老熟妇乱子伦牲交| 深夜精品福利| 国产一区二区激情短视频| 叶爱在线成人免费视频播放| 夜夜爽天天搞| 亚洲五月婷婷丁香| 成人手机av| 热99re8久久精品国产| 一本色道久久久久久精品综合| 熟女少妇亚洲综合色aaa.| 久久国产精品人妻蜜桃| 啦啦啦视频在线资源免费观看| 搡老熟女国产l中国老女人| 在线 av 中文字幕| 黄色 视频免费看| 丁香六月天网| 久久久精品免费免费高清| 国产精品亚洲av一区麻豆| 国产极品粉嫩免费观看在线| 成人国语在线视频| 亚洲综合色网址| 成人18禁高潮啪啪吃奶动态图| 最近最新中文字幕大全电影3 | 热re99久久精品国产66热6| 12—13女人毛片做爰片一| 国产在线一区二区三区精| av在线播放免费不卡| 国产高清激情床上av| 一二三四社区在线视频社区8| 国产男女超爽视频在线观看| 老鸭窝网址在线观看| 国产av又大| 国产成人免费观看mmmm| 免费久久久久久久精品成人欧美视频| 成人18禁高潮啪啪吃奶动态图| netflix在线观看网站| 黄色视频在线播放观看不卡| 天堂俺去俺来也www色官网| 乱人伦中国视频| 日韩视频在线欧美| 欧美黄色片欧美黄色片| 精品人妻熟女毛片av久久网站| 日韩欧美免费精品| 麻豆av在线久日| 中文字幕人妻熟女乱码| 欧美黄色淫秽网站| www.熟女人妻精品国产| 啦啦啦视频在线资源免费观看| 两性午夜刺激爽爽歪歪视频在线观看 | 国产单亲对白刺激| 欧美久久黑人一区二区| 精品国产国语对白av| 一级黄色大片毛片| 亚洲伊人色综图| 久久精品国产亚洲av香蕉五月 | 精品乱码久久久久久99久播| 国产三级黄色录像| 国产精品久久久久久人妻精品电影 | 日本撒尿小便嘘嘘汇集6| 国产精品久久电影中文字幕 | 男人操女人黄网站| 国产单亲对白刺激| 看免费av毛片| 亚洲专区字幕在线| 丰满饥渴人妻一区二区三| 国产日韩欧美亚洲二区| 亚洲成人国产一区在线观看| 国产av精品麻豆| 中亚洲国语对白在线视频| 99热国产这里只有精品6| 亚洲国产欧美在线一区| 老司机福利观看| 99久久精品国产亚洲精品| 国产极品粉嫩免费观看在线| 亚洲伊人久久精品综合| 黄色丝袜av网址大全| netflix在线观看网站| 亚洲熟女毛片儿| 男女午夜视频在线观看| 日韩欧美一区二区三区在线观看 | 午夜精品国产一区二区电影| 巨乳人妻的诱惑在线观看| 午夜福利免费观看在线| 久久免费观看电影| av片东京热男人的天堂| 国产日韩一区二区三区精品不卡| av电影中文网址| 婷婷成人精品国产| 国产成人系列免费观看| 国产精品九九99| 久久久久久久久久久久大奶| 国产熟女午夜一区二区三区| 可以免费在线观看a视频的电影网站| 精品一区二区三区av网在线观看 | 好男人电影高清在线观看| 午夜福利欧美成人| 极品教师在线免费播放| 欧美乱妇无乱码| 色综合欧美亚洲国产小说| 在线观看免费高清a一片| 少妇 在线观看| 丝袜喷水一区| 深夜精品福利| 18禁美女被吸乳视频| 久久香蕉激情| 欧美精品av麻豆av| 热99国产精品久久久久久7| 老司机靠b影院| 久久毛片免费看一区二区三区| 国产精品一区二区在线观看99| 欧美亚洲 丝袜 人妻 在线| av又黄又爽大尺度在线免费看| 精品国产乱子伦一区二区三区| 一边摸一边抽搐一进一出视频| 黑人巨大精品欧美一区二区蜜桃| 精品乱码久久久久久99久播| 久久这里只有精品19| 午夜激情久久久久久久| 一级毛片精品| 国产亚洲av高清不卡| 精品视频人人做人人爽| 欧美日韩中文字幕国产精品一区二区三区 | 欧美 亚洲 国产 日韩一| 999久久久国产精品视频| kizo精华| 男女下面插进去视频免费观看| 精品久久久久久电影网| 色视频在线一区二区三区| 九色亚洲精品在线播放| 丝袜喷水一区| 女人久久www免费人成看片| 国产黄频视频在线观看| 午夜成年电影在线免费观看| 成人永久免费在线观看视频 | 最新的欧美精品一区二区| 中文字幕色久视频| 脱女人内裤的视频| 少妇 在线观看| 国产xxxxx性猛交| 日韩视频一区二区在线观看| 香蕉国产在线看| 在线播放国产精品三级| 免费在线观看视频国产中文字幕亚洲| 99国产极品粉嫩在线观看| 精品人妻熟女毛片av久久网站| 麻豆av在线久日| 免费久久久久久久精品成人欧美视频| 亚洲va日本ⅴa欧美va伊人久久| 十八禁高潮呻吟视频| 亚洲性夜色夜夜综合| 亚洲一卡2卡3卡4卡5卡精品中文| av天堂久久9| 下体分泌物呈黄色| 波多野结衣av一区二区av| 大陆偷拍与自拍| 菩萨蛮人人尽说江南好唐韦庄| xxxhd国产人妻xxx| 女同久久另类99精品国产91| 精品人妻在线不人妻| 日韩欧美免费精品| 国产麻豆69| 色综合婷婷激情| 久久婷婷成人综合色麻豆| 亚洲第一青青草原| 电影成人av| 色综合欧美亚洲国产小说| 精品久久久精品久久久| 嫩草影视91久久| 欧美日本中文国产一区发布| 国产成人影院久久av| 国产成人啪精品午夜网站| 两性午夜刺激爽爽歪歪视频在线观看 | 亚洲中文日韩欧美视频| 欧美日韩亚洲高清精品| 日本精品一区二区三区蜜桃| 亚洲人成电影观看| 国产高清国产精品国产三级| 麻豆成人av在线观看| 日韩欧美一区二区三区在线观看 | 精品国产国语对白av| 在线观看舔阴道视频| 飞空精品影院首页| 91精品三级在线观看| 免费人妻精品一区二区三区视频| 久久久久精品国产欧美久久久| 欧美性长视频在线观看| 精品亚洲成a人片在线观看| 亚洲九九香蕉| 午夜老司机福利片| 国产免费视频播放在线视频| 亚洲专区中文字幕在线| 大香蕉久久网| 久久香蕉激情| 日日夜夜操网爽| 亚洲国产欧美日韩在线播放| 日本黄色日本黄色录像| 别揉我奶头~嗯~啊~动态视频| 午夜免费鲁丝| 久久久久久免费高清国产稀缺| kizo精华| 久久性视频一级片| 91老司机精品| 欧美日韩亚洲高清精品| 国产一区二区三区综合在线观看| 制服诱惑二区| xxxhd国产人妻xxx| 久久久久久人人人人人| videosex国产| 操出白浆在线播放| 国产成人免费无遮挡视频| 韩国精品一区二区三区| 欧美日韩亚洲综合一区二区三区_| 国产一区二区在线观看av| 香蕉久久夜色| 亚洲第一av免费看| 国产福利在线免费观看视频| 亚洲精品国产一区二区精华液| 老熟妇乱子伦视频在线观看| 午夜福利一区二区在线看| 欧美精品人与动牲交sv欧美| 免费av中文字幕在线| 一边摸一边抽搐一进一出视频| 午夜福利乱码中文字幕| 久久久久国内视频| 黄色a级毛片大全视频| 午夜视频精品福利| 亚洲av电影在线进入| 在线av久久热| 在线 av 中文字幕| 国产精品久久久人人做人人爽| 亚洲五月婷婷丁香| 国产精品 欧美亚洲| 在线观看66精品国产| 亚洲第一青青草原| 欧美人与性动交α欧美精品济南到| 黄色a级毛片大全视频| 日韩免费av在线播放| 91成年电影在线观看| 一本—道久久a久久精品蜜桃钙片| 老司机深夜福利视频在线观看| 精品一区二区三卡| 色精品久久人妻99蜜桃| 国产日韩一区二区三区精品不卡| 丝袜在线中文字幕| 日韩熟女老妇一区二区性免费视频| 国产欧美日韩综合在线一区二区| 欧美乱码精品一区二区三区| 少妇精品久久久久久久| 最近最新中文字幕大全免费视频| 法律面前人人平等表现在哪些方面| 亚洲欧美激情在线| 国产欧美亚洲国产| 亚洲免费av在线视频| 亚洲一区中文字幕在线| 久久天堂一区二区三区四区| 国产欧美日韩一区二区精品| 亚洲午夜理论影院| 国产精品久久久av美女十八| 国产精品久久久久久人妻精品电影 | 真人做人爱边吃奶动态| 99久久人妻综合| 天堂中文最新版在线下载| 法律面前人人平等表现在哪些方面| 欧美成人午夜精品| 黄色视频不卡| 色婷婷av一区二区三区视频| 超碰成人久久| kizo精华| 后天国语完整版免费观看| 国产精品电影一区二区三区 | 大型黄色视频在线免费观看| 99精国产麻豆久久婷婷| 日韩中文字幕视频在线看片| 下体分泌物呈黄色| 亚洲成av片中文字幕在线观看| 人人妻人人澡人人看| 高清欧美精品videossex| 午夜福利乱码中文字幕| 青草久久国产| 99国产精品99久久久久| 少妇精品久久久久久久| 中国美女看黄片| 国产精品电影一区二区三区 | 欧美 日韩 精品 国产| 国产一区二区三区在线臀色熟女 | 午夜福利视频在线观看免费| 色精品久久人妻99蜜桃| 国产成人免费无遮挡视频| 欧美av亚洲av综合av国产av| 中亚洲国语对白在线视频| 日本av手机在线免费观看| 久久久国产欧美日韩av| 麻豆av在线久日| 热re99久久国产66热| 国产在线观看jvid| 久久久久国内视频| 中文字幕高清在线视频| 午夜福利欧美成人| 麻豆国产av国片精品| 中文字幕精品免费在线观看视频| 国产不卡av网站在线观看| 热99re8久久精品国产| 欧美日韩av久久| 久久天堂一区二区三区四区| 国产国语露脸激情在线看| av不卡在线播放| 亚洲男人天堂网一区| 亚洲国产中文字幕在线视频| 国产一区二区三区视频了| 啦啦啦 在线观看视频| 亚洲成a人片在线一区二区| 欧美精品亚洲一区二区| 亚洲熟妇熟女久久| 久久久久久亚洲精品国产蜜桃av| 精品一区二区三区视频在线观看免费 | a在线观看视频网站| 纵有疾风起免费观看全集完整版| 成人18禁高潮啪啪吃奶动态图| 动漫黄色视频在线观看| 久久国产精品大桥未久av| 亚洲精品国产色婷婷电影| 自线自在国产av| 国产精品久久久久久人妻精品电影 | a级毛片黄视频| 久久久久精品人妻al黑| 精品欧美一区二区三区在线| 久久久精品区二区三区| 丝袜美腿诱惑在线| 久久热在线av| 国产日韩欧美在线精品| 十八禁人妻一区二区| 王馨瑶露胸无遮挡在线观看| 在线永久观看黄色视频| 99久久国产精品久久久| 国产精品一区二区精品视频观看| 国产成人精品在线电影| 午夜成年电影在线免费观看| 丝瓜视频免费看黄片| 亚洲精品自拍成人| 少妇猛男粗大的猛烈进出视频| 国产免费av片在线观看野外av| 亚洲国产av影院在线观看| 大片电影免费在线观看免费| 亚洲伊人色综图| 视频在线观看一区二区三区| 在线观看舔阴道视频| 人人妻人人爽人人添夜夜欢视频| 性少妇av在线| 黄网站色视频无遮挡免费观看| 精品国产乱码久久久久久男人| 欧美久久黑人一区二区| 啦啦啦免费观看视频1| 亚洲天堂av无毛| 国产在视频线精品| 久久久久精品人妻al黑| 美女扒开内裤让男人捅视频| 狠狠精品人妻久久久久久综合| 亚洲黑人精品在线| 亚洲中文日韩欧美视频| 精品乱码久久久久久99久播| 可以免费在线观看a视频的电影网站| 黄网站色视频无遮挡免费观看| 狠狠狠狠99中文字幕| 久久 成人 亚洲| 欧美变态另类bdsm刘玥| 国产97色在线日韩免费| 成人特级黄色片久久久久久久 | 免费在线观看黄色视频的| 日韩一区二区三区影片| 国产免费福利视频在线观看| 国产在线精品亚洲第一网站| 久久久精品免费免费高清| aaaaa片日本免费| 老司机靠b影院| a级毛片在线看网站| 欧美人与性动交α欧美软件| 国产在线视频一区二区| 午夜久久久在线观看| 免费看a级黄色片| 波多野结衣一区麻豆| 亚洲五月色婷婷综合| 91av网站免费观看| 亚洲国产欧美网| 久久天堂一区二区三区四区| 人妻久久中文字幕网| 亚洲少妇的诱惑av| 免费观看a级毛片全部| 变态另类成人亚洲欧美熟女 | 久久久久久久久久久久大奶| 精品卡一卡二卡四卡免费| 一二三四社区在线视频社区8| 亚洲av欧美aⅴ国产| 国产成人免费观看mmmm| 久久久久久久国产电影| 人人妻,人人澡人人爽秒播| 成年人免费黄色播放视频| 亚洲五月婷婷丁香| 最近最新中文字幕大全电影3 | 欧美日韩亚洲国产一区二区在线观看 | 久久婷婷成人综合色麻豆| 少妇被粗大的猛进出69影院| 日韩欧美免费精品| 高清av免费在线| 亚洲精品在线美女| 国产成人系列免费观看| 国产欧美亚洲国产| 99香蕉大伊视频| 一边摸一边抽搐一进一小说 | 亚洲精品成人av观看孕妇| 久久精品国产a三级三级三级| 欧美变态另类bdsm刘玥| 欧美日韩中文字幕国产精品一区二区三区 | 国产91精品成人一区二区三区 | 真人做人爱边吃奶动态| 久久人妻熟女aⅴ| 女同久久另类99精品国产91| 99香蕉大伊视频| 日韩成人在线观看一区二区三区| 国产高清videossex| 一二三四社区在线视频社区8| 亚洲中文日韩欧美视频| 亚洲七黄色美女视频| 午夜两性在线视频| 国产精品秋霞免费鲁丝片| 99re在线观看精品视频| 天堂8中文在线网| 成在线人永久免费视频| 亚洲欧美精品综合一区二区三区| 亚洲一卡2卡3卡4卡5卡精品中文| 色在线成人网| 制服人妻中文乱码| 国产精品一区二区在线观看99| 美女国产高潮福利片在线看| 欧美日韩视频精品一区| 丝袜喷水一区| 国产日韩一区二区三区精品不卡| 免费高清在线观看日韩| 欧美精品啪啪一区二区三区| 嫁个100分男人电影在线观看| av线在线观看网站| 午夜成年电影在线免费观看| 久久久久久久国产电影| 欧美日韩中文字幕国产精品一区二区三区 | 热re99久久国产66热| 精品福利观看| 亚洲国产毛片av蜜桃av| 日韩视频一区二区在线观看| 日韩有码中文字幕| 真人做人爱边吃奶动态| 波多野结衣av一区二区av| 国产人伦9x9x在线观看| 一个人免费在线观看的高清视频| 麻豆成人av在线观看| 日本撒尿小便嘘嘘汇集6| 考比视频在线观看| 久久 成人 亚洲| 精品福利永久在线观看| 久久久久久久大尺度免费视频| 大码成人一级视频| 99re在线观看精品视频| 国产黄频视频在线观看| 高潮久久久久久久久久久不卡| 亚洲精品中文字幕在线视频| www.999成人在线观看| 夜夜夜夜夜久久久久| 老熟妇乱子伦视频在线观看| 亚洲国产看品久久| 国产一区二区三区综合在线观看| 日本av免费视频播放| 日日夜夜操网爽| 国产男靠女视频免费网站| 美女主播在线视频| 美女午夜性视频免费| 多毛熟女@视频| 亚洲专区中文字幕在线| 精品人妻在线不人妻| 亚洲久久久国产精品| tocl精华| 王馨瑶露胸无遮挡在线观看| 国产高清激情床上av| 免费在线观看黄色视频的| 啦啦啦中文免费视频观看日本| 黄色 视频免费看| 日韩三级视频一区二区三区| 久9热在线精品视频| 91九色精品人成在线观看| 最近最新免费中文字幕在线| av线在线观看网站| 9色porny在线观看| 18禁美女被吸乳视频| 老司机靠b影院| 亚洲精品国产一区二区精华液| 91成年电影在线观看| 亚洲五月色婷婷综合| 夜夜夜夜夜久久久久| 国产欧美日韩精品亚洲av| 美女国产高潮福利片在线看| 久久这里只有精品19| 深夜精品福利| 久久久久网色| 国产极品粉嫩免费观看在线|