 other galaxies even if it was an isolated system at birth. The predictions in the one-zone model therefore may be significantly different from what happened to galaxies in a hierarchical universe.

White and Frenk (1991) developed a detailed analytic formalism to describe the formation and evolution of galaxies while taking account of the hierarchical merging of dark-matter halos, gas cooling, star formation, and supernova feedback. Subsequent numerical approaches in modeling hierarchical merging of dark halos employ two somewhat different algorithms; one is called the ‘block model’ in which a random-Gaussian density fluctuation field is generated by dividing a hypothetical rectangular box recursively ([@CK88]; [@Cole91]; [@Cole94]). While this algorithm is simple and straightforward, the resulting halo masses are necessarily binned in discrete steps of a factor of two. The other generates a realization of halo merger trees according to a probability distribution function predicted by the extended Press–Schechter theory ([@Bower91]; [@Bond91]; [@KW93]; [@SK99]; [@SL99]). The latter is widely used in studying the cosmological evolution of galaxies in a hierarchical universe ([@KWG93]; [@baugh98]; [@SP99]; [@Cole00]; [@nagashima]). Throughout the present paper, we call the latter method the Monte-Carlo modeling of merger histories (simply, the Monte-Carlo modeling), while it is usually referred to as a semi-analytic model of galaxy formation (SAM).

The most important ingredient in Monte-Carlo modeling is the conditional joint-probability distribution function of a set of *progenitor* halos of mass $M_2^{j}$ at a redshift of $z_2$, which is a part of a *parent* halo of mass $M_1$ at $z_1$, conceptually written as $$\begin{aligned}
\label{eq:jointprob}
{\rm Prob}(M_2^1, M_2^2, \cdots, M_2^N, z_2 |
M_1, z_1) dM_2^1 dM_2^2 \cdots dM_2^N \cr
\qquad (N=1, \cdots, \infty).  \end{aligned}$$ Unfortunately only an analytical expression for the conditional one-point probability distribution function, Prob($M_2^i$, $z_2 |M_1$, $z_1$), is known based
other galaxies even if it was an isolated system at parturition. The prediction in the one - zone model therefore may be significantly unlike from what happened to galaxies in a hierarchical universe. 

 White and Frenk (1991) develop a detailed analytic formalism to describe the formation and development of galaxy while taking score of the hierarchical merging of dark - topic halos, gas cooling, star formation, and supernova feedback. Subsequent numerical approaches in model hierarchical merging of dark halos hire two somewhat different algorithms; one is called the ‘ pulley model ’ in which a random - Gaussian density fluctuation field is generated by divide a hypothetical rectangular box recursively ([ @CK88 ]; [ @Cole91 ]; [ @Cole94 ]). While this algorithm is simple and straightforward, the result halo masses are necessarily binned in discrete steps of a factor of two. The other generates a realization of halo merger tree according to a probability distribution function bode by the extensive Press – Schechter theory ([ @Bower91 ]; [ @Bond91 ]; [ @KW93 ]; [ @SK99 ]; [ @SL99 ]). The latter is widely used in study the cosmologic evolution of galaxies in a hierarchical population ([ @KWG93 ]; [ @baugh98 ]; [ @SP99 ]; [ @Cole00 ]; [ @nagashima ]). Throughout the present paper, we call the latter method the Monte - Carlo modeling of merger histories (simply, the Monte - Carlo modeling), while it is usually referred to as a semi - analytic model of galaxy constitution (SAM). 

 The most significant ingredient in Monte - Carlo modeling is the conditional joint - probability distribution function of a set of * progenitor * halos of mass $ M_2^{j}$ at a red shift of $ z_2 $, which is a part of a * parent * halo of mass $ M_1 $ at $ z_1 $, conceptually written as $ $ \begin{aligned } 
 \label{eq: jointprob } 
 { \rm Prob}(M_2 ^ 1, M_2 ^ 2, \cdots, M_2^N, z_2 | 
 M_1, z_1) dM_2 ^ 1 dM_2 ^ 2 \cdots dM_2^N \cr 
 \qquad (N=1, \cdots, \infty).   \end{aligned}$$ Unfortunately only an analytical expression for the conditional one - point probability distribution function, Prob($M_2^i$, $ z_2 |M_1 $, $ z_1 $), is known based
 otjer galaxies even if it das an isolated system at bidth. The oredictions in the one-zone mldwl thtgefore may be signifizantly divferent drom qhat happeisd to galaxiea in c iierarchical unlverse.

White and Frenk (1991) deeeuo'ed a detailed analytic formalism to descrine the formation and qvolhniin of galaxies while taking adcount mf the hierarvhical merging of dark-mattfr hwlos, gas cooling, shar formatiin, age supernova weedback. Slysequent nujerical approaches in modeling fieraxchical mertibg lx dark halow empjoy two somewhat difxerent slgorithms; one is caoled the ‘block model’ mn which a random-Gaufsian denvicy fluctuation field us genetated by auviainf e hgpothehicel rectangumar box recyrsively ([@CK88]; [@Cole91]; [@Coke94]). Qhile this alforithi ys simple and straightforward, the resuluing galo masses are necessaeily binned in discreje steps os a factor of two. The other generates a realizatimn of fali idegfr trees according to a probability distributykn flnction predicted by the extrnfec Press–Schechtgr theoxg ([@Gower91]; [@Bond91]; [@KW93]; [@SK99]; [@SP99]). The lwtter is widelr usrd in studying the cosmologucal evolutijb of galaxies in a hierarchiccl uniferse ([@KWG93]; [@baugh98]; [@SP99]; [@Cole00]; [@nagarhimz]). Throughouh the preadnt paper, we calu tme natter mtghod the Monte-Carjo modeliig of merger hisjories (fimply, the Montc-Warlo modeling), whipe it iv usually geferred to as a semi-analytic model of galaxy fmrmdtion (SAI).

The kost importane ingredient iu Monte-Ccrlo mudeling is the coiditional joynt-probabilitf distributioi functiog of a swt of *pfugenitor* halos of mass $M_2^{j}$ at a reeshift of $z_2$, which ls a kadt of a *parent* kcli of mass $M_1$ at $a_1$, cuncqpnuaklr written as $$\tegiv{alkbned}
\lxbel{eq:jointiroc}
{\rm Lrob}(M_2^1, M_2^2, \cdots, M_2^N, z_2 |
M_1, z_1) dJ_2^1 dM_2^2 \cdots dM_2^N \cr
\qauwd (N=1, \cdojs, \infty).  \end{aligned}$$ Umfortunately only wn anelyticel exptesfion for the conditional one-pojnt probahillty distributyon nuncjion, Prob($M_2^i$, $z_2 |M_1$, $z_1$), is known based
other galaxies even if it was an at The predictions the one-zone model from happened to galaxies a hierarchical universe. and Frenk (1991) developed a detailed formalism to describe the formation and evolution of galaxies while taking account of hierarchical merging of dark-matter halos, gas cooling, star formation, and supernova feedback. Subsequent approaches modeling merging dark halos employ two somewhat different algorithms; one is called the ‘block model’ in which a density fluctuation field is generated by dividing a rectangular box recursively ([@CK88]; [@Cole94]). While this algorithm is and the resulting masses necessarily in discrete steps a factor of two. The other generates a realization of halo merger trees according to a probability function predicted extended Press–Schechter ([@Bower91]; [@KW93]; [@SL99]). The latter used in studying the cosmological evolution a hierarchical universe ([@KWG93]; [@baugh98]; [@SP99]; [@Cole00]; [@nagashima]). the present we call the latter method the modeling of merger histories (simply, the Monte-Carlo modeling), it is usually referred to as a semi-analytic model of galaxy formation (SAM). The most in Monte-Carlo modeling is conditional joint-probability distribution of set *progenitor* of mass at a redshift of $z_2$, which is a part of a halo of mass $M_1$ at $z_1$, conceptually written as $$\begin{aligned} Prob}(M_2^1, \cdots, M_2^N, z_2 M_1, z_1) dM_2^1 dM_2^2 dM_2^N \qquad (N=1, \cdots, \infty). only analytical conditional probability function, Prob($M_2^i$, $z_2 |M_1$, is known based
 other galaxies even if it was aN isolated sYstem At bIrtH. THe prEdicTions in the one-zONe moDel therefore may be signiFicanTlY DiffEReNt froM what haPPeNED to GaLaXieS iN A hIerarChiCal univErse.

White aNd FReNk (1991) developed a DEtAiled analyTic Formalism to dEscRibe thE fOrmATion aNd eVolutIon of gALaxies While takiNg ACcount OF the hieRARcHicaL merging of dark-matTEr HAlos, gas cooling, Star foRmATiON, And SupErnova feedBaCk. SubSEquent nUMeRICAl aPProaches in modEling hierarCHicAl mergInG of DArk halOs empLoY Two Somewhat difFereNt algoritHms; one IS called THe ‘block Model’ iN whIch A ranDOm-gaUssIaN DenSItY flUCtuAtion fieLd Is GenerAted BY DIVidiNg a HypoThetiCal rectangulaR boX recURsiVely ([@Ck88]; [@Cole91]; [@cole94]). whIle thIs algoRithm Is Simple and straigHtfoRward, the rEsuLtIng HaLo masSEs are nEceSsaRily binNed in diSCreTe STEPs Of a factor of two. The oThER GeNerates a RealizATiOn OF halo merGeR trEes aCCOrdinG to a PRoBability DistriBUtIoN functiOn PredicTeD by The ExtenDEd PrEss–SchEchter thEory ([@BOWer91]; [@Bond91]; [@KW93]; [@SK99]; [@SL99]). THE latter is wideLY uSED iN StudYinG the cosmoloGicaL EvolUtioN Of GalAXies iN a hieRaRChICal universe ([@KWG93]; [@baugh98]; [@sP99]; [@cole00]; [@naGashiMa]). Throughout tHe present pAPER, we call tHe laTTeR Method the Monte-carlo Modeling of MErger hisTorieS (simply, tHe Monte-CaRLO modelinG), whIle It iS usUALlY referred to as A SEmi-aNaLytic moDel Of galaxY foRmaTioN (SAm).

THe most impOrtant inGrEdIeNt In MOnte-CARlo modelInG is ThE coNditiONal joiNt-proBabiLiTy DIstRibutioN FuNCTion Of A sEt of *ProGeNitor* HaloS Of mAss $M_2^{j}$ at A redshift Of $z_2$, WHich Is A pArt of a *pArent* halo of maSs $m_1$ at $z_1$, concepTuAllY writtEN As $$\begin{aLigned}
\label{eq:jointprob}
{\rM prob}(M_2^1, M_2^2, \cDotS, M_2^N, z_2 |
M_1, Z_1) dM_2^1 dm_2^2 \cdots dM_2^N \Cr
\qQuad (N=1, \cDotS, \Infty).  \eNd{aligNed}$$ UnFoRtuNATely oNLY aN anAlYtical exprESSioN for tHe CondItional One-point probabilitY DisTribution funcTioN, ProB($m_2^I$, $z_2 |m_1$, $z_1$), iS KnOWn bAsED
 other galaxies even if it  was an is olate d s yst em  atbirt h. The predict i onsin the one-zone modelthere fo r e ma y  b e sig nifican t ly d iff er en t f ro m  w hat h app ened to  galaxiesinahierarchical un iverse.

W hit e and Frenk(19 91) de ve lop e d a d eta iledanalyt i c form alism tode s cribet he form a t io n an d evolution of ga l ax i es while takin g acco un t  o f the  hi erarchical  m ergin g  of dar k -m a t t erh alos, gas coo ling, starf orm ation,  a nds uperno va fe ed b ack . Subsequen t nu merical a pproac h es in m o delinghierar chi cal  mer g in gofda r k h a lo s e m plo y two so me wh at di ffer e n t algo rit hms;  oneis called the  ‘b lock mod el’ i n whi ch a  r andom -Gauss ian d en sity fluctuatio n fi eld is ge ner at edby  divi d ing ahyp oth eticalrectang u lar  b o x re cursively ([@CK88] ;[ @ Co le91]; [ @Cole9 4 ]) .W hile thi salg orit h m  is s impl e  a nd strai ghtfor w ar d,  the re su ltingha lomas ses a r e ne cessar ily binn ed in discrete steps of a factor o f  t w o .T he o the r generates  a r e aliz atio n  o f h a lo me rgertr e es according to a prob ab ilitydistr ibution funct ion predic t e d  by theexte n de d  Press–Schecht er th eory ([@Bo w er91]; [ @Bond 91]; [@K W93]; [@S K 9 9]; [@SL 99] ).The  la t t er  is widely us e d  inst udyingthe  cosmol ogi cal  ev olu ti on of gal axies in  a  h ie ra rch icalu niverse([ @KW G9 3];  [@ba u gh98];  [@SP 99];  [ @C o le0 0]; [@n a ga s h ima] ).  T hrou gho ut  thepres e ntpaper,we call t hel atte rme thod th e Monte-Carlo  m odeling of  m erg er his t o ries (si mply, the Monte-Carlo m o deling) , w hileit i s usually  re ferred  to as a s emi-an alyti cmod e l  of g a l ax y f or mation (SA M ) .
The m os t im portant  ingredient in Mon t e-C arlo modeling  is  the c on dit i on a l j oi n t-p r o bability distri bution fun ct i on  of a seto f * pr ogenito r* halo s ofm ass $M_ 2^{j}$ at  a redshi ft  of$ z _2$ , which is  a partof a *par e nt* h a lo  of m ass  $M_1$  a t $ z_1$,  conce p tua lly w ritten  a s $$\b egin{ al igned}
\ label{eq:jointprob}
{\r m Prob }(M_2 ^1,  M_2^2, \ cdo t s,M_2^N, z_ 2 |M_1, z_1)dM_ 2^1  dM_2 ^2\ cdots  dM_ 2 ^N  \c r 
\qqu ad ( N =1, \cdot s ,\in f t y) .  \end{ali g n e d}$ $ Unf ort u nately  onl y an analytical e x pression for t he c o n dit ion a l on e- point probabil ity  d i s tributio nfunction, P rob($M_2 ^i $ , $z_ 2 |M_1 $, $z_ 1$), is k no w n base d
 other_galaxies even_if it was an_isolated system_at_birth. The_predictions_in the one-zone_model therefore may_be significantly different from_what happened to_galaxies_in a hierarchical universe.

White and Frenk (1991) developed a detailed analytic formalism to describe_the_formation and_evolution_of_galaxies while taking account of_the hierarchical merging of dark-matter_halos, gas_cooling, star formation, and supernova feedback. Subsequent numerical_approaches_in modeling hierarchical_merging of dark halos employ two somewhat different algorithms;_one is called the ‘block model’_in which a_random-Gaussian_density_fluctuation field is generated_by dividing a hypothetical rectangular box_recursively ([@CK88]; [@Cole91]; [@Cole94]). While this_algorithm is simple and straightforward, the resulting_halo masses are necessarily binned in_discrete steps of a factor_of two._The other generates a realization_of halo merger_trees according_to a probability_distribution function predicted by the extended_Press–Schechter theory ([@Bower91];_[@Bond91]; [@KW93]; [@SK99]; [@SL99]). The latter_is_widely used in_studying_the_cosmological evolution_of galaxies in_a_hierarchical universe_([@KWG93];_[@baugh98]; [@SP99]; [@Cole00]; [@nagashima]). Throughout the_present_paper, we call the latter method the_Monte-Carlo modeling of merger_histories_(simply, the Monte-Carlo modeling),_while it is usually referred_to as a semi-analytic model of_galaxy formation_(SAM).

The most_important ingredient in Monte-Carlo modeling is the conditional joint-probability distribution function_of a set of *progenitor* halos_of mass $M_2^{j}$ at_a redshift_of_$z_2$, which is_a_part of_a *parent* halo of mass $M_1$ at_$z_1$, conceptually_written as $$\begin{aligned}
\label{eq:jointprob}
{\rm Prob}(M_2^1, M_2^2, \cdots,_M_2^N, z_2 |
M_1, z_1)_dM_2^1_dM_2^2 \cdots dM_2^N \cr
\qquad (N=1, \cdots,_\infty).  \end{aligned}$$ Unfortunately only an_analytical expression for the conditional_one-point_probability_distribution function, Prob($M_2^i$, $z_2 |M_1$,_$z_1$), is known based