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: recent observational surveys ( e.g. , cosmos , goods , gems , sings , sdss , * ? ? ? * ; * ? ? ? * ; * ? ? ? * ; * ? ? ? * ; * ? ? ? * ) have allowed us for the first time to statistically explore issues of galaxy formation . one of the immediate challenges to disk formation theory that is presented by these observations comes in the form of evidence that there has been little change in the sizes of disk galaxies since @xmath2=1 , despite an expectation that disks should be growing in size over this time . observations tell us that there exist large disks by @xmath2=1 ( before the universe was even half of its present age ) , and suggesting that these disks must be assembled prior to this epoch ( e.g. , * ? ? ? * ; * ? ? ? * ; * ? ? ? * ; * ? ? ? * ; * ? ? ? * ; * ? ? ? * ; * ? ? ? * ; * ? ? ? * ; * ? ? ? * ; * ? ? ? * ; * ? ? ? recent observations lend strong support to the existence of disks back to even higher redshifts @xcite , and for disk assembly at even higher redshifts in the form of clump - chain galaxies in the hubble ultra deep field @xcite . results from large surveys suggest that there has been only weak evolution in disk sizes since @xmath2=1 . there is no evidence that the size function of disks has evolved back to @xmath2=1 @xcite . @xcite found that the shape of the size function did not evolve with redshift , though the normalization ( or amplitude ) did . this can be interpreted in two ways . first , the normalization will vary if the number density of disks varies with time . however , observations suggest that the number density of disks is constant to @xmath2=1 @xcite . second , if it is assumed that galaxies dim with time ( e.g. , due to a passively evolving stellar population or decline in star formation rate ) , a 1 to 1.5 magnitude dimming since @xmath2=1 could explain the change in normalization , while keeping the size distribution of disks constant @xcite . hence , there is no immediate evidence for a change in the sizes of galaxy disks with time . studies of the evolution in the magnitude  size relationship for disk galaxies back to @xmath2=1 have found similar results . importantly , these studies must use surface brightness evolution to interpret evolution within the magnitude  size plane . this is not straightforward , as an increase in disk scale lengths at a fixed magnitude between @xmath2=1 and @xmath2=0 can be mimicked by a decrease in luminosity at a fixed size @xcite . selection effects are difficult to disentangle , and studies that require high and low @xmath2 samples to adhere to the same selection biases have found that large disks are consistent with no evolution in surface brightness @xcite . later studies concluded that this work was too restrictive , and that a careful treatment of completeness as a function of redshift supports surface brightness dimming over time @xcite . most studies conclude that pure size evolution at a fixed magnitude is ruled out @xcite . instead , the observed surface brightness evolution is best explained by luminosity dimming , with the amount of dimming dependent on galaxy size so that lower mass galaxies have undergone more dimming since @xmath2=1 than massive galaxies @xcite . however , some amount of size evolution ca nt be ruled out , and at least weak size evolution is favored @xcite .    while magnitude is assumed to scale with the mass of a galaxy , it is obvious from the above discussion that magnitude at a given mass is not necessarily constant in time . the more fundamental relation is between size and stellar mass . magnitude at a given stellar mass is likely to increase back to @xmath2=1 due to the fact that there is an evolving star formation rate ( sfr )  mass relation @xcite . @xcite examined the stellar mass  size relation for disk dominated galaxies back to @xmath2=1 , and found weak or no evolution . this result can be interpreted in two ways . first , given the change of @xmath81 magnitude of surface brightness dimming in the same population , the result is consistent with a passively evolving stellar population at a given mass , with no growth of galaxy disks . second , it could imply that galaxies are growing , but that they must grow in such a way as to evolve _ along _ the stellar mass  size relation with time .    on the other hand , galaxy disk formation theory predicts that the sizes of disks , both individual and as a population , should be growing since @xmath2=1 . in the standard picture , gas in a halo conserves its specific angular momentum ( equal to that of the dark matter ) as it cools to form a centrifugally supported disk that grows from the inside out @xcite . in the simplest model , in which the density profile of galaxies is modeled as a singular isothermal sphere ( sis ) , the radius of the resulting disk is proportional to the parent halo virial radius , which grows inversely with the hubble parameter , @xmath9 . for a concordance @xmath0cdm cosmology @xcite , this relation predicts that disks at @xmath2=1 should be nearly a factor of two smaller than their @xmath2=0 counterparts at the same circular velocity , @xmath10 @xcite . @xcite adopted more reasonable assumptions to update the simple theoretical model that suggests that galaxy disks should grow by nearly a factor of two since @xmath2=1 , primarily by adopting nfw density profiles rather than sis profiles . unlike the sis model , the concentrations of nfw profiles were lower in the past . combined with disk stability arguments @xcite , this model predicts a much weaker evolution in disk sizes ( 15%-20% ) at a fixed stellar mass , consistent with the @xcite results . these semi - analytic models ( sams ) do not follow the evolution of individual disk galaxies with time . rather , they examine an instantaneous population of disks . clearly , this is the type of information derived from observations , but it prevents an interpretation of the evolution in individual disk galaxies based on observations of the population as a whole . @xcite instead follow the evolution of individual disks in their sam , and demonstrate that disks tend to grow along the stellar mass  size relation , with only weak evolution , consistent with the results of @xcite . the analytic models discussed above ( e.g. , * ? ? ? ? * ; * ? ? ? * ) have until recently been the only available theoretical tool with which to investigate the evolution of galactic disk sizes ( though see * ? ? ? simulations , with their ability to capture complex gas processes in mergers and subsequent sf , should be an ideal tool that allow a better trace of the distribution of stellar light , while sams must assume that light follows mass . however , simulations of disk galaxy formation in a cdm context have historically produced unrealistic disks that are too compact ( dense ) , too small overall , and rotating too fast at a given radius ( e.g. , * ? ? ? * ; * ? ? ? * ; * ? ? ? * ; * ? ? ? * ; * ? ? ? this failure has been named the `` angular momentum catastrophe . '' this catastrophe has been largely associated with the overcooling problem , in which baryons cool too quickly at early times and become very dense and concentrated at the center of halos before merging . these halos then experience dynamical friction in subsequent mergers , and the resulting disks show the classic signs of the angular momentum catastrophe @xcite . thus , feedback mechanisms at early times have historically been invoked to prevent overcooling with moderate success @xcite but additional artificial exchange of angular momentum can still occur between baryons and dark matter in sph simulations with low numerical resolution @xcite . hence , to avoid the angular momentum catastrophe and produce realistic disks , simulations must have both a physically motivated feedback prescription , and very high numerical ( mass and force ) resolution @xcite . the simulations used in this study are part of a suite of very high numerical resolution disk galaxy simulations that incorporate a star formation ( sf ) and supernova ( sn ) feedback scheme that has been shown to overcome many of the past problems of disk galaxy simulations to successfully match a number of observed properties of galaxies ( as discussed further below , and shown by the results of this paper ) .    in this paper , we use simulations of individual disk galaxies to follow their evolution in the magnitude  size plane and the stellar mass  size plane . we attempt to mimic the observations of these galaxies in a manner similar to the observed population , generating mock surface brightness images to measure disk scale lengths , avoiding the assumption that light must follow mass . we demonstrate that the simulated galaxies have properties similar to real galaxies . in  2 we describe our small sample of very high resolution simulated disk galaxies that span a representative range of masses , merger histories , and spin values . we show that our limited sample is in excellent agreement with the population of galaxies surveyed at both @xmath2=0 (  3 ) and @xmath2=1 (  4 ) . in  5 , we disentangle the evolution of individual galaxies from the evolution seen in the population as a whole . we conclude in  6 . these simulations were run with the n - body + smoothed particle hydrodynamics ( sph ) code gasoline @xcite in a fully cosmological @xmath0cdm context using wmap year 3 parameters ( @xmath11 , @xmath0=0.74 , @xmath12 , @xmath13=0.77 , n=0.96 ) . the galaxies were originally selected from either a 25mpc or 50mpc ( depending on their mass ) n - body simulation with uniform mass resolution throughout , and then resimulated at higher resolution with baryons using the volume renormalization technique @xcite . this technique allows for significantly higher resolution on the central galaxy while also capturing the effect of large scale torques that are thought to deliver angular momentum to the galaxy @xcite . the sf / sn scheme used in these simulations has been shown to effectively regulate star formation efficiency as a function of halo mass , resulting in a stellar mass  metallicity relationship for the simulated galaxies that is in excellent agreement with observations both locally and at high @xmath2 @xcite . this regulation of star formation also leads to realistic trends in gas fractions , with our lowest mass galaxies being the most gas rich @xcite , reproducing the observed incidence rate of damped lyman @xmath14 systems ( qso - dlas ) at @xmath2 = 3 @xcite , and the column densities of both qso - dlas and grb - dlas at @xmath2 = 3 @xcite . of critical importance for the present study is the fact that these simulated disks maintain sufficient angular momentum to match the observed tully - fisher relationship @xcite and produce galaxies with realistic disk sizes , as shown below . the full details of our physically motivated sn feedback implementation were originally presented in @xcite . briefly , the sf prescription ensures that the sfr density is a function of gas density according to the observed slope of the kennicutt - schmidt law , and a sf efficiency parameter , @xmath15 , sets the normalization of this relation . each star particle represents a simple stellar population , born with a kroupa initial mass function @xcite . as massive stars go sn , energy and metals are deposited into the nearest neighbor gas particles . the sn feedback recipe calculates the radius affected , and turns off cooling in those affected neighboring gas particles until the end of the snowplow phase as described by the sedov - taylor solution @xcite . the amount of energy deposited amongst those neighbors is 0.4@xmath1610@xmath17 ergs , as was adopted in all of our previous work mentioned above . additionally , we include a uniform uv cosmic background following an updated model of @xcite . @xcite demonstrated that when force resolutions @xmath18100pc can be achieved , high density peaks that mimic sf regions in giant molecular clouds can be resolved with several hundreds of gas particles in @xmath810@xmath19 m@xmath7  clumps ( see also * ? ? ? * ; * ? ? ? this allows for the adoption of a realistic density threshold for sf ( 100 amu/@xmath20 , * ? ? ? * ; * ? ? ? * ; * ? ? ? * ) . to match the kennicutt - schmidt law , the high density sf threshold must be offset by a slightly higher value of @xmath15 = 0.1 . this prescription leads to enhanced gas outflows that remove low angular momentum gas from the central regions of the galaxy ( * ? ? ? * brook et al . , in prep . ) . at low masses , this creates a bulgeless disk with a linearly rising rotation curve , comparable to those observed ( oh et al . , in prep ; * ? ? ? * ; * ? ? ? the low mass galaxies used in this paper , drawn from a 25mpc volume , can achieve force resolutions @xmath8100pc . hence , these simulations adopt the more realistic feedback prescription with high density sf and @xmath15 = 0.1 . however , the more massive galaxies presented here are drawn from 50mpc volumes , making it computational expensive to achieve similar force resolutions in a reasonable time . hence , these simulated galaxies can not resolve the high density sf peaks easily , and a lower sf density threshold must be adopted , 0.1 amu/@xmath20 , with a lower @xmath15 = 0.05 . these latter values have been adopted in all of our previous work on mw mass galaxies @xcite , and are a compromise that allow realistic disks to form , but are inefficient at driving gas outflows . thus , low angular momentum material that might be lost from the central regions is maintained , helping lead to the creation of large bulges in the mw mass galaxies ( see also * ? ? ? * ; * ? ? ? * ; * ? ? ? * ; * ? ? ? * ; * ? ? ? * ; * ? ? ? * ; * ? ? ? additionally , the creation of these large bulges can be due to missing physics . in particular , agn feedback has not yet been added to these simulations , and is potentially a key mechanism to help create realistic bulges . due to our small sample size , it is difficult to quantify if our bulges are substantially different from observed bulges in the same galaxy mass range . the simulated galaxies used in this work were selected to span a range of spin values and merger histories , with the last major merger redshift as low at 0.8 . properties of these simulations are listed in table  [ simsum ] . lcccccccccr h516 & 3.9@xmath1610@xmath21 & 2.6@xmath1610@xmath22 & 1.6@xmath1610@xmath23 & 3.3@xmath1610@xmath24 & 0.05 & 0.05 & 1.2 & 87 & 3.5@xmath1610@xmath25 + h799 & 2.2@xmath1610@xmath21 & 1.3@xmath1610@xmath22 & 1.6@xmath1610@xmath23 & 3.3@xmath1610@xmath24 & 0.04 & 0.05 & 3.0 & 87 & 1.9@xmath1610@xmath25 + h603 & 3.8@xmath1610@xmath26 & 3.1@xmath1610@xmath21 & 3.0@xmath1610@xmath27 & 6.3@xmath1610@xmath23 & 0.06 & 0.11 & 1.1 & 231 & 3.8@xmath1610@xmath25 + h986 & 2.1@xmath1610@xmath26 & 2.2@xmath1610@xmath21 & 3.0@xmath1610@xmath27 & 6.3@xmath1610@xmath23 & 0.04 & 0.07 & 0.8 & 231 & 2.4@xmath1610@xmath25 + h239 & 8.5@xmath1610@xmath26 & 7.8@xmath1610@xmath21 & 1.2@xmath1610@xmath25 & 2.1@xmath1610@xmath27 & 0.03 & 0.05 & 1.1 & 347 & 2.8@xmath1610@xmath25 + h258 & 8.0@xmath1610@xmath26 & 7.4@xmath1610@xmath21 & 1.2@xmath1610@xmath25 & 2.1@xmath1610@xmath27 & 0.04 & 0.07 & 0.8 & 347 & 2.8@xmath1610@xmath25 + h277 & 7.1@xmath1610@xmath26 & 6.9@xmath1610@xmath21 & 1.2@xmath1610@xmath25 & 2.1@xmath1610@xmath27 & 0.03 & 0.04 & 3.0 & 347 & 2.3@xmath1610@xmath25 + h285 & 8.7@xmath1610@xmath26 & 8.1@xmath1610@xmath21 & 1.2@xmath1610@xmath25 & 2.1@xmath1610@xmath27 & 0.02 & 0.05 & 1.9 & 347 & 3.0@xmath1610@xmath25 +      all of our simulations are of field galaxies that do not undergo major mergers below @xmath280.5 , making them disk dominated at redshift zero . in order to compare our simulated disk scale lengths to observational results in as realistic a way as possible , we wish to fit the _ light _ profile rather than the underlying mass profile . we create artificial surface brightness images of our simulated galaxies using sunrise @xcite , a monte carlo radiative transfer program that produces a spectral energy distribution ( sed ) for each resolution element of the simulations . this is done by identifying the age and metallicity of each star particle , which are then convolved with the starburst99 stellar population synthesis models @xcite to produce an sed . sunrise assumes dust tracks with the metallicity of the gas particles , and performs ray tracing from each star particle to compute the observed sed including absorption and scattering . from this , we generate mock images in chosen filter bands . figure  [ fig1 ] shows artificial edge - on and face - on @xmath1 band images , including dust reprocessing , for two different mass galaxies . most of the observational surveys that have investigated the evolution of disk sizes worked in the rest frame @xmath1 band across redshifts . to compare to these results as a function of redshift , we have generated face on @xmath29 band surface brightness images at @xmath2=0.5 and @xmath2=1 that include the effects of surface brightness dimming , and using rest frame face on @xmath1 band surface brightness images for the @xmath2=0 galaxies . two of our disk galaxies are undergoing a major , disruptive merger at @xmath2=1 . for these galaxies , we created face on images for the most massive progenitor at a time just prior to the merger ( @xmath2=1.25 ) while the disk was still rather undisturbed . the face on sunrise images were fit with a sersic bulge component and an exponential disk component using the publicly available 2d surface brightness fitting code galfit @xcite . additionally , bars exist in both h603 and h986 at @xmath2=1 that were fit by a third component , reducing the b / d ratio of these galaxies at this step . the galfit results were checked against a 1d , radially averaged fit generated from the same images , and found to be in good agreement . to mimic observations , the fits were required to be a good match down to a limiting surface brightness value of @xmath826 mag / arcsec@xmath30 ( hence the importance of using surface brightness dimmed images at @xmath2=1 ) . example fits are shown in fig . [ sb ] for the same galaxies shown in fig . [ fig1 ] , with @xmath31 band fits shown simply for a comparison of the @xmath1 band results to those at a longer wavelength . we verified that using redshifted , surface brightness dimmed sunrise images in the @xmath29 band at higher z yielded similar disk scale lengths as the rest frame @xmath1 band . as can be seen in table  [ simsum ] , a couple of these galaxies have undergone major mergers since @xmath2=1 , while some galaxies have had a very quiescent history . we note that at @xmath2=1 , two of our galaxies ( h239 and h285 ) had no component that could be well fit by an exponential . this may be due to mergers that occur in these galaxies near this time . galaxy h239 has the most active merger history of our simulated galaxies , with continual mergers ( both major and minor ) until @xmath280.5 . galaxy h285 begins to accrete a satellite with 1/7 of the stellar mass of the main halo at @xmath2=1.25 . however , the core of this satellite does not merge with the main halo until @xmath280.8 . in both cases , we searched in 250 myr periods from 0.75 @xmath32 1.25 to identify a step with a possible exponential disk , but none could be found . hence , those steps have been excluded from the high redshift analysis below . @xcite demonstrated that measuring disk scale lengths based on light rather than a kinematic decomposition could dramatically increase the resulting disk - to - total , d / t , ratio of their simulated galaxies , making measurement technique a potentially important reconciler between observations and simulations . in fig .  [ ml ] we plot the @xmath2=0 surface mass density as a function of radius for both the stars ( solid line ) and the gas ( dashed line ) in the simulated galaxies . also shown as the dotted line is the exponential fit derived in the @xmath29 band for these simulations ( listed in table  [ simsum ] ) . in a future paper ( mccleary et al . , in prep ) , we will explore in detail the differences between mass / kinematic results and light results ( e.g. , scale lengths , bulge - to - disk ratios , the role of dust and inclination , etc . ) . here , however , we note that an initial comparison of mass versus light exponential scale lengths demonstrates a less dramatic difference in our simulations than in @xcite . as seen in fig . [ ml ] , the @xmath29 band scale length is generally a good match to the stellar surface mass density . for a quantitative comparison , scale lengths were derived for the stellar mass over the same region that the disk dominates the light profile . the two methods were found to agree to within 10% for 75% of the galaxies . variation in the techniques appears to be mostly attributable to varying mass - to - light ratios ( mccleary et al . , in prep ) . the bulge / disk decompositions were performed on the dust extinguished sunrise results . however , disk scale lengths and central surface brightnesses were found to be identical in both the face on sunrise images with and without dust reprocessing , while the central bulge component was dust extinguished . thus , using the dust extinguished versus non dust extinguished decomposition has no effect on the results for the disk analysis presented below . it is generally expected that dust should have some effect on both the central surface brightness and the scale length of the disk component , even for face on orientations @xcite . we speculate that our more massive galaxies have little / no effect from dust due to the fact that they are also more gas poor than observed galaxies ( because they use up gas forming too many stars , discussed in more detail in section  [ comp ] ) . our low mass galaxies are gas rich , in agreement with observed galaxies in this mass range @xcite , but metal poor @xcite , making dust negligible . it is important to note that the following analysis purposely focuses on the evolution of the _ disk _ alone . that is , the discussion below of how properties evolve with time neglects the bulge component of the galaxy . this is done for three reasons . first , as discussed above , our mw - mass galaxies tend to have larger bulges than observed . this is a problem that has historically plagued mw - mass simulated disks ( e.g. , * ? ? ? * ; * ? ? ? * ) , and is likely due to the inability to correctly model sf and feedback at the necessary high resolutions to drive loss of low angular momentum gas , and due to missing physics such as agn feedback . hence , while the disks of our simulated galaxies appear to be in good agreement with observed disks ( as presented below ) , the growth of the bulge in these massive galaxies is probably not modeled correctly . second , observations of the evolution of disk galaxies in the magnitude  size plane have either a ) selected galaxies that can be fit with a small pure sersic index , @xmath33 , so that they are disk dominated , or b ) done a bulge / disk decomposition . in either case , the results we wish to compare to are concerned with disks only . finally , our goal here is to examine the growth of _ disks _ using these simulations , not bulges or spheroidal components . `` observable '' properties of the simulated galaxies are listed in table  [ obssum ] . total magnitudes in each band are derived from sunrise , with each galaxy being integrated out to roughly 50 comoving kpc ( i.e. , essentially all of the flux from the simulated galaxy is contained ) . table  [ obssum ] also lists the @xmath2=0 stellar mass derived for the entire galaxy ( disk + spheroid ) based on sunrise colors using the k_sdss_bell routine in the kcorrect package @xcite . a comparison to table  [ simsum ] shows that the simulated stellar masses are generally @xmath840% larger than those derived from the photometric results . a more detailed comparison will be presented in mccleary et al . , ( in prep ) . in table [ obssum ] we list the bulge - to - disk ratios ( b / d ) despite the caveat mentioned above regarding large bulges in our higher mass galaxies . disk scale lengths in the @xmath29 band at @xmath2=0 are given in table  [ obssum ] , while the @xmath1 band scale lengths used in the following plots are given in table  [ evol ] . as seen in observational studies @xcite , the @xmath1 band scale lengths tend to be longer than the @xmath29 band scale lengths , at least for the six most massive galaxies ( with the two dwarf galaxies exhibiting different behavior ) . it has been debated whether this trend for disk scale lengths to be longer at shorter wavelengths is due to dust effects or age gradients @xcite . as noted above , dust plays no role in our face on _ disk _ fits , meaning that this trend in the simulations is due entirely to age and metallicity gradients . lcccccccc h516 & -16.8 & -15.9 & 0.08 & 0.14 & 0.53 & 3.4@xmath1610@xmath22 & 1.0 & 51 + h799 & -16.2 & -15.6 & 0.00 & 0.00 & 0.37 & 1.1@xmath1610@xmath22 & 0.7 & 43 + h603 & -21.3 & -20.5 & 1.25 & 0.72 & 0.52 & 2.2@xmath1610@xmath21 & 4.0 & 143 + h986 & -21.1 & -20.3 & 0.63 & 0.38 & 0.45 & 1.3@xmath1610@xmath21 & 2.8 & 137 + h239 & -22.5 & -21.8 & 0.35 & 0.16 & 0.41 & 4.2@xmath1610@xmath21 & 2.5 & 246 + h258 & -22.4 & -21.5 & 0.86 & 0.45 & 0.50 & 5.3@xmath1610@xmath21 & 5.3 & 204 + h277 & -22.1 & -21.2 & 0.63 & 0.42 & 0.52 & 4.3@xmath1610@xmath21 & 2.3 & 250 + h285 & -22.3 & -21.3 & 1.33 & 0.78 & 0.55 & 5.6@xmath1610@xmath21 & 4.5 & 203 + size , magnitude ( or mass ) , and rotational velocity make up a fundamental plane for disk galaxies @xcite . * see their figure 5 ) demonstrated that these simulated galaxies lie on the observed magnitude  velocity relation for disk galaxies ( also known as the tully - fisher relation ) . in this paper we derive disk scale lengths in order to extend this analysis , and show that our galaxies are also a good match to the observed magnitude  size and velocity  size relations as well . all of our galaxies have sizes and velocities ( table  [ obssum ] ) that agree well with the observed velocity  size trend @xcite . two of the more massive galaxies have slightly high velocities for their size ( h239 and h277 ) , but still lie within the 2@xmath34 scatter of observed galaxies .    the @xmath1 band disk magnitudes , m@xmath35 , and scale lengths , h@xmath35 , of our simulated galaxies at @xmath2=0 are shown in fig . [ fig : z0 ] . the simulated galaxies ( red circles ) are compared to three observational samples , @xcite and @xcite at the massive end , and @xcite at the low mass end . the @xcite data have adopted a dust correction based on @xcite , while the @xcite data are dust corrected based on inclination , following equations 1 , 2 , and 7 listed in @xcite . again , possibly due to low gas ( and therefore dust ) content , the six more massive galaxies shown here are free of dust effects . as such , they should be compared to dust corrected observational data . at the low mass end , however , @xcite makes no internal dust corrections , as the low metallicities of the dwarf galaxies in her sample are expected to lead to little effect from dust . clearly , the disk sizes of the simulated galaxies are in good agreement with observed disk sizes . we use the @xmath1 band scale lengths of the observational samples in fig . [ fig : z0 ] , which can be compared to the @xmath29 band results at @xmath2=1 in the next section since the @xmath1 band is redshifted into the @xmath29 band at @xmath2=1 . however , the magnitude  size ( and stellar mass  size ) relation has been measured in multiple bands for a large sample of sdss galaxies @xcite . the sdss results are consistent with the results plotted in fig . [ fig : z0 ] , though we purposely select observational samples that have done a bulge / disk decomposition and fit a pure exponential to the disk surface brightness fit , in order to isolate the evolution of the disk and eliminate any contamination from a central spheroid . as discussed in the introduction , historically it has been a challenge for cosmological disk galaxy simulations to match the observed fundamental relations for galaxy disks ( e.g. , * ? ? ? the failure to reproduce observed trends is the result of the angular momentum catastrophe , exacerbated by the fact that previous studies followed mass / kinematic results rather than light profiles . our success at matching the observed disk scaling relations is a result of the fact that our disk baryons maintain angular momentum through time , due to an increased resolution and a physically motivated feedback mechanism ( see * ? ? ? * for a review of how each of these processes contribute to the formation of realistic disks ) . a number of works have attempted to isolate the separate role of resolution and feedback ( e.g. , * ? ? ? * ; * ? ? ? * ; * ? ? ? * ; * ? ? ? * ) . in n - body + sph simulations , dark matter particles are typically an order of magnitude more massive than the gas particles , leading to an exchange of kinetic energy in two body interactions that kinematically heats the disk , randomizing velocity vectors , and potentially turning a disk component into a spheroid @xcite . this effect can be reduced at higher mass resolutions that lower the discrepancy in particle masses . @xcite used controlled ( non - cosmological ) n - body+sph simulations to evaluate the amount of angular momentum that can be lost purely due to resolution . they concluded that 10@xmath19 particles within the virial radius are necessary for a disk galaxy to maintain roughly 90% of its original angular momentum .    yet even if resolution can be increased to the point of maximum angular momentum retention , disks will still suffer dramatic angular momentum loss compared to observed disk galaxies if energy feedback is neglected . feedback prevents rapid , early cooling of gas particles ( the overcooling problem discussed in the introduction ; see references therein ) . heating and expansion of the gas creates a hot reservoir , allowing the gas to cool at later times after the era of rapid mergers , and preventing angular momentum loss via dynamical friction in mergers ( e.g. , * ? ? ? additionally , by preventing overcooling , feedback prevents gas from rapidly turning into collisionless star particles . without feedback , halos of all masses are equally efficient at converting gas into stars @xcite , producing galaxies that are too gas poor compared to @xmath2 = 0 disk galaxies . because feedback regulates star formation , it allows for gas reservoirs to develop that allow disks to survive to the present day @xcite . the existence of large , thin disks at the present day thus requires feedback . not all sn feedback schemes lead to disk galaxies that satisfy observed constraints such as the tully - fisher relation , size  velocity or size  luminosity relations , or the stellar mass  metallicity relation . because the nearest neighbor gas particles surrounding sne are dense and cold , a simple energy deposition will quickly be radiated away and not affect the simulation ( e.g. , * ? ? ? * ; * ? ? ? two main schemes have been adopted to overcome this problem . in the first , a multiphase model of the ism is implemented with a hot gas reservoir inside each gas particle @xcite , overcoming the problem of thermal energy being spread over the entire particle mass , and preventing the hot gas particles from being artificially influenced by their cold gas nearest neighbors . in the second , cooling is turned off in the gas particles near a sn explosion in order to mimic the sub - resolution adiabatic expansion of the sn @xcite . the `` blastwave '' scheme adopted here also turns off cooling in nearest neighbor particles , but attempts to model this based on what is known about actual sne , determining the radius of each sn remnant based on the analytic blastwave solution for a sn remnant @xcite , and cooling is only turned off for those particles within the blast radius . because many sne typically contribute feedback within a dense star forming region , the thermal energy from all of these sne can combine to create a larger blast radius . the differences in the resulting disk between a simple energy deposition ( `` thermal '' feedback ) and the blastwave model have been examined in @xcite and @xcite . as examined by these previous works , the adoption of the `` blastwave '' feedback model , combined with high resolution ( all simulations presented here have more than 10@xmath19 particles within the virial radius at @xmath2 = 0 , see table  [ simsum ] ) , overcome past problems with cosmological disk galaxy simulations and allow for the present study . having established that these simulated disk galaxies match observed disk scaling relations at @xmath2=0 , we now use them to investigate the degeneracy in the magnitude  size evolution with time . first , we establish that the simulations also match the available observational data at @xmath2=1 . the evolution of the magnitude  size relation for these simulated galaxies is shown in fig . [ fig : magsize ] . the top panel shows the simulation results at @xmath2=0 , 0.5 , and 1 . it is evident that the magnitude  size relation for these galaxies is evolving in time , dimming in surface brightness since @xmath2=1 . despite the fact that these are individual galaxies being followed in time , the snapshot of the population of disks at each redshift is consistent with the observations . this is evident in the bottom two panels . the middle panel compares the simulated disk scale lengths at @xmath2=1 to the decomposed @xmath29 band disk scale lengths of galaxies at 0.9 @xmath36 1.2 from @xcite and miller et al.(_in prep _ ) . no dust corrections are applied to the observational @xmath2=1 data , as these galaxies are expected to be low metallicity where corrections are uncertain and often ignored @xcite . any realistic dust corrections are likely to be tiny , and not affect the results . the bottom panel is a reproduction of fig . [ fig : z0 ] , shown again for easy comparison of the results across redshifts . while observational biases allow for a direct comparison only at the massive end at @xmath2=1 , the disk sizes of the simulated galaxies are in good agreement with observed disk sizes , as a function of redshift . we note that observational data that used a single sersic component fit ( e.g. , * ? ? ? * ) are also in good agreement with the @xmath2=1 data presented here ( and lie in the same magnitude range ) , despite the lack of a bulge / disk decomposition . for clarity , we show only the decomposed data in fig . [ fig : magsize ] . the top panel of fig . [ fig : magsize ] shows that there has been apparent dimming of the surface brightness of the simulated population since @xmath2=1 . [ deltamu ] quantifies this evolution , showing the change in the surface brightness , @xmath37 , of the simulated galaxy disks . we have defined @xmath37 such that @xmath38 data points in fig . [ deltamu ] at @xmath2=0 , 0.5 and 1 are for the mean @xmath37 for the population at that time , and error bars reflect the standard deviation of the sample . because our two dwarf galaxies are too faint to be generally observable at @xmath2=1 , fig . [ deltamu ] shows the evolution of these two galaxies separately from the higher mass galaxies . however , the evolution is roughly similar in both mass ranges , with about 1.5 magnitudes in surface brightness dimming between @xmath2=1 and @xmath2=0 . a dimming of 1.5 magnitudes is consistent with observational studies that have examined the surface brightness evolution of disk dominated galaxies ( e.g. , * ? ? ? * ; * ? ? ? * ; * ? ? ? * ; * ? ? ? * ; * ? ? ? * ; * ? ? ? lcccccccc h516 & 22.7 & 1.5 & -15.7 & 2.0 & 0.95 & 1.0 & 4.3@xmath1610@xmath39 & 1.69 + h799 & 22.1 & 0.7 & -15.6 & 0.0 & 0.6 & 1.3 & 1.8@xmath1610@xmath39 & 1.32 + h603 & 22.5 & 2.8 & -19.8 & 1.0 & 4.2 & 2.4 & 1.3@xmath1610@xmath21 & 2.78 + h986 & 21.2 & 1.0 & -19.9 & 0.2 & 3.0 & 1.9 & 1.0@xmath1610@xmath21 & 3.22 + h239 & 19.3 & n / a & -21.5 & n / a & 3.2 & n / a & 9.0@xmath1610@xmath6 & 10.17 + h258 & 21.3 & 1.6 & -21.1 & -0.2 & 6.3 & 2.7 & 3.2@xmath1610@xmath21 & 2.53 + h277 & 19.7 & 0.6 & -20.8 & 0.4 & 2.4 & 1.4 & 1.9@xmath1610@xmath21 & 1.95 + h285 & 21.4 & n / a & -20.6 & n / a & 4.6 & n / a & 2.0@xmath1610@xmath21 & 3.07 +    the previous section examined the population of simulated disks in terms of surface brightness evolution . we now decompose the surface brightness trends into luminosity and size evolution , and demonstrate how individual galaxies evolve in each property with time . table  [ evol ] quantifies the changes in the disk @xmath1 band central surface brightness ( @xmath40 , as opposed to @xmath41 plotted in fig . [ deltamu ] ) , size , and magnitude for each of our simulated galaxies . it is widely expected that galaxies should undergo dimming since @xmath2=1 , due to declining star formation rates . the key question is how much galaxies dim . it is already apparent from fig . [ fig : magsize ] that our more massive galaxies are undergoing growth in size between @xmath2=1 and @xmath2=0 , meaning that the change in surface brightness of @xmath81.5 magnitudes over this time can not be entirely due to luminosity dimming at a fixed size . the changes in size and magnitude that are quantified in table  [ evol ] are also visualized in fig . [ dmagsize ] . now it can be seen that there is in fact a dependency on mass in the evolution of the galaxies . our more massive galaxies are growing in size and generally undergoing only a small amount of dimming in total magnitude ; some are even getting brighter with time . meanwhile , there is little change in the scale length of our lowest mass galaxies , but their evolution in luminosity is erratic ( for this sample of two galaxies ) due to their bursting star formation histories .    the luminosity evolution of the low mass dwarf galaxies is tied up with their particular star formation histories . while h799 is quiescent between @xmath2=1 and @xmath2=0.5 so that it dims , it undergoes another burst of sf at low @xmath2 , increasing its magnitude . on the other hand , h516 stays relatively quiescent between @xmath2=1 and @xmath2=0 , so that it undergoes significant luminosity dimming ( of 2.0 magnitudes ) . this large amount of dimming may at first seem to be at odds with the fact that low mass galaxies are generally undergoing more sf today than they were in the past ( this is discussed further in section  [ starmass ] , and is one form of galaxy `` downsizing '' ) . however , this extreme dimming for low mass galaxies has been seen in observational results as well @xcite . @xcite concluded that small galaxies ( those with half light radii @xmath18 2 kpc ) have undergone roughly 2.5 magnitudes of dimming since @xmath2=1 . given the bursty sf nature of the two low mass dwarfs we have in our sample , it is not possible to compare statistically to the observations , but there is nothing to suggest that our simulated dwarfs are inconsistent with the @xcite results .    at the higher mass end , @xcite find @xmath81.5 mag of dimming since @xmath2=1 for galaxies with half light radii @xmath42 3 kpc ( disk scale lengths @xmath42 1.8 kpc ) , if the strictest interpretation is made that all surface brightness evolution is in luminosity rather than size . clearly , this is not true for our high mass galaxies , which are growing with time . if galaxies are growing with time , then comparing galaxies that are a given size ( e.g. , 3 kpc ) at @xmath2=0 with galaxies that are of similar size at @xmath2=1 means that similar galaxies are _ not _ being compared . those @xmath2=0 galaxies will have been smaller at @xmath2=1 , and it is unknown from our simulation results what size galaxies with scale lengths larger than @xmath82.5 kpc at @xmath2=1 will have evolved to by @xmath2=0 . in summary , if galaxies are growing since @xmath2=1 , this will mimic a larger change in luminosity on the magnitude  size relation than actually occurs @xcite . our simulated disk galaxies _ as a population _ appear to match the observational results for the evolution in the size - magnitude plane very well ( fig . [ fig : magsize ] ) . yet [ dmagsize ] demonstrates what the observations can not witness : individual galaxies are growing with time , so that the evolution in the magnitude  size plane can not be due to dimming alone . disk scale length evolution for each galaxy is plotted in fig . [ evolsize ] , normalized to their size at @xmath2=1 . these values are also quantified in table  [ evol ] . only a few of these galaxies show negligible size evolution , with the two lowest mass galaxies undergoing the smallest size changes . the remaining galaxies show clear growth , and about half evolve by the amount predicted by the simple sis model , ( nearly a factor of two back to @xmath2=1 ) , or more . the dotted line in fig . [ evolsize ] shows the growth predicted for the simple sis model @xcite , while the dashed line shows the growth determined by @xcite using 16 galaxies at @xmath2=1 . however , it is not necessarily expected that any given galaxy should follow the @xmath9 analytic growth predicted by the sis model . @xcite demonstrated that more sophisticated models adopting nfw profiles could bring the observations and theory closer in line . the prediction that sizes scale with @xmath9 is based on the assumption that @xmath43 where @xmath10 is the circular velocity of the halo and @xmath44 is the radius at which the mean density is equal to 200 times the critical density . then @xmath45 where @xmath46 is the dimensionless spin parameter @xcite . if the disk size follows @xmath9 , both @xmath46 and @xmath10 must stay roughly constant with time . while @xmath10 does stay roughly constant back to @xmath2=1 in these galaxies , @xmath46 does not . it is known that @xmath46 may vary for individual halos due to mergers or smooth accretion @xcite . the previous work of analytic models do not require @xmath46 to stay constant , as they simply use a snapshot in time and the spin values of halos at that snapshot . that is , they study a population of galaxies at a given time , and make no attempt to follow an individual galaxy with time as we do here . hence , it is not surprising that some of the simulated galaxies presented here deviate from growing precisely as @xmath9 .    a better question , if we want to compare our disk sizes to analytic models , is to ask if the scale lengths measured based on our light profiles are similar to those predicted by analytic models . for an nfw density distribution , the disk scale length , r@xmath47 is @xmath48 where @xmath49 measures the halo concentration , and @xmath50 is nfw the scale radius . @xmath51 uses @xmath52 to determine the energy of the nfw halo compared to the sis model , and @xmath53 is a factor that accounts for the gravitational effect of the disk ( adiabatic contraction ) . @xmath54 and @xmath55 are the fraction of disk angular momentum and mass , respectively , to the total halo . for sizes to scale as @xmath9 , it is assumed that @xmath56 = 1 . this is likely to not be true , though , as demonstrated by @xcite . we have determined the values of @xmath51 , @xmath57 , @xmath58 , and @xmath59 for each of our galaxy halos at @xmath2=0 , directly from the simulation results . we compared the @xmath60 expected from these quantities , as derived by equation [ eq:4 ] , to the @xmath29 band scale lengths listed in table  [ obssum ] . the results agree to within 10% for only the two lowest mass galaxies . for the more massive galaxies , the predicted nfw @xmath60 is consistently shorter than that found for the @xmath29 band , and for four of the six massive galaxies , the measured @xmath29 band results are a factor of 2 - 4 larger than the nfw prediction . the main source of this discrepancy appears to be in the value obtained for the spin value , @xmath46 , for the entire halo . if the spin of the gas , @xmath61 , within the halo is used instead , the discrepancy shrinks so that six of the eight galaxies have nfw @xmath60 and @xmath29 band results within 10% of each other . as listed in table  [ simsum ] , the values of @xmath46 are generally lower than @xmath61 for all halos except the two lowest mass galaxies . this trend for @xmath61 to be larger than @xmath46 has been seen previously in simulations @xcite , and has been suggested as one reason why disks may have a smaller fraction of low angular momentum than their dark matter halos @xcite . recently , @xcite used a similar simulation run with gasoline to show that gas that enters the virial radius of a halo and cools toward the disk is torqued by the hot halo gas , so that the angular momentum of the disk gas becomes aligned with the hot halo . this is true even for gas that is initially counter - rotating with the disk , and hence is a powerful method to remove negative and low angular momentum material from the disk . the lack of low angular momentum material for the disk gas will lead to larger values of @xmath61 compared to @xmath46 for the total halo . the disk stars , whose light we trace in the results of this paper , form from this cold gas with larger spin values , and thus @xmath61 will be a better predictor for the nfw @xmath60 than @xmath46 . the result that the disk baryons are lacking a low angular momentum component compared to the dm is not trivial . this result needs to be examined in detail . however , this requires a full evaluation of the history of angular momentum in these halos , which is beyond the scope of the current paper . we reserve such a study for future work . as discussed above , numerous observational studies have found that there has been a decline in the surface brightness of disk galaxies , by a magnitude or more , since @xmath2=1 . due to the fact that there is little evidence for a change in the size function of disk galaxies over this time period ( e.g. , * ? ? ? * ) , the observations can only be interpreted as an upper limit to luminosity evolution . that is , if disks do not change in size , then all of the evolution must be in luminosity . however , due to the degeneracy in the magnitude  size plane , the option of possible growth is left open . @xcite investigated this disk growth . after showing that they , too , agreed with @xmath81 magnitude of surface brightness dimming back to @xmath2=1 , they then used galaxy colors to derive stellar masses and considered the stellar mass  size relation for their galaxies . they found that this relation , unlike the magnitude  size relation , showed little or no evolution back to to @xmath2=1 . disk growth may still occur , but must occur in such a way that galaxies remain on the same stellar mass  size relation with time . they concluded that this was evidence for weak inside - out growth of galaxy disks , as did additional later studies @xcite . the right panel of fig . [ dmagsize ] shows the evolution in the stellar mass  size plane for our simulated galaxies . in agreement with previous results @xcite , these galaxies are growing along a path that keeps them on the same stellar mass  size relation with redshift .    the stellar mass shown in the right panel of fig . [ dmagsize ] is for the disk only . to separate the disk growth from total stellar growth of these galaxies ( i.e. , from the bulge and halo spheroidal stellar components ) , a kinematic decomposition was done to identify disk stars at @xmath2=0 . to identify disk particles , the galaxies are first aligned so that the disk angular momentum vector lies along the _ z_-axis . j@xmath62/j@xmath63 is calculated for each star particle in the galaxy , where j@xmath62 is the angular momentum in the _ x - y _ plane , and j@xmath63 is the momentum that a particle would have in a circular orbit with the same orbital energy . disk stars are identified as those having near circular orbits , so that j@xmath62/j@xmath63 @xmath42 0.8 . 0.2 , which matches the eccentricities observed in the milky way disk @xcite . ] these @xmath2=0 disk stars were then searched for in their most massive progenitor at @xmath2=0.5 and @xmath2=1 to find the mass of these stars that had formed at each step , with the results plotted in fig . [ dmagsize ] . if galaxy disks are evolving approximately along the same stellar mass  size relation with time , then the surface mass density of disks should show little evolution back to @xmath2=1 @xcite . that is , a galaxy with stellar mass of 10@xmath21 m@xmath7  at @xmath2=1 is roughly the same size ( though maybe just slightly smaller ) as a galaxy at @xmath2=0 with 10@xmath21 m@xmath7 , yielding roughly the same surface densities . [ barden ] tentatively confirms this result , showing that the surface mass density of our population of disks evolves little with time , with surface mass density , @xmath64 defined as @xmath65 the red points with error bars , connected by the solid black line , show the mean and standard deviation for the whole sample at @xmath2 = 0 , 0.5 , and 1 . the colored lines are the results for individual simulated disks with time . we note , however , that the dashed line shows the prediction for the surface mass density evolution for the simple sis model @xcite . the standard deviation within our small sample is large , so that the simulated galaxies are fully consistent with the predicted growth . while our sample is clearly too small to derive statistical results , the change in the mean surface density of the population between @xmath2=1 and @xmath2=0 is only 0.1 dex , half of the evolution predicted by the sis model . we note that an updated model by @xcite that incorporates the evolution of nfw halos derived in n - body simulations predicts only @xmath80.2 dex change in the stellar surface mass density back to @xmath2=1 . this small evolution in surface mass density is due to evolution along a stellar mass  size relation that changes little with time .      due to the fact that the stellar mass  size relation shows little change with time , examining the evolution of galaxy disks at a fixed _ stellar mass _ is similar to examining the evolution at a fixed _ size_. however , because galaxies are growing with time , a galaxy that is , e.g. , 3 kpc in size and 10@xmath21 m@xmath7  in stellar mass at @xmath2=1 has increased in mass and size by @xmath2=0 and another galaxy has moved from a lower mass and smaller size into the same bin . observations , and the results presented here , tell us that these two galaxies of same size at @xmath2=1 and @xmath2=0 have a different luminosity : a galaxy at @xmath2=1 with similar stellar mass and size to a @xmath2=0 galaxy will be brighter . note that this is _ not _ equivalent to saying that an individual galaxy was brighter in the past . instead , galaxies at fixed stellar mass were brighter in the past . because they are the same stellar mass , the brighter @xmath2=1 galaxy must then be forming more stars . two of our galaxies can be compared in this way . galaxy h603 has a stellar disk mass of 1.28@xmath1610@xmath21 m@xmath7  at @xmath2=0 , while h258 has a stellar disk mass of 1.25@xmath1610@xmath21 m@xmath7  at @xmath2=1 . the sfr within the last 100myr for h258 at @xmath2=1 is 3.1 m@xmath7/yr , while it is 0.6 m@xmath7/yr for h603 at @xmath2=0 . the lower sfr for h603 is due to the fact that it resides in a less massive halo than h258 ( see table  [ simsum ] ) . even at @xmath2=1 , h603 had a lower sfr than h258 , of 2.1 m@xmath7/yr . this is because these galaxies follow a sfr  stellar mass relation ( which also leads to their reproduction of the observed stellar mass  metallicity relation , * ? ? ? [ fig9 ] shows the sfr  stellar mass and specific sfr ( ssfr )  stellar mass relations for our simulated disks . the left panels of fig . [ fig9 ] shows the ssfr ( top ) and sfr within the last 100 myr ( bottom ) versus stellar disk mass . these left panels confirm the values listed in the paragraph above . the dwarf galaxies are included in this plot for full disclosure , but lie outside of the mass range observable at high @xmath2 , unlike the higher mass galaxies . two trends among the higher mass galaxies are noteworthy . first , the sfr of any individual galaxy tends to decrease slightly with time , but second , the entire population as a whole can be observed to shift to lower sfrs and ssfrs @xcite . the decline in sfr and ssfr in has been observed extensively in populations of galaxies @xcite , and is one form of `` downsizing '' in galaxy evolution . the right panels of fig . [ fig9 ] also show sfr and ssfr for the simulated galaxy disks , but now as a function of total halo mass . this is included to demonstrate that although , e.g. , h603 at @xmath2=0 and h258 at @xmath2=1 have similar disk stellar masses , this does not correspond to similar virial masses for those redshifts ( colors yellow and purple in the figures , respectively ) . the deeper potential well of the more massive galaxy ( h258 ) leads to a higher sfr.=1.25 due to the fact that it is undergoing a merger at @xmath2=1 . galaxy h258 is 50% more massive than h603 at @xmath2=1 . ]    in conclusion , attempting to compare galaxies that are a similar stellar mass at @xmath2=1 and @xmath2=0 ( and by the tight stellar mass  size relation , similar size ) leads to inherently comparing galaxies of two different virial masses . this then leads to the observed difference in luminosity for galaxies of a similar size at @xmath2=1 and @xmath2=0 , caused by the sfr  stellar mass relation , and the corresponding stellar mass  halo mass relation @xcite .      in fig . [ fig10 ] we reproduce some of the simulated scaling relations , this time with observed scaling relations directly overplotted for comparison . the left panel of fig . [ fig10 ] shows the low redshift relations for sdss galaxies ( solid line , with 2@xmath34 scatter in dotted lines ) from @xcite . it is worth noting that @xcite performed bulge / disk decompositions with the sdss data , so that the sizes plotted are fits to purely exponential disk scale lengths ( converted from the half light radius presented in their paper ) . @xcite also use higher @xmath2 data from the deep2 survey to measure the change in size since @xmath2=0 , finding @xmath80.1 dex of size evolution at a fixed mass back to @xmath2=1 . in the mass range of overlap ( m@xmath4 @xmath42 10@xmath6 m@xmath7 ) , our simulated galaxies are entirely consistent with observed galaxies within the scatter . the central panel of fig . [ fig10 ] reproduces the stellar mass  sfr relation , now with observational data from @xcite plotted for @xmath2=0 ( solid line ) and @xmath2=1 ( dashed line ) , corrected from a salpeter imf to a kroupa imf . at both redshifts , the minimum mass observed is roughly 10@xmath6 m@xmath7 . the right panel shows the same data , but now for the ssfr . the data point with error bars in the lower right corner of the right panel provides an estimate of the typical scatter in the @xcite data . while our galaxies are certainly consistent with the data within the scatter , and the general slope at each @xmath2 is consistent with the observations , there is overall a tendency for the simulated galaxies to be more massive than the observations . we believe the discrepancy is in stellar mass rather than sfr , because these simulated galaxies have both @xmath2=1 and @xmath2=0 disk sfrs in good agreement with observed galaxies in the same halo mass range @xcite . we note that the @xcite data is for total stellar mass , while our simulated results are for the disk only ( bulge excluded ) . however , including the bulge mass will only exacerbate the discrepancy . for the massive galaxies , the discrepancy is probably due to the fact that these simulations are known to form too many stars @xcite . we can learn something about _ when _ this overproduction of stars occurs in the simulations by also examining the stellar mass  ssfr plot ( right panel of fig . [ fig10 ] ) . while within the observed scatter , is appears that the ssfrs of the simulated disks are consistently lower than observational results at both @xmath2=1 and @xmath2=0 @xcite . these plots demonstrate that these simulated galaxies must be forming too many stars at early times , at redshifts higher than plotted here ( @xmath66 1 ) . because the sfrs of the galaxies are in good agreement with observed galaxies in the same halo mass range at both @xmath2=1 and @xmath2=0 @xcite , the fact that the ssfrs are lower than observed values requires that the stellar mass of these disks be too large by @xmath2=1 . the story is slightly different for the two dwarf disk simulations . these two dwarf galaxies have been published previously in @xcite , where it was demonstrated that their properties are a good match to observed dwarf disk galaxies , including stellar mass @xcite . however , the observed ssfrs of low mass galaxies is on average higher than that of massive galaxies @xcite , indicating that low mass galaxies have been very inefficient in the past at turning gas into stars . clearly , that is not the case with our two simulated dwarf disk galaxies , which have too low ssfrs compared to most observed dwarf galaxies due to the fact that they form the bulk of their stars prior to @xmath2=1 . hence , while approximately the correct mass in stars is formed , the star formation is biased toward @xmath66 1 . we have verified that these two simulations , when run with stronger sn feedback , have burstier sf histories with their sf spread out more evenly down to the present time . this leads to approximately the same mass in stars formed by @xmath2=0 , but with a more uniform sf history . this means that the stronger feedback dwarf runs have higher ssfrs than currently presented here , and ssfr values larger than the massive galaxy simulations shown here . future work will include these new dwarf disks , which show similar stellar masses , rotational velocities , and sizes to those presented in this paper . thus , in all of the simulated galaxies the star formation history is biased toward too much star formation at high @xmath2 . by @xmath2=0 , though , the dwarf galaxies have formed approximately the observed amount of stars for their halo mass , while the more massive galaxies have overproduced the total amount of stars . it is suspected that this overproduction of stars in the massive galaxy disks may be the result of the adopted low sf density threshold , which allows for most of the disk to be capable of sf at any given time , unlike the high density threshold which limits sf only to rare density peaks . however , as discussed in section 2.1 , it is currently computationally expensive to generate the more massive galaxies with the force resolution required to resolve high density peaks where sf occurs , and no simulations of sufficiently high resolution are currently available to test this hypothesis . the adopted sf threshold is a compromise that allows stars to form at radii observed in real galaxy disks , and we have demonstrated that it leads to scale lengths in good agreement with observed galaxies . however , we note that even if our massive galaxies are forming too many stars , it is by roughly a factor of 3 at most , and at early times . this means that the more massive galaxies may shift up to 0.5 dex to lower stellar masses in the stellar mass  size plane , but will still grow along the stellar mass  size relation , so that our conclusions are unaffected . likewise , the @xmath1 band magnitudes we present here are unlikely to be affected , as the bulk of overproduced stars must occur at @xmath66 1 . since @xmath1 band light is dominated by young stars , and our low @xmath2 sfrs are in good agreement with observed galaxies in the same halo mass range , the magnitude  size relations presented in figs . [ fig : z0 ] and [ fig : magsize ] are unlikely to change . hence , the conclusions presented here about the origin of the evolution in the size  luminosity relation remain valid . for the first time , a sample of very high resolution galaxy disks simulated within a fully cosmological context have been used to determine the evolution of the galaxy disk size  magnitude relation . unlike previous simulations , we determine disk scale lengths and magnitudes by fitting the stellar light profile rather than decomposing our simulations based on kinematics ( and hence mass ) . we generate artificial surface brightness images using sunrise @xcite , in order to derive disk properties using methods used by observers . exponential disk scale lengths were derived using galfit @xcite with redshifted , surface brightness dimmed @xmath29 band images at @xmath2=1 and 0.5 , and rest frame @xmath1 band images at @xmath2=0 . this allows for a direct comparison between simulated and observational results in the size  magnitude plane . [ fig : z0 ] demonstrates that these galaxies overcome past inabilities for cosmological disk simulations to retain angular momentum . artificial loss of angular momentum is minimized due to high resolution and a physically motivated sn feedback recipe , and hence we produce disks with realistic sizes compared to observed disks . fig . [ fig : magsize ] shows that the sizes are in agreement with observed disk sizes back to @xmath2=1 . this success , and our previous successes in matching the observed stellar mass  metallicity relation with time @xcite , and the tully - fisher relationship @xcite , indicate that these simulated disks reproduce the fundamental scaling relations for disk galaxies . having established this , we can now use these simulated disks to investigate and interpret observational findings . we use the simulated disk properties to investigate the evolution in the magnitude  size relation with time . observations have generally concluded that there has been 1 to 1.5 magnitudes of surface brightness dimming since @xmath2=1 @xcite . we find similar evolution for our galaxies ( see fig . [ deltamu ] ) . combined with a lack of evidence for an evolving size function for galaxy disks @xcite , and little to no growth in galaxy disks sizes back to @xmath2=1 @xcite , the observed surface brightness evolution can only be interpreted as an upper limit to luminosity evolution if there has been no change in galaxy disk sizes . we have shown that we 1 ) are in good agreement with the observed magnitude  size relation with time and the evolution in surface brightness , but 2 ) a number of the simulated disks are clearly undergoing a large change in size ( figs . [ dmagsize ] and [ evolsize ] ) . our ability to follow individual galaxies with time , unlike the observations , allows us to interpret the evolution in the magnitude  size plane . we find that the evolution is dependent on mass , with our massive galaxies undergoing larger changes in size than magnitude . our dwarf galaxies undergo the least change in size , though there is no immediate theoretical expectation that this should be the case . recent results @xcite show that the rate of growth of halos has little dependence on mass , and thus low mass galaxies should grow just as much in size as their more massive counterparts . we conclude that a larger sample is needed to study if low mass galaxies truly grow less than more massive galaxies . on the other hand , the evolution in magnitude for dwarf disk galaxies is dependent on their bursty sf history , but an individual dwarf can undergo significant dimming of at least 2 magnitudes . a dichotomy of evolution with mass has been seen in observations , with low mass galaxies exhibiting significantly more dimming since @xmath2=1 than more massive galaxies ( e.g. , * ? ? ? * ; * ? ? ? the halo properties for these galaxies ( @xmath67 ) predict a scale length that is generally shorter than the scale lengths derived from the simulated light profiles . fig . [ ml ] demonstrates that the light and mass distributions yield scale lengths in good agreement , so the discrepancy between @xmath60 predicted for nfw halos and from our light profiles can not be due to a discrepancy between light and mass profiles . rather , using @xmath61 , the spin of the gas in the halo , brings the predicted and measured scale lengths into decent agreement . thus , the angular momentum distribution of the gas from which the stars are forming has been modified from that of the dark matter ( see also * ? ? ? * ; * ? ? ? @xcite have demonstrated that accreted gas is torqued by the hot halo after entering the virial radius , potentially preventing counter - rotating and low angular momentum gas from reaching the disk . a full examination of how this process occurs is beyond the scope of this paper and left for future work . having determined that these galaxies are growing in size with time , we can investigate this growth in stellar mass . like previous theoretical models @xcite , we find that these galaxies grow in such a way as to stay along approximately the same stellar mass  size relation with time . hence , at a fixed stellar mass , there is little change back to @xmath2=1 in the size of galaxies .    due to the fact that our larger galaxies have grown since @xmath2=1 , we do not have the ability to directly compare results at a fixed size . in particular , we do not have a population of simulated disks at @xmath2=1 with scale length greater than 3 kpc , though these larger disks do exist observationally . presumably we might generate these larger disks if we simulate higher mass halos . however , it is not clear from this study how these larger disks would evolve to @xmath2=0 . observationally , the number density of disk galaxies shows no evidence for change back to @xmath2=1 @xcite . in conjunction with the apparent lack of evolution in the size function of disks , this fact has been interpreted to mean that either galaxies are done growing by @xmath2=1 , or that the rate of destruction of galaxies at given size must be equally matched by the rate of growth of galaxies into that size . the fact that our galaxies are growing , while still matching the size - magnitude evolution , suggests that galaxies are moving into and out of a given bin in size during this time interval . it is beyond the scope of this work to answer what happens to these larger galaxies , though we may speculate . there are two possibilities that arise due to mergers ; perhaps the larger disks undergo mergers and become early type galaxies @xcite , or perhaps mergers cause them to become bulge dominated disks . as the majority of size - magnitude evolution studies have focused on pure disks or disks with small bulges , this would remove them from the sample being considered @xcite . finally , we demonstrate that the simulated galaxy disks follow sfr  mass and ssfr  mass relations similar to observations . while the growth of dark matter halos is nearly self - similar with mass , observations demonstrate that some process must break this self - similarity in the baryons ( e.g. , * ? ? ? * ; * ? ? ? * ; * ? ? ? * ) . in the mass range of galaxy halos presented here , sn feedback is the process that regulates sf as a function of halo mass , and likely drives gas outflows that vary as a function of mass ( though outflow rates are left for future work ) . the regulation of sf as a function of mass leads to the reproduction of the observed stellar mass  metallicity relationship @xcite , and the sfr  mass relation shown here that is responsible for the reproduction of the size  luminosity relation back to @xmath2=1 . a galaxy at a fixed m@xmath4 at @xmath2=1 will have a higher sfr than a @xmath2=0 counterpart , and thus a higher luminosity as well . in conjunction with the weak evolution of the stellar mass  size relation , this tells us that the difference in luminosity observed at a fixed size between @xmath2=1 and @xmath2=0 is due to the fact that the @xmath2=1 galaxy has a higher sfr and luminosity . we would like to thank s. miller and r. ellis for use of their @xmath2=1 disk scale length data . ab acknowledges support from the sherman fairchild foundation . as was supported by the gordon & betty moore foundation . fg and tq were funded by nsf ast-0908499 . fg acknowledges support from a theodore dunham grant , hst go-1125 , nsf grant ast-0607819 and nasa atp nnx08ag84 g . pj acknowledges support from the w. m. keck foundation . simulations were run at tacc , arsc , and nas .
Abstract:

Answer: a sample of very high resolution cosmological disk galaxy simulations is used to investigate the evolution of galaxy disk sizes back to redshift 1 within the @xmath0cdm cosmology . artificial images in the rest frame @xmath1 band are generated , allowing for a measurement of disk scale lengths using surface brightness profiles as observations would , and avoiding any assumption that light must follow mass as previous models have assumed . we demonstrate that these simulated disks are an excellent match to the observed magnitude size relation for both local disks , and for disks at @xmath2=1 in the magnitude / mass range of overlap . we disentangle the evolution seen in the population as a whole from the evolution of individual disk galaxies . in agreement with observations , our simulated disks undergo roughly 1.5 magnitudes / arcsec@xmath3 of surface brightness dimming since @xmath2=1 . we find evidence that evolution in the magnitude  size plane varies by mass , such that galaxies with m@xmath4 @xmath5 10@xmath6 m@xmath7  undergo more evolution in size than luminosity , while dwarf galaxies tend to evolve potentially more in luminosity . the disks grow in such a way as to stay on roughly the same stellar mass  size relation with time . finally , due to an evolving stellar mass  sfr relation , a galaxy at a given stellar mass ( or size ) at @xmath2=1 will reside in a more massive halo and have a higher sfr , and thus a higher luminosity , than a counterpart of the same stellar mass at @xmath2=0 .