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: quantum hall effects are characteristic of two - dimensional ( 2d ) electron gases in strong magnetic fields . the integral quantum hall effect ( iqhe ) @xcite is explained in terms of weakly - interacting electrons filling quantized landau levels ( ll ) produced by application of a strong magnetic field perpendicular to the 2d gas @xcite . in contrast , the fractional quantum hall effect ( fqhe ) @xcite can occur only as a result of strong electronic correlations in partially - filled landau levels @xcite . because of its massless chiral charge carriers and atomic - monolayer honeycomb lattice exhibiting sublattice , valley , and spin degeneracies , quantum hall effects in graphene could be much richer than in the semiconductor 2d electron gas , where there is no crystal structure and the only degeneracies are those of the ( non - relativistic ) landau levels and spin . for graphene in strong magnetic fields an integral quantum hall effect @xcite and a fractional quantum hall effect @xcite have been observed , with anomalous filling factors that reflect the unique degeneracies of the graphene electronic structure and the dirac - like nature of its electrons . the valley isospin and spin degrees of freedom imply that graphene in a strong magnetic field is described well by a low - energy hamiltonian that commutes approximately with the generators of an su(4 ) lie algebra . this su(4 ) symmetry allows classification of states in graphene , and can serve as the basis for computing explicit breaking of the su(4 ) symmetry by small non - symmetric terms in the effective hamiltonian . however , there is growing evidence that many states observed in modern experiments can not be described by explicit breaking of su(4 ) . for example , the ground state of graphene in a magnetic field at low temperature exhibits a rapid increase of the longitudinal resistance @xmath0 above a critical magnetic field @xmath1 @xcite . the value of @xmath1 decreases for increasingly cleaner samples , indicating that the resistance is not caused by impurity scattering but instead is intrinsic to the state itself @xcite . in quantum hall systems the currents are carried by edge states , so this insulating ground state must correspond to an emergent state that does not support edge currents produced by _ spontaneous _ ( not explicit ) breaking of the su(4 ) symmetry . thus the approximate su(4 ) symmetry of graphene can suggest possible low - energy collective modes exhibiting spontaneously broken symmetry that are important for the properties of graphene in a magnetic field , but the su(4 ) symmetry alone can not determine which of these modes is the ground state . until now , those broken - symmetry modes have been addressed quantitatively by numerical simulations employing limited numbers of states and particles , or by effective low - energy field - theory approximations . these calculations find various possible low - energy collective states resulting from spontaneous breaking of the su(4 ) symmetry with very similar energies . thus they have been unable to give a definitive answer to the nature of the insulating ground state . let us note that other approaches to the problems addressed here have been proposed ( for example , refs . @xcite ) . however , the present discussion will concentrate on methods based on approximate su(4 ) symmetry of the hamiltonian . an alternative and potentially more powerful application of symmetries has been employed extensively in both nuclear structure and condensed matter physics @xcite . fermion dynamical symmetry method _ truncates the hilbert space to a tractable collective subspace by positing a highest symmetry associated with the physical operators for the system , and constructing effective hamiltonians from polynomials in the casimir invariants of the highest symmetry s subgroup chains . in this approach it is possible not only to classify low - energy collective modes , but to solve analytically for the properties of these modes and to determine which lie lowest in energy , either exactly in particular symmetry limits , or approximately using generalized coherent state methods .    in an earlier letter the first application of fermion dynamical symmetry methods to graphene was introduced and applied to determining the ground state in strong magnetic fields @xcite . this paper develops the full dynamical symmetry formalism upon which ref . @xcite rests . it will be shown that the highest symmetry is so(8 ) , with its generators identifiable with particle  hole and pairing degrees of freedom that have been discussed previously in the physics of graphene . this symmetry will be shown to be isomorphic to an so(8 ) symmetry used extensively in nuclear structure physics , which permits already - developed mathematics to be appropriated for the graphene problem , and suggests instructive physical analogies between two very different physical systems . a generalized coherent state approximation will be introduced that corresponds to a hartree  fock - bogoliubov ( hfb ) formalism subject to so(8 ) symmetry constraints . this permits quantitative evaluation of energy surfaces associated with spontaneously - broken symmetry . the so(8 ) highest symmetry will be shown to have an su(4 ) subgroup that recovers the known physics of su(4 ) quantum hall ferromagnetism as a special case , but implies in the more general case new low - energy physics that transcends su(4 ) quantum hall ferromagnetism . hence , a solvable and physically - illuminating approach to the rich low - energy structure of undoped graphene in strong magnetic fields will be proposed that reproduces known physics , but also suggests testable new physics in this complex system . comprehensive reviews of graphene physics may be found in refs . the presentation here will recall only a select set of features that will be relevant for subsequent discussion . undoped graphene is a 2-dimensional semiconductor with zero bandgap . it has a bipartite honeycomb lattice structure illustrated in fig . [ fig : a - b_sublattices ]        that corresponds to two interlocking triangular sublattices , labeled a and b. the two - fold degree of freedom specifying whether an electron is on the a sublattice or b sublattice is a spin - like quantum number termed the _ sublattice pseudospin_. the sublattice pseudospin behaves mathematically like the actual spin of the electron ( which will be introduced later ) , but it is a separate degree of freedom . the dispersion of energy with momentum for undoped graphene in the absence of a magnetic field is illustrated in fig . [ fig : dispersiongraphene ] . 0pt 0.44 electronic dispersion of graphene calculated in a tight - binding model with no magnetic field . two inequivalent points in the brillouin zone are labeled @xmath2 and @xmath3 . near these k - points the dispersion becomes linear , leading to _ dirac cones , _ as shown in the expanded view . for undoped graphene the fermi surface lies at the apex of the cones , where the level density vanishes and the effective electronic mass tends to zero . thus low - energy electrons are governed approximately by a dirac equation for massless electrons . the two inequivalent points @xmath2 and @xmath3 are not connected by reciprocal lattice vectors . the corresponding two - fold k degree of freedom may be viewed mathematically as a 2-valued spin - like quantity , commonly termed the _ valley isospin _ because of the valley - like structure in the dispersion of fig . [ fig : dispersiongraphene ] around the k - points . for brevity the valley isospin will sometimes be termed simply _ isospin . _ obviously this two - fold `` isospin '' degree of freedom is a pseudospin that should be distinguished physically from the actual spin of the electron , and from standard usage of isospin quantum numbers in nuclear and particle physics .    near the _ dirac cones _ ( inset to fig .  [ fig : dispersiongraphene ] ) the dispersion is linear , the density of electronic states tends to zero , and the electrons are described by a massless dirac equation in which the fermi velocity plays the role that the speed of light would play in an actual relativistic system . thus the low - energy electrons for undoped graphene in zero magnetic field behave to good approximation as _ massless chiral fermions , _ with the chirality representing the projection of the sublattice pseudospin on the direction of motion , not the projection of the actual spin . the vanishing of the density of states at the dirac point ( fermi surface ) implies that the transport properties of graphene are different from either a metal or a semiconductor . our interest will be primarily low - energy states in a strong magnetic field . for non - interacting electrons , the quantized levels may be found by solving the dirac equation for massless fermions with a vector potential appropriate to the applied magnetic field . the dispersion of energy with magnetic field strength for massless dirac electrons is illustrated in fig . [ fig : llgraphene_n0](a ) . consider the @xmath4 state , which corresponds to half filling of the fourfold - degenerate @xmath5 ll in graphene , as illustrated in fig . [ fig : llgraphene_n0](b ) . the graphene honeycomb lattice is bipartite , as illustrated in fig . [ fig : a - b_sublattices ] . the @xmath5 ll is located exactly at the dirac point corresponding to @xmath6 . for low - energy excitations in each valley labeled by @xmath2 or @xmath3 , the inter - valley tunneling may be ignored and the electrons in the valley reside entirely on either the a or b sublattice . hence , for the @xmath5 ll labeling with the valley isospin ( indicating whether the electron is in a @xmath2 or @xmath3 valley ) is equivalent to labeling with the sublattice pseudospin ( indicating whether the electron is on the a or b sublattice ) . this is reminiscent of a nel state with spins on two different sublattices , with a nel order defined by the difference in spins on the a and b sublattices . a quantum hall effect is signaled by a plateau in the hall conductance @xmath7 having quantized values @xmath8 where the _ filling factor _ @xmath9 is defined by @xmath10 with @xmath11 the charge - carrier density , @xmath12 the magnetic field strength , and @xmath13 the magnetic flux density in units of the fundamental flux quantum @xmath14 . these plateaus indicate formation of an _ incompressible quantum liquid . _ this is a compact way to say that the ground state is separated from excited states by an _ energy gap , _ which inhibits compression because of the energy required for excitation across the gap . the integral quantum hall effect ( iqhe ) in graphene is similar to the integral quantum hall effect for non - relativistic electrons in that it corresponds to the formation of incompressible states resulting from the complete filling of landau levels by weakly - interacting electrons . however , there are two important differences between the iqhe states observed in graphene and those observed in conventional 2d semiconductor heterostructures :    ( 1 )  in addition to the two - fold spin degeneracy ( in the absence of zeeman splitting ) , there is a two - fold valley degeneracy associated with the distinct @xmath2 and @xmath15 points in the first brillouin zone . thus the filling factor changes in steps of four between plateaus in the hall resistance for graphene . ( 2 )  for graphene the filling factor @xmath9 defined in eq . ( [ iqhegraph1.1 ] ) vanishes at the dirac point for particle  hole symmetric half filling of the graphene lattice , since the electron density @xmath11 tends to zero there . hence , in the absence of a zeeman effect or strong electronic correlations , there is no integral quantum hall effect in graphene for @xmath16 . ( a )  relativistic dispersion for massless dirac electrons as a function of magnetic field strength @xmath12 . states are labeled by a principle quantum number @xmath17 and a quantum number @xmath18 indicating particle - like ( @xmath19 ) and hole - like ( @xmath20 ) states . each landau level labeled by @xmath17 has a high orbital degeneracy @xmath21 ( see eq . ( [ degen1.1 ] ) ) , and an additional 4-fold degeneracy associated with spin and valley isospin . ( b )  one configuration for occupation of the @xmath5 landau level in monolayer graphene . the splitting and occupation are schematic only . the ground state would be a superposition of such configurations and in the su(4 ) symmetry limit ( which obtains for coulomb - only interactions ) the 4 levels shown labeled by valley @xmath2 or @xmath3 and spin up or down would be degenerate . in realistic cases the splitting often is small , suggesting approximate su(4 ) symmetry . ]    in graphene the analog of the integral quantum hall effect was first observed at filling factors @xmath22 by sweeping the field or the carrier density through the landau levels @xcite . this implies hall resistance quantization for filling factors in the sequence @xmath23 where @xmath11 is the charge carrier density and @xmath17 is the landau level index @xcite . this sequence is quite different from the integral quantum hall effect sequence observed in other 2d electron gases . however , it likely results from the same basic physics as the normal iqhe , modified by the 4-fold spin  valley degeneracies of non - interacting , massless dirac electrons .    the period @xmath24 in eq . ( [ iqhegraph1.2 ] ) is is a consequence of the approximate four - fold degeneracy of the graphene landau level , and the added @xmath25 ( which is not present in non - relativistic 2d systems ) is a berry phase effect that results from the special status of the @xmath5 state for massless dirac fermions : a quantum phase arises at the band degeneracy point associated with precession of the pseudospinor describing the 2-fold sublattice degree of freedom , which modifies the quantization condition for electronic orbits @xcite . the four - fold near - degeneracy of the landau level follows because the zeeman energy is small compared with the interaction energy , and the pseudospin degree of freedom representing the two inequivalent dirac cones at the corners of the brillouin zone ( k and k@xmath26 ) does not couple to external fields if the two sublattices are equivalent .      because of the dirac - cone dispersion with the fermi level located at the apex of the cones ( see fig . [ fig : dispersiongraphene ] ) , low - energy excitations in graphene occur in regions of reduced electron density , which disfavors electron correlations . but by placing a strong perpendicular magnetic field on the system the resulting landau quantization ( corresponding semiclassically to requiring that an integral number of debroglie wavelengths wrap around a cyclotron orbit , and implying that an integral number of magnetic flux quanta pass through the area bounded by the orbit ) leads to a bunching of levels into regions of locally high degeneracy separated by gaps . these regions of locally high level density may exhibit conditions more favorable for the development of strong correlations between electrons .    landau levels ( ll ) become strongly correlated when inter - ll excitations are of sufficiently high energy that they may be neglected and the low - energy excitations involve only intra - ll transitions . then the kinetic energy of the ll is a constant that can be omitted . this limit of strong electronic correlations has two important physical implications : ( 1 )  the approximate 4-fold spin  valley degeneracy of the graphene landau levels leads to quantum hall ferromagnetic states that will be discussed further below . ( 2 )  the strong correlations can produce incompressible states at partial filling of the ll that are reminiscent of the fractional quantum hall effect ( fqhe ) in semiconductor devices .    after the discovery of the integral quantum hall effect in graphene , experiments performed at higher magnetic field strengths observed fragile quantum hall states at filling factors @xmath27 . the @xmath28 states are thought to be the result of single - particle zeeman splitting of the landau levels but the states at filling factors @xmath29 are thought to be caused by electron  electron interactions breaking degeneracies of the @xmath5 landau level @xcite . quantum hall incompressible states having actual fractional filling were later observed at filling factors such as @xmath30 @xcite . the fractional states follow the standard composite fermion model @xcite sequence for filling factors @xmath31 , but only even - numerator fractions are seen for @xmath32 . these sequences , and the energy gaps for the corresponding incompressible states , are thought to reflect the interplay of strongly - correlated chiral electrons and the characteristic internal symmetries of graphene that will be discussed further below @xcite . as has been seen , in the quantum hall effect for 2d electron gases produced in semiconductor devices , two fundamental classes of incompressible states are found :    1 .   those where the requisite gaps are produced by complete filling of landau levels ( fermi energy lying between landau levels ) that are explained by weakly - interacting electrons subject to impurity - scattering localization . 2 .   those where the gaps are produced by strong electron correlations within a partially - filled landau level .    for the normal semiconductor quantum hall effect , the first case is commonly termed the integral quantum hall effect ( iqhe ) , because it leads uniquely to quantum hall plateaus at integral filling factors , and the second case is termed the fractional quantum hall effect ( fqhe ) , because it is characterized uniquely by quantum hall plateaus at fractional values of the filling factor .    in graphene , the anomalous counting implied by the four - fold spin  valley degeneracy for massless chiral fermions modifies this correspondence . one may again assume incompressible states divided into those that have gaps produced by weakly - interacting electrons filling landau levels completely and those that have gaps produced by correlations within a partially filled level . however ( 1 )  states of the first class occur at integral fillings but the integers of eq . ( [ iqhegraph1.2 ] ) are not those of the normal iqhe because of the anomalous counting implied by graphene s internal structure . ( 2 )  states of the second class may have either integral or fractional filling factors because of the anomalous counting . thus iqhe and fqhe may be used as shorthand labels in graphene , but it should be understood that what is meant by `` iqhe '' in graphene is the formation of incompressible states by weakly - interacting electrons completely filling landau levels , and by `` fqhe '' the formation of incompressible states by strongly - correlated electrons in a partially - filled level , irrespective of whether the observed filling factors are integers or rational fractions . our primary interest in this paper lies in those incompressible states formed by electron  electron and electron  lattice correlations in partially - filled landau levels , and thus in graphene `` fqhe '' states . in the normal two - dimensional electron gas produced in semiconductor devices the landau levels can contain @xmath33 states , where @xmath34 is the electronic charge and @xmath12 is the magnetic field . as has been seen , in graphene there is an additional 4-fold degeneracy associated with the spin and valley degrees of freedom . it is common to unite these four degrees of freedom in terms of an su(4 ) symmetry that is termed _ quantum hall ferromagnetism ( qhfm ) . _ su(4 ) symmetry for graphene in a strong magnetic field is expected when all four spin and valley levels are degenerate . two conditions must be satisfied for this condition to be fulfilled .    1 . landau - level mixing caused by inter - ll electronic transitions must be negligible . perturbations within a single ll that break the 4-fold spin  isospin symmetry must be small .    the resulting theory predicts quantum hall states that have no analog in semiconductor heterostructures . let us study these states by introducing an effective low - energy hamiltonian exhibiting approximate su(4 ) symmetry . the two - dimensional electronic spin degree of freedom and the two - dimensional electronic valley ( k ) degree of freedom are most elegantly expressed in terms of independent spin and valley isospin states . for the spin , introduce the pauli matrix vector @xmath35 , with the standard representation in terms of the @xmath36 matrices @xmath37 which obey the su(2 ) lie algebra @xmath38 } = i\epsilon_{ijk } \frac{\sigma_k}{2 } , \label{su2comm}\ ] ] where @xmath39 is the completely antisymmetric 3rd - rank tensor . the matrices ( [ algebra1.5 ] ) are assumed to operate on a spinor basis of spin - up and spin - down electrons @xmath40    a set of equations in the valley isospin space completely analogous to eqs . ( [ algebra1.5])([algebra1.6 ] ) for the spin space results if one defines the su(2 ) pauli - matrix representation for the valley isospin operators @xmath41 , @xmath42 which operate on the valley isospinor basis @xmath43 the operators @xmath44 and @xmath45 may now be used to define an effective low - energy hamiltonian having landau - level and internal spin and valley isospin degrees of freedom . the two largest energy scales for graphene in a strong magnetic field are the landau - level separation and the coulomb energy . at the charge neutral point ( fermi energy for undoped graphene ) the ll separation is approximately three times larger than the coulomb energy , which is in turn considerably larger than any additional terms in the interaction . therefore , a strategy is adopted here of ignoring excitations between landau levels and projecting onto the @xmath5 ll . at a quantitative level the landau level mixing can not be ignored ( see the discussion in ref . @xcite ) , but such an approximation gives the correct qualitative physics and the effect of excluded landau levels can be included to some degree by parameter renormalization , which will be sufficient for our purposes .    within this single landau level the hamiltonian is assumed to be dominated by a long - range coulomb interaction that is su(4 ) symmetric , with shorter - range interactions in spin and valley degrees of freedom ( originating in both screened electron  electron interactions and electron  phonon interactions ) causing su(4 ) symmetry breaking . to implement this a graphene hamiltonian projected onto the @xmath5 landau level is adopted that was proposed in ref . @xcite and employed further in ref . @xcite,@xmath46 where the valley - independent coulomb interaction @xmath47 may be expressed as @xmath48 @xmath49 is the short - range , valley - dependent interaction , @xmath50   \delta(\bm r_i - \bm r_j ) ,   \label{so5_1.1c}\ ] ] where the pauli matrices @xmath51 operate on the valley isospin and @xmath52 and @xmath53 are coupling constants , and the zeeman energy @xmath54 is given by @xmath55 where @xmath56 is the bohr magneton , @xmath12 is the magnetic field strength , the pauli matrices @xmath57 operate on the electronic spin degrees of freedom , and the @xmath58 direction for the spin space is chosen to be aligned with @xmath12 .      letting @xmath59 and @xmath60 , the set of 15 operators    @xmath61    [ algebra1.4 ]    is closed under commutation , defining an su(4 ) lie algebra that commutes with the coulomb interaction @xmath47 @xcite . thus , if @xmath49 and @xmath54 are small compared with @xmath47 in eq . ( [ so5_1.1a ] ) , the hamiltonian will have an approximate su(4 ) invariance . the operator @xmath62 represents the total spin and the operator @xmath63 represents the total valley pseudospin . in the @xmath5 landau level for graphene there is an equivalence between valley and sublattice degrees of freedom , so @xmath64 can be viewed as a nel vector in the @xmath5 landau level measuring the difference in spins on the a and b sublattices . the operators @xmath65 coupling spin and valley isospin will be discussed further below . the occurrence of su(4 ) symmetry and its explicit symmetry - breaking pattern depend on the values of the effective coupling parameters @xmath52 and @xmath53 . for the lattice spacings found in graphene , each can be estimated to be considerably smaller than the su(4)-symmetric coulomb term , so one may expect su(4 ) to be broken only weakly by explicit terms in realistic systems . four basic explicit symmetry - breaking patterns have been identified @xcite :    1 .  for arbitrary non - zero values of @xmath52 and @xmath53 , the symmetry is broken to @xmath66 where @xmath67 is associated with global conservation of spin and @xmath68 with conservation of its @xmath58 component , and @xmath69 is associated with conservation of the @xmath70 component of the valley isospin . ( conservation of @xmath70 implies physically that the difference in electronic densities between the @xmath2 and @xmath3 sites is invariant , which might be expected to be true for low - energy states having minimal scattering between valleys . ) in the absence of zeeman splitting spin is conserved , but only the @xmath58 component of the valley isospin is conserved . the full hamiltonian ( [ so5_1.1a ] ) including the zeeman term conserves only the @xmath58 components of the spin and valley isospin . if @xmath71 but @xmath72 , the symmetry is broken to @xmath73 in the absence of zeeman splitting this corresponds to conserving independent spin rotations for the wavefunction in each valley labeled by @xmath2 and @xmath3 , but breaking of the valley isospin to u(1)@xmath74 . the full hamiltonian ( [ so5_1.1a ] ) including the zeeman term conserves only the @xmath58 components of the spin in each valley separately , and the @xmath58 component of the valley isospin .    3 . if @xmath75 , the symmetry is broken to @xmath76 corresponding to full spin and valley isospin rotational symmetry in the absence of zeeman splitting . the complete hamiltonian ( [ so5_1.1a ] ) including the zeeman term conserves the su(2 ) isospin symmetry but only the @xmath58 component of the spin .    4 . if @xmath77 , the hamiltonian commutes with @xmath65 , @xmath62 , and @xmath70 , and these 10 operators generate the lie group so(5 ) , so @xcite @xmath78 where the @xmath79 symmetry is generated by @xmath80 . thus , in the absence of zeeman splitting the system exhibits an so(5 ) symmetry involving both spin and valley isospin . the full hamiltonian ( [ so5_1.1a ] ) including the zeeman term conserves the @xmath58 component of spin and the @xmath79 symmetry . the so(5 ) subgroup plays the role of a transitional symmetry connecting the nel - like states associated with @xmath81 and the states associated with valley degrees of freedom @xmath82 and @xmath83 . the subgroup structure for these four patterns of explicit su(4 ) symmetry breaking is illustrated in fig . [ fig : dynamicalchains_qhf].dynamicalchains_qhf figures / dynamicalchains_qhf_vert_compact_nobkg.eps 0pt 0.40 explicit symmetry - breaking structure for su(4 ) quantum hall ferromagnetism described by the hamiltonian ( [ so5_1.1a ] ) . these symmetries and explicitly - broken symmetries have proven extremely useful in understanding the states of graphene in a strong magnetic field @xcite . let us now consider using symmetries in an even more powerful way than that discussed in the preceding section . specifically , let us attempt to describe the quantum hall ferromagnetic behavior of graphene using _ fermion dynamical symmetries _ of an effective hamiltonian operating in a highly - truncated collective subspace . this has the potential to prescribe dynamics as well as taxonomy and conservation laws , within a many - body model having analytical solutions that illuminate the physics of quantum hall states in graphene .      for clarity of discussion a valence space will be assumed corresponding to a single landau level . to avoid cluttered notation the index @xmath17 labeling the landau level is suppressed and the fermion creation operator @xmath84 and the corresponding hermitian conjugate @xmath85 are introduced . the index @xmath86 takes one of two values ( @xmath87 ) labeling the valley isospin projections corresponding to valleys @xmath2 or @xmath3 , the electron spin polarization @xmath88 takes one of two values ( @xmath89 ) labeling spin - up or spin - down , and @xmath90 is a quantum number distinguishing degenerate states within a given landau level ( typically an angular momentum in symmetric gauge or a linear momentum in landau gauge ) . the fermionic operators @xmath91 and @xmath92 are assumed to obey @xmath93 where @xmath94 . by virtue of this anticommutator , the @xmath95 create a fermion in the state labeled by @xmath96 , the @xmath97 annihilate a fermion in the same state , and @xmath98 counts the number of fermions in the state labeled by @xmath96 . the four states representing possible combinations of @xmath86 and @xmath88 are displayed in table  [ tb : quantumnumbermapping ] ,    ccccc    ' '' ''    valley & @xmath86 & @xmath88 & @xmath99 & @xmath100     +    ' '' ''    @xmath2 & @xmath19 & @xmath101 & @xmath102 & @xmath103     + [ ]    @xmath2 & @xmath19 & @xmath104 & @xmath105 & @xmath106     + [ ]    @xmath3 & @xmath20 & @xmath101 & @xmath107 & @xmath108     + [ ]    @xmath3 & @xmath20 & @xmath104 & @xmath109 & @xmath110     +    and correspond physically to the four possible combinations of the electron being in either the @xmath2 or @xmath3 valley with either spin up or spin down ( see fig . [ fig : llgraphene_n0]b ) . for brevity the label @xmath111 displayed in the last column of the table will often be used to distinguish these states . the four basis states labeled by @xmath100 are illustrated graphically in fig . [ fig : graphenebasis ] . figures / graphenebasis.eps 0pt 0.16 the four isospin  spin basis vectors @xmath112 of table  [ tb : quantumnumbermapping ] . table  [ tb : quantumnumbermapping ] also displays a unique mapping of these four states to a label @xmath99 that takes the four possible projection quantum numbers @xmath113 of a fictitious angular momentum @xmath114 ; the motivation for this mapping will become apparent later . let us now introduce an operator @xmath115 that creates a pair of electrons , one in the @xmath116 level and one in the @xmath117 level , with the total @xmath90 of each pair coupled to zero term by term , @xmath118 and its hermitian conjugate @xmath119 , which annihilates a corresponding electron pair . each index @xmath100 or @xmath120 ranges over four values , implying 16 components in eq . ( [ algebra1.1 ] ) . however , the pair wavefunction must satisfy the pauli principle which , upon expanding the indices @xmath100 and @xmath120 using table  [ tb : quantumnumbermapping ] , eliminates the four diagonal @xmath121 possibilities . furthermore , because of the antisymmetry requirement the pair creation operators are constrained by @xmath122 implying that only half of the remaining 12 operators @xmath123 are independent . thus eq . ( [ algebra1.1 ] ) defines six independent operators @xmath123 , with six independent hermitian conjugates @xmath124 . let us introduce in addition to these pairing operators the 16 particle  hole operators @xmath125 through @xmath126 where @xmath127 is the kronecker delta and @xmath128 is the total degeneracy the single landau level ( see eq . ( [ degen1.2 ] ) ) . the commutators for the 28 operators @xmath124 , @xmath123 , and @xmath12 are found to be @xcite    @xmath129 } & = -b_{db } \delta_{ac } - b_{ca}\delta_{bd } + b_{cb } \delta_{ad } + b_{da } \delta_{bc } \label{algebra1.3a } \\ { \left[b_{ab } , b_{cd}\right ] } & = \delta_{bc } b_{ad } - \delta_{ad}b_{cb } \label{algebra1.3b } \\ { \left[b_{ab } , a^\dagger_{cd}\right ] } & = \delta_{bc } a^\dagger_{ad } + \delta_{bd } a^\dagger_{ca } \label{algebra1.3c } \\ { \left[b_{ab } , a_{cd}\right ] } & = -\delta_{ac } a_{bd } - \delta_{ad } a _ { cb } , \label{algebra1.3d}\end{aligned}\ ] ]    [ algebra1.3 ]    which is isomorphic to an so(8 ) lie algebra . thus , the 28 members of the operator set @xmath130 exhibit so(8 ) symmetry under commutation . the ultimate justification for introducing this so(8 ) lie algebra will lie in the results to which it will lead . however , a general discussion of why the known physics of graphene suggests the efficacy of such an algebra in understanding the possible collective modes of the system may be found in appendix [ h : appendixa ] .      the su(4 ) generators ( [ algebra1.4 ] ) may be expressed in terms of the operators ( [ algebra1.2 ] ) by employing eqs . ( [ algebra1.5])([algebra1.7e ] ) in eq . ( [ algebra1.4 ] ) and utilizing the equivalences in table  [ tb : quantumnumbermapping ] to expand the indices for @xmath125 from eq . ( [ algebra1.2 ] ) . for example , consider the spin operator @xmath131 . from eq . ( [ algebra1.4a ] ) one may write @xmath132 & \left . \qquad\qquad\qquad -i c^\dagger_{-\uparrow m_k } c_{-\downarrow   m_k}{^{\vphantom{\dagger}}}+ i c^\dagger_{-\downarrow m_k } c_{-\uparrow m_k}{^{\vphantom{\dagger}}}\right ) \\ & = -ib_{12 } + i b_{21 } -i b_{34 } + i b_{43}.\end{aligned}\ ] ] by such methods one finds that the spin operators ( [ algebra1.4a ] ) may be expressed in terms of the @xmath125 generators of so(8 ) as    @xmath133    [ algebra1.8 ]    the valley isospin operators ( [ algebra1.4b ] ) as    @xmath134    [ algebra1.9 ]    the nel vector of eq . ( [ algebra1.4c ] ) as    @xmath135    [ algebra1.10 ]    and the operators @xmath65 of eq . ( [ algebra1.4d ] ) as    @xmath136    [ algebra1.11 ]    the inverse transformations expressing the @xmath125 in terms of the @xmath137 are given in appendix [ h : appendixb ] . hence the su(4 ) algebra generated by the operators in eq . ( [ algebra1.4 ] ) is a subalgebra of the so(8 ) algebra , with its generators corresponding to particular linear combinations of the subset of so(8 ) generators defined by the particle  hole operators @xmath125 in eq . ( [ algebra1.2 ] ) . in fact , the present formalism also could be constructed by starting with the su(4 ) algebra generated by the @xmath125 operators and expanding that to the so(8 ) algebra of eq . ( [ algebra1.3 ] ) by adding pairing operators , using the motivation discussed in appendix [ h : appendixa ] . one wishes to investigate possible collective modes for graphene electrons undergoing strong correlations within a single landau level . pairs of fermions often afford a convenient basis for discussing collective states , so let us consider some possible configurations involving pairs of electrons in graphene . let us first consider the degeneracies of undoped graphene placed in a strong magnetic field , confining the discussion to the case of a single landau level for simplicity . the single - particle states within the landau level will be assumed labeled by the quantum numbers @xmath138 , where @xmath17 is the principle quantum number labeling the landau level and @xmath139 is a quantum number distinguishing the degenerate states within the landau level . in the absence of spin and valley degrees of freedom , the states @xmath138 of the landau level are assumed to hold a maximum of @xmath140 electrons ( for consistency with earlier applications of the so(8 ) algebra in nuclear physics , @xmath141 will be defined to be the electron pair degeneracy , so that @xmath140 is the electron degeneracy ) . from the solution of the dirac equation in a magnetic field @xmath142 where @xmath12 is the strength of the magnetic field , @xmath143 is the area of the two - dimensional sample and @xmath144 defines the magnetic flux quantum . but graphene has in addition @xmath145 internal degrees of freedom associated with the @xmath146 space . thus there are four copies of each landau level in graphene and the total electron degeneracy @xmath147 is given by @xmath148 some pair degeneracies calculated from eq . ( [ degen1.2 ] ) as a function of domain size and magnetic field strength are displayed in table  [ tb : degeneracytable ] ,    cccc    ' '' ''    domain size & @xmath149 & @xmath150 & @xmath151     +    ' '' '' @xmath152 & @xmath153 & @xmath154 & @xmath155     + [ ]    @xmath156 & @xmath157 & @xmath158 & @xmath159     +    where it is assumed that the collective wavefunction is delocalized over the entire domain size .    the _ fractional occupation _ @xmath160 of the single landau level [ not to be confused with the filling factor @xmath9 given in eqs . ( [ iqhegraph1.1 ] ) and ( [ fillingfactor ] ) ] may be defined as @xmath161 where @xmath17 is the number of electrons and @xmath162 is the number of electron pairs . for half - filling of the @xmath5 landau level located at the fermi surface ( corresponding to the ground state of undoped graphene ) the electron number @xmath163 is then @xmath164 these degeneracies and occupation numbers are just the standard results for relativistic landau levels in a 2d electron gas subject to a strong perpendicular magnetic field , but modified by the graphene spin and valley degeneracies . graphene exhibits both integral and fractional quantum hall effects but the filling factors are anomalous relative to those for standard quantum hall effects in semiconductor heterostructures . this is because of    1 .   the 4-fold degeneracy associated with the spin and valley degrees of freedom , which introduces factors of four in the counting . 2 .   the nature of the dirac solution , illustrated in fig . [ fig : llgraphene_n0 ] , for which the negative - energy solutions may be interpreted as electron holes , the positive - energy solutions as electrons , and the @xmath5 level is unique , being half - filled in the neutral ground state ( equivalently , it may be thought of as being shared equally by particles and holes ) .    because of the particle - hole symmetry , the charge carriers change sign near the dirac points and the hall conductivity vanishes at charge neutrality ( the electron number density tends to zero at a dirac point ) . for this reason , the filling factor for graphene must be defined _ relative to the charge - neutral state . _ at charge neutrality the @xmath5 landau level is half filled and when the @xmath5 ll is completely full the filling factor is @xmath165 , from eq . ( [ iqhegraph1.2 ] ) . the quantum hall filling factor @xmath9 may be related to the landau level fractional occupation @xmath160 employed in the present formalism by @xmath166 therefore , half filling of the @xmath5 landau level corresponds to a fractional occupation @xmath167 but to a filling factor @xmath4 , @xmath168 corresponds to @xmath169 ( completely empty ) , @xmath170 corresponds to @xmath171 filling , @xmath172 corresponds to @xmath173 filling , and @xmath174 corresponds to @xmath175 ( completely full ) . consider the states created by repeated application of the pair creation operator @xmath115 defined by eq . ( [ algebra1.1 ] ) to the pair vacuum . it is useful to classify states according to a seniority - like quantum number @xmath176 defined to be the number of particles in the system not coupled to one of the pairs defined in eq . ( [ algebra1.1 ] ) . the @xmath177 subspace will be of particular interest since it will contain states of maximal collectivity with respect to the pairs ( [ algebra1.1 ] ) . an @xmath178-pair state in the @xmath177 subspace is given by @xmath179 where the total pair number @xmath178 is @xmath180 with @xmath17 giving the total number of electrons and @xmath181 giving the number of electron pairs created by @xmath182 operating on the vacuum state . for our discussion here it will always be assumed that one is dealing with @xmath177 states , corresponding physically to no broken pairs .      from eq . ( [ algebra1.3 ] ) , the 16 operators @xmath125 are closed under commutation and form a @xmath184 subalgebra of the so(8 ) algebra ( [ algebra1.3 ] ) . let us now investigate the irreducible representations ( irreps ) that are associated with the @xmath185 subgroup chain in the @xmath177 representations .      for @xmath177 at half filling , the number of pairs is @xmath186 and the highest - weight u(4 ) representation is given by @xmath187 . let us define a highest - weight ( hw ) state in the @xmath177 space , and choose it to correspond to the pair state with maximal value of @xmath99 from table  [ tb : quantumnumbermapping ] , which results from placing one electron in the @xmath188 state and one electron in the @xmath189 state . thus , for @xmath190 pairs the highest weight state is given by @xmath191 where the sum runs over the @xmath192 states in the landau level labeled by the @xmath193 quantum number . the other states of the irreducible representation may then be created by the cartan  dynkin algorithm , which consists of using raising and lowering operators in the weight space to construct successively all the other states beginning with the highest weight state @xcite . the state in eq . ( [ hwstate ] ) appears to have a quite complex form , involving a sum with number of terms equal to the pair degeneracy of the landau level raised to a power equal to the pair degeneracy ( with the pair degeneracy typically a large number ) . however , the actual structure of this state is considerably simpler than eq . ( [ hwstate ] ) would suggest because of the pauli principle . as an illustrative example of this assertion , let s construct explicitly the highest - weight state for the case @xmath194 , corresponding to @xmath195 pairs in a single landau level . writing the sum over @xmath196 in eq . ( [ hwstate ] ) out term by term gives @xmath197 where in raising the sum of operators inside the parentheses to the @xmath198 power , all products containing two or more creation operators with the same index vanish because of the pauli principle . similar considerations apply for arbitrary values of @xmath199 and in general the highest - weight state is given by @xmath200 where the simplification in going from the second to the third line is a consequence of the antisymmetry of the fermion creation operators ( the pauli principle ) implied by eq . ( [ anticommutator ] ) . thus the highest - weight state is a product state of pairs , one pair for each of the @xmath190 levels labeled by @xmath90 in the landau level .      by the cartan  dynkin algorithm , other states in the @xmath177 subspace can be constructed by applying successively to the highest - weight state appropriate lowering and raising operators . these will be functions of the generators @xmath125 , so for an arbitrary state @xmath201 in the weight space one has schematically @xmath202 where the function @xmath203 is specified by the cartan  dynkin procedure . as an example , consider the action of the valley isospin lowering operator @xmath204 on the highest - weight state . from eqs . ( [ algebra1.9 ] ) and ( [ algebra1.2 ] ) , @xmath205 thus the state @xmath201 created by applying @xmath204 to @xmath206 is @xmath207 c_{{\scriptscriptstyle k\uparrow } m_k}^\dagger c_{{\scriptscriptstyle k\downarrow } m_k}^\dagger \ , { |0\rangle } \nonumber \\ & = \prod_{m_k }   \left ( c_{3 m_k}^\dagger c_{2 m_k}^\dagger + c_{4 m_k}^\dagger c_{1 m_k}^\dagger \right ) { |0\rangle }   , \label{u0pairs1.7}\end{aligned}\ ] ] where the simplifications are because the only terms that survive correspond to those where an annihilation operator in a factor inside the square brackets is exactly balanced by a creation operator from the factor outside the square brackets . likewise , the other states of the @xmath177 representation can be constructed by using successive applications of raising and lowering operators fashioned from the generators defined in in eqs . ( [ algebra1.8])([algebra1.11 ] ) .      from the preceding discussion , for @xmath208 the states may be written as @xmath209 { |0\rangle } , \label{upairs1.8}\end{aligned}\ ] ] where @xmath210 denote valley isospin projection quantum numbers and @xmath211 denote spin projection quantum numbers . this is the same form as the most general collective pair state used by kharitonov @xcite in his classification of possible broken symmetry states for the @xmath5 landau level in graphene ( see eq . ( [ basis1.0 ] ) in the appendix ) . thus , for undoped graphene the general pairing wavefunction ( [ u0pairs1.0 ] ) characteristic of the so(8 ) @xmath183 su(4 ) dynamical symmetry is in fact _ equivalent to the product form ( [ upairs1.8 ] ) _ employed in standard discussion of quantum hall ferromagnetism , for which the summations are over the internal @xmath212 rather than landau @xmath213 degrees of freedom . the equivalence of eqs . ( [ upairs1.8 ] ) and ( [ u0pairs1.0 ] ) , despite their superficially very different forms , is a fundamental consequence of the pauli principle acting in the collective fermionic pair subspace , which greatly restricts allowed pair configurations . the equivalence established above implies that the present so(8 ) formalism can be used to derive the framework used in ref . @xcite to describe possible collective states in graphene , which establishes an essential connection between the standard discussion of su(4 ) quantum hall magnetism in graphene and the present more general so(8 ) dynamical symmetry formalism . as will be demonstrated further below , this permits the so(8 ) formalism to ( 1 )  encompass the established physics of su(4 ) quantum hall ferromagnetism , ( 2 )  extend quantum hall ferromagnetism to suggest possible additional collective graphene physics beyond the su(4 ) limit , and ( 3 )  provide _ analytical solutions _ for the states suggested by spontaneous breaking of su(4 ) symmetry that can not be obtained using su(4 ) symmetry and must be addressed numerically , or with effective field theory approximations , within standard su(4 ) quantum hall ferromagnetism . furthermore , it will be seen below that the equivalence of product and paired forms for the wavefunction implies a deep formal connection between the collective states resulting from strong electron correlations in graphene landau levels and the collective states produced by strong nucleon correlations in nuclear structure physics , and a suggestive formal connection to the properties of the strongly - correlated electrons responsible for high - temperature superconductivity . the preceding discussion has established that the fermion dynamical symmetry method applied to undoped graphene in a strong magnetic field has one dynamical symmetry chain @xmath185 that recovers exactly su(4)-symmetric quantum hall ferromagnetism . since the @xmath125 operators introduced in eq . ( [ algebra1.2 ] ) form an su(4 ) subgroup of so(8 ) that is in one - to - one correspondence with the operators used to formulate the effective low - energy hamiltonian ( [ so5_1.1a ] ) ,  [ h : fdsgraphene ] implies that all of the physics associated with this effective hamiltonian that has been discussed in the prior literature ( see @xcite and references cited therein ) is implicit in the present formalism . furthermore , the discussion of  [ h : pairrepresentations ] shows that the pair basis ( [ u0pairs1.0 ] ) of the truncated collective subspace for the so(8 ) fermion dynamical symmetry is in fact identical to the most general wavefunction ( [ upairs1.8 ] ) that has been proposed @xcite for collective states breaking the su(4 ) symmetry , despite its superficially very different form . thus the hilbert - space truncation implied by the collectively - paired so(8 ) subspace ( [ u0pairs1.0 ] ) recovers the understanding in the existing literature of the classes of states to be expected from spontaneous breaking of the su(4 ) symmetry by valley - dependent correlations .    however , the existing discussions of these collective states in terms of broken su(4 ) symmetry have been largely qualitative , and have turned to numerical simulations to discuss the actual structure and energy of the states . it will now be demonstrated that the present formalism is capable not only of classifying , but also of addressing the quantitative nature of those collective states _ in analytical fashion . _ furthermore , it will be shown that the so(8 ) highest symmetry implies subgroup chains in addition to @xmath214 that are associated with spontaneous breaking of the symmetry by correlations and have not been discussed in the previous literature and that may play a role in graphene . let us begin that discussion by first transforming to a more convenient representation of the so(8 ) generators . this new representation will be physically equivalent to the original representation , but will offer some advantages in interpretation , and will expose an unexpected relationship between graphene physics and that of a very different field , nuclear structure physics . for the pair creation operators defined in eq . ( [ algebra1.1 ] ) , each electron creation operator @xmath91 carries both spin and valley isospin ; hence the products @xmath215 correspond to a clebsch  gordan series representing sums of terms having different values of total spin and total isospin . likewise , in the particle  hole operators of eq . ( [ algebra1.2 ] ) each creation operator @xmath91 and annihilation operator @xmath92 carries spin and isospin , so the product @xmath216 in eq . ( [ algebra1.2 ] ) represents a superposition of states carrying different total spin and total valley isospin . these representations with indefinite spin and isospin will be termed _ uncoupled representations . _    on physical grounds , the spin is expected to be conserved ( if the zeeman term in the hamiltonian is neglected ) and the valley isospin is expected to be approximately conserved for low - energy excitations . thus , it is desirable to use the uncoupled representation of the pairing and particle  hole operators to construct new _ coupled representations _ that have good total spin and good total valley isospin quantum numbers for bilinear operators .      using standard angular momentum coupling theory @xcite , an electron pair creation operator coupled to good spin and valley isospin may be defined by @xmath217 where @xmath143 is the total spin of the pair with @xmath218 its projection , @xmath219 is the total valley isospin of the pair with @xmath220 its projection , and @xmath221 are clebsch  gordan coefficients for the angular momentum sum @xmath222 that couple the pair to good total spin or total valley isospin . antisymmetry implies that the the pair wavefunction can have only @xmath223 , @xmath224 ; or @xmath225 , @xmath226 ( spin - triplet , isospin - singlet ; or spin - singlet , isospin - triplet pairs ) . explicitly the possibilities are @xmath227 { a_{\scriptscriptstyle 00}^{\raisebox{0.2em}{$\scriptstyle\dagger$ } { \scriptscriptstyle   0 1 } } } = a_{14}^\dagger + a_{23}^\dagger \quad { a_{\scriptscriptstyle 01}^{\raisebox{0.2em}{$\scriptstyle\dagger$ } { \scriptscriptstyle   0 1 } } } = \sqrt 2 a_{13}^\dagger \qquad { a_{\scriptscriptstyle 0 - 1}^{\raisebox{0.2em}{$\scriptstyle\dagger$ } { \scriptscriptstyle   0 1 } } } = \sqrt2 a_{24}^\dagger , \end{array } \label{coupled1.3}\ ] ] with the hermitian conjugates of eq .  ( [ coupled1.3 ] ) giving the six corresponding pair annihilation operators in coupled representation . it is useful to define an alternative set of six coupled pairing operators @xmath228 and @xmath229 according to @xmath230 d_0^\dagger = \tfrac{1}{\sqrt2 } \,{a_{\scriptscriptstyle 00}^{\raisebox{0.2em}{$\scriptstyle\dagger$ } { \scriptscriptstyle   0 1 } } } = \tfrac{1}{\sqrt2 } \left ( a_{14}^\dagger + a_{23}^\dagger \right ) \\[7pt ] d_1^\dagger = \tfrac{1}{\sqrt2 } \ , { a_{\scriptscriptstyle 01}^{\raisebox{0.2em}{$\scriptstyle\dagger$ } { \scriptscriptstyle   0 1 } } } = a_{13}^\dagger \quad d_{-1}^\dagger = \tfrac{1}{\sqrt2 } \,{a_{\scriptscriptstyle 0 - 1}^{\raisebox{0.2em}{$\scriptstyle\dagger$ } { \scriptscriptstyle   0 1 } } } = a_{24}^\dagger \\[7pt ] d_2^\dagger = \tfrac{1}{\sqrt2 } \ , { a_{\scriptscriptstyle 10}^{\raisebox{0.2em}{$\scriptstyle\dagger$ } { \scriptscriptstyle   1 0 } } }   = a^\dagger_{12 } \quad d_{-2}^\dagger = \tfrac{1}{\sqrt2 } \,{a_{\scriptscriptstyle -10}^{\raisebox{0.2em}{$\scriptstyle\dagger$ } { \scriptscriptstyle   1 0 } } } = a_{34}^\dagger   \end{array }   \label{coupled1.4}\ ] ] and the six corresponding hermitian conjugates @xmath143 and @xmath231 . the physical meaning of these pairs may be deduced by constructing the corresponding electronic configurations . consider @xmath232 . from eqs . ( [ coupled1.4 ] ) and ( [ algebra1.1 ] ) @xmath233 where in the last step the correspondence between the index @xmath234 and the valley ( @xmath2 or @xmath3 ) and spin ( @xmath235 ) labels in table  [ tb : quantumnumbermapping ] has been invoked . this implies that @xmath236 creates a state with one spin - up and one spin - down electron on each equivalent site @xmath2 in the brillouin zone , as illustrated schematically in fig . [ fig : actionddagger2 ] . figures / actionddagger2.eps 0pt 0.35 the action of the pair creation operator @xmath237 on the vacuum state is to create a charge density wave with a spin - singlet pair on each site @xmath2 and no electrons on the @xmath3 sites . this is a component of a lattice - scale charge density wave , since the charge differs by two electronic units between adjacent sites . likewise , one finds that @xmath238 creates a charge density wave as in fig . [ fig : actionddagger2 ] , but with the spin - singlet pairs on the @xmath15 sites . the pair configurations produced by all generators of eq . ( [ coupled1.4 ] ) operating on the pair vacuum @xmath239 are summarized in fig . [ fig : s_d_pairs_brillouin ] . figures / s_d_pairs_brillouin_withq.eps 0pt 0.30 configurations created by the @xmath228 and @xmath240 pair creation operators of eq . ( [ coupled1.4 ] ) , and for the linear combinations @xmath241 of eq . ( [ qdef ] ) , operating on the pair vacuum @xmath239 . the upper left box illustrates the valley isospin labeling . a total of six valence electrons as assumed to be distributed on sites in the brillouin zone . in each pair configuration , location of the dots ( @xmath2 or @xmath3 site ) indicates the valley isospin , spin - up electrons are indicated by up arrows , and spin - down electrons are indicated by down arrows . also shown are the configurations generated by the linear combinations @xmath242 which will be useful in later discussion . kharitonov has given a general classification of low - lying collective modes for the @xmath5 landau level of graphene in terms of collective pairs @xcite . the collective pairs created by the so(8 ) pair generators in fig .  [ fig : s_d_pairs_brillouin ] are similar physically to the pairs identified by kharitonov , as will now be described . ( 1 )  the configuration generated by @xmath243 is to the difference of two terms , each with alternating spin - up and spin - down on adjacent sites , implying that all spins on the a sublattice ( identified with valley @xmath2 ) point in one direction and all spins on the b sublattice ( identified with valley @xmath15 ) point in the opposite direction . each term corresponds to a spin density wave ( af order ) , with a nel vector defined by the difference in total spins on the two sublattices serving as an order parameter , but because of the difference of the two terms the net af order for this configuration is zero ( see  [ sh : orderparameters ] and table  [ tb : orderparms ] ) . ( 2 )  the configurations generated by @xmath244 are spin - singlet charge density waves , with alternating charges of two and zero units on adjacent sites . an appropriate order parameter is the difference in charge between the a and b sublattices . ( 3 )  the configurations generated by @xmath245 have one spin on each site , all pointing in the same direction ; this is a ferromagnetic state , with the net spin as an order parameter . ( 4 )  the configuration generated by @xmath246 is the same as that generated by @xmath243 , except for a positive sign for the second term . this also implies alternating spins on adjacent sites and af order for each term , but the total af order vanishes because of the contribution of the two terms . ( 5 )  the configurations corresponding to @xmath247 are states with af order characterized by the difference in spins on the two sublattices labeled by @xmath2 and @xmath3 .    thus the coupled - representation pairs carrying good spin and valley isospin quantum numbers in eq . ( [ coupled1.4 ] ) represent physical degrees of freedom already discussed in the literature as candidate collective modes representing spontaneous breaking of the su(4 ) graphene symmetry by interactions in a single partially - filled landau level . it is desirable to express the particle  hole generators of eq . ( [ algebra1.2 ] ) in coupled representation . let us begin by introducing a set of operators @xmath248 with the definition @xmath249 where @xmath250 and @xmath251 take the values of the fictitious angular momentum projection @xmath99 in table  [ tb : quantumnumbermapping ] , providing a labeling equivalent to that of @xmath100 and @xmath120 in @xmath125 , with @xmath250 or @xmath251 values @xmath252 mapping to @xmath100 or @xmath120 values @xmath253 , respectively . for example , from table  [ tb : quantumnumbermapping ] , @xmath254 and @xmath255 label the same quantity , which is defined in eq . ( [ algebra1.2 ] ) . from the selection rules for angular momentum coupling in eq . ( [ multipole1.1 ] ) , the index @xmath256 can take the values @xmath257 , with @xmath258 projections @xmath259 for each possibility , which gives a total of 16 operators @xmath260 . by inserting the explicit values of the clebsch  gordan coefficients the 16 independent @xmath261 may be evaluated in terms of the @xmath125 . @xmath262 p^1_0 & =   \sqrt{\tfrac{9}{20}}(b_{11 } - b_{44 } ) + \sqrt{\tfrac{1}{20}}(b_{22 } -b_{33 } )   \nonumber   \\   & =   \sqrt{\tfrac{9}{20}}(n_1 - n_4 ) + \sqrt{\tfrac{1}{20}}(n_2 -n_3 ) \nonumber   \\[2pt ] p^1_1 & = -\sqrt{\tfrac{3}{10 } } \,b_{12 } - \sqrt{\tfrac{4}{10 } } \,b_{23 } -\sqrt{\tfrac{3}{10}}\ , b_{34 } \nonumber \\[2pt ]   p^1_{-1 } & =   \sqrt{\tfrac{3}{10 } } \,b_{21 } + \sqrt{\tfrac{4}{10 } } \,b_{32 } + \sqrt{\tfrac{3}{10}}\ , b_{43 } \nonumber   \\[2pt ]   p^2_0 & =   \tfrac{1}{2 } ( b_{11 } - b_{22 } + b_{44 } - b_{33 } )   \nonumber   \\   & = \tfrac12(n_1 - n_2 + n_4 - n_3 ) \nonumber   \\[2pt ]     p^2_1 & =   \tfrac{1}{\sqrt{2 } } ( b_{34 } - b_{12 } )   \qquad   p^2_{-1 } = \tfrac{1}{\sqrt2 } ( b_{21 } - b_{43 } )   \label{pofb }   \\[2pt ]   p^2_2 & =   -\tfrac{1}{\sqrt{2 } } ( b_{13 } + b_{24 } )   \qquad   p^2_{-2 } =   -\tfrac{1}{\sqrt2 } ( b_{31 } + b_{42 } ) \nonumber   \\[2pt ]    p^3_0 & =   \sqrt{\tfrac{1}{20 } } ( b_{11 } - b_{44 } )    + \sqrt{\tfrac{9}{20 } } ( b_{33}-b_{22 } )   \nonumber   \\   & =   \sqrt{\tfrac{1}{20 } } ( n_1 - n_4 )    - \sqrt{\tfrac{9}{20 } } ( n_2 - n_3 ) \nonumber   \\[2pt ]    p^3_1 & = -\sqrt{\tfrac15 } b_{12 } + \sqrt{\tfrac35 } b_{23 } - \sqrt{\tfrac15 } b_{34 } \nonumber   \\[2pt ]   p^3_{-1 } & =   \sqrt{\tfrac15 } b_{21 } - \sqrt{\tfrac35 } b_{32 }   + \sqrt{\tfrac15 } b_{43 } \nonumber   \\[2pt ]    p^3_2 & =   \sqrt{\tfrac12 } ( b_{24 } -   b_{13 } )   \qquad    p^3_{-2}=   \sqrt{\tfrac12 } ( b_{42 } - b_{31 } ) \nonumber   \\[2pt ]   p^3_3 & =   -b_{14 }   \qquad   p^3_{-3 } =   b_{41 } , \nonumber\end{aligned}\ ] ] where the @xmath263 are number operators for each of the four states and the total particle number @xmath17 is the sum over the four states labeled by @xmath100 in table  [ tb : quantumnumbermapping ] , @xmath264 it will be convenient notationally to sometimes replace the operator @xmath265 with the operator @xmath266 , according to the relationship @xmath267 where @xmath147 is the degeneracy of the space for the particles that participate in the so(8 ) symmetry . physically @xmath268 is one half the particle number measured from half filling ( which corresponds to @xmath269 ) .      because the six operators defined by eq . ( [ coupled1.4 ] ) , their six hermitian conjugates , and the 16 operators defined by eq . ( [ multipole1.1 ] ) are independent linear combinations of the so(8 ) generators defined in eqs . ( [ algebra1.1 ] ) and ( [ algebra1.2 ] ) , the 28 operators @xmath270 also close an so(8 ) algebra under commutation . the so(8 ) commutation relations for the coupled representation @xmath271 are given explicitly by @xcite    @xmath272 } & =   -2s_0    \label{commcoupled1.1 }   \\   { \left[d_{\mu ' } , d^\dagger_{\mu}\right ] } & =   -2\delta_{\mu\mu ' } s_0 + \sum_{t { \rm\ odd } } ( -1)^{\mu ' }   \nonumber   \\   & \qquad \times \   { c_{2 , \scriptscriptstyle-\mu ' 2 \mu}^{\,t,\scriptscriptstyle\mu-\mu ' } }   { \left\ {   \begin{matrix } 2&2&t\\#4&\tfrac32&\tfrac32 \end{matrix}\right\ } }   p_{\scriptscriptstyle\mu,-\mu'}^t \label{commcoupled1.2 } \\ { \left[d_\mu^\dagger , s\right ] } & = p_\mu^2 \label{commcoupled1.3 } \\ { \left[p^r_\mu , s^\dagger\right ] } & = 2 \delta_{r2}d^\dagger_\mu + 2 \delta_{r0 } \delta_{\mu0 } s^\dagger \label{commcoupled1.4 } \\ { \left[p_{\mu'}^r , d_\mu^\dagger\right ] } & = 2(-1)^{\mu ' } \delta_{r2 } \delta_{-\mu\mu ' } -4\sqrt{5(2r+1 ) } \nonumber \\ & \qquad\times\   { c_{r \scriptscriptstyle\mu ' 2 \mu}^{2,\scriptscriptstyle\mu+\mu ' } } { \left\ {   \begin{matrix } 2&2&r\\#4&\tfrac32&\tfrac32 \end{matrix}\right\ } } d_{\scriptscriptstyle\mu+\mu'}^\dagger \label{commcoupled1.5 } \\ { \left[p^r_{\mu ' } , p^s_\mu\right ] } & = 2 ( -1)^{r+s } \sqrt{(2r+1)(2s+1 ) }   \,\,\sum_t   { c_{r { \scriptscriptstyle\mu ' } s { \scriptscriptstyle\mu}}^{\,t , { \scriptscriptstyle \mu+\mu ' } } } \nonumber \\ & \times   \left[1-(-1)^{r+s+t}\right ] { \left\ {   \begin{matrix } r&s&t\\#4&\tfrac32&\tfrac32 \end{matrix}\right\ } } p^{\,t}_{\scriptscriptstyle\mu+\mu ' }   \label{commcoupled1.6}\end{aligned}\ ] ]    [ commcoupled ]    where @xmath266 is defined in eq . ( [ p0n ] ) and @xmath273 denotes the wigner 6-@xmath274 symbol @xcite for the recoupling of three angular momenta to good total angular momentum . the action of the so(8 ) pair creation operators on the pair vacuum @xmath178 times creates a @xmath275-particle state @xcite , @xmath276 where the total number of pairs is @xmath277 . the portion of the full hilbert space that is spanned by the states ( [ collective1.1 ] ) will be termed the _ collective subspace . _ it will play an important role in subsequent discussion where it will be shown that the so(8 ) symmetry may be used to construct effective hamiltonians that are diagonal in this space , and that the generators of so(8 ) do not couple the subspace to the remainder of the space . the reason for our alternative labeling of the states in table  [ tb : quantumnumbermapping ] in terms of the index @xmath99 , and our particular choices of phases and normalizations in equations , can now be made clear . with these labelings and choices the six coupled particle  particle operators @xmath228 and @xmath278 defined in eq . ( [ coupled1.4 ] ) , their six hermitian conjugates @xmath143 and @xmath231 , and the 16 coupled particle  hole operators @xmath279 defined in eq . ( [ multipole1.1 ] ) , are mathematically in one - to - one correspondence with the 28 generators for the ginocchio so(8 ) model @xcite and the so(8 ) fermion dynamical symmetry model @xcite . these have found broad application in nuclear structure physics @xcite and may be viewed as a microscopic justification for the interacting boson model ( ibm ) @xcite , which is one of the most commonly used phenomenological models in nuclear structure physics . this correspondence has three important implications : ( 1 )  mathematically , the group - theoretical methodology obtained for so(8 ) already in nuclear physics applications may be appropriated for use in the graphene problem . ( 2 )  physically the nature of the generators for the nuclear physics and graphene so(8 ) symmetries are fundamentally different , but analogs of physical interpretations applied already for nuclear physics so(8 ) symmetries may shed light on the graphene problem . ( 3 )  philosophically , the so(8 ) correspondence between graphene and nuclear structure physics implies a satisfying convergence of mathematical reasoning and physical abstraction in two completely different scientific subfields . this convergence will be elaborated further in  [ h : analogy_su4 ] . the so(8 ) group has various subgroups ( subsets of generators closed under commutation ) and these in turn may have other subgroups . these sequences of subgroups define subgroup chains . these chains will be discussed first in terms of the nuclear physics basis @xmath280 , and then in terms of a new basis that is mathematically equivalent but is physically better suited to describing the physics of graphene .    in the nuclear physics basis @xmath281 is the particle number and generates a group @xmath282 , while @xmath283 is proportional to the total angular momentum and generates a group @xmath284 . in the nuclear physics context the total angular momentum and the particle number are expected to be conserved exactly for all physical states . thus one seeks subgroup chains of so(8 ) that end in the subgroup @xmath285 corresponding to charge and angular momentum conservation . three so(8 ) subgroup chains satisfy these conditions .      from eq . ( [ commcoupled1.1 ] ) the quasispin generators @xmath288 close an @xmath289 algebra that is a subalgebra of so(8 ) , and from eq . ( [ commcoupled1.4 ] ) the operators @xmath260 with @xmath290 close an so(5 ) algebra and commute with these su(2 ) quasispin generators . thus , one subgroup of so(8 ) is @xmath291 where so(5 ) is generated by the 10 operators @xmath292 , where @xmath259 takes the @xmath258 values @xmath293 and the _ quasispin group _ @xmath289 is generated by @xmath294 , where @xmath295 . furthermore , the three generators @xmath296 are components of the total angular momentum @xmath297 and generate an @xmath284 subgroup of so(5 ) , and @xmath266 generates a @xmath282 subgroup of @xmath289 corresponding to conservation of charge . thus one subgroup chain is @xmath298 & \supset \underset{\scriptscriptstyle\ {   p^1\ } } { { { \rm so(3)}}_l}\times\underset{\scriptscriptstyle \ { s_0   \ } } { { { { \rm u(1)}}_{{\scriptstyle\rm c } } } } , \label{chains1.2n}\end{aligned}\ ] ] where the generators of each subgroup are indicated in brackets below the subgroup and the product group on the last line corresponds to conservation of total angular momentum and particle number . the groups su(4 ) and so(6 ) share the same lie algebra . in nuclear physics it is more common to refer to this group as so(6 ) , but to maintain a parallel with the ensuing treatment of graphene it will be labeled su(4 ) in the present discussion . from eq . ( [ commcoupled1.6 ] ) , the 16 operators @xmath300 are closed under commutation , corresponding to @xmath301 where the generator of @xmath282 is @xmath302 and the 15 operators @xmath303 are the generators of su(4 ) . furthermore , the subset of @xmath304 with odd @xmath256 are generators of the so(5 ) symmetry discussed above and so generate an so(5 ) subgroup of this su(4 ) group . hence a second subgroup chain is @xmath305 & \supset \underset{\scriptscriptstyle \ { p^1 , \ , p^3\}}{{{\rm so(5 ) } } } \times \underset{\scriptscriptstyle \ { s_0\}}{{{{\rm u(1)}}_{{\scriptstyle\rm c } } } } \supset \underset{\scriptscriptstyle\ {   p^1\ } } { { { \rm so(3)}}_l}\times\underset{\scriptscriptstyle \ { s_0   \ } } { { { { \rm u(1)}}_{{\scriptstyle\rm c } } } } , \label{chains1.3n}\end{aligned}\ ] ] where eq .  ( [ p0n ] ) has been used to replace @xmath306 with @xmath266 .      from eqs . ( [ commcoupled1.2 ] ) , ( [ commcoupled1.5 ] ) , and ( [ commcoupled1.6 ] ) , the 21 operators @xmath307 close an @xmath308 subalgebra of so(8 ) and the subset @xmath309 closes an @xmath310 subalgebra of so(7 ) . thus a third subgroup chain is given by @xmath311 the relationships of the nuclear analog so(8 ) subgroup chains described above are summarized in fig . [ fig : dynamicalchains ] . nuclear analog so(8 ) subgroup chains with generators in the coupled representation @xmath312 given by eqs . ( [ multipole1.1 ] ) and ( [ coupled1.4 ] ) and obeying the lie algebra of eq . ( [ commcoupled ] ) . generators are indicated below each group factor and @xmath266 and @xmath313 may be interchanged using eq . ( [ p0n ] ) . the subgroup structure expressed in this basis is in one - to - one correspondence with the so(8 ) fermion dynamical symmetry model @xcite of nuclear structure physics . however , for the description of graphene it is more useful to transform this basis to the new one employed in fig .  [ fig : dynamicalchains_graphene ] using eqs . ( [ ngx ] ) , which gives a more direct physical interpretation of quantities important in graphene physics . ] the subgroup chains in  [ h : subgroupchainsnuclear ] were expressed in terms of the nuclear physics basis @xmath314 . this basis demonstrates the deep connection between graphene quantum hall physics and nuclear structure physics , and is suitable mathematically to describe graphene quantum hall effects . however , it is not well suited physically to interpreting the graphene quantum hall effects for three reasons . ( 1 )  the relationship between the generators of the @xmath185 subgroup in the nuclear physics basis and those of the su(4 ) quantum hall ferromagnetism basis defined in eq .  ( [ algebra1.4 ] ) is not clear , which hinders interpretation of the present results in terms of preceding results found in the graphene literature . ( 2 )  charge and electronic spin are expected to be conserved in graphene ( if the zeeman term is neglected ) , but none of the generators @xmath315 can be interpreted as spin in the application to graphene ( instead , spin is a linear combination of these generators ) . for physical applications it is desirable to employ a basis in which the relevant conservation laws are manifest . ( 3 )  in addition to the exact conservation laws expected for charge and spin in graphene , it is expected on physical grounds that for low - energy excitations the scattering between valleys is strongly suppressed and the difference in electron densities between neighboring valleys should be very nearly conserved . this difference is expressed by the @xmath58 component of the valley isospin @xmath70 , and the corresponding approximate invariance is reflected in a @xmath316 symmetry generated by @xmath70 . but @xmath70 is not proportional to any of the @xmath304 generators ( it is a linear combination of these generators ) , so this approximate invariance is not manifest in the nuclear su(4 ) basis . thus a new basis will be employed for the so(8 ) generators in application to graphene for which the particle number ( charge ) operator @xmath17 or @xmath266 and the 12 pairing operators @xmath317 are retained , but the 15 su(4 ) generators @xmath318 in the nuclear representation are replaced with the 15 su(4 ) generators @xmath319 with @xmath320 defined in the graphene representation given in eq . ( [ algebra1.4 ] ) , @xmath321 \ { { \mathscr{s}}_\alpha,\ ,   t_\alpha,\ ,   n_\alpha,\ , { \pi_{\alphax } } , \,{\pi_{\alphay } } ,   \ , s_0 , \ , s , \,s^\dagger , \,d{^{\vphantom{\dagger}}}_{\mu},\ , d^\dagger_\mu \ } { _ { \mbox{\scriptsizegraphene   so(8 ) } } } . \end{array}\ ] ] the transformation from the @xmath322 generators to the @xmath323 generators is given in appendix [ h : appendixb ] . in the new basis it will be convenient to take as order parameters    @xmath324    [ orderp1.1 ]    where @xmath325 is the number operator counting particles in basis state @xmath326 and the expectation value is taken with respect to the collective wavefunction . physically ( 1 )  the net spin is measured by @xmath327 , which characterizes _ ferromagnetic _ order . ( 2 )  the difference in charge between the a and b sublattices is measured by @xmath328 , which characterizes _ charge density wave _ order . ( 3 )  the difference in spins between the a and b sublattices is measured by @xmath329 , which characterizes _ antiferromagnetic _ ( nel or spin density wave ) order . the order parameters evaluated for the states of fig .  [ fig : s_d_pairs_brillouin ] are displayed in table  [ tb : orderparms ] . lccc    ' '' ''    state & @xmath327 & @xmath328 & @xmath329     +    ' '' ''    @xmath330 & @xmath331 & @xmath331 & @xmath331     + [ ]    @xmath332 & @xmath331 & @xmath333 & @xmath331     + [ ]    @xmath334 & @xmath331 & @xmath103 & @xmath331     + [ ]    @xmath335 & @xmath333 & @xmath331 & @xmath331     + [ ]    @xmath336 & @xmath103 & @xmath331 & @xmath331     + [ ]    @xmath337 & @xmath331 & @xmath331 & @xmath331     + [ ]    @xmath338 & @xmath331 & @xmath331 & @xmath103     + [ ]    @xmath339 & @xmath331 & @xmath331 & @xmath333     +      in the new basis it will be assumed that both charge and spin are exactly conserved for the physical states of the model in the absence of the zeeman term @xmath54 , that the charge and the @xmath58-component of spin are exactly conserved if the zeeman term is included in the hamiltonian , and that @xmath70 is conserved , where appropriate . neglecting the zeeman term , the spin  charge symmetry corresponds to a group structure @xmath340 , where @xmath341 is generated by the spin operators and @xmath282 is generated by the particle number operator . thus one seeks subgroup chains of so(8 ) that end in the subgroup @xmath340 corresponding to charge and spin conservation , and in some of these chains a u(1)@xmath74 subgroup implying conservation of @xmath70 will also be required . the group and subgroup structure in the new basis is illustrated in fig . [ fig : dynamicalchains_graphene ] ,        where seven nontrivial subgroup chains may be identified that begin with so(8 ) and end with the symmetry @xmath342 corresponding to conservation of charge and @xmath58-component of the spin in the magnetic field . @xmath343    [ groupchains ]    where for brevity all u(1 ) factors are dropped in the notation and su(2 ) means su(2)@xmath344 corresponding to conservation of spin . each of these corresponds to a different dynamical symmetry that is realized for particular choices of parameters for the so(8 ) hamiltonian , and that yields _ exact manybody solutions _ using the dynamical symmetry methodology . let us now discuss in more detail three subgroup chains of so(8 ) . the quasispin generators @xmath345 close an @xmath289 algebra that is a subalgebra of so(8 ) , and the operators @xmath346 close an so(5 ) algebra and commute with these su(2 ) quasispin generators . putting this together , one subgroup of so(8 ) is @xmath347 furthermore , the three generators @xmath348 are components of the total spin and generate an @xmath341 subgroup of so(5 ) , and @xmath266 generates a @xmath282 subgroup of @xmath289 corresponding to conservation of charge . thus one subgroup chain is @xmath349 & \supset { { { \rm su(2)}}_{\sigma}}\times{{{\rm u(1)}}_{{\scriptstyle\rm c } } } , \label{chains1.2}\end{aligned}\ ] ] where the product group on the last line corresponds to conservation of spin and charge . this subgroup chain with its corresponding generators is illustrated in fig . [ fig : dynamicalchains_graphene ] . alternatively , so(5 ) may be broken according to the pattern @xmath350 & \supset { { { \rm su(2)}}_{\sigma}}\times   { { { \rm u(1)}}_{{\scriptstyle\rm c } } } , \label{chains1.2b}\end{aligned}\ ] ] which also conserves spin and charge , and is illustrated in fig . [ fig : dynamicalchains_graphene ] .      a @xmath351 subgroup of so(8 ) may be obtained by removing the 12 pairing operators from the so(8 ) generator set . the @xmath352 subgroup is generated by the particle number ( charge ) and the su(4 ) subgroup is generated by the 15 remaining operators , which are defined in eq . ( [ algebra1.4 ] ) in the current basis . there are several options for chains corresponding to further subgroups . ( 1 )  the subset @xmath353 defines generators of the so(5 ) symmetry discussed above and so forms an so(5 ) subgroup of this su(4 ) group . hence one su(4 ) subgroup chain is @xmath354 which is displayed in fig . [ fig : dynamicalchains_graphene ] . ( 2 )  physically , the total spin is conserved . if there is little inter - valley scattering one may also assume the spin within each @xmath2 valley and each @xmath3 valley to be separately conserved , corresponding to a @xmath355 symmetry . thus a second su(4 ) subgroup chain corresponds to @xmath356 & \supset { { { \rm su(2)}}_{\sigma}}\times { { { \rm u(1)}}_{{\scriptstyle\rm c } } } , \label{chainy1.1}\end{aligned}\ ] ] where @xmath316 is generated by @xmath70 . this chain also is displayed in fig . [ fig : dynamicalchains_graphene ] . ( 3 )  finally , one can imagine that su(4 ) is broken into subgroups corresponding to simultaneous conservation of both spin and valley isospin , giving a third su(4 ) subgroup chain @xmath357 as illustrated in fig . [ fig : dynamicalchains_graphene ] .    comparing fig . [ fig : dynamicalchains_graphene ] with fig . [ fig : dynamicalchains_qhf ] , it is apparent that the three @xmath358 subgroup chains defined in this section correspond to the three symmetry - breaking patterns described in  [ h : symmetriesfqhegraphene ] and discussed in refs . @xcite for the su(4 ) quantum hall ferromagnetism model . the portion of the so(8 ) subgroup structure leading to su(4 ) quantum hall ferromagnetism is indicated by the shaded boxes with dashed outer boundaries in fig . [ fig : dynamicalchains_graphene ] . thus it is seen explicitly that the special case corresponding to the @xmath358 subgroup chains of the present model imply the results of refs . @xcite . the 21 operators @xmath359 close an @xmath360 subalgebra of so(8 ) , and the subset of generators @xmath361 close an @xmath362 subalgebra of so(7 ) . thus a third subgroup chain is given by @xmath363 as illustrated in fig . [ fig : dynamicalchains_graphene ] . this subgroup chain is of particular interest because it will define a _ critical dynamical symmetry _ that represents an entire phase exhibiting critical behavior and interpolating between two other phases . let us use the subgroup structure of the preceding section to obtain exact solutions of the correlated many - body problem in these dynamical symmetry limits . the basic idea is to use the casimir invariants of the subgroup chains like those described in  [ h : subgroupchainsgraphene ] and illustrated in fig . [ fig : dynamicalchains_graphene ] to label states . then model hamiltonians constructed only from the casimir invariants of a single subgroup chain permit analytical solution of the effective schrdinger equation in that symmetry limit @xcite . specifically , if a hamiltonian @xmath364 can be expressed as a function of the casimir invariants of some subgroup chain @xmath365 , where the @xmath366 represent casimir operators for the groups @xmath367 , then the system is said to possess a dynamical symmetry associated with the subgroup chain . the discussion will be simplified by restricting to the lowest - order casimir invariant for each group , which corresponds physically to omitting @xmath17-body interactions with @xmath368 . elementary properties of lie groups then permit the eigenvalues @xmath369 and eigenfunctions @xmath370 of this hamiltonian to be expressed in closed form as @xmath371 where the @xmath372 stand for the quantum numbers required to specify the irreducible representations ( irreps ) of the groups @xmath367 . the physical properties of the corresponding states can then be elucidated by using the methods of lie groups and lie algebras to evaluate matrix elements for observables . in this way , one generally finds that the dynamical symmetries associated with subgroup chains of some highest symmetry define collective ( emergent ) states that correspond to particular patterns of spontaneous symmetry breaking .      in terms of the generators ( [ coupled1.4 ] ) , ( [ multipole1.1 ] ) , and ( [ algebra1.4 ] ) , the quadratic casimir operator @xmath373 for the so(5 )   subgroup is @xmath374 where @xmath375 . the corresponding casimir operator for the @xmath376 subgroup is @xmath377 the quadratic casimir operator for the so(7 )  subgroup is @xmath378 the casimirs for the su(2 ) subgroups are @xmath379 and for the u(1 ) charge subgroup the invariant operator is trivially the single generator , @xmath380 or some power of it . finally , the quadratic casimir operator for the full so(8 ) group may be expressed as @xmath381 the casimir operators that appear in each of these subgroups chains and the relevant quantum numbers labeling the states for each dynamical symmetry are summarized in table  [ tb : so8subgroups ] .    llllll    ' '' ''    group & dim & generators & quantum numbers@xmath382 & casimir operators @xmath383 & casimir eigenvalues@xmath384     +    ' '' '' @xmath385 & @xmath386 & @xmath387 & @xmath388 & @xmath389 & @xmath390     + [ ]    @xmath391 & @xmath392 & @xmath393 or & @xmath394 & @xmath395 & @xmath396     + [ ]    &  & @xmath397 &  &  & + [ ]    @xmath398 & @xmath399 & @xmath400 or & @xmath401 & @xmath402 & @xmath403     + [ ]     &   & @xmath404 &   & @xmath405 & + [ ]    u(4 ) & @xmath406 & @xmath407 & @xmath408 & @xmath409 & @xmath410     + [ ]    @xmath411 & @xmath412 & @xmath413 or & @xmath414 & @xmath415 & @xmath416     + [ ]    &  & @xmath417 &  & @xmath418 & + [ ]    @xmath289 & @xmath108 & @xmath419 & @xmath9 & @xmath420 & @xmath421     + [ ]    @xmath341 & @xmath108 & @xmath422 & @xmath423 & @xmath424 & @xmath425     + [ ]    @xmath426 & @xmath108 & @xmath63 & @xmath86 & @xmath427 & @xmath428     + [ ]    @xmath282 & @xmath103 & @xmath266 & @xmath17 & @xmath266 & @xmath429     +   +   +   +   +   +   +   +    in the next section these results for the casimir operators will be used to construct the most general hamiltonian permitted in the truncated collective space for specific dynamical symmetries . as has been seen , in a particular dynamical symmetry limit the most general hamiltonian can be constructed from a sum of casimir invariants for the groups contained in the corresponding subgroup chain . for so(8 ) dynamical symmetries the most general hamiltonian in the absence of the zeeman term is represented by the linear combination @xmath430 where @xmath431 is assumed constant in the symmetry limit , the casimir operators @xmath373 have been summarized in the preceding section , and @xmath432 does not appear explicitly because it has been eliminated by the constraint @xcite : @xmath433 hamiltonians representing specific dynamical symmetry limits then correspond to particular choices of the coefficients @xmath434 in this general expression . it may be shown that the most general hamiltonian can also be expressed in the compact form ( see eqs . ( 4.1 ) of ref . @xcite ) @xmath435 where @xmath436 is assumed constant in a symmetry limit and where the different dynamical symmetry limits correspond to specific choices for the values of the parameters @xmath437 , @xmath438 , and @xmath439 . the last term is expressed in terms of the @xmath304 from the nuclear basis . it can be converted to the graphene basis by inverting eqs . ( [ ngx ] ) of appendix [ h : appendixb ] to solve for the @xmath304 . let us now discuss each of the so(8 ) dynamical symmetries and their physical interpretations . for brevity , let us refer to    1 .   the dynamical symmetry structure associated with eqs .  ( [ chains1.2 ] ) and ( [ chains1.2b ] ) as the _ @xmath440 dynamical symmetry , _ 2 .   the dynamical symmetry structure associated with eqs . ( [ chains1.3])([chainy1.2 ] ) as the _ su(4 ) dynamical symmetry , _ and 3 . the dynamical symmetry structure associated with eq . ( [ chains1.4 ] ) as the _ so(7 ) dynamical symmetry . _    initially the role of the zeeman term ( which would break the full spin symmetry down to conservation of its @xmath58 component ) will be ignored and it will be assumed that the chains end in the subgroup @xmath441 corresponding to the physical requirement that spin and charge be conserved exactly . the dynamical symmetry chains given in eqs . ( [ chains1.2])([chains1.2b ] ) and illustrated in fig . [ fig : dynamicalchains_graphene ] correspond to two alternative ways of choosing subgroups of @xmath442 : @xmath443{45}{$\supset$}}}$ } \ { { { \rm su(2)}}_{\sigma}}\times   { { \rm su(2)}}{_{\mbox{\scriptsizep}}}\ \raisebox{-0.8ex}{${\mathbin{\rotatebox[origin = c]{-45}{$\supset$}}}$ } \\   { { \rm so(5)}}\times{{\rm su(2)}}{_{\mbox{\scriptsizep } } } \hphantom{{\mathbin{\rotatebox[origin = c]{45}{$\supset$}}}\ { { { \rm su(2)}}_{\sigma}}\times{{{\rm su(2)}}_{\tau}}{\mathbin{\rotatebox[origin = c]{-45}{$\supset$ } } } } \ \   { { { \rm su(2)}}_{\sigma}}\times   { { { \rm u(1)}}_{{\scriptstyle\rm c}}}\\ \hphantom{{{\rm so(5)}}\times{{{\rm su(2)}}_{\tau}}\ } \raisebox{1.0ex}{${\mathbin{\rotatebox[origin = c]{-45}{$\supset$}}}$ }    \makebox[\widthof{$\ { { { \rm su(2)}}_{\sigma}}\times{{{\rm su(2)}}_{\tau}}\ $ } ] [ c]{${{\rm so(5)}}\ \times \   { { { \rm u(1)}}_{{\scriptstyle\rm c}}}$ } \raisebox{1.0ex}{${\mathbin{\rotatebox[origin = c]{45}{$\supset$}}}$ } . \end{array } \label{so5dyn1.1}\ ] ] in the upper branch of the middle step the so(5 ) symmetry is broken to its @xmath444 spin subgroup , with @xmath289 unbroken . physically this corresponds to conservation of the spin associated with the @xmath445 angular momentum algebra and the pseudospin associated with the @xmath446 pair algebra , but not the full so(5 ) symmetry . in the lower branch of the middle step , the so(5 ) symmetry remains intact and @xmath143-pair pseudospin @xmath444 is broken to @xmath447 charge . in the final subgroup of both chains , only the spin and charge remain as conserved quantities . the hamiltonian in the @xmath448 dynamical symmetry limit corresponds to eq . ( [ mostgeneralham ] ) with the restriction that @xmath449 , @xmath450 where terms that are constant within a given representation have been omitted . the most general so(8 )  state in the @xmath177 collective subspace is given by eq . ( [ collective1.1 ] ) and corresponds to a superposition of @xmath143 and @xmath231 pairs . schematically , @xmath451 with @xmath228 and @xmath452 defined in eq . ( [ coupled1.4 ] ) and @xmath453 , where @xmath178 is the total pair number . on the other hand , the most general states corresponding to the various subgroup chains illustrated in fig .  [ fig : dynamicalchains_graphene ] correspond to pair superpositions having specific constraints on the relative contribution of @xmath143 and @xmath454 pairs . the collective subspace for the @xmath455 subgroup of so(8 ) is of the form ( see eq .  ( 9.13 ) of ref . @xcite ) @xmath456 implying that it is a superposition of @xmath143 pairs . conceptually , the wavefunction of the @xmath455 subgroup chain for @xmath177 is obtained from the most general state in the collective subspace by converting all of its @xmath454 pairs to @xmath143 pairs .    as was seen in  [ sh : coupledpairrep ] and fig . [ fig : s_d_pairs_brillouin ] , the @xmath143 and @xmath454 pairs correspond to coherent superpositions of particular electronic distributions in spin and valley pseudospin . thus , specific @xmath143 and @xmath454 pair content for so(8 ) dynamical symmetry subgroup chains implies specific collective modes associated with coherent distribution of the electrons in spin and valley space . it has been noted above that in the @xmath457 dynamical symmetry limit the ground states correspond to a superposition of @xmath143 pairs . the nature of this collective state may be inferred from fig . [ fig : s_d_pairs_brillouin ] and is illustrated in fig . [ fig : s_d_pairs_brillouin_so5 ] .     the collective state corresponding to the @xmath458 ground state . the state is an antisymmetric combination of a spin density wave with spin - up on the a lattice sites and spin - down on the b sites , and a spin density wave with spin - up on the b sites and spin - down on the a sites . ] the su(4 ) dynamical symmetry corresponds to the three so(8 ) subgroup chains @xmath459 that were introduced in eqs . ( [ chains1.3])([chainy1.2 ] ) and fig . [ fig : dynamicalchains_graphene ] . as already noted , these three dynamical symmetry chains are in one to one correspondence with the explicit symmetry breaking patterns that have been identified for su(4 ) quantum hall ferromagnetism . the most general su(4 ) wavefunction for @xmath178 pairs in the @xmath177 collective so(8 ) subspace is given by @xcite @xmath460^p { |0\rangle } ,   \label{su4pairstate}\ ] ] which corresponds physically to a restriction of the general so(8 ) wavefunction ( [ so8pairstate ] ) to a specific superposition of @xmath143 and @xmath454 pairs . the wavefunction for the @xmath362 subgroup of su(4 ) is given by eq . ( [ so5pairstate ] ) . the wavefunction of the parent su(4 ) group is a superposition of @xmath143 and @xmath454 pairs but the @xmath448 subgroup has a wavefunction containing only @xmath143 pairs . the so(7 ) dynamical symmetry corresponds to the so(8 ) subgroup chain @xmath461 that was introduced in eq . ( [ chains1.4 ] ) and displayed in fig . [ fig : dynamicalchains_graphene ] . the hamiltonian in the so(7 ) dynamical symmetry limit corresponds to eq . ( [ mostgeneralham ] ) with the restriction that @xmath462 , @xmath463 where terms have been dropped that are constant within a given representation .    from the nuclear physics analog so(8 ) symmetry @xcite , one may surmise that so(7 ) will play the role of a _ critical dynamical symmetry _ interpolating smoothly between the collective states corresponding to the @xmath411 dynamical symmetry and the collective states corresponding to @xmath398 dynamical symmetry . such critical dynamical symmetries have been discussed previously in both nuclear physics @xcite and for the strongly - correlated electrons leading to cuprate and iron - based high - temperature superconductivity @xcite . they may be viewed as the generalization of a quantum critical point to an entire _ quantum critical phase , _ and may represent a fundamental organizing principle for quantum critical behavior . the physical implications of this so(7 ) critical dynamical symmetry for graphene quantum hall physics will be discussed further below . the dynamical symmetry limits discussed above represent special solutions resulting from particular choices of the coupling parameters appearing in the hamiltonian . for arbitrary choices of the coupling parameters the solutions will correspond generally to superpositions of the different symmetry - limit solutions and will not have exact analytical forms . in this more general case it is quite feasible to obtain solutions numerically , since the collective subspace is highly truncated relative to the full hilbert space . however , there is a powerful alternative approach : the _ generalized coherent state approximation , _ which permits _ analytical solutions _ for arbitrary choices of the coupling parameters in the hamiltonian .    for the so(8 ) lie algebra introduced in this paper for graphene the gilmore  perelomov algorithm @xcite may be implemented to obtain solutions in terms of a set of generalized coherent states . these solutions represent the most general hartree  fock  bogoliubov theory that can be formulated in the space , subject to a dynamical symmetry constraint @xcite . the solutions of this symmetry - constrained hartree  fock  bogoliubov ( schfb ) theory correspond to determining the stable points of energy surfaces , which represent the coherent - state expectation values of the effective hamiltonian on the coset space . thus the coherent state solutions also represent a microscopically - derived implementation of ginzburg  landau theory . these coherent state solutions are uniquely well suited to study the interplay of competing spontaneous symmetry breaking in determining the ground state of the system and its properties .      the coherent states associated with the full set of subgroup chains in fig . [ fig : dynamicalchains_graphene ] will be discussed in future work . in this paper the power of the method will be illustrated succinctly by restricting to the coherent states associated with the subgroup chains of so(8 ) that contain the so(5 ) subgroup , as illustrated in fig . [ fig : dynamicalchains_coherentso5 ] . 0pt 0.20 subgroup chains included in the present coherent state analysis . thus the corresponding coherent state solutions will represent a superposition of the symmetry - limit solutions for the @xmath464 { { \rm so(8)}}\supset{{\rm su(4)}}\supset{{\rm so(5)}}\qquad { { \rm so(8)}}\supset{{\rm so(7)}}\supset{{\rm so(5)}}\end{array}\ ] ] dynamical symmetries . these coherent state solutions will be seen to have the following properties :    1 . the @xmath391 dynamical symmetry will play the role of a critical dynamical symmetry interpolating between @xmath411 and @xmath398 symmetry - limit solutions . because of fundamental symmetries obeyed by the wavefunction , the coherent state solutions may be parameterized in terms of a single collective parameter @xmath465 that governs the mixture of the @xmath143 and @xmath454 pairs defined in eq .  ( [ coupled1.4 ] ) contributing in the ground state . the collective parameter @xmath465 may also be interpreted physically in terms of the pair configurations displayed in fig . [ fig : s_d_pairs_brillouin ] . the so(8 ) coherent states corresponding to the symmetry structure in fig . [ fig : dynamicalchains_coherentso5 ] have been developed previously in ref . @xcite for nuclear physics applications and will be adapted extensively to development of the present formalism . let us now consider the energy surfaces that may be computed from the coherent states , which link the so(8 ) solutions to ginzburg  landau theory . within the coherent state formalism the ground state energy may be determined through the variational requirement @xmath466 where @xmath467 is the coherent state , @xmath468 is the so(8 ) hamiltonian , @xmath469 and the coefficients @xmath437 and @xmath470 are functions of the effective interaction . for the dynamical symmetry chains having so(5 ) as a subgroup the energies take the general form @xcite @xmath471 , \label{generalenergyformula}\ ] ] where the group - dependent parameters @xmath472 , @xmath473 , @xmath474 , @xmath475 , and @xmath476 are given in table  [ tb : stiffnessparameters ] . cccccc    ' '' ''    g & @xmath472 & @xmath124 & @xmath477 & @xmath478 & @xmath479     +    ' '' ''    su(2 ) & @xmath437 & @xmath480 & @xmath481 & @xmath482 & @xmath483     + [ ]    so(5 ) & @xmath484 & @xmath485 & @xmath486 & @xmath487 & @xmath488     + [ ]    su(4 ) & @xmath489 & @xmath490 & @xmath491 & @xmath492 & @xmath331     + [ ]    so(7 ) & @xmath438 & @xmath493 & @xmath494 & @xmath495 & @xmath496     +    our primary interest in this discussion is in the ground state properties of graphene in a strong magnetic field . the ground states in the coherent state approximation at fixed @xmath497 will be given by those values of @xmath498 that correspond to minima of the energy surface @xmath499 . these are determined by the values of @xmath465 satisfying @xmath500 evaluating these constraints for the energy surfaces ( [ generalenergyformula ] ) , one finds that the minima @xmath501 are given by @xcite @xmath502 the coherent state energy surfaces for the @xmath448 , so(7 ) , and su(4 ) symmetry limits computed from eq . ( [ generalenergyformula ] ) using the entries in table  [ tb : stiffnessparameters ] are shown as functions of @xmath465 for several values of the fractional occupation @xmath160 in fig . [ fig : coherentenergysurfaces_graphene ] .     coherent state energy surfaces as a function of the order parameter @xmath465 for three of the so(8 ) dynamical symmetry limits . the curves are labeled by the fractional occupation @xmath503 defined in eq . ( [ fillingfraction ] ) , where @xmath17 is the particle number and @xmath147 is the maximum number of particles that can be accommodated in the landau level . the formalism is particle  hole symmetric so curves for fractional occupations with @xmath504 are equivalent to those shown above , but with the fractional occupation counted in terms of number of holes . for example , the curves for @xmath505 and @xmath506 are equivalent . for the direct product group @xmath507 the energy surface is the sum of contributions to eq . ( [ generalenergyformula ] ) from su(2 ) and so(5 ) . as will be demonstrated later , the su(2 ) contribution typically dominates that of so(5 ) . ] there one sees that indeed the minima for the @xmath448 and so(7 ) limits are at @xmath508 , and the minima for the su(4 ) limit are at @xmath509 . although the minimum energies for both the @xmath448 and so(7 ) limits are consistent with @xmath510 , fig . [ fig : coherentenergysurfaces_graphene ] shows that these symmetries differ fundamentally in the localization of the minimum . for @xmath448 the energy surface has a deep minimum at @xmath510 but for so(7 ) the energy surface is very flat around @xmath508 , with a broad range of @xmath465 giving essentially the same ground state energy . this highly - degenerate so(7 ) ground state has significant physical implications that will be discussed further below . the so(8 ) generalized coherent state is equivalent to the hartree  fock  bogoliubov ( hfb ) approximation subject to a symmetry constraint . since hfb is a bcs - type approximation married to a hartree  fock mean field , its solutions correspond to symmetry - breaking intrinsic states . in particular , the bcs - like state conserves the physical particle number only on average , and the hartree  fock mean field may break both translational and rotational invariance . let us address these issues for the so(8 ) coherent state .      the fractional uncertainty in electron number @xmath511 for the so(8 ) coherent state is given by @xcite @xmath512 where @xmath513 is the value of @xmath465 at the minimum energy , given by eq . ( [ betaminima ] ) in the symmetry limits . expressing eq . ( [ particlenumberfluctuation ] ) in terms of the fractional occupation @xmath514 , in the @xmath448 and su(4 ) limits , respectively , one obtains @xmath515_{{\rm so}_5\times{{\rm su(2 ) } } } \hspace{-5pt } = \sqrt{\frac{1-{f}}{{f}\omega } } \qquad \left[\frac{\delta n}{n}\right]_{{\rm su}_4 } \hspace{-5pt } =   \sqrt{\frac{1 - 2{f}}{{f}\omega } } , \label{particlenumberfluctuation3}\ ] ] from these results one may notice two important things .    1 . the fluctuation in particle number is large at low degeneracy @xmath128 but decreases with increasing @xmath128 . if su(4 ) symmetry is realized @xmath516 decreases with increasing @xmath160 and _ vanishes identically _ at @xmath167 for any @xmath128 , which corresponds to the fractional occupation for the ground state of undoped graphene . thus it is expected that the current theory applied to graphene has negligible particle number fluctuation @xmath516 in the su(4 ) limit . in the @xmath448 limit the particle number fluctuation @xmath516 remains finite for all @xmath160 but it becomes very small as @xmath128 becomes large , particularly near @xmath517 . thus it too may be neglected in the large-@xmath128 limit . comparison with table [ tb : degeneracytable ] suggests that graphene quantum hall experiments involve sufficient degeneracy that particle number fluctuation in the coherent state solution is not significant . the coherent state approximation represents a mean field localized in spatial position and orientation , so it violates translational and rotational invariance . however , because the crystal is generally macroscopic , the net violation of these symmetries may be expected to be negligible . one concludes that for applications of coherent state methods to graphene , violations of particle number conservation , rotational invariance , and translational invariance are negligible in realistic systems . the generalized coherent state method has been used above to calculate total energy surfaces for quantum hall states in graphene , but one also may use the coherent state wavefunctions and appropriate operators to calculate matrix elements of other relevant observables . this section addresses the nature of the wavefunction and the matrix elements that can serve as order parameters . a significant consequence of the so(8 ) dynamical symmetry structure displayed in fig .  [ fig : dynamicalchains_coherentso5 ] is that the phases may be distinguished in terms of a single order parameter and its fluctuations , which may be taken to be @xmath465 . let us now characterize in more depth the physical meaning of this order parameter . in  [ sh : orderparameters ] an antiferromagnetic order parameter @xmath329 was defined . in the coherent state approximation the onset of af order is signaled by an energy - surface minimum at a finite value of @xmath465 . because @xmath518 [ see eq . ( [ ngx9 ] ) ] , the antiferromagnetic order parameter @xmath329 is related to the coherent state af order parameter @xmath465 by the intrinsic state matrix element of @xmath519 @xcite @xmath520 where @xmath489 is the coupling strength for the @xmath521 term in the hamiltonian . each value of @xmath465 corresponds to a unique value of @xmath329 , so @xmath465 is a measure of antiferromagnetic order . the location of the maxima may be obtained by setting the derivative with respect to @xmath465 of eq . ( [ intrinsic1.1 ] ) equal to zero , which yields that @xmath522 for a given @xmath17 occurs at a @xmath465 of @xmath523 where eq . ( [ betaminima ] ) was used to make the identification in the last step . thus , if @xmath524 the maximum value of @xmath329 maps to a value of @xmath465 that corresponds to a minimum of the energy surface ( a ground state ) in the su(4 ) limit . substituting eq . ( [ betaminima ] ) for @xmath465 in eq . ( [ intrinsic1.1 ] ) , for @xmath524 ground states the af order parameter @xmath329 depends on the electron number @xmath17 as @xmath525 the so(8 ) model is particle  hole symmetric so @xmath17 or @xmath526 count electrons up to half filling and holes for greater than half filling . hence the maximum af collectivity occurs for half filling of the single valence landau level .      as was discussed in  [ h : conservationlaws ] , the coherent state wavefunction corresponding physically to @xmath527 pairs conserves particle number only on average and so is a superposition of terms having different numbers of pairs . in eq . ( 5.27 ) of ref . @xcite the so(8 ) coherent state is decomposed into terms of definite pair number @xmath528 according to    @xmath529    where @xmath530 denotes an intrinsic state with order parameter @xmath465 and closed forms for @xmath531 and @xmath532 are given in ref . @xcite . according to eq . ( 5.28 ) of ref . @xcite , the values of @xmath533 that correspond to the minima of the potential energy surface at @xmath534 for the @xmath448 limit and @xmath535 for the @xmath411 limit [ see eq . ( [ betaminima ] ) ] , respectively , are @xmath536 so the so(8 ) coherent state wavefunction ( [ coherentfn1.1 ] ) in the @xmath448 and su(4 ) limits , respectively , becomes @xmath537 \displaystyle { |{\rm su}_4\rangle } = \sum _ { p } { c_p}\left ( s^\dagger \pm d_0^\dagger \right)^p { |0\rangle } = 2\sum _ { p } { c_p}\left ( q_\pm^\dagger \right)^p { |0\rangle } , \end{array } \label{coherentfn1.4}\ ] ] where eq . ( [ qdef ] ) was used .    as discussed in  [ h : conservationlaws ] , fluctuations in particle number are negligible in the large-@xmath128 limit for so(8 ) coherent states , implying that the summations in eqs . ( [ coherentfn1.4 ] ) become dominated by terms with @xmath538 . thus for large @xmath128 the coherent state wavefunctions are well approximated up to a normalization by    @xmath539    [ largeomegwf ]    as seen from table  [ tb : orderparms ] , the su(4 ) state of eq . ( [ largeomegwf.b ] ) is a coherent superposition of @xmath540 or @xmath541 pairs , each contributing vanishing ferromagnetic order @xmath327 and charge density wave order @xmath328 , but non - zero af order @xmath329 . conversely , the @xmath448 state of eq . ( [ largeomegwf.a ] ) is a coherent superposition of @xmath143 pairs , each with vanishing @xmath327 , @xmath328 , and @xmath329 . thus the @xmath448 and su(4 ) limits of the so(8 ) symmetry are distinguished by the order parameter @xmath329 , which is zero in the @xmath448 state and is non - zero in the su(4 ) state . equivalently , the coherent state order parameter @xmath465 vanishes in the pure @xmath448 limit and is equal to @xmath542 in the pure su(4 ) limit [ see eq . ( [ betaminima ] ) and fig . [ fig : coherentenergysurfaces_graphene ] ] . equation ( [ generalenergyformula ] ) depends only on even powers of @xmath465 so the sign for the two possible spontaneously broken symmetry solutions does not affect the energy .    for undoped graphene the fermi surface corresponds to the @xmath543 curves of fig . [ fig : coherentenergysurfaces_graphene ] . these are shown in fig . [ fig : coherentenergysurfaces_graphenegs ] for the @xmath448 , so(7 ) , and su(4 ) limits . ground - state energy surfaces in coherent state approximation for three of the so(8 ) dynamical symmetry limits . ( a )  the @xmath448 limit . ( b )  the so(7 ) limit . ( c )  the su(4 ) limit . the diagrams indicate schematically the corresponding wavefunctions , as suggested by eqs . ( [ coherentfn1.1])([largeomegwf ] ) and fig . [ fig : s_d_pairs_brillouin ] . ]    * the @xmath448 limit in fig . [ fig : coherentenergysurfaces_graphenegs](a ) has its minimum @xmath513 at @xmath544 . it corresponds to states of the form ( [ largeomegwf.a ] ) , with vanishing expectation values of @xmath545 , @xmath70 , and @xmath546 . * the su(4 ) limit of fig . [ fig : coherentenergysurfaces_graphenegs](c ) corresponds to states of the form ( [ largeomegwf.b ] ) , with minima located at @xmath547 when the symmetry is broken spontaneously by choosing one of these possibilities ( they are equivalent since the energy depends only on even powers of @xmath465 ) , the resulting state has @xmath548 ( no spin or isospin order ) , but @xmath549 ( spin density wave or af order ) . * the so(7 ) limit of fig . [ fig : coherentenergysurfaces_graphenegs](b ) corresponds to a critical dynamical symmetry that interpolates between the @xmath448 and su(4 ) states through critical fluctuations in the antiferromagnetic order . thus the so(8 ) dynamical symmetry limits illustrated in fig . [ fig : coherentenergysurfaces_graphenegs ] represent a rich set of collective states that can be distinguished by the expectation value and fluctuations associated with the order parameter @xmath465 . the so(8 ) coherent state solution can be used to study transitions among the phases defined in fig .  [ fig : coherentenergysurfaces_graphenegs ] . for the @xmath177 space ( no broken pairs ) assumed here , @xmath550 yields a constant that is neglected and eq . ( [ coherent1.12 ] ) may be expressed as @xmath551 the last two terms yield constants when evaluated in a given representation , and @xmath552 is found to contribute negligibly to the total energy compared with @xmath553 . therefore , it will be instructive to set @xmath554 and study the approximate so(8 ) hamiltonian @xmath555 from eqs . ( [ casimir_so5])([casimir_so8 ] ) , one finds that @xmath556 { \langles^\dagger s\rangle } + { \langle p^2{\hspace{-0.1em}\cdot\hspace{-0.1em}}p^2\rangle } \sim { \langle{c_{\scriptscriptstyle{{\rm so(7)}}}}\rangle } \end{array}\ ] ] if constants are neglected . thus the model hamiltonian ( [ modelhamiltonian ] ) may be tuned to favor the @xmath448 , so(7 ) , or su(4 ) phases by varying the ratio of the coupling parameters @xmath437 and @xmath489 .      to study the quantum phase transitions of the so(8 ) model with the approximate hamiltonian ( [ modelhamiltonian ] ) , it is convenient to define a parameter @xmath557 and to rewrite eq . ( [ modelhamiltonian ] ) as @xmath558 thus the value of @xmath559 tunes the hamiltonian ( [ modelhamiltonian2 ] ) between @xmath507 and su(4 ) phases via an intermediate so(7 ) phase exhibiting quantum critical behavior .    1 . if @xmath560 the ground - state energy surface is approximated by fig . [ fig : coherentenergysurfaces_graphenegs](a ) , with a minimum at @xmath544 , no antiferromagnetic order , and @xmath507 symmetry . 2 .   if @xmath561 the ground - state energy surface is approximated by fig . [ fig : coherentenergysurfaces_graphenegs](c ) , with an energy minimum at @xmath562 implying @xmath411 symmetry and antiferromagnetic order . 3 .   if @xmath563 , the ground - state energy surface is approximated by fig . [ fig : coherentenergysurfaces_graphenegs](b ) and the system exhibits so(7 ) critical dynamical symmetry , with large fluctuations in the af order parameter @xmath465 . let us now use the hamiltonian ( [ modelhamiltonian2 ] ) to study the quantum phase transitions and spontaneously broken symmetry of the so(8 ) model .      if terms involving @xmath564 and @xmath565 are ignored ( as justified above ) , eq . ( [ generalenergyformula ] ) with the parameters in table  [ tb : stiffnessparameters ] imply that @xmath566 hence the energy surfaces corresponding to the hamiltonian ( [ modelhamiltonian2 ] ) may be expressed as @xmath567 the variation of the energy surface computed from eq . ( [ coherent1.13d ] ) with the control parameter @xmath568 for half filling ( ground state for undoped graphene ) is shown in fig . [ fig : coherentenergysurfaces_control2d ] ;     quantum phase transitions with coupling strength as control parameter : coherent state energy surface as a function of the control parameter @xmath569 for a degeneracy parameter @xmath570 and a fractional occupation @xmath571 . the solid green curves from @xmath572 correspond to approximate @xmath448 symmetry , the solid blue curves for @xmath573 correspond to approximate su(4 ) symmetry . the dashed red curves near @xmath563 correspond to an approximate so(7 ) symmetry mediating the quantum phase transition from @xmath448 to su(4 ) symmetry . ]    by tuning the control parameter from @xmath574 to @xmath561 , one sees that the system undergoes a quantum phase transition near @xmath575 from an approximate @xmath457 state with energy minimum at @xmath534 [ see fig .  [ fig : coherentenergysurfaces_graphenegs](a ) ] to an approximate su(4 ) state having energy minima at @xmath576 [ see fig . [ fig : coherentenergysurfaces_graphenegs](c ) and eq . ( [ betaminima ] ) ] . for @xmath563 the system has an approximate so(7 ) dynamical symmetry [ see fig .  [ fig : coherentenergysurfaces_graphenegs](b ) ] , with no well - defined minimum for the energy as a function of @xmath465 , but with large fluctuations in @xmath465 implied by a highly - degenerate ground state .    for fixed values of the coupling parameters @xmath437 and @xmath489 , phase transitions may be mediated by changing the particle occupancy . [ fig : coherentenergysurfaces_vs_n_2d ] illustrates for different values of @xmath497 at fixed @xmath577 . one sees that as the particle number is increased the system makes a transition from approximate @xmath448 symmetry with @xmath534 to su(4 ) symmetry with @xmath524 through a critical so(7 ) symmetry for which the energy is highly degenerate in @xmath465 . in the coherent state approximation the dynamical symmetry structure of fig .  [ fig : dynamicalchains_coherentso5 ] has been examined and not the full group structure given in fig .  [ fig : dynamicalchains_graphene ] . for the group chains that contain the so(5 ) subgroup , the only physical implication is to omit the physical effect of zeeman splitting from dynamical symmetry hamiltonian [ which would break the so(5 ) subgroup down into a u(1 ) subgroup generated by the @xmath58 component of the physical spin ] . our primary concern in this discussion is the structure associated with the @xmath5 landau level , for which the effect of the zeeman term is expected to be small ( see  [ h : applicationsgraphene ] ) . thus , one may view the effect of the zeeman term as a perturbation on the results obtained thus far that will act only on the spin part of the wavefunction . as kharitonov @xcite has already discussed , the competition of the zeeman term with the valley interactions will convert the antiferromagnetic solution into a canted antiferromagnetic solution . quantum phase transitions with particle number as control parameter : coherent state energy surfaces for different filling fraction @xmath578 at fixed values for @xmath437 and @xmath489 with @xmath579 . the solid green curves for @xmath580 correspond to approximate @xmath448 symmetry , the solid blue curves for @xmath581 correspond to approximate su(4 ) symmetry . curves near @xmath582 ( dashed red ) correspond to an approximate so(7 ) symmetry mediating the quantum phase transition from @xmath448 to su(4 ) symmetry . ] placing a strong magnetic field on graphene leads to high level degeneracy at energies corresponding to the quantized landau levels . in general , interacting electronic systems with high level degeneracy near the fermi surface can produce ( even for relatively weak interactions ) a rich variety of collective states that differ qualitatively from the non - interacting ground state . these states correspond to a spontaneous breaking of the symmetry and generally can not be obtained through small perturbations of the weakly - interacting ground state since they are emergent in nature . the su(4 ) symmetry of quantum hall ferromagnetism gives rise to a ground state symmetry reflected in the su(4)-symmetric hamiltonian ( [ so5_1.1a ] ) and a possible symmetry breaking structure that has been outlined in  [ h : symmetriesfqhegraphene ] . however , the symmetry - breaking patterns illustrated in fig .  [ fig : dynamicalchains_qhf ] represents _ perturbations around the symmetric ground state ( explicit symmetry breaking ) . _ they can not capture the nature of these possible spontaneously - broken symmetry states , since the broken - symmetry states may differ fundamentally from the possible states in fig . [ fig : dynamicalchains_qhf ] : the states corresponding to the most general linear combination ( [ upairs1.8 ] ) represent a complex superposition of many su(4)-symmetric components and generally can not be classified by pure or any simple linear combination of su(4 ) irreducible representations .    since the nature of the broken symmetry states can not be studied directly within the su(4 ) framework because they are unlikely to be anywhere near eigenstates of an su(4)-symmetric hamiltonian , the broken - symmetry states have typically been studied numerically , or by effective field theory methods . however , as has been shown , the kinds of collective configurations that have been proposed as candidates for low - lying broken - symmetry states in graphene ( see fig .  [ fig : s_d_pairs_brillouin ] and @xcite ) bear strong resemblance to _ eigenstates _ of so(8 ) dynamical symmetry chains . thus , the present so(8 ) symmetry holds the promise of providing _ analytical _ solutions for possible broken - symmetry states in graphene . this is the most important result of the present paper .    at specific filling factors the ground state of graphene will be determined by the competition among the su(4 ) symmetry breaking terms . the most obvious su(4)-anisotropic effect is the zeeman term , which favors a spin - polarized state , but the graphene sublattice structure introduces additional interactions that favor ground states without spin polarization that are characterized by spin density wave or charge density wave order at the lattice scale . the competition between zeeman - term spin polarization and the lattice - scale polarizations can be studied by changing the in - plane component of the magnetic field relative to the perpendicular component , since this changes the zeeman energy but not the orbital energies @xcite . such studies indicate that for the higher - energy landau levels the zeeman term is dominant , producing spin ferromagnets that have skyrmionic excitations at half filling , but in the @xmath5 landau level the lattice - scale interactions dominate the zeeman interaction and drive the system into a spin - unpolarized state @xcite . the remainder of this discussion will concentrate on these spin - unpolarized collective states that are candidates for the ground state in the @xmath5 landau level for charge - neutral graphene , with the zeeman interaction viewed as a perturbation on a collective structure that is dominated by lattice - scale interactions . consideration in this paper will be restricted further to those states that can arise from the dynamical symmetries of fig . [ fig : dynamicalchains_graphene ] that contain the so(5 ) subgroup ( those displayed in fig . [ fig : dynamicalchains_coherentso5 ] ) . it has been shown that so(8 ) describes analytically a number of spontaneously - broken - su(4 ) candidates for the states observed in modern experiments such as those described in refs . these solutions provide a spectrum of excited states as well as ground states . the excited states will not be discussed here , except to note that all ground state solutions have a gap to electronic and collective excitations and so correspond to incompressible states . the general theory to be discussed in forthcoming papers can accommodate fm , cdw , and af states , but for the dynamical symmetries containing so(5 ) that were the focus here , all solutions may be classified by a _ single parameter _ @xmath465 measuring af order : su(4 ) states have finite af order but no cdw or fm order , @xmath448 states have no af , cdw , or fm order , and so(7 ) states correspond to a critical dynamical symmetry interpolating between su(4 ) and @xmath448 with large af fluctuations but no static af order , and with no cdw or fm order .    in a strong magnetic field the zero - energy state of graphene has four - fold spin and valley degeneracy per landau level , and ( neglecting the lattice - scale interactions ) near the sample boundary one might expect the zero landau level to split into one positively - dispersing ( electron - like ) and one negatively dispersing ( hole - like ) mode for each spin projection . this would suggest a ground state having a bulk energy gap at charge neutrality but with electron - like and hole - like states of opposite spin polarization crossing at the edge of the sample ( with the edge - state structure being analogous to the quantum spin hall effect ) @xcite . however , experiments indicate that the ground state of charge - neutral graphene becomes _ strongly insulating _ at high magnetic fields @xcite . the detailed nature of this state remains uncertain , but it is generally expected to correspond to a spontaneously broken symmetry caused by the strong coulomb interactions among the electrons in the zero landau level . transport properties are not manifest in the algebraic solutions presented here but the coherent state approximation is equivalent to symmetry - constrained hartree  fock  bogoliubov ( hfb ) theory @xcite , suggesting that so(8 ) theory can be mapped onto hartree  fock ( hf ) transport calculations . hf calculations for armchair nanoribbons found that af and cdw states similar to ours have no edge currents @xcite . it may be speculated that our af states also are insulating and thus strong candidates for the high - field ground state , but confirmation requires more work . solutions depend on @xmath437 and @xmath489 in eq . ( [ modelhamiltonian ] ) , which define effective interactions in the truncated space [ highly renormalized relative to parameters in eq . ( [ so5_1.1a ] ) ] . they may be fixed by systematic comparison with data , enabling a robust prediction for the nature of the ground and other low - energy states . one expects modest impurity levels to modify the effective interaction parameters but to leave dynamical symmetries intact . a unified model of conventional superconductivity and high - temperature superconductors having an @xmath385 highest symmetry has been proposed @xcite . for sufficiently strong on - site coulomb repulsion ( true for the cuprates and approximately true for the fe - based compounds ) , the most general so(8 ) symmetry is reduced to its su(4 ) subgroup @xcite , which forms the basis of an su(4 ) dynamical symmetry model of high - temperature superconductivity . this su(4 ) model has three dynamical symmetry subchains ending in an @xmath583 subgroup representing conservation of spin and charge . physically , these su(4 ) subgroup chains represent    1 . an su(2 ) pseudospin subgroup chain that describes a @xmath584-wave singlet superconductor ( sc ) . an so(4 ) subgroup chain that describes an antiferromagnetic mott insulator ( af ) . an so(5 ) subgroup chain representing a critical dynamical symmetry that is soft with respect to af and sc fluctuations and interpolates between the af and sc collective modes . thus , in the su(4 ) model of high-@xmath585 superconductivity the so(5 ) subgroup chain plays a similar physical role as the so(7 ) subgroup chain of the graphene so(8 ) dynamical symmetry model and the so(7 ) subgroup chain of the nuclear so(8 ) model . in all three cases the subgroup chains represent the generalization of a quantum critical point to an entire quantum critical phase that exhibits large fluctuations ( in order - parameter space ) connecting collective modes defined by other dynamical symmetries of the problem .        in the high-@xmath585 case the modes connected by the @xmath586 critical dynamical symmetry are antiferromagnetism and superconductivity , in graphene the critical @xmath587 dynamical symmetry connects states with and without to nel ( af ) order , and in the nuclear structure case the critical @xmath587 dynamical symmetry connects collective states that differ in the relative contributions of nucleon pairs carrying total angular momentum 0 and total angular momentum 2 . thus physically these three applications of dynamical symmetry have little in common , but one sees that mathematically they have a deep similarity . although the microscopic physics differs fundamentally , at the level of the observed emergent collective modes in the system one sees that graphene quantum hall states , high - temperature superconductors , and broad classes of nuclear collective states have a unified description in terms of dynamical symmetries associated with very similar compact lie algebras . this remarkable similarity is illustrated in fig . [ fig : nuc_sc_graphene_chains ] the well - known quantum hall ferromagnetic su(4 ) symmetry of graphene in strong magnetic fields has been extended by adding to the particle  hole operators that generate su(4 ) a set of six creation and six annihilation operators that create or destroy fermion pairs in either a total valley isospin triplet , total spin singlet state , or a total valley isospin singlet , total spin triplet state ( the only possibilities allowed by the pauli principle ) . this extended set of operators is shown to close an so(8 ) algebra under commutation , which is formally analogous to the so(8 ) algebra of the fermion dynamical symmetry model of nuclear structure physics . this permits immediate adaptation of mathematical tools developed in nuclear physics to the graphene problem . the previously - known su(4 ) quantum hall ferromagnetism symmetry is recovered as one subgroup , but one finds a richer set of low - energy collective modes associated with the full subgroup structure of so(8 ) . by exploiting the established methodology of fermion dynamical symmetries , it was possible to decouple a collective - pair subspace from the full hilbert space of the problem , permitting exact , analytical , many - body solutions to be obtained in several physically - interesting limits . in addition to exact solutions in specific dynamical symmetry limits , a generalized so(8 ) coherent state approximation has been introduced that permits a broad range of solutions to be obtained even when not in the dynamical symmetry limits . the pairs spanning the collective subspace are shown to be analogous to pairs that have already been discussed at a qualitative level in the graphene literature as defining the possible broken symmetry ground states in the presence of strong electron  electron and electron  phonon correlations in the @xmath5 landau level . the development here places these pairs on a firm , unified mathematical footing and permits analytical solutions to be developed that explore the possible collective states that previously required numerical simulation for their quantitative description . finally , it has been shown that there are uncanny dynamical symmetry analogies among broken symmetry states for graphene in a strong magnetic field , high temperature superconductors , and strongly collective states in atomic nuclei . on the one hand this has the practical utility of allowing technology already developed in one field to be adapted easily to another . on the other hand , it implies a deep and intriguing mathematical affinity among physical problems that are not usually viewed as having more than a superficial connection .    _ acknowledgments._* * this work is supported by the basque government ( grant no .  it472 - 10 ) , the spanish micinn ( project no . fis2012 - 36673-c03 - 03 ) , and the basque country university ufi ( project 11/55 - 01 - 2013 ) . partial support was provided by lightcone interactive llc . the formalism discussed in this paper was introduced by postulating a set of physical operators that were shown explicitly to generate an so(8 ) lie algebra under commutation . it is illuminating to consider a somewhat different perspective on the motivation for introducing an so(8 ) symmetry for the graphene problem . it is well known that in a fermionic space having @xmath275 degrees of freedom the most general set of bilinear products @xmath588 of creation  annihilation operators and their hermitian conjugates generates an su(@xmath275 ) lie algebra under commutation @xcite . ( physically , the restriction to bilinear products corresponds to limiting consideration to two - body interactions . ) likewise , it is well known that adding to this particle  hole operator set the most general pair creation and annihilation operators @xmath589 and @xmath590 extends the su(@xmath275 ) algebra to so(@xmath591 ) . this extension is useful because sometimes more is less . the advantage of expanding the space from su(@xmath275 ) , with @xmath592 generators , to so(@xmath591 ) , with @xmath593 generators , is that the added pair operators permit the definition of a ( collective ) subspace of the full hilbert space spanned by the products of pair creation operators acting on the pair vacuum . if an effective hamiltonian is then constructed by writing the most general polynomial in the casimir invariants of all groups in the subgroup chains of the highest symmetry so(@xmath591 ) ( restricted to quadratic casimir invariants if one considers only two - body interactions ) , it will correspond to the most general hamiltonian that can be written in the collective subspace , and will be diagonal in the collective subspace basis for specific dynamical symmetry subgroup chains . thus , the schrdinger equation can be solved analytically in the symmetry limits defined by each dynamical symmetry subgroup chain , and even away from the symmetry limits it can be solved analytically in coherent - state approximation . the su(@xmath275 ) particle  hole algebra alone can be used to construct a hamiltonian that commutes with its generators , and the corresponding schrdinger equation can be solved analytically . but when one considers the realistic case of adding symmetry - breaking terms to the hamiltonian that do not commute with the generators of su(@xmath275 ) , the best that can be done analytically is to assume that these terms are small and that the physical solutions can be treated as small perturbations around the symmetric solution .    for the non - perturbative case where the added terms lead to spontaneously broken symmetry and new possible ground states that are not connected perturbatively to the symmetric ground state , one has no systematic way to construct the new ground state from the symmetric one except to guess it . but since the su(@xmath275)-symmetric solution is not connected analytically to non - perturbative broken symmetry solutions , the symmetries of the unperturbed ground state could be a poor guide to guessing the nature of the broken - symmetry states and one must rely on numerical solutions or other approximations and not symmetry to determine their properties . let us now apply this general discussion to quantum hall magnetism in graphene . as discussed in  [ h : symmetriesfqhegraphene ] , it is widely accepted that an approximate su(@xmath275 ) particle  hole symmetry with @xmath594 is relevant in the @xmath5 landau level because the dominant long - range coulomb interaction is su(4)-symmetric . however , the short - range terms that break this symmetry in the effective hamiltonian prevent the su(4 ) symmetry from providing a solution for the broken symmetry ground state , unless it is assumed that those terms only perturb the su(4)-symmetric solution ( small explicit symmetry breaking ) . but experimental evidence suggests that the true ground state of graphene in a strong magnetic field at low temperature breaks su(4 ) symmetry _ spontaneously , not explicitly _ ( see the introduction ) , and therefore is produced by a non - perturbative effect that can not be explained in terms of small fluctuations around the su(4)-symmetric solution . the possible ( spontaneously ) broken - symmetry states for the @xmath5 landau level have been described in terms of the most general sets of electron pairs occupying the two valley isospin and two spin degrees of freedom @xcite , @xmath595 { |0\rangle } , \label{basis1.0}\ ] ] where the vacuum state @xmath239 corresponds to completely filled landau levels for @xmath596 and completely empty landau levels for @xmath597 . each factor in the product @xmath598 creates a pair of electrons in the state @xmath599 at orbital @xmath600 of the @xmath5 ll , with @xmath601 equal to sublattice a or b , @xmath211 equal to spin up or down , and with the valley isospin and sublattice pseudospin identified : @xmath602 and @xmath603 . but the pair creation operators @xmath215 in eq . ( [ basis1.0 ] ) and their hermitian conjugates are not generators of su(4 ) and do not commute with the su(4)-symmetric hamiltonian , and the collective states of the form ( [ basis1.0 ] ) that are of interest in the present context are unlikely to represent small fluctuations around the su(4)-symmetric solution . thus the nature of these collective states is not determined by the su(4 ) symmetry and had to be investigated by numerical calculations using a small basis in prior work @xcite .    on the other hand , the so(8 ) pair generators introduced in eq . ( [ coupled1.4 ] ) are included in the most general collective pairs generated by the @xmath604 operators in eq . ( [ basis1.0 ] ) , and include the collective degrees of freedom discussed in ref . @xcite ( see  [ sh : otherstates ] and fig . [ fig : s_d_pairs_brillouin ] ) . thus , the so(8 ) fermion dynamical symmetry permits the nature of possible spontaneously broken symmetries to be investigated in terms of _ symmetry properties _ that permit analytical solutions for the broken - symmetry states . hence , the formalism described in this paper may also be introduced by the following logic . the su(4 ) particle  hole symmetry generated by the operators @xmath125 defined in eq . ( [ algebra1.2 ] ) is known to provide a good starting point for graphene quantum hall states dominated by the long - range coulomb interaction , but does not describe quantitatively the broken - symmetry modes discussed by kharitonov @xcite in terms of collective pairs resulting from short - range correlations . motivated by the preceding discussion in this appendix , the su(4 ) generator set may be extended to include the possible pair creation and pair annihilation operators operating in the space corresponding to the indices in the su(4 ) basis . by the general su(@xmath275 ) @xmath605 so(@xmath591 ) extension discussed above , this gives the so(8 ) lie algebra of eq . ( [ algebra1.3 ] ) . hence , expanding the algebra from su(4 ) to so(8 ) introduces the capability to explore _ analytically _ the possible collective states following from perturbation of the su(4 ) quantum hall ferromagnet by short - range interactions that break su(4 ) both explicitly and spontaneously . this appendix collects some useful transformations among the several bases that have been employed in this paper . for brevity , in the following @xmath606 will be termed the nuclear so(8 ) basis and @xmath607 will be termed the graphene so(8 ) basis .    in transforming from the nuclear so(8 ) basis to the graphene so(8 ) basis the particle number ( charge ) operator @xmath17 or @xmath266 and the 12 pairing operators @xmath608 are retained , but the 15 su(4 ) generators @xmath609 in the nuclear representation are replaced with the 15 su(4 ) generators @xmath610 defined in the graphene representation of eq . ( [ algebra1.4 ] ) . the explicit transformation from the @xmath611 generators to the @xmath612 generators is given by    @xmath613 \label{ngx10 } \\ { \pi_{yx } } & = \tfrac{i}{2}\left [ \sqrt{\tfrac{2}{5 } } p^3_{-3 } + p^3_3 \right . \nonumber \\ & \left . + \sqrt{\tfrac{2}{5}}\left(p^1_{-1 } +   p^1_1\right )   - \sqrt{\tfrac{3}{5}}\left(p^3_1 + p^3_{-1}\right)\right ] \label{ngx11 } \\ { \pi_{zx } } & = -\tfrac{1}{\sqrt2 } \left ( p^3_2 + p^3_{-2}\right ) \label{ngx12 } \\ { \pi_{xy } } & = \tfrac{i}{2}\left[\sqrt{\tfrac{2}{5 } } p^3_{-3 } + p^3_3 \right . \nonumber \\ & \left.- \sqrt{\tfrac{2}{5}}\left(p^1_{-1 } +   p^1_1\right ) + \sqrt{\tfrac{3}{5}}\left(p^3_1 + p^3_{-1}\right)\right ] \label{ngx13 } \\ { \pi_{yy } } & = \tfrac{1}{2}\left [ \vphantom{\sqrt{\tfrac{2}{5 } } } -p^3_{-3 } + p^3_3 \right . \nonumber \\ & \left.-   \sqrt{\tfrac{2}{5}}\left(p^1_1 - p^1_{-1}\right ) +   \sqrt{\tfrac{3}{5}}\left(p^3_1 - p^3_{-1}\right)\right ] \label{ngx14 } \\ { \pi_{zy } } & = -\tfrac{i}{\sqrt2 } \left ( p^3_2 - p^3_{-2}\right ) \label{ngx15}\end{aligned}\ ] ]      in eqs . ( [ algebra1.8])([algebra1.11 ] ) the graphene basis @xmath137 has been expressed in terms of the generators @xmath125 defined in eq . ( [ algebra1.2 ] ) . the inverse transformations giving the @xmath125 in terms of the @xmath614 are        where the unlisted operators may be obtained from @xmath616 and the diagonal operators have been assumed to obey the u(4 ) constraint @xmath617 with @xmath618 the total particle number and @xmath128 the total pair degeneracy given by eq . ( [ degen1.2 ] ) . since from eq . ( [ orderp1.1 ] ) the order parameters for the quantum hall ground states are functions of the expectation values for the number operators @xmath619 specifying the population of the four basis states in table  [ tb : quantumnumbermapping ] and fig .  [ fig : graphenebasis ] labeled by the index @xmath100 , it is useful to have explicit expressions for them in terms of the @xmath260 operators . these are    @xmath620 n_2 & = \tfrac14 n + \tfrac{\sqrt5}{10}p^1_0 - \tfrac12 p^2_0 - \tfrac{3\sqrt5}{10}p^3_0 \label{ptoi2 } \\[2pt ] n_3 & = \tfrac14 n - \tfrac{\sqrt5}{10}p^1_0 -\tfrac12 p^2_0 + \tfrac{3\sqrt5}{10}p^3_0 \label{ptoi3 } \\[2pt ] n_4 & = \tfrac14 n - \tfrac{3\sqrt5}{10}p^1_0 + \tfrac12 p^2_0 -\tfrac{\sqrt5}{10}p^3_0\end{aligned}\ ] ]                                                                                f. t. arecchi , e courtens , r. gilmore , and h. thomas , phys . * a6 * , 2211 ( 1972 ) . r. gilmore , ann . phys .  * 74 * , 391 ( 1972 ) . a. m. perelomov , commun . * 26 * , 222 ( 1972 ) . r. gilmore , rev . de fisica * 23 * , 143 ( 1974 ) . zhang , d. h. feng , and r. gilmore , _ rev . mod.phys._ * 62 * , 867 ( 1990 ) . a. f. young , c. r. dean , l. wang , h. ren , p. cadden - zimansky , k. watanabe , t. taniguchi , j. hone , k. l. shepard and p. kim nature physics * 8 * , 550556 ( 2012 ) .
Abstract:

Answer: a formalism is presented for treating strongly - correlated graphene quantum hall states in terms of an so(8 ) fermion dynamical symmetry that includes pairing as well as particle hole generators . the graphene so(8 ) algebra is isomorphic to an so(8 ) algebra that has found broad application in nuclear physics , albeit with physically very different generators , and exhibits a strong formal similarity to su(4 ) symmetries that have been proposed to describe high - temperature superconductors . the well - known su(4 ) symmetry of quantum hall ferromagnetism for single - layer graphene is recovered as one subgroup of so(8 ) , but the dynamical symmetry structure associated with the full set of so(8 ) subgroup chains extends quantum hall ferromagnetism and allows analytical many - body solutions for a rich set of collective states exhibiting spontaneously - broken symmetry that may be important for the low - energy physics of graphene in strong magnetic fields . the so(8 ) symmetry permits a natural definition of generalized coherent states that correspond to symmetry - constrained hartree  fock  bogoliubov solutions , or equivalently a microscopically - derived ginzburg landau formalism , exhibiting the interplay between competing spontaneously broken symmetries in determining the ground state .