Holographic fermionic system with dipole coupling on Q-lattice

We construct a holographic model for a fermionic system on Q-lattice and compute the spectral function in the presence of a dipole coupling. Both key features of doped Mott insulators, the dynamical generation of a gap and spectral weight transfer, are observed when adjusting the value of the coupling parameter $p$. Of particular interest is that when the background is in a deep insulating phase, the Mott gap opens much easier with a smaller coupling parameter in comparison with a metallic background. The effects of lattice parameters on the width of the gap $\Delta/\mu$ are studied and a turning point is observed near the critical regime of metal-insulator transitions of the background. Furthermore, the temperature dependence of the spectral function is studied. Finally, we also observe that the anisotropic Q-lattice generates anisotropic peaks with different magnitudes, indicating that insulating and metallic phases arise in different directions.


I. INTRODUCTION
To understand and describe Mott metal-insulator transition (MIT) is a long-standing and widely known difficult problem in condensed matter physics because it involves a strongly correlated electron system, in which the conventional theoretical tools prove of little help and non-perturbative techniques are called for. At this stage, holography may provide insights into the associated mechanisms of these strongly correlated electron systems by building a gravitational dual model which is usually solvable in the large N limit. Some excellent examples are the holographic superconductor [1][2][3] and holographic (non)-Fermi liquid [4][5][6][7][8][9][10][11][12][13][14]. In this paper, we shall study the Mott transition induced by a dipole coupling in holographic Q-lattice background.
As early as 1930s, it was reported that many transition-metal oxides (such as NiO) with partially filled bands show insulating behavior. Peierls ascribed that to the strong electron-electron correlation. And then, Mott made a series of pioneering work towards understanding how electron-electron correlations could explain the insulating state [15][16][17][18][19].
Roughly speaking, the main idea of Mott is that the transition from metal to insulator occurs as lattice constant increases. Subsequently, this idea of Mott is formalized in the Hubbard model in which the Mott transition depends on the competition between the kinetic energy (E k ∼ t) and the on-site Coulomb repulsion U . The metallic behavior prevails if the kinetic energy t overcomes the Coulomb energy U while insulating phase is favored for U/t 1 and a gap opens in the single-particle excitation spectrum, resulting in a Mott transition at a critical ratio of U/t.
By adding probe fermions with dipole coupling in RN-AdS black hole [20,21], a Mott gap opens dynamically, which exhibits two key features of doped Mott insulator, i.e., the dynamical generation of a gap and spectral weight transfer. And then, the dipole coupling effects have also been studied in more general geometries in [22][23][24][25][26][27][28][29]. Along this line, here we shall study the holographic fermionic system with dipole coupling in Q-lattice geometry.
Holographic Q-lattice model is firstly proposed in [30], which is similar to the construction of Q-balls [31]. Some extensive studies have been presented in [32][33][34][35]. In this framework, the translational symmetry is broken and MIT is observed through the study of optical conductivity. Different from the holographic scalar lattice and ionic lattice constructed in [36][37][38][39][40][41], which involve solving the PDEs and need a hard numerical work, but here we only need solve ODEs for Q-lattice. In particular, MIT is not observed in scalar lattice or ionic lattice background yet. The difficulties of numerical calculation prevent one from dropping the temperature down to an extremely low level in this context 1 .
This paper is organized as follows. In Section II, based on the holographic model originally presented in [30], which introduces Q-lattice structure in one spatial direction, we generalize it to a two dimensional lattice model which in general can be anisotropic in x − y plane.
And then, the Dirac equations over this Q-lattice background are deduced in Section III.
In Section IV, we present our numerical results for the fermionic spectral function. We conclude in Section V with a summary of main results and suggestions for future research.

II. HOLOGRAPHIC Q-LATTICE GEOMETRY
Here we are interested in a holographic Q-lattice geometry with no translational symmetry in both of the spatial directions. For this purpose, we consider a system containing two complex scalar fields φ 1 and φ 2 plus a Maxwell field A µ . A similar construction can also be seen in [32]. The simplest action may take the form as, where F = dA. The equations of motion can be derived from the above action as follows Consider the following ansatz 1 MIT is also observed in helical lattice model [42]. In addition, a simpler construction which breaks the translational invariance can be found in [43,44], in which the lattice amplitude is absent, and thus no MIT happens in these models. with and substitute them into (2), one has five second order ODEs for V 1 , V 2 , a, ϕ 1 , ϕ 2 and one first order ODE for Q. Note that in above ansatz, k 1 and k 2 are two wave-numbers along x and y directions, respectively, such that φ 1 is periodic in x direction with a lattice constant 2π/k 1 and φ 2 is periodic in y direction with a lattice constant 2π/k 2 . In addition, in our holographic setup, the dual CFT involves two complex scalar operators with scaling dimension ∆ 1,2 = 3/2 ± (9/4 + m 2 1,2 ) 1/2 . When either of the masses falls in the interval [−9/4, −3/2), the AdS 2 BF bound near the horizon will be violated by the scalar field, possibly leading to a different ground state. To avoid this possibility, we will set m 2 1,2 = −5/4 for definiteness in our paper, which corresponds to the scaling dimension of ∆ 1,2+ = 5/2 and ∆ 1,2− = 1/2.
To solve the ODEs numerically, we impose a regular boundary condition at the horizon z = 1 and impose the following conditions on the conformal boundary In this paper, we only focus on the standard quantisation of the scalar field, in which the asymptotic behaviors of ϕ 1 and ϕ 2 look like ϕ 1 = λ 1 z ∆ 1− and ϕ 2 = λ 2 z ∆ 2− , respectively.
The UV behavior of the scalar field corresponds to a Q-lattice deformation with lattice amplitudes λ 1,2 . In addition, the temperature of the black hole is given by As a result, each of our Q-lattice solutions is specified by five dimensionless quantities T /µ, λ 1 /µ ∆ 1− , λ 2 /µ ∆ 2− , k 1 /µ and k 2 /µ. For simpleness, we shall still use T , λ 1,2 and k 1,2 to denote the above five dimensionless quantities in what follows.

III. DIRAC EQUATION
To explore the properties of fermionic spectral function in the Q-lattice geometry, we consider the following action including the dipole coupling with strength p between the fermion and gauge field [20,21,45,46] where D a = ∂ a + 1 4 (ω µν ) a Γ µν − iqA a and / F = 1 2 Γ µν (e µ ) a (e ν ) b F ab with (e µ ) a and (ω µν ) a being a set of orthogonal normal vector bases and the spin connection 1-forms, respectively.
Note that we have set the gauge coupling constant g F = √ 2 here. It is different from the conventions in [20,21], in which they set g F = 1. Consequently the charge q and dipole coupling p here will correspond to q/ √ 2 and p/ √ 2 in [20,21] as the relevant quantities are the products g F q and g F p.
The Dirac equation can be deduced from the above action To cancel off the spin connection, we can make a redefinition of ζ = (g tt g xx g yy ) − 1 4 F. At the same time, by the Fourier expansion, where k = (−ω, k x , k y ), one has Choose the following gamma matrices and split the 4-component spinor into two 2-component spinors as Furthermore, by the decomposition F α ≡ (A α , B α ) T with α = 1, 2, the above Dirac equation can be expressed as At the horizon, we can find that In order to obtain the retarded Green function on the boundary by holography, the independent ingoing boundary condition should be imposed at the horizon, i.e., Near the AdS boundary, the Dirac field reduces to And so by holography, the retarded Green function can be read off Note that since the four components of the Dirac fields couple to one another, we need to construct a basis of finite solutions, (A I α , B I α ) and (A II α , B II α ), to obtain the boundary Green function. We are mainly interested in the measurable spectral function, which is

IV. MOTT TRANSITION ON Q-LATTICE
In this section, we study the properties of spectral function on Q-lattice. We shall firstly address Mott transition in an isotropic Q-lattice, in which we set λ 1 = λ 2 and k 1 = k 2 . So, in this paper, we shall denote λ = λ 1 = λ 2 and k = k 1 = k 2 except in Subsection IV D, where we give a brief discussion on anisotropic spectral function. For definiteness, we work exclusively with the massless fermion and fix q = 1. In addition, we work at a very low but non-zero temperature of T 0.00398 for this paper except the subsection IV C, in which we explore the dynamics at different temperatures.

A. Free fermionic spectral function
In this subsection, we present the results for free fermionic spectral function on Q-lattice.
Because the notion of Fermi surface is only well defined at zero temperature, for our Q-lattice system of low but non-zero temperature, we need an operational definition, which has been proposed in [47] and adopted in holographic models [40,48]. It is argued that the Fermi momentum k F can be determined by searching the peak of A(ω, k) with a tiny frequency  coupling, the holographic Q-lattice background itself is not able to drive fermionic probes to undergo a Mott transition. To model Mott physics, we shall introduce the dipole coupling term as proposed in [20].

B. The Mott transition
A Mott insulator is signaled by a vanishing spectral weight, i.e., a suppression of A(ω) around ω = 0. The 3D plot of spectral function A(ω, k x ) with p = 4.5 for λ = 0.5 and k = 0.8 in FIG.3 shows such a gap, which has previously been observed in RN-AdS geometry [20,21]  and other backgrounds [22][23][24][25][26][27]. In addition, two bands are located in the positive frequency and negative frequency regions, respectively. Furthermore, we also show the spectral function To gain a better understanding of the spectral measure with dipole coupling on Q-lattice, we study the density of states (DOS) for our fermionic system. The DOS A(ω) signals the total weight of the spectral measure and is defined by the integral of A(ω, k x , k y ) over k x and k y . Since we have set k y = 0, the DOS only probes the effects in k x direction. is to open the gap dynamically. In numerical calculations, the onset of gap can be identified as that the DOS A(ω) at the Fermi level drops below some small number (here, we take 10 −3 ). With this in mind, we find that the critical value p c 3.98 for λ = 0.5 and k = 0.8, which corresponds to a metallic phase in [30]. While for an insulating phase with λ = 2 and k = 1/2 3/2 , p c = 3.74. It implies that the gap opens easier in a deep insulating phase than metallic phase. This observation can be further confirmed by evaluating the critical value for various wave number k and lattice amplitude λ, as we list in Table I and II. Nevertheless, we notice that an anomalous behavior occurs for k ≥ k c ( where k c is the critical value for MIT at zero temperature, with k c 0.8 for λ = 1), in which the critical value p c diminishes slightly with the increase of k, as opposite to that for k < k c .  Subsequently, we shall study the features of Mott insulating phase over a Q-lattice background. A characteristic measure is the width of the gap ∆/µ. FIG.6 shows the width ∆/µ as a function of the dipole coupling p for λ = 0.5 and k = 0.8. As p increases, the width of the gap becomes larger. It is just the effects of on-site interaction strength U of Hubbard model, which p plays.
We also study the effects on the width of the gap of wave-number k and lattice amplitude λ, which can be summarized as follows.
• First of all, a turning point is observed in the plot of gap width ∆/µ versus k,which is shown in FIG.7. Specifically, as k increases from small one, the width of the gap ∆/µ will decrease at first, which coincides with the idea of Mott that the strength of electron-electron correlations increases with the decrease of lattice constant [15][16][17][18][19].
However, when the wave-number exceeds a critical value k c , ∆/µ rises up slightly.
Interestingly enough, we find the values of k c falls into the critical regime of MIT in the phase diagram of the background, which is quite general. We expect to understand this anomaly with an analytical treatment in near future.
• The critical value k c is almost independent of the dipole coupling p but becomes larger with the increase of the lattice amplitude, as shown in FIG.7 and Table III. In addition, it is interesting to notice that in the right plot of Fig.7 all the lines almost intersect at one point of k 0.38.
• The width of gap as a function of the lattice amplitude λ is plotted in Fig.8. For k < 0.38, with the augmentation of the lattice amplitude λ, the width of the gap ∆/µ rises up at first and then tends to a constant or very slightly decreases, depending on the dipole coupling p. But for k > 0.38 the cases are just opposite.
• As a consistent check, one may find that in the right plot in FIG.8, all   . puzzles with V O 2 . It indicates that the driven force to form the gap in Mott transition is strongly correlated. By holography, the authors in [21] show the ratio ∆/T * in RN-AdS black hole is approximately 10 for p = 6 (or p = 7), which suggests holography can provide a good description on strong correlation induced transition. Here we shall address this issue on Q-lattice by holography.  Before working out the results of the ratio ∆/T * , we first make a qualitative description.
FIG .9 shows a 3D plot of the spectral function at T 0.115 and the DOS for different T with λ = 1 and k = 0.3. It is obvious that the gap closes when the temperature exceeds some certain value. Now, we work out the ratio ∆/T * for different lattice amplitudes λ and wave-numbers k for p = 6 (  system with dipole coupling on Q-lattice, the temperature dynamics possesses non-trivial behavior, as revealed in RN-AdS black hole [21].
In this subsection, we briefly discuss the Fermi surface and gap on anisotropic Q-lattice.
As shown in FIG.10, if we set q = 1, λ 1 = 2, λ 2 = 0.1 and k 1 = k 2 = 0.8 for p = 0, we find that the Fermi peaks along k x direction develop into some bumps. It exhibits an anisotropic peaks with different magnitudes, which suggests that insulating phase arises in one direction while metallic in the other. Furthermore, when we switch on the dipole coupling p, the gap along k x direction gradually opens but the bump still remain along k y direction. Therefore, the anisotropic geometry with insulating phase in one direction but metallic in the other provides us more space to model real materials with anisotropy by holography.

V. CONCLUSION AND DISCUSSION
Our main results in this paper are: • Two key features of doped Mott physics, the dynamical emergence of a gap and the spectral weight transfer, are observed in the Q-lattice background with dipole coupling.
Such features have been observed in many of holographic fermionic systems including dipole coupling, e.g., RN-AdS geometry and other geometry. It confirms the robustness of the generations of Mott gap induced by dipole coupling, which play the double roles of doping as well as the interaction strength U in the Hubbard model.
• The most important thing is that the fermionic system with dipole coupling on Qlattice can exhibit abundant Mott physics due to the introduction of wave-number k and lattice amplitude λ.
• The evolution of the spectral function as a function of temperature reveals the nontrivial temperature dynamical behavior, which indicates the Mott transition induced by dipole coupling on Q-lattice is due to strong correlations but not the spontaneous breaking of some symmetry.
• For the free fermionic system, the wave-number k and the lattice amplitude λ can not generate the Mott gap. When dipole coupling p exceeds some critical value, the gap opens and, usually, the gap opens much easily in deep insulating phase.
• The anisotropic peaks with different magnitudes occurs in our holographic fermionic system on Q-lattice. It indicates that insulating and metallic phases arise in different directions, respectively.
Still, a number of problems deserve further exploration.
• Rather than the dipole coupling, we may introduce other sorts of couplings to induce the Mott transition. Especially the interaction of electron-phonon asζ(η 1 φ 1 +η 2 φ 2 )ζ + h.c may have a manifest periodic structure and generate Brillouin zones. Our work on this subject is under progress.
• By the poles and zeros duality through the detG R , the pseudo-gap phase can be observed in the holographic fermions with dipole coupling in RN-AdS geometry [49].
It is valuable to address this problem in the holographic Q-lattice model.
• It is interesting to study the Fermi arc phenomenology in the anisotropic Q-lattice model.