Holographic Floquet states II:
Floquet condensation of vector mesons in nonequilibrium phase diagram
Abstract
With the aim to reveal universal features of hadronic matter and correlated Dirac insulators in strong ACelectric fields, we study the supersymmetric QCD with a finite quark mass driven by a rotating electric field . The analysis is done in the holographically dual D3/D7 system in the corotating frame, effectively. The nonequilibrium phase diagram is determined from the threshold electric field at which the insulator phase breaks down to a conductive phase due to the AC version of the Schwinger mechanism. The external field induces a rotating current originating from vacuum polarization and dissipative current in the insulating and conductive phases respectively. Intriguing features are observed as the frequency approaches resonance with the meson excitation energy . There, the threshold minimizes and a condensate of vector mesons with oscillating current exists even in the zero driving field limit. This state, which we call Floquet condensate of vector mesons, is expected to be dynamically stable realizing a nonthermal fixed point that breaks time translational and reversal symmetries. Our finding has many similarities with exciton BEC discussed in solid state systems, where the semiconductor is to be replaced by materials hosting gapped Dirac electrons, e.g. 3D topological insulators or bismuth. Vector meson Floquet condensate may also have implications in the prethermalized dynamics in heavy ion collision experiments.
1]Shunichiro Kinoshita, 2]Keiju Murata 3]and Takashi Oka \affiliation[1]Department of Physics, Chuo University, Tokyo 1128551, Japan \affiliation[2]Keio University, 411 Hiyoshi, Yokohama 2238521, Japan \affiliation[3] MaxPlanckInstitut für Physik komplexer Systeme (MPIPKS), Nöthnitzer Straße 38, Dresden 01187, Germany \affiliation[3]MaxPlanckInstitut für Chemische Physik fester Stoffe (MPICPfS), Nöthnitzer Straße 40, Dresden 01187, Germany \emailAddkinoshita@phys.chuou.ac.jp \emailAddkeiju@physh.keio.ac.jp \emailAddoka@pks.mpg.de
1 Introduction and overview
Nonequilibrium dynamics after abrupt phase transitions in strongly coupled gauge theories is of fundamental interest. Several processes sequentially take place after the transition. Imagine that we initiate the transition by destroying the groundstate vacuum using a pulsed electric field. Being driven by the external field, the system evolves quantum coherently for some time. A prominent process here is the Schwinger effect of quarkantiquark pair production [1, 2, 3] which can trigger a vacuum breakdown. If the pulse is long enough, a nonequilibrium steady state may be realized as a result of a balance between the excitation and relaxation processes, e.g. pair annihilation. The interest in this paper is to study the phase transition that is realized within the nonequilibrium steady state in the presence of a driving field.
Our analysis is done by combining the framework of AdS/CFT correspondence [4, 5, 6] with the Floquet theory for periodically driven systems, subsequently to the previous work [7]. The Floquet theory is a framework to study steady state solutions that are time periodic. In this paper, we study a specific case where the state is realized in a rotating electric field
(1) 
with a field in the plane. This simplifies the analysis considerably since we can obtain a homogeneous and stationary solution by moving to the corotating frame (Sec. 2.2). Experimentally, a strong circularly polarized laser can approximately realize a similar electric field.^{1}^{1}1 Note that this “rotating electric field” is distinct from a “circularly polarized electromagnetic wave”, a vacuum solution of the Maxwell equation, since its field does not have a spatial modulation along direction such as .
Among the strongly coupled gauge theories that can be treated holographically, we study a minimal model that shows phase transitions. The model is the supersymmetric QCD (SQCD) consisting of a single flavor quark and gluons that mediates interaction with gauge group SU at large and at large ’t Hooft coupling . The holographic dual is the D3/D7 system in the probe limit: We study the classical motion of a single probe D7brane in the AdS spacetime [8]. In the single flavor case , the back reaction from the Dbrane to AdS is negligible. Physically, this means that the gluon degrees of freedom is large and act as a bosonic heat bath enabling us to have a stable nonequilibrium steady state. In this model we can consider quarks with mass . When the quarks are massive (), the quark and antiquark in this system are bound to each other by a potential that are linear in short range, and becoming Coulomb like in the long range limit [9]. In the vacuum, there are no asymptotically free quarks, and the excitation spectrum consists of a series of mesons (quarkantiquark bound states), as well as massless gluons. These mesonic excitations, such as coherent excitations of , , and so on, correspond to fluctuations for various fields on the Dbrane. A phase transition can be induced by external electric fields or finite temperature changing the quark and antiquark excitations from gapped to gapless (dissociated). Geometrically, this corresponds to a change of the Dbrane embedding in the bulk spacetime and to an emergence of an effective “horizon”; The D fluctuations can now dissipate into the horizon. The D3/D7 system and its phase transition under nonrotating electric fields, namely static (DC) electric fields, has been studied in Refs. [10, 11, 12, 14, 15, 17, 13, 16].
The electric field can make the vacuum unstable against a creation of quarkantiquark () pairs (Fig. 1(a)), which is known as the Schwinger effect [1] (see Refs. [18, 19] for reviews on the QFT treatment).^{2}^{2}2 The holographic analogues of the Schwinger mechanism have been studied recently, some focusing on the elementary tunneling process described by fundamental strings [20, 21, 22, 23, 24, 25, 26, 27], and others studying the phase transition of the vacuum via Dbrane dynamics [28, 29]. Relations with the generation of entanglement via the EPR pairs were also discussed [30]. The rotating electric fields can also trigger the Schwinger mechanism (Fig. 1(b)). In our previous work, we have studied the case when the quarks are massless [7]. In this case,  pairs can be created without any energetic cost and, thus, the system is always in a gapless conductive phase [10] (Fig. 2(b)). A holographic version of the Floquet Weyl semimetal that shows anomalous quantum Hall response, i.e., a current can be induced that flows perpendicular to the applied electric field, was studied in details. In this paper, we holographically study the nonequilibrium steady state of strongly interacting massive quarks in rotating electric fields. In the remaining of this section we summarize our findings.
1.1 Floquet condensation and nonequilibrium phase diagram
When quarks are massive, a singleparticle excitation of the quarks or antiquarks is forbidden and only mesons can be excited under sufficiently weak electric fields. It is the insulator phase of the D3/D7 system. As illustrated in Fig. 2(a), an electric field can induce “vacuum polarization” [1] by aligning the vacuum charge fluctuation. In the case of QCDs, vacuum polarization is nothing but a coherent excitation of vector mesons . In the insulator phase, by applying a rotating electric field, the polarization will follow it as . The insulator phase of D3/D7 is a perfect insulator and the DCconductivity is zero. However, in a time dependent field, a polarization current can exist given by the time derivative of the polarization
(2) 
This is perpendicular to the electric field, i.e., . Here, is introduced to generically parameterize the phase delay of the current relative to the electric field. In the insulator phase, since the current is always perpendicular to the electric field, Joule heating is absent and the system is dissipationless. A time dependent oscillating condensate such as Eq. (2) is generically found in periodically driven systems and is called the Floquet condensate [31]. A prominent example studied theoretically in condensed matter physics is the superconducting order in periodically driven lattice models [32, 33, 34]. Thus, in this article, we call the time periodic condensate expressed as Eq. (2) the electric field induced Floquet condensate of vector mesons. However, this state in the insulator phase is nothing special nor surprising; It is a state where the order parameter is induced by an external field. A similar example would be a paramagnet in an oscillating magnetic field where the magnetization will simply follow the field.
When the field strength exceeds a threshold critical value, an insulatortoconductor transition (i.e., dielectric breakdown) takes place. In the conductor phase, as illustrated in Fig. 2(b), gapless excitation of quarks and antiquarks become possible. Charge carriers are now separated and dissipative current that are not associated to polarization can flow allowing Joule heating to become finite. The critical field is known as the Schwinger limit ; Since it is now dependent on the frequency , we coin the breakdown process induced by ACelectric fields as the ACSchwinger effect (Fig. 1(b)). The nonequilibrium phase diagram in ACelectric fields is separated into regions with conductive and insulating states by the ACSchwinger limit. This is given in Fig. 3 (a) recovering the results from section 4. Starting from the DC case (), the critical Schwinger limit becomes smaller as the frequency is increased, and vanishes when resonance with meson modes occurs. Since SQCD has an infinite tower of meson mass spectrum, the threshold touches zero at each resonances, and thus the phase structure becomes “lobeshaped” in the plane.
Figure 3(b) is the nonequilibrium phase diagrams focused around the first resonance. The purple dots represent the insulating state and the green dots represent conductive states. Near resonance, we find an interesting state which we call a zero field Floquet condensate of vector mesons, or simply a vector meson Floquet condensate, represented by an orange line in the diagram. It is known that multiple phases may coexist near the Schwinger limit [13, 14], which resembles first order phase transitions in equilibrium. This degeneracy of states can be lifted by displaying the phase diagram as a three dimensional plot using the amplitude as the depth. At a critical frequency, which we define , the lower boundary of the conductive phase touches the zero field () plane. From there, on the zero field plane, a continuous set of solutions (orange line) stretches up to the meson mass frequency . The vector meson Floquet condensate, as illustrated in Fig. 2(c), is a subset of the insulator phase solutions, which supports a rotating current even when the external driving field is zero. The current breaks time translational symmetry as well as the time reversal symmetry (). This state exists not only for the lowest meson resonance but for all resonances in the vector meson mass spectrum. It is a state with a coherent oscillation of “vacuum polarization” [1].
The vector meson Floquet condensate is expected to be dynamically stable within the present theory since mesonic excitations cannot dissipate in the probe limit (). If the system is isolated from an energy bath, there is no relaxation channel within the model that changes this state, in other words, it is a fixed point of time evolution, i.e., a nonthermal fixed point [36]. We note the relation of the Floquet condensate discussed here with more exotic states known as “time crystals”. The ground state time crystal proposed by F. Wilczek is, in our language, a Floquet condensate that is thermodynamically stable existing as a lowest energy state [35]. Another class of time crystal is the “Floquet time crystal” [37, 38], a state in a periodically driven closed quantum system showing oscillations with a periodicity longer than the external driving force that was experimentally studied in artificial matters [39, 40]. The vector meson Floquet condensate is not included in neither of these categories since it is an excited state, and the periodicity is the same as the driving.
1.2 Possible experiments: Bose Einstein condensation of chiral excitons in gapped Dirac materials
The mesons in the current D3/D7 system shows many similarities with excitons discussed in solid state physics. Quarkantiquark pairs and the interaction mediated by gluons replace the electronhole pairs and Coulomb interaction there. The vector meson Floquet condensate is then the counterpart of an exciton Bose Einstein condensate (BEC). Exciton BEC is already experimentally realized in semiconductors [41]. However, in contrast to the nonrelativistic electrons in semiconductors, the quarks and antiquarks in the D3/D7 system are Dirac particles that have spin and angular momentum degrees of freedom. Thus, if we translate the states in our theory to the condensed matter language, it should be excitons in gapped Dirac materials in the strong coupling limit, i.e., BEC limit (Fig. 4). It is possible to find materials that host gapped 3D Dirac electrons: bulk Bismuth [42], and three dimensional topological insulators [45, 46] such as BiSe and its family materials are good candidates. In such materials, by applying circularly polarized laser, one can excite excitons with an orbital angular momentum. The chirality, i.e., the sign of the angular momentum, is selected by the helicity of the laser. The vector meson Floquet condensate in SQCD is translated to the Bose Einstein condensate (BEC) of chiral excitons in gapped Dirac materials. We note that related theoretical models which shows oscillating exciton condensation are currently being studied [43, 44].
The formation of Floquet states and its quasienergy spectrum were observed in the surface two dimensional Dirac band a time resolved angleresolved photoemission spectroscopy (ARPES) experiment done in a three dimensional topological insulator [47]. This was done with circularly polarized laser whose photon energy is lower than the bulk gap. In order to realize BEC of chiral excitons, a laser with photon energy close to the bulk gap should be used. If BEC is realized, our theory predicts that even after the pump laser pulse has disappeared a Floquet state will remain that is maintained by the bulk rotating current.
How can we efficiently create a chiral exciton BEC with a large rotating current? The phase diagram Fig. 3 gives us some hints. First, using a pulse field, it is possible to directly create the BEC if the peak strength is large enough. This process is illustrated in Fig. 5, which is a constant slice of the phase diagram Fig. 3. Imagine that the field is given by where the envelope changes slowly. In a pulse field, increases from zero up to a maximum and then decrease. If this process is slow enough, the state will adiabatically follow the insulator state, as indicated by red arrows, until the state no longer exists. Then, a abrupt jump will take place to the other solution, which is conductive for and insulating for . As the field weakens, the state will follow the solution down to the chiral exciton BEC. This process is similar a hysteresis loop of magnetization in ferromagnets.
Second, the chirping technique, i.e., time dependent shift of frequency may be useful in following the chiral BEC solution up to large rotating current. As Figure 3(b) suggests, the resonant frequency of the condensate (orange line) redshifts from to as the condensation density increases. We can adiabatically follow the orange line by using a chirped pulse in which the laser frequency slowly redshifts. The merit of this approach is that there is no threshold in the field strength and can be done with relatively weak fields. There is an upper limit in the magnitude of the rotating current in the insulator phase since a insulator to conductive transition takes place at . This transition is nothing but the exciton (meson) Mott transition [49, 50, 48, 51, 16] which happens when the exciton density exceeds a threshold at which the Coulomb attraction that binds the excitons becomes strongly screened and the excitons dissolve into plasma of electrons and holes.
Another important realization of a rotating electric field is heavy ion collision. Calculations suggest that strong electromagnetic fields can be created when heavy ions fly by each other [52, 53, 54, 55, 56]. The field that is created in such processes mimics a halfcycle rotating electric field. If the frequency of the pulse is resonant with the meson excitation energy, it is possible to produce vector meson Floquet condensates by such fields efficiently.
The rest of the paper is organized as follows. In Sec. 2, we will see that the equations of motion for the Dbrane can be reduced to ordinary differential equations effectively by moving to the corotating frame [7], and in Sec. 3 we will explain three types of their solutions with different brane embeddings. On the basis of solutions obtained numerically, we show a nonequilibrium phase diagram and physical quantities of the dual gauge theory near the phase boundary in Secs. 4 and 5. In Sec. 6, we will give an analytical derivation of some features in the phase diagram for the high frequency limit. In Sec. 7, we provide a simple understanding of the phase diagram using a toy model.
2 Formalism for a Dbrane in external rotating electric fields
2.1 Timedependent setup for Dbrane
We consider supersymmetric QCD as the field theory. Its gravity dual is described by the dynamics of a probe D7brane in the AdS spacetime. The metric of is written as
(3) 
where is the AdS radius, is the AdS boundary. Hereafter, we take the unit . We study the D7brane in this geometry. Dynamics of the D7brane is described by the DiracBornInfeld (DBI) action,
(4) 
where is the tension of the brane, is the induced metric and is the gauge field strength on the brane worldvolume. In the AdS spacetime, the directions are filled with the D7brane. We assume spherical symmetry in and translational symmetry in the space. We also set for simplicity, which means zero baryon number density in the boundary theory. We write the worldvolume gauge field as
(5) 
For later convenience, we introduce a complex expression for the gauge field,
(6) 
For the brane position, we can consistently assume^{3}^{3}3 Because of the rotational symmetry in the plane, the schematic form of the equation of motion for becomes , where are invariant under rotation. It follows that is a consistent solution.
(7) 
Although we take into account the full timedependence of the gauge fields and the brane position, we will see that we can eliminate the timedependence by a change of variables in next subsection. In terms of the brane embedding and the complex field , the DBI action is written as
(8) 
where and . From the equations of motion derived from the above action, we obtain an asymptotic solution near the AdS boundary () as
(9)  
(10) 
where is a constant with the dimension of length and the dot denotes a derivative. and are expansion coefficients, which depend on time in general.^{4}^{4}4 For the rotating electric field, we will see that can be timeindependent. Therefore, we have dropped the argument for and . They are related to the electric field , electric current , quark mass and quark condensate in the boundary theory by
(11) 
Note that we use the complex notation for the electric field and current: real and imaginary parts of and correspond to  and components of the electric field and current, such as and .
2.2 Corotating frame and reduction to an 1dimensional problem
We focus on the rotating external electric field that is described by
(12) 
where is a complex constant whose argument represents the initial direction of the electric field at . The electric field rotates anticlockwise with a frequency in the plane. From Eq. (9), the boundary condition for the gauge field at the AdS boundary becomes
(13) 
We define two real variables and as
(14) 
where and represent the amplitude and the phase for , respectively. For new variables and , the boundary condition becomes timeindependent as
(15) 
Note that, since at the AdS boundary describes the direction of the electric field at a certain time, one can fix it to an arbitrary value.
Using new variables and , the DBI action (8) can be rewritten as
(16) 
This action does not depend on explicitly. The boundary condition (15) is also timeindependent. Thus, we can consistently assume that the variables , and do not depend on : , and . By the reduction to 1D problem, the DBI action (16) is simply written as
(17) 
where . Obviously, this action is invariant under the constant shift of the phase variable : , where is an arbitrary real constant. As mentioned above, this arbitrariness of a constant phase rotation of corresponds to changing the initial direction of the electric field in the plane in the boundary theory. Thus, we have a conserved charge defined by
(18) 
As we will see later, the above conserved charge corresponds to the Joule heating in the boundary theory. This implies that should be satisfied to make the Joule heating nonnegative. The ansatz for the gauge field in Eq. (14) is similar to that of an ordinary circularly polarized electromagnetic wave in the Maxwell theory. In fact, is nothing but the condition for the “circularly polarized” gauge field to propagate inward along the direction from the AdS boundary on the Dbrane. The above conserved quantity is the energy flux of this “circularly polarized” gauge field in the bulk theory.
From Eq. (17), equations of motion for , and are given by^{5}^{5}5 For numerical calculations, and are not suitable variables when a solution approaches the origin of the complex plane: although they are convenient for the simple expressions. (Note that Eq. (20) is singular at because the phase becomes illdefined .) To avoid this “coordinate singularity”, we use Cartesianlike variables in our actual numerical calculations. We can easily cast the equations of motion or boundary conditions for into those for .
(19)  
(20)  
(21) 
It turns out that is a trivial solution of Eq. (20). In this case, the conserved quantity is zero and the gauge field becomes a standing wave without a net energy flux rather than a propagating wave. As we will see later, whether such solutions can be admitted physically depends on boundary conditions for the other fields and .
It is worth noting that we can easily recover the DC limit from the expression of the action (17). This is done by fixing and while taking the limit, where denotes a direction determined by the alignment at , and then we recover the expression for static electric fields [28]. This is because this system still holds a translational invariance with respect to instead of time translational invariance with respect to . Therefore, the rotating electric field we are studying now is a natural generalization of static electric fields. Since the vector represents a corotating frame with angular velocity , the above translational invariance implies that the system should be stationary in the corotating frame. This enables us to treat the system in a similar way to static systems and provides its phase structure sharply as we will see later.
2.3 Observables in the boundary quantum field theory
Once we solve the equations of motion, observables in the boundary theory can be read out from asymptotic behavior of the solution near the AdS boundary . From the asymptotic form of , we can read off quark mass and condensate as in Eq. (11). Multiplying to Eq. (9) and using the rotating electric field ansatz (12), we obtain the asymptotic form of as
(22) 
We have defined by
(23) 
where is a timeindependent constant representing the electric current at . In particular, we have from Eq. (22). By using the symmetry of selecting the initial field direction, i.e., shift of , we can fix to be a real and nonnegative value.
Substituting the above expression into Eq. (18), which is the conserved charge in the bulk, we can rewrite the conserved charge as
(24) 
This explicitly confirms that the conserved quantity in the bulk theory corresponds to the Joule heating in the boundary theory. Besides the expression (22), the complex electric current can be written as
(25) 
where the last term in the bracket is arising from a counterterm in holographic renormalization [10]. It is known that the electric current has an ambiguity corresponding to a finite counterterm if the electric field is timedependent, and in fact the ambiguity is reflected in a constant of (25). However, it is worth noting that the Joule heating (24) is unaffected by this ambiguity as is clear from (25). In general, an electric current in a dielectric medium can be decomposed into the active current contributing to the Joule heating via resistive dissipation and the reactive current (or polarization current) originating in polarizations. It turns out that a part of proportional to represents the active current, while the remnant is the reactive current.
There is a scaling symmetry in our setup:
(26) 
We will nondimensionalize physical quantities using the quark mass .
3 Black hole, Minkowski and critical embeddings
The equations of motion (19)(21) are singular at and where satisfies
(27) 
The singular points or give the inner boundary of the integration of equations of motion. In section 3.4, we will see that is the effective horizon with respect to the effective metric on the Dbrane. There are three possibilities of the brane solutions depending on the position of the inner boundary:

The brane reaches , the axis of wrapped by the Dbrane, and the effective horizon does not exist. (i.e. no satisfying Eq. (27) exists.)

The brane intersects with the effective horizon at before it reaches the axis.

The brane reaches the axis which coincides with the locus of the effective horizon, namely .
The brane embeddings in the cases 1, 2 and 3 are called “Minkowski embeddings”, “black hole embeddings” and “critical embeddings”, respectively [57, 58, 59].
Figure 6 shows typical profiles of the Dbrane for the black hole, Minkowski and critical embeddings in the plane. We fix the electric field and its angular velocity as and . Asymptotic position of the brane corresponds to the quark mass. (See Eq. (11).) The trajectory of the effective horizon emerging on the Dbrane for fixed is shown by the black curve in the bulk spacetime. (Note that there is no horizon in the bulk spacetime.) For static electric fields , it is known that the trajectory of the effective horizon is given by a circle: for a fixed [11, 12, 16]. In this figure, by the effect of nonzero , the trajectory expands and is slightly deformed from the round circle.
3.1 Minkowski embeddings
Here, we consider the Minkowski embeddings (case 1). In this case the conserved quantity of (18) leads to because the brane reaches . Furthermore, since should be conserved for , we have the trivial solution We will focus on and , and expand the solutions around . The asymptotic solutions become
(28)  
(29) 
where . When the brane does not intersect with the effective horizon, no Joule heating is generated in the boundary theory. As we mentioned, means the worldvolume gauge field becomes standing wave without any net energy flux.
3.2 Black hole embeddings
In the case of the blackhole embeddings (case 2) there exists satisfying (27). The effective horizon emerges at on the Dbrane. Equations (18) and (27) directly give us the conserved quantity as
(30) 
Thus, in contrast to the Minkowski embeddings becomes nontrivial and the energy flux appears on the Dbrane. In the boundary theory this means Joule heating occurs due to the current.
The boundary conditions for , and are given by the regularity at the effective horizon . Expanding the solutions around , we obtain
(31)  
where and is arbitrary. Explicit forms of the first derivatives , and are shown in appendix A.
3.3 Critical embeddings
Solutions of the phase boundary between the Minkowski embeddings and the blackhole embeddings are the critical embeddings (case 3). In this case the locus of the effective horizon coincides with the axis, so that (27) yields because of . We expand the solutions around as
(32)  
(33) 
Since , solution of is arbitrary constant.
3.4 Effective metric and horizon
The effective metric is defined by [60, 61, 62, 63, 64]
(34) 
For the current system, the effective metric is written as
(35) 
where
(36) 
We have also defined 1forms as
(37) 
Note that becomes zero at . This implies that is the event horizon with respective to the brane effective metric . In general, dynamics of the brane itself and various fields on the brane is governed by the effective (induced) metric. The effective horizon is a causal boundary for fluctuations on the brane, so that the fluctuations across the effective horizon can never go back and will dissipate.
The effective Hawking temperature is given by
(38) 
where explicit expressions of , and at the effective horizon are shown in appendix A.
4 Nonequilibrium SQCD phase diagram in rotating electric fields
Here we explain the structure of the nonequilibrium phase diagram that we obtain by solving the equation of motions (19)(21) for the probe D7brane. There are two phases. The insulator phase is given by the Minkowski embedding solutions. In this phase fluctuations on the Dbrane, correspond to mesonic excitations, can never dissipate, so that it is a gapped phase where stable mesons can exist. Since we have , the Joule heating is absent although the rotating current is finite. This means that this current is a reactive current with no resistive dissipation. The conductive phase is given by the black hole embeddings. Fluctuations will dissipate into the effective horizon emerging on the Dbrane. This phase is gapless where the mesons have been dissociated and there are no stable mesons. Now, becomes a nontrivial function and an active current with the Joule heating emerges. This phase emerges after the dielectric breakdown of the insulator phase. The Joule heating plays a role as an order parameter to distinguish the phase transition between insulator and conductive phases.
We should describe the phase diagram in a physical parameter space. The brane solutions are originally specified by the boundary conditions with three parameters , , and at the inner boundary; for the Minkowski embeddings we specify and () at the axis , and for the black hole embedding and at the effective horizon , where . Using the scaling symmetry (26), we can set one of the parameters to unity without loss of generality ^{6}^{6}6In our numerical calculation, we set .. This means that the phase diagram of our system can be described by a twoparameters family of solutions, specified by the frequency and strength of the rotating electric field.
Figure 3(a) shows a phase diagram in the parameter space . The phase boundary between the insulator and the conductive phases is given by the critical embeddings (solid curve) giving the AC extension of the Schwinger limit for dielectric breakdown [1, 2, 3], that is, the ACSchwinger limit. The region above the curve indicates the conductive phase while that below indicates the insulator phase, basically. (As will see, actually, multiple phases coexist in the vicinity of the critical curve.) We see a lobeshaped phase structure, i.e., the insulator phase is separated into infinite but discrete regions analogous to a lobe (leaf). We can see that “lobes” become thinner as the frequency becomes larger. Its height decreases following a powerlaw (Fig. 7), which can be analytically verified in Sec. 6. The red points on the horizontal axis are mass spectrum of mesons in SQCD: [9]. The critical embedding solution (solid curve) touches the horizontal axis at discrete points smaller than but close to the meson mass frequencies. We refer to this frequency as the critical frequency and denote it as . The numerical values of the critical frequencies and meson mass spectrum are given by
(39) 
for the first four modes. In Sec. 6, we will examine more deeply how the critical and meson mass frequencies appear by employing a perturbative approach in a highfrequency regime. We find that the ACSchwinger limit periodically becomes infinitesimal near, but not equal to, frequencies corresponding to the meson mass spectrum. This indicates that the rotating electric fields can resonantly excite the vector mesons close to the meson mass frequencies, so the dielectric breakdown can easily occurs even though the rotating electric field is too weak.
We have overviewed the global structure of the phase diagram obtained from the critical solution. Next, let us focus on the detailed structure of the phase diagram near the phase boundary. Discrete selfsimilarity is known to emerge near the phase boundary, so that multiple Minkowski and black hole embedding solutions with a spiral structure take place near the phase boundary resulting in a characteristic first order phase transition. Such behavior has been studied in various probe brane systems [57, 58, 59]. The selfsimilarity leads to oscillatory behavior of physical quantities as one approaches the transition. Figure 3(b) shows fine structure of the phase diagram around the lowest meson excitation energy . We take as the third axis to distinguish degenerated states for a same (). We plot Minkowski and black hole embedding solutions by purple and green dots that are obtained by the changing parameter set ()^{7}^{7}7 For the Minkowski embedding, we set parameters at inner boundary as , and , where and (, ). For the black hole embedding, we set , and , where and (, ). . We see that both solutions coexist for the same parameters () near the critical electric field. Such multiplicity of the solutions was observed in the case of the static electric fields (), so that the above feature is a natural AC extension of that around the Schwinger limit. It is especially notable that the multiplicity extends to the horizontal axis () within a range of the frequency: (see Fig. 3 also). Among the two states at , one is the thermal ground state, while the other is a nontrivial excited state. This nontrivial state in the absence of external electric fields can be interpreted as a nonlinear excitation of the vector mesons: Floquet condensate of the vector mesons. It is known that the vector mesons can stably exist in the insulator phase. In general, coherent excitations of these vector mesons appear as oscillations of the expectation values with frequencies corresponding to each masses. However, if the vector mesons are circularly polarized, its amplitude does not oscillate while its phase oscillates (the direction of the polarization is rotating). When the excitation is small enough to be described by linear perturbations, the frequency is given by the meson mass . As the amplitude becomes larger, nonlinearity of the excitation cannot be negligible and its frequency will shift from the meson mass frequency.
When multiple states coexist, in the long time limit, a state that has lesser energy will be selected as a thermodynamical stable state, while the other may exist as a transient metastable state. Within the present model, although Floquet condensate of the vector mesons is thermodynamically unstable, this metastable state is expected to be a dynamically stable regarding timeevolution and defines a “nonthermal fixed point”. This is because the mesonic excitations can be stable without dissipating in the insulator phase.
If we finetune the frequency of the rotating electric field , we can control the insulatortoconductor phase transition by a weak electric field. To see this explicitly, let us focus on in Fig. 3(b). At this point, we have . Above this point (, ), only the conductive phase exist. This indicates that, for this frequency , we can cause the phase transition just by at most. Note that, for DC electric field , we need for the phase transition. (See the vertical axis in Fig. 3(a).) This is 100 times larger than the case of the finetuned frequency .
5 Physical quantities near the phase boundary
In order to understand the meaning of the multiple solutions for given (), we focus on their physical quantities in the dual gauge theory. From the asymptotic brane solution near the AdS boundary, we can read off two vacuum expectation values in the dual gauge theory. One of them is the rotating electric current , which appears as the nonlinear response against the external rotating electric field. Note that is a complex number, whose amplitude and phase represent the magnitude and the initial direction of the electric current, respectively. The other is the quark condensate . As the subsequent quantity, we also consider the Joule heating playing the role of the order parameter of the phase transition between the insulator and conductive phases. The Joule heating becomes nonzero only in the conductive phase, while the electric current can exist in the both conductive and insulator phase. It is a contrast to the case of the static electric fields () in which the insulator phase has neither the electric current nor the Joule heating.
Figure 8 shows the electric currents as functions of the electric field for several fixed frequencies . Near the phase boundary, the electric current becomes multivalued. In the multivalued region, the most stable state is expected to be realized. If we start increasing the external field from , the system goes through an insulatortoconductor transition at the critical electric field. Then, there will be a jump in the electric current as well in the Joule heating at the phase transition. (There is also a jump in the Joule heating. See Fig. 10.) From this perspective, the phase transition from insulator to conductor is “first order”. Such phase transitions between nonequilibrium states with the electric current have been often observed in a similar kind of holographic systems described by D3/D7branes, especially around the critical embedding solutions.
Spiral structure of the solutions can be seen taking place around the critical embedding solution as shown in the inset of the left panel of Fig. 8. As we increase the frequency from to (left panel), the center of the spiral approaches the vertical axis and touches at the critical frequency . Between (right panel), the insulator phase continues to touch and reflect at the vertical axis . At reflecting points, the electric field is zero but the current is nonzero. This is the Floquet condensate of the vector mesons, which is the nonthermal fixed points solutions (existing between ). As explained in section 2.3, the definition of the electric current contains an ambiguity with a term proportional to the electric field. (The ambiguous parameter in Eq. (25) is set to .) However, only the current at is determined without this ambiguity. This means that the expectation value of the Floquet condensate obtained here is definite at .
Figure 9 shows the quark condensates as functions of the electric field for several fixed frequencies . Note that, for the definition of the quark condensate, there is no ambiguity as in the electric current. The basic behavior of the condensate is very similar to that of the electric current Fig. 8 and we find it quite remarkable that the quark condensate is also nonzero as well as the Floquet condensate of the vector mesons.
Compared with meson condensate realized by linearly polarized fields ( and ), the condensate realized by rotating electric field is distinct as follows. The rotating electric field breaks time translational and reversal symmetries and makes the system a time periodic steady state, i.e., a Floquet state. Indeed its vector quantities such as the directions of the electric current and the polarization oscillate, but its scalar quantities such as the condensation of a scalar meson, and the magnitude of the polarization associated with vector mesons, do not oscillate.
Figure 10 shows the Joule heating as functions of the electric field . This quantity is only nonzero in the conductive phase. We see that the Joule heating is also multivalued in terms of the electric field. It is remarkable that, within a small region of the frequency , the multiplicity of the Joule heating disappears and the nonequilibrium phase transition between conductive phases does not occur [14]. In general, the multiplicity around the critical solutions may be infinitely large, i.e., infinite number of solutions can coexist for the same value of the electric field. In the current numerical analysis, it is still unclear whether all of the infinite multiplicity can disappear around the critical frequency.
So far, we have considered physical quantities for fixed . In Fig. 11, we show the Joule heatings as functions of for several fixed . When we take the slice of in Fig. 3(a), the slice can cross several lobelike regions in the insulator phase. In the lobelike regions, the Joule heating is zero. Therefore, for small such as the curve for in the figure, the Joule heating becomes zero in several discrete regions. For fields above the critical field, e.g., , we still find a “dip” in the Joule heating around reflecting the lobeshaped phase structure.
6 Analytic derivation for features of the phase diagram in the highfrequency limit
When the frequency of the rotating electric fields is sufficiently high, , the gauge field is welldescribed by linear perturbations because (for example, see Eq. (15)). In this section, we study the highfrequency regime analytically.
We suppose that unperturbed background solutions are the vacuum state without any electric field, that is, the trivial solutions and . Perturbed solutions are given by and . For convenience, will be treated as a complex function including its phase and amplitude. Note that we can set for the background solution without loss of generality because of the scaling symmetry (26) of the current systems.
The nonlinear equations of motion (19) and (21) can be linearized, so that the linearized equations of motion become for ,
(40) 
and for ,
(41) 
Since is real, the complex conjugate satisfies the same equation of motion as (40).
General solution of Eq. (40) is given by
(42) 
where
(43)  
and . Each particular solution is a real function and has the following asymptotic behavior:
(44) 
and
(45) 
Thus, in terms of asymptotic behavior near the AdS boundary , and describe socalled nonnormalizable modes and normalizable modes, respectively. At we find is regular while is singular. This fact is closely related to whether the solution describes the insulator or conductive phase (namely, the Minkowski or black hole embeddings).
For later convenience, using perturbative quantities we write down the conserved quantity defined by Eq. (18), which corresponds to the Joule heating in the boundary theory, as
(46) 
Since it is obviously proportional to the Wronskian of solutions of the linear differential equation (40), is conserved at this order as well as is conserved. Substituting the general solution of Eq. (48), we explicitly obtain
(47) 
General solution of Eq. (41) is
(48) 
where is an integration constant equivalent to a boundary condition .
6.1 The insulator phase
As seen previously, the Minkowski embedding solution describing the insulator phase can reach to