given the diversity of mass and length scales , redshifts , initial conditions and physical processes that take part in the hierarchical formation of galaxies , the fact that most galaxies exhibit some common key properties , such as an exponential disk profile e.g. @xcite , and follow common scaling relations such as star formation rate - stellar mass relation ( e.g. @xcite ) is remarkable .
such simple relations may imply that despite the diversity of initial and end states , there is a  subset " of universal processes that effectively dominate these relations .
even though many constituents are necessary to model correctly galaxy formation to its last detail , a toy model that involves as few parameters and processes as possible , and can nevertheless yield those key properties , is powerful in its ability to offer intuition towards their understanding .
this approach has motivated the _ ` equilibrium - model ' _ ( or ` bathtub - model ' ) @xcite that has been used to investigate the content of galaxies , both gaseous and stellar , the associated star formation rate ( sfr ) , metallicities , and the relations between them , especially for the ( blue )  main sequence " of galaxies .
these models typically assume an equilibrium relation ( see for example @xcite ) , @xmath9 , where @xmath10 is the sfr , @xmath11 is the gas mass of the disk and @xmath12 is the gas - to - star depletion time ( or  turn - over " time , @xcite ) .
in addition , the outflowing baryonic mass out of the galaxy , back to the circum- or intergalactic medium ( cgm / igm ) in these models is proportional to the sfr , @xmath13 , where @xmath14 is the mass - loading factor , set by energy and momentum injection back into the interstellar medium ( ism ) .    to date ,
all of the equilibrium models have been applied to galaxies as a whole , ignoring the internal disk structures .
in contrast , in this paper we formulate a radially resolved model , in which the equilibrium relations are  local " and functions of ( planar ) distance from the disk centres , and assuming that the gas dynamics is set by the mass - continuity equation and global angular momentum conservation .
our model is self - consistent , as star formation , disk rotation , outflows and the disk size are all set dynamically with the disk evolution .
we study in detail how systems that follow these simple assumptions evolve , reach a ( quasi-)steady state , accumulate stellar mass , and analyze dynamically their surface density and toomre-@xmath15 profiles . in our presentation , we show that many of our numerical results can be derived analytically , by solving simplified cases of the full ( partial differential ) equations , for our single phase ism .
the paper is organised as follows . in
[ sec : setup ] we write the evolution equations and the assumptions of our model , and discuss a dimensionless representation for it . in
[ sec : cossnap ] we consider a cosmological snapshot  in which we hold the cosmological conditions constant , and analyze the consequent steady state and the convergence to it .
next , in
[ sec : cosevol ] we generalise the cosmological snapshot , by accounting for the evolution of the halos , and varying baryon accretion rates . in  [ sec : results ] we analyze the results of our model for the cosmologically evolving system , and present a parameter study .
we summarise in  [ sec : discussion ] .
we consider an evolving , rotationally supported , star forming disk , within a cosmologically growing dm halo .
the halo accretes baryonic matter from its cosmological neighbourhood , at some given rate @xmath16 , forms stars continuously , and expels matter back to the cgm / igm at a rate proportional to the star formation rate .
our radially resolved equilibrium model is based on the dynamics set by mass - continuity and angular momentum conservation , together with local equilibrium relations .
the mass continuity equation in cylindrical coordinates ( assuming cylindrical symmetry ) includes cosmological accretion , and star formation and outflows as source and sink terms , and follows the gas surface density in the disk , @xmath17 .
it can therefore be written as @xmath18 where @xmath19 and @xmath20 are the
local " sfr and outflow rates respectively , @xmath21 is the radial gas velocity within the disk ( inflow velocity is @xmath22 ) , and where for the second equality we have assumed ` local ' equilibrium relations , @xmath23 to solve the mass - continuity equation , we will have to specify the boundary condition at some outer radius @xmath24 , and the gas radial inflow velocity @xmath25 at each radius .
we do that as follows :    * we assume that the entire ( thin ) baryonic disk , gas+stars , is confined to a radius @xmath26 , so that @xmath27 .
we determine @xmath26 by setting the specific angular momentum of the accreted gas onto the disk , at any time @xmath28 , equal to that of the parent halo with a given ( cosmological ) spin parameter @xmath29 where @xmath30 and @xmath31 are the angular momentum , energy and the virial mass of the halo , respectively , and @xmath32 is newton s gravitational constant .
observations and simulations suggest that the specific angular momentum of the baryonic disks indeed remain globally equal to the parent halos as the disks evolve despite the complex inflow and outflow patterns ( e.g. @xcite ) .
we therefore do not impose cosmological initial conditions for the cumulative mass distribution of specific angular momentum found in dark matter only simulations ( e.g. * ? ? ?
* in realistic disks , the radial gas inflows depend in a complex way on local conditions and global tidal torques and bar instabilities . in our model
we assume an azimuthally averaged gas inflow velocity , @xmath25 , at each radius @xmath33 , and we consider three cases .
the first two are ( _ a _ ) constant radial inflow velocity , @xmath34 , and ( _ b _ ) radial inflow velocity that satisfies @xmath35 , where @xmath36 is the inflow velocity at the outer disk radius @xmath26 .
the third case is motivated by a viscous disk in which the dissipated energy in one rotational timescale , which is of order @xmath37 , where @xmath38 and @xmath39 is the rotational velocity at each radius , is compensated by the gain in potential energy , @xmath40 @xcite , leading to ( _ c _ ) @xmath41 . our case ( _ c _ ) in combination with our mass conservation eqn .
( [ eq : dimen ] ) is similar to the viscous evolution models presented by @xcite and @xcite but with the inclusion of inflow and outflows , and self consistent star formation .    with the gas surface density @xmath42
we can keep track of the gas mass in the disk out to distance @xmath43 as a function of time @xmath44 the total gas mass in the disk is @xmath45 . to obtain the accumulated stellar mass @xmath46 one
has to integrate the sfr over time , and assume some stellar migration profile .
we will assume in this paper that once formed , stars do not migrate radially , but rotate around the disk center together with the inflowing gas , so that @xmath47 where @xmath48 denotes a time before star formation has begun .    solving for @xmath42 will also enable us to examine the local toomre-@xmath15 parameter @xcite , that quantifies the stability of the gas in the disk with respect to star formation , as a function of the galactocentric distance . by definition , @xmath49 where @xmath50 is the gas sound speed ( thermal plus turbulent ) and @xmath51 is the epicyclic frequency @xmath52 we approximate @xmath39 by the keplerian rotational velocity for a spherical mass distribution @xmath53 , that is composed of the gaseous and stellar disk components , plus the extended dark matter distribution @xmath54 we assume that the rotation curves for our two - dimensional disks differ only slightly than for spherical mass distributions ( see @xcite ) .
we ignore pressure support along the disk plane ( c.f .
@xcite ) .
it is useful to rescale to dimensionless form .
we define dimensionless time and distance coordinates , @xmath55 where @xmath56 is the crossing time , @xmath57 the normalised inflow velocity is @xmath58 we then define the dimensionless
accumulation " parameter @xmath59 and the normalised surface density @xmath60 thus in dimensionless form the mass - continuity equation is @xmath61 this dimensionless representation is valid only if @xmath16 and @xmath36 are constant , which is a reasonable assumption only if the disk reaches a steady - state faster than the time it takes @xmath16 to change significantly .
we assume this in  [ sec : cossnap ] , but relax this assumption in  [ sec : cosevol ] , which will force us back to eqn .
( [ eq : dimen ] ) .    in terms of this dimensionless notation ,
the fixed constant inflow velocity is written as @xmath62 , the @xmath63 case is written as @xmath64 and the @xmath65 case can be written as @xmath66 , where @xmath67 is a  dimensionless mass " , and the coefficient @xmath68 is set by the input parameters .
we start by solving the disk evolution equation for a  cosmological snapshot " , for which the halo mass and cosmological inflow rates are constant , and the dimensionless version of our model is valid .
we focus first on steady state solutions with (  [ sec : ysf ] ) and without (  [ sec : nsf ] ) star formation .
we demonstrate that simple mass conservation naturally leads to disks with exponential profiles over many length scales .
finally , in  [ sec : cca ] we examine the timescale for the convergence to steady - state .
quantitatively , the case of no star formation ( i.e. infinite depletion time ) , is considered simply by setting @xmath69 in eqn .
( [ eq : cont ] ) . in steady state ,
the inflow velocity is independent of time , and the continuity equation can be easily solved , yielding @xmath70 this solution is valid for @xmath71 or @xmath72 ( cases ( _ a _ ) and ( _ b _ ) ) discussed above . for @xmath66 ( case ( _ c _ ) ) eqn .
( [ eq : yqss ] ) is also satisfied but depends on @xmath73 , so this is actually a differential equation on its own , rather than a solution for the surface density . ignoring the dm halo , i.e. , assuming that the entire rotation curve , and hence the inflow velocity are determined solely by the baryonic matter in the disk
, this differential equation can also be solved analytically . using eqn .
( [ eq : m ] ) the differential equation for the surface density can be written as @xmath74 solving this differential equation , setting the constant of integration so that @xmath75 , and writing the solution in terms of @xmath76 , we get @xmath77 comparing this solution to eqn .
( [ eq : yqss ] ) shows that in the no - star formation no - dm case , the constant inflow velocity is a better fit to the physically motivated case ( _ c _ ) , in which the inflow velocity changes by @xmath78 per cent across the entire disk .    the dimensionless gas profiles , for all three inflow models , are generic results for @xmath69 , and do not depend on the accretion rate , the size of the disk or the sound speed .
it is instructive to display these solutions dimensionally for a specific ( characteristic ) cosmological accretion rate @xmath16 and disk inflow velocity @xmath36 . in fig .
[ fig : yqss ]        we plot the steady state solution for these three @xmath79-cases , in a dimensionless form ( bottom - left axes ) , and dimensional setting ( top - right axes ) , for @xmath80 , @xmath81 and @xmath82 , a cosmological accretion rate that corresponds to a uniformly distributed cosmological accretion by a @xmath83 median halo at @xmath3 @xcite .
it can be seen , that the surface density in all three cases , fits an exponential profile ( dashed lines , same colours ) quite well for the bulk of the disk . at the inner part of the disk
, however , the exponential profile ceases to fit the steady state solution , and in fact , the surface density diverges towards the center of the disk .
however , the enclosed mass @xmath84 for these cases does not diverge .
the gas profile shapes we find and display in figs .
[ fig : yqss ] and [ fig : y ] resemble those found assuming the `` strong form '' of angular momentum conservation , such that the cumulative distribution of mass versus specific angular momentum , @xmath85 , is conserved and equal to that of a uniform density sphere in solid body rotation ( e.g. * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
in contrast , our profile shapes are essentially set by radial inflow and mass conservation only ( eqn .  ( [ eq : dimen ] ) ) .
we interpret the divergence of the surface density towards the center as material that forms a bulge ( including a central black hole ) . in order to study the bulge to disk mass ratio ( b / d ) , defined as @xmath86 of our solutions ,
we first define the bulge radius , as the distance from the center of the disk at which the gas surface density is larger than the surface density predicted by the best fitting exponential profile by a factor of two .
the calculated disk scale radii , bulge radii and the corresponding b / d for the three @xmath79 test - cases are summarised in table  [ tab : thetable ] ,    .a summary the disk scale radius , bulge radius and the bulge to disk mass ratios for the three @xmath79-test cases discussed in the test , in steady state , with no star formation . [ cols="^,^,^,^",options="header " , ]     where it can be seen that the ( _ a _ ) and ( _ c _ ) disks , are indeed very similar , as predicted by comparing eqns .
( [ eq : yqss ] ) and ( [ eq : wow ] ) . enhanced stellar feedback efficiencies could remove the accumulated `` bulge material '' near the disk centres ( e.g. * ? ? ? * ; * ? ? ? * ; * ? ? ?
* ; * ? ? ?
* ) but we do not include such feedback effects in our computations .
so far we have considered the gas steady state solution without any star formation or outflows .
thus the only parameter that determines the gas surface density profile is the inflow velocity .
we now turn on the star formation circuit , which is also responsible for the outflowing gas from the disk back to the cgm / igm . in this section
we will focus on cases ( _ a _ ) and ( _ b _ ) , and will return to case ( _ c _ ) only when we take the dm component into account .    in the current treatment ,
fixing @xmath79 to the two limiting cases , there is no difference between star formation and outflows . from a physical point of view , though , these are different in essence , since while outflowing gas escapes the gravitational potential of the disk , forming stars stay in the disk , and contribute to the gravitational potential , the rotational velocity and the toomre-@xmath15 parameter .
we will distinguish between those two processes in
[ sec : cosevol ] .    solving the steady - state equation ( @xmath88 ) with @xmath87 for cases ( _
a _ ) and ( _ b _ ) allows us to quantify the effect of star formation on the steady state solutions for the surface density . although tedious , both cases are analytically solvable , and converge continuously to eqn .
( [ eq : yqss ] ) as @xmath89 .
we present the analytic solution for case ( _ a _ ) in appendix  [ sec : sssan0 ] .    by plugging in physically motivated values for @xmath90 and @xmath91
, we can get a feel of reasonable @xmath92-values . considering a @xmath93 disk , inflow velocity @xmath94 , depletion time of @xmath95 at @xmath3 @xcite , and @xmath96
( see for example @xcite ) we get @xmath97 .    in fig .
[ fig : y ] ( _ left _ ) we plot the steady state solutions @xmath98 for @xmath99 and @xmath100 of case ( _ a _ ) .
as @xmath92 increases the gas surface density of the bulk of the disk has a shallower slope .
this is because larger @xmath92 means more efficient star formation and less efficient inflow , or @xmath101 . in this limit
the evolution at every radius @xmath102 is almost independent , thus resulting in a more uniform disk . for small @xmath92 ,
as the depletion time becomes larger , gas flows towards the center rapidly , not staying long enough to form stars , until there is practically no star formation at all , and eqn .
( [ eq : yqss ] ) is relevant .
[ fig : y ] shows that ( quasi-)exponential disks form naturally when the radial gas crossing time is short compared to the star - formation gas depletion times .    the effect of @xmath103 on case ( _ b _ ) is similar .
moreover , even though the steady state profiles of the two @xmath79-cases with @xmath69 are significantly different across the disk , as @xmath92 increases both surface densities diminish , until they are barely distinguishable for @xmath104 .
having analysed the steady state solutions of eqn .
( [ eq : cont ] ) , we turn to solving its ( dimensionless ) time evolution , but still for a cosmological snapshot , i.e. assuming constant cosmological inflow rate and non - varying halo parameters . for simplicity
, we focus on the @xmath69 case . by carrying out this analysis
, we will be able to estimate the dimensionless time @xmath105 it takes the disk to reach a steady state .
as we already noted , the cosmological snapshot solution should be a reasonable approximation , if the steady state convergence timescale is short compared to the time during which the cosmological conditions change significantly .    in fig .
[ fig : y ] ( _ right _ ) we plot the dimensionless surface density @xmath106 , for case ( _ a _ ) , at a few @xmath107 intervals .
the gas surface density @xmath108 reaches a steady state by @xmath109 ( the @xmath110 case reaches a steady state even sooner , within @xmath111 ) , that is within one crossing time , converging to the analytic solution of eqn .
( [ eq : yqss ] ) , starting off with an empty disk . considering , for example ,
a @xmath112 and @xmath113 disk , means that the system reaches a steady state within half a gyr . at @xmath114 , the dm halo ( and therefore the cosmological accretion onto the disk as well )
change significantly within half a gyr . at later times , however , the cosmological accretion rate changes only slightly , meaning that a quasi steady state is maintained .    in addition , we also find that the convergence to the steady state surface density is an  outside - in " process . as a consequence ,
the fact that the radial surface density profile is well approximated by an exponential profile is not true throughout the evolution of the disk .
this approximation becomes better as the evolution approaches its quasi steady state .
so far we have considered halos of constant mass , and thus constant baryonic accretion rate , to which we referred to as a cosmological snapshot .
as we have shown , within this picture the disk evolution equation can be written in a purely dimensionless form and solved analytically .
we now generalise our model by accounting for the evolution of those cosmological conditions as well as other galactic parameters and local star formation , which is a necessary step for the model to be self consistent :          * * dm halo *  the dm halo is a crucial ingredient in disk formation and evolution for a couple of reasons .
first , the halo growth rate , that depends on @xmath115 sets the baryonic accretion rate onto the disk . in our model
we follow @xcite , taking the baryonic accretion rate as @xmath116 , where the halo growth rate is @xmath117 and can be integrated to give the exact halo mass at redshift @xmath118 , assuming that its mass at some @xmath119 is known , as we show in appendix
[ sec : mhofz ] .
+ second , the dm halo also affects the rotation curve of the baryonic disk . as a result not only
the mere total halo mass is important , but also its spatial distribution @xmath120 .
furthermore , due to the dependence of the rotation curve on the halo mass distribution , the inflow velocity and the disk size will also be affected .
therefore , to account for the mass distribution correctly , we must specify the halo concentration @xmath121 as a function of @xmath115 and @xmath118 .
we assume continuously evolving median concentration halos following @xcite .
we neglect adiabatic contraction of the dark matter halos * * sound speed *  since in our model the gas inflow velocity is @xmath122 and @xmath123 , the gas sound speed is an essential parameter .
higher sound speed means gas that flows inwards more efficiently , resulting in low @xmath42-values at large @xmath33 , with the addition of a larger numerator in @xmath124 .
therefore , we expect that the larger the gas sound speed is , the harder it gets for the disk to form stars , and the smaller the stellar disk becomes .
+ we follow the redshift evolution of the sound speed , as in @xcite , i.e. @xmath125 as noted in @xcite , this relation is valid over a wide range of redshifts
. however , due to measurement uncertainties , @xmath50 shows a substantial scatter around this value ( for example @xmath126 at @xmath127 , and @xmath128 at @xmath3 ) .
therefore , we will use eqn .
( [ eq : cs ] ) as our fiducial sound speed in  [ sec : fid ] , and study how our model predictions change when considering different @xmath50 in  [ sec : bfm ] , keeping the @xmath129 dependence on redshift . *
* disk size *  to constrain the disk size , we assume that the specific angular momentum of the accreted matter , @xmath130 , equals that of the parent dm halo , @xmath131 , where the accreted gas is assumed to be in centrifugal equilibrium with the disk , i.e. it starts rotating around the center of the disk at a rotational velocity @xmath132 immediately as it enters the disk . not making this assumption , necessarily means that other assumptions , about the time it takes the newly accreted gas to reach the disks rotational velocity and the amount of energy that is dissipated in this process , have to be made .
it would also cause a situation in which at each radius @xmath33 , the gas rotation is inhomogeneous .
since we assume that the cosmological accretion at each time is uniform across the disk , then the constraint on the conservation of the angular momentum reads @xmath133 therefore for a given @xmath134 , since we know the halo mass at each redshift , we can solve eqn .
( [ eq : amconstraint ] ) for the disk size @xmath135 . with eqn .
( [ eq : amconstraint ] ) we are essentially adopting the results of both observations ( e.g. @xcite ) and cosmological hydrodynamic simulations ( e.g. * ? ? ? * ; * ? ? ? * ) that show that specific angular momentum is conserved from halo to disk .
this despite the different angular momentum histories of the gas and dm .
thus , we only specify the spin parameter @xmath134 , but not the cumulative mass distribution of the specific angular momentum , @xmath136 , found in dark matter only simulations as cosmological initial conditions ( e.g. * ? ? ?
* ; * ? ? ?
* * depletion time *  at the base of our radially resolved equilibrium model lies the equilibrium relation ( [ eq : eq ] ) , that relates the sfr to the gas content of the disk .
the proportionality , @xmath91 , quantifies the efficiency of star formation  the shorter @xmath91 is , the more efficient star formation becomes , which can proceed even if the gas in the disk flows rapidly inwards due to a high inflow velocity .
for the analysis below , we take @xmath137 following @xcite . *
* star formation *  it has already been noted by @xcite that @xmath15 in the disk will self - regulate to maintain @xmath138 .
the intuitive explanation for this self - regulation in @xcite is generally valid in the spatially resolved case as well .
nevertheless , resolving the disk radially , means that @xmath139 and @xmath15 are all calculated as a function of the distance from the center of the disk . wherever the ratio between the epicyclic frequency and @xmath42 is small enough
, @xmath15 will reach its threshold value , and self regulation will proceed . at the outer regions of the disk , however , the low gas surface density will dominate @xmath124 , which will consequently be @xmath140 , and therefore star formation will not be initiated as expected for the outer lower surface density portions of galaxy disks ( e.g. * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
this does not imply that regions in the disk where no star formation is taking place at some redshift will remain quiescent .
however , due to the cosmological evolution of the parent halo , and as a result of the baryonic accretion rate and the size of the disk , star formation can begin at a later time .
+ in order to account for star formation self - consistently in the analysis below , we calculate @xmath124 , and turn the star formation - circuit on only at those regions of the disk in which @xmath1 .
in this section we consider cosmological evolution of the halo and the accretion rate , employ a self consistent star formation criterion , and solve numerically the evolution equations for the gas and the stars in the disk .
first , in  [ sec : fid ] , we define and analyze thoroughly our fiducial model .
next , in
[ sec : bfm ] we examine how our results depend on model parameters , @xmath141 and @xmath14 for a fixed ( at @xmath3 ) halo mass . finally , in
[ sec : bfm2 ] we study how the model results change when the halo mass is modified .          as our fiducial model
we assume a @xmath142cdm cosmology , with @xmath143 , @xmath144 , and baryon fraction @xmath145 @xcite .
we also assume that the dm halo profile follows a navarro - frenk - white profile @xcite with median concentration parameter that varies with redshift , following @xcite .
the fiducial dm halo we consider is of mass @xmath146 at @xmath3 and spin parameter @xmath147 . in addition , we assume @xmath148 , no instantaneous recycling , and @xmath149 at @xmath3 within the disk .
our massive fiducial halo represents the tip of the galaxy main - sequence at @xmath150 .
solving for the cosmologically time - dependent fiducial halo , in which the star formation term is turned on only when and where @xmath1 , we plot in fig .
[ fig : fid ] ( _ top left _ ) the resulting surface density of the two disk components at @xmath3 , and the corresponding @xmath124 in fig .
[ fig : fid ] ( _ top right _ ) . as can be seen in the figure , for the fiducial model both gas and the stellar profiles are well approximated by exponential profiles ( dashed lines ) , for a wide range of radii .
for the stars and gas distributions , exponential profiles with scale lengths equal to @xmath151 and @xmath152 respectively , represent the computed distributions within a factor of two , from 0.5 to 9 @xmath153 for the stars , and from 1 to 11 @xmath153 for the gas .    remarkably , though somewhat fortuitously , the stellar length scale precisely equals the length scale given by @xmath154 for assumed exponential disk in an isothermal halo @xcite .
according to our bulge - criterion , the fiducial model bulge extends to @xmath155 , and the bulge to disk mass ratio is @xmath156 .
in addition , while the gaseous disk extends roughly to @xmath157 , the stellar disk is somewhat more compact , stretching out to @xmath158 . the gas half mass surface density , defined as @xmath159 is @xmath160 , where @xmath161 is the the distance form the center of the disk such that @xmath162 the stellar half mass surface density , defined in a similar way is @xmath163 .
at @xmath3 we find @xmath164 , consistent with observed sizes for this stellar mass and redshift @xcite .
the star formation rate ( at @xmath3 ) for the fiducial model is @xmath165 , and the gas fraction is @xmath166 .    comparing our predicted sfr to the sfr - main sequence , we find that the predicted sfr for our fiducial halo , @xmath167 , is in very good agreement with the sfr according to @xcite , @xmath168 , for the same redshift @xmath3 , and stellar mass @xmath169 .
using our model we can also follow the actual evolution of the disk , as a function of redshift . in fig .
[ fig : fid ] ( _ bottom row _ ) we plot the surface density of the gaseous ( _ left _ ) and stellar disks ( _ right _ ) for @xmath170 and @xmath171 .
it can be seen that even though that the characteristic gas surface density reaches its  final " ( @xmath3 ) value already at @xmath172 , the disk growth continues even at smaller redshifts .
studying the stellar disk surface density evolution , we find that star formation begins in the inner regions of the disk at @xmath172 , and propagates outwards with time , until at @xmath3 the stellar disk is almost double in size with respect to @xmath173 , or @xmath174 later .    in fig .
[ fig : msofz ] ( _ left _ ) we plot the redshift evolution of the total mass for each of the disk components ( gas and stars ) , together with the dm halo and hot gas for our fiducial model , where for the hot gas we sum all the gas that flows through the disk and out from the center , plus the accumulated outflow that is proportional to the total stellar mass .
most of the baryons remain as hot halo gas throughout the evolution . at @xmath172 , the onset of star formation in the fiducial model case
, the disk gas mass growth slows down , as is to be expected from the equilibrium relations .
the effect on the hot gas mass , however , is insignificant since most of the hot gas is ejected from the disk center .    in fig .
[ fig : msofz ] ( _ right _ ) we plot the rotation curve @xmath175 ( black curve ) together with the baryonic and dark matter contributions to the total mass ( solid and dashed red curves ) .
the disk is baryon dominated throughout .
it is also of interest to plot the resulting cumulative mass versus specific angular momentum , @xmath85 , for our fiducial model . in fig .
[ fig : mlessthanj ]     versus @xmath176 for our fiducial model ( _ solid blue _ ) and @xmath177 ( _ dashed black _ ) for a sphere in solid body rotation . ]
we plot @xmath178 versus @xmath176 for our fiducial model ( solid curve ) together with the relation @xmath177 for a sphere in solid body rotation .
the two curves differ only slightly .
this is consistent with our resulting ( quasi)-exponential disks , and the known similarity between the angular momentum distributions of centrifugally supported exponential disks and rotating spheres ( e.g. * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ) .
that allows star formation at @xmath3 for a @xmath83 ( _ blue _ ) and a @xmath179 ( _ dark - blue _ ) halo as a function of the gas sound speed . the fiducial model ( @xmath180 ) is designated by a small black circle
only the regions below the curves are star forming .
no star formation occurs above the blue curves .
_ top - right axes : _ @xmath134 dependence of the total sfr for the @xmath181 case .
it can be seen that the sfr dependence on @xmath134 is weak , changing by less than @xmath182 within the @xmath183 range . ]
we wish to examine how the disk properties are modified when the halo or disk parameters are altered .
we note that a key advantage of our approach  a simple toy - model of disk evolution  is that it enables scanning through a large parameter - space at a low computational cost .
we focus on @xmath83 halos at @xmath3 , as for our fiducial model .
we examine how the redshift at which star formation begins , @xmath184 , depends on @xmath50 and @xmath134 .
since the initial conditions for our disk evolution are of a median halo and an empty disk that is fed by cosmological accretion , star formation can not take place until enough gas has been accreted and @xmath42 is large enough .
the growing disk size , on the other hand , makes @xmath42 smaller , thus delaying the onset of star formation . the impact of having a large @xmath50 is threefold .
first , the large inflow rate empties the gas at each radius more efficiently , resulting in reduced @xmath42 .
second , it delays star formation due to the @xmath15-barrier it induces .
third , as we show below , although counter - intuitive , larger @xmath50 results in larger disks , and hence smaller @xmath42 , thus delaying star formation .
even though larger @xmath134 can also delay star formation , through the disk size , for @xmath134 from @xmath185 to @xmath186 , which we consider , its effect is weaker , as seen in fig .
[ fig : cs_l_first ] ( _ top left _ ) .
interestingly , after star formation has begun , @xmath134 hardly has an effect on sfr , despite its effect on the disk size .
qualitatively , increasing @xmath134 means increasing the angular momentum of the halo , and hence also of the accreted baryons .
therefore , as @xmath134 increases , the disk size @xmath24 increases too .
however , since the baryonic accretion onto the disk @xmath187 is independent of @xmath134 , then to first order the total baryonic mass in the disk will also be independent of @xmath134 . consequently ,
even though large @xmath134 means smaller characteristic @xmath42 , as long as it is high enough to maintain @xmath138 , star formation , which is @xmath188 should proceed at more or less the same rate , as seen in fig .
[ fig : cs_l_first ] ( _ top right _ ) .
nevertheless , there should exist a @xmath189 , that once reached , renders @xmath42 too low and @xmath15 too high , which will not allow any more star formation . in fig
[ fig : lambdamax ] ( _ bottom - left axes _ ) we plot @xmath189 that allows star formation at @xmath3 , as a function of the gas sound speed , for our fiducial halo ( @xmath83 at @xmath3 ) , and for a slightly less massive halo ( @xmath179 ) .
it can be seen that @xmath190 , as preferred by observations and simulations is well within the star forming portion of the parameter space , and becomes borderline only as the gas sound speed increases significantly . in fig  [ fig : lambdamax ] ( _ top - right axes _ )
we plot the total sfr for the fiducial sound speed , as a function of @xmath134 .
it can be seen that for @xmath134 altered by more than an order of magnitude , the total sfr changes by less than @xmath100 per cent , following our qualitative argument .
studying the disk size , one might expect that since @xmath191 , then for higher sound speeds , the gaseous disk will be more compact , as the radial migration of the gas becomes more efficient .
however , the disk size is primarily set by the baryonic accretion onto the disk , whose extent , in our model , is set by conservation of angular momentum .
let us consider a halo of mass @xmath115 , with some specific angular momentum , @xmath131 that is , of course , independent of the sound speed .
a larger @xmath50 at the early stages of the disk evolution , will result in a mass profile that is more concentrated near the center of the disk , and thus lower baryonic specific angular momentum .
therefore , to satisfy the angular momentum constraint , eqn .
( [ eq : amconstraint ] ) , the disk must grow for higher sound speeds .
this is seen in fig .
[ fig : cs_l_first ] ( _ bottom left _ ) , where we plot @xmath26 ( the truncation radius of the baryonic disk ) as a function of the gas sound speed and @xmath134 , for our fiducial halo at @xmath3 .    for the stellar disk size dependence on @xmath50 , which is plotted in fig .
[ fig : cs_l_first ] ( _ bottom right _ ) , we note that there are two competing effects . as @xmath50 increases , eqn .
( [ eq : q ] ) implies that in the outer regions of the disk , where @xmath42 is lower , star formation is suppressed .
but , since the gaseous disk grows as @xmath50 grows , then the absolute size of the stellar disk is potentially larger , at least as long as overall @xmath15 is small enough .    in fig .
[ fig : cs_l_more ] we plot the corresponding @xmath192 ( _ top row _ ) and @xmath193 ( _ bottom row _ ) , the half mass radius and surface density of the gas ( _ left column _ ) and stars ( _ right column _ ) .
it can be seen that regions in parameter space that result in small @xmath26 and @xmath194 have large @xmath193 .
this is because small disks get a larger mass contribution from the central divergence of the surface density .
the opposite is not exactly satisfied  the largest disks do not have the lowest half mass surface densities . for the gas component ,
larger disks contain more baryons , and hence their characteristic surface density is higher , while for the stellar component not all of the parameter space results in disks that form stars , resulting in lower characteristic surface densities .
owing to these effects , the total disk mass , as well as the relative gaseous and stellar contributions to it depend on @xmath50 and @xmath134 in a non - trivial way . as seen in fig .
[ fig : cs_l_msandmore ] ( _ top left _ ) , the total disk mass increases as the sound speed decreases , as one would expect , primarily due to the fact that higher @xmath50 results in a rapidly emptying disk , where the excess gas that inevitably reaches the center of the disk and ekected back to the hot halo .
the hot halo is also fed by the outflowing baryons , whose mass @xmath195 sfr .
the total hot gas mass due to these contributions , is plotted as a function of @xmath196 for our fiducial halo at @xmath3 in fig .  [
fig : cs_l_msandmore ] ( _ top right _ ) .
it can be seen that the hot gas mass depends very weakly on @xmath134 , similarly to the sfr .
this implies that outflows constitute a significant part of the total hot gas mass .
the relative contributions of gas and stars to the total disk mass can be read off fig .
[ fig : cs_l_msandmore ] ( _ bottom left _ and _ bottom right _ ) ,          where we plot the gas fraction @xmath197 , and @xmath198 , respectively .
it can be seen that except for the fact that the gas fraction is lower where sfr is higher , it exhibits a stronger dependence on @xmath134 than the sfr , since as argued above higher @xmath134 generates larger disks , that results in a smaller portion of it that forms stars .
apart from @xmath134 and @xmath50 , another free parameter in our model is the mass loading factor @xmath14 .
considering variation in @xmath14 from @xmath199 to @xmath200 results in insignificant modifications to the results reported above , and we thus do not pursue this @xmath14 analysis further .
we note , however , that analysis of the extreme case @xmath201 , i.e. no outflows at all , shows a stronger dependence of sfr and @xmath202 on @xmath134 , verifying our argument about the importance of outflows to the hot gas mass .
as we have already mentioned , the tight scaling relation between the sfr in disk galaxies and their stellar mass is remarkable .
this scaling relation has been repeatedly observed , in many independent observations and different surveys .
therefore , it is interesting to compare our model predictions to actual observations of sfr , and how well those predictions fall within the galaxy main sequence . in particular , we compare our model predictions to the best - fitting sfr-@xmath203 relation , as found in @xcite , for star forming galaxies between @xmath127 and @xmath171 @xmath204 where @xmath205 since the initial conditions of our simulated disks are of empty disk , and since , as we have shown , in our model star formation is initiated at @xmath172 , it is clear that the correspondence between our model predictions and eqn .
( [ eq : whit ] ) can not hold for any high redshift , and thus we are interested in comparing these sfr s at a lower redshift , so that the disk has already reached an equilibrium . in fig .
[ fig : sfrscatter ]     relation ( or `` main sequence '' , as predicted by our model for various model parameters ( @xmath206 ) at @xmath3 , and the corresponding observed relation ( dashed ) according to @xcite . ]
we plot the whitaker relation , eqn .
( [ eq : whit ] ) , as a function of @xmath203 for @xmath3 ( _ black dashed _ ) compared to our model predictions , all at @xmath3 and for a halo of mass @xmath83 , and various model parameters , @xmath207 and @xmath14 .
it can be seen from this plot that our model predicts very well the sfr-@xmath203 relation , even as @xmath196 and @xmath14 are varied .
so far we have only considered the evolution of disks within a dm halo of mass @xmath83 at @xmath3 , and examined how perturbations to model parameters affect the evolution of those disks .
we now wish to consider less massive halos .
these imply smaller baryonic accretion rates , and hence less massive baryonic disks .
this affects the disk size and sfr , and consequently also the gas fractions in the disks .
reduced halo masses lead to smaller disks as less gas is accreted , but as shown in fig .
[ fig : changem ] ( _ top left _ ) , however , the dependence of the disk size on the halo mass weakens as @xmath50 decreases .    since a less massive halo means reduced baryonic accretion , and hence lighter gaseous disks , it is also to be expected that a smaller @xmath115 implies smaller sfr , for a given @xmath50 .
as we have already shown , as @xmath50 increases , the characteristic gas surface density has to become larger to sustain star formation .
therefore , for each @xmath50 there is the corresponding @xmath115 below which star formation is no longer possible , since @xmath208 at all times for all radii .
this forms a star - formation threshold , which is visible in fig .
[ fig : changem ] ( _ bottom left _ ) of the sfr as a function of @xmath50 and @xmath115 at @xmath3 .
this naturally points out a validity range for our model .    examining the stellar disk size as a function of @xmath115 and @xmath50 , as in fig .
[ fig : changem ] ( _ top right _ ) , we find that the behaviour we found in the analysis of fig .  [
fig : cs_l_first ] ( _ middle right _ ) , of our fiducial halo , is modified for less massive halos . for halos of mass @xmath209
the stellar disk size is monotonically decreasing as @xmath50 increases .
this implies that for the less massive halos the fact that the total size of the gaseous disk increases as @xmath50 increases is subdominant to the @xmath15-barrier to star formation due to the increase of the gas sound speed .    finally , in fig .
[ fig : changem ] ( _ bottom right _ ) we plot the total stellar mass in the disk as a function of @xmath115 and @xmath50 at @xmath3 .
we find that even though the total stellar mass increases monotonically with @xmath115 , as might be expected , it does not follow the general behavior of the stellar disk size discussed above , which means that maximising the star forming region does not guarantee a higher total sfr through that region .
in this paper we have presented a
radially - resolved - equilibrium - model " for the growth of baryonic disks in cosmologically accreting massive halos .
we have solved the time - dependent equations that govern the radially dependent gaseous and stellar contents of disks , as well as their star - formation rates , inflows and outflows from and to the inter- and circum - galactic medium , and inward radial gas flows within the disks , while imposing mass and angular momentum conservation .
we have focused on the evolution up to redshifts @xmath0 . the equilibrium model for galaxy evolution
has been widely employed to account for the interrelated observed parameters of the ( blue ) main - sequence @xcite . in this paper
we have generalised the analytic model to account for the global disk structures  their gaseous and stellar surface densities , radially resolved sfr  and study the impact of the variation of various parameters on this structure .
we have shown that :    * the dynamical time scales required for both components of baryonic disks to reach their equilibrium configuration are much shorter than the those of cosmological variations in the dm halo mass and profile .
this feature guarantees that the disk structures and star formation rates remain close to equilibrium during the cosmological evolution . *
radial ( in plane ) gas inflows together with mass conservation naturally give rise to ( quasi-)exponential disks over many length scales .
in our models the disk sizes , surface densities , and star formation rates are self - consistently computed assuming global conservation of the halo specific angular momentum . *
the stellar disk size is determined by the sfr , which is self consistently calculated according to the condition toomre @xmath1 . for a roughly flat rotation curve
, this condition is not satisfied when the surface density becomes too small , as in the outer parts of the gaseous disk , which naturally leads to stellar disks that are more compact than the corresponding gaseous disks .
* focusing on cosmological evolution up to @xmath3 , and on massive , @xmath2 halos ( at @xmath3 ) , our model reproduces the star formation main sequence , @xmath210 .
in particular , our fiducial model predicts sfrs of @xmath4 , stellar masses @xmath5 , gas contents @xmath6 , half mass sizes of 4.5 and 5.8 kpc for the stars and gas , and characteristic surface densities @xmath211 .
our analytic models may be used to study the interrelationships between spin parameter , halo mass , and redshift , on the disk gas and stellar surface densities and radially dependent star - formation rates .
we thank andreas burkert , reinhard genzel and linda tacconi for many helpful discussions .
we thank the referee for very helpful comments that improved our presentation .
this work was supported by the dfg via german - israel project cooperation grant ste1869/1 - 1/ge625/15 - 1 .
we also acknowledge support by the raymond and beverly sackler tel aviv university - harvard / itc astronomy program , and israel science foundation i - core grant 1829/12 .
for the velocity case ( _ a _ ) , the steady state solution of the @xmath212 ( that is , the star forming ) scenario may be obtained analytically . to do that
, we first rewrite the differential equation , as a classical non - homogeneous first order ordinary differential equation of the form @xmath213 : @xmath214 identifying the coefficients @xmath215 we can immediately find the solution to this equation , with the constant of integration chosen so that @xmath216 , to yield @xmath217 it can be easily verified that for @xmath218 @xmath219 which converges continuously to the no - star formation case .
another notable fact is that having @xmath87 does not change the functional form of the divergence of @xmath108 towards the center of the disk , where the surface density satisfies @xmath220 where @xmath221 can be read off eqn .
( [ eq : yofalpha ] ) and satisfies @xmath222 .
finally , we note that the same is true for case ( _ b _ ) : the dimensionless surface density may be calculated analytically , and shown to converge continuously to the @xmath69 scenario .
let us write eqn .
( [ eq : dekel ] ) in a more general way @xmath224 we note that @xmath225 where @xmath226 is the cosmological scale factor and @xmath227 are the cosmological redshift and time as usual
. therefore eqn .
( [ eq : dekelgeneral ] ) can be written as @xmath228 where @xmath229 is the hubble constant at @xmath118 .
this equation can be solved numerically , as we do for the quantitative analysis in this paper .
it is illustrative though to assume a matter dominated universe ( @xmath230 ) , so that @xmath231 , where @xmath232 is the present day hubble constant , and solve this differential equation analytically . in a matter
dominated universe then @xmath233 where @xmath234 .