Excited states of holographic superconductors from massive gravity

Qian Xiang,Li Zhao,*and Yong-Qiang Wang

1 Institute of Theoretical Physics,Lanzhou University,Lanzhou,Gansu 730000,China

2 Lanzhou Center for Theoretical Physics and Key Laboratory of Theoretical Physics of Gansu Province,Lanzhou University,Lanzhou,Gansu 730000,China

Abstract In this paper,we generalize the study of the model of holographic superconductors in excited states to the framework of massive gravity at the probe limit.By taking into account the effect of a massive graviton,we numerically present a family of solutions for holographic superconductors in excited states and find that the critical temperatures can be higher due to the effect of the massive graviton,in comparison with the superconductor in Einstein gravity.We also investigate the condensates and conductivities in the ground state and the excited states by studying various parameters that determine the framework of gravity background.

Keywords:AdS/CFT,massive gravity,holographic superconductors

1.Introduction

In condensed matter physics,Landau-Ginzburg’s phenomenological theory and BCS theory suffer from some difficulties in fully explaining the mechanism of superconductivity.Due to the strong coupling,electrons can still pair up in some materials with higher critical temperatures,which the BCS theory failed to predict.Fortunately,as one of the greatest achievements of string theory,the AdS/CFT correspondence which establishes the relations between strongly correlated fields on the boundary and weak gravity in a one-dimensional higher bulk spacetime has provided a novel way to study condensed matter physics.In order to model superconductivity,Horowitz et al[1—4]coupled a complex scalar field to a U(1)gauge field.Considering a spontaneous U(1)symmetry breaking under the critical temperature,they successfully modeled the Cooper pair-like superconductor condensate with the scalar field,which provides us with a new way of understanding the mechanism of high temperature superconductivity.Afterward,the holographic superconductor condensates were realized in various ways.By replacing the complex field with a symmetric and traceless second-rank tensor,the holographic d-wave model was constructed in[5—7].The holographic p-wave superconductor condensate can be modeled by a two-form field[8]or a complex,massive vector field with U(1)charge[9,10].In addition,the holographic p-wave superconductor can also be constructed by the SU(2)Yang-Mills field coupling to gravity[11].For reviews of holographic superconductors,see[12—15].

Recently,inspired by the seminal paper of Gubser et al[1],Yong-Qiang Wang et al generalized the studies of holographic superconductors from ground states to excited states.In[16],the authors found a family of solutions of excited states for the s-wave holographic superconductor in the probe limit,where the scalar field possesses n nodes along the radial coordinate corresponding to the n-th excited state.Moreover,in[17],considering the holographic superconductor with the full backreaction,the author also numerically studied the condensate and conductivity of their ground state and excited states.It was found that the highly excited states of〈1O 〉condensates converge to about 4.4 as the temperature approaches zero.It was also found that the condensation values ofO2operators in the excited states are larger than in the ground state.Other studies on the excited state of holographic superconductor can be found in[18,19],where in[18],the author proposed an analytic technique based on a variational method for the Sturm-Liouville eigenvalue problem,to investigate the excited states of the holographic superconductor in the probe limit.The non-equilibrium process of the holographic s-wave superconductor with the excited states was investigated in[19].

In the above discussion,the Cooper pair-like superconductor condensate was modeled by matter fields coupling to Einstein gravity.As a holographically modeled superconductor,it is necessary to involve momentum dissipation to model features of a material in the experiment.In[20],the authors pointed out that there exists a connection between momentum relaxation and graviton mass.Besides,as a natural curiosity,it is always interesting to ask how will such a system behaves,if the graviton is massive.However,this is not an easy task,since the diffeomorphism invariance in GR has already restricted the graviton to be the massless spin-2 boson[21].Simple generalization would usually make the massive gravity unstable because of the problem of the Boulware-Deser ghost[22].To avoid this,C.de Rham and G.Gabadadze(dRGT)introduced polynomial terms in a general action,successfully building a ghost-free nonlinear massive gravity theory[23—25].Such a novel achievement arouses great concerns in the realm of holographic superconductivity.The holographic DC and Hall conductivity was studied in[26]in massive Einstein-Maxwell-Dilaton gravity providing a holographic method to explain strange metals.The nonequilibrium process of a holographic s-wave superconductor was investigated in the dRGT massive gravity in[27].In[28],Cao H.Nam investigated the effect of massive gravity on the p-wave holographic superconductor and found that the critical temperature and condensate depend crucially on the sign of the massive gravity couplings.In addition,attempting to avoid complicated numerical calculation,a framework for translational symmetry breaking and momentum dissipation was set up in[29].In[30],Ya-Peng Hu et al studied a holographic Josephson junction of its critical temperature,tunneling current and coherence length,and found that due to the graviton mass,the transition from a normal state to a superconducting state would be more difficult.The model of holographic superconductor with backreaction in massive gravity can be found at[31].

However,until now,the investigation of the holographic superconductor of its condensate and conductivity in the ground state and the excited states has not yet been studied in the framework of massive gravity.Therefore,in this paper,we will numerically investigate the holographic superconductor of its condensate and conductivity in the ground state and the excited states considering the effect of the dRGT massive gravity.

The paper is arranged as follows:in section 2 we give a brief review of the dRGT massive gravity.The construction of holographic superconductor in(3+1)-dimensional AdS spacetime with massive graviton can be found in section 3.The numerical results and properties of scalar condensate and optical conductivity are in section 4.The conclusion and discussion are in the last section.

2.Review of dRGT Massive Gravity

In order to construct a holographic superconductor in the framework of massive gravity we consider the following action in a 4-dimensional ghost-free dRGT massive gravity,which is the Hilbert-Einstein action including the nonlinear polynomial terms[23]

where

Since we study the holographic superconductor in a 4-dimensional background,we choose the reference metric as:

where c0is a positive constant that,for convenience,is set to one.Such that,a general black hole solution can be found[32,33]

with

Where k is a parameter usually used in topological black holes allowing one to study a hyperbolic,Ricci flat or sphere black hole horizon with its value chosen to be -1,0,1.Here,we choose k=0.The line element hijdxidxjis related to a 2-dimensional flat space.The m0is the mass of the black hole.In(2.7),two independent parameters c1and c2emerging from the ghost-free massive graviton self-interacting potential can be realized as couplings of massive graviton λ in the holographic superconductor.By tuning the product of c1λ2/2 and c2λ2which serve as a linear term of the radius and a constant term in the metric solution,one can study the effect of the dRGT massive gravity on the holographic superconductor.Besides,it is obvious that the solution will go back to a planar Schwarzschild-AdS black hole if one sets c1=c2=0 or simply λ=0.The temperature is introduced by the black hole temperature on the radius rh=1 horizon as below:

In this paper,we will study the holographic condensate and conductivity via tuning the coupling factors c1,c2and the graviton mass λ.

3.Holographic setup

In order to model the scalar condensate at low temperature,we consider the following Maxwell and complex scalar field action in AdS4as

in which A is the U(1)gauge field with the corresponding field strength Fμν=?μAν-?νAμ.α is the mass of the complex scalar field ψ.Note that,we have rescaled the above action with the charge of the scalar field and therefore the backreaction of the matter field on the metric can be negligible if the charge is large enough.

From the above action,the EoMs can be obtained as

On the two-dimensional flat space,an isotropic ansatz for the U(1)gauge field and the scalar field can be written as

With the above ansatz,the equations of motion for the scalar field ψ(r)and electrical scalar potential φ(r)in the background of the Schwarzschild-AdS black hole are

Near the infinite boundary r →∞,the asymptotic behaviors of the functions ψ(r)and φ(r)have the following forms

with

Where the condensate is modeled by the expectation value of scalar operatoriO dual to the field ψ on the r →∞boundary.μ,ρ are the chemical potential and charge density of the holographic superconductor,respectively.Usually,by numerically setting one of the ψ(i)to be zero on the r →∞boundary,the other can be used to describe the corresponding expectation value of the scalar operator with dimension Δ±without being sourced.In this paper,we choose α2=-2 to meet the Breitenlohner-Freedman(BF)bound condition that α2≥-9/4[34].And hence,the dimensions of the condensate operators are Δ+=2 and Δ-=1,respectively.We will study both〈1O 〉 and〈O2〉 condensates and conductivity corresponding to〈O2〉in the following section.

4.Numerical results

In condensed matter physics,when temperature T drops below a critical temperature Tc,electrical resistivity will suddenly drop to zero leading to superconductivity.In the holographic approach,due to the spontaneously broken U(1)gauge symmetry,the scalar condensate will turn on when the chemical potential μ(or ρ)is larger than μc(or ρc).In this paper,the temperature is defined by the ratio between Hawking temperature on the horizon and ρ-1/2,such that Tcis represented by TBHdivided byThe dimensionless relation between temperature and condensate can be presented.In this section,we will numerically study the holographic superconductor in the framework of massive gravity of its condensate and conductivity in ground and excited states.

Firstly,we conduct a simple coordinate transformation,which reads:

We will set the spatial infinity z∞=0 and the black hole radius zh=1 for simplicity.After imposing the above transformation,we conduct all of our numerical calculations based on the spectral method.In the integration region 0 ≤z ≤1,our choice of typical grids have sizes from 50 to 300.

In order to study the scalar condensate of its ground and excited states for both1O andO2operators,the coupled equations(3.13)and(3.14)will be put into an iterative process by means of the Newton-Raphson method.We will give a good initial guess that leads to a solution for a ground or nexcited state for each time.For the ground state of〈iO 〉condensate,the initial guess for ψ(z)is the profile that has no node along the z coordinate.For the n-th excited state,we choose the initial guess that allows ψ(z)to have n nodes on the z axis to reach the excited state.The estimated relative error for our numerical results is below 10-5.

Recall that the effect of massive gravity on the holographic superconductor is involved in the black hole solution(2.7)performed by the coupling coefficients c1,c2and the graviton mass λ.Numerically changing their values,will enable us to study the behaviours of such a holographic system,its scalar condensate and its conductivity under different configurations of massive gravity framework.However,the contribution of massive gravity will disappear and the situation will go back to[2,16]by setting these parameters to be c1=0,c2=0,or simply,λ=0.Therefore,in the study of the effect of c1(or c2)we change c1(or c2)only with c2=0(or c1=0)and λ ≠0.As for studying graviton mass,we change λ only while maintaining c1=1,c2=-0.5.

4.1.Critical Temperature

In figure 1,we present our numerical results of critical temperatures for bothO1(first row)andO2(second row)operators,where from left to right they correspond to the cases of tuning c1with fixed c2=0,λ=1,tuning c2with fixed c1=0,λ=1 and tuning λ with fixed c1=1,c2=-0.5,respectively.

From the figure,we can see that the Tcof excited states(red and blue lines)corresponding to scalar field solutions that possess nodes along radius direction are all lower than ground states,while the Tcof higher excited states is even lower than their former states.This means that by decreasing the temperature,the holographic superconductor will reach superconductivity firstly through the ground state,then new bound states will begin to appear with lower temperatures.Recall that the critical temperature is equal to the ratio of the Black hole temperature andthat both depend on c1,c2,λ.Thus,the value of Tcis no longer proportional toalone as it was in the Einstein-Maxwell framework.An interesting finding from the results is the monotonicity of the ground state of Tcfor〈O1〉(black lines in the first row of figure 1)will change during the process of c1,c2,λ tuned from small to large,individually.However,it is notable that Tcof ground states of〈O1〉 cannot be infinitely high by simply lowering these parameters that determine the value of TBH,because sufficiently small values of c1and c2will cause minus black hole temperature.On the other hand,except for the curves of ground states of〈O1〉,as we can see from the remaining curves,the Tcof each state and each operator increases steadily when the parameters are tuned from small to large.

In tables 1—3,we present our results of critical chemical potentialfrom ground state to fifth excited state.In tables 1 and 2,we investigate the coupling coefficients c1and c2changing from -1 to 3 solely.In table 3,we investigate the graviton mass changing from 0.1 to 0.5 solely.

From the tables 1—3,we find that the μcof ground states,which corresponds to the ground states of ψ that has no node along the radius,is the smallest for fixed coupling factors and graviton mass in both〈O1〉 and〈O2〉 condensates.It is the same for the results of Einstein gravity,the critical chemical potentials of higher excited states are all larger than their former states.By enlarging the effect of massive gravity,μcof each state and each operator will have a larger value.

Table 1.The table refers to critical chemical potential of the holographic superconductor from ground state to fifth excited state for the coupling coefficient c1,where only c1 changes its values and c2 and λ remain fixed.

Table 1.The table refers to critical chemical potential of the holographic superconductor from ground state to fifth excited state for the coupling coefficient c1,where only c1 changes its values and c2 and λ remain fixed.

iO〈 〉 c1 c2 λc 0μc 1μc 2 μc 3 μc4μc 5 μ 1O〈 〉 -1 0 0.2 0.668 6.383 11.585 16.759 21.927 27.092 1 0 0.2 1.429 6.603 11.816 17.037 22.262 27.488 3 0 0.2 1.884 6.820 12.045 17.312 22.593 29.710 2 O〈 〉 -1 0 0.2 4.014 9.113 14.251 19.400 24.554 29.710 1 0 0.2 4.114 9.263 14.462 19.676 24.897 30.120 3 0 0.2 4.213 9.411 14.672 19.950 25.237 30.527

Table 2.The table refers to critical chemical potential of the holographic superconductor from ground state to fifth excited state for the coupling coefficient c2,where only c2 changes its values and c1 and λ remain fixed.

Table 2.The table refers to critical chemical potential of the holographic superconductor from ground state to fifth excited state for the coupling coefficient c2,where only c2 changes its values and c1 and λ remain fixed.

iO〈 〉 c1 c2 λc 0μc 1μc 2 μc 3 μc4μc 5 μ 1O〈 〉 0 -1 0.2 1.098 6.456 11.642 16.817 21.989 27.161 0 1 0.2 1.142 6.530 11.759 16.980 22.199 27.418 0 3 0.2 1.184 6.603 11.875 17.141 22.406 27.672 2 O〈 〉 0 -1 0.2 4.037 9.140 14.287 19.445 24.681 29.775 0 1 0.2 4.091 9.236 14.527 19.632 24.842 30.055 0 3 0.2 4.144 9.329 14.565 19.815 25.072 30.332

and the fittings of coupling coefficient c2,listed in table 2

Similarly,the fittings of graviton mass λ,listed in table 3

A previous study showed that,in the Einstein gravity model,the difference of μcbetween consecutive states is approximately 5 for both〈O1〉 and〈O2〉,indicating a linear relation between μcand n[16].This interesting finding was then found in a semi-analytical study[18].Here,from(4.19)to(4.24),we can see that this property also appears when taking into account the effect of massive gravity.

4.2.Condensate

In this subsection,we study the condensates of our holographic superconductor for bothO1andO2operators in the framework of massive gravity.According to AdS/CFT dictionary,the expectation value of〈iO 〉dual to ψ(i)is

As usual,we study the effects of a massive graviton on scalar condensates by numerically solving equation(3.13),equations(3.14)and(4.25)with changing parameters c1,c2and λ solely.

In figure 2,we present our condensation results in various configurations of a massive gravity background,where the〈O1〉and〈O2〉condensates are presented in the left and the right panel,respectively.In the first row,we fix c2=0,λ=1 and study the coupling coefficient c1in the cases of c1=0.2,c1=1.0,c1=2.0.In the middle row,we also study the coupling coefficient c2in the same way while setting c1=0,λ=1.In the last row,we fix c1=1,c2=-0.5 while tuning the graviton mass from λ=0.1 to λ=1.8.

Table 3.The table refers to critical chemical potential of the holographic superconductor from ground state to fifth excited state for the graviton mass λ,where only λ changes its values and c1 and c2 remain fixed.

Table 3.The table refers to critical chemical potential of the holographic superconductor from ground state to fifth excited state for the graviton mass λ,where only λ changes its values and c1 and c2 remain fixed.

iO〈 〉 c1 c2 λc 0μc 1μc 2 μc 3 μc 4μc 5 μ 1O〈 〉 1 -0.5 0.1 1.203 6.516 11.722 16.923 22.123 27.324 1 -0.5 0.3 1.712 6.698 11.894 17.119 22.353 27.591 1 -0.5 0.5 2.371 7.053 12.236 17.508 22.809 28.123 2 O〈 〉 1 -0.5 0.1 4.073 9.201 14.375 19.562 24.754 29.950 1 -0.5 0.3 4.146 9.303 14.516 19.744 24.980 30.219 1 -0.5 0.5 4.291 9.506 14.796 20.108 25.429 30.756

From the three plots of〈O1〉condensates,we can see that the condensates of ground states rise rapidly at T/Tc=1 and then go to infinity at lower temperatures.Such a singularity indicates that the charge of the scalar field is not large enough at low temperatures and consideration of the backreaction is needed.The excited states also begin to condense around Tc,however,with the condensation values converge to limited values at low temperatures.We can also see that the condensation values of the ground states are larger than the excited states when the coupling coefficients c1,c2and the graviton mass λ are small.However,when the effect of massive gravity is stronger,the condensation values of excited states can be equal or even less to or than ground states.Moreover,from the sub-figures shown inside the first and third plots in the left panel,one can find that,with larger coupling coefficients and graviton mass,the condensation values of the second excited states are higher than the first excited states when the temperature is not sufficiently low.At T →0,the scalar condensates of the second excited states converge to lower values compared to the first excited states.

For the〈O2〉 condensates in the three plots of the right panel,the condensation values of all three plots eventually converge to constants at low temperatures.There exists a gravitational version of the Higgs mechanism which creates connections between the lattice or impurity of a modeled material and graviton mass leading to momentum relaxation rate Γ ～λ2[20].The results presented in figure 2 illustrate that the effect of dRGT massive gravity will lower condensation values of all states when couplings and graviton mass are tuned large.In addition,another feature is that the condensation values of higher excited states are larger than their former states which is similar to the massless gravity case.However,it is worth mentioning that even though the excited states have larger condensates since these are metastable and will last for a considerable time,they will eventually evolve to ground states[19].

It is obvious that our model also possesses a second order transition from normal state to superconducting state near the critical point,which is a square root behaviour predicted by the mean field theory.Therefore,we fit the condensation curves versus temperatures for both O1and O2operators near the critical temperatures according to our results shown in figure 2.

The formats of the fittings are as follows:

where ζ(n)andare the fitting coefficient and critical temperature of the n-th state.The results of ζ(n)are presented in table 4.The corresponding critical temperatures are shown in figure 1 where we have also calculated other critical temperatures in different massive gravity backgrounds to show their tendencies.

From the fitting coefficients of〈O2〉condensate listed in table 4,one can easily find that,similar to the models in Einstein gravity,each ζ(n)of the n-th state is higher than its former state in every case.And,as we enlarge the coupling coefficients or the graviton mass,ζ(n)of each state decreases rapidly for both O1and O2operators.The difference between our model in the massive gravity framework and the model in Einstein gravity is reflected by the〈O1〉condensate,wherein in the latter model the fitting coefficient ζ(n)of the n-th state is smaller than its former state around the critical temperature.Whereas,on the contrary,in our model around the critical temperature,the values of the scalar field condense larger in the n-th excited state than its former state as shown in cases i to iii,the cases viii and ix.One can also directly discover such a tendency by looking at the first and the third plot around Tcin the left panel of figure 2.As for the four remaining cases from iv to vii,we can see that the ζ of excited states are smaller than the ground states with ζ(2)being slightly larger than ζ(1).In addition,as we extend our calculations to a fourth excited state for the cases iv to vii,ζ(n)of the excited state increases slightly with the n-th excited state.

4.3.Conductivity

In this subsection,we study the effects of massive gravity on the conductivity in the ground and the excited states of our holographic superconductor.

By turning on perturbations of the vector potential Axin the bulk geometry of a Schwarzschild-AdS black hole,we consider the Maxwell equation with a time dependence of e-iωt,the linearized equations are given as

When imposing the ingoing boundary conditions at the horizon,the asymptotic behaviour of the Maxwell field on the boundary is

The conductivity can be obtained according to the Ohm’s law

In figures 3 and 4,we plot the real and imaginary parts of optical conductivity as a function of frequency at low temperature T/Tc≈0.100 for the operatorO2of ground states and excited states.In both figures,we use black lines,red lines and blue lines to represent the ground states,first and second excited states.The studied parameters from large to small sequencing are denoted by short dashed lines,long dashed lines and solid lines.In both figure 3 and figure 4,plots.(a)and(b),we study the cases c1=-1,c1=1,c1=3 with c2and λ fixed at c2=0,λ=0.2 and the cases c2=-1,c2=1 c2=3 with c1and λ fixed at c1=0,m=0.2 respectively.In plot.(c),we study the cases λ=0.1,λ=0.3,λ=0.5 with c1and c2fixed at c1=1,c2=-0.5.

Table 4.Fitting coefficients ζ(n) for both〈1O 〉and〈O2 〉condensation curves.

As it is known to all that in the real,dissipative part of conductivity,the horizontal lines at Re(σ)=1 correspond to temperatures higher than the critical temperature Tcwhere there is no scalar condensate.In the three plots of figure 3,the gaps of ground states(black lines)open at ω/T ≈60 until the Re(σ)exponentially converges to 1.Comparing with the Re(σ)in the ground states(black lines),we can see one additional peak for the first excited state(red lines)and two additional peaks(blue lines)for the second excited state.The appearance of additional peaks is the result of excited states which are bound states where ψ has n nodes along the radius.Similarly,in[4]additional peaks of conductivity with the scalar field mass equaling to 3/2 were also found and explained that it is the result of interactions between quasiparticles that form new bound states excited beyond the ground state of BCS theory.Hence,the appearance of new bound states will lead to additional peaks.Besides,since we study the optical conductivity at T/Tc→0.100,the peaks of the ground states develop to almost delta functions as they are generated to excited states.Moreover,the number of the additional peaks of the n-th excited state is equal to n.To explain this,on the one hand,one can ascribe this to the monotonicity’s change of ψ along a radius that affects massive gravity,tuning these parameters will allow us to modify Re[σ]and Im[σ]on the frequency space.in(4.28),where monotonicity of this term alters once for the ground state,three times for the first excited state and five times for the second excited state.On the other hand,by plotting the profile of -Re(Ax)(whose value of the first derivative corresponds to the value ofIm(σ))for excited states around nodes ofIm(σ),one can indeed see the derivative value of -Re(Ax)changes its sign around the nodes.

In both figures,compared with Einstein gravity,where conductivity,of course,develops only one curve in the frequency space with fixed temperature,under the influence of

4.4.Conclusions and discussions

In this paper,we investigated a holographic superconductor that was constructed by a Maxwell field coupled to a scalar field in dRGT massive gravity.The effects of massive gravity are involved in the black hole solution carried by coupling factors c1,c2and graviton mass λ in(2.7).By tuning their values,we numerically studied critical chemical potential,critical temperature,condensate,and conductivity of ground and excited states changing with the dRGT nonlinear massive gravity.

As we tuned c1,c2,λ separately and thus tuned the effect of massive gravity,the change of μc,Tc,condensation value and conductivity of ground and excited states are as follows:

· For〈O1〉 condensate,as couplings c1,c2and graviton mass λ were tuned from small to large in a range that guarantees a positive TBH,the Tcof ground states declined from high values to minimums and then went upward.The Tcof higher excited states,lower than their former states,however,rose steadily with the growth of these parameters.For〈O2〉 condensate,the Tcof ground and excited states all went higher monotonically with larger coupling factors and graviton mass.What’s more,the difference of μcbetween consecutive states is also about 5 for bothO1andO2operators,similar to the discovery in Einstein gravity.

· Under the same configuration of dRGT massive gravity coefficients,condensation values of scalar fields in higher excited states are larger than their former states,for〈O2〉.For〈O1〉,when the effect of massive gravity is not manifest,the condensation value of its ground state is larger than its excited states,while the ground state’s condensation value can be smaller than the excited states with the coupling factors and the graviton mass becoming bigger.Moreover,the stronger effect of massive gravity will reduce the condensation values for both〈O1〉and〈O2〉.

· When studying conductivity,we fixed temperature T/Tc≈0.100 and then tuned c1,c2,λ from small to large.With bigger parameters,the optical conductivities of each state shift to lower frequencies.Besides,similar to the massless gravity case,we also find that the additional peaks(comparing with theRe(σ)andIm(σ)of ground state)of n-th excited state are equal to n.

There could be many interesting extensions of our work.First,as we have noticed in figure 2 the large scalar field charge limit could not suffice for the ground state of〈O1〉,which required us to solve the coupled differential equations with an Einstein equation.Second,around the critical temperature,the condensation of the O1operator performs differently compared with the model in Einstein gravity and the mechanism is unclear.It is interesting to answer this question with a semi-analytical method[18].In the end,we would like to extend our study of excited states in massive gravity to the p-wave and d-wave holographic superconductors in the future.

Acknowledgments

We would like to thank Tong-Tong Hu for the helpful discussions.Parts of computations were performed on the shared memory system at the Institute of Computational Physics and Complex Systems in Lanzhou University.This work was supported by the National Key Research and Development Program of China(Grant No.2020YFC2201503)and the National Natural Science Foundation of China(Grant No.12 047 501).

ORCID iDs