thermochemical models have been used to describe alkali @xcite , titanium and vanadium @xcite , carbon , nitrogen , and oxygen @xcite , sulfur and phosphorus @xcite , condensate ( e.g. , * ? ? ?
* ; * ? ? ?
* ) chemistry in the atmospheres of gas giant planets , brown dwarfs , and low - mass dwarf stars .
here we continue and extend these previous studies by using thermochemical equilibrium calculations to model the chemical behavior of fe , mg , and si in substellar objects .
iron , magnesium , and silicon are the most abundant rock - forming elements in a solar composition gas , and  condensed as iron metal ( fe ) , and forsterite ( mg@xmath0sio@xmath1 ) and enstatite ( mgsio@xmath2 )  will produce the most massive cloud layers in substellar atmospheres .
cloud formation strongly affects the optical and infrared spectra of substellar objects by removing gases from the overlying atmosphere and by introducing solid or liquid cloud particles ( e.g. , * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
the large and growing number of discovered brown dwarfs ( @xmath3750 ) and extrasolar planets ( @xmath3400 ) makes it impractical to model the thermochemistry of all objects individually @xcite , as was done for jupiter and saturn @xcite and gliese 229b @xcite .
instead , we adopt an approach similar to that of our previous papers @xcite and determine the abundance of each chemical species as a function of pressure , temperature , and metallicity , plotted in abundance contour diagrams .
our results are thus independent of any particular pressure - temperature profile , and in principle , the atmospheric profile for any object may be superimposed on the abundance diagrams to determine its equilibrium atmospheric chemistry . in some instances , the behavior of key gases may be diagnostic of atmospheric temperature and/or tracers of weather in substellar objects .
we begin with a brief description of our computational method ( @xmath4[s computational method ] ) , and then present our results for iron chemistry in substellar atmospheres in @xmath4[s iron chemistry ] .
we first give an overview of iron chemistry in a solar composition gas and identify important gases and condensates ( @xmath4[ss overview of iron chemistry ] ) .
this is followed by more detailed discussion of the chemical behavior individual fe - bearing gases as a function of temperature , pressure , and metallicity ( @xmath4[ss chemical behavior of iron - bearing gases ] ) . wherever possible
, we note relevant spectroscopic observations of fe - bearing gases in substellar objects . our results for iron
are then summarized by illustrating fe gas chemistry along the atmospheric profiles of representative substellar objects ( @xmath4[ss iron chemistry in substellar objects ] ) . a similar approach to magnesium and
silicon chemistry follows in @xmath4[s magnesium chemistry ] and @xmath4[s silicon chemistry ] , respectively .
we conclude with a brief summary in @xmath4[s summary ] .
thermochemical equilibrium calculations were performed using a gibbs free energy minimization code , previously used for modeling the atmospheric chemistry of saturn @xcite and sulfur and phosphorus chemistry in substellar objects @xcite .
thermodynamic data for the equilibrium calculations were taken from the compilations of @xcite , @xcite , the fourth edition of the janaf tables @xcite , and the thermodynamic database maintained in the planetary chemistry laboratory @xcite .
this database includes additional thermodynamic data from the literature for compounds absent from the other compilations .
all calculations were conducted using elemental abundances from @xcite for a solar system ( i.e. , protosolar ) composition gas .
the effect of metallicity on sulfur and phosphorus chemistry was examined by running computations at [ fe / h ] = -0.5 dex ( subsolar ) , [ fe / h ] = 0 dex ( solar ) , and [ fe / h ] = + 0.5 dex ( enhanced ) metallicities .
the metallicity factor , @xmath5 , is defined as @xmath6 $ ] .
we assume that the elemental abundance ratios for mg , si , and other elements of interest vary similarly with [ fe / h ] ( e.g. , [ mg / h ] @xmath7 [ si / h ] @xmath7 [ fe / h ] ) over the range of metallicities considered here ( see * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
when considering the chemical behavior of individual gases , we focus on higher temperatures ( 800 k and higher ) , where thermochemical processes are expected to dominate over disequilibrium processes such as photochemistry or atmospheric mixing ( e.g. , see * ? ? ?
we assume that condensates settle gravitationally into a cloud layer and are removed from the cooler , overlying atmosphere .
this equilibrium cloud condensate scenario for the deep atmospheres of giant planets and brown dwarfs is supported by several lines of evidence ( e.g. , see * ? ? ?
* and references therein ) .
first , the presence of germane ( geh@xmath1 ) and the absence of silane ( sih@xmath1 ) in the upper atmospheres of jupiter and saturn ( even though si is expected to be much more abundant than ge ) can be explained by the removal of si from the gas into silicate clouds deeper in the atmosphere , whereas ge remains in the gas phase @xcite .
secondly , the detection of h@xmath0s in jupiter s troposphere by the _ galileo _
entry probe indicates that fe must be sequestered into a cloud layer at deep atmospheric levels , because the formation of fes would otherwise remove h@xmath0s from the gas above the @xmath8 k level @xcite .
third , absorption from monatomic k gas in the spectra of t dwarfs @xcite requires the removal of al and si at deeper atmospheric levels , because k would otherwise be removed from the observable atmosphere by the condensation of orthoclase ( kalsi@xmath2o@xmath9 ) @xcite .
the presence of monatomic na gas in brown dwarfs @xcite also suggests al and si removal , because albite ( naalsi@xmath2o@xmath9 ) condensation would otherwise effectively remove na from the observable atmosphere .
furthermore , the removal of na by na@xmath0s cloud formation is consistent with the observed weakening of na atomic lines throughout the l dwarf spectral sequence and their disappearance in early t dwarfs ( e.g. , * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
finally , as we note below , the disappearance of iron , magnesium , and silicon spectral features in later spectral types is consistent with removal of these elements into cloud layers . in our thermochemical model , the abundances of fe- ,
mg- , and si - bearing gases above the clouds are computed assuming saturation ( equilibrium ) vapor pressure .
if supersaturation occurs , a condensate will form and settle toward the cloud layer to restore equilibrium .
figure [ figure iron chemistry overview ] illustrates model atmospheric profiles for an m dwarf ( @xmath10 = 2600 k , @xmath11 ; * ? ? ?
* ) , an l dwarf ( @xmath10 = 1800 k , @xmath11 ; * ? ? ?
* ) , the hot , close - in ( pegasi ) planet hd209458b ( @xmath10 = 1350 k ; * ? ? ?
* ) , the t dwarf gliese 229b ( @xmath10 = 960 k ; * ? ? ?
* ) , and jupiter ( @xmath10 = 124 k ) , indicated by dashed lines .
we note that jovian atmospheric chemistry differs slightly than that for a solar - metallicity gas because jupiter has a heavy element enrichment comparable to [ fe / h ] @xmath7 + 0.5 dex @xcite .
also shown in figure [ figure iron chemistry overview ] are lines indicating where a(ch@xmath1 ) = a(co ) and a(h@xmath0 ) = a(h ) .
these boundaries are important because carbon and oxygen affect the chemical behavior of many fe , mg , and si - bearing gases .
methane is the dominant carbon - bearing gas in jupiter and t dwarfs ( such as gliese 229b ) whereas co is the dominant carbon - bearing gas in l dwarfs and pegasi planets ( such as hd209458b ) .
molecular hydrogen dissociates into monatomic h at high temperatures and low pressures ( lower right corner of figure [ figure iron chemistry overview ] ) .
however , h@xmath0 is the dominant form of hydrogen in substellar objects , and we therefore take @xmath12 throughout the following .
this approximation holds for metallicities up to [ fe / h ] @xmath7 + 0.5 dex ; at higher metallicities the h@xmath0 mole fraction abundance decreases as the relative abundance of heavy elements increases ( e.g. , @xmath13 at [ fe / h ] @xmath7 + 1.0 dex ) .
the dotted line in figure [ figure iron chemistry overview ] shows the condensation curve for fe metal , with an open circle denoting its normal melting point ( 1809 k ) .
iron condensation occurs via @xmath14 the equilibrium condensation temperature ( @xmath15 ) of fe is approximated by @xmath16,\ ] ] where higher pressures and/or metallicities lead to higher fe condensation temperatures .
iron cloud formation effectively removes nearly all iron from the atmosphere , and the abundances of fe - bearing gases above the clouds rapidly decrease with altitude .
monatomic fe gas is the dominant fe - bearing gas in a protosolar composition gas until it is replaced by fe(oh)@xmath0 at low temperatures and high pressures .
the conversion between fe and fe(oh)@xmath0 is represented by the net thermochemical reaction @xmath17 and the solid line dividing the fe and fe(oh)@xmath0 fields in figure [ figure iron chemistry overview ] indicates where these gases have equal abundances a(fe ) = a(fe(oh)@xmath0 ) @xmath18fe@xmath19 , where @xmath20fe@xmath19 is the total amount of iron in the gas .
the position of this line is given by @xmath21,\ ] ] showing that an increase in metallicity will shift the fe - fe(oh)@xmath0 boundary to higher temperatures and lower pressures .
in other words , an increase in metallicity increases the stability field of fe(oh)@xmath0 for otherwise constant conditions .
the mole fraction abundance of monatomic fe as a function of pressure and temperature is shown in figure [ figure iron gas 1]a .
as mentioned above , monatomic fe is the dominant fe - bearing gas over a wide range of pressures and temperatures in a solar system composition gas .
below the fe clouds , the abundance of fe gas is representative of the total iron content of the atmosphere ( @xmath22 ) , and its abundance is given by @xmath23.\ ] ] upon fe metal condensation , the amount of iron in the gas rapidly decreases and the equilibrium fe gas abundance is governed by its vapor pressure over solid or liquid iron , represented in reaction ( [ reaction fe condensation ] ) .
the mole fraction abundance of monatomic fe above the clouds is given by @xmath24 inversely proportional to @xmath25 .
the fe gas abundance is independent of metallicity in this region because it depends solely on the temperature - dependent vapor pressure of iron . as described below ( see @xmath4[sss iron hydride , feh ] and @xmath4[sss iron monohydroxide , feoh ] ) , this expression can be used with chemical equilibria to determine the abundances of other iron gases in substellar atmospheres .
equations giving the abundance of fe ( and other fe - bearing gases ) as a function of temperature , pressure , and metallicity below and above the fe clouds are listed in table [ table iron reactions ] .
neutral monatomic fe possesses several hundred spectral lines in the @xmath26 and @xmath27 bands @xcite , and fe i features are observed in the spectra of brown dwarfs and low - mass dwarf stars ( e.g. * ? ? ?
the strong fe feature at 1.189 @xmath28 m weakens in mid- to late - type m dwarf spectra and generally disappears in mid - type l dwarfs @xcite .
this trend is consistent with the removal of iron from the gas into an fe metal cloud deck located deeper below the photosphere in objects with low effective temperatures ( e.g. * ? ? ?
the chemical behavior of feh in a protosolar composition gas is illustrated in figure [ figure iron gas 1]b .
the conversion between fe and feh occurs via the net thermochemical reaction @xmath29 here we show how chemical equilibria may used to derive equations giving the abundance of fe - bearing species as a function of temperature , pressure , and metallicity . rearranging the equilibrium constant expression for reaction ( [ reaction fe : feh ] ) yields @xmath30 using the fe abundance from equation ( [ equation fe below clouds ] ) , the temperature dependence of @xmath31 ( @xmath32 from 800 to 2500 k ) , and the hydrogen abundance ( @xmath33 ) , the feh abundance between the h@xmath0-h boundary and the fe cloud deck is given by @xmath34,\ ] ] proportional to @xmath35 and @xmath5 . at high temperatures ( @xmath36 k ) and low pressures ( @xmath37 bar )
as monatomic h becomes increasingly abundant the h@xmath0 abundance begins to decrease near the h@xmath0-h boundary .
the reduced h@xmath0 abundance , in turn , reduces feh formation via reaction ( [ reaction fe : feh ] ) and changes the chemical behavior of feh at high temperatures and low pressures ( i.e. , lower right corner of figure [ figure iron gas 1]a ) .    above the iron clouds ,
the fe abundance in equation ( [ equation fe above clouds ] ) is used in equation ( [ equation feh abundance generic ] ) to give @xmath38 where @xmath39 is proportional to @xmath40 and independent of metallicity .
the differences in chemical behavior of feh below ( equation [ equation feh abundance below clouds ] ) and above ( equation [ equation feh abundance above clouds ] ) the fe cloud deck are illustrated in the shape of the feh abundance contours in figure [ figure iron gas 1]b .
the sharp bends in the contours correspond to the fe condensation curve in figure [ figure iron chemistry overview ] .
absorption bands from feh are common in the near - infrared spectra of brown dwarfs ( e.g. * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
the most prominent feh feature is the band located at 0.9896 @xmath28 m , which weakens throughout the l dwarf spectral sequence and in early t dwarfs @xcite , consistent with the removal of iron into a cloud located deeper and deeper in the atmosphere .
this feature unexpectedly strengthens again in mid - type t dwarf spectra @xcite , prompting different explanations for the shape of the color - magnitude diagram for brown dwarfs near the l - t transition .
@xcite suggested that the strengthening feh bands are caused by upward convective mixing of feh gas from deeper levels where it is more abundant .
however , as pointed out by @xcite and demonstrated @xcite , the fragile fe@xmath41h bond is unlikely to survive convective upwelling in a t dwarf atmosphere . instead ,
the observations are plausibly explained by cloud disruption and clearing ( in a @xmath421 @xmath28 m window ) which allows the observation of feh gas at deep atmospheric levels @xcite .
mole fraction contours of feoh are illustrated in figure [ figure iron gas 1]c .
the equilibrium between fe and feoh is represented by the net thermochemical reaction @xmath43 and expressions giving the feoh abundance as a function of temperature , pressure , and metallicity are listed in table [ table iron reactions ] . above the iron clouds , the curvature in the feoh abundance contours along the ch@xmath1-co boundary results from the effect of carbon chemistry on the h@xmath0o abundance in reaction ( [ reaction fe : feoh ] ) ( e.g. * ? ? ?
for example , the atmospheric water abundance may be written as @xmath44,\ ] ] where @xmath45 is the water abundance in a solar - metallicity gas , which is @xmath46 inside the ch@xmath1 field and @xmath47 inside the co field in figure [ figure iron chemistry overview ] . at temperatures and pressures near the ch@xmath1-co boundary ,
@xmath45 may be derived from ch@xmath1-co equilibria @xcite .
rearranging the equilibrium constant expression for reaction ( [ reaction fe : feoh ] ) gives @xmath48 substituting for the fe abundance from equation ( [ equation fe above clouds ] ) , the h@xmath0o abundance from equation ( [ equation h2o abundance generic ] ) , the h@xmath0 abundance ( @xmath49 ) , and the temperature dependence of @xmath50 ( @xmath51 from 800 to 2500 k ) , equation ( [ equation feoh abundance generic ] ) becomes @xmath52 + \log x_{\textrm{\scriptsize{h}}_{2}\textrm{\scriptsize{o}}}'.\ ] ] this expression gives the feoh abundance above the iron clouds and includes the effect of metallicity on the atmospheric water abundance .
the kinks in the feoh contours denote the position of the fe condensation curve .
as demonstrated in equation ( [ equation feoh abundance above clouds ] ) and shown in figure [ figure iron gas 1]c , the mole fraction abundance of feoh is proportional to @xmath40 throughout this region .
the chemical behavior of fe(oh)@xmath0 as a function of temperature and pressure is illustrated in figure [ figure iron gas 1]d .
the equilibrium conversion between fe and fe(oh)@xmath0 is represented by reaction ( [ reaction fe : fe(oh)2 ] ) .
below the fe clouds ( kinks in abundance contours ) , the fe(oh)@xmath0 abundance increases with total pressure , and has a strong @xmath53 dependence on metallicity ( see table [ table iron reactions ] ) .
the fe(oh)@xmath0 abundance in reaction ( [ reaction fe : fe(oh)2 ] ) is very sensitive to the h@xmath0o abundance since the formation of one mole of iron hydroxide requires two moles of water . as a result
, there is a large shift in the fe(oh)@xmath0 abundance contours above the clouds ( inflections in figure [ figure iron gas 1]d ) when moving between the ch@xmath1 and co fields . within each field ,
the fe(oh)@xmath0 abundance above the clouds is essentially pressure - independent ( see table[table iron reactions ] ) and is therefore , in principle , diagnostic of atmospheric temperature .    as shown in figure [ figure iron chemistry overview ]
, fe(oh)@xmath0 becomes the most abundant fe - bearing gas at low temperatures and high pressures in a solar composition gas .
iron hydroxide is therefore expected to be the most abundant iron gas at @xmath54 on jupiter ( @xmath55\approx+0.5 $ ] ) and @xmath56 on gliese 229b .
however , we emphasize that the amount of iron remaining in the gas is greatly diminished by the condensation of fe metal at deeper atmospheric levels .
for example , even where it is the most abundant fe - bearing gas , the predicted fe(oh)@xmath0 mole fraction is @xmath57 at the 1000 k level in jupiter s atmosphere .
mole fraction abundance contours for feo are shown in figure [ figure iron gas 2]a for a solar composition gas .
monatomic iron reacts with water to form feo via the net reaction @xmath58 between the h@xmath0-h boundary and the fe cloud deck , @xmath59 is proportional to @xmath60 and is independent of pressure ( see table [ table iron reactions ] ) . in principle , the feo abundance in this region is therefore diagnostic of either temperature or metallicity . above the clouds , the feo abundance decreases with pressure and increases with metallicity . as for feoh and fe(oh)@xmath0 ,
the feo abundance is sensitive to the h@xmath0o abundance and the feo abundance contours display a subtle shift along the ch@xmath1-co equal abundance boundary .
the abundance of fes gas as a function of temperature and pressure is illustrated in figure [ figure iron gas 2]b .
the chemical behavior of fes is governed by the reaction between fe and h@xmath0s : @xmath61 above the iron clouds , h@xmath0s is the dominant sulfur - bearing gas ( @xmath62 $ ] ) , and the fes abundance is proportional to @xmath63 and @xmath5 ( see table [ table iron reactions ] ) .
below the fe cloud deck , the fes abundance is @xmath64 ppb and decreases at higher temperatures and lower pressures as h@xmath0s is replaced by sh and s @xcite .
iron sulfide is predicted to be the second or third most abundant iron gas throughout l dwarf atmospheres and in the upper atmospheres ( @xmath65 k ) of pegasi planets .
figures [ figure iron chemistry summary]a-[figure iron chemistry summary]d summarize the iron gas chemistry in four representative substellar objects : jupiter , the t dwarf gliese 229b , the pegasi planet hd209458b , and an l dwarf ( @xmath66 k ) . iron chemistry in substellar objects
is strongly affected by fe metal condensation at deep atmospheric levels .
monatomic fe is the dominant fe - bearing gas in the deep atmospheres of jupiter and gliese 229b , and throughout the atmospheres of hd209458b and the l dwarf .
a number of fe - bearing gases become relatively more abundant in objects with lower effective temperatures . on jupiter and gliese 229b , fe(oh)@xmath0 replaces fe at lower temperatures .
the fe(oh)@xmath0 and feo abundances are pressure - independent and thus potentially diagnostic of atmospheric temperature , respectively , above and below the iron clouds . for all four objects ,
iron hydride ( feh ) is the second most abundant fe - bearing gas at deep atmospheric levels , until it is replaced at lower temperatures by fe(oh)@xmath0 , feoh , or fes .
because of its strong absorption and relatively high abundance near the iron condensation level , feh is a tracer of weather in brown dwarfs ( e.g. * ? ? ?
other fe - bearing gases ( e.g. fe , feoh , fes ) are potential tracers of weather since they typically achieve their maximum abundance near iron cloud base .
figure [ figure magnesium chemistry overview ] gives an overview of mg chemistry as a function of pressure and temperature in a protosolar composition gas .
the dotted lines labeled mg@xmath0sio@xmath1(s , l ) , mgsio@xmath2(s , l ) , and mgo(s ) show the equilibrium condensation curves for forsterite , enstatite , and periclase , respectively , and the open circles denote the normal melting temperatures for forsterite ( 2163 k ) and enstatite ( 1830 k ) .
forsterite ( mg@xmath0sio@xmath1 ) condenses via the net thermochemical reaction @xmath67 and its condensation temperature as a function of pressure and metallicity is approximated by @xmath68,\ ] ] at slightly lower temperatures , enstatite ( mgsio@xmath2 ) condensation occurs via the net reaction @xmath69 the enstatite condensation curve is approximated by @xmath70,\ ] ] where higher pressures and/or metallicities lead to higher condensation temperatures . at very high pressures , forsterite and enstatite condensation temperature
are depressed as sio is replaced by sih@xmath1 ( see @xmath4[ss overview of silicon chemistry ] ) .
periclase ( mgo ) condenses via the net thermochemical reaction @xmath71 at pressures greater than @xmath72 bar .
the condensation curve for akermanite ( ca@xmath0mgsi@xmath0o@xmath73 ) is not shown here because most ca is expected to be removed at deeper atmospheric levels by the condensation of refractory calcium aluminates ( e.g. , * ? ? ?
if no calcium is removed , ca@xmath0mgsi@xmath0o@xmath73 condensation would consume @xmath74 of the total atmospheric mg inventory .
in the same way , olivine ( ( mg , fe)@xmath0sio@xmath1 ) and fayalite ( fe@xmath0sio@xmath1 ) are not expected in substellar atmospheres because nearly all fe is removed from the gas phase by iron metal condensation at higher temperatures ( e.g. , * ? ? ? * ; * ? ? ? * )    magnesium - silicate cloud formation is very effective at removing nearly all ( @xmath399% ) magnesium from the atmosphere , and the abundances of mg - bearing gases rapidly decrease with altitude above the clouds .
this behavior is generally consistent with the disappearance of mg spectral features by early - type l dwarfs ( see  [ sss monatomic magnesium , mg ] ) .
furthermore , @xcite find a si@xmath41o absorption feature at 10 @xmath28 m in mid - type l dwarfs which is consistent with the presence of silicate grains and a weak 9.17 @xmath28 m feature tentatively attributed to crystalline enstatite .
monatomic mg is the dominant mg - bearing gas in substellar atmospheres until it is replaced by mg(oh)@xmath0 at low temperatures and high pressures .
the conversion between mg and mg(oh)@xmath0 is represented by the net thermochemical reaction @xmath75 and the solid line in figure [ figure magnesium chemistry overview ] indicates where mg and mg(oh)@xmath0 have equal abundances
a(mg ) = a(mg(oh)@xmath0 ) @xmath18mg@xmath19 , where @xmath20mg@xmath19 is the total amount of magnesium in the gas .
the position of the equal - abundance line is given by @xmath76,\ ] ] where an increase in metallicity shifts the mg - mg(oh)@xmath0 boundary to higher temperatures and lower pressures .
figure [ figure magnesium gas 1]a shows the chemical behavior of monatomic mg as a function of pressure and temperature in a solar - metallicity gas .
the chemical behavior of mg is strongly affected by silicate cloud formation . below the forsterite clouds ,
the abundance of mg gas is given by @xmath77,\ ] ] and comprises nearly 100% of the total elemental mg content in the atmosphere ( @xmath78 ) .
upon forsterite and enstatite condensation , the magnesium gas abundance is governed by its vapor pressure over rock , represented by reactions ( [ reaction forsterite condensation ] ) and ( [ reaction enstatite condensation ] ) .
curvature in the mg mole fraction contours in figure [ figure magnesium gas 1]a occurs along the ch@xmath1-co boundary , which affects the h@xmath0o abundance in reactions ( [ reaction forsterite condensation ] ) and ( [ reaction enstatite condensation ] ) .
the mg abundance is therefore @xmath79.\ ] ] in ch@xmath1-dominated objects , and @xmath80.\ ] ] in co - dominated objects .
we can use these expressions along with chemical equilibria to determine the abundances of other mg - bearing gases in substellar atmospheres .
equations giving the abundance of mg - bearing gases as a function of pressure , temperature , and metallicity in a protosolar composition gas are listed in table [ table magnesium reactions ] .
several mg absorption lines are present in the near infrared spectra of m dwarfs , including the prominent feature at 1.183 @xmath28 m @xcite .
these features weaken in mid- to late - type m dwarfs and generally disappear by @xmath42 l1 @xcite , consistent with the removal of magnesium into mg@xmath0sio@xmath1 and mgsio@xmath2 clouds located at increasingly greater depths below the observable atmosphere ( e.g. , see * ? ? ?
* ; * ? ? ?
mole fraction contours for mgh are illustrated in figure [ figure magnesium gas 1]b .
the mgh abundance is governed by equilibrium with monatomic mg , via the net thermochemical reaction @xmath81 between the h@xmath0-h boundary and the mg - silicate cloud base , the mgh abundance is proportional to @xmath35 and @xmath82 ( see table [ table magnesium reactions ] ) .
in contrast , the mgh abundance above the magnesium - silicate clouds is proportional to @xmath40 and @xmath83 .
the sharp bends in the mgh abundance contours occur at the condensation temperature of forsterite .
slight inflections in the mgh abundance occur along the ch@xmath1-co boundary , because the mg abundance in reaction ( [ reaction mg : mgh ] ) is sensitive to the h@xmath0o abundance in reactions ( [ reaction forsterite condensation ] ) and ( [ reaction enstatite condensation ] ) .
magnesium hydride is an important opacity source from 0.44 to 0.56 @xmath28 m @xcite , and mgh bands at 0.48 and 0.52 @xmath28 m have been found in the optical spectra of bright l dwarfs and extreme subdwarfs @xcite .      the chemical behavior of mgoh as a function of pressure and temperature is shown in figure [ figure magnesium gas 1]c .
the equilibrium abundance of mgoh is governed by the reaction @xmath84 decreasing with temperature and proportional to @xmath35 and @xmath60 below the mg - silicate clouds .
the kinks in the mgoh contours occur where forsterite condenses . above the clouds
, @xmath85 is proportional to @xmath40 ( see table [ table magnesium reactions ] ) .
interestingly , the mgoh abundance in this region is independent of metallicity and shows no shift at the ch@xmath1-co boundary because these effects cancel out in reaction ( [ reaction mg : mgoh ] ) .
for example , the h@xmath0o abundance slightly decreases when moving from the ch@xmath1 to the co field ( e.g. * ? ? ?
however , as can be seen by comparing equations ( [ equation mg above clouds ch4 ] ) and ( [ equation mg above clouds co ] ) , there is a corresponding increase in the mg abundance , so the resulting mgoh abundance in reaction ( [ reaction mg : mgoh ] ) remains unaffected by carbon chemistry .
magnesium monohydroxide is typically the third most abundant mg - bearing gas in substellar atmospheres .
mole fraction abundance contours for mg(oh)@xmath0 are shown in figure [ figure magnesium gas 1]d for a solar - metallicity gas . the equilibrium conversion between mg and mg(oh)@xmath0
is represented by reaction ( [ reaction mg : mg(oh)2 ] ) .
below the mg - silicate clouds , the mg(oh)@xmath0 abundance decreases with temperature and pressure , and has a very strong ( @xmath86 ) dependence on metallicity .    above the magnesium - silicate clouds ( the kinks in figure [ figure magnesium gas 1]d ) , the mg(oh)@xmath0 abundance contours show inflections along the ch@xmath1-co boundary , because mg(oh)@xmath0 is sensitive to the water abundance in reaction ( [ reaction mg : mg(oh)2 ] ) . within either the ch@xmath1 or co field
, the mg(oh)@xmath0 abundance is pressure - independent . in principle , the abundance of mg(oh)@xmath0 is therefore diagnostic of atmospheric temperature if the object s metallicity is known .    at low temperatures and high pressures , mg(oh)@xmath0 becomes the most abundant mg - bearing gas , as illustrated in figure [ figure magnesium chemistry overview ] .
magnesium hydroxide is thus expected to be the dominant magnesium gas below @xmath421550 k in the atmosphere of jupiter ( cf .
* ) and below @xmath42980 k on gliese 229b .
the abundance of mgo as a function of pressure and temperature is illustrated in figure [ figure magnesium gas 2]a .
magnesium monoxide forms via the net thermochemical reaction @xmath87 below the mg - silicate cloud deck , at pressures less than @xmath64 bar , the mgo abundance is effectively pressure - independent and is therefore potentially diagnostic of atmospheric temperature if the metallicity is known .    for similar reasons as for mgoh ( see  [ sss magnesium monohydroxide , mgoh ] ) ,
the mgo abundance above the magnesium - silicate clouds is independent of metallicity and is unaffected by the ch@xmath1-co boundary .
for example , the mg abundance decreases with metallicity , as shown by equations ( [ equation mg above clouds ch4 ] ) and ( [ equation mg above clouds co ] ) , whereas the h@xmath0o abundance increases with metallicity .
these effects cancel out in reaction ( [ reaction mg : mgo ] ) , with the result that mgo is unaffected by changes in metallicity or the prevailing carbon chemistry .
figure [ figure magnesium gas 2]b displays mole fraction contours for mgs in protosolar composition gas .
the equilibrium abundance of mgs is governed by the reaction @xmath88 near the silicate cloud base , the mgs abundance is @xmath64 ppb in a solar - metallicity gas and decreases at higher temperatures and lower pressures as h@xmath0s is replaced by monatomic s. above the mg - silicate clouds , @xmath89 is proportional to @xmath63 and is metallicity - independent ( see table [ table magnesium reactions ] ) .
magnesium sulfide is typically among the more abundant mg - bearing gases in brown dwarfs and extrasolar giant planets , and becomes the second most abundant magnesium gas at temperatures below @xmath421600 k in the atmospheres of l dwarfs and pegasi planets .
magnesium gas chemistry along the pressure - temperature profiles of jupiter , gliese 229b , hd209458b , and an l dwarf ( @xmath66 k ) are illustrated in figures [ figure magnesium chemistry summary]a-[figure magnesium chemistry summary]d .
monatomic mg is the dominant mg - bearing gas throughout the atmospheres of pegasi planets and l dwarfs , and in the deep atmospheres of giant planets and t dwarfs .
furthermore , the mg abundance below the mg - silicate cloud deck is essentially constant and representative of the total mg abundance in the atmosphere . upon condensation
, the abundances of magnesium - bearing gases rapidly decrease with decreasing temperature above the magnesium - silicate clouds . in objects with lower effective temperatures ,
a number of other mg - bearing gases become relatively abundant and mg(oh)@xmath0 replaces mg as the most abundant magnesium gas at the @xmath90 k level on jupiter and the @xmath91 k level on gliese 229b .
magnesium hydride ( mgh ) is the second most abundant magnesium gas in the deep atmospheres of substellar objects until it is replaced at lower temperatures by mg(oh)@xmath0 and mgoh ( in giant planets and t dwarfs ) or mgs ( in l dwarfs and pegasi planets ) .
magnesium hydroxide ( mg(oh)@xmath0 ) and mgo are potential atmospheric temperature probes , respectively , above and below the magnesium - silicate clouds .
an overview of silicon chemistry as a function of pressure and temperature is illustrated in figure [ figure silicon chemistry overview ] .
the dotted lines labeled mg@xmath0sio@xmath1(s , l ) and mgsio@xmath2(s , l ) show the condensation temperatures of forsterite and enstatite , which together remove nearly all silicon from the overlying atmosphere .
@xcite found that silicate absorption features near 10 @xmath28 m in mid - type l dwarf spectra are consistent with the presence of these magnesium silicates , but noted the possibility of additional absorption by quartz ( sio@xmath0 ) grains based upon the predictions of @xcite . in contrast to the models of @xcite , we find that quartz will not condense in the atmospheres of substellar objects unless enstatite condensation is suppressed .
this is demonstrated in figure 10 , which shows the elemental distribution of si in condensed phases at 1 bar total pressure in a solar - metallicity gas with ( figure [ figure enstatite]a ) and without ( figure [ figure enstatite]b ) enstatite condensation .
as shown in figure [ figure enstatite ] , mg@xmath0sio@xmath1 formation consumes nearly half of the total si abundance because the solar elemental abundances of mg and si are approximately equal .
enstatite formation plausibly proceeds via reactions between sio gas and pre - existing forsterite grains and continues until nearly all silicon is consumed .
thus , quartz ( @xmath92 k ) can only form in the absence of enstatite ( @xmath93 k ) , because mgsio@xmath2 otherwise efficiently removes silicon from the gas phase . even in the absence of gas - grain reactions between sio and mg@xmath0sio@xmath1 ,
the vapor pressures of mg and sio above forsterite ( reaction [ reaction forsterite condensation ] ) remain high enough to drive enstatite cloud formation via the net thermochemical reaction ( [ reaction enstatite condensation ] ) , so that mgsio@xmath2 condenses instead of sio@xmath0 .
we therefore conclude that sio@xmath0 will not condense within the silicate cloud .    the most abundant si - bearing gas over a wide range of pressures and temperatures is sio , until it is replaced at higher pressures by silane , sih@xmath1 ( see figure [ figure silicon chemistry overview ] ) .
the equilibrium conversion between sio and sih@xmath1 is @xmath94 the position of the line where sio and sih@xmath1 have equal abundances a(sio ) = a(sih@xmath1 ) @xmath95 ( where @xmath96 is the total amount of silicon in the gas ) is given by @xmath97.\ ] ] as temperatures decrease , sih@xmath1 is replaced by sih@xmath2f via the net thermochemical reaction @xmath98 this reaction is independent of pressure and the position of the sih@xmath1-sih@xmath2f boundary is given by @xmath99),\ ] ] and occurs at @xmath100 960 k in a solar - metallicity gas .
meanwhile , the equilibrium between sio and sih@xmath2f is @xmath101 where the position of the sio - sih@xmath2f line is given by @xmath102 independent of the metallicity .
the sih@xmath1-sio , sih@xmath1-sih@xmath2f , and sio - sih@xmath2f equal abundance lines intersect to form a `` triple point '' at @xmath100 960 k and @xmath103 bar in solar - metallicity gas , where all three gases have equal abundances [ a(sio ) = a(sih@xmath1 ) = a(sih@xmath2f ) @xmath104 .    at lower temperatures
, sih@xmath0f@xmath0 replaces sih@xmath2f via the net thermochemical reaction @xmath105 this reaction is also independent of pressure and the sih@xmath2f - sih@xmath0f@xmath0 boundary is given by @xmath106),\ ] ] and is located at @xmath100 917 k in a solar - metallicity gas .
the conversion between sio and sih@xmath0f@xmath0 takes place by the reaction @xmath107 where the position of the sio - sih@xmath0f@xmath0 boundary is given by @xmath108.\ ] ] in a solar - metallicity gas , equations ( [ boundary line sio : sih3f ] ) and ( [ boundary line sio : sih2f2 ] ) intersect to form the sio - sih@xmath2f - sih@xmath0f@xmath0 `` triple point '' at @xmath109 k and @xmath110 bar , where all three gases have equal abundances ( @xmath111 for [ fe / h ] = 0 ) .
however , the abundances of sih@xmath2f and sih@xmath0f@xmath0 in this region are extremely low because most silicon is removed from the atmosphere by cloud formation at deeper levels .        the chemical behavior of sio in a protosolar composition gas is illustrated in figure [ figure sio ] as a function of pressure and temperature . within the sio field ,
the sio abundance below the magnesium - silicate clouds is given by @xmath112,\ ] ] and sio contains @xmath42100% of the atmospheric silicon inventory ( @xmath113 ) . upon mg - silicate condensation , silicon is efficiently removed from the gas phase and the sio abundance rapidly decreases with decreasing temperature . above the clouds ,
the sio abundance is governed by its vapor pressure over rock .
curvature in the sio contour lines along the ch@xmath1-co equal abundance boundary results from the effect of carbon chemistry on the h@xmath0o abundance in reactions ( [ reaction forsterite condensation ] ) and ( [ reaction enstatite condensation ] ) .
the sio abundance is thus @xmath114,\ ] ] in ch@xmath1-dominated objects and @xmath115,\ ] ] in co - dominated objects , inversely proportional to pressure and metallicity .
these expressions , along with chemical equilibria , are used to determine the equilibrium abundances of other si - bearing gases .
expressions giving the abundances of silicon species as a function of temperature , pressure , and metallicity in a protosolar composition gas are listed in table [ table silicon reactions ] .
we note that the types of cloud condensates present will affect the gas chemistry of sio and subsequently all other si - bearing gases .
for example , as shown in figure [ figure sio condense ] , sio mole fraction abundances above the clouds are @xmath420.5 dex higher if enstatite formation is suppressed and replaced by sio@xmath0 condensation ( see  [ ss overview of silicon chemistry ] above ) .
the sio abundance is therefore potentially diagnostic of weather and cloud composition in brown dwarf atmospheres .
silicon monoxide has not yet been detected in the atmospheres of brown dwarfs or giant planets , but has been observed in numerous objects including molecular clouds , circumstellar envelopes , the photospheres of late - type stars , and sunspots ( e.g. , * ? ? ? * and references therein ) .
abundant sio gas was recently detected in the circumstellar disk of the @xmath116 pic analog hd172555 @xcite .
mole fraction contours for sis are shown in figure [ figure silicon gas 1]a for a solar - metallicity gas .
silicon monosulfide is formed by the reaction between sio and h@xmath0s : @xmath117 below the mg - silicate cloud deck , @xmath118 - 10 ppm in a solar - metallicity gas .
the sis abundance decreases at higher temperatures and lower pressures as h@xmath0s is replaced by sh and monatomic s ( e.g. * ? ? ?
* ) .    above the clouds ,
the sis abundance is inversely proportional to metallicity and total pressure ( see table [ table silicon reactions ] ) .
the sis abundance contours in figure [ figure silicon gas 1]a also display curvature along the ch@xmath1-co boundary .
this shift is more pronounced for sis than for sio because the sis abundance in reaction ( [ reaction sio : sis ] ) depends on the sio and h@xmath0o abundances , both of which are affected by carbon chemistry .
for example , when reaction ( [ reaction sio : sis ] ) is at equilibrium , lechtelier s principle shows that more sis is produced either by adding sio ( or h@xmath0s ) or removing h@xmath0o . when moving from the ch@xmath1 field to the co field in a protosolar composition gas , the sio abundance increases and the h@xmath0o abundance decreases , yielding a correspondingly large increase in the sis abundance .
silicon sulfide is expected to be the second most abundant silicon - bearing gas ( after sio ) in the atmospheres of brown dwarfs and extrasolar giant planets .
furthermore , it is a potential tracer of weather in these objects because its maximum abundance is typically achieved near the magnesium - silicate cloud base .
figure [ figure silicon gas 1]b illustrates the chemical behavior of monatomic si gas as a function of pressure and temperature .
the abundance of si is governed by the net thermochemical reaction @xmath119 below the mg - silicate cloud deck and at pressures less than @xmath64 bar , the si abundance is effectively independent of pressure and metallicity and thus potentially diagnostic of atmospheric temperature . above the clouds , the si abundance rapidly decreases with decreasing temperature and
is proportional to @xmath63 and @xmath120 ( see table [ table silicon reactions ] ) .
a number of si absorption bands are observed in the near infrared spectra of low - mass dwarf stars .
these features generally weaken and disappear in late - type m dwarfs @xcite , consistent with the removal of silicon into mg@xmath0sio@xmath1 and mgsio@xmath2 clouds deeper in the atmosphere .
the abundance of sih as a function of pressure and temperature is illustrated in figure [ figure silicon gas 2]a .
the sih abundance is governed by the net thermochemical reaction @xmath121 below the mg - silicate clouds , the sih abundance is proportional to @xmath35 and is independent of metallicity . above the clouds , the sih abundance is proportional to @xmath40 and @xmath120 ( see table [ table silicon reactions ] ) . when moving from the ch@xmath1 to the co field , the sio abundance increases and the h@xmath0o abundance decreases in reaction ( [ reaction sio : sih ] ) .
both effects serve to increase the sih abundance , resulting in a shift in the sih contour lines along the ch@xmath1-co boundary .
the @xmath120 dependence on metallicity means that for every @xmath55=+1 $ ] dex increase in metallicity , the sih abundance decreases by a factor of 100 at a given pressure and temperature .
mole fraction contours for sih@xmath0 are shown in figure [ figure silicon gas 2]b .
the equilibrium between sio and sih@xmath0 is given by the reaction @xmath122 beneath the magnesium - silicate cloud deck , the metallicity - independent silylene abundance increases with total pressure . upon rock condensation
, the amount of sih@xmath0 remaining in the gas rapidly decreases with decreasing temperature . as for sih ,
the sih@xmath0 abundance is sensitive to the sio and h@xmath0o abundances in reaction and thus a shift occurs in the sih@xmath0 contour lines when moving between the ch@xmath1 and co fields . within each field
, the sih@xmath0 abundance is pressure - independent ( see figure [ figure silicon gas 2 ] and table [ table silicon reactions ] ) , and thus potentially diagnostic of temperature if the metallicity is known and sih@xmath0 is thermochemical in origin .
the chemical behavior of the sih@xmath2 radical in a solar - metallicity gas is illustrated in figure [ figure silicon gas 2]c .
the abundance of sih@xmath2 is governed by the reaction @xmath123 between the h@xmath0-h boundary and the mg - silicate cloud base , the sih@xmath2 is proportional to @xmath124 . above the clouds , the sih@xmath2 abundance is proportional to @xmath124 and @xmath120 ( see table [ table silicon reactions ] ) .
there is curvature in the sih@xmath2 contour lines along the ch@xmath1-co boundary because the prevailing carbon chemistry affects both the sio and h@xmath0o abundances in reaction ( [ reaction sio : sih3 ] )      mole fraction abundance contours for sih@xmath1 are illustrated in figure [ figure silicon gas 2]d .
the equilibrium conversion between sio and sih@xmath1 is represented by reaction ( [ reaction sih4:sio ] ) , and expressions giving the sih@xmath1 abundance as a function of pressure , temperature , and metallicity are listed in table [ table silicon reactions ] . below the clouds , the silane abundance has a strong ( @xmath125 ) dependence on total pressure .
upon mg - silicate cloud formation , there is an inflection in the sih@xmath1 contour lines , and the sih@xmath1 abundance above the clouds is proportional to @xmath25 and @xmath120 .
curvature in the sih@xmath1 contour lines along the ch@xmath1-co boundary results from the effect of carbon chemistry on the sio and h@xmath0o abundances in reaction ( [ reaction sih4:sio ] ) .    at high pressures and temperatures ( @xmath126 k for @xmath55=0 $ ] ) ,
sih@xmath1 becomes the most abundant silicon - bearing gas .
we thus expect sih@xmath1 to be the dominant silicon gas below the 1031 k level in the deep atmosphere of jupiter ( @xmath55\approx+0.5 $ ] ) .
however , we again emphasize the efficiency with which rock condensation removes silicon from the gas , and we expect an abundance of @xmath127 at the 1031 k level on jupiter . at lower temperatures , sih@xmath1 is surpassed by sih@xmath2f and sih@xmath0f@xmath0 .
silane is expected to be the third most abundant si - bearing gas in the deep atmosphere of gliese 229b .
we also point out the trend in pressure dependence for the silicon hydrides below ( @xmath128 , @xmath129 , @xmath130 , @xmath131 ) and above the silicate clouds ( @xmath132 , @xmath133 , @xmath134 , @xmath135 ) , which is evident in the shapes of the contour lines in figure [ figure silicon gas 2 ] .
the abundance of each of the silicon hydrides is metallicity - independent below the mg - silicate clouds and proportional to @xmath120 above the clouds .
figures [ figure silicon chemistry summary]a-[figure silicon chemistry summary]d summarize the chemical behavior of silicon gases along the pressure temperature profiles of jupiter , the t dwarf gliese 229b , the pegasi planet hd209458b , and an l dwarf ( @xmath66 k ) .
there is a clear trend in silicon chemistry as a function of effective temperature . on jupiter
, sih@xmath1 is the dominant si - bearing gas throughout the deep atmosphere , and is replaced by sih@xmath2f and sih@xmath0f@xmath0 at lower temperatures ( see figure [ figure silicon chemistry overview ] ) .
the second most abundant si - bearing gas is sio , followed by sis and a number of other silicon gases . in the warmer atmosphere of gliese 229b
, sio is the dominant si - bearing gas , followed by sis and sih@xmath1 .
the relative importance of sih@xmath1 decreases with increasing effective temperature , and sio and sis are the most important silicon gases throughout the atmospheres of pegasi planets and l dwarfs .
the abundances of all the silicon gases shown in figure [ figure silicon chemistry summary ] rapidly decrease with decreasing temperature above the silicate clouds , which explains the non - detection of sih@xmath1 and other silicon species in the atmospheres of jupiter and saturn . below the cloud base ,
the important silicon gases sio and sis reach their maximum abundance and are therefore potential tracers of weather in brown dwarfs and pegasi planets .
silylene ( sih@xmath0 ) and monatomic silicon ( si ) are potentially diagnostic of atmospheric temperature , respectively , above and below the magnesium silicate clouds .
the abundant silicon gases sio and sis typically achieve their maximum abundance near the cloud base and are potential tracers of weather in brown dwarfs and pegasi planets .
the chemical behavior of iron species in substellar atmospheres is strongly affected by fe metal condensation , which efficiently removes most iron from the gas phase .
similarly , most magnesium and silicon is removed from the gas by forsterite ( mg@xmath0sio@xmath1 ) and enstatite ( mgsio@xmath2 ) cloud formation .
the equilibrium abundances of fe- , mg- , and si - bearing gases rapidly decrease with increasing altitude ( and decreasing temperature ) above the clouds .
monatomic iron is the dominant fe - bearing gas throughout the atmospheres of l dwarfs and pegasi planets .
other less abundant iron gases become increasingly important in objects with lower effective temperatures , and fe(oh)@xmath0 replaces fe at low temperatures in t dwarfs and giant planets .
magnesium gas chemistry is similar to that of iron .
monatomic mg is the most abundant magnesium gas throughout the atmospheres of l dwarfs and pegasi planets and in the deep atmospheres of giant planets and t dwarfs , where mg is replaced by mg(oh)@xmath0 at lower temperatures .
a number of mg - bearing gases become relatively abundant with decreasing effective temperature .
silicon monoxide ( sio ) is the most abundant si - bearing gas , followed by sis , throughout the atmospheres of l dwarfs and pegasi planets and in the deep atmospheres of t dwarfs . in objects with lower effective temperatures ,
a number of other silicon gases become increasingly important and sih@xmath1 is the dominant silicon gas in the deep atmosphere of jupiter . at high pressures and low temperatures sih@xmath1 and
sio are replaced by sih@xmath2f and/or sih@xmath0f@xmath0 .
the abundances of several fe- , mg- , and si - bearing gases are pressure - independent and thus , in principle , diagnostic of atmospheric temperature .
these include fe(oh)@xmath0 , mg(oh)@xmath0 , and si below the clouds and feo , mgo , and sih@xmath0 above the clouds .
in addition , a number of gases ( e.g. fe , feh , feoh , fes , mg , mgh , mgoh , mgs , sio , sis ) may serve as indicators of weather since they generally reach their maximum abundance just below the iron metal or magnesium - silicate cloud decks .
this may be particularly useful for late m dwarfs and early l dwarfs in which the metal and silicate clouds are located at relatively shallow depths below the photosphere .
this research was conducted at washington university in st .
louis and was supported by the nasa planetary atmospheres program ( nng06gc26 g ) .
support for k.  lodders was also provided by nsf grant ast-0707377 .
final preparation of the manuscript was supported by the lunar and planetary institute / usra ( nasa cooperative agreement ncc5 - 679 ) .
lpi contribution no .
xxxx .
llc fe & @xmath137 $ ] & + feh & @xmath138 $ ] & [ reaction fe : feh ] + feoh & @xmath139 $ ] & [ reaction fe : feoh ] + fe(oh)@xmath0 & @xmath140 $ ] & [ boundary line fe : fe(oh)2 ] + feo & @xmath141 $ ] & [ reaction fe : feo ] + fes & @xmath142 $ ] & [ reaction fe : fes ] + fe & @xmath143 & [ reaction fe condensation ] + feh & @xmath144 & [ reaction fe : feh ] + feoh & @xmath145 + \log x_{\textrm{\scriptsize{h}}_{2}\textrm{\scriptsize{o}}}'$ ] & [ reaction fe : feoh ] + fe(oh)@xmath0 & @xmath146 + 2\log x_{\textrm{\scriptsize{h}}_{2}\textrm{\scriptsize{o}}}'$ ] & [ reaction fe : fe(oh)2 ] + feo & @xmath147 + \log x_{\textrm{\scriptsize{h}}_{2}\textrm{\scriptsize{o}}}'$ ] & [ reaction fe : feo ] + fes & @xmath148 $ ] & [ reaction fe : fes ] + [ table iron reactions ]    llc mg & @xmath149 $ ] & + mgh & @xmath150 $ ] & [ reaction mg : mgh ] + mgoh & @xmath151 $ ] & [ reaction mg : mgoh ] + mg(oh)@xmath0 & @xmath152 $ ] & [ reaction mg : mg(oh)2 ] + mgo & @xmath153 $ ] & [ reaction mg : mgo ] + mgs & @xmath154 $ ] & [ reaction mg : mgs ] + mg & @xmath155- \log x_{\textrm{\scriptsize{h}}_{2}\textrm{\scriptsize{o}}}'$ ] & [ reaction forsterite condensation ] , [ reaction enstatite condensation ] + mgh & @xmath156- \log x_{\textrm{\scriptsize{h}}_{2}\textrm{\scriptsize{o}}}'$ ] & [ reaction mg : mgh ] + mgoh & @xmath157 & [ reaction mg : mgoh ] + mg(oh)@xmath0 & @xmath158 + \log x_{\textrm{\scriptsize{h}}_{2}\textrm{\scriptsize{o}}}'$ ] & [ reaction mg : mg(oh)2 ] + mgo & @xmath159 & [ reaction mg : mgo ] + mgs & @xmath160 & [ reaction mg : mgs ] + [ table magnesium reactions ]    llc sio & @xmath161 $ ] & + sis & @xmath162 $ ] & [ reaction sio : sis ] + si & @xmath163 & [ reaction sio : si ] + sih & @xmath164 & [ reaction sio : sih ] + sih@xmath0 & @xmath165 & [ reaction sio : sih2 ] + sih@xmath2 & @xmath166 & [ reaction sio : sih3 ] + sih@xmath1 & @xmath167 & [ reaction sih4:sio ] + sih@xmath2f & @xmath168 $ ] & [ reaction sio : sih3f ] + sih@xmath0f@xmath0 & @xmath169 $ ] & [ reaction sio : sih2f2 ] + sio & @xmath170 - \log x_{\textrm{\scriptsize{h}}_{2}\textrm{\scriptsize{o}}}'$ ] & [ reaction forsterite condensation ] , [ reaction enstatite condensation ] + sis & @xmath171 - 2\log x_{\textrm{\scriptsize{h}}_{2}\textrm{\scriptsize{o}}}'$ ] & [ reaction sio : sis ] + si & @xmath172 - 2\log x_{\textrm{\scriptsize{h}}_{2}\textrm{\scriptsize{o}}}'$ ] & [ reaction sio : si ] + sih & @xmath173 - 2\log x_{\textrm{\scriptsize{h}}_{2}\textrm{\scriptsize{o}}}'$ ] & [ reaction sio : sih ] + sih@xmath0 & @xmath174 - 2\log x_{\textrm{\scriptsize{h}}_{2}\textrm{\scriptsize{o}}}'$ ] & [ reaction sio : sih2 ] + sih@xmath2 & @xmath175- 2\log x_{\textrm{\scriptsize{h}}_{2}\textrm{\scriptsize{o}}}'$ ] & [ reaction sio : sih3 ] + sih@xmath1 & @xmath176- 2\log x_{\textrm{\scriptsize{h}}_{2}\textrm{\scriptsize{o}}}'$ ] & [ reaction sih4:sio ] + sih@xmath2f & @xmath177 - 2\log x_{\textrm{\scriptsize{h}}_{2}\textrm{\scriptsize{o}}}'$ ] & [ reaction sio : sih3f ] + sih@xmath0f@xmath0 & @xmath178 & [ reaction sio : sih2f2 ] + [ table silicon reactions ]


the antarctic neutrino telescope amanda  @xcite and the future km@xmath0 icecube  @xcite are designed to observe high - energy neutrinos from astrophysical sources .
the ice is instrumented with photomultipliers to pick up the cherenkov light from secondary charged particles . in order to reach the large volume needed to detect the expected small fluxes at high energies ,
the density of optical modules is far too sparse to measure , for example , solar neutrinos .
however , it has been recognized for a long time that these instruments can detect a supernova ( sn ) neutrino burst because the cherenkov glow of the ice can be identified as time - correlated noise among all phototubes  @xcite .
this approach has been used by amanda to exclude the occurrence of a galactic sn over a recent observation period  @xcite .    for amanda
the physics potential of a possible sn observation is essentially limited to its detection , notably in the context of the supernova early warning system ( snews ) that would alert the astronomical community several hours before the optical explosion  @xcite . for the future icecube with 4800 optical modules
, however , the number of detected cherenkov photons would be of order @xmath2 and thus so large that several interesting physics questions could be addressed in earnest .
the observed quantity is the number of cherenkov photons caused by the sn neutrinos as a function of time , i.e.  a measure of the energy deposited by the neutrinos in the ice .
therefore , the information about the sn signal is far more limited than what can be extracted from a high - statistics observation in super - kamiokande or other low - energy experiments that detect individual events .
however , galactic sne are so rare , perhaps a few per century , that the chances of observing one depend crucially on the long - term stability of the neutrino observatories .
once icecube has been built it may well operate for several decades , backing up the low - energy experiments . besides the detection and associated early warning one could measure important details of the neutrino light - curve , for example the existence and duration of the initial sn accretion phase , the overall duration of the cooling phase , and so forth . such an observation would provide a plethora of astrophysically valuable information
.    however , from the perspective of neutrino physics a _ simultaneous _ observation in both icecube and another large detector such as super - kamiokande or hyper - kamiokande would be especially useful .
assuming that the neutrinos have traversed significantly different paths through the earth , the two signals could well show measurable differences caused by neutrino oscillations in matter  .
as this earth effect shows up only for certain combinations of neutrino mixing parameters , a dual observation may well distinguish , for example , between the normal and inverted neutrino mass hierarchy .
it is well known that observing sn neutrinos with two or more detectors with different earth - crossing lengths is extremely useful , but icecube s potential has not been explored in this context . with all the low - energy observatories being in the northern hemisphere , icecube s location in antarctica is uniquely complementary for this purpose .
any oscillation signature depends on the small flavor - dependent differences between the fluxes and spectra at the source .
if these differences were as large as had been assumed until recently there would be little question about icecube s usefulness for co - detecting the earth effect . however , a more systematic study of the flavor - dependence of the sn neutrino fluxes and spectra reveals that these differences are more subtle , although by no means negligible  @xcite
. we evaluate icecube s potential as a co - detector from the perspective of these `` pessimistic '' assumptions about the primary fluxes and spectra .
this paper is organized as follows . in sec .
[ at - icecube ] , we show that the neutrino signal from a galactic supernova can be measured at icecube with a sub - percent statistical precision . in sec .
[ earth - effects ] , we calculate the earth matter effects on this signal and illustrate that it is possible to detect them in conjunction with another high statistics experiment .
sec .  [ concl ] concludes .
the sn neutrinos streaming through the antarctic ice interact according to @xmath3 and some other less important reactions .
the positrons , in turn , emit cherenkov light that is picked up by the optical modules ( oms ) frozen into the ice . while the expected number of detected photons per om was calculated in refs .
@xcite , we revisit their estimate for two reasons .
first , the sn signal was directly scaled to the historical sn  1987a observation in kamiokande  ii so that the exact assumptions about the neutrino flux are not directly apparent .
second , the expected number of cherenkov photons detected by one om was based on estimating an effective ice volume seen by one  om . however , it is much simpler to work in the opposite direction and start with the homogeneous and isotropic cherenkov glow of the ice caused by the sn neutrinos .
the om is immersed in this diffuse bath of photons and picks up a number corresponding to its angular acceptance and quantum efficiency .
as a first simplification we limit ourselves to the signal caused by the inverse @xmath4 reaction @xmath3 .
the @xmath5 fluence ( time - integrated flux ) at earth is @xmath6 we define the `` sn fudge factor '' as @xmath7 where @xmath8 is the total energy leaving the sn in the form of @xmath5 after flavor oscillations have been included , @xmath9 is the average @xmath5 energy , and @xmath10 the distance .
the energy deposited in the ice per target proton is @xmath11 . for the inverse @xmath4 cross section we ignore weak - magnetism and
recoil corrections and also the difference between @xmath5 and positron energy so that  @xcite @xmath12 for the neutrino flux of each neutrino and anti - neutrino species we assume a distribution of the form  @xcite @xmath13\,,\ ] ] where @xmath14 is the average energy , @xmath15 a parameter that typically takes on values 2.55 depending on the flavor and the phase of neutrino emission , and @xmath16 the overall flux at the detector in units of @xmath17 .
this distribution implies @xmath18 altogether we thus find @xmath19 with @xmath20 this fudge factor can also be taken to include deviations from the simplified energy dependence of the cross section and deviations from the assumed spectral shape .
the cherenkov angle for photon emission by a charged particle is @xmath21 where @xmath22 is the medium s refractive index and @xmath4 the particle s velocity . with @xmath23 for ice , neglecting the @xmath24-dependence , and @xmath25 we have @xmath26 .
a particle with unit charge produces cherenkov photons per unit path length and per unit wavelength band according to @xmath27 where @xmath28 is the fine - structure constant .
assuming that @xmath22 and thus @xmath29 are independent of wavelength we integrate over @xmath24 and find @xmath30 taking the useful wavelength range to be 300600  nm this translates into 319 photons per cm pathlength . taking the positron mean free path to be 12  cm for an energy of 20  mev , and taking it to be proportional to its energy , the number of useful cherenkov photons per deposited neutrino energy is @xmath31 with yet another fudge factor @xmath32 .
the density of ice is @xmath33 , corresponding to about @xmath34 proton targets .
therefore , the sn neutrinos produce @xmath35 useful cherenkov photons per unit volume of ice .
multiplying this number with the speed of light and dividing by @xmath36 gives us the resulting diffuse photon flux in units of @xmath37 .
however , the average lifetime of these photons is @xmath38 with @xmath39 the absorption length .
therefore , the neutrino - induced photon fluence is found by multiplying the flux with @xmath38 , @xmath40 where @xmath41 .
the number of events produced by this fluence in a given om depends on the average quantum efficiency taken to be @xmath42 .
in addition , it depends on the angular acceptance , i.e.  the effective photo cathode detection area @xmath43 times the angular acceptance range @xmath44 . therefore , in one om we expect @xmath45 with @xmath46 this result is independent of the presence of bubbles in the ice that scatter the photons .
the cherenkov glow of the ice represents an isotropic and homogeneous distribution that is not changed by elastic scattering .      in order to compare our result with the one derived in ref .
@xcite we need to translate their assumptions into our fudge factors .
the @xmath5 distribution was taken to follow a fermi - dirac spectrum with @xmath47  mev , implying @xmath48  mev .
the distance of the sn was taken to be 10  kpc , and the total energy release was scaled to the kamiokande  ii signal for sn  1987a . with our choice of the @xmath4
cross section these assumptions correspond to @xmath49  erg , i.e.  to @xmath50 .
these authors also used a quadratic energy dependence of the cross section . integrating over their fermi - dirac spectrum they effectively used @xmath51 .
further , they assumed 3000 useful cherenkov photons for 20  mev deposited energy , i.e.  effectively @xmath52 .
for the absorption length they used 300  m , i.e.  @xmath53 . finally , they assumed a quantum efficiency of 25% , a cathode area of 280  @xmath54 , and an acceptance range of @xmath55 , i.e.@xmath56 .
altogether , we find for these assumptions @xmath57 per om .
this compares with 273 in ref .
@xcite , i.e.  our result is larger by a factor  1.5 .
the result in ref .
@xcite was backed up by a detailed monte carlo treatment of the production and propagation of cherenkov photons in the amanda detector .
therefore , the difference may well relate to details of the om acceptance and wavelength - dependent quantum efficiency and photon propagation .
many of these details will be different in icecube where 10  inch photomultiplier tubes and different regions of ice will be used .
detailed values for the detector - dependent fudge factors must be determined specifically for icecube once it has been built .
the main difference between the assumptions in ref .
@xcite and our estimate is the absorption length . when using amanda as a sn observatory a realistic value was taken to be around 100  m  @xcite .
the vast difference between these estimates is that the former was based on the measured absorption length in a dust - free region of the ice . for our further estimates
we stick to 100  m as a conservative assumption .      in our derivation
we have used the time - integrated neutrino flux , amounting to the assumption of a stationary situation .
the absorption time for photons is very small , @xmath58 .
the sn signal will vary on time scales exceeding 10  ms .
therefore , the cherenkov glow of the ice follows the time - variation of the sn signal without discernible inertia .
hence one may replace the neutrino fluence with a time - dependent flux and @xmath59 with an event rate @xmath60 .
moreover , for our further discussion it will be useful to consolidate our fudge - factors into one describing the detector response , and others characterizing the neutrino flux .
therefore , we summarize our prediction for the event rate per om in the form @xmath61 where @xmath62 is the @xmath5 luminosity after flavor oscillations and @xmath63 here , @xmath64 also includes corrections for the energy dependence of the @xmath4 cross section .
we stress that our simple estimate of the counting rate primarily serves the purpose of determining its magnitude relative to the background .
the important feature is that the signal relative to the background can be determined with a good statistical precision . of course , for an absolute detector calibration a detailed modeling would be necessary . for our present purpose , however , even an uncertainty of several 10% in our estimated counting rate is irrelevant .
icecube will have 4800  oms so that one expects a total event number of @xmath65 , taking all fudge factors to be unity .
assuming a background counting rate of 300  hz per om over as much as 10  s this compares with a background rate of @xmath66 .
assuming poisson fluctuations , the uncertainty of this number is @xmath67 , i.e.  0.25% of the sn signal .
therefore , one can determine the sn signal with a statistical sub - percent precision , ignoring for now problems of absolute detector calibration .    in order to illustrate the statistical power of icecube to observe a sn signal we use two different numerical sn simulations .
the first was performed by the livermore group  @xcite that involves traditional input physics for mu- and tau - neutrino interactions and a flux - limited diffusion scheme for treating neutrino transport .
the great advantage of this simulation is that it covers the full evolution from infall over the explosion to the kelvin - helmholtz cooling phase of the newly formed neutron star .
we show the livermore @xmath5 and @xmath68 lightcurves in fig  [ fig : simulations ] ( left panels ) . here and in the following
we take @xmath68 to stand for either @xmath69 or @xmath70 .
apart from very small differences the sn fluxes and spectra are thought to be equal for @xmath71 , @xmath72 , @xmath69 , and @xmath70 .
our second simulation was performed with the garching code  @xcite .
it includes all relevant neutrino interaction rates , including nucleon bremsstrahlung , neutrino pair processes , weak magnetism , nucleon recoils , and nuclear correlation effects .
the neutrino transport part is based on a boltzmann solver .
the neutrino - radiation hydrodynamics program allows one to perform spherically symmetric as well as multi - dimensional simulations .
the progenitor model is a 15@xmath73 star with a @xmath74 iron core .
the period from shock formation to 468  ms after bounce was evolved in two dimensions .
the subsequent evolution of the model is simulated in spherical symmetry . at 150
ms the explosion sets in , although a small modification of the boltzmann transport was necessary to allow this to happen  @xcite .
unmanipulated full - scale models with an accurate treatment of the microphysics currently do not obtain explosions  @xcite .
this run will be continued beyond the current epoch of 750  ms post bounce ; we here use the preliminary results currently available  @xcite .
we show the garching @xmath5 and @xmath68 lightcurves in fig  [ fig : simulations ] ( right panels ) .
we take the livermore simulation to represent traditional predictions for flavor - dependent sn neutrino fluxes and spectra that were used in many previous discussions of sn neutrino oscillations .
the garching simulation is taken to represent a situation when the @xmath68 interactions are more systematically included so that the flavor - dependent spectra and fluxes are more similar than had been assumed previously  @xcite .
we think it is useful to juxtapose the icecube response for both cases .
another difference is that in livermore the accretion phase lasts longer .
since the explosion mechanism is not finally settled , it is not obvious which case is more realistic .
moreover , there could be differences between different sne .
the overall features are certainly comparable between the two simulations .    in fig .
[ fig : simulationsicecube ] we show the expected counting rates in icecube on the basis of eq .
( [ eq : rateprediction ] ) for an assumed distance of 10  kpc and 4800 oms for the livermore ( left ) and garching ( right ) simulations .
we also show this signal in 50  ms bins where we have added noise from a background of 300  hz per om .
the baseline is at the average background rate so that negative counts correspond to downward background fluctuations .
one could easily identify the existence and duration of the accretion phase and thus test the standard delayed - explosion scenario .
one could also measure the overall duration of the cooling phase and thus exclude the presence of significant exotic energy losses .
therefore , many of the particle - physics limits based on the sn  1987a neutrinos  @xcite could be supported with a statistically serious signal .
if the sn core were to collapse to a black hole after some time , the sudden turn - off of the neutrino flux could be identified . in short ,
when a galactic sn occurs , icecube is a powerful stand - alone neutrino detector , providing us with a plethora of information that is of fundamental astrophysical and particle - physics interest .
in addition , icecube is extremely useful as a co - detector with another high - statistics observatory to measure neutrino oscillation effects , a topic that we now explore .
neutrino oscillations are now firmly established by measurements of solar and atmospheric neutrinos and the kamland and k2k long - baseline experiments  .
evidently the weak interaction eigenstates @xmath75 , @xmath71 and  @xmath72 are non - trivial superpositions of three mass eigenstates @xmath76 , @xmath77 and  @xmath78 , @xmath79 where @xmath80 is the leptonic mixing matrix that can be written in the canonical form @xmath81 here @xmath82 and @xmath83 etc . , and @xmath84 is a phase that can lead to cp - violating effects , that are , however , irrelevant for sn neutrinos .
the mass squared differences relevant for the atmospheric and solar neutrino oscillations obey a hierarchy @xmath85 .
this hierarchy , combined with the observed smallness of the angle @xmath86 at chooz @xcite implies that the atmospheric neutrino oscillations essentially decouple from the solar ones and each of these is dominated by only one of the mixing angles .
the atmospheric neutrino oscillations are controlled by @xmath87 that may well be maximal ( 45@xmath88 ) .
the solar case is dominated by @xmath89 , that is large but not maximal . from a global 3-flavor analysis of all data
one finds the 3@xmath90 ranges for the mass differences @xmath91 and mixing angles summarized in table  [ tab : mixingparameters ] .
@lll observation&mixing angle&@xmath92 [ mev@xmath93 + sun , kamland & @xmath94 27@xmath8842@xmath88 & @xmath95 55190 + atmosphere , k2k & @xmath96 32@xmath8860@xmath88 & @xmath97 14006000 + chooz & @xmath98 & @xmath99 +    a sn core is essentially a neutrino blackbody source , but small flavor - dependent differences of the fluxes and spectra remain .
we denote the fluxes of @xmath100 and @xmath101 at earth that would be observable in the absence of oscillations by @xmath102 and @xmath103 , respectively . in the presence of oscillations
a @xmath5 detector actually observes @xmath104 f_x^0\ , , \label{fedbar}\ ] ] where @xmath105 is the @xmath5 survival probability after propagation through the sn mantle and perhaps part of the earth before reaching the detector .
a significant modification of the survival probability due to the propagation through the earth appears only for those combinations of neutrino mixing parameters shown in table  [ tab : earthcases ] .
the earth matter effect depends strongly on two parameters , the sign of @xmath106 and the value of @xmath107  @xcite .
the `` normal hierarchy '' corresponds to @xmath108 , i.e.  @xmath109 , whereas the `` inverted hierarchy '' corresponds to @xmath110 , i.e.@xmath111 .
note that the presence or absence of the earth effect discriminates between values of @xmath112 less or greater than @xmath113 , i.e.  @xmath86 less or larger than about @xmath114 .
thus , the earth effect is sensitive to values of @xmath86 that are much smaller than the current limit .
@lll 13-mixing&normal hierarchy&inverted hierarchy + @xmath115&@xmath75 and @xmath5 & @xmath75 and @xmath5 + @xmath116&@xmath5&@xmath75 +    let us consider those scenarios where the mass hierarchy and the value of @xmath86 are such that the earth effect appears for @xmath5 . in such cases the @xmath5 survival probability @xmath105 is given by @xmath117 where the energy dependence of all quantities will always be implicit .
here @xmath118 is the mixing angle between @xmath100 and @xmath119 in earth matter while @xmath120 is the mass squared difference between the two anti - neutrino mass eigenstates @xmath121 and @xmath119 in units of @xmath122ev@xmath123 , @xmath124  is the distance traveled through the earth in units of 1000  km , and @xmath125 is the neutrino energy in mev .
we have assumed a constant matter density inside the earth , which is a good approximation for @xmath126 , i.e.  as long as the neutrinos do not pass through the core of the earth .      in order to calculate the extent of the earth effect for icecube
, we will assume that the relevant mixing parameters are @xmath127 and @xmath128 .
we further assume that the source spectra are given by the functional form eq .
( [ eq : spectralform ] ) .
the values of the parameters @xmath15 and @xmath129 for both the @xmath100 and @xmath68 spectra are in general time dependent .    in fig .
[ fig : exampleeartheffect ] we show the variation of the expected icecube signal with earth - crossing length @xmath124 for the two sets of parameters detailed in table  [ tab : examplecases ] .
the first could be representative of the accretion phase , the second of the cooling signal .
we use the two - density approximation for the earth density profile , where the core has a density of @xmath130 and a radius of 3500  km , while the density of the earth mantle was taken to be @xmath131 .
we observe that for short distances , corresponding to near - horizontal neutrino trajectories , the signal varies strongly with @xmath124 . between about 3,000 and 10,500  km
it reaches an asymptotic value that we call the `` asymptotic mantle value . '' for case  ( a ) , this value corresponds to about 1.5% depletion of the signal , whereas for  ( b ) it corresponds to about 6.5% depletion .
@llllllll example&phase&@xmath9 & @xmath132 & @xmath133&@xmath134 & @xmath135 & asymptotic + & & [ mev]&[mev]&&&&earth effect + ( a)&accretion&15&17&4&3&1.5&@xmath136% + ( b)&cooling & 15&18&3&3&0.8&@xmath137% +    beyond an earth - crossing length of @xmath13810,500  km , the neutrinos have to cross the earth core with another large jump in density .
the core effects change the asymptotic mantle value by @xmath139% as can be seen in fig .
[ fig : exampleeartheffect ] .
we neglect the core effects in the following analysis , and the `` asymptotic value '' always refers to the asymptotic mantle value .
for the largest part of the sky the earth effect either appears with this asymptotic value ( `` neutrinos coming from below '' ) , or it does not appear at all ( `` neutrinos from above '' ) .
therefore , we now focus on the asymptotic value and study how the signal modification depends on the assumed flux parameters . in table
[ tab : asymptotic1 ] we show the signal modification for @xmath140 , @xmath141 , and @xmath142 as a function of @xmath132 and the flux ratio @xmath135 . in table
[ tab : asymptotic2 ] we show the same with @xmath143 .
the results are shown in the form of contour plots in fig .
[ fig : cont ] .
@lllllll flux ratio & + @xmath135 & 15&16&17&18&19&20 + 2.0 & 1.026 & 1.014 & 1.002 & 0.988 & 0.974 & 0.960 + 1.9 & 1.023 & 1.011 & 0.999 & 0.985 & 0.971 & 0.956 + 1.8 & 1.021 & 1.009 & 0.995 & 0.982 & 0.967 & 0.952 + 1.7 & 1.018 & 1.005 & 0.992 & 0.978 & 0.963 & 0.948 + 1.6 & 1.015 & 1.002 & 0.988 & 0.974 & 0.959 & 0.944 + 1.5 & 1.012 & 0.998 & 0.984 & 0.969 & 0.954 & 0.939 + 1.4 & 1.008 & 0.994 & 0.980 & 0.965 & 0.949 & 0.934 + 1.3 & 1.004 & 0.990 & 0.975 & 0.960 & 0.944 & 0.928 + 1.2 & 1.000 & 0.985 & 0.970 & 0.954 & 0.938 & 0.922 + 1.1 & 0.995 & 0.980 & 0.964 & 0.948 & 0.932 & 0.915 + 1.0 & 0.989 & 0.974 & 0.957 & 0.941 & 0.925 & 0.908 + 0.9 & 0.983 & 0.967 & 0.950 & 0.934 & 0.917 & 0.901 + 0.8 & 0.976 & 0.959 & 0.942 & 0.925 & 0.909 & 0.892 + 0.7 & 0.967 & 0.950 & 0.933 & 0.916 & 0.899 & 0.883 + 0.6 & 0.958 & 0.940 & 0.923 & 0.906 & 0.889 & 0.873 + 0.5 & 0.946 & 0.928 & 0.911 & 0.894 & 0.877 & 0.862 +    @lllllll flux ratio & + @xmath135 & 15&16&17&18&19&20 + 2.0 & 1.036 & 1.024 & 1.012 & 1.000 & 0.986 & 0.972 + 1.9 & 1.033 & 1.022 & 1.010 & 0.996 & 0.983 & 0.968 + 1.8 & 1.031 & 1.019 & 1.006 & 0.993 & 0.979 & 0.964 + 1.7 & 1.028 & 1.016 & 1.003 & 0.989 & 0.975 & 0.960 + 1.6 & 1.025 & 1.013 & 0.999 & 0.985 & 0.971 & 0.955 + 1.5 & 1.022 & 1.009 & 0.995 & 0.981 & 0.966 & 0.951 + 1.4 & 1.019 & 1.005 & 0.991 & 0.976 & 0.961 & 0.945 + 1.3 & 1.015 & 1.001 & 0.986 & 0.971 & 0.955 & 0.940 + 1.2 & 1.010 & 0.996 & 0.981 & 0.965 & 0.949 & 0.933 + 1.1 & 1.006 & 0.991 & 0.975 & 0.959 & 0.943 & 0.927 + 1.0 & 1.000 & 0.985 & 0.969 & 0.952 & 0.936 & 0.919 + 0.9 & 0.994 & 0.978 & 0.961 & 0.945 & 0.928 & 0.911 + 0.8 & 0.986 & 0.970 & 0.953 & 0.936 & 0.919 & 0.903 + 0.7 & 0.978 & 0.961 & 0.944 & 0.926 & 0.910 & 0.893 + 0.6 & 0.968 & 0.950 & 0.933 & 0.916 & 0.899 & 0.882 + 0.5 & 0.956 & 0.938 & 0.920 & 0.903 & 0.886 & 0.870 +    even for mildly different fluxes or spectra the signal modification is several percent , by far exceeding the statistical uncertainty of the icecube signal , although the _ absolute _ calibration of icecube may remain uncertain to within several percent .
however , the signal modification will vary with time during the sn burst . during
the early accretion phase that is expected to last for a few 100  ms and corresponds to a significant fraction of the overall signal , the @xmath68 flux may be almost a factor of 2 smaller than the @xmath5 flux , but it will be slightly hotter and less pinched  @xcite .
this corresponds to case  ( a ) above ; it is evident from fig .
[ fig : exampleeartheffect ] and table  [ tab : asymptotic1 ] that this implies that the earth effect is very small . during the kelvin - helmholtz cooling phase the flux ratio is reversed with more @xmath68 being emitted than @xmath5 , but still with the same hierarchy of energies .
this corresponds to case  ( b ) ; in this case the earth effect could be about 6% .
this time dependence may allow one to detect the earth effect without a precise absolute detector calibration .    in order to illustrate the time dependence of the earth effect we show in fig .
[ fig : percenteffect ] the expected counting rate in icecube for both the livermore ( left panels ) and garching ( right panels ) simulations . in the upper panels we show the expected counting rate with flavor oscillations in the sn mantle , but
no earth effect ( solid lines ) , or with the asymptotic earth effect ( dashed lines ) that obtains for a large earth - crossing path .
naturally the differences are very small so that we show in the lower panels the ratio of these curves , i.e.  the expected counting rate with / without earth effect as a function of time for both livermore and garching .
while for the livermore simulation there is a large earth effect even at early times , the change from early to late times in both cases is around 45% .
therefore , the most model - independent signature is a time variation of the earth effect during the sn neutrino signal .    in order to demonstrate the statistical significance of these effects
we integrate the expected signal for both simulations separately for the accretion phase and the subsequent cooling phase ; the results are shown in table  [ tab : integratedrates ] .
for both simulations the earth effect itself and its change with time is statistically highly significant .
based on the livermore simulation , the earth effect is much more pronounced than in garching , the latter involving more up - to - date input physics for neutrino transport .
however , the _ difference _ between the earth effect during accretion and cooling is not vastly different between the two simulations . recalling that the absolute detector calibration may be very uncertain so that one has to rely on the temporal variation of the earth effect , the difference between livermore and garching becomes much smaller .
we expect that it is quite generic that the temporal change of the earth effect is a few percent of the overall counting rate .
@lrrrr & & + & & & & + integration time [ s]&00.500&0.5003&00.250&0.2500.700 + sn signal [ counts ] + no earth effect & 519,080&818,043&173,085&407,715 + asymptotic earth effect&488,093&751,137&171,310&390,252 + difference & 30,987 & 66,906 & 1,775 & 17,463 + fractional difference&@xmath144%&@xmath145%&@xmath146%&@xmath147% + background [ counts]&720,000&4,320,000&360,000&648,000 + @xmath148/signal&0.16%&0.25%&0.35%&0.20% +      one can measure the earth effect in icecube only in conjunction with another high - statistics detector . we do not attempt to simulate in detail the sn signal in this other detector but simply assume that it can be measured with a precision at least as good as in icecube .
one candidate is super - kamiokande , a water cherenkov detector that would measure around @xmath149 events from a galactic sn at a distance of 10  kpc .
therefore , the statistical precision for the total neutrino energy deposition in the water is around 1% and thus worse than in icecube .
even though super - kamiokande will measure a larger number of cherenkov photons than icecube , a single neutrino event will cause an entire cherenkov ring to be measured , i.e.  the photons are highly correlated .
therefore , in the estimated statistical @xmath150 fluctuation of the signal , the fluctuating number @xmath151 is that of the detected neutrinos . if the future hyper - kamiokande is built , its fiducial volume would be about 30 times that of super - kamiokande . in this case
the statistical signal precision exceeds that of icecube for the equivalent observable .
we denote the equivalent icecube signal measured by super- or hyper - kamiokande as @xmath152 and the icecube signal as @xmath153 .
if the distances traveled by the neutrinos before reaching these two detectors are different , the earth effect on the neutrino spectra may be different , which will reflect in the ratio @xmath154 .
of course , in the absence of the earth effect this ratio equals unity by definition .
the geographical position of icecube with respect to super- or hyper - kamiokande at a latitude of @xmath155 is well - suited for the detection of the earth effect through a combination of the signals . using fig .
[ fig : exampleeartheffect ] we can already draw some qualitative conclusions about the ratio @xmath156 .
clearly , @xmath157 if neutrinos do not travel through the earth before reaching either detector . if the distance traveled by neutrinos through the earth is more that 3000  km for both detectors , the earth effects on both @xmath152 and @xmath158 are nearly equal and their ratio stays around unity .
if the neutrinos come `` from above '' for sk and `` from below '' for icecube , or vice versa , the earth matter effect will shift this ratio from unity .    in fig .
[ fig : earthmap ] , we show contours of @xmath159 for the sn position in terms of the location on earth where the sn is at the zenith .
the map is an area preserving hammer - aitoff projection so that the sizes of different regions in the figure gives a realistic idea of the `` good '' and `` bad '' regions of the sky . in order to generate the contours we use the parameters of case  ( b ) in table  [ tab : examplecases ] so that the asymptotic suppression of the signal is about 6.5% .
the sky falls into four distinct regions depending on the direction of the neutrinos relative to either detector as described in table  [ tab : skymapregions ] .
when the neutrinos come from above for both detectors ( region  d ) there is no earth effects .
if they come from below in both ( region  c ) , the earth effect is large in both .
depending on the exact distance traveled through the earth , the event ratio can be large , but generally fluctuates around  1 . in the other regions where the neutrinos come from above for one detector and from below for the other ( regions a and  b ) the relative effect is large .
@lllll region&sky fraction & & @xmath156 + & & icecube&super - k & + a&0.35&below&above&1.070 + b&0.35&above&below&0.935 + c&0.15&below&below&fluctuations around 1 + d&0.15&above&above&1 +
for assumptions about the flavor - dependent sn neutrino fluxes and spectra that agree with state - of - the - art studies , the earth matter effect on neutrino oscillations shows up in the icecube signal of a future galactic sn on the level of a few percent .
if the icecube signal can be compared with another high - statistics signal , notably in super - kamiokande or hyper - kamiokande , the earth effect becomes clearly visible as a difference between the detectors .
as one is looking for a signal modification in the range of a few percent , the absolute detector calibration may not be good enough in one or both of the instruments .
however , for typical numerical sn simulations the effect is time dependent and most notably differs between the early accretion phase and the subsequent neutron star cooling phase .
therefore , one would have to search for a temporal variation of the relative detector signals of a few percent .
the large number of optical modules in icecube renders this task statistically possible .
in fact depending on the differences in flavor - dependent fluxes , the statistical accuracy of super - kamiokande may turn out to be the limiting factor .
this limitation is not significant for hyper - kamiokande .
the unique location of icecube in antarctica implies that for a large portion of the sky this detector sees the sn through the earth when super- and hyper - kamiokande sees it from above , or the other way round , i.e.  the chances of a relative signal difference between the detectors are large . if both detectors were to see the sn from above there would be no earth effect to detect .    assuming that the magnitude of the mixing angle @xmath86 can be established to be large in the sense of @xmath160 by a long - baseline experiment  @xcite
, observing the earth effect in sn anti - neutrinos implies the normal mass hierarchy . on the other hand ,
if @xmath161 has been established , the earth effect is unavoidable . not observing
it would imply that the primary sn neutrino fluxes and spectra are more similar than indicated by state - of - the - art numerical simulations .    if @xmath160 is known , and we do not observe the earth effect , it still does not prove the inverted mass hierarchy .
it could also mean that we do not properly understand the flavor - dependent source fluxes and spectra .
therefore , even if @xmath160 is known , our method only allows one to detect the normal mass hierarchy , it does not strictly allow one to exclude it . as far as neutrino parameters are concerned , only a positive detection of the earth effect would count for much .
of course , a normal mass hierarchy and @xmath160 is certainly a plausible scenario so that expecting a positive identification of the earth effect is not a far - fetched possibility .    in summary , even though galactic sne are rare , the anticipated longevity of icecube and the long - term neutrino program at super- or hyper - kamiokande imply that detecting the earth effect in a sn neutrino burst is certainly a distinct possibility .
this could identify the normal neutrino mass hierarchy , a daunting task at long - baseline experiments  @xcite . given the difficulty of pinning down the mass hierarchy at long - baseline experiments
, both icecube and super- or hyper - kamiokande should take all instrumental and experimental steps required to ensure the feasibility of a high - statistics simultaneous sn observation .
this work was supported , in part , by the deutsche forschungsgemeinschaft under grant no .
sfb-375 and by the european science foundation ( esf ) under the network grant no .  86 neutrino astrophysics .
we thank francis halzen and robert buras for helpful comments on an early version of this manuscript .
99        f.  halzen , j.  e.  jacobsen and e.  zas , `` possibility that high - energy neutrino telescopes could detect supernovae , '' _ phys .
d _ * 49 * , 1758 ( 1994 ) .
f.  halzen , j.  e.  jacobsen and e.  zas , `` ultra - transparent antarctic ice as a supernova detector , '' _ phys .
d _ * 53 * , 7359 ( 1996 ) [ astro - ph/9512080 ] .
j.  ahrens _ et al .
_ [ amanda collaboration ] , `` search for supernova neutrino - bursts with the amanda detector , '' _ astropart .
* 16 * , 345 ( 2002 ) [ astro - ph/0105460 ] .
k.  scholberg , `` snews : the supernova early warning system , '' astro - ph/9911359 .
snews : the supernova early warning system , http://hep.bu.edu/@xmath138snnet/    a.  s.  dighe and a.  y.  smirnov , `` identifying the neutrino mass spectrum from the neutrino burst from a supernova , '' _ phys .
d _ * 62 * , 033007 ( 2000 ) [ hep - ph/9907423 ] .
a.  s.  dighe , `` earth matter effects on the supernova neutrino spectra , '' hep - ph/0106325 . c.  lunardini and a.  y.  smirnov , `` supernova neutrinos : earth matter effects and neutrino mass spectrum , '' nucl .
b * 616 * , 307 ( 2001 ) [ hep - ph/0106149 ] . c.  lunardini and a.  y.  smirnov , `` probing the neutrino mass hierarchy and the 13-mixing with supernovae , '' hep - ph/0302033 .
k.  takahashi and k.  sato , `` earth effects on supernova neutrinos and their implications for neutrino parameters , '' _ phys .
d _ * 66 * , 033006 ( 2002 ) [ hep - ph/0110105 ] .
g.  g.  raffelt , `` mu- and tau - neutrino spectra formation in supernovae , '' _ astrophys . j. _ * 561 * , 890 ( 2001 ) [ astro - ph/0105250 ] .
r.  buras , h.  t.  janka , m.  t.  keil , g.  g.  raffelt and m.  rampp , `` electron - neutrino pair annihilation : a new source for muon and tau neutrinos in supernovae , '' _ astrophys . j. _ * 587 * , 320 ( 2003 ) [ astro - ph/0205006 ]
. m.  t.  keil , g.  g.  raffelt and h.  t.  janka , `` monte carlo study of supernova neutrino spectra formation , '' _ astrophys .
, in press ( 2003 ) [ astro - ph/0208035 ] . g.  g.  raffelt , m.  th .
keil , r.  buras , h .- t .
janka and m.  rampp , `` supernova neutrinos : flavor - dependent fluxes and spectra , '' proc
.  noon  03 ( 1014 february 2003 , kanazawa , japan ) , astro - ph/0303226 .
p.  vogel and j.  f.  beacom , `` angular distribution of neutron inverse beta decay , @xmath162 , '' _ phys .
* 60 * , 053003 ( 1999 ) [ hep - ph/9903554 ] .
t.  totani , k.  sato , h.  e.  dalhed and j.  r.  wilson , `` future detection of supernova neutrino burst and explosion mechanism , '' _ astrophys .
j. _ * 496 * , 216 ( 1998 ) [ astro - ph/9710203 ] .
m.  rampp and h .-
janka , `` radiation hydrodynamics with neutrinos : variable eddington factor method for core - collapse supernova simulations '' , _ astron .
_ * 396 * , 361 ( 2002 ) [ astro - ph/0203101 ] .
janka , r.  buras , k.  kifonidis , t.  plewa and m.  rampp , `` explosion mechanisms of massive stars '' , in : _ core collapse of massive stars _ , edited by c.l .
fryer ( kluwer academic publ . ,
dordrecht , 2003 ) [ astro - ph/0212314 ] .
r.  buras , m.  rampp , h .- t .  janka and k.  kifonidis , `` improved models of stellar core collapse and still no explosions : what is missing ? '' , _ physical review letters _ , submitted ( 2003 ) [ astro - ph/0303171 ]
. g.  g.  raffelt , `` particle physics from stars , '' _ ann .
nucl . part .
sci . _ * 49 * , 163 ( 1999 ) [ hep - ph/9903472 ] .
m.  c.  gonzalez - garcia and y.  nir , `` developments in neutrino physics , '' hep - ph/0202058 .
g.  l.  fogli , e.  lisi , a.  marrone , d.  montanino , a.  palazzo and a.  m.  rotunno , `` solar neutrino oscillation parameters after first kamland results , '' _ phys .
d _ * 67 * , 073002 ( 2003 ) [ hep - ph/0212127 ] .
j.  n.  bahcall , m.  c.  gonzalez - garcia and c.  pea - garay , `` solar neutrinos before and after kamland , '' _ jhep _ * 02 * , 009 ( 2003 ) [ hep - ph/0212147 ] .
p.  c.  de holanda and a.  y.  smirnov , `` lma msw solution of the solar neutrino problem and first kamland results , '' _ jcap _ * 02 * , 001 ( 2003 ) [ hep - ph/0212270 ] . m.  maltoni , t.  schwetz and j.  w.  valle , `` combining first kamland results with solar neutrino data , '' _ phys .
d _ * 67 * , 093003 ( 2003 ) [ hep - ph/0212129 ] .
m.  maltoni , t.  schwetz , m.  a.  tortola and j.  w.  valle , `` constraining neutrino oscillation parameters with current solar and atmospheric data , '' _ phys .
d _ * 67 * , 013011 ( 2003 ) [ hep - ph/0207227 ] .
m.  apollonio _ et al . _
[ chooz collaboration ] , `` limits on neutrino oscillations from the chooz experiment , '' _ phys .
b _ * 466 * , 415 ( 1999 ) [ hep - ex/9907037 ] . v.  d.  barger , s.  geer , r.  raja and k.  whisnant , `` exploring neutrino oscillations with superbeams , '' _ phys .
d _ * 63 * , 113011 ( 2001 ) [ hep - ph/0012017 ] .
a.  cervera , a.  donini , m.  b.  gavela , j.  j.  gomez cadenas , p.  hernandez , o.  mena and s.  rigolin , `` golden measurements at a neutrino factory , '' _ nucl .
b _ * 579 * , 17 ( 2000 ) ; erratum _ ibid . _ * 593 * , 731 ( 2001 ) [ hep - ph/0002108 ] . m.  freund , p.  huber and m.  lindner , `` systematic exploration of the neutrino factory parameter space including errors and correlations , '' _ nucl . phys .
b _ , * 615 * , 331 ( 2001 ) [ hep - ph/0105071 ] .