Instruction: Write an abstract for the following acadamic paper.

Input:
Paper: we have recorded the flight of a fly during take off and landing using digital high speed photography . it is shown that the dynamics of flexible wings are different for these procedures . during this observation fly flew freely in a big box and it was not tethered .
Abstract: in this fluid dynamics video , we demonstrated take off and landing of a fly . the deformation of wings is in focus in this video .

Paper: since its prediction in 1954 @xcite , the effect of superradiance has attracted a lot of theoretical and experimental attention @xcite . its basic principle  the fact that a coherent interaction of an unexcited gas with the vacuum of a common field can create a spontaneous macroscopic excitation of both matter and field subsystems  appears in various incarnations in diverse branches of physics @xcite . the dicke model @xcite of superradiance resorts to a maximal simplification of the problem to capture the main features of the superradiant transition in the most transparent way . the model shows a thermal phase transition , analyzed and discussed in refs.@xcite , as well as a zero - temperature ( ground - state ) quantum phase transition ( qpt ) , which was realized experimentally with the aid of a superfluid gas in a cavity @xcite . recent theoretical analyses showed that the model exhibits also a novel type of criticality  the one observed in the spectrum of excited states @xcite . these so - called excited - state quantum phase transitions ( esqpts ) affect both the density of quantum levels as a function of energy and their flow with varying control parameters , and are present in a wide variety of quantum models with low numbers of degrees of freedom @xcite .    in this work , we analyze properties of a simple superradiance model interpolating between the familiar dicke @xcite and tavis - cummings @xcite hamiltonians . a smooth crossover between both limiting cases is achieved by a continuous variation of a model parameter , which allows one to observe the system s metamorphosis on the way from the fully integrable ( hence at least partly understandable ) to a prevailingly chaotic ( so entirely numerical ) regime . one of the aims of our work is to survey the phase - transitional properties of the extended model and to investigate the nature of its _ quantum phases_the domains of the excited spectrum in between the esqpt critical borderlines . a usual approach to characterize different phases related to the ground state of a quantum system makes use of suitable order parameters , i.e. , expectation values of some observables . we show that an unmistakable characterization of phases involving excited states is not achieved through the expectation values alone but rather through their different smoothed energy dependences ( trends ) .    the main emphasis of our analysis is nevertheless put on the entanglement properties of excited states across the whole spectrum and their potential links to esqpts and quantum phases of the model . it is known that a continuous ground - state qpt in many models ( including the present one ) is characterized by a singular growth of entanglement within the system , which can be seen as a quantum counterpart of the diverging correlation length in continuous thermal phase transitions @xcite . a question therefore appears whether there exist any entanglement - related signatures of esqpts . the extended dicke model is rather suitable for a case study of this type since it allows one to analyze at once various types of entanglement  that between the field and all atoms , and that between individual atoms . it should be stressed that robust , systematic studies of entanglement in excited spectra of any systems are rather scarce . we hope that results of our analysis , which is unavoidably model - specific , will initiate similar studies in other relevant systems . the plan of the paper is as follows : basic quantum and classical features of the model are described in sec.[edm ] . thermal and quantum critical properties and a classification of thermodynamic and quantum phases are presented in sec.[pt ] . the atom  field and atom  atom entanglement properties are investigated in sec.[ent ] . conclusions come in sec.[sum ] . consider single - mode electromagnetic field with photon energy @xmath0 ( polarization neglected ) interacting with an ensemble of @xmath1 two - level atoms , all with the same level energies @xmath2 . the free hamiltonian reads as @xmath3 where the first term describes the field of photons , created and annihilated by bosonic operators @xmath4 and @xmath5 , respectively , and the second term represents the ensemble of atoms , with @xmath6 standing for the third pauli matrix acting in the space of the @xmath7th atom . in the following we use all pauli matrices @xmath8 with the subscripts @xmath9 or @xmath10 . the interaction term is taken in the dipole approximation , with the dipole operator for the @xmath7th atom proportional to @xmath11 . the size of the atomic ensemble is assumed to be much smaller than the wavelength of photons ( cavity size ) . then @xmath12 where the coupling parameter ( in units of energy ) can be written as @xmath13 , with @xmath14 standing for the electric dipole moment matrix element , @xmath15 for the vacuum permitivity , and @xmath16 for the cavity volume . the hamiltonian @xmath17 represents the well known dicke model @xcite . it is commonly invoked to explain basic features of the transition of matter to superradiant phase @xcite . indeed , as shown below , if the coupling strength @xmath18 is sufficiently increased , the superradiant phase with @xmath19 and @xmath20 is created at low temperatures . note that the diamagnetic term @xmath21 is neglected here though it turns out to be very important in the discussion of superradiance in realistic situations connected to cavity electrodynamics @xcite . we emphasize that recent experimental realizations of the dicke hamiltonian @xcite do not implement the above electrodynamic content , nevertheless in this presentation we keep this original formulation as a widely comprehensible background . introducing collective quasi - spin operators for the whole atomic ensemble ( which is possible due to the long wavelength approximation ) , @xmath22 one immediately sees that the dicke hamiltonian conserves the squared quasi - spin @xmath23 with eigenvalues @xmath24 , where @xmath25 is integer for @xmath1 even or half - integer for @xmath1 odd . the full atomic hilbert space is the span of all @xmath26 possible configurations of atoms , thus @xmath27 , where @xmath28 is the space of the @xmath7th atom . however , due to the conservation of @xmath29 , the dynamics can be investigated separately in any of the single-@xmath25 subspaces @xmath30 with dimension @xmath31 . the decomposition reads as follows @xcite : @xmath32 where @xmath33 enumerates replicas ( their number is @xmath34 ) of the space with given @xmath25 differing by the exchange symmetry of atomic components . in the following , we will investigate thermodynamic properties of the model in the full atomic space @xmath35 , as well as quantum properties in a single-@xmath25 space @xmath36 . the most natural choice in the latter case is the unique ( @xmath34=1 ) subspace with maximal value @xmath37 , which is fully symmetric under the exchange of atoms and therefore emphasizes the collective character of the superradiance phenomenon . a general-@xmath25 subspace has a mixed exchange symmetry such that only a number of @xmath38 atoms can be excited independently , while excitations of the remaining @xmath39 atoms have to compensate each other ( in ref.@xcite the quantum number @xmath25 is called a cooperation number of the atomic gas ) @xcite . the reduced single-@xmath25 model has only two effective degrees of freedom , one associated with the bosonic field , the other with the su(2 ) algebra of collective quasi - spin operators ( [ jbul ] ) , hence @xmath40 . in contrast , in the full ( all-@xmath25 ) model each su(2)@xmath41 algebra of atom @xmath7 brings an independent degree of freedom , so the whole atom  field system has @xmath42 . the dicke model is sometimes simplified by omitting the terms @xmath43 and @xmath44 , that for @xmath45 and very small @xmath18 yield negligible contributions to the transition amplitudes @xcite . the reduced model , in case of @xmath46 atoms called the tavis - cummings model @xcite , has an advantage of being completely integrable because of the existence of an additional integral of motion : @xmath47 where @xmath48 is the number of bosons and @xmath49 the number of excited atoms ( taking values @xmath50 $ ] , where @xmath51 is the @xmath52 quantum number ) . this conservation law follows from a u(1 ) symmetry under gauge  transformations consisting of a phase shift @xmath53 and a simultaneous rotation @xmath54 , where @xmath55 is the rotation matrix by angle @xmath56 around axis @xmath57 . the quantity ( [ m ] ) is not conserved in the dicke model , although the full hamiltonian has a discrete z(2 ) symmetry with conserved parity  @xmath58 .    in this work , we analyze properties of the single - mode superradiance model across the transition between the tavis - cummings and dicke regimes . following refs.@xcite , we introduce a second interaction hamiltonian , @xmath59 where @xmath60 $ ] is an additional control parameter , and compose the total hamiltonian as @xmath61 , so : @xmath62 the tavis - cummings and dicke limits of the model are obtained for @xmath63 and 1 , respectively . we see that an increase of @xmath64 from zero breaks the u(1 ) symmetry while preserving the z(2 ) symmetry of the @xmath65 model . extended dicke model as a function of @xmath18 for @xmath66 and @xmath67 . the panels correspond to ( a ) the tavis - cummings limit @xmath68 , ( b ) intermediate regime @xmath69 and ( c ) the dicke limit @xmath70 . ] figure  [ f_spe ] shows quantum spectra of hamiltonian with @xmath37 depending on the interaction parameter @xmath18 for three different values of @xmath64 . the evolution of the ground state , namely its sudden drop to lower energies above a certain value of @xmath18 , indicates the superradiant transition at zero temperature . we note that the eigensolutions of the model for a given @xmath25 are in general obtained by a numerical diagonalization of hamiltonian in the tensor product of atomic and field hilbert spaces , @xmath71 , where @xmath72 is the finite atomic space from eq . spanned by basis vectors @xmath73 with quasi - spin projections @xmath74 and @xmath75 is the infinite - dimensional space of the boson field spanned by vectors @xmath76 with boson numbers @xmath77 . since the diagonalization procedure requires a truncation of @xmath75 , the eigensolutions must be carefully checked for convergence . a general state vector of the model for any fixed @xmath25 is expressed as @xmath78 where @xmath79 are expansion coefficients satisfying the normalization condition @xmath80 . note that for even- and odd - parity states , respectively , the sums in eq . go either over even or odd values of @xmath81 . any vector can be visualized as a wave function @xmath82 , where the angle @xmath83 comes from the quasi - spin representation of the atomic subsystem and variable @xmath84 from the oscillator representation of the bosonic field . for integer @xmath25 ( even @xmath1 ) , the wave function can be obtained by substitutions ( up to normalization constants ) @xmath85 and @xmath86 , with @xmath87 standing for the hermite polynomial . for half - integer @xmath25 ( odd @xmath1 ) , however , the wave function @xmath82 acquires a spinorial character . a possible representation can be obtained by coupling an integer angular momentum @xmath88 ( the even core of @xmath89 atoms ) with @xmath90 ( the odd atom ) to the total @xmath25 as in the case of spinor spherical harmonics . this leads to a mapping of @xmath73 to a 2-valued function @xmath91 , which is @xmath92-periodic , with @xmath93 denoting the corresponding clebsch - gordan coefficients . examples of squared wave functions of the @xmath37 hamiltonian eigenstates for an even @xmath1 in the medium-@xmath64 regime are depicted in fig.[f_wav ] . they are taken at the parameter and energy values indicated by dots in fig.[f_qpd]d ( @xmath94 ) . all eigenstates have a good parity , which means that their wave functions are symmetric or antisymmetric under the transformation @xmath95 and @xmath96 .     of selected eigenvectors of the hamiltonian with @xmath97 , @xmath66 for @xmath69 taken at the @xmath98 values indicated by dots in fig.[f_qpd]d . panels ( a)(f ) , respectively , correspond to the vertical row of points at the right of fig.[f_qpd]d from bottom to top [ panel ( a ) shows the ground state ] , panels ( g),(h ) to the left pair of points . note that variable @xmath99 is @xmath92 periodic . ] finding eigensolutions of the model is much simpler in the @xmath63 limit . in that case , the conserved quantity from eq . splits the hilbert space into a sum @xmath100 of dynamically invariant subspaces @xmath101 with different values of quantum number @xmath102 . each of these subspaces is spanned by states @xmath103 where @xmath51 goes from @xmath104 to @xmath105 and @xmath106 takes values from @xmath107 to @xmath108 , so its dimension is : @xmath109 note that the saturation of the @xmath101 dimension at @xmath110 for @xmath111 results from the finiteness of the atomic space @xmath36 . the full spectrum for @xmath63 is comprised of a number of mutually non - interacting spectra with different values of @xmath108 , each of them being obtained by truncation - free diagonalization in a finite dimension . as seen in panel ( a ) of fig.[f_spe ] , the @xmath63 spectrum indeed shows numerous crossings of levels with different @xmath108 . in contrast , an increase of @xmath64 in panels ( b ) and ( c ) leads to strengthening of level repulsion , which is a quantum signature of chaos . spectrum ( for @xmath97 ) in the @xmath112 subspace for the @xmath68 model with @xmath66 ( panel a ) , @xmath113 ( panel b ) and @xmath114 ( panel c ) . the wave - function entropies @xmath115 of individual eigenstates in the unperturbed basis @xmath103 are encoded in the varying shade of the corresponding lines . ] the structure of the spectrum for @xmath63 can be further elucidated by considering the hamiltonian with @xmath116 ( the tuned case ) . in a single-@xmath108 subspace it reads as @xmath117 where @xmath118 is the energy in absence of interaction and @xmath119 are the interaction matrix elements in the unperturbed basis @xmath120 enumerated by @xmath121 . a simple calculation yields : @xmath122 . the eigenvectors of this hamiltonian obviously do not depend on @xmath18 since both terms in eq . the eigenvalues have the form @xmath123 , where @xmath124 are eigenvalues of the matrix in the second term including the @xmath125 factor . since the characteristic polynomial of this matrix contains either even ( for @xmath126 even ) or odd ( for @xmath126 odd ) powers of @xmath127 , the eigenvalues come in pairs with opposite signs . for @xmath126 odd , there is an additional unpaired eigenvalue @xmath128 . therefore , the spectrum linearly expands with increasing @xmath18 around the non - interacting energy @xmath118 , while the eigenvectors remain unchanged . both these conclusions concerning the tuned system with @xmath63 are verified in panel ( a ) of fig.[f_cary ] , which besides the energy spectrum of a single-@xmath108 subspace ( the linearly diverging bundle of levels ) shows also the wave - function entropy @xmath129 of individual eigenstates ( encoded in the shade of the respective line ) . this quantity measures the localization of instantaneous eigenstates in the unperturbed basis @xmath103 , ranging from @xmath130 for any unperturbed state to @xmath131 for a uniformly distributed superposition @xcite . in sec.[ent ] we will see that for @xmath63 the wave - function entropy quantifies the atom  field entanglement contained in the given state , but here we use it just as an indicator of structural evolution of eigenstates with varying @xmath18 . in fig.[f_cary ] we see  in accord with the above analysis  that the wave - function entropies of all states remain constant . it is clear that the full spectrum in fig.[f_spe]a results just from a pile up of the trivial spectra similar to that in fig.[f_cary]a with all values of @xmath108 . in contrast , the spectrum of a detuned hamiltonian ( with @xmath132 ) is less trivial . examples of level dynamics in a single-@xmath108 subspace and the corresponding wave - function entropies for @xmath133 and @xmath134 are shown in panels ( b ) and ( c ) of fig.[f_cary ] , where we see that the eigenvalues are no more straight lines and that the eigenvectors change with @xmath18 . the peculiar shapes of these spectra will be subject of our interest in sec.[qpt ] . here we only point out that the tuned situation is restored in an approximate sense for very large @xmath18 , when the last term in the detuned hamiltonian @xmath135 , with @xmath136 , becomes negligible . therefore , investigating properties of the model for large coupling strengths , one can make an assumption that @xmath45 . this conclusion , which holds independently of @xmath64 , will be seen in action in sec.[afe ] . a straightforward way to the classical limit of the extended dicke model leads via the following operator to c - number mappings : @xmath137 the spherical angles @xmath138 and @xmath139 $ ] determine the orientation of the classical quasi - spin vector @xmath140 . note that angle @xmath141 is measured here from the south to the north pole , thus a fully de - excited state of atoms ( @xmath142 ) corresponds to @xmath143 while a maximally excited state of atoms ( @xmath144 ) is associated with @xmath145 . it can be shown that the pair of variables @xmath99 and @xmath146 , respectively , represents canonically conjugate coordinate and momentum associated with the collective states of the atomic subsystem . similarly , quantities @xmath147 and @xmath148 in eq.([map2 ] ) are suitable coordinate and momentum of the field subsystem . let us point out that while the atomic phase space for each @xmath1 is finite ( it is a ball with radius @xmath25 ) , the field phase space covers the whole plane , expressing the fact that the number of bosons is unlimited . the classical hamiltonian resulting from substitution of eqs.([map1 ] ) and ( [ map2 ] ) to reads as ( see also ref.@xcite ) : @xmath149 \label{hcl } \,,\ ] ] since the classical limit of the model is reached in the limit @xmath150 , the above expression needs an appropriate scaling . it is achieved by introducing quantities @xmath151 , @xmath152 , @xmath153 and @xmath154 that for any value of the size parameter take values of order @xmath155 . this leads to the following size - invariant expression : @xmath156 \label{hclscal } \,.\ ] ] note that eq . contains an effective frequency @xmath157 and an effective coupling strength @xmath158 , which both increase linearly with @xmath25 and reach values of the bare parameters @xmath159 and @xmath18 for the maximal @xmath160 . stationary points of the classical hamiltonian will be analyzed below to determine critical values of control parameters and energy separating distinct quantum phases of the model . we know that for @xmath63 the hamiltonian is integrable as it preserves the quantity @xmath108 from eq .. this introduces a specific correlation between the degrees of freedom associated with both atomic and field subsystems . as a consequence , the system with @xmath40 degrees of freedom is for any fixed value of @xmath108 reduced to a system with just a single effective degree of freedom . to make this explicit , we employ a transformation : @xmath161 where @xmath162 is conserved . it can be checked that this transformation is canonical , so @xmath163 and @xmath164 are new pairs of conjugate coordinates and momenta . the @xmath163 variables for each @xmath165 satisfy a constraint @xmath166 $ ] , so they form a disc ( for @xmath167 ) or an annulus ( for @xmath168 ) . the classical tavis - cummings ( @xmath63 ) hamiltonian in terms of the new variables becomes : @xmath169 } \ x ' \label{hcl1d } \,,\ ] ] which for a constant @xmath165 depends only on @xmath163 , so has effectively @xmath170 degree of freedom . this not only guarantees the full integrability of the @xmath63 model , but also results in some emergent critical phenomena which will be discussed below . maps of the single-@xmath25 model with @xmath66 for six values of parameter @xmath171 $ ] . the tones of blue indicate the classical regular fraction . the curves correspond the ground state energy @xmath172 ( full line ) and the critical esqpt energies @xmath173 ( dotted ) , @xmath174 ( dashed ) , and @xmath175 ( dash - dotted ) from eqs.([ec1])([ec3 ] ) . acronyms d ( dicke ) , tc ( tavis - cummings ) , n ( normal ) and s ( saturated ) indicate different quantum phases  of the model ( see sec.[qpt ] ) . dots in panel ( d ) mark places where wave functions in fig.[f_wav ] are taken . ] variation of parameter @xmath64 from 0 to 1 breaks the integrability of hamiltonian and leads to increasingly chaotic dynamics . the evolution of a classical measure of chaos with @xmath64 is demonstrated in fig.[f_qpd ] . it displays maps of so - called regular fraction @xmath176 in the plane @xmath98 for various @xmath64 . the regular fraction is defined as the volume of the phase space domain filled with regular trajectories relative to the total phase space volume available at given energy @xmath177 : @xmath178 here , @xmath179 stands for a characteristic function defining the regular part of the phase space , i.e. , the domain where classical dynamics is regular in the sense of vanishing lyapunov exponents ( @xmath180 in the regular part and @xmath181 in the chaotic part ) . the way how the regular fraction is obtained from numerical simulation of classical dynamics is described in ref.@xcite . the value of the regular fraction in fig.[f_qpd ] is encoded into the tones of blue : the full color indicates perfectly regular areas ( @xmath182 ) and white completely chaotic ones ( @xmath183 ) . we observe that the degree of chaos exhibits an overall increase with @xmath64 . however , even in the dicke limit the model is not entirely chaotic . the most chaotic domain in all @xmath184 panels is that with parameter @xmath18 above the superradiant transition and energy @xmath177 exceeding a certain value above the ground - state @xmath172 . besides the main trends in the dependence of @xmath176 we also observe some surprising fine structures  for instance the ribs  in panels ( d)(f ) . note that various aspects of chaos in the extended dicke model were recently studied in refs.@xcite . the extended dicke hamiltonian ( [ h ] ) has a complex phase structure @xcite . in this section , we analyze thermal phase transition of the full model with atomic hilbert space ( [ hil ] ) and the number of degrees of freedom @xmath42 . following the standard approach described in ref.@xcite , we use the glauber coherent states @xmath185 , where @xmath186 is an eigenvalue of the @xmath5 operator , @xmath187 is a normalization factor , and @xmath188 denotes the field vacuum . an average number of bosons @xmath189 in the state @xmath190 represents an indicator of the phase transition from normal to superradiant phase : it is zero in the normal phase , while in the superradiant phase we expect @xmath191 . on the way to the thermodynamic limit @xmath192 , it is convenient to perform the scaling transformation @xmath193 ( cf.sec.[cl ] ) , so that the scaled coherent - state variable @xmath194 taking values between 0 and @xmath195 becomes a suitable complex order parameter of the superradiant phase transition . we start by introducing two special values of the atom  field coupling strength : @xmath196 while the first one will be shown to represent a critical strength for the occurrence of the superradiant phase , the second value will turn out to be a strength at which a specific form of the superradiant phase occurs for @xmath197 . to get an equilibrium field configuration at a given temperature @xmath198 ( measured in energy units ) , we minimize in variable @xmath199 a scaled free energy @xmath200 , where @xmath201 is the partition function , calculated for a given @xmath199 from a partial trace of the hamiltonian exponential over the full atomic hilbert space ( [ hil ] ) . for hamiltonian ( [ h ] ) this procedure yields : @xmath202 \\ \lambda\left[(1\+\delta){\bar\alpha'}\!\-i(1\-\delta){\bar\alpha''}\right ] & \frac{1}{2}\omega_0 \end{array}\right ) \right\ } \label{z } \,.\ ] ] the two eigenvalues of the matrix in the exponential are : @xmath203}\equiv\pm e_{{\bar\alpha } } \label{e } \,,\ ] ] which leads to a simple expression of the scaled free energy : @xmath204    stationary points of @xmath205 can be obtained from a simple analysis . the first derivative in variable @xmath206 or @xmath207 , respectively , vanishes at the points satisfying the first or second line of the following array : @xmath208 the point @xmath199=0 is always a trivial solution of both lines , but additional solutions appear for @xmath209 . this coupling strength represents a critical point where the normal phase @xmath210 becomes unstable at the lowest temperatures and a new , superradiant equilibrium is created at some non - zero values of @xmath199 . the critical temperature for the superradiant transition , i.e. , the upper temperature limit for the existence of the @xmath211 solution in the region @xmath209 , is : @xmath212 for @xmath213 , the stable equilibrium of the system is again only the normal solution @xmath210 . model ( @xmath66 , @xmath69 ) and sample landscapes of the scaled free energy at the @xmath214 cut of the phase diagram ( vertical line ) . the full curve in the phase diagram marks the critical temperature @xmath215 for the transition to the normal ( n ) phase , while the dotted curve indicates the temperature @xmath216 separating the tc and d superradiant phases ( with or without saddles of free energy ) . the free energy landscapes in the complex-@xmath217 plane are , for selected temperatures , visualized by contour - shade plots ( darker areas indicate lower values ) and by the cuts along real and imaginary axes . ]    for @xmath63 ( the tavis - cummings limit ) , both right - side conditions in eqs.([con1 ] ) and ( [ con2 ] ) become identical and yield a solution @xmath218 that grows from zero with increasing difference @xmath219 . therefore , the superradiant minimum of the free energy forms a circle around @xmath210 . for @xmath184 , however , a simultaneous solution of both right - side conditions in ( [ con1 ] ) and ( [ con2 ] ) is no more possible . the superradiant equilibrium is then represented by a pair of points on the real axis , @xmath220 , where @xmath221 solves the right - side equation in ( [ con1 ] ) . if @xmath222 and @xmath223 , a pair of saddle points appears for low temperatures on the imaginary axis , @xmath224 , where @xmath225 solves the right - side equation in ( [ con2 ] ) . these unstable solutions exist for temperatures below @xmath226 a possibility of thermodynamic quasi - equilibrium states associated with the saddle points of free energy was recently discussed in ref.@xcite .    a thermal phase diagram in the @xmath227 plane for @xmath94 is shown in fig.[f_tpd ] along with samples of the free energy landscapes in various phases and at @xmath228 . the critical temperature @xmath215 from eq.([tc ] ) determines the phase transition between the superradiant and normal ( acronym n ) phases of the model . the superradiant phase exists in two forms : the tavis - cummings phase ( acronym tc ) with the saddles in the free energy landscape and the dicke phase ( acronym d ) without the saddles . while the tc phase is the only type of superradiant phase in the tavis - cummings limit ( where @xmath229 and @xmath230 ) , the d phase is exclusive in the dicke limit ( in which @xmath231 ) . for intermediate @xmath64 both phases coexist , being separated by temperature @xmath216 from eq.([t0 ] ) . let us stress that the tc@xmath232d transition ( in contrast to d@xmath232n ) is not a phase transition in the standard sense since it does not affect the global minimum of the free energy . we return now to the analysis of the extended dicke model with the restriction to a single-@xmath25 collective subspace of atomic states . we assume in particular @xmath37 , but in the following we keep @xmath25 as a changeable parameter . the number of relevant degrees of freedom of the single-@xmath25 model is @xmath40 , independently of the size parameter @xmath1 , which implies that the infinite - size limit , @xmath192 , coincides with the classical limit @xcite . both ground - state and excited - state quantum phase transitions can be predicted from the classical version of the model , namely from the behavior of stationary points of the classical hamiltonian . in particular , the esqpts result from singularities ( non - analyticities ) in the semiclassical density of states @xmath233 where @xmath234 is the step function ( @xmath235 for @xmath236 and @xmath237 for @xmath238 ) , that appear at the points satisfying @xmath239 ( with @xmath240 standing for the gradient in the phase space ) .    it has been shown @xcite that the esqpts caused by non - degenerate stationary points  those which are locally quadratic as their hessian matrix of second derivatives has only non - zero eigenvalues  can be classified by a pair of numbers @xmath241 , where a so - called index of the stationary point @xmath242 is a number of negative eigenvalues of the hessian matrix . in particular , for @xmath243 , the first derivative of the level density in a vicinity of a stationary - point energy @xmath244 behaves as : @xmath245 hence exhibits either a jump ( for @xmath242 even ) or a logarithmic divergence ( for @xmath242 odd ) at the critical energy @xmath244 .     from eq.([trans ] ) in the phase space of the atomic subsystem defined by spherical angles @xmath141 and @xmath99 for @xmath66 and @xmath69 . various panels correspond to the indicated values of the coupling strength @xmath18 . each contour - shade polar plot ( with @xmath246angle and @xmath247radius ) is accompanied by the corresponding horizontal and vertical cuts . ]    in the present case , the evaluation of the level density according to eq . can be simplified with the aid of the volume - preserving substitution @xmath248 which transforms the classical hamiltonian to the form @xmath249 with fully separated variables of the field and atomic subsystems . although the mapping does _ not _ represent a canonical transformation ( so @xmath250 is not a new coordinate - momentum pair  otherwise the model would be separable and thus fully integrable for any @xmath64 ) , it simplifies the analysis of the level density . in particular , the integration in eq . over variables @xmath250 , on which the transformed hamiltonian depends quadratically , can be performed explicitly ( cf . ref.@xcite , where an analogous calculation is performed for a hamiltonian with a quadratic kinetic term ) . the result is a simplified expression @xmath251 in which the integration , involving function @xmath252 defined in eq . , goes only over the 2-dimensional collective phase space of the atomic subsystem  a ball with radius @xmath25 . eq . is proportional to an area of the ball region where @xmath253 takes values less than ( or equal to ) the chosen energy @xmath177 , hence it can be visualized as flooding of a landscape with profile @xmath253 on a globe . in analogy to eq . , the @xmath253 function can be expressed in a scaled form , @xmath254 , with parameters @xmath159 and @xmath18 replaced by their effective values @xmath255 and @xmath256 depending on @xmath25 . the @xmath253 function is for selected values of control parameters @xmath18 and @xmath64 depicted in fig.[f_impo ] . , for the single-@xmath25 model with @xmath66 and @xmath69 as a function of @xmath18 and @xmath257 . the shade diagram in the left panel was obtained from eq . ; darker areas represent larger values and vice versa . panels ( a ) and ( b ) show cuts ( full curves ) at ( a ) @xmath258 and ( b ) @xmath214 in comparison with finite - size results ( dashed curves ) based on numerical diagonalization of the hamiltonian for @xmath97 and gaussian smoothening of the spectrum ( @xmath259 and @xmath260 for @xmath258 and @xmath261 , respectively ) . scales in panels ( a),(b ) are arbitrary and not the same . ] it is clear from the expression that the stationary points of @xmath262 in all four variables correspond to @xmath263 and @xmath264 determined as stationary points of the function @xmath253 . since the quadratic minimum in variables @xmath250 has a null index , any stationary point of @xmath262 has an index @xmath242 equal to that of the corresponding stationary point of @xmath253 . for the determination of esqpts it is therefore sufficient to find and classify stationary points of the function @xmath253 on a ball . non - degenerate stationary points of @xmath253 with index @xmath265 , @xmath266 , and @xmath267 cause , respectively , an upward jump , logarithmic divergence and a downward jump of @xmath268 .    inserting the effective parameters @xmath256 and @xmath255 from the scaled hamiltonian into eq . , we obtain the following @xmath25-dependent values of the critical couplings : @xmath269 these will play the roles of @xmath270 and @xmath271 in individual subspaces of states with fixed values of @xmath25 . note that in the highest subspace with @xmath160 both expressions in eq . yield the bare values of critical couplings . a straightforward analysis leads to the following conclusions :    \(i ) stationary points with @xmath265 represent the global minimum of both @xmath253 and @xmath262 functions , demarcating the ground state of the @xmath192 system . the minimum appears at @xmath272 ( @xmath99 arbitrary ) for @xmath273 , and at @xmath274 , @xmath275 and @xmath276 ( a pair of degenerate minima ) for @xmath277 . the ground - state energy is given by the following formula : @xmath278 & { \rm for\ } \lambda\in[\lambda'_c,\infty)\ , , \end{array } \right . \label{e0}\ ] ] which exhibits a second - order qpt from normal to superradiant ground - state phase at the critical coupling @xmath279 , where @xmath280 has a discontinuity . \(ii ) stationary points with @xmath281 represent saddles of @xmath253 . they are located at @xmath272 ( @xmath99 arbitrary ) for @xmath282 , and at @xmath283 , @xmath284 and @xmath285 for @xmath286 . these stationary points correspond to an esqpt ( a logarithmic divergence of @xmath268 ) at energy : @xmath287 & { \rm for\ } \lambda\in[\lambda'_0,\infty)\ , . \end{array } \right . \label{ec1}\ ] ] for @xmath288 we have @xmath289 and eq.([ec1 ] ) is reduced to its first line . \(iii ) stationary points with @xmath290 are maxima of @xmath253 at @xmath272 ( @xmath99 arbitrary ) for @xmath286 and at @xmath291 ( @xmath99 arbitrary ) for @xmath292 . related esqpts ( downward jumps of @xmath268 ) appear at energies @xmath293 the second maximum of @xmath253 at energy is the global one , so for @xmath294 the formula yields a constant ( saturated ) value of the level density equal to @xmath295 . subspace of the @xmath296 tavis - cummings model ( @xmath68 ) with @xmath113 as a function of @xmath18 and @xmath297 ( shade plot on the left , with dark areas indicating larger values and vice versa ) and three dependences of the corresponding hamiltonian function on @xmath298 for @xmath299 ( on the right ) . an esqpt due to an inflection point of the classical hamiltonian above the qpt critical point results in a logarithmic divergence of the semiclassical level density ( the dark band in the shade plot ) , cf.fig.[f_cary]b . ] all esqpt critical borderlines @xmath173 , @xmath174 and @xmath175 from eqs. for various values of @xmath64 are demarcated in fig.[f_qpd ] above . their existence is numerically verified in fig.[f_dender ] , which depicts the @xmath98 dependence of @xmath300 for @xmath37 somewhere in between the dicke and tavis - cummings limits . the shade plot was obtained through the phase - space integration in eq . and a comparison with results obtained by a numerical diagonalization is shown in the two panels on the right . we stress that due to the @xmath40 character of our model , the esqpt singularities occur in the _ first derivative _ of the level density with respect to energy . however , the conservation of quantity in the @xmath63 limit and the corresponding reduction of the number of effective degrees of freedom to @xmath170 ( see sec.[cl ] ) leads to a possibility to generate an esqpt singularity in the level density itself . such an effect can really be identified within a particular _ single @xmath108-subspace _ of the full model . the lowest - energy state of any of such subspaces for @xmath192 coincides with the global minimum of classical hamiltonian within the available ( for a given @xmath108 ) domain of the phase space . it turns out that the most interesting @xmath108-subspace is the one with @xmath301 . this particular subset of states exhibits for @xmath133 a second - order ground - state qpt at the coupling strength @xmath18 equal to a critical value : @xmath302 at this coupling , the main minimum of function moves from @xmath303 away to @xmath304 while @xmath305 becomes an inflection point . the inflection point is present for @xmath306 and generates an esqpt ( a logarithmic divergence of the level density in the @xmath301 subspace ) at energy @xmath175 from eq .. note that for @xmath134 the esqpt appears for @xmath307 at the upper energy of the unperturbed spectrum . this effect , described already in ref.@xcite , is demonstrated in fig.[f_tcdens ] . it displays a shade plot of the semiclassical density of levels with @xmath301 ( obtained by the @xmath170 phase  space integration ) and the underlying forms of the hamiltonian . a finite - size sample of the predicted singularity in the level density was seen in panel ( b ) of fig.[f_cary ] above . a more detailed analysis of eq . shows that the present type of criticality is absent in the @xmath308 subspaces . this can be intuitively understood from the evaluation of matrix elements @xmath309 of the interaction hamiltonian @xmath310 with @xmath63 in the unperturbed basis ; see eq . and the text below it . the interaction matrix elements quantify the mixing induced by @xmath310 in the unperturbed eigenbasis of @xmath311 . neglecting the trivial state with @xmath312 which does not mix at all , the matrix elements @xmath309 are particularly small in the @xmath313 subspace ( dimension 2 ) and for the @xmath314 state of the @xmath301 subspace ( dimension @xmath31 ) . only the latter state can develop a singularity in the @xmath150 limit . it is the state @xmath315 , which for @xmath133 represents the lowest state of the @xmath301 subspace , while for @xmath134 it is the highest state ; cf.panels ( b ) and ( c ) of fig.[f_cary ] . although the @xmath308 subspaces show no quantum critical effects , it is interesting to realize that a pile up of all @xmath108-subspaces produces the downward jump of @xmath300 ( with @xmath316 standing for the total semiclassical level density ) , as observed in the @xmath317 esqpt of the @xmath63 system . the critical borderlines @xmath173 , @xmath174 and @xmath175 in the @xmath98 plane separate spectral domains that we consider to constitute different quantum phases  of the model ( in analogy to thermodynamic phases ) . they are denoted by acronyms d ( dicke ) , tc ( tavis - cummings ) , n ( normal ) and s ( saturated ) , see fig.[f_qpd ] . both phases d and tc , which coexist for @xmath64 in between 0 and 1 , contain quantum states of a superradiant nature because their energy is lowered with respect to the minimal non - radiant value @xmath318 . in the limits of @xmath64 , the model shows only one type of the superradiant phase : tc for the tavis - cummings limit @xmath63 and d for the dicke limit @xmath319 . in contrast , both phases n and s contain states that resemble excitations in the fully non - radiant regime at @xmath320 . the s phase above @xmath321 yields a constant , saturated value of level density @xmath316 . ( panel a ) and @xmath322 ( panel b ) in individual eigenstates of the @xmath296 hamiltonian ( [ h ] ) for @xmath97 , @xmath66 , @xmath214 . the upper panels depict lattices for the limiting values of @xmath68 and 1 , whereas the lower panels correspond to @xmath69 . the insets of the @xmath69 panels show smoothed dependences of both observables on energy ( averages over 20 neighboring eigenstates and their piecewise fits in various quantum phases ) . ] while the definition of esqpts is obvious from the behavior of quantities like level density at the critical energies , the meaning of quantum phases in between the critical energies is not a priori clear . it should be looked for in the structure of the energy eigenstates in the corresponding energy domains . however , sample wave functions in fig.[f_wav ] [ where panels ( a ) , ( b ) , ( d ) , ( f ) and ( h ) correspond to phases d , d , tc , n and n , respectively ] indicate that individual eigenstates do not show sufficient clues for the identification of phases . we have to resort to a more efficient visualization tool allowing us a bulk inspection of the eigenstate properties . the approach employed in this work is based on the method of so - called _ peres lattices _ @xcite . it was proven useful in various models ( see ref.@xcite and references therein ) including the dicke model @xcite . a general peres lattice shows expectation values @xmath323 of a selected observable @xmath324 in individual energy eigenstates @xmath325 ( enumerated by integer @xmath326 ) arranged into lattices with energy @xmath327 on one of the axes . examples of peres lattices for observables ( a ) @xmath328 and ( b ) @xmath329 are shown in the respective panels of fig.[f_pere ] . a comparison of both @xmath94 lattices ( in the main panels ) with the map of chaos in fig.[f_qpd]d ( for @xmath18 demarcated by the second vertical line ) yields an overall agreement with the original conjecture concerning the correlation of orderly arranged domains of a generic lattice with more regular regions of classical dynamics and vice versa @xcite . also the limiting @xmath63 lattices in the upper panels of fig.[f_pere ] are perfectly regular , while those for @xmath319 above a certain energy are almost fully chaotic . the @xmath94 peres lattices in fig.[f_pere ] cross all four quantum phases of the model . there exist apparent similarities between parts of both @xmath94 lattices located in the tc and d phases , respectively , and the corresponding tavis - cummings and dicke lattices displayed in the @xmath63 and 1 upper panels . however , a more specific distinction of quantum phases results from an averaging of both lattices over the neighboring eigenstates , as shown in the insets of both lower panels of fig.[f_pere ] . points in these _ smoothed lattices _ represent averages @xmath330 and @xmath331 of the respective expectation values @xmath332 and @xmath333 over 20 neighboring eigenstates . smoothed dependences of the averaged quantities on energy are given by lines , resulting from piecewise fits within the four quantum phases . we observe that various quantum phases are recognized by different characters of the energy dependences , namely : ( i ) the phase d by slowly descending dependences of both @xmath330 and @xmath331 averages , ( ii ) the phase tc by roughly constant dependences , ( iii ) the phase n by linearly increasing dependences , and ( iv ) the phase s by saturated @xmath334 and linearly increasing @xmath331 . hamiltonian ( [ h ] ) with @xmath97 , @xmath66 and @xmath69 at @xmath214 . quantum phases d , tc , n and s are distinguished by different smoothed energy dependences of the average slopes  see the piecewise fits indicated by full black lines . ] peres lattices for both quantities @xmath52 and @xmath322 are connected with the lattice of the interaction hamiltonian @xmath335 via the energy conservation relation . using the hellmann - feynman formula , @xmath336 , one derives the following expression for the slope of individual energy levels : @xmath337 the slopes averaged over bunches of neighboring 20 levels , as in the insets of fig.[f_pere ] , are presented in fig.[f_slo ] for the same values of control parameters . piecewise fits of the energy dependences in individual quantum phases are again indicated by lines . we see that any of the d , tc , n and s quantum phases is characterized by a specific , roughly invariant energy dependence of the averaged level slope @xmath338 within the corresponding energy domain . abrupt changes of the character of these dependences coincide with critical esqpt energies . this is in accord with a general relation between the density and flow signatures of esqpts , namely with the fact that a typical esqpt generates the same type of non - analyticity in the energy dependences of both quantities @xmath316 and @xmath338 @xcite . indeed , the discontinuities of @xmath339 at the critical energies @xmath174 and @xmath175 , as observed in fig.[f_slo ] , are consistent with the analogous behavior of the @xmath300 ( fig.[f_dender ] ) . on the other hand , the anticipated point of a singular growth of the level slope ( logarithmic divergence of its energy derivative ) at @xmath173 is smoothed out in the finite spectrum for a moderate system s size . we turn to the study of quantum entanglement properties of individual eigenstate of the model hamiltonian and their links to the esqpts . our analysis includes two types of bipartite entanglement : ( a ) that between the bosonic field and the set of all atoms ( atom  field entanglement ) , and ( b ) that between any pair of individual atoms ( atom  atom entanglement ) . the entanglement of type ( a ) is an important ingredient of superradiance since the interaction term of the hamiltonian carries a direct coupling between the atomic and field subsystems . in contrast , the entanglement of type ( b ) appears only due to an indirect coupling of individual atoms via the bosonic field , so it may be expected to be just a higher - order  effect . we start with a brief description of the measures used to quantify both types of entanglement .    in the following , the atom  field and atom  atom entanglement will be evaluated in individual eigenstates @xmath340 of hamiltonian . some examples of the eigenstate wave functions were shown in fig.[f_wav ] . since the wave - function arguments @xmath99 and @xmath341 , respectively , correspond directly to the atomic and field coordinates , a compound state @xmath342 is factorized with respect to the atom  field partitioning of the system if it has a product wave function @xmath343 , where @xmath344 and @xmath345 are arbitrary atomic and field wave functions . the method to quantify a departure of a given pure state @xmath342 from exact factorization , i.e. , an amount of atom  field entanglement involved in @xmath342 , makes use of the von neumann entropy corresponding to the reduced density operators @xmath346 and @xmath347 of the atomic and field subsystem , respectively . these operators are obtained by partial tracing of the total density operator @xmath348 over the irrelevant part of the total hilbert space , that is @xmath349 ( trace goes over @xmath350 ) for the atomic subsystem and @xmath351 ( trace goes over @xmath352 ) for the field subsystem . von neumann entropy of the pure compound state @xmath353 is by definition zero , but entropies of both reduced density operators @xmath346 and @xmath347 satisfy : @xmath354 . the case @xmath355 implies a separable compound state , while the case @xmath356 indicates that the reduction of @xmath353 to a single subsystem leads to a loss of information on mutual entanglement between subsystems . a maximally entangled state yields @xmath357 , where @xmath31 is a dimension of @xmath358 , the smaller of both subspaces . we therefore define a normalized atom  field entanglement entropy in a compound state @xmath342 as : @xmath359}{\ln(2j+1 ) } = -\frac{\textnormal{tr}\left[\rho_{\textnormal{f}}\,\textnormal{ln}\rho_{\textnormal{f}}\right]}{\ln(2j+1 ) } \label{s } \,.\ ] ] it changes between @xmath360 for separable states and @xmath361 for maximally entangled states .    quantifying the atom  atom entanglement , i.e. , quantum correlations between a randomly chosen pair @xmath362 of atoms , is a more complex problem . the use of the above entropic approach is disabled by the fact that for any state @xmath342 of the whole atom  field system , an arbitrary pair of atoms generically occurs in a _ mixed _ quantum state . it is known that mutual entanglement of a pair of objects in a mixed compound state can not be recognized by a non - zero entropy of the reduced density operators @xcite . indeed , a mixed compound state of the pair generates mixed reduced states of the objects even in absence of entanglement . in such cases , the evaluation of an entropy - based measure of entanglement ( so - called entanglement of formation ) has to be performed with respect to all possible decompositions of the compound density operator into statistical mixtures of pure states @xcite , which is from the computational viewpoint a difficult task .    a way to bypass this obstacle for two - qubit systems was proposed in refs.@xcite in terms of a quantity called concurrence . the idea was applied and further elaborated @xcite for multi - qubit systems in fully symmetric states , like our ensemble of @xmath1 two - level atoms with @xmath160 . the atom  atom entanglement in this case is characterized by a scaled concurrence : @xmath363 where @xmath364 are eigenvalues ( real and non - negative ) of a non - hermitian matrix @xmath365 , with @xmath366 denoting a reduced density matrix of the pair of atoms @xmath362 and @xmath367 the corresponding spin - flipped  ( conjugate under the time reversal ) state ( with @xmath368 denoting the tensor product of operators for the selected atoms and @xmath369 a complex conjugated density matrix in the @xmath370 eigenbasis ) . the argument @xmath371 in the definition again reminds the compound state @xmath342 in which the concurrence is calculated . it has to be stressed that the full symmetry of @xmath342 under the exchange of atoms , which is guaranteed in the @xmath160 subspace , ensures that the reduced density matrix @xmath366 , as well as @xmath372 and @xmath373 , are the same for _ any _ pair of atoms , hence do not in fact depend on @xmath374 . model in the plane of control parameters @xmath18 and @xmath64 ( for @xmath66 , @xmath375 ) . the red curve indicates the qpt critical coupling @xmath279 . left : atom  field entanglement measured by the entropy . right : atom  atom entanglement measured by the concurrence . the darker areas indicated larger entanglement . ] matrix elements of the reduced two - atom density matrix @xmath366 for any symmetric state of an @xmath1-qubit system can be expressed through expectation values @xmath376 , @xmath377 and @xmath378 of collective quasi - spin operators in a given state  see formulas ( 4 ) and ( 11 ) in ref.@xcite . this makes it possible to perform a straightforward numerical computation of @xmath373 in the @xmath160 subspace of our model . moreover , as shown in ref.@xcite , the scaled concurrence in this subspace is related to the so - called kitagawa  ueda spin - squeezing parameter @xmath379 though the relation @xmath380 , so it varies between @xmath381 for separable states and @xmath382 for the states that exhibit the maximal entanglement allowed for a given @xmath1 . an entropy corresponding to the entanglement of formation of a single atomic pair reads as @xcite : @xmath383 \label{saa } \,.\ ] ] for @xmath381 we have @xmath384 , while for @xmath382 we obtain @xmath385 , which decreases from @xmath386 for @xmath387 to zero in the classical limit @xmath388 . with @xmath389 , the maximal entropy behaves as @xmath390 . note that the @xmath389 maximal atom  field entanglement entropy _ per atom _ decreases as @xmath391 , see eq . , so it exceeds the atom  atom entanglement entropy about @xmath1 times . these scaling properties are verified theoretically and by large-@xmath1 numerical calculations @xcite . the atom  field and atom  atom entanglement properties have been studied for the ground state of the dicke model @xcite . as in other models of similar nature , see e.g. refs.@xcite , the second - order qpt was shown to induce a singularity in both entanglement measures and . figure  [ f_eng ] verifies this conclusion in our extended model with variable parameter @xmath64 . we observe that both @xmath392 and @xmath373 measures exhibit an increase at about the critical coupling @xmath279 from eq .. for @xmath184 , the atom  atom entanglement drops back to nearly zero values with @xmath393 . the atom  field entanglement saturates at a value @xmath394 for @xmath395 , which is due to an irreducible atom  field coupling in the lowest positive parity state in the strong coupling limit @xcite . in the @xmath396 region , close to the integrable tavis - cummings limit , both entanglement measures show roughly a step - like increase at @xmath279 . this is due to a specific mechanism , in which the ground state at each @xmath18 is formed via unavoided crossings of levels with different values of the conserved quantum number @xmath108 , see fig.[f_spe]a . we will analyze the integrable case in more detail below .      in the following , we present results of a numerical study of the atom  field entanglement in individual eigenstates of hamiltonian with @xmath37 . we start by analyzing the integrable tavis - cummings limit @xmath63 . this simple setting will allow us to obtain some insight into the entanglement properties from quasi - analytic solutions , which will serve as a useful reference for the less trivial @xmath184 case . we know that for @xmath63 the hilbert space splits into the subspaces @xmath101 with fixed values of quantum number @xmath108 , see eq . , with the corresponding dimensions @xmath126 given by . the reduced density matrix of the atomic subsystem within each @xmath101 reads as @xmath397 where @xmath398 is a number of bosons associated with quasi - spin projection @xmath51 for a given value of @xmath108 . this expression implies that for @xmath63 the entanglement entropy of any eigenstate is equal to the wave - function entropy corresponding to its expansion in the non - interacting basis .     model with @xmath66 for @xmath296 and @xmath97 ( left panel ) , and the atom  field entanglement entropies @xmath392 in individual eigenstates corresponding to @xmath399 cut of the spectrum ( right panel ) . the lowest v - shaped chain of points ( connected by the green curve ) in the entropic spectrum corresponds to the lowest states from various @xmath108-subspaces , the second chain to the second states etc . the points corresponding to two selected @xmath108-subspaces ( @xmath400 and 8) are highlighted . ] let us focus first on the tuned case @xmath116 . the wave - function entropies for the @xmath301 spectrum of a tuned @xmath63 hamiltonian were shown in panel ( a ) of fig.[f_cary ] . we notice that the entropy of individual eigenstates does not change with @xmath18 and that it has an apparent symmetry with respect to the vertical reflection of the spectrum around the centroid energy @xmath118 . the former feature was explained by constancy of eigenstates of the simple hamiltonian , the latter follows from a recursive relation for the eigenstate components in eq . which yields the same distributions @xmath401 , hence the same entropies , for the pair of levels with opposite slopes . the full spectrum of the tavis - cummings model is obtained by combining the spectra for all values of @xmath108 . this is for @xmath116 seen in fig.[f_pileup ] . it shows the energy spectrum as a function of @xmath18 ( in the left panel ) and the atom  field entanglement entropy of individual eigenstates for one particular value of the coupling strength ( the right panel ) . as follows from the above theoretical analysis , the energy spectrum in the left panel of fig.[f_pileup ] is formed by straight lines . the qpt critical coupling @xmath402 coincides with the point where the energy of the @xmath312 ground state is crossed by the lower state from the @xmath313 space . a further increase of @xmath18 above @xmath402 leads to a sequence of consecutive level crossings in which the lowest states from subspaces with increasing @xmath108 become instantaneous ground states of the system . a similar mechanism applies also in the spectrum of excited states , where we observe a sequence of separated caustic structures formed by states with increasing ordinal numbers within each @xmath108-subspace ( the @xmath322th caustic structure in the vertical direction represents an envelope of lines corresponding to the @xmath322th states from individual @xmath108-subspaces ) . the entropic spectrum on the right of fig.[f_pileup ] arises in a similar way . it is composed of several mutually shifted v - shaped chains of points , each of them containing a collection of states from different @xmath108-subspaces . for example , the lowest v - chain of points is formed by the lowest - energy states from all @xmath108-subspaces , the second lowest chain by second lowest - energy states and so on . to indicate a contribution of a single @xmath108-subspace to the entropic spectrum , we highlighted the points corresponding to two values of @xmath108 . for @xmath403 , the contribution has a seagull - like form ( see the @xmath404 example in the figure ) , whose reflection symmetry around @xmath118 results from the above - discussed properties of eq .. the instantaneous ground state at the given @xmath405 , which is well above the critical point @xmath402 , has a high value of @xmath108 and therefore carries a considerable amount atom  field entanglement ( as in general expected in superradiance ) . on the other hand , the trivial one - dimensional subspace with @xmath312 formed by the state @xmath406 , which was the ground state of the @xmath63 hamiltonian at @xmath407 and becomes an excited state with energy @xmath408 for @xmath409 , has apparently zero atom  field entanglement . we will see below that some remnants of this state remain present in entropic spectra also in the non - integrable @xmath184 regime . model with @xmath296 and @xmath97 in the detuned case ( @xmath113 ) for several values of the coupling strength @xmath18 . ] let us consider now the detuned case @xmath132 . field entropic spectra of a hamiltonian with @xmath63 and @xmath410 are displayed in fig.[f_detun ] for increasing values of the coupling strength @xmath18 ( the critical coupling is @xmath411 ) . we know ( see sec.[ha ] ) that in the strong coupling limit , @xmath412 , the detuned situation becomes very similar to the tuned one . this is obvious from a comparison of the @xmath413 entropic spectrum in fig.[f_detun ] with that in fig.[f_pileup ] . on the other hand , for small coupling strengths , @xmath414 , the detuning @xmath415 plays an important role . a detuned system shows generally lower values of the atom  field entanglement than the tuned one , which can be attributed to a smaller perturbation efficiency of the interaction hamiltonian @xmath310 due to larger energy gaps between the unperturbed ( factorized ) states ( the energy spread of a single-@xmath108 spectrum at @xmath320 increases with @xmath415 ) . this overall lowering of entanglement entropies at small @xmath18 ( more apparent for @xmath133 than @xmath134 ) can be verified by simple numerical calculations , which we do not present here . it is clear that the factorized state @xmath406 with @xmath312 at critical energy @xmath416 is an invariant eigenstates of any @xmath63 hamiltonian , irrespective of @xmath18 and @xmath0 s . so the @xmath360 point at the place of the unperturbed ground state is present in all entropic spectra of fig.[f_detun ] . however , for an @xmath132 system , a local decrease ( converging to @xmath360 with @xmath150 ) of the atom  field entanglement is present also at the highest critical energy @xmath417 . this is due to quantum critical properties of the @xmath301 subspace discussed in sec.[qpt ] . we consider here the @xmath133 case ( see fig.[f_cary]b ) , but the same phenomenon takes place also in the reversed detuning hierarchy @xmath134 ( fig.[f_cary]c ) . for @xmath418 , where @xmath419 is the critical coupling , the lowest state of the @xmath301 subspace at energy @xmath175 keeps the factorized form @xmath315 ( mind the sign difference in @xmath51 from the other factorized state , which means that now a maximal possible , for a given @xmath25 , number of atoms is in the _ upper _ level ) . distinct lowering of the atom  field entanglement entropy at @xmath317 is observed in the @xmath420 and 0.4 panels of fig.[f_detun ] ( in which @xmath421 ) . this qpt - related effect becomes sharper with increasing @xmath1 and @xmath25 , and in the @xmath422 limit it results in a similar @xmath360 cusp as that at energy @xmath174 . for @xmath134 , the above factorized state , causing the same effect , is the highest state of the @xmath301 subspace at @xmath320 . subspace of the @xmath68 detuned ( @xmath113 ) model with @xmath296 : the ground state ( panel a ) and excited states which at @xmath423 have 15% ( panel  b ) and 24% ( panel  c ) of the maximal excitation energy . the dashed curves are for @xmath97 , the full ones for @xmath424 . the energy of the corresponding level ( for @xmath97 ) is shown at the bottom of each panel . the log - log plot in panel ( b ) inset shows the minimum value of @xmath392 in the dip as a function of @xmath25 for the level at 15% of the spectrum ( for this level , the dependence is roughly @xmath425 ) . ] an analogous lowering of the atom  field entanglement entropy at the upper critical energy @xmath175 appears also for the coupling strengths _ above _ the qpt critical value @xmath419 . the effect in this domain is caused by the @xmath170 type of esqpt in the @xmath301 subspace . indeed , when excited states within this subspace cross the critical energy , they get temporarily localized in the coordinate region around the point @xmath426 , i.e. , at the position of the factorized state . this is due to dwelling  of the @xmath317 classical trajectory and a semiclassical wave function at an inflection point of the classical hamiltonian . an analogous phenomenon was observed also in other quasi @xmath170 models , see e.g. ref.@xcite . the lowering of the eigenstate entropy in the @xmath427 domain was visible already in panels ( b ) and ( c ) of fig.[f_cary ] , but it is not apparent in fig.[f_detun ] . to see the effect in the pile up of all-@xmath108 spectra , one would have to go to higher values of @xmath1 and @xmath25 , when the localization becomes stronger and the entropy decreases to its lower limit @xmath360 . entropic spectra of the atom  field entanglement with increasing parameter @xmath60 $ ] at @xmath214 , @xmath428 , @xmath97 . the esqpt energies are indicated by vertical lines . the spectrum contains @xmath429 well converged eigenstates . ] the esqpt - based localization effect is nevertheless well documented in fig.[f_entrotc ] , which depicts the evolution of the atom  field entanglement for three selected levels ( ground state and two excited states ) from the @xmath430 subspace for @xmath431 and 1000 . we see that the passage of the selected excited states through the esqpt energy @xmath175 is correlated with a local decrease of the entanglement entropy [ sharp dips in the dependences in panels ( b ) and ( c ) ] . as demonstrated by the curves for two values of @xmath1 , the drop of entropy becomes sharper and deeper with increasing size of the system . the decrease of @xmath392 at the lowest point of the dip with @xmath37 for one of the excited states is shown in the inset of fig.[f_entrotc]b . the depicted log - log plot indicates a very slow power - low decrease according to @xmath432 . the exponent @xmath433 is rather small and depends on the selected level ( e.g. , for the displayed level at 15% of the spectrum we obtain @xmath434 , while for the level at 20% we would have @xmath435 ) note that the entanglement entropy for excited states for large @xmath18 saturates at a value close to @xmath436 , where @xmath437 is the dimension for @xmath301 . this expression follows from the use of the random matrix theory for estimating an average wave - function entropy for @xmath126 orthonormal states in the limit of strong mixing @xcite . note also that the critical coupling @xmath419 of the ground - state qpt coincides with the sharpest increase of the ground - state entanglement entropy ( fig.[f_entrotc]a ) . now we are ready to turn our attention to the non - integrable @xmath184 model . figure  [ f_sdelta ] shows the evolution of the entropic spectrum with parameter @xmath171 $ ] for @xmath116 and @xmath438 , well above the critical coupling @xmath402 for any @xmath64 . the vertical lines demarcate the esqpt energies , the zones between them corresponding to the d , tc , n and s quantum phases . one can see that the departure from integrability ( an increase of @xmath64 ) gradually destroys regular patterns in the entropic spectra , in analogy with peres lattices . orderly patterns in the entropic spectrum occur in the domains of regular dynamics , while disorderly scattered points indicate chaotic domains . note that the accumulation of points in a narrow band close to ( slightly below ) the maximal entropy @xmath439 for high energies is a consequence of ergodicity of chaotic dynamics @xcite : a common eigenstate in the chaotic part of the spectrum contains virtually all atom  field configurations , whose dominant part is strongly entangled . on the other hand , the regular low - energy part of the entropic spectrum tends to a minimal value @xmath440 , corresponding to states with a good parity quantum number @xcite .    , @xmath69 ( cf.fig.[f_wav]e ) and ( b ) the factorized @xmath441 eigenstate for @xmath68 . the other parameters are the same as in fig.[f_sdelta ] . ] an interesting aspect of the entropic spectra in fig.[f_sdelta ] is the evolution of the atom  field entanglement in transition between the superradiant and normal phases at the critical energy @xmath318 . for @xmath184 , this energy corresponds to @xmath173 for @xmath442 $ ] and to @xmath174 for @xmath443 . we know that in the @xmath63 limit , there always exists the fully factorized trivial eigenstate with @xmath444 located right at this energy and a group of weakly entangled eigenstates with low values of @xmath108 located nearby ( see figs.[f_pileup ] and [ f_detun ] ) . to what extent do these structures survive an increase of parameter @xmath64 ? the answer can be read out from fig.[f_sdelta ] . we see that the eigenstates with lowered entropy and even a single state with @xmath445 are preserved in the entropic spectra as far as the tc phase is present for the chosen value of the coupling strength @xmath18 . this is so if @xmath446 , hence @xmath447 , see eq .. if @xmath448 $ ] , that is for @xmath64 large enough to avoid the existence of the tc phase for a given @xmath18 , the decrease of entropy no longer takes place . these observations allow us to say that the occurrence of states with lowered atom  field entanglement entropy is a signature of the esqpt between the n and tc phases , but not of that between the n and d phases . the wave function corresponding to the @xmath449 eigenstate with the lowest entanglement entropy for @xmath94 can be seen in panel ( e ) of fig.[f_wav ] . its detail is depicted in panel ( a ) of fig.[f_wavdet ] , in comparison with the fully factorized state with @xmath312 shown in panel ( b ) . although the fine structures of the wave function for @xmath184 ( panel a ) apparently disable its full factorization ( in contrast to the @xmath63 case ) , a distant view in fig.[f_wav]e makes clear that this state is indeed very similar to the @xmath312 state . following our findings for @xmath63 , one might expect an analogous decrease of the atom  field entanglement entropy also at the critical energy @xmath175 ( in transition between the s and n phases ) in a detuned ( @xmath132 ) system with @xmath184 . a numerical evidence of this phenomenon is however hindered by the above - presented ( cf . the inset in fig.[f_entrotc]b ) slowness of the decrease of the entanglement entropy dip to zero . to see an effect in the @xmath184 system , one would have to perform a truncated - space diagonalization of the hamiltonian for a very large value of @xmath1 , which is computationally demanding . model with @xmath66 , @xmath214 , @xmath296 , @xmath97 . the contribution of selected @xmath108-subspaces ( the u - shaped chains of points ) is highlighted . ]     for @xmath214 , @xmath66 , @xmath296 , @xmath97 . the concurrence spectra contain 5000 well converged states . ] we turn now to a brief numerical analysis of the atom  atom entanglement , which is limited to the symmetric @xmath37 subspace of the atomic hilbert space . as in the previous case , we focus at first on the integrable tavis - cummings model . the scaled concurrence for individual eigenstates ( the concurrence spectrum ) of the tuned hamiltonian with @xmath63 and @xmath409 is plotted in fig.[f_concint ] . the highlighted points denote states belonging to selected @xmath108-subspaces . like in the entropic spectrum ( fig.[f_pileup ] ) , the states corresponding the same @xmath108-subspace form a pattern which is always symmetric with respect to the reflection around the central energy @xmath118 . trivially , the @xmath312 state @xmath406 has zero concurrence as the atomic state with minimal quasi - spin projection is fully factorized ( all atoms in the lower level ) . so this state shows no entanglement in either atom  field or atom  atom sense . in contrast , some states with low values of @xmath108 , which all have only weak ( increasing with @xmath108 ) atom  field entanglement , show relatively large ( decreasing with @xmath108 ) atom  atom entanglement . consider as an example a doublet of states with @xmath313 . the scaled concurrence for these states is @xmath450 , which changes between @xmath451 for @xmath387 and the maximal value @xmath382 for @xmath192 ( see the upper pair of points in fig.[f_concint ] ) . as follows from eq . , the @xmath313 eigenstates are expressed by the superpositions @xmath452 , whose first term is fully factorized in the atomic subspace ( as in the @xmath312 case ) while the second term is maximally entangled ( a symmetric state with one atom in the upper level and the rest of atoms in the lower level ) . it turns out that for increasing @xmath453 , only the states with the largest positive and negative slopes within the given @xmath108-subspace yield a non - vanishing atom  atom entanglement . these states form the feelers  of the u - shaped dependences of @xmath373 in each @xmath108-subspace . the scaled concurrence of these states decreases with @xmath108 , forming together the left and right chains of @xmath454 points in fig.[f_concint ] . the rest of the @xmath63 eigenstates has @xmath381 . figure  [ f_cdelta ] shows an evolution of the concurrence spectrum with increasing @xmath64 for a fixed @xmath455 . apparently , an overall trend of the atom  atom entanglement in individual eigenstates is a decrease with increasing @xmath64 . we stress that the concurrence of an absolute majority of states for any @xmath64 is zero ( the spectrum contains 5000 states ) and both the number of states with @xmath454 and the values of @xmath373 in these states further decrease with @xmath64 . the only states in which the atom  atom entanglement remains significant even for a relatively large values of @xmath64 appear in a narrow energy interval around the esqpt critical energy @xmath174 demarcating the transition between the tc and n phases . these states partly coincide with those in which we previously observed a reduced atom  field entanglement , see fig.[f_sdelta ] . for @xmath64 large enough to avoid the existence of the tc phase , the atom  atom entanglement of all states vanishes . model with @xmath456 for two values of @xmath64 . the number of principal components of an actual eigenstate in the unperturbed basis is encoded in the shades of gray of the respective line ( black means perfect localization , @xmath457 , white complete delocalization , @xmath458 ) . revival of some @xmath459 eigenstates at @xmath460 is seen in both spectra . ] we analyzed properties of a generalized dicke model of single - mode superradiance allowing for a continuous ( governed by parameter @xmath64 ) crossover between the integrable tavis - cummings limit and the prevailingly chaotic dicke limit . we considered three constraints on the model hilbert space @xmath461 , namely : ( i ) the all-@xmath25 model with the full @xmath462 containing the complete @xmath26-dimensional atomic hilbert space @xmath35 , ( ii ) a single-@xmath25 model whose atomic hilbert space is restricted to a subspace @xmath463 with fixed quasi - spin quantum number @xmath25 ( usually @xmath37 ) and ( iii ) a single-@xmath464 model with only a single subspace @xmath465 of quantity @xmath108 , which is conserved in the integrable limit @xmath63 . we determined the classical limit and critical properties of the model . the all-@xmath25 model has @xmath42 degrees of freedom , the single-@xmath25 model @xmath40 and the single-@xmath464 model @xmath170 . the thermal phase diagram of the all-@xmath25 model for @xmath197 shows a coexistence of the tc and d types of superradiant phase ( with and without saddles in the free energy landscape ) , which contract to a single type tc or d in the limits @xmath63 or 1 , respectively . the quantum phase diagram of a single-@xmath25 model contains the ground - state superradiant quantum phase transition and three types of excited - state quantum phase transitions , characterized by singularities ( jumps and a logarithmic divergence ) in the first derivative of the semiclassical level density as a function of energy . the quantum phase diagram of a critical single-@xmath464 model with @xmath301 shows in the detuned case ( @xmath132 ) an excited - state quantum phase transition characterized by a singularity ( logarithmic divergence ) directly in the level density . we defined the spectral domains in between the esqpt critical borderlines of the single-@xmath25 model as the tc , d , n and s quantum phases . these phases were characterized by different shapes of a smoothed energy dependence of some expectation values , which show abrupt changes at the esqpt energies . a natural choice of the phase - defining quantity is the interaction hamiltonian , whose expectation value @xmath466 determines the slope @xmath467 of the given level @xmath326 . according to the density  flow relation of the esqpt signatures @xcite , a smoothed level slope should exhibit the same type of non - analyticity as the semiclassical level density . therefore , quantum phases of systems with @xmath40 ( like our single-@xmath25 model ) are distinguished by @xmath468 , while phases in other systems are expected to differ in the @xmath469th energy derivative of @xmath470 . we analyzed entanglement properties of the single-@xmath25 and single-@xmath464 models . in general , an absolute majority of eigenstates in the spectrum for any choice of model parameters has large atom  field entanglement but vanishing atom  atom entanglement . exceptional in this sense for @xmath63 are the eigenstates with low values of @xmath108 , which yield a weak atom  field entanglement and simultaneously an increased atom  atom entanglement . we showed that remnants of these states exist relatively far in the non - integrable regime with @xmath184 . this concerns also the @xmath63 state with @xmath312 located at the esqpt energy @xmath174 , which has zero atom  field entanglement . for @xmath132 , a state with reduced atom  field entanglement entropy , originating in the critical @xmath301 subspace , exists also at the esqpt energy @xmath175 . a numerical evidence of its partial survival for @xmath184 is however difficult because of the observed slow decrease of the entropy dip with the size parameter @xmath1 . it needs to be stressed that the above states with exceptional entanglement properties form only a small subset of all states ( its relative size goes to zero in the @xmath388 limit ) . even so is their occurrence a significant phenomenon . the remnants of the @xmath63 states with low @xmath108 s appear in the range of parameter @xmath64 where the tc quantum phase still exists in the spectrum . so if @xmath64 is fixed at a non - zero value and the coupling strength @xmath18 increases from 0 to values above @xmath402 and @xmath471 , these states exhibit a surprising evolution . they first ( not later than at @xmath402 ) disappear from the spectrum , but later above the critical coupling ( somewhere in the @xmath472 region ) some of them reappear again . indeed , a systematic analysis of the spectrum discloses approximate revivals of a number of weak - coupling ( @xmath414 ) eigenstates in the strong - coupling ( @xmath473 ) regime . this is illustrated in fig.[f_fog ] where we again show spectra of the tuned model with @xmath37 for two values of @xmath64 . the fullness ( shade of gray ) of each level encodes the number of principal components @xmath474 of the corresponding hamiltonian eigenstate @xmath475 in the eigenbasis @xmath476 taken at @xmath477 infinitesimally displaced from the degeneracy point @xmath320 . the number of principal components varies from @xmath478 for a perfectly localized state ( identical with one of the @xmath479 eigenstates ) to @xmath480 for totally delocalized states . the revival of several @xmath479 eigenstates is clearly observed in fig.[f_fog ] as reappearance of certain dark lines in the spectrum for large @xmath18 in the tc and n phases . one can only speculate that the mechanism responsible for this surprising phenomenon can be perhaps related to a specific partial dynamical symmetry @xcite . the authors thank to t.brandes , n.lambert , j.hirsch and m.a.bastarrachea-magnani for illuminating discussions . this work was supported by the czech science foundation under the project no.p203 - 13 - 07117s . dicke , phys . rev . * 93 * , 99 ( 1954 ) . m. gross and s. haroche , phys . rep . * 93 * , 301 ( 1982 ) . benedict ( ed . ) , _ super - radiance : multiatomic coherent emission _ , taylor & francis , new york , 1996 . t. brandes , phys . rep . * 408 * , 315 ( 2005 ) . bekenstein and m. schiffer , phys . d * 58 * , 064014 ( 1998 ) . n. auerbach and v. zelevinsky , _ rev . phys . _ * 74 * 106301 ( 2011 ) . jaynes and f.w . cummings , proc . ieee * 51 * , 89 ( 1963 ) . m. tavis and f.w . cummings , phys . rev . * 170 * , 379 ( 1968 ) . wang and f.t . hioe , phys . rev . a * 7 * , 831 ( 1973 ) . k. hepp and e.h . lieb , phys . rev . a * 8 * , 2517 ( 1973 ) . k. hepp and e.h . lieb , ann . * 76 * , 360 ( 1973 ) . k. rzazewski , k. wodkiewicz and w. zakowicz , phys . 35 * , 432 ( 1975 ) . k. baumann , c. guerlin , f. brennecke and t. esslinger , nature * 464 * , 1301 ( 2010 ) k. baumann , r. mottl , f. brennecke and t. esslinger , phys . lett . * 107 * , 140402 ( 2011 ) p. prez - fernndez , p. cejnar , j.m . arias , j. dukelsky , j. e. garca - ramos and a. relao , phys . rev . a * 83 * , 033802 ( 2011 ) . p. prez - fernndez , a. relao , j.m . arias , p. cejnar , j. dukelsky and j. e. garca - ramos , phys . e * 83 * , 046208 ( 2011 ) . t. brandes , phys . e * 88 * , 032133 ( 2013 ) . bastarrachea - magnani , s. lerma - hernndez and j.g . hirsch , phys . rev . a * 89 * , 032101 , 032102 ( 2014 ) . p. cejnar , m. macek , s. heinze , j. jolie and j. dobe , j. phys . a * 39 * , l515 ( 2006 ) . caprio , p. cejnar and f. iachello , ann . ( n.y . ) * 323*,1106 ( 2008 ) . p. cejnar and p. strnsk , phys . e * 78 * , 031130 ( 2008 ) . p. strnsk , m. macek and p. cejnar , ann . phys . * 345 * , 73 ( 2014 ) . p. strnsk , m. macek , a. leviatan and p. cejnar , ann . phys . * 356 * , 57 ( 2015 ) . p. strnsk and p. cejnar , phys . a * 380 * , 2637 ( 2016 ) . p. cejnar and p. strnsk , phys . * 91 * ( 2016 ) 083006 . a. osterloh , l. amico , g. falci and r. fazio , nature * 416 * , 608 ( 2002 ) . j. vidal , g. palacois and r. mosseri , phys . a * 69 * , 022107 ( 2004 ) . j. vidal , r. mosseri and j. dukelsky , phys . a * 69 * , 054101 ( 2004 ) . n. lambert , c. emary and t. brandes , phys . lett . * 92 * , 073602 ( 2004 ) . n. lambert , c. emary , and t. brandes , phys . a * 71 * , 053804 ( 2005 ) . p. cejnar and j. jolie , phys . e * 58 * , 387 ( 1998 ) . p. cejnar and p. strnsk , aip conf . proc . * 1575 * , 23 ( 2014 ) . m.a . bastarrachea - magnani , b. lpez - del - carpio , s. lerma - hernndez and j.g . hirsch , phys . a * 90 * , 068015 ( 2015 ) . bastarrachea - magnani , b. lpez - del - carpio , j. chvez - carlos , s. lerma - hernndez and j.g . hirsch , phys . e * 93 * , 022215 ( 2016 ) . j. chvez - carlos , m.a . bastarrachea - magnani , s. lerma - hernndez and j.g . hirsch ( to be published ) ; see arxiv:1604.00725 [ nlin.cd ] ( 2016 ) . bastarrachea - magnani , s. lerma - hernndez and j.g . hirsch ( to be published ) ; see arxiv:1605.05357 [ quant - ph ] ( 2016 ) . a. peres , phys . rev . lett . * 53 * , 1711 ( 1984 ) . p. strnsk , p. hruka and p. cejnar , phys . e * 79 * , 066201 ( 2009 ) . a. peres , phys . * 77 * , 1413 ( 1996 ) . bennett , d. divincenzo , j. smolin and w.k . wooters , phys . a * 53 * , 3824 ( 1996 ) . s. hill and w.k . wootters , phys . * 78 * ( 1997 ) 5022 . wootters , phys . * 80 * ( 1998 ) 2245 . x. wang and k. mlmer , eur . j. d * 18 * 385 ( 2002 ) , x. wang and b.c . sanders , phys . a * 68 * , 012101 ( 2003 ) . s. heinze , p. cejnar , j. jolie and m. macek , phys . c * 73 * 014306 ( 2006 ) . y. alhassid and a. leviatan , j. phys . a : math . gen . * 25 * , l1265 ( 1992 ) . a. leviatan , prog . . phys . * 66 * , 93 ( 2011 ) .
Abstract:

Answer: we study a simple model describing superradiance in a system of two - level atoms interacting with a single - mode bosonic field . the model permits a continuous crossover between integrable and partially chaotic regimes and shows a complex thermodynamic and quantum phase structure . several types of excited - state quantum phase transitions separate quantum phases that are characterized by specific energy dependences of various observables and by different atom  field and atom  atom entanglement properties . we observe an approximate revival of some states from the weak atom field coupling limit in the strong coupling regime .    single - mode superradiance model , excited - state quantum phase transitions , thermal and quantum phases , entanglement properties of excited states