Supersymmetric SYK model and random matrix theory
Abstract
In this paper, we investigate the effect of supersymmetry on the symmetry classification of random matrix theory ensembles. We mainly consider the random matrix behaviors in the supersymmetric generalization of SachdevYeKitaev (SYK) model, a toy model for twodimensional quantum black hole with supersymmetric constraint. Some analytical arguments and numerical results are given to show that the statistics of the supersymmetric SYK model could be interpreted as random matrix theory ensembles, with a different eightfold classification from the original SYK model and some new features. The timedependent evolution of the spectral form factor is also investigated, where predictions from random matrix theory are governing the late time behavior of the chaotic hamiltonian with supersymmetry.
1400 R Street, Lincoln, Nebraska 68588, USA Walter Burke Institute for Theoretical Physics, California Institute of Technology,
1200 east California Boulevard, Pasadena, California 91125, USA Department of Physics, Boston University,
590 Commonwealth Avenue, Boston, MA 02215, USA Perimeter Institute for Theoretical Physics,
31 Caroline Street North, Waterloo, Ontario N2L 2Y5, Canada
1 Introduction
Physical systems with some stochastic or chaotic properties have some randomness in the setup of the fundamental hamiltonian, which could be effectively simulated in the context of random matrix theory. When choosing an ensemble from random matrix theory for a chaotic hamiltonian, we often need to consider the symmetries in the dynamics of the related physical system. The choice of standard matrix ensembles from symmetries, historically comes from the invention of Dyson testbook (), which is called threefold way when classifying Gaussian Unitary Ensemble (GUE), the Gaussian Orthogonal Ensemble (GOE), and the Gaussian Symplectic Ensemble (GSE). For more general symmetry discussion of interaction systems, the AltlandZirnbauer theory gives a more complete description as a tenfold classification Zirnbauer1996 (); AZ1997 (). In the practical usage, one of the most celebrated works would be the classification of interaction inside topological insulators and topological phases in a tenfold way ludwig (); ki ().
In the recent study, the rising interests of studies on SachdevYeKitaev (SYK) model gives another profound application in the random matrix theory classification. SYK model kitaev (); Sachdev:1992fk () is a
microscopic quantum hamiltonian with random Gaussian nonlocal
couplings among majonara fermions. As is maximally chaotic and nearly
conformal, this model could be treated as a holographic dual of
quantum black hole with horizon through the (near)
AdS/CFT correspondence
Almheiri:2014cka (); Cvetic:2016eiv (); Fu:2016yrv (); Polchinski:2016xgd (); Jevicki:2016bwu (); Maldacena:2016hyu (); Jensen:2016pah (); Jevicki:2016ito (); Bagrets:2016cdf (); Maldacena:2016upp (). In
the recent research people have also discussed several generalizations
of the SYK model
Gu:2016oyy (); Gross:2016kjj (); Berkooz:2016cvq (); Fu:2016vas (), such as
higher dimensional generalizations and supersymmetric
constraints. Some other related issues and similar models are
discussed in
old (); Sachdev2 (); Sannomiya:2016mnj (); Sachdev:2010um (); GarciaAlvarez:2016wem (); Hayden:2007cs (); Anninos:2013nra (); Sachdev:2015efa (); Perlmutter:2016pkf (); Anninos:2016szt (); Danshita:2016xbo (); Roberts:2016hpo (); Betzios:2016yaq (); Witten:2016iux (); Patel:2016wdy (); Klebanov:2016xxf (); Blake:2016jnn (); Nishinaka:2016nxg (); Davison:2016ngz (); Anninos:2016klf (); Liu:2016rdi (); Magan:2016ehs (); Peng:2016mxj (); Krishnan:2016bvg (); Turiaci:2017zwd (); Ferrari:2017ryl (); GarciaGarcia:2017pzl (); Bi:2017yvx (); Ho:2017nyc (). In the recent discussions, people have discovered that the SYK
hamiltonian has a clear correspondence with the categories of the three fold standard
Dyson ensembles, unitary, orthogonal and sympletic ensembles, in the random matrix theory You:2016ldz (); GarciaGarcia:2016mno (); Dyer:2016pou (); Cotler:2016fpe (). In the recent work, Dyer:2016pou (); Cotler:2016fpe (), it is understood that the timedependent quantum dynamics of the temperaturedependent spectral form factor, namely, the combinations of partition functions with a special analytic continuation in SYK model, is computable in the late time by form factors in the random matrix theory with the same analytic continuation, as a probe of the discrete nature of the energy spectrum in a quantum black hole, and also a solid confirmation on the threefold classification Cotler:2016fpe ().
In the route towards Dyson’s classification, one only considers the set of simple unitary or antiunitary operators as symmetries when commuting or anticomuting with the hamiltonian. An interesting question would be, what is the influence of supersymmetry, the symmetry between fermions and bosons in the spectrum, in the classification of symmetry class?
As is illuminated by research in the past, supersymmetry Sohnius:1985qm () has several crucial influences in the study of disorder system and statistical physics Efetov:1997fw (), and could be emergent from condensed matter theory models Lee:2010fy (). Originating from particle physics, supersymmetry will enlarge the global symmetry group in the theory, has fruitful algebras and strong mathematical power used in several models in quantum mechanics and quantum field theory, and is extremely useful to simplify and clarify classical or quantized theories. In the recent study of SYK model, the supersymmetric generalization for the original SYK has been discussed in detail in Fu:2016vas (), which shows several different behaviors through supersymmetric extensions. This model might give some implications in the quantum gravity structure of black hole in two dimension in a supersymmetric theory, and also a related conjecture in Cotler:2016fpe () for spectral form factor and correlation functions in super YangMills theory.
In order to explore the supersymmetric constraints on the random matrix theory classification, in this paper we will study the symmetry classification and random matrix behavior of the supersymmetric extension of SYK model by FuGaiottoMaldacenaSachdev’s prescription Fu:2016vas (). The effect of supersymmetry in the symmetry classification could be summarized in the following aspects,

Supersymmetry will cause the hamiltonian to show a quadratic expression. Namely, we could write as the square of . This condition will greatly change the distribution of the eigenvalues. From random matrix language statistical (), if is a Gaussian random matrix, then should be in a WishartLaguerre random matrix, with the eigenvalue distribution changing from Wigner’s semicircle to the MarchenkoPastur distribution. In another sense, the quadratic structure will fold the eigenvalues of and cause a positivity condition for all eigenvalues. Namely, if has the eigenvalue distribution that eigenvalues come in pair with positive and negative signs, the squaring will cause larger degeneracies and a folded structure in eigenvalues of energy. More over, the coupling degree might be changed when considering instead of . For instance, in the extended SYK model, is a nonlocal three point coupling, which is not even. This will change the previous classification in the hamiltonian based on the representation of Clifford algebra from mathematical point of view.

We find the Witten index or Witten parity operator , which is wellknown as a criterion for supersymmetry breaking Sohnius:1985qm (); Witten:1981nf (); Cooper:1982dm (); Cooper:1994eh (), is crucial in classifying the symmetry class for supercharge . Some evidence of this point also could be found in some other models or setups. For instance, Witten parity is the Klein operator which separates the bosonic and fermionic sectors in the supersymmetric systems oikonomou:2014kba (); Oikonomou:2014jea (). MateosGuilarte:2016mxm () provides a more nontrivial example, where the odd parity operators are used to move states along a chain of different fermion sectors. Reversely, in some systems where one can define a graded algebra, Klein operator serves as a key factor in realizing supersymmetry, which is helpful in models of bosonization and higher spin theories, etc.Brink:1993sz (); Plyushchay:1994re (); Plyushchay:1996ry (); Bekaert:2005vh (). For example, Plyushchay:1996ry () constructs the bosonized Witten supersymmetric quantum mechanics by realizing the Klein operator as a parity operator. Plyushchay:1994re () realize a BoseFermi transformation with the help of the deformed Heisenberg algebra which involves a Klein operator. Another interesting application of Witten operator is VanHove:1982cc (), where the author argues that incorporating the Witten operator is crucial in some computation in supersymmetric systems with finite temperature. In the supersymmetric SYK model we are considering, Witten parity and the antiunitary operator together become a new antiunitary operator, which will significantly enlarge the set of symmetries in the hamiltonian, and change the eightfold story for supercharge and hamiltonian .
These aspects will be investigated in a clearer and more detailed way in the paper.
This paper will be organized as the following. In Section 2 we will review the model construction and thermodynamics of SYK model and its
supersymmetric extensions. In Section 3 we will discuss the random matrix classification for models, especially supersymmetric extensions of the SYK model. In Section 4 we will present our numerical confirmation for symmetry classifications from the exact diagonalization, including the computation of the density of states and spectral form factors. In Section 5, we will arrive at a conclusion and discuss the directions for future works. In the appendix, we will review some knowledge to make this paper selfcontained, including basics on AltlandZirnbauer theory and a calculation on the random matrix theory measure.
2 Introduction on models
In this paper, we will mostly focus on SYK models and their extensions. Thus before the main investigation, we will provide a simple introduction on the necessary knowledge of related models to be selfcontained.
2.1 The SYK model
In this part, we will simply review the SYK model mainly following Maldacena:2016hyu (). The SYK model is a microscopic model with some properties of quantum black hole. The hamiltonian^{1}^{1}1One could also generalize the SYK model to general point nonlocal interactions where are even numbers larger than four. The hamiltonian should be (9) where (10) Sometimes we will discuss the general in this paper but we will mainly focus on the case. is given by
(1) 
where are Majorana fermions and they are coupled by the four point random coupling with Gaussian distribution
(2) 
where and are positive constants, and . The large partition function is given by
(3) 
where is the total ground state energy proportional to and it is roughly Cotler:2016fpe (). is the ground state entropy contributed from one fermion, and one can estimate it theoretically Maldacena:2016hyu (),
(4) 
where is the Catalan number. is the specific heat, which could be computed by
(5) 
and is a positive constant. This contribution is from the Schwarzian, the quantum fluctuation near the saddle point of the effective action in the SYK model. The Schwarzian partition function is
(6) 
where the path integral is taken for all possible reparametrizations of the thermal circle in different equivalent classes of the symmetry. The Schwarzian corresponds to the broken reparametrization symmetry of the SYK model. One can compute the oneloop correction from the soft mode of the broken symmetry,
(7) 
As a result, one can consider the correction from the soft mode if we consider an external oneloop factor . The density of states could be also predicted by the contour integral of the partition function as
(8) 
2.2 supersymmetric extension
Following Fu:2016vas (), in the supersymmetric extension of SYK model, firstly we define the supercharge^{2}^{2}2For the generic positive integer we can also define the supersymmetric extension with nonlocal interaction of fermions. The supercharge should be (20) where (21) And will recover the case in the main text.
(11) 
for Majonara fermions . is a random tensor with the Gaussian distribution as the coupling,
(12) 
where is also a constant with mass dimension one. The square of the supercharge will give the hamiltonian of the model
(13) 
where
(14) 
where is the summation of all possible antisymmetric permutations. Besides the shifted constant , the distribution of is different from the original SYK model because it is not a free variable of Gaussian distribution, which changes the large behavior of this model. In the large limit, the model has an unbroken supersymmetry with a bosonic superpartner . The Lagrangian of this model is given by
(15) 
In this model, the Schwarzian is different from the original SYK model. We also have the expansion for the large partition function
(16) 
But the results of and are different (while the specific heat is the same for these two models). In the large limit, the supersymmetry is preserved, thus we have the ground state energy . The zero temperature entropy is given by
(17) 
Moreover, the oneloop correction from Schwarzian action is different. As a result of supersymmetry constraint, the oneloop factor is
(18) 
which predicts a different behavior for the density of states
(19) 
3 Random matrix classification
It is established that SYK model is classified by random matrix theory in that the random interacting SYK hamiltonian fall into one of the three standard Dyson ensembles in the eightfold way You:2016ldz (); GarciaGarcia:2016mno (); Dyer:2016pou (); Cotler:2016fpe (). It is natural to believe that the supersymmetric extension can also be described by random matrix theory. To sharpen the argument, we derive the exact correspondence between each SYK hamiltonian and some random matrix ensembles, in other words, the eightfold rule for supersymmetric case. A priori, the supersymmetric SYK hamiltonian should lead to a different random matrix theory description than the original case. Superficially, the original SYK theory and its supersymmetric cousin are different have two major differences, which have been also mentioned in the previous discussions.

The degeneracy of the two hamiltonian matrices are different. The degeneracy of supersymmetric SYK model is also investigated by Fu:2016vas (), which we derive again using some different discussion in Section 3.2.2. The degeneracy space is enlarged by supersymmetry. Generally, the energy level distribution of random matrices is sensitive to the degeneracy and is thus sensitive to the supersymmetric extension.

Another difference is the apparent positive semidefiniteness of the hamiltonian being the square of the supercharge. We will see later that the positive constraint leads to a new eigenvalue distribution different from those of Gaussian ensembles.
Symmetry analysis is crucial in classifying the random matrix statistics of hamiltonian matrices. You:2016ldz (); Cotler:2016fpe () argue that the particlehole symmetry operator determines the class of random matrix theory statistics. The random matrix classification dictionary is determined by the degeneracy and the special relations required by having the symmetry. The systematic method of random matrix classification is established as the AtlandZirnbauer theory Zirnbauer1996 (); AZ1997 (), reviewed in appendix A. The antiunitary operators play a central role in the classifications. The AtlandZirnbauer also applies to extended ensembles different from the three standard Dyson ensembles, which we find useful in classifying the supersymmetric SYK theory. In Section 3.1 we derive again the eightfold way classification of original SYK hamiltonian using AtlandZirnbauer theory and find unambiguously the matrix representations of hamiltonian in each mod eight sectors. We notice that the matrix representation of hamiltonian takes block diagonal form with each block being a random matrix from a certain ensemble. This block diagonal form is also found by You:2016ldz () in a different version.
Naively one would apply the same argument to the supersymmetric hamiltonian, since it also enjoys the particlehole symmetry. But this is not the full picture. First, one need to take into account of hamiltonian being the square of the supercharge and is thus not totally random. In Section 3.2.1 we argue that the supercharge has a random matrix description which falls into one of the extended ensembles. Using the AtlandZirnbauer theory on we obtain its matrix representation in block diagonal form and use it to determine the matrix representation of the hamiltonian in Section 3.2.2. Second, in order to obtain the correct classification one needs to consider the full set of symmetry operators. Apparently particlehole is not enough since supersymmetry enlarges the SYK degeneracy space. We argue that the Witten index operator, , is crucial in the symmetry analysis of any system with supersymmetry. Incorporating we obtain the full set of symmetry operators. Finally, the squaring operation, will change the properties of the random matrix theory distribution of supercharge , from Gaussian to WishartLaguerre. The quantum mechanics and statistics in supersymmetric SYK models, based on the main investigation in this paper, might be a nontrivial and compelling example of supersymmetric symmetry class.
3.1 Syk
Now we apply the AltlandZirnbauer classification theory (see appendix A for some necessary knowledges) to the original SYK model You:2016ldz (); GarciaGarcia:2016mno (); Dyer:2016pou (); Cotler:2016fpe (). This is accomplished by finding the symmetry of the theory (and has been already discussed in other works, see You:2016ldz (); Cotler:2016fpe ()). First, one can change the majonara fermion operators to creation annihilation operators and by
(22) 
where . The fermionic number operator divides the total Hilbert space with two different charge parities. One can define the particlehole operator
(23) 
where is the complex conjugate operator ( and are real). The operation of on fermionic operators is given by
(24) 
where
(25) 
From these commutation relation we can show that
(26) 
To compare with the AltlandZirnbauer classification, we need to know the square of and this is done by direct calculation
(27) 
Now we discover that can be treated as a operator and it completely determines the class of the hamiltonian. Before we list the result, it should be mentioned that the degeneracy of hamiltonian can be seen from the properties of :

:
The symmetry exchanges the parity sector of a state, so there is a twofold degeneracy. However, there is no further symmetries caused by in each block, Thus it is given as a combination of two GUEs, where two copies of GUEs are degenerated. 
:
The symmetry is a parityinvariant mapping and , so there is a twofold degeneracy. There is no further independent symmetries. From AltlandZirnbauer theory we know that in each parity block there is a GSE matrix. Also, where two copies of GSEs are independent. 
:
The symmetry is a parityinvariant mapping and . There is no further symmetries so the degeneracy is one. From AltlandZirnbauer theory we know that in each parity block there is a GOE matrix. Also, two copies of GOEs are independent.
We summarize these information in the following table as a summary of SYK model,
Deg.  RMT  Block  Type  Level stat.  

0  1  GOE  GOE  
2  2  GUE  GUE  
4  2  GSE  GSE  
6  2  GUE  GUE 
where the level statistics means some typical numerical evidence of random matrix, for instance, Wigner surmise, number variance, or statistics, etc. Although the SYK hamiltonian can be decomposed as two different parity sectors, we can treat them as standard Dyson random matrix as a whole because these two sectors are either independent or degenerated (The only subtleties will be investigating the level statistics when considering two independent sectors, where two mixed sectors will show a manybody localized phase statistics instead of a chaotic phase statistics, which has been discussed originally in You:2016ldz ().) In the following we will also numerically test the random matrix behavior, and based on the numerical testing range of we can summarize the following table for practical usage.
10  12  14  16  18  20  22  24  26  28  

Ensemble  GUE  GSE  GUE  GOE  GUE  GSE  GUE  GOE  GUE  GSE 
3.2 supersymmetric classification
Supersymmetry algebra is a graded algebra, where states and operators are subdivided into two distinct parity sectors. In such an algebra there may exist a Klein operator JunkerGeorg:2012 () which anticommutes with any operators with odd parity and commutes with any operators with even parity. The Klein operator of supersymmetry algebra is naturally the Witten index operator.
Witten index might plays a role in the symmetric structure and block decomposition in the supersymmetric quantum mechanics. A simple example is JunkerGeorg:2012 (), in supersymmetry algebra, Define be the Witten operator. The Witten operator has eigenvalue and separates the Hilbert space into two parity sectors
(28) 
We can also define projection operators . In the parity representation the operators take block diagonal form
(29) 
Because of and the complex supercharges are necessarily of the form
(30) 
which imply
(31) 
In the above equation, takes and its adjoint takes . The supersymmetric hamiltonian becomes diagonal in this representation
(32) 
In this construction, the Hilbert space is divided by Witten parity operator. The hamiltonian is shown to take the block diagonal positive semidefinite form without even referring to the explicit construction of the hamiltonian. It is remarkable that the above computation is very similar to our work from Section 3.2.1 to 3.2.2. Applications of this property can be found inoikonomou:2014kba (); Oikonomou:2014jea (). They describe a supersymmetric Quantum Mechanics system where fermions scatter off domain walls. The supercharges are defined as a differential operator and its adjoint. From (32) the number of ground states of each sector is simply the kernel of the differential operator and the Witten index is computed. A more nontrivial example is provided by MateosGuilarte:2016mxm (). In this work, the Hilbert space is divided into an fermions Fock space. Thus the Hamiltonian can be expressed as the direct sum of all fermion sectors. The ladder operators and are odd operators and move states between different sectors.
The argument can also work in reversive way. Hidden supersymmetry can be found in a bosonic system such as a Calogerolike model
Calogero:1969xj (), a system of one dimensional Harmonic oscillators with inverse square interactions and extensions. What makes supersymmetry manifest is the Klein operator. The model and its various extensions are studied in Brink:1993sz (); Plyushchay:1994re (); Plyushchay:1996ry (); Bekaert:2005vh (); Brink:1992xr (). A trivial simple Harmonic operator has algebra . The algebra describes a bosonic system. grading is realized by introducing an operator . The new operator anticommutes with and thus is a Klein operator. Based on the Klein operator one can construct the projection operators on both sectors and also the supercharge. In this way the simple harmonic oscillator is “promoted” to have supersymmetry. A generalization to simple hamornic oscillator is the deformed Heisenberg algebra, . The corresponding system is an supersymmetric extension of the 2body Calogero. The model is also used in considerably simplifying Calogero model.
These evidences strongly support the argument that supersymmetry will change the classification of symmetry class in quantum mechanical models. In the following work, we will show that supersymmetric SYK model symmetry class can be explicitly constructed and change the classification of random matrix theory ensembles.
3.2.1 Supercharge in Syk
In the supersymmetric model, it should be more convenient to consider the spectrum of instead of , because is the square of . Although is not a hamiltonian, since we only care about its matrix type, and the AltlandZirnbauer theory is purly mathematical, can be treated as a hamiltonian. Similiar to the original SYK model, we are concerned about the symmtry of the theory. We notice that the Witten index is
(33) 
which is the fermionic parity operator up to a sign . Witten index and particlehole symmetry have the following commutation relation:
(34) 
Now we define a new operator, . It has a compact form
(35) 
and are both antiunitary symmetries of , with commutation relations:
mod 8  

0  
2  
4  
6 
and squares
(36) 
Thus, in different values of , the two operators and behave different and replace the role in and in the AltlandZirnbauer theory. Now we can list the classification for the matrix ensemble of supersymmetric SYK model
mod 8  Cartan Label  Type  

0  1  BDI (chGOE)  
2  1  DIII (BdG)  
4  1  CII (chGSE)  
6  1  CI (BdG) 
One can also write down the block representation of . Notice that the basis of block decomposition is based on the eigenspaces of antiunitary operators, namely, it is decomposed based on the parity.
3.2.2 Hamiltonian in theory
Now we already obtain the random matrix type of the supercharge. Thus the structure of the square of could be considered case by case. Before that, we can notice one general property, that unlike the GOE or GSE group in SYK, in the supersymmetric model there is a supercharge contains odd number of Dirac fermions as a symmetry of , thus it always changes the parity. Thus the spectrum of is always decomposed to two degenerated blocks. Another general property is that the spectrum of is always positive because is Hermitian and . Thus the random matrix class of will be some classes up to positivity constraint.

: In this case is a BDI (chGOE) matrix. Thus we can write down the block decomposition as
(37) where is a real matrix. Thus the hamiltonian is obtained by
(38) Since and share the same eigenvalues ( thus flips the sign of eigenvalues of , but after squaring these two eigenvalues with opposite signatures become the same), and there is no internal structure in (in this case is a symmetry of , , but , thus cannot provide any further degeneracy), we obtain that has a twofold degeneracy. Moreover, because and are both real positivedefinite symmetric matrix without any further structure, it is nothing but the subset of GOE symmetry class with positivity condition. These two sectors will be exactly degenerated.

: In this case is a CII (chGSE) matrix. Thus we can write down the block decomposition as
(39) where is a quaternion Hermitian matrix. Thus after squaring we obtain
(40) Since and share the same eigenvalues, and each block has a natural twofold degeneracy by the property of quaternion (Physically it is because thus flips the sign of eigenvalues of , but after squaring these two eigenvalues with opposite signatures become the same. Also, in this case is a symmetry of , , and ), we get a fourfold degeneracy in the spectrum of . Because and are quaternion Hermitian matrices when is quaternion Hermitian^{3}^{3}3We say a matrix is a quaternion Hermitian matrix if and only if for some real in a basis, and is symmetric while is skewsymmetric. There is an equivalent definition that, defining thus is a quaternion Hermitian matrix if and only if and . Thus it is shown directly that if is quaternion Hermitian then and , thus is still a quaternion Hermitian matrix. , are both quaternion Hermitian positivedefinite matrix without any further structure. As a result, it is nothing but the subset of GSE symmetry class with positivity condition. These two sectors will be exactly degenerated.

: In this case is a DIII (BdG) matrix. Thus we can write down the block decomposition as
(41) where is a complex, skewsymmetric matrix. Thus after squaring we obtain
(42) Firstly let us take a look at the degeneracy. Since and share the same eigenvalues and each block has a natural twofold degeneracy because in skewsymmetric matrix the eigenvalues come in pair and after squaring pairs coincide (Physically it is because thus flips the sign of eigenvalues of , but after squaring these two eigenvalues with opposite signatures become the same. Also, in this case is a symmetry of , , and ), we obtain a fourfold spectrum of .
Now take the operator as a whole, from the previous discussion, we may note that it is quaternion Hermitian because it could be easily verified that and . Thus must be a quaternion Hermitian matrix (there is another way to see that, which is taking the block decomposition by another definition of quarternion Hermitian, squaring it and check the definition again). Moreover, has a twofold degenerated parity decomposition thus in each part it is also a quarternion Hermitian matrix. Because in the total matrix it is a subset of GSE symmetry class (with positivity constraint), in each degenerated parity sector it is also in a subset of positive definite GSE symmetry class (one can see this by applying the total measure in the two different, degenerated part). 
: In this case is a CI (BdG) matrix. Thus we can write down the block decomposition as
(43) where is a complex symmetric matrix. Thus after squaring we obtain
(44) Since and share the same eigenvalues ( thus flips the sign of eigenvalues of , but after squaring these two eigenvalues with opposite signatures become the same), and there is no internal structure in (in this case is a symmetry of , , but , thus cannot provide any further degeneracy), we obtain that has a twofold degeneracy.
Similar with the previous case, we can take the operator and as the whole matrices instead of blocks. For we notice that the transposing operations make the exchange of these two sectors. However, the symmetric matrix statement is basisdependent. Formally, similar with the quarternion Hermitian case, we can extend the definition of symmetric matrix by the following. Define(45) and we could see that a matrix is symmetric real (or symmetric Hermitian) if and only if and (where means the basis changing over two sectors). We can check easily that satisfies this condition, thus must satisfy. Thus we conclude that the total matrix in a subset of GOE symmetry class (with positivity constraint).
Although from symmetric point of view, the hamiltonian of model should be classified in the subsets of standard Dyson ensembles. But what the subset exactly is? In fact, the special structure of the squaring from to will change the distribution of the eigenvalues from Gaussian to WishartLaguerre Wishart (); statistical (); MIT () (Although there are some differences in the powers of terms in the eigenvalue distributions.) We will roughly called them as LOE/LUE/LSE, as has been used in the random matrix theory research. Some more details will be summarized in the appendix B.
However, the difference in the details of the distribution, beyond numerical tests of the distribution function of the one pointeigenvalues, will not be important in some physical tests, such as spectral form factors and level statistics (eg. Wigner surmise). The reason could be given as follows. From the supercharge point of view, because is in the AltlandZirnbauer distribution with nontrivial (see appendix B), the squaring operation will not change the level statistics such as Wigner surmise and spectral form factors (which could also be verified by numerics later). From the physical point, as is explained in You:2016ldz (), the details of distribution (even if not Gaussian), cannot change the universal properties of symmetries.
Finally, we can summarize these statements in the following classification table (the degeneracies have been already calculated in Fu:2016vas ()),
Deg.  RMT  Block  Type  Level stat.  

0  2  LOE  GOE  
2  4  LSE  GSE  
4  4  LSE  GSE  
6  2  LOE  GOE 
For our further practical computational usage, we may summarize the following table for different s in the supersymmetric SYK random matrix correspondence. As we show in the next section, for , these theoretical consideration perfectly fits the level statistics.
10  12  14  16  18  20  22  24  26  28  

RMT  LSE  LSE  LOE  LOE  LSE  LSE  LOE  LOE  LSE  LSE 
Universal Stat.  GSE  GSE  GOE  GOE  GSE  GSE  GOE  GOE  GSE  GSE 
4 Exact Diagonalization
In this part, we will present the main results from numerics to test the random matrix theory classification in the previous investigations. One can diagonalize the hamiltonian exactly with the representation of the Clifford algebra by the following. For operators acting on qubits, one can define
(46) 
where means standard Pauli matrices acting on the th qubit, tensor producting the identity matrix on the other parts, and . This construction is a representation of the Clifford algebra
(47) 
And one can exactly diagonalize the hamiltonian by replacing the majonara fermions with gamma matrices to find the energy eigenvalues. Thus, all quantities are computable by brute force in the energy eigenstate basis.
The main results of the following investigation would be the following. In the density of supercharge eigenstates and energy eigenstates in the supersymmetric SYK model, the behavior is quite different, but coincides with our estimations from the random matrix theory classification: the spectral density of supercharge shows clearly the information about extended ensembles from AltlandZirnbauer theory, and the spectral density of energy shows a clear MarchenkoPastur distribution from the statistics of WishartLaguerre. Moreover, because both and both belongs to the universal level statistical class for GOE, GUE and GSE, the numerics from Wigner surmise and spectral form factor will show directly these eightfold features.
4.1 Density of states
The plots for density of states in SYK model and its supersymmetric
extension are shown in Figure 1 for comparison. For each
realization of random hamiltonian, we compute all eigenvalues. After
collecting large number of samples one can plot the histograms for all
samples as the function . For density of states in SYK model,
in small tiny vibrations are contained, while in the large the
distribution will converge to a Gaussian distribution besides the
small tails. However, in the supersymmetric SYK model the energy eigenvalue structure is
totally different. All energy eigenvalues are larger than zero because
. Because of supersymmetry the lowest energy eigenvalues
will approach zero for large , and the figure will come to a
convergent distribution. The shape of this distribution matches the eigenvalue distribution of WishartLaguerre, which is
the MarchenkoPastur distribution lec () in the large limit. For the supercharge matrices, as becomes
larger the curve acquires a dip at zero, which is a clear feature for extended ensembles and could match the averaged
density of eigenvalues of random matrices in CI, DIII AZ1997 ()
and chiral Jackson:1996xt () ensembles at large .
For numerical details, we compute (40000 samples), (25600 samples),
(12800 samples), (6400 samples), (3200 samples),
(1600 samples), (800 samples), (400 samples),
(200 samples), and (100 samples). The results for
original SYK model perfectly match the density of states obtained in
previous works (eg. Maldacena:2016hyu (); Cotler:2016fpe ()).
4.2 Wigner surmise
There exists a practical way to test if random matrices from a theory are from some specific ensembles. For a random realization of the hamiltonian, we have a collection of energy eigenvalues . If we arrange them in ascending order , we define, to be the level spacing, and we compute the ratio for the nearest neighbourhood spacing as . For matrices from the standard Dyson ensemble, the distribution of level spacing ratio satisfies the WignerDyson statisticswignerSurmise ()) (which is called the Wigner surmise
(48) 
for GOE universal class, , ; for GUE universal class, , ; for GSE universal class, , (In fact, these are labels for the field of representation. See appendices for more details). Practically we often change to , and the new distribution after the transformation is . Standard Wigner surmises are shown in the Figure 2. You:2016ldz () has computed the nearestneighbor level spacing distribution of the SYK model, which perfectly matches the prediction from the eightfold classification.
What is the story for the supersymmetric SYK model? A numerical investigation shows a different correspondence for the eightfold classification, which is given by Figure LABEL:wignersusy. One can clearly see the new correspondence in the eightfold classification for supersymmetric SYK models, as has been predicted in the previous discussions.
Some comments should be given in this prediction. Firstly, one have
some subtleties in obtaining correct s. Considering there are two
different parities in the SYK hamiltonian (), each
group of parity should only appear once in the statistics of
. For in SYK, the particle hole operator
maps each sector to itself, thus if we take all the
distribution will be ruined, serving as a manybodylocalized
distribution (the Poisson distribution). For
in SYK, the particle hole operator maps even and odd parities to
each other, and one can take all possible s in the distribution
because all fermionic parity sectors are degenerated. Similar things
are observed for all even in the supersymmetric SYK model. As we
mentioned before, the reason is that the supercharge is a symmetry
of , which always changes the particle number because it is an
oddpoint coupling term. Moreover, the standard ensemble behavior is
only observed for , and for small enough s we have no
clear correspondence. Similar things happen for original SYK model,
where the correspondence works only for , because there is no
thermalization if is too small You:2016ldz (). However, the
threshold for obtaining a standard random matrix from
supersymmetric extension is much larger.
In Section 3.2.1, we argued that the supercharge operator in supersymmetric SYK theory are also random matrices in some extended ensembles AZ1997 (); Zirnbauer1996 (). We compute the level statistics of and compare it with the Wigner surmises of three standard Dyson ensembles in cases with different . The result is presented in Figure LABEL:fig:qstat. We see the level statistics of matrices match the same ensembles as the corresponding hamiltonian. This result confirms the relationship between ’s random matrix ensemble and that of the corresponding . That we do not see extended ensemble in the ’s level statistics because the level statistic does not see all the information in the ensembles.
4.3 Spectral form factors
Before presenting the numeric results of spectral form factors, we will review the discreteness of spectrum and the spectral form factor following Cotler:2016fpe (). For a quantum mechanical system, the partition function
(49) 
could be continued as
(50) 
The analytically continued partition function is an important quantity to understand a discrete energy spectrum. Typically, people will compute the time average to understand the late time behavior, but for , it vibrates near zero at late time and the time average should be zero. Thus, we often compute . For a discrete energy eigenvalue spectrum, we have
(51) 
It’s hard to say anything general directly for a general spectrum, but one can use the longterm average
(52) 
for large enough ( means the degeneracy). For a nondegenerated spectrum, it should have a simple formula
(53) 
However, for a continuous spectrum, the quantity has vanishing longterm average. Thus, the quantity should be an important criterion to detect the discreteness. In this paper, we will use a similar quantity, which is called the spectral form factor
(54) 
In the SYK model, these quantities will have similar predictions with the hamiltonian replaced by random matrix from some specific given Dyson ensembles. For example, for a given realization from a random matrix ensemble with large , we have the analytically continued partition function
(55) 
where
(56) 
The properties of spectral form factors given by random matrix theory, , have been studied in Cotler:2016fpe (). There are three specific periods in . In the first period, the spectral form factor will quickly decay to a minimal until dip time . Then after a short increasing (the ramp) towards a plataeu time , will arrive at a constant plataeu. This pattern is extremely similar with SYK model. Theoretically, in the early time (before ), should not obtained by because of different initial dependence on energy, while in the late time these two systems are conjectured to be coincide Cotler:2016fpe ().
With the data of energy eigenvalues one could compute the spectral form factors, which have been shown in Figure 5 for supersymmetric SYK model. We perform the calculation for three different functions , and with and several s. Clear patterns similar with random matrix theory predictions are shown in these numerical simulations. One could directly see the dip, ramp and plateau periods. For small s there exist some small vibrations in the early time, while for large this effect disappears. The function is strongly vibrating because we have only finite number of samples. One could believe that the infinite number of samples will cancel the noisy randomness of the curves.
A clear eightfold correspondence has been shown in the spectral form factor. Near the plateau time of one should expect roughly a smooth corner for GOEtype, a kink for GUEtype, and a sharp peak for GSEtype. Thus, we observe roughly the smooth corners for , while the sharp peaks for (although the peaks look not very clear because of finite sample size). For , as shown in Figure LABEL:wignersusy there is no clear random matrix correspondence because is too small, thus we only observe some vibrations near the plateau time.
We also perform a similar test on the supercharge , plotted in Figure 6.
In Section 4.2, we numerically tested the nearest neighbour level statistics of which matches perfectly the statistics of the corresponding . The spectral form factors of are slightly different from those of , yet they show exactly the same eight fold behavior.
4.4 Dip time, plateau time and plateau height
More quantitative data could be read off from the spectral form factors. In Figure 7, Figure 8 and Figure 9 we present our numerical results for dip time of , plateau time of , and plateau height of respectively. For numerical technics, we choose the linear fitting in the ramp period, and the plateau is fitted by a straight line parallel to the time axis. The dip time is read off as the averaged minimal point at the end of the dip period, and the error bar could be computed as the standard deviation.
It is claimed in Cotler:2016fpe () that polynomial and exponential fitting could be used to interpret the dip time as a function of with fixed temperature. We apply the same method to the supersymmetric extension. However, we find that in the supersymmetric extension, the fitting is much better if we fit the GOEtype group () and the GSEtype group () separately. On the other hand, although we cannot rule out the polynomial fitting, the fitting effect of exponential function is relatively better. On the exponential fittings with respect to different degeneracy groups, the coefficients before are roughly the same ( for ) while the overall constants are different. That indicates that the eightfold degeneracy class or random matrix class might influence the overall factors in the dip time exponential expressions.
One could also read off the plateau time and exponentially fit the data. Similar with dip time, we could also separately fit the plateau time with respect to two different random matrix classes, and one could find a difference in the overall factors of these two groups, while the coefficients before are the same. There is a nontrivial check here. Theoretically from random matrix theory one can predict that the plateau time scales like Cotler:2016fpe (). In the large limit, the entropy should be roughly the ground state entropy. Analytically, the entropy is predicted by . Now check the largest we take (), we can read off the entropy by (GSEtype), (GOEtype), or (two groups together), which perfectly matches our expectation.
For the plateau height, one can clearly see an eightfold structure. From the previous discussion we obtain that the plateau height should equals to times a contribution from the degeneracy, which is clearly shown in the figure. For (GOEtype), the degeneracy is two thus points should be on the lower line, while for (GSEtype), the degeneracy is four thus points should be on the upper line. These observations match the prediction from random matrix theories.
5 Conclusion and outlook
In this paper, we use analytic arguments and numerical evidence to explore the supersymmetric constraints on the random matrix theory symmetry class. We focus on the supersymmetric SYK model, a supersymmetric generalization of nonlocalcoupled majonara fermions with similar chaotic behavior for a two dimensional quantum black hole.
Use the direct classification from random matrix theory, we show that for supersymmetric SYK model has a different behavior for structure. These arguments might be made to be more general: supersymmetry could directly change the universal class of Hamiltonian (GOE/GUE/GSE) by classifying the symmetry class of supercharge, where combinations of Witten index and antiunitary operators will make some new antiunitaries; On the other hand, the quadratic structure of the Hamiltonian will change the original type of distribution from Gaussian to WishartLaguerre. These points may happen for generic supersymmetric statistical physics models.
We also use numerical method, exact diagonalization to confirm the random matrix theory classification on the Hamiltonian and the supercharge of the supersymmetric SYK model. It is clear that if we check the spectrum density, the supercharge shows a clear feature from onepoint function of extended random matrix theory ensembles, while the Hamiltonian shows a feature of quadratic semicircle (MarchenkoPastur). However, for level statistics (eg. Wigner surmise and spectral form factor), the universal class GSE/GOE could capture important physical features, and the new eightfold rule could be verified.
Several future directions could be investigated. Firstly, one may consider higher supersymmetry constraints on the SYK model, such as generalization. Many thermodynamical and field theory properties of higher supersymmetric SYK theory are nontrivial, and it might be interesting to connect these properties to random matrix theory. Moreover, to understand the spectral form factor with supersymmetric constraints, one could also try to study superconformal field theory partition functions at late time. Finally, introducing supersymmetries in the symmetry classification of phases in the condensed matter theory will bring more understanding at the boundary of condensed matter and high energy physics. We leave these interesting possibilities to future works.
Acknowledgments
We thank Xie Chen, Kevin Costello, Liam Fitzpatrick, Davide Gaiotto, Yingfei Gu, Nicholas HunterJones, Alexei Kitaev, Andreas Ludwig, Evgeny Mozgunov, Alexandre Streicher for valuable discussions. We thank Takuya Kanazawa for comments on the draft. JL is deeply grateful to Guy GurAri for communications on the symmetry of the original and supersymmetric SYK models. TL, JL, YX and YZ are supported by graduate student programs of University of Nebraska, Caltech, Boston University and Perimeter Institute.
Appendix A Review on AltlandZirnbauer theory
In this appendix we make a brief review the AltlandZirnbauer theory
(eg., see Zirnbauer1996 (); AZ1997 ()) that brings hamiltonians to ten
different random matrix classes. In a physical system, symmtries can
appear and they consist a group , then the space of physical states
is a projective representation of the symmetry group. A fundamental
question we can ask is, what is the most general type of hamiltonian
the system can have.
We may visit the simplest example to get some intuitions. The action
of an element of on the Hilbert space can be either unitary or
antiunitary, thus there is a homomorphism from group to
which labels unitarity of operators. Let be the
subgroup of unitary operators, then splitts into irreps of :
(57) 
where are irreps and are their multiplicities in . If there is no antiunitary operators then followed by Schur’s lemma, the most general Hamitonians are those belong to the set
(58) 
plus Hermicity. This is called Type A in the AltlandZirnbauer theory, without any antiunitary operators. The case with the presence of antiunitary operators is more complicated. Let be an antiunitary operator, then the conjugation by , i.e. , is an automorphism of , thus maps a component to another . A simple case is when , which is easy to see that the most general hamiltonian is of form Zirnbauer1996 (); AZ1997 ()
(59) 
where is an Hermitian operator in component and acts on component . Thus it’s also of Type A.
The Type A is the simplest structure without any further symmetries. However, if we consider , and consider more antiunitary operators, the situation is much more technical. It turns out that possible hamiltonians with specific symmetric structures can be classified into ten classes. Here we skip the detailed analysis and directly present the final results. These classes are classified by the following three different operators,

, antiunitary, commutes with hamiltonian, and

, antiunitary, anticommutes with hamiltonian, and

, unitary, anticommutes with hamiltonian, and
If two of these three operators exist, the third will be determined by the following identity,
(60) 
The properties of these three operators can classify the hamiltonian in the following ten classes,
Cartan label  Block  Type  

A (GUE)  
1  AI (GOE)  
AII (GSE)  
1  AIII (chGUE)  
C (BdG)  
1  D (BdG)  
1  1  1  BDI (chGOE)  
1  1  CI (BdG)  
1  1  DIII (BdG)  
1  CII (chGSE) 
where there are no values in some corresponding operators we mean that there is no such a symmetry in the system. We also present the block representation in this table, where blocks are classified by the eigenspace of antiunitary operators. The first three ensembles in this table are original Dyson ensembles, while other extended ensembles are their subsets with higher symmetries. These classifications are widely used in theoretical physics, for example, the symmetry classifications of topological insulators and topological phases ludwig (); ki ().
Appendix B Eigenvalue distribution
This appendix is a simple introduction on the random matrix theory eigenvalue distribution (for instance, see Wishart (); MIT ()), the measure in the eigenvalue basis. For WignerDyson ensemble, this is given by the formula
(61) 
where is the set of eigenvalues, is the Vandermont determinant defined by
(62) 
and is a normalization constant depending on and . For different ensembles, is defined as
RMT  

AI(GOE)  1 
A(GUE)  2 
AII(GSE)  4 
For the remaining ensembles, the eigenvalues occur in pairs (because the operator introduced in the last appendix anticommutes with ), and the eigenvalues probability distribution is given by
(63) 
where we only take the positive one from a pair of eigenvalues, and is defined also as the corresponding normalization constant. In the AltlandZirnbauer classification, constants and are set as (considering the real model of us, we have set the flavor number and the topological index in chiral ensembles)
RMT  

BDI(chGOE)  1  0 
AIII(chGUE)  2  1 
CII(chGSE)  4  3 
CI(BdG)  1  1 
D(BdG)  2  0 
C(BdG)  2  2 
DIII(BdG)  4  1 
We will also need the eigenvalue distribution of the hamiltonian which is the square of , so we can take the square distribution of 63, which will change Gaussian distribution to WishartLaguerre, which is
(64) 
here are nonnegative and is a new normalization constant which is one half of . We could also write
(65) 
where . The following table summarize the related index for supersymmetric SYK model
0  BDI (chGOE)  0  1  
2  DIII (BdG)  1  4  
4  CII (chGSE)  3  4  
6  CI (BdG)  1  1 
In , the index precisely matches Wishart matrix. Moreover, Although the result has dependence for , which does not precisely match Wishart matrix from Dyson Gaussian ensemble by index , we could also use the terminology LOE/LSE to refer the universal class from squaring of Gaussian matrix, similar with AltlandZirnbauer classification as a subset of Dyson, regardless multiple antiunitary symmetries. Thus, we call as LOE/LSE/LSE/LOE respectively,
References
 (1) F. J. Dyson, J. Math. Phys. 3, 140 (1962).
 (2) M. R. Zirnbauer, J. Math. Phys. 37, 4986 (1996) [mathph/9808012].
 (3) A. Altland and M. R. Zirnbauer, Phys. Rev. B 55, 1142 (1997).
 (4) S. Ryu, A. P. Schnyder, A. Furusaki, and A. W. W. Ludwig, New Journal of Physics 12, 065010 (2010), arXiv:0912.2157 [condmat.meshall].
 (5) L. Fidkowski and A. Kitaev, Phys. Rev. B 83, 075103 (2011), arXiv: 1008.4138 [condmat.strel].
 (6) A. Kitaev, Talks at KITP, April 7, 2015, May 27, 2015, and Feb. 12, 2015; A. Kitaev, Talk given at the Fundamental Physics Prize Symposium, Nov. 10, 2014.
 (7) S. Sachdev and J. Ye, Phys. Rev. Lett. 70, 3339 (1993) [condmat/9212030].
 (8) A. Almheiri and J. Polchinski, JHEP 1511, 014 (2015) [arXiv:1402.6334 [hepth]].
 (9) W. Fu and S. Sachdev, Phys. Rev. B 94, no. 3, 035135 (2016) [arXiv:1603.05246 [condmat.strel]].
 (10) J. Polchinski and V. Rosenhaus, JHEP 1604, 001 (2016) [arXiv:1601.06768 [hepth]].
 (11) A. Jevicki, K. Suzuki and J. Yoon, JHEP 1607, 007 (2016) [arXiv:1603.06246 [hepth]].
 (12) J. Maldacena and D. Stanford, Phys. Rev. D 94, no. 10, 106002 (2016) [arXiv:1604.07818 [hepth]].
 (13) K. Jensen, Phys. Rev. Lett. 117, no. 11, 111601 (2016) [arXiv:1605.06098 [hepth]].
 (14) M. Cvetic and I. Papadimitriou, JHEP 1612, 008 (2016) Erratum: [JHEP 1701, 120 (2017)] [arXiv:1608.07018 [hepth]].
 (15) A. Jevicki and K. Suzuki, JHEP 1611, 046 (2016) [arXiv:1608.07567 [hepth]].
 (16) D. Bagrets, A. Altland and A. Kamenev, Nucl. Phys. B 911, 191 (2016) [arXiv:1607.00694 [condmat.strel]].
 (17) J. Maldacena, D. Stanford and Z. Yang, PTEP 2016, no. 12, 12C104 (2016) [arXiv:1606.01857 [hepth]].
 (18) Y. Gu, X. L. Qi and D. Stanford, arXiv:1609.07832 [hepth].
 (19) D. J. Gross and V. Rosenhaus, arXiv:1610.01569 [hepth].
 (20) M. Berkooz, P. Narayan, M. Rozali and J. Simon, JHEP 1701, 138 (2017) [arXiv:1610.02422 [hepth]].
 (21) W. Fu, D. Gaiotto, J. Maldacena and S. Sachdev, Phys. Rev. D 95, no. 2, 026009 (2017) [arXiv:1610.08917 [hepth]].
 (22) O. Parcollet and A. Georges, Phys. Rev. B, 59, 5341 (1999) [condmat/9806119].
 (23) A. Georges, O. Parcollet, and S. Sachdev, Phys. Rev. Lett. 85, 840 (2000), condmat/9909239.
 (24) P. Hayden and J. Preskill, JHEP 0709, 120 (2007) [arXiv:0708.4025 [hepth]].
 (25) S. Sachdev, Phys. Rev. Lett. 105, 151602 (2010) [arXiv:1006.3794 [hepth]].
 (26) D. Anninos, T. Anous, P. de Lange and G. Konstantinidis, JHEP 1503, 066 (2015) [arXiv:1310.7929 [hepth]].
 (27) S. Sachdev, Phys. Rev. X 5 (2015) no.4, 041025 [arXiv:1506.05111 [hepth]].
 (28) E. Perlmutter, JHEP 1610, 069 (2016) [arXiv:1602.08272 [hepth]].
 (29) D. Anninos, T. Anous and F. Denef, JHEP 1612, 071 (2016) [arXiv:1603.00453 [hepth]].
 (30) I. Danshita, M. Hanada and M. Tezuka, arXiv:1606.02454 [condmat.quantgas].
 (31) P. Betzios, N. Gaddam and O. Papadoulaki, JHEP 1611, 131 (2016) [arXiv:1607.07885 [hepth]].
 (32) L. Garcï¿½ï¿½aAlvarez, I. L. Egusquiza, L. Lamata, A. del Campo, J. Sonner and E. Solano, arXiv:1607.08560 [quantph].
 (33) D. A. Roberts and B. Yoshida, JHEP 1704, 121 (2017) [arXiv:1610.04903 [quantph]].
 (34) E. Witten, arXiv:1610.09758 [hepth].
 (35) A. A. Patel and S. Sachdev, arXiv:1611.00003 [condmat.strel].
 (36) I. R. Klebanov and G. Tarnopolsky, arXiv:1611.08915 [hepth].
 (37) M. Blake and A. Donos, arXiv:1611.09380 [hepth].
 (38) T. Nishinaka and S. Terashima, arXiv:1611.10290 [hepth].
 (39) R. A. Davison, W. Fu, A. Georges, Y. Gu, K. Jensen and S. Sachdev, arXiv:1612.00849 [condmat.strel].
 (40) N. Sannomiya, H. Katsura and Y. Nakayama, arXiv:1612.02285 [condmat.strel].
 (41) D. Anninos and G. A. Silva, arXiv:1612.03795 [hepth].
 (42) C. Peng, M. Spradlin and A. Volovich, arXiv:1612.03851 [hepth].
 (43) Y. Liu, M. A. Nowak and I. Zahed, arXiv:1612.05233 [hepth].
 (44) C. Krishnan, S. Sanyal and P. N. B. Subramanian, arXiv:1612.06330 [hepth].
 (45) J. M. Magan, arXiv:1612.06765 [hepth].
 (46) G. Turiaci and H. Verlinde, arXiv:1701.00528 [hepth].
 (47) F. Ferrari, arXiv:1701.01171 [hepth].
 (48) A. M. Garcï¿½ï¿½aGarcï¿½ï¿½a and J. J. M. Verbaarschot, arXiv:1701.06593 [hepth].
 (49) Z. Bi, C. M. Jian, Y. Z. You, K. A. Pawlak and C. Xu, arXiv:1701.07081 [condmat.strel].
 (50) W. W. Ho and D. Radicevic, arXiv:1701.08777 [quantph].
 (51) Y. Z. You, A. W. W. Ludwig and C. Xu, arXiv:1602.06964 [condmat.strel].
 (52) A. M. GarciaGarcia and J. J. M. Verbaarschot, Phys. Rev. D 94, no. 12, 126010 (2016) [arXiv:1610.03816 [hepth]].
 (53) E. Dyer and G. GurAri, arXiv:1611.04592 [hepth].
 (54) J. S. Cotler et al., arXiv:1611.04650 [hepth].
 (55) M. F. Sohnius, Phys. Rept. 128, 39 (1985).