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: bistable systems are widely utilized to model the biological processes which exhibit distinct phenotypes under homogeneous conditions @xcite . switching between phenotypes ( stable states ) is facilitated by the stochasticity arising from molecular noise @xcite . the paths of switching have been studied in various systems , such as the @xmath0-phage lysis - lysogeny decision @xcite and cellular development and differentiation @xcite , to gain insights into the molecular processes of biological decision making .    in nature , we sometimes observe situations where multiple bistable systems are coupled . for example in bacterial quorum sensing , every cell produces small signaling molecules whose production is regulated by a positive feedback to synchronize the population . this positive feedback may induce bistability between the high and low concentrations of the signaling molecules . the cells are further coupled by secreting and sensing the signaling molecules in the medium@xcite . another example is the toxin - antitoxin ( ta ) loci in _ escherichia coli _ , where there are 10 known mrnase toxins and every pair may act as a bistable system allowing the cells to switch between the normal growing state and the dormant state which exhibits antibiotic persistence @xcite . the ta systems may interact each other by interfering the protein synthesis and the cellular growth . motivated by these systems , in this paper we analyze a coupled bistable system that consists of interacting identical bistable motifs . we consider the case where coupling is such that the coupled system itself is also bistable and for each stable state the individual motifs are in the same steady state . in other word , the coupling is positive to allow all the bistable motifs to jointly switch from one state to the other state . while the individual bistable systems without coupling show the same pattern of transition paths by definition , the switching properties of the coupled system remain unclear . the individual systems may transit synchronously at the same pace , resulting in one transition path . however , it is also possible that individual systems lead the switching process . as a result , the transition paths of the coupled system are split into multiple ones , a phenomenon called `` bifurcation of transition path '' . since many properties of the coupled system , such as the transition rates between the steady states , are dependent on the transition paths , it is interesting to study whether the transition paths of a given coupled bistable system bifurcate and how this bifurcation relates to the individual bistable systems and their coupling . bifurcation of transition paths was first demonstrated in the maier - stein model @xcite , but its relevance and implication to biochemical systems remain to be explored . in this work , we first construct a model of coupled bistable gene circuits in section [ modeling ] . then we demonstrate that the transition paths bifurcate with appropriate parameter sets . in section [ theory ] , we construct a general formulation of coupled bistable systems . we consider the transition of this general model between its steady states as a noise - induced exit process from a metastable state and propose a criterion for the bifurcation of transition paths by extending the previous works on the maier - stein model @xcite . finally in section [ application ] , we apply our criterion to the model of coupled bistable gene circuits . we confirm the theory numerically and discuss the transition rates . consider a model of coupled bistable gene circuits ( fig . [ model]a ) . first we restrict our attention to one gene . the promoter of the gene is weak and the proteins of this gene bind to the promoter in the form of tetramers to activate the gene expression . we may model the proteins of this gene using @xmath1 where @xmath2 is the concentration of the protein , @xmath3 refers to the basal synthesis rate of the protein , the hill term describes the activation of gene expression by the tetramers and the last term models the linear degradation . with appropriate parameter values , the positive feedback on gene expression allows eq . [ example1d ] to show bistability . we now couple @xmath4 such genes in one cell and we assume that all these genes ( and their promoters and proteins ) have the same kinetic properties . the genes are coupled in a way that the proteins of the genes are well - mixed in the cells and the mixture activates a cell by binding to a promoter in a tetramer . multiple coupling strategies may be used . for example , if the proteins of the @xmath4 genes are identical , any four monomers may bind a promoter and we may model the coupled bistable gene circuits using @xmath5 where @xmath6 is the concentration of the products of the @xmath7-th gene and @xmath8 is the hill constant for the coupled system . meanwhile , if the genes are equipped with identical promoters but encode different proteins , and the tetramer activating the gene expression consists of four monomers from the same gene , we model the coupled system using @xmath9    we then generalize these examples and we propose the following model of coupled bistable gene circuits @xmath10 where the parameter @xmath11 , called the configuration parameter , governs the general coupling strategy . it is straightforward that @xmath12 corresponds to eq . [ alpha1 ] and @xmath13 corresponds to eq . [ alpha4 ] . here we allow @xmath11 to be arbitrary positive values , though not all coupling strategies are biologically plausible . we also choose the values of the hill constants ( @xmath8 in eq . [ alpha1 ] and @xmath14 in eq . [ alpha4 ] ) such that the steady states of every individual gene are not affected by the coupling . if the model for one gene . [ example1d ] ) is bistable , it is straightforward that the general model of coupled bistable gene circuits ( eq . [ example ] ) is also bistable and contains three steady states : two stable ones and one saddle . furthermore , one can show that @xmath15 holds at every steady state . for convenience , throughout this work we call the stable state where the concentrations of all gene products are low the `` lower stable steady state '' @xmath16 and the other stable steady state the `` higher stable steady state '' @xmath17 .    if we set the volume of the cell to be @xmath18 , the concentrations of proteins @xmath6 can be converted to the absolute numbers of molecules ( @xmath19 , should be an integer ) and the coupled system is governed by the chemical master equation . one may then sample the transition paths between the two steady states using the gillespie algorithm @xcite . here we restrict our attention to @xmath20 and the transition from the lower stable steady state to the higher one . as illustrated in fig . [ model]b , when @xmath12 , the transition paths are narrowly distributed around the diagonal ( the line satisfying @xmath21 ) , suggesting that the switching of the two genes is synchronized . the distribution becomes wider as the value of @xmath11 increases , and at @xmath22 , the transition paths exhibit a visible bifurcation . therefore , coupled bistable systems are capable of exhibiting both bifurcating and non - bifurcating transition paths , and we can modulate the bifurcation pattern for the model of coupled bistable gene circuits with the parameter @xmath11 . ( color online ) model of coupled bistable gene circuits . * a*. a schematic illustration of the model . a cell contains @xmath4 genes with weak promoters and their products mix in the cell and activate the gene expression by binding to the promoters in the form of tetramers . * b*. distribution of transition paths from the lower stable steady state ( @xmath16 ) to the higher one ( @xmath17 ) . we set @xmath23 and we carry out 100 simulations using gillespie algorithm . the distribution of the last instanton trajectories , i.e. the trajectories associated with the successful escapes , is calculated . we present the frequencies in the logarithmic scale in the form of heat plots . red indicates high frequency and blue indicates low frequency . the black lines represent the most probable escape paths computed in the zero - noise limit . the parameter values are @xmath24 , @xmath25 , @xmath26 and @xmath27 . in appendix d , we verify that the gillespie simulation is carried out in the low - noise limit . ] in this section we develop a criterion for coupled bistable systems to predict whether the transition paths bifurcate or not . to formulate a general model for coupled bistable systems , we notice that in the model of coupled bistable gene circuits ( eq . [ example ] ) , every gene is governed by the concentration of its own protein ( @xmath6 , in the degradation term ) and the average concentration of all proteins ( @xmath28 , in the production term ) . here for a general coupled bistable system , we may model the deterministic drifts of the individual bistable systems in the same fashion and describe the effect of noises using @xmath29 which is interpreted as an ito - type stochastic differential equation . here @xmath6 is the state of the @xmath7-th bistable system and @xmath4 is the number of systems to be coupled . the function @xmath30 , defined as @xmath31 , computes the average state . we choose this formulation because it allows modulating of @xmath11 and @xmath4 without changing the steady states . obviously , the state of every individual bistable system is governed by its own state and the average state of all bistable systems , as modeled by the deterministic drift @xmath32 . the second term of eq . [ model_defn ] arises from expanding the chemical master equation in the continuous limit and serves as the noise term for the coupled bistable systems . here @xmath33 are independent gaussian white noise sources ( @xmath34 ) and we assume the multiplicative noise @xmath35 for the @xmath7-th system is also a function of its own state and the average state . the overall noise is modulated by a small parameter @xmath36 in order to keep the validity of the continuous limit . we assume that the deterministic drift with @xmath37 gives three steady states - two stable ones and an unstable one . we further assume the coupled system also allows three steady states - two stable ones and a saddle point . at each steady state , the states of all bistable systems are assumed to be equal . a wide range of coupled bistable systems may be modeled in the fashion of eq . [ model_defn ] . for example , the model of coupled bistable system , as mentioned in section [ modeling ] , corresponds to @xmath38 and a noise amplitude of the form @xmath39 @xcite . the parameter @xmath11 in @xmath40 corresponds to the coupling mode of the genes and the small @xmath41 corresponds to a cell with a large volume @xmath18 ( in the large volume limit , @xmath42 ) . we derive our criterion by performing linear perturbation analysis in the zero - noise limit . when the noise level of a system reaches zero ( @xmath43 ) , the transition between stable steady states is dominated by the path(s ) associated with the lowest energy cost ( or `` action '' in standard terminology ) . this path is called the most probable escape path ( mpep)@xcite . the mpep can be quantified analytically since it is governed by the freidlin - wentzell hamiltonian @xcite , allowing ones to use the tools of analytical mechanics for derivation . one property of mpep in conventional cases is that the mpep of transition passes through the saddle of the system ( refer to @xcite ) and after that the system follows the deterministic flow to the final stable steady state @xcite . in our formulation of coupled bistable systems eq . [ model_defn ] , the individual bistable systems are always driven by the same deterministic drift after passing the saddle . we then conclude that no bifurcation of transition paths will occur between saddle and the final steady state . therefore , the transition capable of bifurcation is the one from the initial stable steady state to the saddle , and we may restrict our attention to this exit process . the exit process may be studied with the approach of linear perturbation . to be precise , we first assume that the mpep does not bifurcate , i.e. all the bistable systems take the same state during the transition . we then perturb the shape of mpep by allowing some bistable systems to be in different states from the others during transition and we examine how the associated actions for the perturbed mpep change . if the actions always increase regardless of the perturbation , the non - bifurcating path is locally energetically favorable . for the model of coupled bistable gene circuits , we confirm by numerical simulation that this path is indeed the mpep . meanwhile , if any perturbation leads to a decreased action , the non - bifurcating path is energetically unfavorable . the mpep should then be some other paths . due to the symmetry of the coupled system , multiple paths must exist . therefore , the coupled system exhibits a bifurcation of transition paths . we now follow this approach and we begin our analysis by deriving the equation governing the mpep , closely following the procedure found in previous work@xcite . the model of the coupled bistable system eq . [ model_defn ] , viewed as a stochastic differential equation , gives the fokker - planck equation @xcite @xmath44\ ] ] where @xmath45 denotes the probability of the coupled system at state @xmath46 and time @xmath47 . the drift vector @xmath48'$ ] satisfies @xmath49 and the covariance @xmath50 is a diagonal matrix satisfying @xmath51 . as the initial condition , we use a delta function at the initial stable steady state . the separatrix of the basins gives an absorbing boundary condition . in the zero - noise limit , @xmath52 and we replace @xmath45 by the quasi - steady state distribution @xmath53 governed by @xcite @xmath54 = 0\ ] ] we further assume the quasi - steady state distribution to take an arrhenius form @xmath55 , where @xmath56 is the quasi - potential@xcite . we perform a wkb expansion by plugging the form into eq . [ fp_quasi ] and keeping the terms of the lowest order of @xmath41 . one can show that the mpep of the transition is a classical zero - energy trajectory of the freidlin - wentzell hamiltonian @xmath57 where the momentum vector can be computed by @xmath58@xcite . the equation @xmath59 can be viewed as the equation governing the mpep , but it is worth noting that not all trajectories satisfying @xmath59 are the mpeps , as one has to examine the associated actions . we now determine whether the mpep of the transition bifurcates or not by linearly perturbing the non - bifurcating path and examining the actions . we can show that the non - bifurcating path is a zero - energy trajectory of the freidlin - wentzell hamiltonian ( appendix a ) . the action associated with the non - bifurcating path @xmath60^t$ ] can then be computed by the quasi - potential @xmath56 and one may quantify how the actions change with perturbations to this path by considering the expansion @xmath61 . due to symmetry , the first - order term vanishes if the system is perturbed in the directions perpendicular to the non - bifurcating path . note that these directions are of our primary interest , and we then discuss the second - order term and examine the eigenvalues of the hessian matrix @xmath62 . the equation governing the hessian matrix , computed by differentiating the @xmath59 twice over @xmath46 and making use of the hamiltonian equation @xmath63 ( also refer to ref . @xcite ) , is @xmath64 where @xmath65 is a linearization of the drift vector . the linearization of the variances gives @xmath66 . we may simplify eq . [ hessian_full ] significantly ( details in appendix b ) . by utilizing the fact that all the bistable systems share the same state along the non - bifurcating path , we may express the hessian matrix @xmath67 , where @xmath68 is an @xmath69 identity matrix and @xmath70 is an @xmath69 one matrix . we can further express the remaining terms , namely @xmath71 , @xmath72 , @xmath73 and @xmath74 , in a similar fashion . by plugging the new forms into eq . [ hessian_full ] , we obtain the equations governing @xmath75 and @xmath76 : @xmath75 is shown to follow @xmath77 and @xmath76 is shown to follow eq . [ z2_dt ] . here @xmath78 , @xmath79 and @xmath80 and similar for @xmath81 . in eq . [ z1_dt ] , @xmath2 and @xmath82 refer to the state and the momentum of every bistable system . the two arguments to the function @xmath32 and @xmath81 are @xmath2 .    here we analyze the eigenvalues of the hessian matrix @xmath83 , namely @xmath75 ( repeat @xmath84 times ) and @xmath85 , and we need to determine the eigenvalues corresponding to perturbations that break the non - bifurcating assumption . in the appendix c we show that the eigenvalue @xmath85 is only governed by one individual bistable system and it represents the change in action along the non - bifurcating path . so this eigenvalue is not relevant in our context and we focus on the directions perpendicular to the non - bifurcating path , which is given by the eigenvalue @xmath75 . we should examine the sign of @xmath75 along the non - bifurcating path from the initial stable steady state to the saddle . if @xmath75 is always positive , the perturbation is energetically unfavorable and the non - perturbing path is locally stable . if @xmath75 is negative somewhere , the mpep bifurcates into multiple paths . sometimes it is convenient to express @xmath75 as a function of @xmath2 rather than time . based on the freidlin - wentzell hamiltonian for the coupled bistable system , we have that the momentum vector may be computed by @xmath86 . in addition , the transition for the individual bistable system satisfies that @xmath87 since it is moving against the deterministic drift @xcite . by plugging these relations into eq . [ z1_dt ] , we show that @xmath88 usually we assume that the system is equipped with a constant , non - zero noise at the initial stable steady states @xcite . we may then solve the equation @xmath89 at the initial stable steady states and we set the nontrivial solution to be the initial condition , which is @xmath90 . we can integrate eq . [ z1 ] from the initial stable steady state to the saddle and determine whether the mpep bifurcates or not by looking at the sign of @xmath75 . [ z1 ] generates several insights . the first two terms of eq . [ z1 ] illustrate that the nonlinearity induces the bifurcation of mpep : @xmath91 and @xmath92 gives information about the nonlinearity for a bistable system without coupling , while @xmath93 gives the nonlinearity created by the coupling of the bistable system . we also notice that eq . [ z1 ] contains no terms of @xmath4 , indicating that the bifurcation of mpep is not dependent on the number of bistable systems to be coupled . in this section we apply our theoretical criterion to the model of coupled bistable gene circuits ( eq . [ fuv]-[guv ] ) , and we first consider the mpep of the model with two coupled genes . for the forward transition from the lower stable steady state @xmath16 to the higher one @xmath17 , if we choose @xmath12 , i.e. no nonlinearity is generated from the coupling of the genes , and integrate eq . [ z1 ] from @xmath16 to the saddle , we find that the eigenvalue @xmath75 monotonously decreases along the path ( fig . [ 2d_result]a , left panel ) . @xmath75 remains positive throughout the path , indicating that at @xmath12 the mpep does not bifurcate . for @xmath94 , the nonlinearity from the coupling of the genes drives @xmath75 to a lower value , but @xmath75 remains positive and the mpep does not bifurcate as well . for @xmath95 and 4 , the coupling is highly nonlinear and @xmath75 reaches negative values and diverges , suggesting that the mpep bifurcates . meanwhile , for the backward transition from the higher stable steady state @xmath17 to the lower one @xmath16 , we may integrate eq . [ z1 ] from @xmath17 to the saddle and we find that the eigenvalue @xmath75 remains positive regardless of the values of @xmath11 ( fig . [ 2d_result]b , left panel ) , suggesting that the mpep never bifurcates .    to verify our prediction , we use the geometric minimum action method ( gmam ) @xcite to compute the mpep numerically ( fig . [ 2d_result]ab , right panels)@xcite . as predicted , the mpep for the forward transition with @xmath12 and @xmath96 do not bifurcate , the one with @xmath95 and @xmath97 bifurcate to two symmetric ones with equal actions and no bifurcation can be observed for the backward transition . therefore , the predictive power of the theoretical criterion is confirmed and this bifurcation pattern can explain the distribution of transition paths illustrated in fig . [ model]b . we then examine the transition rates of two coupled genes . the action associated with the mpep ( @xmath98 ) is defined as the difference in the quasi - potential between the initial and final stable steady states , and one can show that the transition rate is proportional to @xmath99 @xcite . we sample the configuration @xmath11 in a wide range and compute the actions . in accordance with the maier - stein model @xcite , the actions for the model of coupled bistable gene circuits are not modulated by the value of @xmath11 if no bifurcation of mpep exists ( fig . [ nd_result]a , @xmath20 ) . meanwhile , in the presence of bifurcation , the actions are driven to a lower value by this configuration parameter and the transition rates increase . in short , the bifurcation of transition paths facilitates transition . ( color online ) analysis of the coupling of two genes . the eigenvalues induced by the coupling of genes were integrated along the path with no bifurcation from the initial stable steady state to the saddle and plotted in the left panels . the right panels illustrate the most probable escape paths which are numerically computed with gmam method . the dashed lines represent the separatrix of the basins surrounding the two stable steady states . the arrows represent the directions of transition . * a*. forward transition from the lower stable steady state ( @xmath16 ) to the higher one ( @xmath17 ) . * b*. backward transition from the higher stable steady state ( @xmath17 ) to the lower one ( @xmath16 ) . ] finally we examine the coupling of more genes . as predicted by the theory , the bifurcation in the transition paths is not dependent on the number of genes and the actions @xmath98 exhibit a similar dependency on @xmath11 ( fig . [ nd_result]a , upper panels ) . numerical simulation reveals that bifurcated mpep has an interesting shape . as illustrated in fig . [ nd_result]b for @xmath100 and @xmath13 , we plot how the states of individual genes differ from the arithmetic average of all states , and it shows that only one gene is in the leading position ( blue curve ) and the remaining ones follow with the same states ( other curves ) . we numerically verify that the pattern of one leading / the rest following is the optimal one and other patterns are associated with a higher action . this finding is consistent with the probabilistic view that transient alternation in the state of one gene against the drift is more probable than alternating multiple simultaneously . numerical analysis also shows that the actions associated with non - bifurcated mpep scale linearly with the number of genes @xmath4 , suggesting an exponential scaling law in the transition rate . the actions for bifurcated mpep scales sublinearly with @xmath4 ( fig . [ nd_result]c ) . this pinpoints that the bifurcation in mpep leads to a softer dependency between the transition rate and the size of the system . ( color online ) analysis of the coupling of multiple genes . * a*. the actions @xmath98 ( defined as the difference in the quasi - potential between the initial stable steady state and the saddle ) are computed numerically for both forward transitions ( left panels ) and backward transitions ( right panels ) . we consider the coupling of two to five genes and sample the configuration @xmath11 in a wide range ( the sampled values are shown in dots ) . the dash line represents the predicted critical @xmath11 for bifurcation . the theory predicts that the transition paths for systems with @xmath101 are bifurcated and the action becomes smaller . with the present numerical accuracy , we are able to show that the transition paths are bifurcated , but not able to find the correct action . * b*. the most probable escape path of the forward transition for systems with @xmath100 and @xmath13 . the horizontal axis represents the arithmetic average of the genes states ( @xmath102 ) and the vertical axis represents the difference between the state of individual genes and the average state ( @xmath103 ) . * c*. the actions @xmath98 of forward transitions for multi - dimensional systems with @xmath12 ( no bifurcation ) and 4 ( bifurcation ) . ] in this work , taking the model of coupled bistable gene circuits as an example , we consider a general model of coupled bistable systems and derive a criterion to determine whether the most probable escape paths of the transition bifurcate or not . we show that in the present setup where the steady states are not modulated by @xmath11 and @xmath4 , this criterion is independent on the number of individual bistable systems . furthermore , we apply our criterion to the model of coupled bistable gene circuits and verify the theory s predictive power numerically . numerical analysis reveals that only one bistable system takes the leading position in the bifurcated mpep . we also show that the transition rates associated with non - bifurcated mpep scale exponentially with the number of bistable systems while the coupled systems with bifurcated mpep exhibit a softer scaling behavior . the theoretical criterion in this work is developed for a restricted class of systems satisfying the conditions that every individual system is bistable and its steady states are not affected by the coupling and the mpep passes through the saddle of the coupled system . one may study bifurcation of a wide range of systems by following our procedures and making appropriate changes . for example , it is possible to study the coupling of identical systems where every individual system is not bistable but the coupled system is , as well as to study the coupled systems whose the steady states are dependent on the coupling . the saddle point could also be replaced by the `` global maximum along the dominant path '' defined by the point along the mpep where the drift along the path changes its sign ( refer to @xcite ) . finally , one may extend the theory to study the coupling of multidimensional bistable systems , though one needs to numerically search for the optimal transition path for one bistable system and integrate @xmath75 ( now a matrix rather than a scalar ) along this path . in summary , our work provides with a convenient method to study the bifurcation of transition paths of coupled bistable systems and may lead to wide range of practical insights . c.t . and n.m . acknowledge kim sneppen and erik aurell for helpful discussions . this study is funded by the danish national research foundation through the center for models of life ( c.t . and n.m . ) and center for bacterial stress response and persistence ( n.m . here we show that the non - bifurcating path of a coupled bistable system is a zero - energy trajectory of the corresponding freidlin - wentzell hamiltonian . first we restrict our attention to one bistable system . since it is one - dimensional , there exists only one transition path in the zero - noise limit and it is obviously a zero - energy path of the corresponding freidlin - wentzell hamiltonian @xmath104 where we write the drift and variance explicitly . we then discuss the non - bifurcating path of a coupled bistable system . note that at every point along the non - bifurcating path , @xmath105 and @xmath106 hold , and we have that the drifts for every individual bistable systems are the same and the diagonal elements of the covariance matrix @xmath50 are equal as well . one may then express the corresponding freidlin - wentzell hamiltonian in the form that @xmath107 where we use @xmath2 and @xmath82 to represent the state and momentum of every bistable system . in other words , the coupling vanishes effectively and the hamiltonian is only dependent on the transition of every individual bistable system . since the transition path of every individual system is zero - energy , we can claim that the non - bifurcating path of the coupled system should satisfy @xmath59 . we compute the eigenvalues of the hessian matrix @xmath108 along the non - bifurcating path . the element @xmath109 in the matrix @xmath83 contains information about the interaction between the @xmath7-th and the @xmath110-th bistable systems . recall that all the bistable systems take the same state along the non - bifurcating path , and we have that all the diagonal elements should take the same value ( denoted as @xmath111 ) and same for the off - diagonal elements ( denoted as @xmath76 ) , though we should notice that the diagonal elements may not equal to the off - diagonal ones . the hessian matrix @xmath83 then takes the form of @xmath112 and it is straightforward to show that the eigenvalues are @xmath75 ( repeat @xmath84 times ) and @xmath85 . the non - bifurcating path also allows us to simplify the computation of the drifts and the noise levels in eq . [ hessian_full ] . take the drifts as an example . since all the bistable systems take the same state , for any bistable systems @xmath7 , @xmath110 , @xmath113 and @xmath114 ( @xmath115 , @xmath116 ) , we have @xmath117 and therefore @xmath118 . the linearization of the drifts is then @xmath119 . note that the function @xmath30 is symmetric for all individual bistable systems , and we have @xmath120 . we then show that @xmath121 . in other words , the off - diagonal elements of the linearization matrix @xmath71 take the same value . similarly , one can prove that the diagonal elements also take the same value and we can express @xmath71 in the form of @xmath122 where the arguments to the functions are @xmath2 and @xmath123 is defined as @xmath124 . meanwhile , we may also simplify the hessian of the drift @xmath73 . by noticing that @xmath106 , we have that @xmath125 where @xmath82 is the momentum of one bistable system , @xmath126 is defined as @xmath127 and @xmath128 is defined as @xmath129 for @xmath115 . following the same procedure , one can show that @xmath130 and @xmath131    we are readily to obtain the eigenvalues of the hessian matrix @xmath83 by plugging the equations above into eq . [ hessian_full ] . by making use the facts that @xmath132 , @xmath133 and @xmath134 along the non - bifurcating path , we have eq . [ z1_dt ] and        we analyze the eigenvalue @xmath85 along the non - bifurcating path . by making use of eq . [ z1_dt ] and eq . [ z2_dt ] , we show that this eigenvalue follows @xmath136 in eq . [ z1nz2_dt ] , the term @xmath137 is equivalent to @xmath138 , @xmath139 is equivalent to @xmath140 and similar for the terms of @xmath32 . since @xmath141 and @xmath142 are the drift and noise induced by one bistable system , we claim that the eigenvalue @xmath85 corresponds to the direction of perturbation where the coupling of multiple bistable systems effectively vanishes . the only possible direction is the one along the non - bifurcating path . therefore , the eigenvalue @xmath85 should not be considered in our context as it does not break the non - bifurcating assumption . we should then focus on the eigenvalue @xmath75 , as discussed in section [ theory ] . we verify that the gillespie simulation in fig . [ model]b is carried out in the low - noise limit . we perform simulation with three volumes ( @xmath143 ) and estimate the actions by fitting the simulation data to @xmath144 where @xmath145 is the mean transition time and @xmath146 . we show that the estimated actions are in good agreement with the values associated with the mpep which are computed in the low - noise limit ( fig . [ nd_result ] , table [ comparison_dw ] ) . furthermore , we plot the mpep on top of the heat plots in fig . [ model]b and we show that the transition paths sampled from gillespie simulation match the mpep . therefore , we claim that our gillespie simulation is in the low - noise limit .       5 j.e . ferrell , w. xiong , chaos . * 11 * , 227 ( 2001 ) . ferrell , curr . opin . cell . * 14 * , 140 ( 2002 ) . veening , w.k . smits , o.p . kuipers , annu . micro . * 62 * , 193 ( 2008 ) . g. balazsi , a. van oudenaarden , j.j . collins , cell . * 144 * , 910 ( 2011 ) e. aurell , k. sneppen , phys . rev . lett . * 88 * , 048101 ( 2002 ) j. wang , k. zhang , l. xu , e. wang , proc . 108 * , 8257 ( 2011 ) m.b . miller , b.l . bassler , annu . microbiol . * 55 * , 165 ( 2001 ) j. muller , c. kuttler , b.a . hense , m. rothballer , a. hartmann , j. math . biol . * 53 * , 672 ( 2006 ) e. maisonneuve , k. gerdes , cell * 157 * , 539 ( 2014 ) i. cataudella , a. trusina , k. sneppen , k. gerdes , n. mitarai , nuc . acids res . * 40 * , 6424 ( 2012 ) . fasani , m.a . savageau , proc . 110 * , e2528 ( 2013 ) . maier , d.l . stein , j. stat . phys . * 83 * , 291 ( 1996 ) d.t . gillespie , j. phys . chem . * 81 * , 2340 ( 1977 ) . if we expand the chemical master equation for the model of coupled bistable gene circuits using kramers - moyal expansion and preserve up to the second order , we obtain a fokker - planck equation ( eq . [ fokker - planck ] ) with a drift vector equal to the production rates minus the degradation rates and a diagonal diffusion matrix where the diagonal elements equal to the sum of the production and degradation rates . by defining the langevin equation using eq . [ fuv ] and [ guv ] , we will obtain the same fokker - planck equation and thus the definition is physical . alternatively , one may derive chemical langevin equation from the master equation and it contains the sum of two noise terms : one for production @xmath147 and one for degradation @xmath148 @xcite . within the scope of this manuscript , these two formulations are equivalent as they give the same fokker - planck equation . maier , d.l . stein , phys . e. * 48 * , 931 ( 1993 ) m.i . freidlin , a.d . wentzell , random perturbations of dynamical systems , 3rd ed . ( springer - verlag , berlin , 2012 ) r.s . maier , d.l . stein , siam j. appl . math . * 57 * , 752 ( 1997 ) h. risken , the fokker - planck equation : methods of solution and applications , 2nd ed . ( springer - verlag , berlin , 1996 ) t. naeh , m.m . klosek , b.j . matkowsky , z. schuss , siam j. appl . 50 * , 595 ( 1990 ) r. kupferman , m. kaiser , z. schuss , e. ben - jacob , phys . a. * 45 * , 745 ( 1992 ) eq . [ fw_hamiltonian ] corresponds to the system constructed by preserving the kramers - moyal expansion of the chemical master equation up to the second - order term ( cf . m. heymann , e. vanden - eijnden , phys . * 100 * , 140601 ( 2008 ) . m. heymann , e. vanden - eijnden , commun . pure appl . math . * 61 * , 1052 ( 2008 ) . the numerical simulation is carried out with eq . [ fw_hamiltonian ] which is a quadratic approximation to @xmath149 $ ] where @xmath150 is the conjugate momentum for @xmath6 @xcite . to verify the validity of the quadratic approximation , we perform numerical simulation with @xmath151 . with the present numerical accuracy , we do not observe any difference in the mpep and the associated actions @xmath98 differ by less than 2% for the two regimes . all the conclusions presented in section [ application ] remain valid for the non - approximating regime . we therefore claim that our work based on eq . [ fw_hamiltonian ] is a valid analysis for the chemical master equation . gillespie , j. chem . phys . * 113 * , 297 ( 2000 ) . h. feng , k. zhang , j. wang , chem . sci . * 5 * , 3761 ( 2014 ) .
Abstract:

Answer: we discuss the transition paths in a coupled bistable system consisting of interacting multiple identical bistable motifs . we propose a simple model of coupled bistable gene circuits as an example , and show that its transition paths are bifurcating . we then derive a criterion to predict the bifurcation of transition paths in a generalized coupled bistable system . we confirm the validity of the theory for the example system by numerical simulation . we also demonstrate in the example system that , if the steady states of individual gene circuits are not changed by the coupling , the bifurcation pattern is not dependent on the number of gene circuits . we further show that the transition rate exponentially decreases with the number of gene circuits when the transition path does not bifurcate , while a bifurcation facilitates the transition by lowering the quasi - potential energy barrier .


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: it is well accepted today that the fractional quantum hall effect results from the formation of a strongly correlated spin polarized electron liquid in the lowest landau level , that crystallizes below filling 1/7 @xcite . it is intriguing , however , that a mean field charge density wave solution exists at all fillings , called a hall crystal ( hc ) , which yields in a natural way the odd - denominator rule that characterizes the effect.@xcite the state is characterized by a unit cell that contains a fractional charge of even denominator , such as e/2 . the detailed investigation of this class of mean field states was stopped largely due to their higher energy than that of the liquid , and the very small correlation correction obtained by yoshioka and lee ( yl ) for the wigner solid ( ws ) , which accomodates one whole electron per plaquette.@xcite however , while in the latter case there is little overlap of electrons from neighboring cells , in the former the charge density forms continuous ridges between cells opening the way for an appreciable increase of the cohesive correlation energy.@xcite this idea was first explored in the simplest case of @xmath0 filling.@xcite the results confirmed the effect yielding a second order correlation correction an order of magnitude larger than for the best ws state . the same line of thinking was also recently considered by mikhailov.(@xcite ) .    in this work we report results for the correlation energy of the hc at @xmath1 filling . the work rests on an analytic solution of the mean field problem constructde using symmetry considerations . we show that the electron - electron interaction breaks the single particle degeneracy of the lowest landau level ( lll ) , splitting it into three remarkably flat bands . all bands contain the same number of states so that only one is filled while the other two are empty , separated by a large gap . while the hc unit cell traps 3/2 flux quanta of the original field , it is only pierced by one flux quanta of the reduced field @xmath2 where @xmath3 and @xmath4 are the fields at @xmath1 and @xmath0 filling factor , respectively . the narrowness of the occupied band and this latter fact suggest that the magnetic wannier states are approximate solutions of the problem and behave as nearly free quasiparticles filling the lll of the effective field @xmath5 much as the composite fermion model predicts.@xcite    we compute the correlation energy using the yl method , that is , evaluating the second order perturbation correction to the hartree - fock ( hf ) solution . two methods are used , one employing a monte carlo approach for the space integrals over a large sample , and the other by computing in a discrete momentum space using a set defined by periodic boundary conditions applied to a comparatively small sample . as for the half filling case , both methods yield an energy correction of about 10% , an order of magnitude larger than that for the ws state@xcite , confirming the earlier conjecture that at all fillings the high electron overlap of the hc enhances significantly the cohesive energy corrections@xcite , thus re - opening the case for this state as the proper precursor in a first principles perturbative approach to the true gound state of the system in the thermodynamic limit . one further outcome of our calculations is that over 50 particles would be required in order to approach convergence , way above current permissible sizes used in first principles numerical evaluations normally performed with up to 12 particles , number that has been extended to treat at most 27 particles @xcite . this may be understood again looking at the large inter cell overlap through hexagonal ridges , which seems to point out the importance of cooperative ring exchange effects in the system.@xcite to be operative , these need large sample sizes and would be missed by such numerical diagonalizations of the hamiltonian . if the true ground state retains the essential features of its hf approximation , then one would expect the hall crystal to be the proper mean field precursor due to its higher correlation energy than other mean field solutions.@xcite    in sec . 2 the single particle hartree - fock hamiltonian is diagonalized analytically . section 3 is devoted to the second order correction os the mean field solution . finally , in section 4 we present our conclusions . details of the calculation are given in appendices a - d . in particular , in appendix b a formula is derived expressing the single particle hartree - fock hamiltonian in terms of the magnetic translations associated to an arbitrary periodicity lattice of the  hf  problem . this formula , up to our knowledge , is absent in the literature . we consider a 2deg in a strong perpendicular magnetic field . we are interested in an analytic treatment of the hartree - fock problem at filling 1/3 . we start out by writing the hartree - fock hamiltonian in the compact form    @xmath7    where the coefficients @xmath8 are given by @xmath9 where @xmath10 is the background dielectric constant . the magnetic translation operators t are defined in appendix a , while form ( 1 ) is derived in appendix b. the particle density in real space @xmath11 is assumed to be periodic under displacements covering the triangular lattice defined by the vectors @xmath12 the fourier components of the density@xmath13are defined as @xmath14 where @xmath15 is the unit cell area @xmath16 here @xmath17 is a unit vector normal to the electron gas plane . through this cell traverse a flux @xmath18 in units of the magnetic flux quantum @xmath19 . associated with the above real space lattice is the reciprocal lattice @xmath20      we next show that it is possible to find a basis in which the matrix representation of ( 1 ) has a diagonal form composed of simple 3-dimensional blocks . in addition , the functions have such a structure that they automatically furnish the translation symmetry of the total density over the lattice ( 2 ) . for this purpose we consider the basis functions @xmath21defined in appendix a , constructed over the lattice with primitive vectors @xmath22 the magnetic translations then have the form @xmath23 since the flux piercing the unit cell of the lattice ( 2 ) is not an integral number of flux quanta , the set of translation operators obtained by varying the integers @xmath24and @xmath25 in ( [ trans ] ) do not commute and it is not possible to find common eigenfunctions to all of them . the basis we shall construct decomposes in a set of three dimensional subspaces , closing each of them under the action of translations for all values of @xmath26 . a first step in finding the basis is to define a set of eigenfunctions of a translation in the vector @xmath27 which is a period of the density . for a given value of the momentum @xmath28 such orbitals may be written as @xmath29 using the formulas in appendix a it can be readily proven that these functions satisfy the eigenvalue relation @xmath30 the two values of @xmath31 appearing in these equations will play an important role in what follows . they will allow us to impose the periodicity of the density under the shifts in @xmath32 and @xmath27 in spite of the impossibility of obtaining a basis of eigenvectors of all the magnetic translations in the lattice . it should be stressed that the range of values of @xmath28 defining independent functions in the new basis have been reduced in half , the two values of @xmath31 compensating for this reduction . the restriction comes from the singular property of the basis @xmath33 that a magnetic translation is fully equivalent to a shifting of the momenta argument as implied by the relation @xmath34 where @xmath35 is a pure phase factor . see appendix a and ref ( 16 ) for the justification of this relation . it directly implies that the shift done in @xmath36 in constructing the new basis precisely changes the momentum @xmath28 of the particle in @xmath37 . therefore , what has been basically done is to form linear combinations of the original functions associated with different values of the momentum . in the new basis the magnetic translation in @xmath32 has the simple effect @xmath38 that is , it merely changes the sign of @xmath31 and multiplies it by a phase factor .    as the next step let us employ the fact that , although the functions are not eigenvectors of translations in @xmath32 , they are eigenfunctions of the double sized translations in @xmath39 . this is because its effect , when considered as two consecutive shifts in @xmath40 have the simple result of making two consecutive changes of sign of @xmath31 that reproduce the original function . therefore , if for @xmath28 and @xmath41 fixed we construct the triplet of states formed by the function @xmath42 and the other two obtained by a pair of successive translations in the vector -@xmath43 , the operation of performing a translation in an arbitrary multiple of these vectors leaves the triplets invariant . a specific linear combination within each triplet which is also an eigenfunction of the translation in -@xmath43 can be obtained by constructing the new basis @xmath44 after solving the linear equations for the constants @xmath45 obtained by imposing the condition that these functions be solutions of the eigenvalue problem @xmath46 one finds @xmath47 substituting in ( [ functions ] ) yields @xmath48 the last relation expresses that , modulo a phase factor , the states of the new basis are equivalent upon a shift of their momenta @xmath28 in any linear combination with integer coefficients of half the unit cell vectors of the reciprocal lattice corresponding to the periodicity of the density . the periodicity of the states under the shifts in @xmath49 was discussed above , and the one related with @xmath50 similarly follows from the relation ( [ phase ] ) in appendix a , expressing the equivalence of a magnetic translation with a shift in momentum . the functions just defined have an alternative and more compact form given by @xmath51 where the effective momenta and the new elementary lattice of vectors @xmath52 over which the sum is performed are given by @xmath53 the double sum ( [ altern ] ) can be evaluated to obtain an explicit formula in terms of the elliptic theta functions as follows @xmath54 where @xmath55    let us inspect now the action of a magnetic translation by @xmath56 on the new functions . if such a transformation leaves the triplets invariant , then the matrix reduction of the hartree - fock hamiltonian will follow . one has , @xmath57 but after using ( [ commutation ] ) for changing the order of the two operators within the sum , it follows that @xmath58,\sigma ) } ( % \mathbf{x}),\ ] ] where the square bracket defines the number among the set \{-1,0,1 } that is equivalent , modulo 3 , to the integer in the argument .    for fixed @xmath28 and @xmath31 the matrix elements of the hamiltonian ( [ hamiltonian ] ) in the new basis can readily be found to have the form @xmath59}. \label{heigen}\ ] ] the problem has thus been reduced to the self - consistent diagonalization of a three dimensional matrix for each value of momenta @xmath28 and parameter @xmath31 . the basis can be checked to have the following set of transformations properties , @xmath60,\sigma ) } ( \mathbf{x } ) , \nonumber \\ t_{3\mathbf{c}_{1}}\chi _ { \mathbf{p}}^{(r,\sigma ) } ( \mathbf{x } ) & = & \ \exp ( -3i\ \mathbf{p.c}_{1}\,)\chi _ { \mathbf{p}}^{(r,-\sigma ) } ( \mathbf{x } ) , \nonumber \\ t_{3\mathbf{c}_{2}}\chi _ { \mathbf{p}}^{(r,\sigma ) } ( \mathbf{x } ) & = & \ \sigma \exp ( -3i\ \mathbf{p.c}_{2}\,)\chi _ { \mathbf{p}}^{(r,\sigma ) } ( \mathbf{x } ) , \nonumber \\ i\ \chi _ { \mathbf{p}}^{(r,\sigma ) } ( \mathbf{x } ) & = & \chi _ { -\mathbf{p}% } ^{(-r,\sigma ) } ( \mathbf{x } ) ,   \nonumber\end{aligned}\ ] ] where the parity transformation @xmath61 is defined as usual by @xmath62  as shown in appendix a , from these symmetry properties it follows that the density associated with any slater determinant constructed by selecting one orbital within each triplet has exact periodicity under shifts in all vectors @xmath63 .    in order to find the solution of the mean field problem by an iterative technique it is sufficient to make an ansatz for the density in the first step , and then diagonalize numerically the matrices for a sufficiently high partition of the reduced brillouin cell momenta @xmath28 defined by ( [ brillouin ] ) . by selecting the normalized lowest energy state within each three dimensional quantum mechanical problem , the fourier components of the density corresponding to the new step should be constructed . following its definition , it can be done by means of the following expression , @xmath64}% { \large . } \nonumber\end{aligned}\ ] ] this formula can be obtained by evaluating the gaussian integrals appearing after substituting the expansions defining the functions @xmath65the coefficients @xmath66 , @xmath67 and @xmath68 , define the components of the eigenvectors of the single particle hf hamiltonian in the basis of states @xmath69 . they fix the wave - functions of the filled band and the empty bands as @xmath70 here @xmath71 labels the filled band in each triplet and @xmath72 label the two empty bands .      in order to proceed within an analytical context let us consider the observation from former numerical studies , that the particle density when the number of flux quanta piercing the unit cell is a half integer essentially vanishes at all lattice points.@xcite then , let us first assume that the density rigorously vanishes at this set of points . if so is the case , the wave - function of any of the filled states should then also vanish at those points . this requirement follows from the fact that the hartree - fock particle density is a sum over the individual densities of all occupied orbitals ,    @xmath73    we can then use this property to fix the coefficients of the wave functions within each triplet .    after imposing the vanishing conditions at the origin , the coefficients defining the functions ( [ solution ] ) are fully determined and take the form @xmath74 note that the coefficients @xmath75 are all independent of @xmath31 . this completes the definition of our functions . that they are true solutions of the hartree - fock problem is proven in appendix c.    the particle density may now be computed replacing these functions in eq . ( [ charge ] ) . the real space particle density thus obtained is shown in fig.[fig : densi ] . a main property to be noticed in this figure is the formation of sharp hexagonal channels surrounding the low density regions at the center of which the vanishing density occurs . these structures mark the difference with the wigner solid whose unit cell encloses one flux quanta . the charge density in this latter instance is made up essentially of well localized gaussians centered at each lattice point . in our case there is strong overlap , suggesting that cooperative ring exchange involving many unit cells is a large contribution to the correlation energy . further , the insertion of the calculated density in the eigenvalue equation associated with the matrix representation of the hamiltonian in each triplet , eqs . ( [ poten ] ) and ( [ heigen ] ) , allows for the evaluation of the one particle spectrum of the system . as it was expected , three energy bands appear , each associated with a value of the index r and covering the full range of @xmath28 within each triplet . we also note that states associated with @xmath76 turn out to be degenerate . the bands dispersion relations are illustrated in fig . [ fig : bands ] . note the narrowness of the bands as compared with the gaps separating them . this fact leads to the idea that in this mean field approximation the coulomb interaction , although breaking the degeneracy of the first landau level , it reorganizes the states in three equally populated separate sets that again are aproximately degenerate , as if they were landau levels of a renormalized problem . since at 1/3 filling one band is full and the other two empty , one expects the electrons to behave dynamically similarly to filling one , except for a different effective mass , as the composite fermion model suggests.@xcite the magnetic wannier states are expected to be approximate solutions , playing the role of the angular momentum states in the non - interacting problem at filling one . these properties seem to corroborate the possibility of tracing a link between weakly interacting composite fermions and the bloch or wannier orbitals in the mean field solutions considered earlier by one of us.@xcite to finish this section we note that the mean field energy per particle obtained from the above formalism for the solution we have constructed confirms the value reported in an earlier numerical computation,@xcite @xmath77 in order to obtain the energy correction to second order of perturbation theory we proceed similarly as in ref ( 8) , starting with the expression@xcite    @xmath78    here @xmath79 , @xmath80 are the mean field slater determinant and total hartree - fock energy associated with the ground state , respectively , and @xmath81 is the projection of the exact many particle hamiltonian onto the first landau level . the many particle excited states @xmath82 are slater determinants constructed with the basis ( [ solution ] ) , mixing states in the filled band with states in the empty bands . it follows that @xmath83 , a property that allowed to write the last equality in ( [ foryosh ] ) . in the second quantized representation the hamiltonian @xmath81 will have non - vanishing matrix elements linking the hf state and excited states of the form @xmath84 only , where @xmath85 creates an electron of quantum numbers @xmath86 , etc . the index @xmath87 is a shorthand notation for the set of two pairs of filled @xmath88 and empty ( @xmath89 ) electron states , where @xmath90 and @xmath91are the set of all states in the landau level and the filled ones , respectively . the indices @xmath92 @xmath93 and @xmath94 for @xmath72 denote the quantum numbers of filled and empty states , respectively@xmath95 the total energy of an excited state is @xmath96 . then , the second order correction can be rewritten in the form @xcite @xmath97    where the total projected hamiltonian is @xmath98 the matrix elements of the coulomb interaction are given by @xmath99 where we have used the shorthand notation @xmath100 @xmath101 . by using the anti - commutation relations @xmath102=\delta _ { \alpha \mathbf{,\alpha } ^{\prime } } $ ] , formula ( [ foryosh ] ) can be expressed as    @xmath103    where the two particle states @xmath104 are defined by @xmath105 the pairs @xmath106and @xmath107 are considered as unordered .      as pointed out above the single particle bands are remarkably flat . we can use this property to simplify the calculation of the energy correction . first , we approximate the filled band energies appearing in the denominator of ( [ e2 ] ) by their mean value , @xmath108 in addition , and slightly more crudely , we substitute the energies in the excited bands by a common energy equal to half the sum of the mean energies of the two bands , @xmath109 this last approximation is taken in view of the small relative gap separating these two bands .    with this simplifying substitution ( [ e2 ] ) can be expressed in the simpler form @xmath110 where @xmath111 is the projection operator on the subspace of states of the filled band and @xmath112 the projector associated to the subspace of states formed with the union of the empty bands . the projectors have the following expression in terms of the corresponding densities @xmath113 where @xmath114 is the ordinary kronecker delta and @xmath115is the projection operator onto the first landau level , defined in appendix a. these expressions can be obtained from formula ( [ onepart ] ) in appendix b. after evaluating a few spatial integrals , the following formula for the correlation energy per particle @xmath116 is obtained @xmath117    further progress in the evaluation of the energy correction is aided by noting the presence of an exponential factor in the squared fourier wave - vectors , allowing to take just a few of them only , for sufficient convergence . we thus assume the 36 fourier components associated with the shortest values of @xmath26 are nonzero , only . of these , those associated with the longest @xmath26 are slightly altered to assure that the sum rule@xcite @xmath118 is satisfied . this condition is used to cancel various fictitious divergences in the formula defining the correlation energy @xcite . then , we obtain    @xmath119    where the order parameters of the filled and empty bands have been defined as usual , @xmath120        the integrals in eq . ( [ correc1 ] ) were calculated using the monte - carlo method with a variable number of sample points in the range @xmath121 . the region of integration was taken as a square of size @xmath122 in units of the magnetic length @xmath123 centered at the origin . at 1/3 filling this sample contains 132 particles and traps 396 flux quanta . the results for the correlation energy as a function of @xmath124 are shown in fig . [ fig : data1 ] . in units @xmath125 the estimated value for the correlation energy is @xmath126 adding this correction to the hartree - fock energy ( [ energy ] ) one obtains the value @xmath127 , to be compared with the result of yoshioka and lee with its own correction @xmath128 , and that of laughlin @xmath129 , all in the same units . the lack of a definite convergence pattern provided by the method just described prompted us to try a different procedure . using the formulas derived in appendix d it is possible to recast eq . ( [ correc1 ] ) to give it the form ,    @xmath130    where the function @xmath131 is defined by    @xmath132 , \end{aligned}\ ] ]    the form of the pure phase factors @xmath133 and @xmath134 , the momenta @xmath135 and the special reciprocal lattice vectors @xmath136 are all specified in appendix d. as before , the functions @xmath137 are the coefficients determining the single particle hf excitations . we remark that once the coefficients for the filled band are known the other set of coefficients for the empty bands can be evaluated analytically as the two eigenvectors of the ( 3x3 ) matrix representation of the single particle hf hamiltonian , orthogonal to the vector @xmath138 .    in order to evaluate numerically expression ( [ e2alter ] ) we restrict the hilbert space of the single particle hf  problem by defining a cell of sides @xmath139 , @xmath140 with @xmath141 , @xmath142 integers , and imposing periodic boundary conditions over its borders . using the symmetry properties ( [ symmetry ] ) @xmath143    these constraints can be translated into the relations    @xmath144    they restrict the values of the quasimomentum @xmath28 to a discrete set , as expressed in compact form by the condition @xmath145 where @xmath146 and @xmath147 are defined in appendix a. the cell chosen traps @xmath148 flux quanta and contains @xmath149 particles . the integrals over the continuum in the brillouin zone of eq . ( [ e2alter ] ) are now restricted to summations over a finite number of values of @xmath28 . formula ( [ e2alter ] ) was first evaluated by for an elongated sample with  @xmath150 and @xmath151 varying from @xmath152 to @xmath153 .  in order to keep the computing time within our allowed practical limits , the values of @xmath154 and @xmath155 in the expression defining the summation argument @xmath156 were restricted to the regions @xmath157 and @xmath158 the results of this calculation are shown as empty circles in fig . ( [ fig : data2 ] ) . note that the correlation energy per particle grows with n , approaching rather slowly the thermodynamic limit whose value @xmath159 exceeds the result ( [ e2num1 ] ) of the previous method . the energy correction was also calculated for samples with an aspect ratio close to one , taking @xmath160 , up to @xmath161 an upper bound set by computing limitations . the results are shown in fig . ( [ fig : data2 ] ) by filled circles .  as can be observed ,  the energy values for those samples grow similarly with @xmath162 than those associated with elongated samples . this property indicates that these more appropriate samples also need values of @xmath162 greater than five in order to approach convergence . however , already at @xmath163 the number of particles involved is @xmath164 , considerably higher than the numbers handled in exact diagonalization . the most relevant finding of this last method of evaluation is the need for very high sample sizes in order to approach the thermodynamical limit . this outcome strongly supports the relevance of  long distance correlations in the attainment of the true character of the ground state through numerical computations .    . the circles indicate the calculated values for elongated samples of the form @xmath165 . the black dots signal the values corresponding more appropriate samples with form defined by @xmath166,width=336 ] an analytic solution of the hartree - fock problem for a 2deg at filling @xmath1 and 1/2 electron per unit cell is found by employing symmetry considerations and a special complete set of common eigenfunctions of the magnetic translations . a triplet of bands arises that turn out surprisingly flat as functions of the two dimensional quasimomentum . the coulomb interaction breaks the first landau level in three narrow sub - bands resembling effective landau levels of composite fermions , one of which is fully occupied and the other two , empty . the energy per particle reproduces an earlier numerical result for this quantity @xcite . we find that the charge density forms hexagonal rings throughout the lattice , suggesting strong long range correlations owing to cooperative ring exchange . an evaluation of the correlation energy using second order perturbation theory yields a correction an order of magnitude larger than that for the localized single - particle features of the usual wigner crystal state,@xcite and lowers the energy to make it comparable with the lowest values obtained using other models . our results also suggest that the thermodynamic limit in a cell with aspect ratio near one is not reached if less than 50 particles are included . this work was supported in part by fondecyt grants 1020829 and 7020829 , the ictp associateship program , and the third world academy of sciences . helpful discussions with a. gonzalez , g. baskaran , n. h march and k. esfarjani are also acknowledged . consider a 2deg constrained to move in a plane of area @xmath167 under a perpendicular magnetic field @xmath168 . a useful basis set of single particle bloch - like states in the lowest landau level can be defined in terms of linear combinations of the normalized zero angular momentum eigenfunction    @xmath169    in the compact form @xcite ,    @xmath170    @xmath171 due to its role in the above definition the function @xmath172 is called the  seed  function . the sum runs over all integers @xmath173 defining a planar lattice @xmath174 , through @xmath175 where the unit cell intercepts one flux quantum , so that @xmath176 the magnetic translation operators @xmath177 acting on any function @xmath178 are defined by @xmath179 where the vector potential is assumed in the axial gauge @xmath180 and the electron charge @xmath181 is taken with its negative sign . this basis was employed before to obtain exact mean field solutions of the related problem at filling 1/2.@xcite for arbitrary vectors @xmath182 and @xmath183 the translation operators satisfy the commutation relation @xmath184 as it may be easily verified , the functions @xmath185 satisfy the eigenvalue equation    @xmath186    arranged in an arbitrary slater determinant these functions are exact solutions of the hartree - fock problem.@xcite this strong property happens because the hf single particle hamiltonian commutes with all translations leaving @xmath174 invariant . @xcite the functions ( [ fi ] ) are common eigenfunctions of the commuting magnetic translations.moreover , the set of eigenvalues ( [ eigen2 ] ) uniquely determines them . therefore , the hf hamiltonian associated with the slater determinant can not change those eigenvalues and the @xmath187 should be eigenfunctions .    finally , let us show that the effect of an arbitrary translation on the basis functions is equivalent to a shift in the momentum label , modulo a phase factor.@xcite operating twice with the translation operator involving an arbitrary vector @xmath188 and a vector in the lattice @xmath189 and using eqs . ( [ commutation ] ) and ( [ eigen1 ] ) one readily gets ,    @xmath190    again using ( [ eigen1 ] ) we have    @xmath191 then , taking into account that the set of eigenvalues defines uniquely the wave - functions modulo a phase , it follows that    @xmath192    from which also follows ,    @xmath193    that is , a magnetic translation is equivalent to a shift in the quasi - momentum . any function@xmath195 in the first landau level , and its inverse , can be represented as    @xmath196    where @xmath197 @xmath198 with @xmath199 with @xmath200 the number of flux quanta in the system area , @xmath201 is the brillouin zone defined by the unit cell vectors @xmath202    the orthogonality and completeness relations in the first landau level of the modified functions take the forms    @xmath203      lets us now verify that the functions @xmath204 constructed to diagonalize in blocks the hf hamiltonian in sec . 2 are simply constant phase factors multiplied by a particular kind of basis functions @xmath206 . these differ from the orbitals considered in sec . 2 and defined by the unit cell vectors ( [ bas32 ] ) in that their basis vectors are given by @xmath207 to start out it will be useful to consider the @xmath194-transform of  a translated  seed  function@xmath208 for this purpose we notice that it is possible to fix the second argument of the projector operator to be equal to the translation vector @xmath209  it follows that the projector operator in the first landau level can be rewritten as @xmath210 where the last equality is merely the definition of the projector operator in coordinate space . therefore , these relations determine the following @xmath206-transform of the translated  seed  function @xmath211 next we will consider the definition @xmath212 and represent the vectors @xmath213  in terms of alternative vectors @xmath189 and the indices @xmath214 and @xmath215 defined through the expressions    @xmath216    the symbols @xmath217 and @xmath218 are specified by the relations    @xmath219 = { \huge \ { } \begin{tabular}{l } $ -1 $ \ \ \ \ $ \ if$ $ \ \ m_{2}=-1 $ $ \ modulo(3)$ \\ $ \ \ 0 $ \ \ \ \ \ $ if$ $ \ \ $ \ \ \ 0 $ $ \ modulo(3)$ \\ $ \ \ 1 $ \ \ \ \ \ $ if$ $ \ \ $ \ \ \ 1 $ $ \ modulo(3)$% \end{tabular } , \\ _ { 1 } & = & \frac{m_{1}-\left\langle m_{1}\right\rangle } { 2 } , \\ \mathbf{\ell } _ { 2 } & = & \frac{m_{2}-\left [ m_{2}\right ] } { 3\ } .\end{aligned}\ ] ]    the use of these alternative definitions allows to write ( [ csinew ] ) in the form    @xmath220    now , it can be noticed that the above expression differs form the definition ( [ fi ] ) only in that the  seed  function is changed by a superposition of itself translated to the points @xmath221 however , it was shown in ref . ( 16 ) that the basis functions are completely independent of an arbitrary change in the  seed  for any other function in the first landau level . an alternative proof of this curious property is possible .  to see it let us represent the  seed  functions in ( [ indep ] ) by their @xmath222-transform ( [ transform ] ) , yielding @xmath223 after employing the character of eigenfunctions of the operators @xmath224 that the functions @xmath225 have . but now the following identity can be employed    @xmath226 where the @xmath227 function is non - vanishing for equal arguments modulo a vector of the reciprocal lattice @xmath228 this relation allows to show the desired connection between the basis functions @xmath204and @xmath229 ,    @xmath230 where the phase factor is given by @xmath231      here we show that the single particle density associated with a slater determinant formed with a set of functions , each one corresponding to an arbitrary linear combinations of the triplet of orbitals characterized by the indices @xmath232 is periodic over the lattice @xmath233 to consider this question , let us use the basic property that the density of a many electron slater determinant constructed with orthogonal orbitals is the sum of the individual densities of each orbital . writing the selected linear combinations as @xmath234 where the coefficients @xmath235 are arbitary , and considering the symmetry properties ( [ symmetry ] ) , it follows that    @xmath236    one can show in a similar way that @xmath237 here we shall derive the general formula expressing the single particle hartree - fock hamiltonian as a sum of  magnetic translation operations . the arguments of these operators are spatial vectors determined by the reciprocal lattice associated with the fourier components of the density . it is worth remarking here that the lattice considered in this appendix is not restricted in any way . therefore , the representation of the single particle hf hamiltonian is valid for a general though periodic mean field problem .    consider the action of the kernels of the direct and exchange interactions defining the hf hamiltonian  over any of the elements @xmath238 of the complete basis functions as follows    @xmath239    where @xmath240 is the jellium background charge density making the overall system neutral . the representation we shall derive is valid for both direct and exchange kernels and will be discussed below for each case separately . the potential term representing the direct coulomb interaction can be written in its fourier transform representation as follows    @xmath241    then , using the following representation of the kernel for the projection operator in the first landau level    @xmath242 we obtain , @xmath243 recalling definition ( [ fi ] ) we obtain for each translation term @xmath244 this last relation can be further simplified by employing relation ( [ transl1 ] ) and the similar property @xmath245 leading to the formula    @xmath246 after adding corresponding the integrals over all translations @xmath247 one obtains @xmath248 since the equality is valid for any element of the complete basis @xmath249 it then follows , @xmath250      the derivation of the analogous representation for the exchange interaction kernel is more involved and needs for some special properties of the one particle density matrix . then , let us initially consider these properties . the definition of the one - particle density matrix in terms of the slater determinant @xmath251 can be transformed by performing a simultaneous magnetic translation operation in a vector @xmath63 of the periodic lattice over all the particle coordinates . since this is assumed to be a symmetry transformation of the system , this map should leave the many particle state invariant . thus by assumption    @xmath252    an expression that furnishes the transformation law of the one particle density matrix under spatial shifts in the vectors of the periodic lattice . below we will show that the fourier components of the density completely determine the whole one - particle density matrix .  for this purpose , let us use the new variable @xmath253 in the one - particle density matrix , so that @xmath254 where the tilde over @xmath255 underlines the different functional expression of the new definition . it now transforms under spatial shifts as @xmath256 therefore , the function    @xmath257 is fully periodic in the variable @xmath258 under the lattice shifts . that is @xmath259 then , fourier expanding @xmath260  leads to the following expression for @xmath261    @xmath262 where    @xmath263    further , let us consider the density written as the sum @xmath264 and the fact that each @xmath265 in the lower landau level can be expanded as a linear combination of the complete basis functions@xmath266 defined over a lattice with one flux quantum per unit cell.@xcite then , a generic term in the integral in ( [ elele ] ) has the form @xmath267 since the function @xmath172 is a gaussian , its magnetic translations ( [ magntra ] ) have also this character , and thus the integrals in ( [ elele ] ) can be calculated explicitly . we obtain , @xmath268 but since the phase factor here is independent of the indices @xmath269 ,  the summation of all the contributions in ( [ elele ] ) leads to the relation    @xmath270    which expresses the interesting result that the full one - particle density matrix is completely determined by the particle density . we now consider the exchange term ( [ exchange ] ) . using the above results the integrals in    @xmath271 can be calculated  explicitly . for this purpose , let us consider again the integral of a generic term @xmath272 in the sum over @xmath189 defining the functions @xmath194 . the integral for this term is then reduced to two simple integrals , which after some algebra can be explicitly evaluated to obtain @xmath273 henceforth , after adding all terms for different values of @xmath274 follows @xmath275 using the independence and completeness of the basis formed by the @xmath276 s in the first landau level , we thus obtain @xmath277 combining this results and the previous one ( [ directterm ] ) the following representation for the one - particle hartree - fock hamiltonian in the first landau level follows ,    @xmath278    with the coefficients defined by    @xmath279 in this section we prove that the eigenfunctions that vanish at all latice points @xmath63 are solutions of the hf equations@xmath280 for this purpose , let us consider the fourier components @xmath281 of the configuration space periodically extended density operator    @xmath282    we will examine the projection of this operator in the first landau level . for this aim it is useful to consider the set of operators @xmath283 where the indices @xmath284  take the values 1,2 . these quantities satisfy the commutation relations    @xmath285 & = & 0 , \\ \left [ p_{i},p_{j}\right ] & = & -\frac{\epsilon _ { ij}}{r_{0}^{2 } } , \\ \left [ g_{i},g_{j}\right ] & = & \frac{\epsilon _ { ij}}{r_{0}^{2}}.\end{aligned}\ ] ]    with their use , the spatial fourier modes can be identically represented as    @xmath286    where the 3d - vectors @xmath287 and @xmath288 are defined as @xmath289 and @xmath290 . but the two components @xmath291 are linear combinations of the rising and lowering operators which transform functions of a landau level in functions of the higher and lower contiguous levels . therefore , the projection onto the first landau level of the fourier component can be obtained by simply omitting the operator @xmath287 in the above formula . then , the projection of the density operator will be @xmath292 where a new set of spatial lattice vectors was defined in terms of the fourier wavevectors , @xmath293 let us consider now the identity @xmath294 where the three coefficients @xmath295 are fixed by the condition that the linear combination of the functions of the two triplets of  functions @xmath296for @xmath76 vanish at the origin @xmath297 projecting the above identity onto the first landau level by acting with @xmath298 at the left it follows @xmath299 implying that the projected periodic density has the defined functions as eigenvectors with zero eigenvalues . let us inspect further the commutation properties of the projected density operator . using the representation of the single particle hamiltonian in terms of the magnetic translations found in appendix b and the relations @xmath300 it follows that    @xmath301   & = & \sum_{\mathbf{x}% ^{*}}\sum_{\mathbf{y}^{*}}u\text{(}\mathbf{x}^{*})\left [ t_{\mathbf{x}% ^{*}},t_{\mathbf{y}^{*}}\right ]   \\ & = & \sum_{\mathbf{x}^{*}}\sum_{\mathbf{z}^{*}}u\text{(}\mathbf{x}^{*})(\exp ( % \frac{ie}{\hbar c}\mathbf{a(x^{*}).z}^{*})-\exp ( -\frac{ie}{\hbar c}\mathbf{% a(x^{*}).z}^{*}))t_{\mathbf{z}^{*}}.   \nonumber\end{aligned}\ ] ]    considering that @xmath302 is the coefficient of the superposition of translations giving the single particle hf hamiltonian which is a function of @xmath303 , and using the inversion symmetry of the solution under consideration @xmath302(@xmath304(-@xmath305 one finds @xmath306 = 0.\ ] ] hence , acting with the hamiltonian on ( [ vanish ] ) and employing the definition of the projection operators and the periodic density operator it follows that    @xmath307    finally , taking into account that the set of functions @xmath308 for all values of the lattice ( [ perlatt ] ) are linearly independent @xcite , it follows that    @xmath309    in addition , since the coefficients of the linear combinations of the triplet functions for the two signs of @xmath31 are uniquely defined by the vanishing condition of the function at the origin ( or any point of the lattice @xmath63 ) , we have , @xmath310 where @xmath311 is here a normalization constant . finally , we obtain @xmath312 that is , the special functions that vanish at all latice points are eigenfunctions of the one particle hf hamiltonian . in this appendix we sketch the derivation of the formula ( [ e2alter ] ) employed in the evaluation of the sample size dependence of the correction to the energy per particle . as before , we use the shorthand notation @xmath313 the relations expressing the functions @xmath314 in terms of @xmath315times a phase factor has also been employed in order to define the new coefficients @xmath316  by representing the coulomb potential by its fourier transform , and substituting all the expansions of the functions @xmath317 in terms of  magnetic translations of the  seed  in the matrix elements ( [ matrix ] ) , the resulting gaussian integrals for each term can be evaluated . after some algebra it is possible to obtain the expresion      where the summation vectors @xmath136 and the norm functions are defined as follows , @xmath319 and , as before , @xmath320 equals one when the arguments are same modulo vectors of the class @xmath136 defined above , and zero otherwise@xmath95 it is possible to show now that the ratio multiplying this latter function a simply a phase factor . for this purpose let us consider the definition of the @xmath321        @xmath323 further , let us recall the formula connecting a translated function of the basis @xmath194 and the same function for a shifted momenta ( [ phase ] ) . after considering this relation evaluated at @xmath324 , the following expression can be obtained
Abstract:

Answer: an analytic solution of the hartree - fock problem for a 2deg at filling 1/3 and half an electron per unit cell is presented . the coulomb interaction dynamically breaks the first landau level in three narrow sub - bands , one of which is fully occupied and the other empty , as in the composite fermion model . the localized orbitals associated to the bloch like single electron wavefunctions are nearly static , resembling the angular momentum eigenstates within a landau level for non - interacting fermions . strong correlations are expected owing to the large charge density overlap between neighboring plaquettes . a numerical evaluation brings the cohesive energy close to that of the best present day models . it is also found that correlations are long range , requiring over 50 particles spread over a finite sample to approach convergence . since presently allowed exact calculations are far from this number , the question of how relevant the considered wave - function is for the description of the ground state of the 2deg system remains open . pacs numbers : 73.43.cd , 73.43.-f       * on leave of absence from _ theoretical physics group , instituto de ciberntica , matematica y fsica , calle e , no . 309 , vedado , la habana , cuba . _     preprint ictp , ic/2003/60(2003 )