the goal of cosmology is to describe the evolution history of the universe as a whole , tracing it as far back as possible from the present epoch .
this seemingly daunting task is made tractable by the fact that the observed matter density in the universe appears to be approximately homogeneous when viewed on the largest scales .
further , observations of the cosmic microwave background ( cmb ) radiation reveal it to be isotropic ( after subtracting a dipole contribution due to our local motion ) down to the level of one part in @xmath0 .
the twin conditions of homogeneity and isotropy , if valid on all length scales , would immensely simplify the form of the metric and matter tensors used in the general relativistic description of cosmology .
it can be shown @xcite that on purely geometric grounds , the most general metric describing a homogeneous and isotropic universe takes on the form of the friedmann - lema^ itre - robertson - walker ( flrw ) line element .
the standard method of proceeding is to assume that the geometry of the universe can be described by a background metric with small perturbations , where the background is exactly homogeneous and isotropic and is hence of the flrw form .
( for a detailed discussion of the assumptions involved in the standard cosmology , see ref .
@xcite . )    this approach is adopted not only for the early universe , for example in describing big bang nucleosynthesis and the dynamics of the cmb anisotropies , but is also assumed to be valid at recent epochs . in particular ,
the luminosity distances constructed to compare data from type ia supernovae ( sne ) with theory , assume a metric of the flrw form with its dynamics given by the einstein equations for a homogeneous and isotropic universe .
we know , however , that the matter content of the universe today is _ not _ homogeneous on small scales - one finds inhomogeneity at the scale of individual galaxies for example ; and there is evidence @xcite that voids of order @xmath1 mpc is the conventional parameter which appears in the definition of the hubble constant @xmath2 kms@xmath3mpc@xmath3 . ] in diameter account for @xmath4-@xmath5 of the observed universe .
there is also evidence for voids 3 - 5 times this size @xcite , as well as local voids on smaller scales @xcite .
the evidence for homogeneity of the matter distribution in the universe only appears at scales larger than about @xmath6 mpc or so @xcite .
it is then reasonable to assume that the metric describing our universe must also reflect the inhomogeneity of the matter on the corresponding scales .
one may now ask whether it is at all justified to use the luminosity distances constructed using the flrw metric with its associated dynamics to compare theory with observational data from sne . in this paper we argue that the use of such luminosity distances is indeed justified , provided we are only interested in using such constructions in probing length scales larger than a certain minimum scale which is fixed by requiring that the assumption of homogeneity be valid ( a lower bound on such a scale could be @xmath6 mpc for example ) .
more importantly , we propose that the _ dynamics _ of this ( effective or template ) macroscopic metric must be governed not by the standard einstein equations applied to the flrw metric , but by the effective equations obtained after spatial averaging on constant time slices , as developed by buchert @xcite .
an outstanding open issue has been the following : how does one compare the predictions of the inhomogeneous , spatially averaged cosmology with observations , and should one be attempting an interpretation in terms of spatially averaged variables in the first place ?
for instance , one could study light propagation and determine the luminosity distance versus redshift ( @xmath7 ) relation which best fits the observations , within an _ inhomogeneous _ model of the universe .
following the important and interesting early work of dyer and roeder , where the angular - size and luminosity distances were obtained for beams propagating in `` lumpy '' flrw universes @xcite , detailed investigations have been carried out in @xcite .
such a program has also received considerable attention in the context of the spherically symmetric inhomogeneous lematre - tolman - bondi ( ltb ) models @xcite . however , it is not clear to us that this is necessarily the best way to proceed , at least as far as distant sne are concerned .
a realistic computation of light propagation in an inhomogeneous geometry seems difficult at best , and while the use of simplified toy models such as the ltb solutions in this endeavor may be justified on the grounds that we may inhabit a locally underdense region @xcite , it is not clear that the assumptions justifying the use of the toy model remain valid at large scales ( see also ref .
@xcite for a recent review of the ltb approach ) . and
after all , the concordance model infers the presence of a dark energy by comparing the predictions of the standard homogeneous and isotropic flrw cosmology against observations @xcite . if one decides to interpret the flrw cosmology as a template against which observations
may be compared , then in order to decide whether cosmic inhomogeneities could be an alternative to the concordance model , one must construct a homogeneous and isotropic cosmology , after a suitable smoothing procedure , and compare its predictions with observations such as the @xmath7  relation for type ia sne .
one should note that the effects of averaging within the standard perturbed flrw framework of cosmology , have been studied by several authors ( see , e.g. , refs .
@xcite ) , and it has been argued that these effects must necessarily be small .
the papers in refs .
@xcite however , do not necessarily address the entire problem .
for example , ref .
@xcite deals with only super - horizon fluctuations , whereas work done by others ( see , e.g , ref .
@xcite ) indicates that there may be nontrivial contributions due to averaging from _
_ sub__horizon fluctuations where the matter perturbations have gone nonlinear . and
while the authors of ref .
@xcite construct a cosmological model with negligible effects of averaging , there do exist toy models of inhomogeneous spacetimes in which these effects are large ( see , e.g. , ref .
@xcite ) .    furthermore ,
the possibility of _ nonperturbative _ effects of averaging has not , to our knowledge , been satisfactorily addressed , much less ruled out .
it is important to realize that current @xmath8-body simulations for example , which are based on a newtonian approximation of gravity , consistently miss possible relativistic effects of averaging .
this is because the modifications due to averaging in a newtonian setting can be shown to reduce to boundary terms in general @xcite , and to vanish in the case of spherical symmetry @xcite .
these terms therefore do not contribute in simulations which either set periodic boundary conditions on the simulation box , or use a spherically symmetric simulation region . as we will see , the `` backreaction '' effects in buchert s framework ( i.e. , the modifications to the cosmological equations due to averaging ) are coupled to the average of the spatial curvature , this being a completely relativistic effect .
nontrivial consequences of this coupling deserve more attention than they have received in the context of @xmath8-body simulations ( see ref .
@xcite ) .    we also note that an argument for describing the universe using flrw plus newtonian perturbations ( e.g. , as given in ref .
@xcite ) , is based on the fact that it is newtonian perturbations _ and _ a negative pressure component , which describe the universe very well .
in other words , one argues that the effect of averaging inhomogeneities must be small , but at the cost of adding an ill - understood dark energy component in the universe . on the other hand
one should allow for the following possibility ( as of now neither confirmed nor ruled out ) : at an initial epoch , say at the end of inflation , the perturbations in the homogeneous metric would give a small `` backreaction '' as a result of averaging . as a result of its coupling with the curvature ,
this backreaction would grow and become significant around the time of structure formation , when matter inhomogeneities on small scales become nonlinear .
the backreaction would also nontrivially affect the growth of the metric perturbations .
such a possibility could only be tested within a fully nonperturbative treatment of the evolution of inhomogeneities .
the goal of this paper is modest .
we argue that the scale factor , defined in the spatially averaged cosmology via the volume of a domain , should correspond to an effective flrw metric , provided the size of the averaging domain coincides with the scale at which cosmological homogeneity arises .
this allows us , in principle , to compare predictions of a spatially averaged cosmology with observations , in the standard manner , for instance by computing the @xmath7  relation .
( however , see ref .
@xcite for an alternative approach to matching the spatially averaged dynamics with observations . )
the predictions of the spatially averaged cosmology would in general differ from standard flrw cosmology , because the scale - factor now obeys modified einstein equations .
we will demonstrate that consistent models of spatially averaged flrw templates exist which can , at least in principle , explain observations such as those of the @xmath7  relation for supernovae .
[ intro ]    in a nutshell , our new results in this paper are that we propose a template metric which takes into account the kinematical effects of inhomogeneities in the universe .
to implement this idea , we first suggest that in buchert s spatial averaging scheme , the size of the averaging domain should coincide with the cosmological scale at which homogeneity sets in .
we then argue that , as a first step towards predicting the relation between inhomogeneities and observations , the expansion factor @xmath9 defined by buchert should be identified with the scale - factor @xmath10 in standard flrw cosmology .
the template metric is then the flrw metric , with @xmath9 playing the role of @xmath10 .
we then show , with the help of a few toy examples , that the luminosity distance - redshift relation in a model universe with an accelerating @xmath9 can be computed in the standard manner . to the best of our knowledge
, such a template metric has not been proposed in earlier work by other researchers .
at the length scales of cosmological interest , individual galaxies are small enough to be treated as point - like objects . given the large number ( @xmath11 ) of galaxies estimated to exist in the observed universe , it is also assumed that the matter content of the universe can be described to a good approximation as a continuum fluid . at the level of the matter distribution ,
this picture is comparable to that of a gas composed of a large number of molecules , wherein the internal structure of the molecules can be ignored when dealing with the macroscopic properties of the gas which is treated as a continuum fluid .
further , the complex intermolecular interactions in the gas at the microscopic level can be dealt with in a mean field approximation by considering volume averages of the quantities of interest .
this picture is made precise by introducing the concept of _ physically infinitesimal _
volumes @xcite or _ coarse - graining cells_. the fluid approximation is assumed to be valid provided the ( effectively ) infinitesimal volumes used in defining continuum quantities such as the matter density of the fluid , correspond to _ physical _ volumes which are large enough to contain a very large number of the particles constituting the fluid . at the same time , since these volumes are to be effectively infinitesimal , their physical size must also be much smaller than the scale of the fluid as a whole .
we can then say that the length scales @xmath12 being probed by these physically infinitesimal volumes must lie in a certain range @xmath13 which is determined by the details of the fluid ( see also @xcite ) .
clearly we must satisfy @xmath14 where @xmath15 is the size of the system . in the cosmological context
, @xmath15 would be the hubble scale , and in order to deal with homogeneous and isotropic effective descriptions we would also have to satisfy @xmath16 where the homogeneity scale @xmath17 as mentioned earlier can be as large as @xmath6 mpc .
the physical quantities we wish to deal with , such as the matter density , pressure , fluid velocity field , etc . are defined by suitably averaging or _
coarse graining _ over the physically infinitesimal volumes ( hence the terminology of coarse - graining cells ( cgcs ) ) .
for example , the matter density of a given cgc is defined as the mass contained in it divided by its physical volume .
the ( effectively ) continuum matter density is built up by repeating this for every cgc in the fluid .
since we deal with the notion of a continuum fluid in cosmology , clearly we need a working description of coarse graining appropriate in this context . however , the underlying theory being general relativity , there are problems .
the crucial issue is that along with the matter distribution ( which is now tensorial ) , the _ metric _ for the spacetime also needs to be coarse grained . as discussed by ellis @xcite , in defining a smoothing operation to go from the scale of stars ( scale 1 in ellis terminology , _ microscopic _ in ours ) to cosmological scales where homogeneity can be assumed to have set in ( scale 5 for ellis , _ macroscopic _ for us )
, one would employ two logically different smoothing operators , one each on the metric and the matter tensor .
we denote the smoothing operator acting on the metric as @xmath18 and the one acting on the matter tensor as @xmath19 ( our notation differs from that used in ref .
we assume that einstein s equations are valid on the microscopic scale of stars , which is reasonable since it is at these scales that solar system tests and binary pulsar studies probe general relativity .
symbolically then ( see fig .
[ fig1 ] ) , the metric @xmath20 at this scale is used to construct the einstein tensor @xmath21 which is equated to the matter tensor @xmath22 , giving einstein s equations @xmath23 .        to obtain the corresponding description at the macroscopic scale , we use the smoothing operators to get a new metric tensor @xmath24 and a matter tensor @xmath25 .
the macroscopic metric @xmath26 is used to construct the einstein tensor @xmath27 at the macroscopic scale .
formally , one can think of the map @xmath18 as inducing a map @xmath28 which acts on the microscopic einstein tensor @xmath21 to yield @xmath29 .
since the operators @xmath28 and @xmath19 are in general different , we have @xmath30 . hence ,
if the einstein equations are valid on a given length scale , then they are in general not valid on some other length scale connected to the first by appropriately defined smoothing operators .
if we can re - express the new einstein tensor @xmath27 in terms of quantities obtained by acting the operator @xmath19 on _ both sides _ of the einstein equations at the microscopic scale , then it should be possible to write @xmath31 where @xmath32 represents corrections to the einstein equations arising from an explicit smoothing . in the standard approach
, it is assumed that these corrections are small enough to be neglected .
( it should be emphasized that given the recent evidence for voids in the universe , one or more additional scales may eventually have to be included in the averaging scheme . )
added to all this is the complication that smoothing or averaging operations are notoriously difficult to define for tensor quantities in a gauge covariant manner .
some examples of averaging techniques employed in general relativity ( and in cosmology in particular ) can be found in ref .
@xcite . here
we restrict ourselves to the formalism developed by buchert @xcite for volume averages of scalar quantities in a given choice of time - slicing . in the following sections we recall this formalism and argue that it is naturally adapted to the notion of the smoothing operator for the _ matter _ sector .
we will reserve the term `` averaging '' for the matter smoothing defined in this way . as for the smoothing of the metric ( which we will refer to as `` coarse graining '' ) , we will take recourse to cosmological observations and argue that the result of such an operation on the inhomogeneous metric must yield the flrw metric on a suitable coarse graining scale .
for simplicity we restrict our attention to a pressure - less fluid or `` dust '' , although a generalization to fluids with non - zero pressure is also possible ( see ref .
@xcite ) . for a general spacetime containing irrotational dust ,
the metric can be written in synchronous and comoving gauge .
] , @xmath33 the expansion tensor @xmath34 is given by @xmath35 where the dot refers to a derivative with respect to ( proper ) time @xmath36 .
the traceless symmetric shear tensor is defined as @xmath37 where @xmath38 is the expansion scalar .
the einstein equations can be split @xcite into a set of scalar equations and a set of vector and traceless tensor equations .
the scalar equations are the hamiltonian constraint and the evolution equation for @xmath39 ,    @xmath40    @xmath41    [ avg2 ]    where @xmath42 is the ricci scalar of the 3-dimensional hypersurface of constant @xmath36 , @xmath43 is the rate of shear defined by @xmath44 and @xmath45 is the matter density of the dust .
note that all quantities in eqns . generically depend on both position @xmath46 and time @xmath36 .
eqns . and can be combined to give raychaudhuri s equation @xmath47 the continuity equation @xmath48 which gives the evolution of @xmath45 , is consistent with eqns . , .
we only consider the scalar equations , since the spatial average of a scalar quantity can be defined in a gauge covariant manner within a given foliation of space - time . for the space - time described by , the spatial average of a scalar @xmath49 over a _ comoving _
domain @xmath50 at time @xmath36 is defined by @xmath51 where @xmath52 is the determinant of the 3-metric @xmath53 and @xmath54 is the volume of the comoving domain given by @xmath55 .
spatial averaging is , by definition , not generally covariant .
thus the choice of foliation is relevant , and should be motivated on physical grounds . in the context of cosmology , averaging over freely - falling observers is a natural choice , especially when one intends to compare the results with standard flrw cosmology .
following the definition the following commutation relation then holds @xcite @xmath56 which yields for the expansion scalar @xmath39 @xmath57 introducing the dimensionless scale factor @xmath58 normalized by the volume of the domain @xmath50  at some initial time @xmath59 , we can average the scalar einstein equations , and the continuity equation to obtain @xcite    @xmath60    @xmath61    @xmath62    [ avg7 ]    here @xmath63 , the average of the spatial ricci scalar @xmath42 , is a domain dependent spatial constant .
the ` kinematical backreaction ' @xmath64 is given by @xmath65 and is also a spatial constant over the domain @xmath50 .
a necessary condition for to integrate to takes the form of the following differential equation involving @xmath64 and @xmath63 , @xmath66 [ avg ]
the averager defined in eqn . is a natural candidate for the smoothing operator @xmath67 which must act on the matter tensor , although we have not yet specified the domain size .
note that this averaging is defined on an inhomogeneous underlying geometry , i.e. it does not smoothen out the geometry itself .
it remains to determine the operator @xmath18 that will coarse grain the metric .
although we do not have an explicit definition for this operator , we can define a notional smoothing operator acting on the metric , by using physical inputs from the universe we observe .
namely , we propose that since the universe has a matter content which is homogeneous on large scales , and since the cmb is highly isotropic around us , this notional operator must reflect these observations by acting on the general metric of eqn . to yield a flrw metric .
crucially , the _ scale factor _ of this flrw metric must ( within our assumptions ) be given by buchert s functional @xmath68 , where the domain size to be used will be specified shortly .
the argument leading to this proposal is as follows .
we begin by making precise the notion of `` homogeneity on large scales '' , by requiring that in the range @xmath69 of cgc sizes allowed by the fluid approximation validity criterion , there exist a subset @xmath70 satisfying the following condition : when the microscopic matter density @xmath71 is averaged over a length scale @xmath72 , the resulting smoothed density _ must lose all information _ about the particular cgc over which the smoothing was performed .
in other words , for @xmath72 , at any given time @xmath36 , the smoothed matter density must have the _ same _ value in all cgcs .
it is in this sense that the smoothed matter density is homogeneous ( see however , the discussion at the end of this section concerning an explicit smoothing of the metric ) .
this definition can be related to the concept of _ statistical _ homogeneity if we assume that the total volume of the observable universe is large enough that each cgc effectively samples one of an ensemble of configurations of matter distribution ( see ref .
@xcite for a careful discussion of this point ) .
it is possible that homogeneity sets in at a scale too large to be considered as the size of the cgcs , namely that the range @xmath73 is not contained in the range @xmath74 and that @xmath75 is not satisfied .
we will assume that this does not happen ( for a more detailed discussion on this assumption , see sec .
2.1 of ref .
we also make the crucial assumption that @xmath73 does not change with time . this will later allow us to consistently analyze light propagation using the effective flrw metric .
this assumption needs some justification and we will return to this issue in the discussion .
now define a coarse - graining cell to be an averaging domain @xmath50as used in buchert s formalism , with a characteristic comoving size @xmath72 .
the averaged density is then @xmath76 , independent of the location of the particular cgc being considered .
it also follows that the averaged expansion scalar @xmath77 is homogeneous in the above sense ( using eqn . ) and describes the evolution of the physical volume of any given cgc .
thus , @xmath78 plays the role of the expansion scalar for a geometry whose infinitesimal spatial volume elements probe length scales of the size of our cgcs ( see however the discussion at the end of this section ) .
it is therefore natural to demand that @xmath78 be the expansion scalar for our notionally smoothed out macroscopic metric .
this is equivalent to demanding that the scale factor of this homogeneous and isotropic effective metric be identified with the functional @xmath68 ( up to a constant factor which we choose to be unity by an appropriate choice of units for the effective comoving coordinates ) .
we are thus left with an effective metric whose line element can be written as @xmath79 which probes comoving length scales of size @xmath12 as chosen above .
the 3-ricci scalar of this effective metric is given by @xmath80 and should be thought of as arising from the coarse graining of the geometry .
following ellis ideas from sec .
[ cgcpic ] then , @xmath81 is in general different from the averaged 3-ricci scalar @xmath63 which arises from applying the matter averager to the inhomogeneous scalar curvature ( @xmath82 ) . the constant @xmath83 has the usual interpretation as a parameter characterizing all possible homogeneous and isotropic 3-metrics
. we can write @xmath84  in terms of @xmath85 by defining @xmath86 .
so we are merely shifting focus from ( the _ a priori _ unknown ) @xmath84  to another unknown @xmath87 .
as the notation suggests , we think of @xmath87  as being a correction ( a part of @xmath32 in the notation of sec .
[ cgcpic ] , see eqns . below ) to the 3-ricci scalar of the flrw metric , although this correction need not be small compared to @xmath85 and need not evolve proportional to @xmath88 .
we choose not to normalize the magnitude of the constant @xmath89 , retaining instead the normalization of the functional @xmath68 at the initial time @xmath59 .
the evolution equations appropriate at this scale are buchert s equations and together with the integrability condition . for completeness , and in the spirit of ellis arguments we can rewrite equations and as follows    @xmath90    [ smop2 ]    where we have symbolically denoted the specific components of the general modified equations to which the various terms of buchert s equations belong .
the integrability condition which supplements these equations is purely a consequence of having non - zero corrections @xmath32 to the standard flrw equations . as a quick check
, we note that requiring the corrections @xmath32 to vanish implies both @xmath91 ( from eqn . ) and @xmath92 , i.e. @xmath93 ( from eqn . ) ; the integrability condition is identically satisfied and we recover the usual flrw solution .    in the arguments leading to eqns .
, we have ignored the effects of a virtual _ change _ in the averaging length scale @xmath12 .
considerable effort has been spent by buchert and carfora @xcite in studying such effects . in particular , in the spirit of the real space renormalization group formalism
, these authors derive a novel _ curvature backreaction _
effect analogous to the kinematical backreaction @xmath94 , together with a volume scaling effect .
the inclusion of these effects would change the evolution equations satisfied by the scale factor @xmath95 .
we also note that buchert and carfora s techniques ( involving ricci flow of the inhomogeneous geometry into a smooth geometry ) form one candidate for an explicit realization of the metric coarse graining operator @xmath18 .
seen in this light , it is apparent that by accounting for the above effects , we would also have to appropriately modify our notion of large scale homogeneity ( to include the volume scaling effect and the residual information from the curvature fluctuations ) . in the arguments leading to the effective metric , the scale factor @xmath95
would be replaced by an appropriately `` dressed '' quantity characterizing the evolution of effectively infinitesimal volumes in the smoothed geometry .
however , the main arguments leading to an effective metric which is flrw , will not be altered given the existence of the range @xmath96 .
for simplicity , we ignore these additional effects in this paper , and hope to investigate them in future work .
it should be noted that a notion of coarse - graining is always implicit in standard flrw cosmology .
it is assumed that a coarse - graining ( while satisfying the conditions of the fluid approximation ) is possible on a length scale which coincides with the scale at which homogeneity is achieved ( in other words , the existence of @xmath96 is assumed ) .
standard cosmology however ignores the possible modifications to flrw cosmology brought about by the process of explicitly smoothing over the inhomogeneities - an issue which buchert s approach addresses .
our key proposal in this paper is that the size of the averaging domain should be taken to coincide with the scale of homogeneity , and that ( up to some corrections ) the quantity @xmath68 corresponds to the scale factor of a template flrw metric .
we now illustrate how our proposal can facilitate the comparison of the predictions of the spatially averaged cosmology with observations , for example by computing the @xmath7  relation using the scale factor @xmath68 , which has a modified time evolution compared to the corresponding scale factor in standard cosmology .
in our approach , we take the effective or template metric to be flrw , and since we are holding the size @xmath12 of the averaging domain fixed , there is nothing conceptually new to be done while studying its null geodesics . the important difference , of course , is that the evolution of the scale factor is now given by the modified friedmann equations ( subject to the caveats listed at the end of the previous section ) .
we note that the use of the flrw metric is , strictly speaking `` unrealistic '' in the sense that in the real universe , light travels mostly in vacuum which has a non - zero weyl tensor whereas in the flrw cosmologies which are conformally flat , light travels through a geometry with vanishing weyl tensor ( see also refs . @xcite and @xcite ) .
this issue is clearly relevant in the standard approach as well .
while we have not addressed this issue in the context of averaging , we emphasize that our approach must be thought of as a first step in gauging the effects of inhomogeneities , by using an _ explicitly _ averaged construction of the template flrw model to be used .
the arguments relating the redshift of a source to the scale factor @xmath95  and to its comoving distance @xmath97 from the observer , are exactly the same as in the standard approach . in particular
, we have for the redshift @xmath98 of a source that emits light at time @xmath36 which is received by an observer at time will indicate the value of the quantity at the present epoch @xmath99 , and the subscript @xmath100 , the value at the initial epoch @xmath59 . ]
@xmath99 , @xmath101 and the comoving distance @xmath97 to the source is determined by solving the equation @xmath102 a caveat to the above relations is that the `` step - size '' used in computing the integrals involved must be understood to be no smaller than the size @xmath12 of the cgcs .
further , the entire construction must be interpreted as a `` fitting template '' that allows us to compare observations with theoretical predictions in a straight forward manner . with our choice of conventions ,
the magnitude of @xmath89 is a dynamically determined parameter of the model and determines the form of the luminosity distance @xmath7  as follows .
solving eqn . for @xmath103 gives @xmath104 where @xmath105 is the hubble parameter and the function @xmath106is defined as @xmath107 in order to complete the picture we need the functional dependence of the hubble parameter @xmath108  on the scale factor @xmath95  .
this requirement is complicated by the fact that the system of equations , , and is only consistent , it does not close . in order to obtain the required relation
then , it is necessary to make certain assumptions about the evolution of the kinematical backreaction @xmath64 and the averaged 3-ricci scalar @xmath84 .
_ @xcite have analyzed a class of scaling solutions of the form @xmath109 the case @xmath110 , @xmath111 is the only scaling solution with @xmath112 . for the cases
@xmath113 , the backreaction is proportional to the averaged 3-curvature should not be confused with the comoving coordinate in the effective metric . ] , @xmath114 the forms of the integrability condition and the solutions indicate that the scaling solutions can be superposed to obtain new solutions .
we will use two such superpositions to construct models which contain an accelerating scale factor @xmath95  and which also yield analytic expressions for the luminosity distance @xmath7 .
the solutions we consider are of the form @xmath115 in keeping with our earlier discussion of the quantity @xmath87 .
namely , we model @xmath87  by a single scaling solution . here
@xmath97 and @xmath116 are related as in eqn . , and @xmath89 and @xmath117 are constants with dimensions of ( length)@xmath118 .
we are interested in models which admit an accelerating scale factor @xmath95 , and will hence assume @xmath119 and @xmath120 ( the reason for which will be apparent shortly ) . inserting the relations in eqns .
and we find    @xmath121    @xmath122      although both models considered above must be regarded as unphysical ( since the luminosity distance of a source is a measurable quantity ) , consistent models also exist .
for example , the case @xmath156 , @xmath157 corresponding to a constant backreaction , mimics the cosmological constant of the standard cosmology .
further , it may be argued that current observations seem to indicate that we live in a universe with @xmath158 and hence the above examples are not observationally relevant in any case . to show that the formalism does admit consistent ( accelerating ) models with @xmath158 as well , we now give an example of such a model
. +   + _ case 3 _ : @xmath159 , @xmath158 . in this case , although the _ effective _ spatial curvature @xmath160 is zero , the average of the _ physical _ spatial curvature @xmath63 _ does not _ vanish , and is given by ( see eqn . ) @xmath161 .
if we take @xmath119 , which guarantees that the scale factor will accelerate for @xmath162 ( see eqn . ) , then we have @xmath163 . the quantity @xmath136 is now given by @xmath164 since @xmath165 and @xmath166 are positive , this integral is well behaved for all @xmath98 .
further , since @xmath158 , we have @xmath167 . defining @xmath168 and noting that @xmath169 and @xmath170 , we get @xmath171^{1/2 } } }    \ , .
\label{lumin16 - 2}\ ] ] in this notation , which is reminiscent of that in the standard lcdm model , if we use eqn .
evaluated at present time we also see that @xmath172 , and that acceleration begins at a redshift given by @xmath173    as an example , in fig . [ fig2 ] we have plotted the behaviour of @xmath174 in units of @xmath175 in this model ( for redshifts between @xmath176 and @xmath177 ) , setting @xmath178 , @xmath179 . in this case
the acceleration redshift is @xmath180 . for comparison
we also show the @xmath7  curve in the standard lcdm model with @xmath181 and @xmath182 , for which the acceleration redshift is @xmath183 .
[ interestingly , a model in which the backreaction @xmath64 behaves as @xmath184 at early times , has recently been studied by li and schwarz @xcite in the context of spatial averaging of a perturbed flrw cosmological model . ]
+   + we have therefore demonstrated the existence of models of the backreaction which can admit an accelerating scale factor .
it is important to keep in mind though , that the models considered above were not physically well - motivated , they served only to demonstrate the construction of luminosity distances using the modified dynamics .
a more realistic scenario would probably contain a backreaction which becomes significant only at late times , when structure formation has advanced sufficiently far .
it will be an interesting exercise to numerically construct luminosity distances for such models and compare them with those of the concordance model .
we have argued that for the purposes of cosmological study , an explicit coarse graining of the matter density in the universe must be performed using the formalism of spatial averaging developed by buchert . due to the observed large scale homogeneity of the matter density and the isotropy of the cmb radiation ,
the notionally smoothed out template metric obtained after this coarse graining must be described by a flrw line element , with the scale factor being given by buchert s functional @xmath68 .
the procedure for constructing luminosity distances is therefore identical to that in the standard approach , except that the _ evolution _ of the effective scale factor must be given by the modified flrw equations supplemented by the integrability condition .
( the scale factor and its evolution will be further modified by including perfect fluid sources with non - zero pressure @xcite and the curvature backreaction and volume effects @xcite . )
such a construction of luminosity distances must be interpreted as a convenient way of comparing observations with theoretical models , while retaining the essence of the physical inhomogeneous geometry .
it will be an interesting project to compare such luminosity distances in realistic backreaction scenarios with current supernovae observations in an attempt to constrain the magnitude and evolution of the backreaction .
we end by discussing some additional issues concerning an explicit smoothing of the metric which we avoided in the present treatment , and also some potential problems with our construction .
firstly , we return to the issue of holding the comoving range @xmath73 , and in particular the smoothing scale @xmath12 , constant .
this is justified only so long as the condition @xmath185 is satisfied , where @xmath186 is the appropriate comoving scale associated with the observable universe at redshift @xmath98 . for example , if we choose a scale @xmath12 which satisfies this condition at the present epoch , the same scale will not in general be valid at the last scattering epoch of the cmb .
however , it seems reasonable that for a small enough redshift range , the condition @xmath187 can be satisfied with a single constant @xmath12 .
for epochs sufficiently far back in the past , the scale @xmath12 would have to be reduced , and the construction of the @xmath7  relation would no longer be valid .
further , a change in the size of the averaging domain would in general affect both the scale factor @xmath95  and the `` constant '' @xmath89  in a non - trivial way .
in particular , any change in @xmath89  would render the effective flrw metric meaningless .
also , a change in @xmath89  need not necessarily be caused by a change in the averaging domain size .
it seems likely @xcite that an explicit smoothing of the geometry using buchert and carfora s technique will yield a 3-space of different spatially constant curvatures at different times , this difference being independent of any change in the averaging domain size with time and being possibly incommensurate with our assumption of @xmath188 .
this could be interpreted in terms of a change in @xmath89  with time , and in such a case , one would be dealing with a _ different _ flrw template at different times .
the effects of such a construction are not clear at present
. clearly this issue needs more careful study , and we hope to address this in the future .
of course , if @xmath89  is made a function of time , care must be taken that the matter content is now so chosen that einstein s equations continue to be satisfied , as pointed out by rsnen @xcite in response to @xcite .
also , as pointed out by rsnen @xcite , the use of a homogeneous and isotropic flrw template will necessarily ignore physical effects such as a non - trivial shear in the underlying geometry .
a potential problem with our template construction is that effects of the non - trivial scalar curvature @xmath189 on light propagation are also simplified to a proportionality to @xmath88 .
although the full scalar curvature will show up in the evolution of @xmath68 via its spatial average @xmath84 , it is possible that this construction may inaccurately model the real universe .
a study of the evolution of the effective spatial curvature as mentioned in the preceding paragraph may help resolve this issue .
there is also the related issue of degeneracy , namely that several underlying inhomogeneous models may reproduce the same backreaction effects when averaged .
for example , it was shown @xcite in the context of the ltb solutions , that models with a vanishing spatial curvature ( the so called marginally bound models ) in the inhomogeneous geometry also have a vanishing backreaction @xmath94 after averaging .
one is now faced with the situation wherein interpreting say luminosity distance observations within the inhomogeneous geometry may distinguish two inhomogeneous models , but interpreting the data _ after _ averaging would not make this distinction and would hence perhaps rule out _ both _ models .
this issue has been recently highlighted by enqvist and mattsson ( the last paper in ref .
@xcite ) who show that a marginally bound ltb model can fit the data from type ia sne without any spatial averaging . since the standard flrw matter dominated model is also ( a special , homogeneous case of ) a marginal ltb model which _ does nt _ fit the sne data @xcite , the formalism after spatial averaging would rule out both models by demanding a non - trivial backreaction instead . however , since the ltb ( and in fact any inhomogeneous ) models themselves face a similar degeneracy issue ( for example several different ltb models can currently claim to fit the sne data @xcite ) , this matter can not be considered as settled , and also deserves further attention .
lastly , we note that buchert s averaging scheme itself has been criticized in the literature ( see , e.g. , ref .
@xcite ) on the grounds that the effects of averaging in a noncovariant manner as done in this scheme , are likely to be gauge artifacts .
it should be noted that within second - order perturbation theory , gauge invariance of the averaged equations has been demonstrated in @xcite and @xcite .
it is also possible to approach the problem of determining modifications to the standard cosmological equations via a covariant method such as zalaletdinov s macroscopic gravity @xcite for example .
such an attempt has in fact been made recently and has led to effective cosmological equations very similar to those of buchert ( see ref .
the proposal in the present paper however , is attractive to us due to its simplicity in implementation , and may be treated as a phenomenological model of the more rigorous equations derived in ref .
[ discuss ]
we have greatly benefited from our correspondence with thomas buchert and syksy rsnen , and it is a pleasure to acknowledge their insightful comments and criticism of an earlier draft .
ap also thanks rakesh tibrewala for useful and stimulating discussions .    1 s weinberg , _ gravitation and cosmology _ , john wiley and sons , inc . , 1972 ; + t padmanabhan , _ theoretical astrophysics , vol.iii_ , cambridge univ .
press , 2002 ; + m trodden and s carroll , _ tasi lectures , introduction to cosmology _ , [ astro - ph/0401547 ] , ( 2004 ) .
f hoyle and m s vogeley , _ astrophys . j. _ , * 566 * , 641 , ( 2002 ) [ astro - ph/0109357 ] ; _ astrophys . j. _ , * 607 * , 751 , ( 2004 ) [ astro - ph/0312533 ] ; + m s vogeley _
et al . _ , in _ proceedings iau colloquium no .
195 _ , ( 2004 ) , a. diaferio , ed . + s g patiri , _ et al .
_ , _ mon . not
_ , * 369 * , 335 , ( 2006 ) [ astro - ph/0506668 ] .
_ , _ astrophys .
* 642 * , 54 , 2005 [ astro - ph/0411197 ] ; +
pietronero l and labini f s , _
proc . _ * 822 * , 294 , 2006 [ astro - ph/0406202 ] ; + joyce m , _ et al .
_ , _ astron . &
_ , * 443 * , 11 , 2005 [ astro - ph/0501583 ] ; + kerscher m , schmalzing j , buchert t and wagner h , _ astron . &
_ , * 333 * , 1 , 1998 [ astro - ph/9704028 ] ; + kerscher m , _ et al .
_ , _ astron . &
. _ , * 373 * , 1 , 2001 [ astro - ph/0101238 ] ; + hikage c _ et al . _
( sdss collaboration ) , _ publ .
_ , * 55 * , 911 , 2003 [ astro - ph/0304455 ] .
+ see also ref .
@xcite .
c c dyer and r c roeder , _ astrophys .
j. _ , * 189 * , 167 , ( 1974 ) ; + r c roeder , _ astrophys . j. _ , * 196 * , 671 , ( 1975 ) ; + m sasaki , _ mon . not .
_ , * 228 * , 653 , ( 1987 ) . + an attempt to use this formalism in the spatial averaging context has been made by + d palle , _
nuovo cim .
_ , * 117b * , 687 , ( 2002 ) [ astro - ph/0205462 ] .
n mustapha , c hellaby , and g f r ellis , _ mon . not
r. astron .
_ , * 292 * , 817 , ( 1997 ) [ gr - qc/9808079 ] ; + m
n celerir , _ astron . and astrophys .
_ , * 353 * , 63 , ( 2000 ) [ astro - ph/9907206 ] ; + h iguchi , t nakamura and k nakao , _ prog . of theo
_ , * 108 * , 809 , ( 2002 ) [ astro - ph/0112419 ] ; + j w moffat , _ jcap _ , * 10 * 012 , ( 2005 ) [ astro - ph/0502110 ] ; + h alnes , m amazguioui and o gron , _ phys .
d _ * 73 * , 083519 , ( 2006 ) [ astro - ph/0512006 ] ; + r mansouri , [ astro - ph/0512605 ] , ( 2005 ) ; + r a vanderveld , e e flanagan and i wasserman , _ phys .
d _ * 74 * , 023506 , ( 2006 ) [ astro - ph/0602476 ] ; + d garfinkle , _ class .
* 23 * , 4811 , ( 2006 ) [ gr - qc/0605088 ] ; + k enqvist and t mattsson , [ astro - ph/0609120 ] , ( 2006 ) .
t futamase , _ phys .
_ , * 61 * , 2175 , ( 1988 ) ; + m kasai , _ phys .
_ , * 69 * , 2330 , ( 1992 ) ; + t w noonan , _ gen .
rel . and grav .
_ , * 16 * , 1103 , ( 1984 ) ; + r m zalaletdinov , _ gen . rel . and grav . _ ,
* 24 * , 1015 , ( 1992 ) ; _ gen .
rel . and grav .
_ , * 25 * , 673 , ( 1993 ) ; + j p boersma _ phys .
_ * d57 * , 798 , ( 1998 ) [ gr - qc/9711057 ] ; + a a coley , n pelavas and r m zalaletdinov , _ phys . rev
_ , * 95 * , 151102 , ( 2005 ) [ gr - qc/0504115 ] ;   a a coley and n pelavas , _ phys . rev .
_ * d74 * , 087301 , ( 2006 ) [ astro - ph/0606535 ] .
+ several other references and a discussion on the issue of averaging in cosmology can be found in - a krasiski , _ inhomogeneous cosmological models _ , cambridge univ .
press , 1997 .
t buchert and m carfora , _ class .
_ , * 19 * , 6109 , ( 2002 ) [ gr - qc/0210037 ] ; + t buchert and m carfora , _ phys . rev .
_ , * 90 * , 031101 ( 2003 ) [ gr - qc/0210045 ] ; + t buchert and m carfora , _ proc .
of 12th workshop on general relativity and gravitation ( jgrg 12 ) , tokyo , japan _ , ( 2002 ) , [ astro - ph/0312621 ] ; + m carfora and k piotrkowska , _ phys .
d _ , * 52 * , 4393 , ( 1995 ) [ gr - qc/9502021 ] ;


the management of heat in dielectric beams and membranes is an important part of the design in micro - electro - mechanical systems ( mems )  @xcite .
while the physics of thermal transport in dielectric materials is well understood , the effect of fabrication processes on the surface physics is less clear and widely treated as a hidden variable in the evaluation of the device performance .
this uncertainty can prolong the design - test cycle , where a parameter such as beam geometry is iterated until the target thermal conductance @xmath3 is obtained  @xcite .    in our application ,
the mems is a cryogenic detector known as the transition - edge sensor ( tes )  @xcite .
figure  [ fig : q - band - detector ] shows a tes designed for measurements of the polarization state of the cosmic microwave background  @xcite .
the tes is comprised of a superconducting moau bilayer deposited on a single - crystal silicon substrate .
the silicon is etched around the bilayer and other electrical components to form a membrane that is supported by dielectric beams .
the total thermal conductance @xmath4 of the beams is a critical parameter that determines the sensitivity ( @xmath5 ) , response time ( @xmath6 ) , and saturation power of the tes detector ( @xmath7 ) ( saturation energy in the case of a tes micro - calorimeter )  @xcite .
the response time is also a function of the total heat capacity @xmath8 of the detector , which is the sum of the heat capacity of the silicon membrane and normal metal films .
the thermal conductance from the membrane to the frame is effectively set by the ballistic - dominated thermal conductance of the short beam ( magnified ) .
the long silicon beams are @xmath10 @xmath1m@xmath9 , and merely support tes bias and rf leads .
the thermal conductance of the long silicon beams is diffusive - ballistic and sub - dominant to the conductance of the ballistic beam .
the pd layer deposited on the membrane sets the detector heat capacity .
, width=321 ]    at temperatures far below the debye temperature of the lattice ( @xmath11 k ) , anharmonic phonon collisions in the bulk of the dielectric are negligible , and the phonon - mean - free path @xmath0 is restricted by the cross sectional dimensions of the beam  @xcite . in state - of - the - art far - infrared and millimeter tes detectors operating at sub - kelvin temperatures  @xcite , the beam length @xmath12 is much larger than @xmath0 , and the thermal conductance is often described by the fourier law for phonon diffusion , @xmath13 , where @xmath14 and @xmath15 are the width and thickness of the beam , and @xmath16 is the bulk thermal conductivity of the dielectric .
@xmath16 is determined empirically with the fabrication and testing of numerous beams of varying width and length , and the aspect ratio @xmath17 is tuned to achieve the desired conductance . however
, this approach rarely provides a complete view of the thermal physics , and higher order effects such as surface roughness of the beams have been identified as the likely cause of the reported non - uniformity in @xmath3 across detector arrays  @xcite .
a simple analysis can readily show the sensitivity of phonon propagation in long beams ( @xmath18 ) to the details of the surface physics .
for a random gaussian surface , the roughness can be characterized by the ratio of the standard deviation @xmath19 of the surface features to the peak thermal wavelength @xmath20 of the incident phonon field . when @xmath21 , phonons are reflected specularly from the surface  @xcite . as
the wavelength becomes much smaller than the scale size of the surface , @xmath22 , phonons are scattered over all angles , or diffused .
assuming isotropic diffusion , it can be shown that for a rectangular beam , the phonon mean - free - path @xmath0 scales logarithmically with @xmath23 ( @xmath24 for @xmath25  @xcite ) .
hence , in a beam where @xmath18 , phonons experience a large number of scattering events , @xmath26 , and the heat flux entering the beam at one end @xmath27 is reduced to @xmath28 at the exit port , where @xmath29 is the diffuse fraction .
the derivative @xmath30 is a measure of the sensitivity of @xmath3 to changes in @xmath31 and proportional to the power - law function @xmath32 .    in practice
, the statistics of the surface roughness can be heavily biased by the fabrication and is likely non - gaussian .
dielectric beams can be roughened in various forms depending on the compatibility and properties of the etchants with the dielectric material .
preliminary tests can shed light on the roughening steps , however , fabrication processes can be difficult to control in practice .
details such as base pressure and process chamber conditioning during dry - etch steps can vary between batches of devices , especially in shared facilities . as a result ,
the characteristics of the roughness of a single surface can change significantly over spatial scales shorter than the beam length , and between wafers .
furthermore , not all surfaces of a beam are rough to the same degree .
these effects add to the unpredictable variability in @xmath31 and therefore @xmath3 in beams dominated by diffusive thermal transport .
the work in this paper is focused on understanding the effects of fabrication on the thermal conductance in single - crystal silicon beams .
the beam geometries explored span several orders of magnitude in @xmath33 .
the results demonstrate that @xmath3 in short ( @xmath34 @xmath1 m , @xmath35 1 ) uniform beams with ballistic - dominated phonon transport is insensitive to the fabrication conditions and can be realized with precision in devices on different wafers .
this approach is suitable for tes detector arrays in ground - based and air - borne telescopes , where the detector noise equivalent power requirements is greater than 10@xmath36 w/@xmath37 .    throughout this paper , the terms _
diffuse _ and _ ballistic _ are used to describe the scale size of the phonon transport in a silicon beam . in the presence of surface roughness , the thermal transport in a long beam with @xmath38
is dominated by diffuse reflections off the beam walls , and @xmath39 .
for a short beam , ballistic transport dominates and @xmath40 . we emphasize that these limits bound the diffusive - ballistic conductance observed in practice .    in sec .
[ sec : method ] , we briefly describe the device fabrication and readout hardware employed for the low - temperature measurements of @xmath3 and @xmath41 . in sec .
[ sec : model ] , we describe the theory of radiative phonon transport in diffusive - ballistic beams . in sec .
[ sec : diff - ball ] , we present results that illustrate the effect of fabrication on the diffusive - ballistic nature of the thermal conductance . in sec .
[ sec : controlg ] , we describe the superior performance of ballistic beams for the control of @xmath3 of tes detectors .
finally , we demonstrate that because of the low specific heat capacity of single - crystal silicon , the heat capacity of the detectors fabricated on this substrate can be effectively controlled using a normal metal film with high specific heat .
the fabrication of the tes detectors reported in this paper is described in detail in ref .
@xcite , but for completeness , we provide a brief description here .
fabrication starts with a float - zone - refined p - type ( 100 ) boron doped silicon - on - insulator wafer .
heat flow is primarily in the [ 110 ] direction .
the single - crystal silicon layer is 5 @xmath1 m thick , with a room temperature resistivity of 5 k@xmath42-cm and a carrier density of 5@xmath4310@xmath44 @xmath45 .
a moau bilayer is deposited by a combination of dc plasma sputter deposition and electron beam evaporation .
the mo and au are photolithographically patterned to define the superconducting sensor .
the au layer is etched by ion milling and the mo layer is etched by reactive ion etching ( rie ) .
further metallic components are deposited on the device layer to enable the required operation of the pixel at rf frequencies , and to enable on - chip diagnostic capabilities such johnson noise thermometry as needed .
au resistors are deposited by electron beam evaporation and patterned by lift - off . following an in - situ rf ar plasma cleaning of the wafer
, a nb microstrip layer is sputter deposited , photolithographically patterned , and etched in fluorine - based rie chemistry .
the silicon membrane is defined and reactive ion etched down to the buried oxide layer .
the silicon wafer is then bonded to a pyrex handle wafer using a sacrificial wax .
the backside of the silicon wafer is patterned by photolithography and the bulk silicon is etched by inductively - coupled plasma deep reactive - ion etching ( drie ) .
the wax is then dissolved in solvent to release the devices .
measurements were performed in a mechanically - cooled system that reaches 3 k at the coldest stage .
further cooling of the detector stage down to 100 mk was achieved using an adiabatic demagnetization refrigerator .
the devices were mounted and thermally anchored to a gold - plated copper package , and wirebonded to a superconducting quantum interference device ( squid ) cold stage amplifier for johnson noise current measurements .
further amplification was achieved with a squid series array mounted at the 3 k stage .
the temperature of the 100 mk cold stage was read with a calibrated ruo@xmath46 thermometer  @xcite with calibration traceable to the national institute of standards and technology .
the precision of the temperature control was 32 @xmath1k rms at 150 mk .
the room temperature electronics was comprised of a custom low - noise pre - amplifier , feedback electronics for flux - locked - loop operation , commercial low - noise high - gain amplifier and low - pass filter ( srs-650 ) , and commercial 16 bit digital - to - analog electronics ( ni pci-6120 ) for digitization of the amplified feedback signal .
the method of johnson noise ( jn ) thermometry is frequently utilized for the characterization of the thermal properties of thin - films below 1 k  @xcite . at audio frequencies ( 0.5 - 20 khz ) and sub - kelvin temperatures ,
the current noise spectral density of a resistor @xmath47 in series with the squid input coil takes the form @xmath48^{-1}/r$ ] , where @xmath49 is the electron temperature , @xmath50 is the characteristic frequency of the circuit , and @xmath12 is the total inductance of the input coil and stray inductance in the input circuit loop .
this equation directly relates measurements of @xmath51 to @xmath49 . for the test devices that are described in ref .
@xcite , the jn thermometer and heater are 240@xmath43230@xmath430.27 @xmath1m@xmath9 au resistors deposited on top of interdigitated nb leads .
the kapitza conductance at the au - silicon interface and the electron - phonon conductance in the au are approximately 10 nw / k each at 100 mk  @xcite .
these conductances are a factor of 50 greater than the highest silicon beam conductance measured in our devices .
we therefore assume the electron system in the thermometer is in thermal equilibrium with the phonons in the silicon membrane .
the maximum expected temperature gradient across the membranes is less 0.5 mk across the 0.1 - 1 k temperature range .
the random error in a temperature measurement from the jn current can be estimated from the radiometer equation @xmath52  @xcite , where @xmath53 khz is the bandwidth over which the current noise power spectral density was measured and @xmath54 s is the measurement time .
hence , the fractional error throughout the temperature range of our measurements is less than 0.5% .
measurements of the jn current as a function of bath temperature were used to fully calibrate the cryogenic and warm readout amplifiers .
typical readout noise of 15 pa/@xmath37 referenced to the input circuit was achieved .
the jn signal was 30 - 100 pa/@xmath37 in the range 0.1 - 1 k.    to measure @xmath3 of the silicon beams , a small signal dc heater power was applied to the membrane and the increase in temperature was recorded .
the relationship between heater power and membrane temperature at a given bath temperature @xmath55 can be parameterized by the power law function @xmath56 , where @xmath57 is the heater power , @xmath58 is the membrane temperature , @xmath59 is a constant that is proportional to the number of phonon modes , and @xmath60 is the temperature exponent that is indicative of the temperature dependence of the phonon mean - free - path . the first order derivative yields the thermal conductance @xmath61 . in our measurements ,
the change in temperature was small , @xmath6210% , and the approximation @xmath63 was valid .
we have also utilized tes current - voltage curves to determine the thermal conductance of detectors in which normal metal films for jn thermometry were omitted . when a tes is voltage - biased on the normal - to - superconducting phase transition , the silicon membrane temperature is fixed at the critical temperature of the moau bilayer , @xmath64 mk .
the thermal conductance is then determined from a fit to the joule power dissipated in the tes as a function of bath temperature . in this case ,
@xmath65 and @xmath66 , where the parameters @xmath59 and @xmath60 are equivalent to those derived from jn thermometry in the same temperature range .
in building a model of heat flow in a silicon beam , it is important to first identify the thermal baths that define the boundary conditions in the model .
a thermal bath is ideally a blackbody that provides a matched termination for phonons at all wavelengths . in a detector such as that shown in fig .
[ fig : q - band - detector ] , each silicon beam is connected to two thermal baths , the frame and the membrane . the electron system in the pd film is strongly coupled to the phonon system in the membrane . at 150 mk ,
the electron - phonon conductance in the pd is 30 nw / k  @xcite , and the kaptiza conductance between the pd and the silicon membrane is @xmath62 130 nw / k  @xcite . for comparison ,
the highest beam conductance measured is 240 pw / k at 150 mk .
the silicon frame is effectively a blackbody since its lateral dimensions are much larger than the peak thermal wavelength . at 100 mk
, the peak thermal wavelength is @xmath67 @xmath1 m in silicon , where @xmath68 and @xmath69 have the usual definitions and @xmath70 m s@xmath71 is the speed of sound for transverse waves .
a phonon scattered from a beam into the silicon frame has a low probability of coherently backscattering into the beam .
given the 5 @xmath1 m thickness of the silicon layer , and the thermal wavelengths above 100 mk ( @xmath72 @xmath1 m ) , the membrane and the frame can be modeled as three - dimensional multi - moded thermal baths , where the number of modes is given by the planck distribution and is much larger than unity .
the beams behave as three - dimensional multi - moded acoustic waveguides since the widths ( @xmath73 @xmath1 m ) and 5 @xmath1 m thickness are much larger than @xmath74 at 0.1 - 1 k. phonons emitted by the baths couple radiatively to the silicon beams , propagate unimpeded in the bulk of the beams , and scatter at the beam surfaces . when the surfaces of the beams are perfectly smooth , phonons coupled to the beams are scattered elastically without any loss of coherence in the beams .
the conductance in this case is determined by the beam areas and is referred to as ballistic . in the case of rough surfaces ,
phonons scatter diffusively with a large probability of backscatter .
the scattering is elastic but the phonon wavevectors are randomized at the surface .
if the surface features form cavities much larger than @xmath74 , the surface can be thought of as a blackbody .
the conductance is thus boundary - limited and was first described theoretically by casimir  @xcite .
an important assumption in casimir s formalism is the isotropic ( lambertian ) reemission of phonons by the beam surfaces .    in practice ,
the scattered pattern of phonons from a rough surface is strongly dependent on the surface characteristics .
when the roughness can be described statistically , the reflection pattern is averaged over the ensemble of the surface realizations . for a surface with a gaussian correlation function ,
the roughness can be characterized by the correlation length @xmath75 and standard deviation @xmath19 of the surface features above a mean height .
figure  [ fig : beampattern ] illustrates the scattered pattern as a function of @xmath76  @xcite . for a nearly smooth surface , @xmath77 , phonons are largely reflected in the specular direction . when @xmath78 , phonons incident on the surface are diffused in all directions , however , the flux is not lambertian . nonetheless , in numerous theoretical studies of boundary - limited scattering  @xcite , the reflected pattern is decomposed into a diffuse lambertian fraction @xmath29 and a specular fraction @xmath79 , as illustrated in fig .
[ fig : beampattern ] .
the diffuse reflection is modeled as inelastic scattering that redistributes the phonon wavevector , although strictly speaking all reflections are elastic at the surface facets that define the roughness .
the surface can then be thought of as a blackbody with emittance @xmath29 .
hence , a completely lossless but rough surface can be modeled as a lossy surface  @xcite .
this approach is numerically less intensive than calculations of the reflection pattern based on the exact surface profile , or the statistics of the surface  @xcite .    ) and a very rough surface ( @xmath78 ) .
the diffuse fraction of the realistic beam patterns are generally approximated as isotropic ( lambertian ) , which is relatively simple to integrate into models of thermal transport .
, width=321 ]    our model of heat flow in a silicon beam has the following assumptions : ( i ) the ends of the beam constitute thermal baths with constant temperature and unity diffuse fraction .
the baths continue to distances much larger than the beam length as shown in fig .
[ fig : modeldiagram](a ) .
this construct is numerically convenient from the perspective of the heat transport model presented , and is similar to that adopted in ref .
( ii ) the surfaces of the beam reflect phonons with a diffuse fraction @xmath29 that is isotropic , and a specular fraction given by @xmath79 .
( iii ) the gradient of the temperature profile across a beam is linearly dependent on @xmath29 and constant as a function of length .
this assumption is valid in multi - moded beams , where the number of modes is much larger than unity and the diffuse phonon reflection can be treated as diffuse emission  @xcite .
models that iteratively solve for the temperature profile in multi - moded beams show negligible deviations from a constant temperature gradient  @xcite .
when the phonon mean - free - path is much smaller than the beam length , or @xmath29 is close to unity , phonon diffusion is the dominant term of the total heat flux , and the constant temperature gradient is an excellent approximation .
the approximation also holds when @xmath29 approaches zero , or the mean - free - path is comparable to the beam length , since the temperature gradient is linearly dependent on @xmath29 ( the mathematical description of this detail is in the appendix ) .
( iv ) the phonon modes that propagate are two orthogonal transverse modes with a sound speed of @xmath80 m / s , and one dilatational mode with a sound speed of @xmath81 m / s .
each mode has a thermal spectrum .
( v ) phonons propagate unimpeded in the crystalline bulk of the silicon beams , and only scatter at the beam surfaces .
( vi ) the diffuse flux is reemitted according to lambert s law of thermal radiation .
( vii ) all modes scatter diffusively with the same probability @xmath29 , and there is no spectral dependence for @xmath29 .
it is possible to map @xmath82 by fitting measurements of @xmath4 and calculating the peak thermal wavelength as a function of @xmath58  @xcite . for a gaussian surface  @xcite , when @xmath77 , the specular power reflection coefficient is proportional to @xmath83 and the diffuse coefficient to @xmath84 ; when @xmath78 the reflection coefficient is independent of wavelength .
hence , the thermal spectrum can be scattered very differently as temperature changes , and in recognizing the limitations of our model , we do not extrapolate the results to @xmath82 .    )
that extend to distances much larger than @xmath12 .
( b ) each surface element of the beam radiates phonons into all solid angles .
the radiated flux from the infinitesimal area @xmath85 that reaches the shaded beam port is reduced depending on the diffuse fraction of the beam surfaces .
( c ) using the method of images , the path of specularly scattered phonons can be re - imaged to one vector .
the components of the vector map the number of scattering events at each surface of the beam . , width=321 ]    figure  [ fig : modeldiagram](b ) illustrates the radiative model .
the beam is divided into infinitesimal areas @xmath86 that radiate heat into all solid angles .
the phonons are scattered several times before reaching the beam exit port and thermalizing in a bath . during each scattering event ,
the phonon flux is reduced by a factor @xmath87 that is the diffuse fraction of the surface @xmath88 .
the nature of the problem is then to simultaneously calculate the flux that is viewed by @xmath89 from @xmath85 and count the number of scattering events from each surface of the beam . for the latter ,
the complex path of the phonons can be unfolded , or re - imaged , to a single vector that reaches an equivalent beam exit port as shown in fig .
[ fig : modeldiagram](c ) .
the components of the vector then equate to the number of scattering events from each surface .
the details of the radiative model are described in the appendix  [ sec : app : a ] .
we define an equivalent view factor for a rectangular waveguide @xmath90 , \label{eqn : fx - final}\end{aligned}\ ] ] where @xmath91 and @xmath92 are equivalent view factors between parallel and perpendicular surfaces of the beam respectively , and @xmath93 is the contribution of the direct ( unscattered ) flux from the surfaces of the beam to the cross sectional area at @xmath94 ( reference plane ) .
@xmath95 is the diffuse fraction of the surface emitting the flux , and the index @xmath96 runs over all surfaces of the beam .
the index @xmath97 runs over the surfaces perpendicular to @xmath96 . included in the terms @xmath91 and @xmath92 is a factor @xmath98 , where @xmath99 is the total number of scattering events from each surface as described in fig .
[ fig : modeldiagram](c ) . for the baths ,
@xmath100 is obtained by setting @xmath101 in eq .
[ eqn : fx - final ] .
the total heat flux through the beam is @xmath102 , @xmath103dx_i ,   \label{eqn : px0}\end{aligned}\ ] ] where @xmath104 , and @xmath105 is an average group velocity . given the linear temperature profile @xmath106 , the approximation @xmath107 where @xmath58 is the average of the bath temperatures , and the relation @xmath108 , eq .
[ eqn : px0 ] reduces to @xmath109 where @xmath110 is the net diffuse power from the beam surfaces , and @xmath111 is the net specular power that is exchanged between the baths .
the conductance of the beam as a function of temperature when scaled by the maximum ( ballistic ) conductance @xmath112 is @xmath113 in the short beam limit , @xmath114 , and @xmath115 . for non - zero diffuse fractions and as the beam length becomes much larger than the cross - sectional dimensions , @xmath116 , and using @xmath117 where @xmath118 is the volume specific heat of the material , the mean free path takes the form @xmath119 .    when the diffuse fractions of the surfaces of the rectangular beam are different , all terms in the summation over the index @xmath96 in eq .
[ eqn : fx - final ] must be evaluated .
the numerical evaluation of the problem is reduced by a factor of two if there is symmetry in the diffuse fraction of the surfaces , @xmath120 and @xmath121 . for the silicon beams discussed in this paper , the top surface and sidewalls are diffuse , and the bottom surface is specular .
it is important to note that the diffuse fraction of the top and sidewall surfaces are uncorrelated .
the terms @xmath91 and @xmath92 in eq .
[ eqn : fx - final ] are evaluated using monte - carlo integration with 10@xmath122 uniformly distributed random points averaged 10@xmath123 times .
@xmath93 is evaluated with an adaptive quadrature routine .
figure  [ fig : modelsimulation ] illustrates the accuracy of the model compared to an analytic expression for the mean - free - path @xmath0 in a rectangular beam  @xcite . at the diffuse limit ( @xmath39 ) ,
the conductance scales as @xmath124 , and at the ballistic limit , the conductance is independent of length .
the expression @xmath125 is remarkably accurate over the range of beam length explored in fig .
[ fig : modelsimulation ] .
this fact highlights the importance of the phonon mean - free - path , that a single parameter can capture a significant fraction of the heat transport in the ballistic - diffusive limit .
figure  [ fig : g - aol ] shows the thermal conductance of silicon beams as a function of beam area - to - length ratio .
the cause of the variability in conductance ( a factor of @xmath625 for the same beam geometry ) is the roughening of the silicon beam surfaces during fabrication .
the top surface of the beams are primarily roughened during a dry etch step that patterns a nb layer . the plasma chemistry , 75% cf4 , 15% o@xmath126 , etches silicon at a higher etch rate than nb .
the exact processes that roughen the silicon are unclear , but the likely candidates are micro - masking due to non - volatile species in the plasma , and slower etch of native and in - situ nbo@xmath46 in the fluorine chemistry .
the use of an etch stop such as 100 nm thick alo@xmath46 was successfully employed in early devices , however , these detectors exhibited very large uncontrolled heat capacity due to the amorphous nature of the etch stop  @xcite .
test devices were also fabricated with a lift - off process that left the beam surfaces smooth .
the conductance of these beams were higher , suggesting boundary - limited phonon propagation as the cause of the variability in conductance of the beams with rough surfaces .
( excluding the lift - off points ) , and is the simplest interpretation of the data that is intended to merely guide the eye . details of this data are summarized in table  [ tbl : data ] .
, width=321 ]    we have explored the relationship between thermal conductance and surface roughness with an extensive dataset of atomic force microscopy ( afm ) images .
non - destructive afm of a silicon beam surface a millimeter in length is not practical . instead ,
the silicon device layer was imaged on the handle wafer , in close proximity to each device . for a silicon surface with a gaussian roughness , @xmath127 nm and @xmath128 nm depending on the processing conditions .
however , we have found that a large fraction of the roughened surfaces did not satisfy the mathematical assumptions required to apply the theory of phonon scattering by nearly - smooth gaussian surfaces , and the measured thermal conductance differed from expectations of such a theory , even when @xmath19 was much less than the typical phonon wavelength at the temperatures explored ( 0.1 - 1 k ) .
we observed power spectral densities of surface height that were approximately power laws over the range of observable wavevector values . in this case
, we found an empirical relationship between the power - spectral - density @xmath129 of the surface profile and the measured magnitude and temperature dependence of the phonon mean - free - path in our silicon beams , @xmath130 , where @xmath15 is the silicon beam thickness and @xmath131 is the phonon wavevector .
as expected from beckmann s theoretical treatment of scattering from non - gaussian surfaces  @xcite , the observed relationship involves the correlation function of the local tilt angle of the surface facets , @xmath132 .
it is challenging to precisely predict the phonon mean - free - path from information contained in the afm data due to practical sampling considerations .
this occurs for several reasons : the statistics of the surface roughness are observed to vary significantly across the wafer and are sparsely sampled in @xmath133 @xmath1m@xmath123 patches set by the resolution and stability of the afm . in addition , for a multi - moded beam with thickness of order the width , details such as surface roughness of the beam sidewalls are significant , as is shown in sec .
[ sec : controlg ] .
non - destructive afm imaging of the sidewall geometry is extremely difficult to perform , and such a task becomes impractical when considering the number of silicon beams that constitutes a detector imaging array .
however , we have found that afm images can provide qualitative insight to the thermal conductance in beams where the surface roughness relative to the thermal wavelength has a significant amplitude and is non - gaussian .
suggests resonant phonon scattering from the rough surface .
( b ) afm images of 6.4@xmath436.4 @xmath1m@xmath123 patches representative of the top surface of the beams described in ( a ) .
the scale size estimated from the peaks in @xmath134 is consistent with the roughness scales measured from these images .
the rough surface is heavily pitted , with an average pit depth of 30 nm .
the large distribution in pit size and shape is the likely cause of the observed broad resonance .
since the scale size is much smaller for the smooth surfaces ( @xmath135 nm ) , the onset of the peaks occur at a much higher temperature ( shorter thermal wavelengths ) .
the striping in the afm image of the smooth surface is due to residue and accumulation of particulates on the afm tip .
beam dimension are @xmath136 mm , and @xmath137 @xmath1 m for devices labelled with squares , diamonds , and triangles respectively.,width=321 ]    figure  [ fig : cas - lim](a ) shows the conductance of smooth ( lift - off ) and roughened ( dry etch ) silicon beams when compared to the boundary - limited conductance in each case  @xcite .
the phonon propagation transitions from diffuse to ballistic at 300 mk , where the peak thermal wavelength is @xmath138 nm . the minimum in conductance suggests the transition is resonant for the rough beam  @xcite .
the inset in fig .
[ fig : cas - lim](a ) shows the inverse phonon mean free path @xmath134 .
the scale height of the scatterers @xmath139 can be estimated from the temperature at which @xmath134 peaks . using the rayleigh criterion , @xmath140  @xcite , where @xmath74 is the phonon wavelength at the peak temperature , and @xmath141 is the phonon angle of incidence to the surface . assuming an average over all angles , @xmath142
, the rayleigh criterion suggests @xmath143 nm for the rough beam and @xmath144 nm for the smooth beams .
these estimates are consistent with the roughness measured from the afm images shown in figure  [ fig : cas - lim](b ) .
the rough ( dry etch ) surface is heavily pitted , with an average pit depth of 30 nm and average size of 250 nm that is comparable to @xmath74 at 300 mk . the scale size is much smaller for the smooth ( lift - off ) surfaces , @xmath135 nm , which leads to the onset of resonance at a much higher temperature .
in addition to phonon scattering by rough surfaces of a beam , organic residue such as photoresist on the surfaces can enhance diffuse scattering .
thin layers of metals and non - metals hundreds of angstroms thick have been shown to induce diffuse transport in silicon  @xcite . in our devices ,
the reduction and variability in @xmath3 in a few cases were at least qualitatively attributed to residue accumulation on beams of order a millimeter in length .
care must be taken when devices are cleaned since the cleaning process may introduce further residue , or roughness ( see example of silicon dipped in hf in ref .
@xcite ) .    table  [ tbl : data ] summarizes the silicon beam geometries and conductances shown in fig .
[ fig : g - aol ] .
the exponent @xmath145 is inversely proportional to @xmath3 , indicating a transition from diffuse to ballistic transport .
low @xmath3 devices have boundary - limited transport , and @xmath146 .
high @xmath3 devices are approaching ballistic transport , and @xmath147 in the transition region .
.conductance fit parameters @xmath148 and @xmath145 as a function of beam geometry , @xmath149 .
the entries marked with an asterisk indicate beams with smooth top surfaces ( lift - off devices described in the text ) . [ cols="^,^,^,^,^ " , ]     [ tbl : data ]
in the design of prototype tes detectors , the choice of the beam length and width was guided by the data in fig .
[ fig : g - aol ] .
several strategies were adopted in order to control the conductance and produce uniform detectors for the focal plane array described in ref .
@xcite . to minimize the top surface area of the silicon beams roughened by the dry etch ,
the width of the nb leads were extended to cover a larger fraction of this surface , and the beam widths were defined with a combination of drie using the standard bosch process  @xcite followed by rie to remove the silicon down to a nb ground plane .
the two processes are separated to reduce potential metallic contamination in the drie chamber while minimizing lateral silicon undercut from the isotropic rie .
the bottom surface of the beams in contact with the nb ground plane are smooth , and scanning electron microscopy images of the sidewalls showed no significant features from the rie etching of the beams employed in test structures , fig .
[ fig : sidewall - sem](a ) .
the beam sidewalls were predicted to be nearly specular , with conductances comparable to the lift - off devices shown in fig .
[ fig : g - aol ] .
the top surface of the beams was expected to be slightly rough , @xmath150 , since the area exposed during the dry etch step was reduced by a factor of two .
four beams 785 @xmath1 m in length and 13 @xmath1 m in width were integrated into the tes design .
post fabrication , the beams had well - defined top surfaces , but the sidewalls were roughened by the drie steps , fig .
[ fig : sidewall - sem](b ) .
the total conductance was a factor of 5 below the predicted target range . in principle , the scale length and height of the sidewall features that include the drie scalloping as well as passivation residue present a very diffuse surface for phonons above 100 mk .
figure  [ fig : longbeam - model ] compares the total conductance of the four beams to simulations based on the diffusive - ballistic model described in sec .
[ sec : model ] .
the diffuse fractions of the sidewall and top surface of the beams vary across a wafer and between wafers .
the inferred diffuse fraction of the top surface of the beams is @xmath151 .
a possible source for the moderate roughness is the 200 nm thick nb lead sputter - deposited on each beam .
the lead and beam are generally considered as decoupled systems , and the nb conductance is estimated to be far below that of the silicon beam  @xcite .
however , it is the diffusion process of the phonons in the silicon by the nb layer that is of importance and generally unknown .
it has been shown that even in the presence of large acoustic mismatch , phonons in silicon can be diffused at surfaces with metallic layers tens to thousands of angstroms thick ( see fig .  8 in ref .
@xcite for examples of phonon diffusion caused by al and ti films below 1 k ) .
, devices marked with crosses , squares , diamonds , and circles have similar sidewall profiles .
( b ) to minimize the top surface area roughened by the dry etch of the nb lead layer , drie was used to define the beam width
. the ridges at the top of the beam are scalloping from the drie cycles .
possible sources for the bright features on the sidewall are leftover passivation layer , and charged silicon features . ,
width=321 ]    it is evident from the results in fig .
[ fig : longbeam - model ] that a solution to the problem of conductance control is a design that is least sensitive to the effect of fabrication conditions such as surface roughening and residue .
our approach was to reduce the length @xmath12 of a beam , which increases the phonon mean - free - path @xmath0 by reducing the interaction probability of the phonons with the rough sidewalls of the beam .
the conductance is less sensitive to @xmath0 in the limit @xmath152 .
figure  [ fig : shortbeam ] shows the simulated and measured total conductance when a ballistic beam 10 @xmath1 m long and 13 @xmath1 m wide is integrated into the detector as shown in fig .  [ fig : q - band - detector ] . given our knowledge of the surface roughness after the dry etch step ( fig .
[ fig : longbeam - model ] ) , a diffuse fraction @xmath153 was expected and subsequently obtained in practice .
m long beams shown in fig .  [
fig : q - band - detector ] .
the measured conductance was a factor of @xmath625 below the design conductance ( hatched box ) .
the cause of the variability is boundary - limited scattering .
( right ) simulations of the conductance as a function of beam wall diffuse fractions .
the inferred diffuse fraction of the samples are moderately rough top surfaces , @xmath154 , and very rough sidewalls , @xmath155 .
, width=321 ]    .
the shaded and hatched regions are defined by the measured range of conductance of the long diffuse beams shown in fig .
[ fig : longbeam - model ] .
( right ) the measured conductances are within the predicted range of the simulations ( light grey ) .
, width=321 ]    the total conductance , which is largely set by the ballistic beam , was predicted prior to fabrication .
the diffuse beams merely support microstrip and bias leads necessary for the operation of the detector .
however , they contribute the largest fraction to the variance in the total conductance across a wafer and between wafers .
the magnitude and variance in @xmath3 of these beams were estimated from the data in fig .
[ fig : longbeam - model ] .
the fractional standard deviation in total conductance is @xmath28% and within the requirement of the target instrument  @xcite .
the range of conductances that can be accessed with a ballistic beam is suitable for low - noise bolometric detectors operating at ground - based or air - borne telescopes where the signal power is of order 1 - 3 pw / k . in the ideal design , a single ballistic beam sets the conductance of the detector and carries all signal and bias leads .
this approach is very well suited to the fabrication of kilo - pixel arrays , especially when the signal is coupled to the detector with free - space absorbers . for space applications where the signal power is less than 1 pw
, the ballistic beam can be a phononic structure with phonon bandgaps that can be engineered to reduce the conductance  @xcite .
the heat capacity of the dielectric also plays an important role in the performance of thermal detectors . for a bolometer and calorimeter
, @xmath41 affects the detector time constant , and in some cases anomalous heat capacity associated with the dielectric has been reported to increase thermodynamic noise of these detectors  @xcite .    the heat capacity of the silicon membranes were determined from @xmath156 , where the thermal time constant @xmath157 was measured from the decay in temperature of a silicon membrane to a pulse in heater power .
compared to the debye specific heat  @xcite , @xmath158 at 200 mk and 0.7 at 1 k. the systematic difference between the measurements and the debye approximation is due to the large contribution to the total heat capacity of the devices from the normal metal au heater / thermometer and pd - au termination on the silicon membranes .
nonetheless , the measured specific heat are extremely small . for the detectors shown in fig .
[ fig : q - band - detector ] , the estimated silicon heat capacity is @xmath622 fj / k at 150 mk , which is advantageous since the overall detector heat capacity can be controlled with a normal metal that can be reliably deposited on the membrane .
we have successfully used pd to achieve a heat capacity of 2 pj / k  @xcite .
measurements indicate that the specific heat of pd below 1 k is 5.8 mj mol@xmath71 k@xmath159 , which is lower than the room temperature value of 9.4 mol@xmath71 k@xmath159  @xcite .
we have successfully deposited 400 nm pd layers with a lift - off process .
it is important to emphasize that there is no evidence for anomalous specific heat in the single - crystal silicon device layers .
this is in contrast to the large excess specific heat of amorphous sin@xmath46 that is widely used in the fabrication of thermal detectors  @xcite .
the thermal conductance of single - crystal silicon devices can be precisely controlled using a short beam with ballistic - dominated phonon transport .
the conductance of the short beam is largely set by its area and is insensitive to the details of the fabrication conditions that can vary over time .
this approach is superior to the conventional use of long beams for the control of thermal conductance .
phonon transport in long beams is diffusive - ballistic and very sensitive to the detailed surface physics of the beams .
the conductance is thus boundary - limited , and has been observed in beams with modest surface roughness ( 5 - 10 nm rms in surface height ) .
resonant phonon scattering has also been observed at 300 mk in beams with highly pitted surfaces with pit depths of order 30 nm .
boundary - limited scattering has contributed to a large variance , up to a factor of 5 , in the conductance of devices fabricated and tested during this work . in contrast ,
when a beam with ballistic - dominated thermal transport is integrated into the design , the uniformity in conductance is reduced to a fractional deviation of @xmath28% in devices fabricated across two wafers , and this variability is largely set by the conductance of diffuse beams that support electrical and microwave leads to the device . in addition
, we have found no evidence of excess specific heat in single - crystal silicon membranes .
hence , the total heat capacity of the detector can be effectively controlled with a normal metal film . for the transition - edge sensors ( tess )
described in this work , the heat capacity was determined by a pd layer 400 nm thick .
hence , the strategies outlined for the precision control of thermal conductance and heat capacity are well suited to the fabrication of uniform large - format arrays of tess with sensitivities approaching 10@xmath36 w / hz@xmath160 .
we gratefully acknowledge financial support from the nasa roses / apra program .
k. rostem was supported by the nasa postdoctoral program at the goddard space flight center during the initial stages of this work .
we thank samelys rodriguez for wirebonding and cryo - cable support .
we thank our collaborators a. ali , j. appel , and t. a. marriage for their contribution to the data acquisition .
we thank the national institute of standards and technology for providing the squid series array and time division multiplexer chips .
we start from the usual definition of the infinitesimal view factor @xmath161  @xcite , @xmath162 where @xmath163 is the vector between the infinitesimal surfaces @xmath164 and @xmath165 , and @xmath166 are the angles between the normal vectors to @xmath167 and @xmath163 ( see fig .  [
fig : modeldiagram ] ) .
we define an equivalent view factor for convenience , @xmath168 and note that the heat exchanged between the infinitesimal surfaces is @xmath169\,f_{da_i - da_j},\end{aligned}\ ] ] where @xmath170 are the temperatures of the surface elements .
the total heat flux is the integral of @xmath171 over the beam length , @xmath172 where @xmath173 represents the surfaces of the waveguide , which in our case is a rectangular beam .
the equivalent view factor between any two perpendicular surface elements of the beam is @xmath174 ^ 2 } \left [ \,f_m\prod_k ( 1-f_k)^{n_k}\right ] \right\ } \,da_i\,da_j , \\          & = & \mathcal{f}^\perp ( m ; x_i , y_i , y_j , z_j ) \,dx_i\,dy_i\,dy_j\,dz_j , \nonumber   \end{aligned}\ ] ] where the indices @xmath96 and @xmath88 sum over all surfaces of the beam . for parallel surface elements , @xmath175 ^ 2 } \left [ \,f_m\prod_k ( 1-f_k)^{n_k}\right ] \right\ } \,da_i\,da_j , \\          & = & \mathcal{f}^\parallel ( f_m ; x_i , y_i , y_j , z_j ) \,dx_i\,dy_i\,dx_j\,dy_j . \nonumber\end{aligned}\ ] ] @xmath99 is the number of specular reflections experienced by phonons reflecting off the element @xmath89 first , and then all other surfaces until a reference plane ( e.g. @xmath94 in fig .
[ fig : modeldiagram](a ) ) is reached .
@xmath99 can be determined from the geometry of the specular scattering as described in fig .
[ fig : modeldiagram](c ) , and is the sum of scattering events from each surface . for example , if the radiating element is on the top surface , @xmath99 can be found for phonons reflecting off the top and bottom surfaces , and reflecting off the sidewalls only . in this way
, a different diffuse fraction can be assigned to each surface .
it is also useful to note that the linear dependence of the equivalent view factors on @xmath95 can be translated to a linear dependence of the temperature gradient on the diffuse fraction when eq .
[ eqn : perp ] and  [ eqn : para ] are finally substituted into eq .
[ eqn : px01 ] .
the integral equations of the equivalent view factors over the surfaces of the beam are @xmath176 where the limits of the integrals depend on the surface @xmath96 from which the phonons were originally radiated .
since a temperature gradient is only present in the principal direction of heat flow ( x - axis ) , all surface elements at the plane of constant @xmath177 are at the same temperature and radiate the same flux .
hence , the equivalent view factor for the direct ( non - scattered ) heat flux @xmath93 from the surface to the reference plane can be calculated analytically from eq .
[ eqn : vf ] by considering the view factor of a infinitesimal strip of width @xmath14 to the reference plane at @xmath94  @xcite .