Diffusion in higher dimensional SYK model with complex fermions

We construct a new higher dimensional SYK model with complex fermions on bipartite lattices. As an extension of the original zero-dimensional SYK model, we focus on the one-dimension case, and similar Hamiltonian can be obtained in higher dimensions. This model has a conserved U(1) fermion number Q and a conjugate chemical potential \mu. We evaluate the thermal and charge diffusion constants via large q expansion at low temperature limit. The results show that the diffusivity depends on the ratio of free Majorana fermions to Majorana fermions with SYK interactions. The transport properties and the butterfly velocity are accordingly calculated at low temperature. The specific heat and the thermal conductivity are proportional to the temperature. The electrical resistivity also has a linear temperature dependence term.


Introduction
The Sachdev-Ye-Kitaev model is a strongly interacting quantum system at low energy [1,2]. As the analysis of the recent works [3,4,5,6,7,8,9,10], this model have some interesting properties. It is solvable at large N limit with an emergent conformal symmetry [11]. The reparameterization and gauge symmetries provide connection of SYK to the AdS 2 horizon, and the effective action could be obtained from both SYK model and gravity theory on the boundary [12,13]. It consists of Majorana fermions with Gauss distribution random at a time. A fermion only moves by entangling with another fermion [1,14]. The q-body interaction (q is even) is (i) q/2 J i1,...,iq χ i1 ...χ iq with < J 2 i1,...,iq >= J 2 (q − 1)!/N q−1 . This model has maximal chaos [15], and the Lyapunov time τ L = 1/λ L describes how long a many-body quantum system becomes chaotic. There is a upper bound on the Lyapunov exponent defined in out-of-time correlations in thermal quantum systems [16]. The diffusion constants from holography are investigated in previous works, such as [17,18,19]. Naturally, the transport and diffusivity properties in SYK model could be an interesting aspect. The Figure 1: The 1-dimensional model: Each unit cell consists not only two kinds of SYK fermions (creation f † , annihilation f ) on A-sublattice but also two kinds of free fermions (creation g † , annihilation g) on B-sublattice. fermion can be written in Majorana fermion As illustrated in Figure 1, our 1-dimensional model comprises of L unit cells, and there are two kinds of sublattices in each unit cell. The coupling J ijkl,x denotes the interaction among fermions on the same sub-lattice, and the coupling t iα,x denotes the interaction among fermions on different sub-lattices. The coupling t iα,x denotes interaction for adjacent cell.
In order to describe a grand-canonical ensemble, we replace the Hamiltonian by where x labels the lattice site, L is the length of the chain. Note that another choice for the interaction [54,55] for details). These Gauss distribution randoms are all anti-symmetric tensor with zero mean obeying We define the fermion number In the large N limit, we introduce the bilocal fields with O(N ) symmetry Different replica indices m, m do not interact.
The green function G f in our model is given as, and G g is analogous.
After integrating the fermions, the partition function can be written as a path integral with collective modes G f /g (τ 1 , τ 2 ) and Σ f /g (τ 1 , τ 2 ), with a collective action based on the Luttinger-Ward analysis in [56], Σ is the self energy which contains only irreducible graphs. We define r ≡ M/N , it changes from Wigner-Dyson distribution (r < r c ) to Possion distribution (r > r c ). r c denotes the critical value. The equation for the self energy in the grand canonical in the large N limit is In the IR limit, we drop the ∂ τ term in the effective action, and obtain the self-consistency equations (i.e. the saddle point equation) by taking δS col δG f /g = 0, δS col δΣ f /g = 0 in the large N limit The above Schwinger-Dyson equations also can be derived by summing up the one particle irreducible diagrams [30].

Fluctuations and diffusion constants
The reparametrization symmetry maintains the emergent conformal symmetry of the original SYK model, which is spontaneously broken to SL(2, R) leading to zero modes [12]. SL(2, R) is the isometry group of AdS 2 . In the complex SYK model, an additional U(1) phase field φ is needed [10]. The conserved U(1) density is related to the fermion number constraint [14] and the chemical potential [57]. So under reparametrization of time τ → h(τ ), we havẽ where a = f, g.
Note that the second part in (26) is Schwarzian, which corresponds to the action of black hole with AdS 2 × R 2 horizon. When r = 1, t = 0 , t = 0 and the adjacent cell has been decoupled. Since the isolated SYK islands are not connected, the diffusion constant D 2 vanishes. When r 1, the SYK fermions dominate such that a finite thermal diffusion constant exist. When r 1 the free fermions dominate such that the diffusion constant is imaginary. So there is a phase transition, and |r − 1| could be considered as the order parameter.
Since the phenomenological coupling K, γ could not be evaluated analytically, we begin the large q expansion expand by small ∆ = 1/q at low T with three universal thermodynamics quantities Q, S, E. The T → 0 limit of the entropy is given by When Q = 0, (30) return to Appendix C of [50]. To solve the inverse function 2πE = dS(Q) dQ , we make an ansatz Taking E independent, we obtain The grand potential is given by the express 4 Then, we follow the analysis in [10] and obtain where υ satisfy πυ cos(πυ/2) = J T , J 2 = q 2 J 2 2(2 + 2 cosh µ T )) q/2−1 . Considering the thermodynamic relation, we have the free energy, the chemical potential and the entropy The entropy at low temperature could be obtained as, Additionally, we collect the relationship between thermodynamic parameters and diffusion constants according to the effective action derivation directly in Appendix H of [10].
Combined with the O(∆ 2 ) result, we obtain When r → 0, there are almost entirely SYK fermions. As the previous analysis, the system is composed of isolated stacks of SYK Majorana fermions. Therefore, the diffusion constant D 1 vanish. When r → 1, SYK fermions are as many as free fermions. In this case, the diffusion constants also vanish, and the behavior ∝ (r c − r) 5/2 implies a phase transition at r = r c = 1.

Transport properties and the butterfly velocity
In order to investigate the diffusivity and conductivity properties, we characterize transport by two-point correlations of the conserved number density firstly [59] and obtain transport coefficient by Green-Kubo relation. Similar to the analysis in [10,60,61], we obtain the dynamic susceptibility matrix where χ s = lim k→0,ω→0 χ(k, ω), Q , E is the conserved charge and the energy density in per site lattice, which depend on wavevector k and frequency ω The conductivities matrix is where σ is the DC electric conductivity, α is the thermoelectric conductivity, κ = κ−(T α 2 )/σ is the thermal conductivity. Then, the Wiedemann-Franz ration could also be given as From the analysis above, we observe that D 2 is really the thermal diffusion constant, as D 2 is related to the thermal conductivity. In many cases [20,60,63], there is a simple relation between the butterfly velocity and the thermal diffusive and the previous work [50] show the relationship remains by calculating the growth of out-oftime-ordered four-point function(see Appendix for the relation between the butterfly velocity and the thermal diffusive). We continue to find the relation between the butterfly velocity and the charge diffusive which depends on the interaction J of our generalized SYK model. Apparently, the relationship between υ B and the charge diffusion constant, unlike the thermal diffusion constant, does not depend on temperature T. Furthermore, (44) also shows that the butterfly velocity υ B vanishes at the critical point r → 1, which indicates an MBL phase. Next, we study the conductivity properties at low temperature. By considering the leading and next-to-leading contribution in the large q expansion, we obtain the DC electric conductivity σ and the thermal conductivity κ in the small T expansion, The first term in the electrical conductivity (45) is constant, and the second term in (45) is proportional to 1/T . This result is close to the dependence of the electric resistivity ρ = 1/σ on the low temperature in the normal phase of high temperature superconductors [64]. The thermal conductivity κ (46) is linear in the temperature . This result coincides with the cuprate strange metal as in [65,66]. The approximate behavior in (45)(46) indicate not only a dynamical phase transition but also an interesting temperature-dependence. Besides, the heat capacity could be written as which shows that the heat capacity is linearly with the temperature. Because of the experiment on optimally doped YBCO gives c ∼ T above the critical temperature [67], our result shows a good agreement qualitatively in the normal phase of cuprates.

Conclusion and discussion
In this paper, the diffusive property of an extended SYK model with complex fermions is investigated in the large N limit. We derived the collective action and the Schwinger-Dyson equations. With our saddle point solutions in the IR limit, we studied the fluctuations in the effective action. It consists P SL(2, R) and U (1) symmetries. Then, we obtained the quantum transport when the temperature is taken to zero. With the diffusion constants in our model, we calculate the DC electric conductivity, the thermal conductivity,the heat capacity and the relation between the butterfly velocity and the diffusion constants. As the analysis shows in Section 4, our novel results are qualitatively similar to strange metal. As noted in (29) and (39), the thermal diffusion constant D 2 relates to the coupling J, and the charge diffusion constant D 1 relates to the temperature T . Moreover, we explore the relationship between υ B and D 1 , D 2 . We notice that the thermal diffusion constant in the complex model is similar to the one in the real model [50]. When M = 0 (i.e. r = 0) which means that there are no free fermions, the relevant components of the Hamiltonian L x=1 vanishes. Thus, it is natural that the effective action in our model could return to the Gaussian action for the zero-dimensional complex SYK model as presented in [10] consistently.
Besides, as an extension of the original zero-dimensional SYK model, we focus on the one-dimension case. The higher dimension case is straightly replaced with vectors x → x. For example, the vectors connecting neighboring unit in the two-dimensional model are (0, 1), (1, 0), (1, 1). However, there are still some subtle uncertainties in our calculations. First, our calculation is based on the large q expand. Therefore, a future analytical study would be interesting and natural. Second, our model is not applied at r = 0, 1. The issue also occurs in the real case as in [50], we leave the critical theory for future study.