solar flares are sudden , brief , and powerful outbursts , releasing an energy of @xmath0 ergs on time scales of a second to several tens of minutes . in solar flares , both protons and electrons are accelerated to non - thermal energies .
although these energetic charged particles are strongly affected by the solar magnetic fields and can not usually be detected , we can utilize neutral secondary particles , neutrons and photons , to probe them . among them , hard x - ray to gamma - ray continua , produced via bremsstrahlung process , provide the best diagnostics of the flare - accelerated electrons .
although high - spatial - resolution images of flares in the hard x - ray and gamma - ray energies have become available by the yohkoh hxt @xcite and rhessi @xcite respectively , we still depend heavily on the spectral data when attempting to study the non - thermal electrons . actually , to calculate back
the electron spectrum from the observed gamma - ray spectra has been an important issue in solar physics .
as we understand , flare - accelerated electrons having a power - law distribution penetrate the solar chromosphere , and produce bremsstrahlung photons via `` thick - target '' emission .
the thick - target emission is a process wherein non - thermal electrons impinge on a thick matter , and lose almost all their energies through repeated coulomb collisions while emitting x - ray / gamma - ray photons via bremsstrahlung .    in the thick - target condition , the spectrum of x - ray emitting electrons is formed by an equilibrium between the injection of newly accelerated electrons and the loss of their energies predominantly through collisions . if the injected electron number spectrum is given by a power - law as @xmath1 ( electrons kev@xmath2 s@xmath2 ) , with @xmath3 the electron kinetic energy and @xmath4 a constant called electron index , then we expect the x - ray photon - number spectrum also to take a power - law from as @xmath5 where @xmath6 is a quantity called photon index . in 1974
, brown found an approximate relation as @xmath7 in non - relativistic regime , assuming that all emitted photons are collected .
the photon spectrum thus becomes relatively flat , because of `` loss - flattening '' effect working on the electrons ; as they penetrate deeper into the target material , the spectrum gradually hardens because lower - energy electrons have shorter time scales of energy loss through coulomb collisions .    besides the thick - target condition , another extreme case called thin - target condition ( brown 1971 ) was studied extensively .
this represents a situation wherein a bunch of energetic electrons pass through a thin medium while emitting bremsstrahlung photons , but they leave the emission region with their energies almost unchanged . in the thin - target condition ,
the energy distribution of the incident electron is not significantly affected by collisions with the target .
the emergent photon spectrum is steeper , because the electrons are free from the loss flattening effect and is characterized as @xcite @xmath8    since then , there have been many attempts to improve equations ( [ eq.brown ] ) and ( [ eq.brownthin ] ) , particularly toward relativistic regime . some utilized purely analytic calculations of various probability distributions , while others employed numerical solutions to the fokker - planck equation , or monte carlo techniques . in electron energies up
to 10 mev , @xcite investigated angular dependence of the bremsstrahlung in the hard x - ray range , and discussed the associated electron transport . @xcite
dealt with the particle transport by employing a simple thick target emission model , assuming an isotropic electron distribution .
@xcite conducted a detailed study of the directional bremsstrahlung emission by anisotropically accelerated electrons .
@xcite made a calculation including pitch angle scattering of electrons by plasma waves .
these previous works have indicated that the gamma - ray continua of solar flares can be generally explained in terms of chromospheric thick - target emission , produced by electrons having a power - law energy distribution with @xmath9 .    in spite of these extensive studies ,
the estimation of electron spectra from flare gamma - ray data would not be self - contained unless we take into account another important physical process ; namely , compton scattering .
the relativistic bremsstrahlung from precipitating electrons must be significantly forward - beamed , so that a majority of gamma - ray photons must be compton back - scattered by the solar materials in order to reach us .
this process will strongly degrade the gamma - ray energies , and affect the photon spectrum to be observed . to treat this kind of physical condition which involve multiple compton scatterings ,
we clearly need to use monte carlo techniques .
the effect of compton scattering has been studied by several authors , beginning , e.g. , with @xcite who used monte - carlo simulations in such early days .
kontar et al .
( 2006 ) considered this issue using a semi - analytic green s function by magdzias & zdziarski ( 1995 ) , and calculated `` primary '' photon spectra from observed ones .
these works have clarified that the compton back - scattering can significantly modify the bremsstrahlung spectra . however , these pioneering works have treated electrons with sub - mev kinetic energies , leaving higher energies less studied .    in the present paper
, we calculate not only the electron transport and bremsstrahlung emission , but also the photon propagation via compton process . for this purpose
, we employ a monte carlo simulation toolkit named geant4 , which is widely used in experimental high energy physics .
our results imply that the compton scattering has significant effects on the observed gamma - ray spectra .
together with full - relativistic effects , this is expected to require revision to the relation such as equation([eq.brown ] ) and equation([eq.brownthin ] ) .
we show the relevance of these results to actual flares , by comparing the simulated gamma - ray spectra with those observed with yohkoh from a sample of solar flares .
the geant4 monte carlo simulation toolkit @xcite , developed in experimental high energy physicists , has four important features .
firstly , it allows a very flexible construction of `` geometry '' in which particles interact .
secondly , it can `` track '' individual particles while monitoring their physical quantities such as the energy , momentum , and position ; the name geant in fact comes from `` geometry and tracking '' .
thirdly , it allows us to incorporate the desired interactions ; in the present case , coulomb scattering and bremsstrahlung for electrons , and compton scattering for photons .
finally , the toolkit has been tested and calibrated extensively by many investigators under different conditions .      in our monte carlo simulation using geant4
, we basically shoot a parallel beam of a large number of electrons , either vertically or slantly , into our `` solar atmosphere '' to be explained later .
the electrons are assumed either to be mono - energetic , or to have a power - law energy distribution with an index @xmath4 . in the latter case ,
the initial electron energies are assumed to be in the range of 1100 mev , so as to simulate gamma - ray flares .
we implemented the following physics processes which we require our particles to obey in our simulation : ionization , multiple scattering , and bremsstrahlung for electrons , while photo - absorption , compton scattering and electron - positron pair creation for photons .
we suppressed the production of secondary electrons for simplicity .
because those secondary electrons ( typically less than 100 kev ) have much lower energies than the primary ones , they can be considered negligible .
one of the essential features of geant4 is step - by - step tracking of each individual particle , considering all the implemented elementary processes .
the particle tracking algorithm roughly consists of the following 4 steps .
1 ) the initial particle velocity is calculated .
2 ) each physical process calculates the mean - free - path of the specified particle under the relevant interaction , and generates a random number around the mean .
then it proposes this random number as a step length .
3 ) among various competing processes , the interaction which proposes the shortest physical length is adopted .
4 ) track properties , such as the kinetic energy and momentum , of the current particle are updated .    among the implemented processes ,
the coulomb scattering of electrons need special consideration , because its mean free path is much smaller than those of the other processes ( e.g. , @xmath10 of that of bremsstrahlung for a 10 mev electron ) . as a result , its full monte carlo treatment would demand too much computing times . therefore , the expected total energy loss of an electron due to multiple coulomb scatterings , in each step which is determined by the processes other than the coulomb scattering itself , is computed based on semi - empirical formulae by @xcite .
fluctuations around the mean energy loss is expressed by generating random numbers after @xcite .
similarly , the total spatial displacement caused by the multiple scattering in each step is calculated based on lewis method @xcite , which solves a diffusion equation of electrons in the matter .
fluctuations around the expected mean is again represented by random numbers @xcite .
thus , geant4 simulates the effect of coulomb interactions after a given step as a statistical sum of a large number of scatterings , instead of treating them one by one .    in this way ,
the injected electrons and the bremsstrahlung - generated photons are tracked , until they escape out of the model boundary , or they become less energetic than a certain threshold energy which we set at 1 kev .
we collect these photons under a specified condition , and make their spectra .
the simulation of flare - accelerated electrons , precipitating onto the chromosphere , requires yet another physical process to be implemented : electron gyration around the solar magnetic fields .
its exact treatment , however , would make the simulation extremely time - consuming , because the gyration radius of a @xmath11 mev electron in a @xmath11 gauss field is as small as @xmath12 m. therefore , we have decided not to take gyration into our simulation , but instead , to consider slant injections ; as illustrated in figure [ fig.coordinate ] , injecting electrons at a polar angle @xmath13 , and collecting the emergent photons at another polar angle @xmath14 but all over the azimuth angle @xmath15 , are equivalent to simulating electrons gyrating with a constant pitch angle @xmath13 , as long as the field lines are close to normal to the photosphere .
the case of tilted magnetic fields can be reproduced by a superposition of different injection angles , @xmath13 .      before actually simulating the electron transport and gamma - ray emission employing a realistic model for the solar atmosphere
, we carried out simpler simulations with two purposes in mind ; to validate our simulation code by examining whether each of the relevant physical processes is correctly reproduced , and to grasp the essence of the physics to be investigated .
we accordingly approximated the sun as a hydrogen gas box of @xmath16 , having a uniform density of 3.2@xmath17 g @xmath18 which is the density at the solar photosphere .
hereafter we call this model `` simple solar atmosphere model '' ; it is one of three uniform target models employed in the present paper .
we vertically injected @xmath19 mono - energetic electrons with an initial kinetic energy of 50 mev into our `` simple solar atmosphere '' . in order to monitor how the coulomb interaction affects the spectrum and angular distribution of the electrons
, we prepared imaginary boundaries at different depths , 100 , 200 , and 300 km from the injection surface , and collected electrons at each boundary .
figure [ fig.monoelectronspectrum ] a shows the electron spectra obtained in this way , while figure [ fig.monoelectronspectrum]b the associated angular distributions . from these , we can see that the maximum energy of the spectrum decreases as the electrons penetrate deeper .
since a 100 km thick slab in our `` solar atmosphere '' has a column density of 3.2 g @xmath18 , an electron of energy @xmath20 50 mev passing through it is expected to lose @xmath20 15 mev , as we can easily calculate from the bethe - bloch formula .
the results in figure [ fig.monoelectronspectrum]a indeed meet this expectation .
the geometrical thickness ( 10,000 km ) of the model is thus large enough to make it physically thick to coulomb loss , for electrons with initial energies up to @xmath201.5 gev .
furthermore , the electron energy distribution is seen to broaden as the energy loss proceeds .
figure [ fig.monoelectronspectrum]b reveals that the injected electrons , initially forming a parallel beam , are gradually deflected due to coulomb interactions with ambient atoms .
the simulated angular distributions approximately consist of two components : a gaussian component caused by small - angle multiple scatterings , and a power - law tail due to large - angle rutherford scatterings .    in a similar way
, we simulated the transport of energetic electrons which has a power - law distribution over the 1 - 100 mev range with @xmath21 .
the spectra and angular distributions , measured at depths of 100 , 200 , 300 , 400 , and 500 km , are presented in figure [ fig.powelectronspectrum ] .
the maximum electron energy ( initially at 100 mev ) thus decreases according to the bethe - bloch formula .
furthermore , the low - energy part of the spectrum gradually flattens , because the energy loss per unit length is roughly independent of the electron energy , and hence lower - energy electrons suffer a larger _ fractional _ energy loss than higher - energy ones .
this is essentially the same as the `` loss flattening '' effect mentioned in section [ sec.introduction ] .
the electron index we employed here , @xmath21 , is significantly smaller ( flatter ) than would be achieved in the standard diffusive shock acceleration process @xcite .
this is motivated by two reasons .
one is that the effects of compton scattering are more clearly observed as @xmath4 gets smaller .
the other is that such a flat electron distribution could be realized in actual solar flares , via , e.g. , direct electric acceleration .
hereafter , we hence consider both flat ( @xmath22 ) and steep ( @xmath23 ) electron distributions .
we investigated photon spectra created by the bremsstrahlung process under the two representative conditions , namely the thick - target and thin - target conditions explained in section [ sec.introduction ] .
the geant4-simulated gamma - ray spectra are compared with analytic formulae , which are fully relativistic unlike the non - relativistic calculations of eq.([eq.brown ] ) and eq.([eq.brownthin ] ) .      to validate the bremsstrahlung process under the thin - target condition which is simpler of the two ,
we let @xmath24 mono - energetic electrons vertically penetrate a second simplified model atmosphere , namely a uniform hydrogen gas with a size of ( 1000 km)@xmath25 , having a density of @xmath26 g @xmath27 ; this is henceforth called `` thin atmosphere model '' .
the model has a column density of 1.0 mg @xmath18 , so that the injected electrons lose their energies by no larger than 1 % , in agreement with the thin - target condition .
we collected the information of each bremsstrahlung photon when it was created , and obtained the spectra shown in figure [ fig.bremsmonoenergyspectrum ] .
also shown in figure [ fig.bremsmonoenergyspectrum ] are analytically calculated spectra , using a fully relativistic formula @xcite , as @xmath28 \right\ } \label{km3bse}\end{aligned}\ ] ] with @xmath29 here , @xmath30 is a cross section for an electron of total energy @xmath31 ( including the rest mass energy ) to interact with a target atom of charge @xmath32 , and emit a bremsstrahlung photon of energy @xmath33 , thus achieving a final energy of @xmath34 .
the classical electron radius is denoted @xmath35 .    as seen in figure [ fig.bremsmonoenergyspectrum ] ,
the difference between our simulation and the analytical calculation is less than 10 % , indicating that the bremsstrahlung process has been correctly implemented . in the fully relativistic regime(@xmath36,@xmath37 ) , in particular , the emitted gamma - ray spectrum has a power - law shape with a photon index @xmath38 , or an energy index of @xmath39 , up to the maximum photon energy that is identical to the initial photon energy .
this is reasonable , because eq.([km3bse ] ) predicts @xmath40 when @xmath41 and hence @xmath42 .
figure [ fig.powerthinspectrum ] shows the bremsstrahlung photon spectra calculated under the same thin - target condition , but when the electrons are assumed to have a power - law spectrum distributed from 1 to 100 mev , with @xmath43 or 1.5 .
the simulated spectra again exhibit approximately power - law like shapes , with a steeper slope than in the case of mono - energetic electrons because softer electrons can obviously emit only lower - energy photons .    in figure [ fig.powerthinspectrum ]
, we also show photon spectra calculated analytically using equation ( [ km3bse ] ) as @xmath44 where @xmath45 is the photon number spectrum , @xmath46 is the injected electron spectrum ( @xmath47 in the present case ) , @xmath48 is the target ion density , and @xmath49 is the velocity of an electron having an energy @xmath50 .
the simulation agrees with the analytic results within @xmath51 .
thus , the geant4 simulation and the fully - relativistic analytic formalism consistently indicate @xmath52 for @xmath53 and @xmath54 for @xmath55 , both measured in a typical gamma - ray energy range of 1 to 10 mev .
these photon indicies are significantly smaller ( flatter ) , particularly for @xmath55 , when compared to the conventional non - relativistic relation of equation ( [ eq.brownthin ] ) .    in this way , we repeated the geant4 simulation with difference values of the electron index , and calculated the bremsstrahlung photon index over an energy of 1 - 10 mev .
the obtained relation between the electron index and the photon index is shown in figure [ fig.thinthick ] , in comparison with the prediction by equation ( [ eq.brownthin ] ) .
we see that the formula of equation ( [ eq.brownthin ] ) in the non - relativistic regime , namely @xmath56 , is no longer valid in the present relativistic regime , with the discrepancy increasing toward larger electron indicies .
the cross section of bremsstrahlung as a function of emitted photon energy and angle is given by @xcite and @xcite as    @xmath57\ln m(y)\right\ } , \label{km2bs}\ ] ]    where `` reduced photon angle '' @xmath58 , together with its function @xmath59 , is given by @xmath60 and @xmath61 is the initial angle of an emitted photon .
figure [ fig.bremsangdistvalid ] shows equation ( [ km2bs ] ) as a function of @xmath62 , where the results of the simulation using the thin atmosphere model are superposed .
thus , the whole emission becomes concentrated roughly within @xmath63 ; the emission is strongly forward - peaked to within @xmath64 .
the thick - target emission is significantly more complex than the thin - target case , because it results from electrons while their energies change continuously due to coulomb scattering .
moreover , it has a higher practical importance , because the flare hard x - rays and gamma - rays are thought to be produced mainly through this process in the solar chromospheric and denser regions .    under the thick - target condition , we first investigated the emission from mono - energetic electrons , as we did for the thin - target case .
we shot @xmath65 electrons with an initial kinetic energy of 10 mev or 50 mev , into a hydrogen target of which the size is ( 10@xmath66 km)@xmath25 ; this is our third simplified model atmosphere called `` thick atmosphere model '' .
figure [ fig.monothicksim ] shows the spectra produced by collecting all photons emitted by each electron while its energy gradually decreases down to @xmath67 kev .
compared with the thin - target emission from mono - energetic electrons , the photon spectra are somewhat steeper , and roll over more prominently toward the maximum energy .
obviously , this is because the electrons keep losing their energies while radiating .    in figure
[ fig.monothicksim ] , we also show the results of analytic calculations , which are formulated as @xmath68 where @xmath69 is the time , @xmath70 is the initial total energy of an electron , @xmath71 is its energy loss rate , @xmath72 is the instantaneous electron velocity , @xmath48 is the ion density of the target plasma , and @xmath73 refers to equation ( [ fig.bremsmonoenergyspectrum ] ) .
again , the simulation and analytic calculation agrees within @xmath74% .
we investigated the thick - target emission from power - law distributed electrons , with @xmath43 or @xmath75 , over the 1 - 100 mev energy range . employing @xmath76 electrons , we collected information of the photons in the same way as in the preceding simulations , and obtained figure [ fig.powthicksim ] . designating the electron injection rate ( electrons @xmath77 ) by @xmath78 and the target area by @xmath79 , and employing equation ( [ fig.bremsmonoenergyspectrum ] ) ,
the photon production rate is analytically given by @xmath80 this prediction is again show on the figure .
the simulation agrees within @xmath81% with the analytic calculation via eq.([eq.powbremsthick ] ) .
as can be seen from figure [ fig.powthicksim ] , the photon spectra emitted by power - law distributed electrons under the thick - target condition again exhibit power - law shape up to @xmath12 mev , like in the case of the thin - target emission ( figure [ fig.powerthinspectrum ] ) .
the obtained photon index is @xmath82 for @xmath83 and @xmath84 for @xmath85 .
in contrast , the hardest portion of the spectra exhibits rather opposite trend , and can not become as hard as implied by equation.([eq.brown ] ) , since the maximum energy of electrons decreases as they propagate deeper through the thick target .
figure [ fig.thinthick ] ( right ) is the same comparison as figure [ fig.thinthick ] ( left ) , but for the thick - target condition .
photon indicies in this case are generally smaller compared with thin target emission , because of the loss - flattening effect ( section [ sec.introduction ] ) . in addition
, @xmath6 depends more weakly on @xmath4 than in the thin - target case .
in the previous section , we treated the solar atmosphere as a uniform gas slab .
this is basically reasonable , because the energy losses of electrons and photons are primary determined by the column density of matter along its path . the actual solar atmosphere , with the strong vertical density gradient ,
could be transformed into such a uniform zone , as long as all the particles are traveling vertically .
however , particle trajectories with significant transverse components can no longer be conserved by such a transform , and hence we need to construct a more realistic solar model taking into account its vertical density gradients .
we here refer to `` harvard - smithsonian reference atmosphere '' @xcite for the solar mass - density profile , and approximate this numerical table analytically by an exponential function as @xmath86 here , @xmath87 is the vertical coordinate measured from the photosphere , @xmath88 is the mass density in the unit of g @xmath27 , and @xmath89 is the scale height which is @xmath90 km for @xmath91 ( the solar interior ) and @xmath92 km for @xmath93 ( the coronal region ) .
we construct our new solar atmosphere model as a series of boxes ( or `` zones '' ) stacked in the vertical direction .
individual zones are defined to have different densities but to share approximately the same column density , and the density within each of them is assumed to be uniform .
each zone is assumed to have a lateral extent of @xmath94 km , to simulate the cross section of a typical magnetic loop . on this transverse scale ,
the curvature of the solar surface is negligible .
as illustrated in figure [ fig.zones ] , we divide the coronal ( @xmath93 ) region into 5 zones , each having approximately the same column density of 0.70 g @xmath18 .    according to the harvard smithonian model ,
the total column density above @xmath95 is 3.5 g @xmath18 , which corresponds to a depth of @xmath11 km in our simple atmosphere model .
also , the sum of the 5 coronal zones in our new model gives an overall column density of 3.5 g @xmath18 , thus faithfully representing the actual solar atmosphere for @xmath93 . while traversing each zone
, a 100 mev electron loses @xmath96 mev in total , of which only @xmath81 % is in radiation whle the rest is in coulombic loss .    in simulating the solar interior ( @xmath91 ) ,
we use 10 thicker zones having 3.5 g @xmath18 each ; again , this corresponds to a depth of @xmath11 km in our simple atmosphere model . with the overall column density of 35 g @xmath18 below @xmath97 ,
our new model goes down to a depth of @xmath98 km , below the solar photosphere in comparison with equation ( [ eq.harvarddensity ] ) .
this depth is considered sufficient , because a 100 mev electron , injected at the top of our model , would lose @xmath11 % of its initial energy by the time it reaches the bottom of the deepest zone .
furthermore , the thicker 10 zones , when summed up , has a compton optical depth of @xmath12 , for a vertically precipitating photons ; few photons would penetrate deeper .
although our `` solar atmosphere '' is comprised of pure hydrogen , the actual solar atmosphere contains helium to @xmath99 by number , or @xmath100 by mass .
since the ionization loss of electrons is proportional to the electron number density , the helium would decrease the ionization loss per unit mass column density by @xmath99 .
bremsstrahlung per unit mass should remain the same , because a helium atom has 4 times higher mass than a proton but 4 times higher bremsstrahlung cross section at the same time .
the compton scattering depends only on the number of electrons , and hence would decrease by @xmath99 if helium was included .      to examine the flare photon spectra
, we vertically injected @xmath76 electrons into our solar atmospheric model from a height of @xmath101 km , and tracked them from that height to @xmath102 km , as they lose energy via coulomb scattering and emit bremsstrahlung photons .
the electrons are assumed to be distributed initially in the 1 - 100 mev range with a power - law index @xmath21 .
the electrons , initially moving in the vertical direction , gets gradually deflected as indicated in figure [ fig.powelectronspectrum ] ; in actual configuration , this corresponds to pitch - angle scattering .
nevertheless , as represented by figure [ fig.powelectronspectrum ] in red , still 98 % of the electrons are contained within @xmath103 of the vertical direction when they have reached the photosphere .
therefore , the vertical injection employed here is valid as long as the magnetic fields are sufficiently normal to the solar photosphere , the electrons pitch angle is initially concentrated near @xmath104 , and non - coulombic ( e.g. , magnetohydro - dynamic , or due to magnetic mirroring ) pitch - angle scatterings can be neglected .
the bremsstrahlung photons emitted by the almost vertically moving electrons are strongly forward collimated , as indicated by figure [ fig.bremsangdistvalid ] . therefore ,
most of them are expected to enter the atmosphere and undergo compton scattering or photoabsorption .
although some of these photons will die , others will be scattered back to escape out of the photosphere .
we collected these outcoming photons at different viewing angles @xmath14 ( figure [ fig.coordinate ] ) , and averaged the results over @xmath105 .
figure [ fig.comptongamma0deg ] shows the photon spectra obtained in this way for equal intervals in @xmath106 ( hence over an equal solid angle ) .    in figure [ fig.comptongamma0deg ] , spectra with small @xmath14 , simulating disk - center flares , exhibit a clear `` knee '' , which represents the effect of single compton scatterings . in the case of @xmath107 ( blue in figure [ fig.comptongamma0deg ] ) , the knee energy becomes @xmath108 , because of the strong energy degradation in the single compton back - scatterings .
photons with energies above the knee result either directly from the bremsstrahlung source , or through multiple compton scatterings with a relatively small energy loss in each step .
the spectrum in this region has a much steeper slope than in the region below the knee .    as @xmath14 gets larger in figure [ fig.comptongamma0deg ] , the knee energy increases , because the angle needed in the single compton scattering decreases
nevertheless , the spectrum remains much softer ( e.g. , @xmath109 in the 0.2 - 0.6 mev range ) than is expected when the compton effects are not taken into account ( @xmath110 ; figure [ fig.thinthick]b ) .
the observed photon flux below the knee energy decreases as @xmath14 increases and approach @xmath111 , due to the same limb darkening effect as in visible light .
as we have confirmed in figure [ fig.comptongamma0deg ] , bremsstrahlung photons back - scattered from the solar photosphere exhibit very steep spectra as a result of the compton degradation . in order to explain the origin of the observed flare spectra often extending beyond 1 mev , we must then considered photons which reach us after compton scatterings with much smaller angles . for this purpose , we have to consider electrons penetrating into the atmosphere with shallow angles ( i.e. , large @xmath13 ) , namely `` slant injection '' .
we then examined the case of slant injection , by shooting @xmath76 power - law electrons ( in the 1 - 100 mev energy range ) , into our realistic solar model .
specifically , we repeated the geant4 simulation under the same condition , but assuming this time an injection angle of @xmath112 ( figure [ fig.coordinate ] ) , instead of the vertical injection ( @xmath113 ) considered in the previous subsection .
we now inject the electrons from a height of 100 km ( i.e. , neglecting the uppermost two zones in figure [ fig.zones ] ) , to avoid them escaping sideways before arriving at thicker regions .
figure [ fig.comptongamma80deg ] shows the results of this simulation for the same three intervals of @xmath14 as figure [ fig.comptongamma0deg ] , again averaged over the observer s azimuth angle @xmath15 .
the spectrum becomes considerably harder than in the case of vertical electron injection , and extends much beyond @xmath114 mev when @xmath14 is rather large ( @xmath115 ) . in this case , electrons make grazing incidence on the atmosphere , and some of their bremsstrahlung photons get nearly forward scattered to leave the zone with a small energy loss , thus producing the hard power - law continuum .    in this way , we have monte - carlo simulated slant injections for a family of electron spectra with various indicies @xmath4 ( again distributed over 1100 mev ) , assuming two representative injection angles of @xmath116 and @xmath117 .
figure [ fig.deltagamma ] summarizes the photon index @xmath6 of 110 mev gamma - rays , collected over the most favorable viewing angle of @xmath118 ( i.e. , rather grazing to the solar photosphere ) and integrated over the observer s azimuth .
thus , even with the extreme assumption of @xmath119 , the simulated gamma - ray spectra become significantly softer ( typically by @xmath114 in @xmath6 ) than the original bremsstrahlung spectra before the photons are compton scattered ( reproduced in figure [ fig.deltagamma ] in blue ) .
moreover , @xmath6 further increases by @xmath120 as @xmath13 decreases from @xmath121 to @xmath122 .
these results imply that the compton scattering significantly soften the emergent photon spectra , of which the effect depends sensitively on @xmath13 and @xmath14 .
an interesting prediction of figure [ fig.deltagamma ] is that the observed gamma - ray photon index should appear in a relatively narrow range of @xmath123 , regardless of the electron index @xmath4 .
of course , we would observe considerably larger values of @xmath6 when either @xmath13 or @xmath14 is rather small .
however , such a spectral component would not be easily detected in the mev energy region .
although our simulations of the slant injection neglected the uppermost two coronal zones , their effects may not be necessarily negligible , because the large injection angle @xmath124 ( or @xmath122 ) effectively increases the electron path length in each zone by @xmath125 [ or @xmath126 times .
since the electron energy losses in these two regions , if properly taken into account , would reduce the maximum electron energies while flatten the low - energy electron slope ( figure 3 ) , the emergent photon flux above @xmath20 1 mev would further be suppressed if the electron injection in actual flares takes place at a height much exceeding the assumed 100 km . in other word ,
our simplification employed in simulating the slant injection makes the case more conservative .
in order to investigate the effects of compton scatterings on the solar gamma - ray spectra , we have conducted monte carlo simulations of energetic electrons , impinging on the solar atmosphere , using geant4 as our basic toolkit .
we modeled the solar atmosphere as vertically stacked parallel hydrogen zones , each assumed to be uniform .
the density gradient was expressed by assigning higher densities to lower zones .
we injected electrons which are power - law distributed over the 1 - 100 mev range , changing their spectral index @xmath4 and incident angle @xmath13 .
the results of the present simulations can be summarized as follows .    1 .
the observed photons hardly reach mev energies if the electrons are injected vertically , because the bremsstrahlung - produced gamma - rays must be compton back - scattered to reach the observer , and hence their energies are subject to strong compton degradation .
the emergent gamma - ray spectrum , integrated over the observer s azimuth , gradually hardens as the injection angle @xmath13 increases ( becoming more `` slant '' ) , and also as the observer s polar angle @xmath14 gets larger .
3 .   as long as @xmath13 and @xmath14 are both relatively large
, the azimuthally - averaged 1 - 10 mev gamma - ray spectra exhibit @xmath127 over a relatively wide range of @xmath4 .
these results generally agree with , and further extend , the previous works by , e.g. , @xcite and @xcite .
there are some effects which are not considered in our simulation . among them , those of the secondary electrons and helium were already described in section 2.2 and section 3.1 , respectively . as the electron energy increases
, their synchrotron energy loss becomes significant .
however , in the chromospheric region in which we are interested , the magnetic field intensity is typically @xmath11 g , and hence the rate of synchrotron loss of 100 mev electrons is only @xmath114 % of their coulomb loss , and only @xmath128 % of their bremsstrahlung loss .
this is because the coulomb and bremsstrahlung losses are rather large due to the high matter density .
thus , the synchrotron process can be neglected .
we also ignored the inverse compton process between flare - accelerated electrons and ambient photons , because electrons with the maximum energy assumed , 100 mev , can produce at most hard x - rays up to @xmath11 kev by scattering off the solar visible photons .
however , if the electrons have a very flat spectrum well extending to energies above 350 mev , the inverse compton process may become a key process in creating gamma - ray photons .
the size of the world used in our simulation , 10000 km @xmath129 10000 km @xmath129 20000 km , simulates a typical size of a foot point of a magnetic tube .
the number of photons which escaped from the boundaries below the photosphere is less than @xmath130 of that of produced photons , and energies of these escaped photons are less than 100 kev each .
hence we ignored those photons in the simulation .
possible effects of the relatively low heights of our slant injection was already considered in subsection 3.3 .    in the present simulation , we truncated the power - law electron spectrum at 100 mev . in order to examine whether this abrupt cutoff produces any artifact , we repeated some of the simulations by changing the maximum cutoff energy to 1 gev .
q then we have found that @xmath6 decreases only slightly , typically 0.1 or less , if @xmath4 is larger than 2 .
finally , a separate gamma - ray component could be produced without the compoton effect , simply if the electron spectrum consists of two power - law components with different slopes .
even so , the good agreement between the predicted range of power - law indices and the _ yohkoh _ data suggests that the compton scattering is playing a significant role ( see the next section ) .      in order to evaluate the significance of the results of our simulation ,
let us compare them with actual observations of hard x - ray and gamma - ray emission from solar flares .
for that purpose , we have also to relate the simulation parameters @xmath13 and @xmath14 to observational parameters , including pitch - angle distributions of electrons , the tilt of magnetic field lines , and the flare longitude on the solar disk . as the observational data , we utilize 40 flares that show significant gamma - ray emissions , selected out of 2788 x - ray flares observed from 1991 october to 2001 december by yohkoh @xcite .
yohkoh , the japanese solar satellite launched in 1991 , was able to take four - color hard x - ray flare images with the hard x - ray telescope(hxt ) @xcite and perform spectro - photometry over the @xmath131 mev range with the gamma - ray spectrometer ( grs ) @xcite .
as the simplest case , we may consider a condition where the magnetic field lines are perpendicular to the photosphere . in this case , we can identify @xmath13 with the pitch - angle of an electron , and the effect of electron gyration around the magnetic field lines is fully represented by our photon collection method which takes an average over @xmath15 . if the electrons have a narrow pitch - angle distribution around 0 , the strong compton degradation , described in section 3.2 and shown in figure [ fig.comptongamma0deg ] , will make the observed photon spectrum very steep , and will prevent it from reaching mev energies . if , on the other hand , the pitch - angle distribution is rather broad , electrons with larger values of @xmath13 will emit harder spectra , as represented by figure [ fig.comptongamma80deg ] .
conversely , as long as the field lines are close to vertical , flares with significant gamma - ray emission must involve broad pitch - angle distributions of electrons .
analyzing the 40 gamma - ray flares detected by yohkoh , @xcite found that their gamma - ray ( typically 1 mev ) to hard x - ray ( @xmath20 70 kev ) flux ratios scatter largely ( by an order of magnitude ) , with no correlations to their hard x - ray spectral slopes .
hence they conclude that there is a hidden parameter which causes the gamma - ray flux to vary significantly relative to the hard x - ray extrapolation , in such a way that the gamma - ray component often appears as a separate spectral hard tail ( see their figure [ fig.thinthick ] ) .
given the present results , this parameter is very likely to be the electron pitch - angle distribution , with broader distributions yielding higher gamma - ray to hard x - ray flux ratios .
in addition to the electron injection angle @xmath13 , the viewing angle @xmath14 has been confirmed to strongly affect the flare gamma - ray spectra ( section [ sec.solarsimulation ] ) .
assuming again the vertical field geometry , we can identify @xmath14 with the angle @xmath132 which our line - of - sight to the flare makes relative to the normal to the photosphere .
hereafter , we simply call @xmath132 `` flare - viewing angle '' . using the flare longitude @xmath133 and latitude @xmath134 , it can be calculated as @xmath135    to examine the effect of @xmath14 , we plotted in figure [ fig.longitude ] the gamma - ray to hard x - ray hardness ratio of the 40 yohkoh flares , as a function of @xmath132 calculated via equation ( [ eq.theta_v ] ) .
it is similar to figure 7 of @xcite , but differs from it in that @xmath132 is used instead of @xmath136 .
as predicted by our simulation , the hardness ratio clearly increases with @xmath132 ; this effect was already known as `` gamma - ray limb brightening '' in solar flares using _
data @xcite .
@xcite interpreted this effect in terms of the radiation anisotropy : supposing that the electrons have a pitch - angle distribution which increases with angle from the outward normal , then higher - energy continua are expected to be more strongly beamed parallel to the photosphere , because harder bremsstrahlung photons are more strongly forward - beamed than softer ones even when emitted by electrons of the same energy .
we presume that this mechanism can explain at least part of the observed limb - hardening effects , while the compton back - scattering process must be enhancing them .
in particular , our scenario can explain limb hardening even when the electrons are mostly directed downwards .    for a more quantitative comparison of our simulation with the yohkoh observations ,
let us limit ourselves to limb flares which have absolute longitudes larger than @xmath137 , because gamma - ray data of disk flares are of poor photon statistics than those of limb flares as is clear from figure [ fig.longitude ] , and because these disk flares are more contaminated by nuclear gamma - ray lines @xcite which are unrelated to the accelerated electrons .
figure [ fig.color_color2 ] shows a scatter plot between two hardness ratios of the 40 sample flares , one calculated in the hard x - ray range ( 5793 kev vs 3357 kev ) while the other in the gamma - ray range ( 1.436.21 mev vs 0.221.43 mev ) . in this `` color - color '' plot
, the dashed line represents the condition that the two colors imply the same photon index @xmath6 , of which values are given in the figure . the hard x - ray ( or gamma - ray ) photon index of a flare
can be found by drawing a vertical ( or horizontal ) line from the data point , and read its intersection with the dashed line .    from figure [ fig.color_color2 ]
, we can see that the gamma - ray spectral index of the 40 flares is distributed in a typical range of 1.7 - 2.5 .
this agrees very well with the prediction by our monte - carlo simulations ( figure 13 ; subsection 3.3 ) , on condition that the electrons in these flares have broad pitch - angle distributions beyond @xmath138 .
further , in figure [ fig.color_color2 ] , a majority of flares are distributed to the left side of the single power - law locus .
such a concave - shaped wide - band spectrum can be naturally explained as a consequence of a broad pitch - angle distribution of electrons , because electrons with small pitch angles are expected to emit steep - spectrum photons , while those with large pitch angles emit much harder spectra .
a superposition of these spectral components with different slopes will result naturally in a concave spectrum just as is observed .
thus , as long as the magnetic lines are assumed perpendicular to the photosphere and the flare viewing angle @xmath132 is large , flares with significant gamma - ray signals are inferred to have rather broad pitch - angle distributions .      as a more general case ,
we may briefly consider tilted magnetic - field conditions . in this case , it is somewhat difficult to directly compare the results of our simulation with the observational data , because the electron gyration can no longer be faithfully emulated by our method of taking an average over @xmath15 .
nevertheless , we may employ the present results to predict some aspects of that case in a qualitative way .
let us first consider limb flares as an extreme case .
in this case , a field - line tilt on the sky plane is considered to have little effects , because the most efficient gamma - ray emission ( as seen from the observer ) occurs over the gyration phase wherein the electrons are approaching us , and this phase is the same even the filed lines are tilted in that way . in contrast , a field - line tilt along our line - of - site is expected to have a prominent effect : if the field lines are tilted away from us , the observed spectra tend to reach higher energies , e.g. a few mev , even if the electrons have a rather narrow pitch - angle distribution .
if the lines are tilted toward us , the gamma - ray spectra are expected to soften , on the contrary .
then , the observed scatter in the gamma - ray to hard x - ray intensity ratio for limb flares ( figure [ fig.longitude ] , with @xmath139 ) may be contributed by the field - line tilts , in addition to the pitch - angle distribution .
the other extreme case is a disk - centered flare , where the field - line tilt is specified only by its angle @xmath140 with respect to the photosphere normal ( i.e. , our line - of - sight ) without depending on the tilt direction .
then , an electron with a fixed pitch angle , @xmath141 is expected to sweep a range of injection angle , from @xmath142 to @xmath143 .
therefore , a larger value of @xmath140 is favorable in delivering strong gamma - ray signals to the observer , even though @xmath14 is still limited to the unfavorable value of @xmath144 . if @xmath141 and @xmath140 are both rather large so that @xmath145 , the electron is implied to be approaching us over a certain phase of its gyration , thus allowing a considerable fraction of the bremsstrahlung photons to directly reach us without being compton scattered .
therefore , a disk flare may become a strong gamma - ray source under a limited condition , that the electrons have a broad pitch angle and the magnetic field lines are significantly tilted .
now that broad pitch - angle distributions of accelerated electrons are considered to be an important element of observing intense gamma - ray emission , it may be imperative to simulate such a case , and compare the results with gamma - ray spectra from actual solar flares .
therefore , we repeated our solar simulation , assuming that the electrons have a continuous incident - angle distribution as @xmath146 which corresponds to @xmath147 in @xcite except for covering only the @xmath148 range . when the magnetic fields are vertical , equation ( [ eq.dermerramaty1986 ] ) becomes identical to the pitch - angle distribution of the accelerated electrons
. under this angular distribution
, we injected @xmath65 electrons in the same way as in section  [ sec.solarsimulation ] , assuming the incident electron spectrum to have a canonical index of @xmath149 , and to extend again over 1100 mev .
collecting the produced gamma - rays over a viewing angle range of @xmath150 , and convolving the results with an approximate response of the yohkoh grs @xcite , we have obtained the spectrum shown in figure  [ fig.grs_sim_comp ] in red .    to compare with this simulation ,
we have chosen a goes x3.7 class flare of 1998
november 22 , 06h37 m ut , which is one of the gamma - ray brightest events among the 40 yohkoh flares @xcite .
since this is an extreme limb flare located at s32w90 , our viewing angle @xmath14 should be close to @xmath151 , in agreement with the simulation condition .
actually , the hxt highest - band images of this flare reveal intense hard x - ray emission from a pair of magnetic loop foot points , both located right on the solar west limb @xcite .
the grs data , accumulated over a period of 06:37:2206:40:22 , are shown in figure [ fig.grs_sim_comp ] in blue , after subtracting background ( accumulated for 06:56:1406:59:02 ) , but without removing the detector response .
thus , the spectrum is featureless , and the 2.23 mev neutron capture line is undetectable with its equivalent width being @xmath152 kev @xcite .
therefore , this flare is considered to be a typical `` electron dominated '' event .    in figure
[ fig.grs_sim_comp ] , the simulation and the actual data agree reasonably well , at least over the 110 mev band .
when the grs response is considered , these spectra can be approximated both by @xmath153 @xcite , which is typical in figure  [ fig.color_color2 ] .
incidentally , we would obtain from figure  [ fig.deltagamma ] @xmath154 for the assumed @xmath149 , if we were observing direct thick - target emission .
we should instead find @xmath155 ( for @xmath117 ) to @xmath156 ( for @xmath116 ) , if the compton back - scattering dominates .
therefore , the values of @xmath157 found in figure  [ fig.grs_sim_comp ] are in between these two cases .
indeed , an inspection of the monte - carlo photons revealed that @xmath158 of them are compton - scattered gamma - rays , while the remaining @xmath159 are direct bremsstrahlung events .
given the good agreement between the two data sets , it is possible that the actual gamma - rays from the november 22 flare are also a mixture of the direct and scattered photons .
in order to quantitatively estimate the spectra of energetic electrons generated in solar flares based on the observed gamma - ray spectra , we numerically studied the electron transport and gamma - ray emission in the solar atmosphere .
as elementary processes , we mainly considered coulomb scattering , bremsstrahlung and compton scattering .
we modeled the solar atmosphere with a vertical stack of parallel zones with different densities , and neglected the magnetic fields .    as a result of our simulation
, we have found that the thick - target gamma - ray spectra emitted from magnetic loop foot - point can hardly reach mev energies , because of heavy compton degradation , unless the electrons make rather grazing angles to the photosphere and we are observing forward - scattered gamma - rays . when this compton effect is taken into account , the traditional @xmath4 vs. @xmath4 relation
is drastically modified .
the pitch - angle distribution of flare - accelerated electrons are suggested to be broad and extend to large angles for those flares of which the emitted photon spectra reach mev energies .
actually , a gamma - ray spectrum simulated under electron injection with a broad pitch angle distribution reproduced reasonably well the yohkoh spectrum of the 1998 november 22 flare .      1 .
when the magnetic field lines are perpendicular to the photosphere , and the electrons have a narrow pitch - angle distribution around 0 , the gamma - ray spectra are very soft and can not reach mev energies .
2 .   even in the same field configuration
, the gamma - ray spectra become harder if the electrons have larger pitch angles .
therefore , the electrons producing gamma - rays are inferred to have relatively large pitch angles .
limb flares are expected to have harder spectra than disk flares , because of their smaller angles of compton scattering .
this agrees with the observed fact that the gamma - ray to hard x - ray intensity ratio increases as the flare longitude gets closer to the solar limb .
4 .   in case of a broad pitch angle distribution , electrons with small pitch angles
emit steep hard x - rays , while those with large angles emit photons with flatter spectra .
this can explain another yohkoh result , that the gamma - ray slope tends to be flatter than the contemporaneous hard x - ray slope .
5 .   under the conditions of @xmath160 and @xmath161 ,
the values of @xmath6 obtained in our simulations in the 110 mev range are concentrated over 1.72.5 .
this agrees again with the yohkoh observations .
a disk flare may become a strong gamma - ray source under a limited condition , that the electrons have a broad pitch - angle and the magnetic field lines are significantly tilted .
the gamma - ray spectra of actual solar flares are expected to be an appropriate mixture of the direct and compton - scattered photons .      ,
s. , allison , j. , amako , k. , apostolakis , j. , araujo , h. , arce , p. , asai , m. , axen , d. , banerjee , s. , & barrand , g. 2003 , nuclear instruments and methods in physics research section a : accelerators , spectrometers , detectors and associated equipment , 506 , 250                              , r.  p. , dennis , b. , hurford , g. , smith , d.  m. , & zehnder , a. 2004 , in telescopes and instrumentation for solar astrophysics .
edited by fineschi , silvano ; gummin , mark a. proceedings of the spie , volume 5171 , pp .
38 - 52 ( 2004 ) . , 3852


there is a strange omission in the traditional version of the second law of thermodynamics ( slt ) .
it does not take gravity into account .
however the existence of the arrow of time is usually explained with reference to slt . and most discussions of the origin of the arrow of time appeal ultimately to the initial condition and the evolution of the universe .
gravity plays an essential role for this evolution .    as pointed out by p.c.w .
davies @xcite there seems to be a paradox that the material contents of the universe began in a condition of thermodynamic equilibrium , whereas the universe today is far from equilibrium . hence the thermodynamic entropy has been reduced in conflict with slt which says that the entropy of the universe is an increasing function of cosmic time .
this is due to the tendency of self - gravitating systems irreversibly to grow inhomogeneous .    in order to include this gravitational effect into a generalized version of slt , one has to define a gravitational entropy .
several tentative definitions have been given based on penrose s weyl curvature hypothesis @xcite .
the weyl curvature tensor vanishes identically in the homogeneous and isotropic friedmann - robertson - walker universe models .
it was suggested by penrose @xcite to use this tensor as a measure of inhomogeneities of the universe models .
the weyl curvature scalar is expressed by scalars constructed from the riemann curvature tensor as follows c^c_=r^r_-2r^r_+13
r^2    penrose has given several formulations of the conjecture .
a version that emphasizes the entropy aspect was given in ref .
@xcite p.178 .
it seems that in some way the weyl tensor gives a measure of the entropy in the space- time geometry .
the initial curvature singularity would then be one with large ricci tensor and vanishing weyl tensor ( zero entropy in the geometry ) ; the final curvature singularity would have weyl tensor much larger than ricci tensor ( large entropy in the geometry ) .    in order to quantify penrose s conjecture wainwright and collaborators
@xcite have suggested that the quantity p^2= [ defp]may represent a `` gravitational entropy '' , at least in a cosmological context with a non - vanishing ricci - tensor .
we could of course ask ourself whether the entropic behavior of the weyl curvature invariant ( or an invariant composed thereof ) results because the weyl tensor is _ directly related _ to the gravitational entropy , or whether it is only a side - effect because the universe evolves towards a state of maximal entropy . hence ,
even though the results may show that the weyl curvature invariant has an entropic behaviour , it is by no means a proof that the weyl tensor should be identified as the `` gravitational entropy '' .
in an earlier article @xcite we discussed the weyl curvature tensor in a homogeneous but anisotropic model as well as in an inhomogeneous model . our survey will be summarized in this and the next sections , as well as bringing new arguments into the discussion .    as a basis for our study , we used two different models which induce weyl curvature effects in two conceptually different ways .
one was the anisotropic bianchi type i model , the other was the inhomogeneous lematre - tolman model .
the weyl curvature conjecture has been investigated in the szekeres cosmological model , that generalize the lematre - tolman model , by w.b .
bonnor @xcite .
inhomogeneous modes are in general local modes for general relativity while anisotropic modes are global modes .
the global topology and geometry of our universe has significant consequences for the possibilities for anisotropic modes @xcite .
the role of inhomogeneous modes can be considered as more local in origin , but the effect on the global geometry and topology are not known in detail .
hence , it is important to study both an inhomogeneous model as well as a homogeneous model to get a more complete picture of the behaviour of the weyl tensor in generic cosmological models .
the bianchi type i model is the simplest anisotropic generalisation of the frw models .
it has flat spatial sections , and hence , can be a good candidate for the universe we live in .
the metric for this model can be written as ds^2=-dt^2+e^2_ijdx^idx^j where @xmath0 .
in @xcite we solved the einstein field equations for dust and a cosmological constant @xmath1 . by introducing the volume element @xmath2 , the einstein field equations gives the following equation for @xmath3 : ^2=3v^2 + 3mv+a^2 where @xmath4 is the total mass of the dust , and @xmath5 is an anisotropy parameter .
the equations for @xmath6 are in terms of the volume element _ = where the constants @xmath7 are related to @xmath5 via @xmath8 .
it is useful to define an angular variable @xmath9 by @xmath10 and @xmath11 .
an important special case of the bianchi type i is the kasner vacuum solutions .
the kasner solutions are characterised only by the angular variable @xmath9 .
these solutions have a weyl scalar ( c^c_)_i= ( 1 - 2z3+z^2 ) [ weylbti]where @xmath12 and @xmath13 for the kasner solutions which have vanishing ricci tensor .
so in a sense , the kasner solutions are the anisotropic counterpart to the inhomogeneous schwarzschild solution
.    the weyl tensor will decay as the volume expands , even at late times .
one should expect that this entity would increase if it represents gravitational entropy .
but is decreases monotonically , and hence , it is doubtful that it is the correct measure in these models .
inserting dust and a cosmological constant will not have a significant effect on this decreasing behaviour . if @xmath4 is the mass of the dust inside the volume @xmath3 , then eq .
( [ weylbti ] ) still holds but with z=. hence , except for the special case is somewhat interesting . for
@xmath14 it can be shown that the kasner solution is just a special part of minkowski spacetime , which has no singularities .
it seems a bit odd that for any @xmath15 we will have an initial singularity , while for @xmath14 we do not have one .
however , if we compactify the spatial sections this oddity disappears .
the model has a singularity at @xmath16 even for the case @xmath14 ! ]
@xmath14 , the weyl tensor diverges as @xmath17 as @xmath18 .    for the bianchi type i universe models the entity @xmath19 defined in eq .
( [ defp ] ) , turns out to be ( p^2)_i = .
also this measure for the weyl entropy diverges as @xmath20 .
let us now consider the inhomogeneous lematre - tolman ( lt ) models .
the line element for the lt models can be written as [ ltmetric ] ds^2=-dt^2+q^2dr^2+r^2(d^2+^2d^2 ) where @xmath21 and @xmath22 .
in @xcite we used this model and investigated the solutions of the einstein field equations where the spacetime contains a cosmological constant @xmath1 and dust .
the equations then turn into @xmath23 where @xmath24 is an arbitrary function , and [ meq ] r^2+(1 -f^2)r - r^3=m .
[ masseq]the function @xmath25 is given by the integral m(r)=^r_0 4r^2 rdr , where prime denotes derivative with respect to @xmath26 and @xmath27 is the dust density .
hence , @xmath28 can be interpreted as the total mass of the dust inside the spherical shell of coordinate radius @xmath26 . inverting the equation for @xmath28
, the dust density can be written in terms of @xmath29 and @xmath28 : 4= [ rhoeq ]    it is also useful to define the _
mean dust density function _ @xmath30 by the relation m(r)= |r^3 .
[ rhobareq]interestingly , the weyl curvature scalar can now be written quite elegantly as ( c^c_)_lt= ^2(|-)^2 .
[ rhoweyl]this relation provides us with a physical interpretation of the weyl scalar in the lt models .
it is just the difference between the mean dust density and the actual dust density .
also , the weyl tensor is everywhere zero , if and only if @xmath31 . in this case
the lt models turn into the frw models with homogeneous spatial sections .
the frw models are conformally flat , hence they have a zero weyl tensor .
let us return to the equations of motion for the lt model .
the mass equation ( [ masseq ] ) can be written as an `` energy '' equation : ^2+v(r , r)=e(r ) [ energyeq]where the `` potential '' @xmath32 and the `` energy '' @xmath33 are given by v&= & --r^2 + e&=&-(1-f^2 )    this energy equation may be integrated and solved exactly .
a summary of the results and a qualitative description of the physical meaning of the solutions is given in  @xcite .
the solutions may be written in terms of the weierstrass elliptic functions  @xcite .
the actual expressions are not very informative unless the reader has massive knowledge of these elliptic functions . however a lot of qualitative information can be extracted from simple classical considerations .
the classical solutions will move on level curves of the energy function @xmath34 , since the total energy @xmath33 is independent of @xmath35 . in figure
[ levelcurves ] the level curves of a typical energy function are drawn . if @xmath36 there will exist a saddle - point of the energy function .
this saddle point will be at @xmath37 where the energy function will have the value @xmath38 .
the saddle point solution is static , and is as a matter of fact the einstein static universe .
this can be done since we have to assume that @xmath39 on physical grounds . in the case
@xmath40 , the metric become degenerate . ] .
if @xmath41 , the solutions fall into two distinct classes :    1
.   schwarzschild - like solutions : these solutions expands , but they do not possess enough energy to escape the gravitational collapse , so they end as black holes . if @xmath42 is constant these solutions are those of a schwarzschild black hole in lematre coordinates .
de sitter - like solutions : solutions where the universe evolves approximately as that of de sitter solutions with positively curved hypersurfaces .    if @xmath43 the ( test ) matter has enough energy to escape the gravitational collapse ( expanding solutions ) or enough energy to prevent the gravitational repulsion from the cosmological constant ( contracting solutions ) .
near the initial singularity , both @xmath44 and @xmath27 diverge . unless they are identically equal
, the weyl tensor will diverge near the initial singularity .
explicitly we have near the initial singularity 4(|-)= where we have used that near the initial singularity , we can approximate the solutions with @xmath45 the free function @xmath46 is the big bang time . to avoid intersecting world - lines we have to assume @xmath47@xcite , thus @xmath48
. the weyl scalar will diverge as @xmath49 unless @xmath50 .
also the entity @xmath51 will diverge in general near the initial singularity .
@xmath51 is given by @xmath52 which near the initial singularity can be approximated by ( p^2)_lt()^2 .
hence , @xmath51 diverges as @xmath53 unless @xmath50 , i.e. unless the big bang is homogeneous .
this does not prove to be a very promising behaviour for the wcc .
the curvature scalars diverge near the initial singularity , and hence , come in conflict with the wcc . but let us analyse the situation in the lt a bit more carefully . as the universe expands , both @xmath44 and @xmath27 will decrease .
they both decrease from an infinite value at the initial singularity . as the universe expands the value of @xmath44 decrease as @xmath54 while @xmath27 decreases as @xmath55 close to the initial singularity .
hence , the value of @xmath44 approaches the value of @xmath27 by a factor or @xmath55 , but since they both diverge , the weyl tensor and @xmath51 diverge as well .
however , even though r. penrose and s. hawking @xcite showed that according to the general theory of relativity the big bang must have started in a singularity , it has later been emphasized that the physical universe must obey not only the general relativistic laws of nature , but also the quantum mechanical laws .
hence , the initial singularity is a fiction which not corresponds to physical reality .
the classical laws can be applied only after the planck time .
we should therefore study the behaviour of @xmath19 not only in the limit @xmath56 , but rather at a very small cosmic time .
choosing the origin of cosmic time @xmath57 at @xmath58 we introduce @xmath59 . hence close to @xmath60
@xmath61_{t\lra 0}<0 $ ] .
bonnor @xcite considered the model with @xmath62 and a vanishing cosmological constant , and found the opposite result .
however he restricted his investigation to models with homogeneous initial singularity , i.e. @xmath63 . as we will show this is an exceptional case .
the behaviour of more general models with @xmath64 is different . to see what is actually happening
if we push @xmath65 towards zero , we utilize that close to @xmath60 the models with @xmath62 and @xmath66 behave according to r^3=92 mt^2(1-(r)t^23 ) [ smallr3]where @xmath67 is assumed to be a small function of @xmath26 . from eqs .
( [ rhoeq ] ) , ( [ rhobareq ] ) and ( [ rhoweyl ] ) with @xmath66 we get |p|=(3 - 1 ) . using ( [ smallr3 ] )
we obtain |p|=. this function is plotted in fig .
[ pperturb ] .
note that this function diverges for @xmath68 , obtains a minimum for some small @xmath57 and increases thereafter .
as @xmath69 is driven towards zero this minimum goes towards @xmath60 .
let us consider a finite 3-volume @xmath32 comoving with the cosmic gas in our space time . according to the first law of thermodynamics the matter entropy @xmath70 and the internal energy @xmath71 will evolve as : tds_m = du+pdv .
if the matter content is dust then @xmath72 .
there are a couple of things to note .
firstly , in a dense dust cloud , we expect the internal energy of the dust to be large , thus we expect the entropy to be large . secondly , the entropy is an increasing function of the volume .
let us therefore consider a co - moving volume @xmath32 in our spacetime . what could the expression for a _ gravitational entropy _ be ? from the previous sections we noticed that an lt model with @xmath31 is homogeneous .
the quantity @xmath73 is an inhomogeneity measure in the lt models . in the absence of
a cosmological constant we notice that p= ( ) .
[ eq : pdiff]the sign of @xmath51 is here chosen so that the configuration @xmath74 is associated with @xmath75 while the more realistic configuration @xmath76 has @xmath77 corresponding to eq .
( [ eq : pdiff ] ) positive .
let us consider the entity defined by : = _ v p dv .
introducing co - moving coordinates @xmath78 we write @xmath79 where @xmath80 is the 3-volume element .
the integration range is now constant as a function of time and if we integrate over a unit coordinate volume which is so small that the integrand is approximately constant , we may write = _ v p dvp [ defs ]    this formulation of the gravitational entropy is an intuitive one , but it is not a covariant formulation . to find a covariant formulation
we define the entropy current vector by @xcite = s*u*+ where @xmath81 and @xmath82 are orthogonal , @xmath83 is the entropy density , @xmath81 is the material flow vector ( @xmath84 ) and @xmath82 is the entropy flux . the second law of thermodynamics can now be expressed as _ ; 0 . writing this in a local coordinate system ,
the divergence is _ ; = _ ( ^ ) . in our case
we take the entropy density proportional to @xmath51 and assume a vanishing entropy flux .
hence , p*u*. in comoving coordinates @xmath85 so that _ ; _ t(p)=. for the `` second law to hold '' , we have to check whether the entity @xmath86 is increasing .
if we take the corresponding one - form of the entropy vector , and take the dual of this form ( hodge dual ) , we obtain the entropy three - form which is given by the contraction of the space - time volume form @xmath87 with the entropy current current vector = i_. in our case , this three - form is simply = p*dx * hence , if the coordinates are comoving , then it is the component of this three - form that is increasing as long at the matter and fields obeys the sec .
this explains why the entropy should scale as the volume .
it is usually this entropy three - form we are thinking of .
it has a more intuitive behaviour than the entropy current vector .
we will therefore use @xmath88 in the further to study the wcc .
the square of the weyl tensor and the entity @xmath51 do not , as we have seen , capture the entropic behaviour properly .
they both diverge at the initial singularity and thus , can not be a proper measure of the gravitational entropy .
the entity @xmath88 on the other hand , seems to be a more promising candidate for the weyl entropy .
let us start with the lt model .
the motivation for studying the entity @xmath88 given in ( [ defs ] ) came from the lt models .
therefore it is natural to start with this model . in all the lt models , we can approximate the behaviour near the initial singularity , @xmath49 with s_lt=2(m(t - t_0)-2t_0 ) .
[ eqsmallt ] for physically realistic spacetimes , @xmath89 , so as @xmath90 , @xmath88 is positive and finite .
note also that @xmath88 is increasing , as suggested by the wcc . in the absence of a cosmological constant and
if @xmath91 , this will be the exact expression for @xmath88 . if a cosmological constant is present , the universe will eventually go into a de sitter phase if the universe is allowed to expand for ever .    for the sake of illustration
, it is useful to investigate a specific case
. let us choose @xmath92 .
the energy equation eq .
( [ energyeq ] ) , can now be solved to yield r=()^13^ 23[eqexactlt]where @xmath93 . in general , the late time behaviour of the lt model ( if the universe has one ) is that of a de sitter universe .
this exact solution provides us with a solution which connects the initial singularity with the late time de sitter era .
the initial behaviour of @xmath94 is the same as in eq .
( [ smallimit ] ) .
the late time behaviour of @xmath29 is r=()^13e^h(t - t_0(r ) ) .
hence , at late times the entity @xmath88 is approximately constant : s=. [ eq : latetimedesitter]at late times , @xmath88 approaches a constant value which is inversely proportional to @xmath95 .
the larger the value of the cosmological constant is , the smaller the value of the final value of @xmath88 . in fig .
[ fig : sevol ] we have plotted , using the exact solution ( [ eqexactlt ] ) , the evolution of @xmath88 for three different values of @xmath1 .
the volume under consideration expands at late times exponentially as @xmath96 , where @xmath95 is the hubble parameter .
hence , the entropy per unit volume will decrease exponentially : e^-3ht .
the entropy in a unit volume will therefore decrease rapidly during a de sitter stage .
this is a very important consequence of inflation . at the exit of the inflationary period ,
the entropy in a unit volume is very low , there are very little inhomogeneities left of the primordial ones .
let us summarize the evolution of the entity @xmath88 in the lt models :    1 .   * large @xmath1 and ever expanding : * in the initial epoch the dust dominates and @xmath88 is increasing linearly in @xmath35 .
the universe is becoming more and more inhomogeneous .
after the universe has grown considerably , the cosmological constant becomes dominant , @xmath88 stops growing and if the cosmological constant is large enough , it evolves asymptotically towards a constant value .
the universe is smoothened out .
this scenario corresponds approximately to the lowest two graphs in fig .
[ fig : sevol ] .
* small @xmath1 and ever expanding : * again the dust dominates initially .
the cosmological constant is too small to make @xmath88 decreasing .
the entity @xmath88 is ever increasing but is bounded from above by a relatively large constant value .
* zero @xmath1 and ever expanding : * the @xmath88 will again be ever increasing and will asymptotically approach a function @xmath97 where @xmath98 and @xmath99 are constants and @xmath100 iff @xmath101 and @xmath102 iff @xmath92 .
* recollapsing universe : * due to the dust term , the final singularity will not be similar to the initial singularity .
hence this entity is asymmetric in time for a recollapsing universe .    in the lt models we see that @xmath88 behaves in agreement with the wcc .
let us investigate more carefully the schwarzschild spacetime , where the whole motivation of gravitational entropy comes from @xcite .
the schwarzschild spacetime is a special case of the general lt models ; it has @xmath103 and a vanishing ricci tensor .
if we look at the entity @xmath88 in the region outside the schwarzschild singularity , @xmath88 will diverge .
this is in some sense the maximal possible value of @xmath88 , the weyl tensor is as large as possible and the ricci tensor is the smallest as possible .
thus at this classical level it seems that the schwarzschild spacetime has the largest possible @xmath88 which is a good thing if one wants to connect @xmath88 with the entropy of the gravitational field . outside a real black hole
, @xmath104 is probably impossible . even though the classical vacuum has vanishing ricci tensor , the quantum vacuum will probably not have @xmath104
the quantum fields will fluctuate and cause the expectation value of the square of the ricci tensor to be non - zero : @xmath105 .
hence , there will probably be an upper bound of how large @xmath88 can be , even in a vacuum .      in the bianchi type i with dust ,
the explicit expression for the entity @xmath88 is s_i = p=()^ where @xmath106 . in a neighbourhood of the singularity we can make a taylor expansion to first order in @xmath3 : @xmath107 we interpret @xmath108 as the dust density in coordinate space @xcite . hence
, as we go towards the initial singularity @xmath88 will be finite .
for the vacuum case , @xmath109 and @xmath66 , and the solutions are the kasner solutions .
the kasner solutions have a vanishing ricci tensor , but a non - zero weyl tensor . hence ,
as in the schwarzschild case , the entity @xmath88 will diverge .    from the taylor expansion we see that , if @xmath110 then @xmath88 will increase immediately after the initial singularity . if @xmath111 then @xmath88 will decrease .
this can be understood as follows .
if the cosmological constant is to large , the universe will increase too rapidly initially to allow the dust to contribute significantly to the anisotropy .
initially the universe will become more and more isotropic . in the presence of a cosmological constant
, the universe will eventually enter a de sitter phase . again
@xmath88 will evolve towards a constant value given by @xmath112    if the cosmological constant vanishes
@xmath88 is monotonically increasing : @xmath113 note that for large @xmath35 we have @xmath114 , which is according to the same power law as the @xmath92 and @xmath66 case of the lt model ( compare with eq .
( [ eqsmallt ] ) ) .    to summarize our investigation of @xmath88 in the bianchi type i model
we can say that the entity @xmath88 at the initial singularity is a constant determined by the inverse of the dust density . for @xmath110 it will increase immediately after the initial singularity . in most cosmological considerations
it is assumed that the cosmological constant is small .
the exception is in the inflationary era in which the vacuum energy dominates over all other matter degrees of freedom .
the inflationary era will smooth out anisotropies as well as inhomogeneities , and the behavior of @xmath88 in this case is therefore expected . in the de sitter limit @xmath88 will asymptotically move towards a constant .
large @xmath1 means small value , while small @xmath1 corresponds to a large value .
this is in full agreement with the lt models .
it is also interesting that the entity @xmath88 is very sensitive to different matter configurations .
this makes it a lot easier to check whether the @xmath88 has the right entropic behaviour .
the question now arises : how generic is this behaviour ? does the entity behave in the correct way for all physically realistic models ?
we know that our universe today is close to homogeneous on a scale larger than a billion light years .
as mentioned earlier , the frw models are conformally flat and hence , they have a zero weyl tensor .
so a question would be , at late times when the universe is close to isotropic and homogeneous , does the entity @xmath88 still behave in the correct manner ?
does it still increase ?
we saw that at late times , both the lt model and the bianchi type i model , did behave correctly for @xmath66 even though they both isotropise and evolve towards homogeneity .
let us choose a more general model , a flat model that allows for inhomogeneities and anisotropies .
we will assume that the model at late times asymptotically evolve towards a frw model . in @xcite , barrow and maartens
investigated the general equations of motion for such models with @xmath66 . in the velocity dependent regime
they derived some approximate solutions to the field equation close to a frw model , having both inhomogeneities and anisotropy .
also an anisotropic stress of the form @xmath115 where @xmath116 is a constant matrix , was included .
their result was that the late time behaviour of such a model is consistent with a perfect fluid where the equation of state parameter @xmath9 for the perfect fluid obeys @xmath117 .
for @xmath118 the frw model is unstable at late times . in the latter case
one can get for instance different stable anisotropic magnetic solutions @xcite .
let us therefore assume that @xmath119 ( so that the matter obeys the sec ) .
the ratio of the weyl tensor squared and the ricci tensor squared is in this case p^2=o(t^-2n)+o(t^-3n)+o(t^-4n ) where @xmath120 for @xmath121 . in the radiation case , @xmath122 , we get a logarithmic decay of @xmath19 : p^2_=4/3 .
thus in the range of validity of the assumption @xmath123 , this entity will decrease in the future .
however , the entity @xmath88 shows more promising behaviour s = p = a^3 pt^-n [ eq : bmsol]for @xmath124 , while the @xmath122 case yields s(t)^-1t^32 . hence , @xmath88 increases as long as 23 .
this entity increases in the future ( as any entity describing entropy should do ) and tells us that even if the universe itself asymptotically goes towards isotropy , the entropy of the gravitational field actually increases ( if we should believe the wcc ) .
note that this result is slightly different than the case where no anisotropic stress is present .
if the anisotropic stress is not present , we get @xmath125 at late times for all @xmath9 .
we should also mention a work done by hervik @xcite which investigates the evolution of the weyl curvature invariant for generic solutions of spatially homogeneous models containing a @xmath9-law perfect fluid with @xmath126 .
the conclusion was that all spatially homogeneous models , except for sets of measure zero , had an increasing @xmath88 at late times .
as long as the matter obeys the strong energy condition ( sec ) , matter will behave more or less attractive .
hence , when the sec is fulfilled , we should expect the gravitational entropy to increase .    during inflation the sec
is violated and gravity is not necessarily attractive . we have already seen how the inclusion of a cosmological constant could alter the behaviour of @xmath88 .
the late time behaviour of @xmath88 in the presence of a cosmological constant is that @xmath88 evolves approximately as a constant .
the constant itself is a decreasing function of @xmath1 ( see eq .
( [ eq : latetimedesitter ] ) ) .    what happens for a more general inflationary fluid ?
consider a perfect fluid with equation of state p=(-1)for @xmath127 .
this perfect fluid will violate the strong energy condition and will in general cause a power law inflation .
note that the case @xmath14 can be considered the same as including a cosmological constant . in the frw cases ,
all the models will now have a late time behaviour similar to the flat case .
the frw models will become inflationary and will be dominated by this fluid at some stage in the future ( provided that the universe is ever - expanding ) . to simplify
, we will therefore consider a flat universe , and perturb the frw flat universe model with @xmath128 .
following @xcite , the anisotropic stresses will not dominate the shear modes at late times .
the shear will under these assumptions decrease as _ t^- while the hubble parameter is the same as in the frw case ( to lowest order ) h = t^-1 .
the late time behaviour of @xmath19 is now p^2=
o(^2 t^- ) which is decreasing for all @xmath129 .
the late time behaviour of @xmath88 is s = p = o(t ) which is increasing linearly in @xmath35 at late times .
for @xmath14 we recover the cosmological constant case where @xmath88 is approximately a constant at late times .
note also that this result coincide with eq .
( [ eq : bmsol ] ) for @xmath130 .
hence , during the inflationary period , @xmath88 increases much slower than for ordinary matter ( @xmath131 for dust and @xmath132 for radiation ) .
maybe the closest to a quantum theory of gravitation that have been obtained , is what we call quantum cosmology ( qc ) . we will in this section
show how one might be able to determine the likelihood of a certain initial state to occur .
qc is perhaps best described as a `` theory of initial conditions '' .
we will again use the results of the previous paper to try determine which of the initial states are more probable .    in qc
, the wave function of the universe satisfies the _ wheeler - dewitt ( wd ) equation _ : ( -g_ijkl+ v(h_ij))=0 where @xmath133 is called dewitt s supermetric , and @xmath134 is the metric on the 3-dimensional spatial hypersurfaces .
the potential term @xmath135 consists of the ricci scalar of the three - dimensional hypersurfaces and possibly matter potentials and a cosmological constant .
we want now to calculate the expectation value for the weyl scalar , in these models .
all of the curvature invariants have to go over to their respective curvature operators .
especially , the entity @xmath88 goes over to the curvature operator @xmath136 .
the question is now , what is more likely , creation of a universe with large expectation value of @xmath136 or a universe with a low value ?      in the paper @xcite we considered semi - classical tunneling wave functions which are solutions of the wd equation .
the universe was tunneling from a matter - dominated universe classically confined to a finite size , into a de sitter like universe .
after tunneling across the classically forbidden region the universe became @xmath1-dominated , similarly to an inflationary model of the universe .
thus we have a `` flow '' from matter - dominated universes towards de sitter like universes .
the actual expectation values of entities like @xmath137 , @xmath19 and @xmath88 are not explicitly obtained for these models , because the actual calculations suffer from `` endless '' expressions and highly time - consuming quantities . we will therefore give a more general description of the evolution of the weyl tensor for the lt models .
it would be useful first to recapitulate some of the discussion done in @xcite .
first of all we discussed tunneling wave functions in the wkb approximation . in the wkb approximation
we assume that the wave function has the form @xmath138 , where @xmath88 will to the lowest order satisfy the hamilton - jacobi equation : @xmath139=0\ ] ] here @xmath140 denotes the functional derivative with respect to the function @xmath29 .
if we assume that @xmath141 , the resulting equation will be the hamilton - jacobi equation for a point particle with action @xmath142 ( @xmath26 is only a parameter and the functional derivatives turn into ordinary partial derivatives ) . in the hamilton - jacobi equation
the functional @xmath143 turns out to be the action at the classical level .
since the classical action can be written as an integral over @xmath26 the assumption @xmath141 is therefore reasonable at the lowest order wkb level .
we can interpret the action @xmath144 as the action of a point particle moving in a potential @xmath145 $ ] with zero energy .
the wkb wave function @xmath146 for the point particle can then be written @xmath147 .
the two wkb wave functions can therefore be related by @xmath148 .
finding first the wave function @xmath149 we can then relate its wkb approximation to @xmath150 through @xmath148 .
let us now ask the question : _ is it more likely for a universe with small weyl tensor to tunnel through the classical barrier than a universe with a large weyl tensor ? _ the question is difficult to answer in general but we shall make some simple considerations in order to shed some light upon it .
we assume that the dust density near the origin of the coordinates is larger than further out .
we define the homogeneous mass function for a closed universe ( @xmath151 ) as @xmath152 where @xmath153 is a constant .
the constant @xmath153 is determined by demanding @xmath154 . if the dust density is larger near the origin of the @xmath26-coordinate than for larger values of @xmath26 then @xmath155
. this will not in general change the size of the universe so we can look at the effects from @xmath28 alone .
since @xmath28 is greater in general for an inhomogeneous universe than for a homogeneous universe , we see that the potential barrier will be smaller for an inhomogeneous universe than for a homogeneous universe .
thus an inhomogeneous universe will tunnel more easily through the classical barrier than the homogeneous universe . since an inhomogeneous universe will have a larger weyl tensor than an almost homogeneous one
, we see that universes with large weyl tensor tunnel more easily than those with a small weyl tensor .
if we look at the tunneling amplitude concerning effects from the @xmath1 term , it is evident that larger @xmath1 will yield a larger tunneling probability . in the initial era inhomogeneities
will increase the value of @xmath88 .
we saw that an inhomogeneous state will tunnel more easily through the potential barrier than a homogeneous state .
the largest tunneling probability amplitude thus occurs for universes with a large cosmological constant and large local inhomogeneities . from a classical point of view
the value of @xmath88 initially was large ( but increasing thereafter ) , but as the universe entered the inflationary era the cosmological constant had a large value , hence the value of @xmath88 at the end of the inflationary era was relatively small .
in the initial epoch the universe is not believed to be dust dominated .
the dust does not exert any pressure and dust particles do therefore not interact with each other . a more probable matter content is matter which has internal pressure . even though gravitation tends to make the space inhomogeneous , internal pressure from the matter will try to homogenise the space . since dust is the only matter source in our model , the model only indicates the tendency for gravitation itself to create inhomogeneities .
since the universes tunnel into a de sitter - like state , the cosmological constant will rapidly dominate the evolution .
the larger the cosmological constant the lower will the entity @xmath88 be after the inflationary era ends .
we can write the general solution of the wd - equation for the bianchi type i models as : @xmath156\ ] ] where @xmath157 is a particular solution of the wd equation ( with dust ) , @xmath158 is a `` normalizing constant '' , and @xmath159 is a distribution function in momentum space .
this distribution function satisfies the equation : @xmath160 the function @xmath161 is the _
whittaker function_.    there is still a factoring problem in turning the classical entities to operators .
let us first investigate the expectation value of @xmath137 . in the expression for the weyl tensor squared , there is a term ( p_+^2+p_-^2)^2 where @xmath162 is the conjugated momenta to @xmath6 .
upon quantization , the above term is replaced with the operator : ( p_+^2+p_-^2)^2(+ ) ^2 in the expectation value of the square of the weyl tensor , the above term will contribute with @xmath163 ; @xmath164 , the contribution from this term to the weyl curvature invariant will diverge as @xmath17 for small @xmath3 .
thus , we have to conclude that , in the small @xmath3 limit the expectation value of @xmath165 goes as : @xmath166 just as in the classical case .
investigating the invariant @xmath167 , we notice that things are not so easy .
the ricci square also has a term which presumably would contribute with a @xmath168 term .
however , looking at the classical expression we see that the ricci square is independent of the anisotropy parameter @xmath5 .
this indicates that at the classical level all terms involving the anisotropy parameters , have to cancel exactly .
this is not the case quantum mechanically . in the quantum case
operators do not necessarily commute . hence there may be contributions from terms which classically would cancel each other . in other words ,
the fact that the classical vacuum has @xmath104 , does not mean that the quantum vacuum has @xmath169 .
we assume that @xmath170 is a set of factor - ordering parameters which represents the `` true '' quantum mechanical system in such a way that @xmath171 represents the classical system . with this parameterization of the factor ordering we would expect the ricci square expectation value to be : @xmath172 where @xmath173 is some function of @xmath3 which has the property : @xmath174 .
thus for small @xmath3 and @xmath175 the expectation value would behave as @xmath176    and @xmath177 the weyl square divided by the ricci square is in general finite as @xmath178 for a quantum system . in some sense ,
the quantum mechanical effects renormalizes the infinity that the classical system possesses at @xmath16 .
the expectation value at @xmath16 is , however , strongly dependent on the factor - ordering .
as the factor ordering parameters approach zero , the value will diverge .
quantum effects in the early epoch are essential for the behaviour of this entity near the initial singularity .
as @xmath178 we expect the quantum effect to be considerable , thus expecting the factor - ordering parameters @xmath175 to be large .    as indicated in the above discussion , the quantum mechanical expectation value of @xmath179 will be lower and presumably finite at the initial singularity .
therefore the expectation value of @xmath180 is also presumed to be considerably lower in the initial stages than its classical counterpart .
comparing different tunnelling amplitudes in the bianchi type i model is difficult and more speculative because the bianchi type i universe has no classically forbidden region for @xmath181 .
this causes the lowest order wkb approximation to be purely oscillatory . the lowest order wkb wave function will therefore be approximately constant
. in the paper
@xcite we did however construct under some assumptions a wave function which clearly peaked at small values of the anisotropy parameter .
thus these wave functions predicts universes that have a relatively low value of @xmath88 .
let us now recapitulate how the evolution in the context of the wcc _ might _ have been .
the universe was created in a rather arbitrary state . as the time ticked past the planck time @xmath182 s a rather inhomogeneous universe appeared .
qc suggests that this state was rather inhomogeneous , but it is more or less a guess how inhomogeneous the universe was at that state .
nevertheless , as the universe grew larger , the expectation value of the weyl entropy increased initially . at
what rate the weyl entropy increased is highly uncertain , it depends very much on the true nature of our universe .
it depends on what matter fields that were present , the topology of the universe , whether it was anisotropic or not , quantum effects etc .
nevertheless , at some time very short after the big bang , an enormous effective cosmological constant appeared .
the universe was driven unconditionally into an inflationary period .
the weyl entropy stopped increasing and began instead to evolve asymptotically towards a constant .
whether or not the weyl entropy had earlier a higher value than at the exit of the inflationary regime , is difficult to say . during the inflationary period ,
the universe was more or less in an adiabatic expanding state , the weyl entropy was more or less constant .
since the universe was expanding exponentially during this period , the entropy per unit volume dropped exponentially .
if the scale factor increased by a factor of 60 @xmath183-foldings during the inflationary epoch , then the entropy per unit volume of space would have decreased by a factor of ~e^-36010 ^ -78 .
this is quite a drastic decrease , and shows how powerful inflation is when it comes to smoothing out the inhomogeneities and anisotropies of our universe .    at the exit of the inflationary regime ,
the weyl entropy had dropped by an enormous factor compared to what is would have been if no inflation had occurred .
the weyl entropy was very small , compared to the maximal possible value .
the universe was more or less uniform and homogeneous .
this homogeneity can be seen in the cmb radiation today .
however , small fluctuations in the spectrum can also be seen , reflecting the state of the universe 300 000 years after the big bang .
after the inflation , small seeds of inhomogeneities from the quantum fluctuations of quantum fields were the only thing left of the initial inhomogeneities .
nevertheless , these seeds were large enough to gradually clump together and form galaxies and stars .
the weyl entropy began to increase again after the inflationary regime was over .
the radiation- , and later the matter - dominated universe caused the weyl entropy to grow steadily and firmly for almost 15 billion years . during the radiation era , there were only small inhomogeneities and anisotropies left of the primordial fluctuations .
using the calculations from the earlier section , the weyl entropy increases as @xmath132 in the radiation era .    at
about @xmath18410 000 years the radiation become sub - dominant .
the universe evolved into a matter dominated era and the weyl entropy increased as @xmath185 . today
it is still growing .
recent observations suggest that the universe has entered a new era with accelerated expansion .
the surprising fact that the universe appears to be in an accelerating state today , can be explained with the presence of a vacuum energy .
if this is true , the weyl entropy will increase steadily and asymptotically towards a constant .
however , the late time behaviour of our universe is still quite speculative , and whether or not this vacuum dominated period will persist , is very uncertain .
it might happen that the period ends like inflation did , and perhaps the vacuum period is followed by a curvature dominated period .
if so , the weyl entropy might again rise to new heights , increasing towards a value where all the matter is collected in black holes .
sh was founded by the research council of norway .
r. penrose , in _ proc .
first marcel grossmann meet .
( ictp trieste ) _ , ed .
r. ruffini , north - holland ( 1977 ) r. penrose , in _ general relativity , an einstein centenary survey _ ,
s.w . hawking and w. israel , cambridge univ . press ( 1979 ) .
r. penrose , _ j. stat .
phys , _ , * 77 * ( 1994 ) 217    j. wainwright and p.j .
anderson , _ gen .
* 16 * ( 1984 ) 609 .
goode and j. wainwright , _ class .
quantum grav . _
* 2 * ( 1984 ) 99 .
grn and s. hervik , _ class .
quantum grav . _ * 18 * ( 2001 ) 601 .
bonnor , _ class .
grav . _ * 3 * ( 1986 ) 495 .
h. kodama , _ prog .
* 99 * ( 1998 ) 173 j.d . barrow and h. kodama , _ class .
quantum grav . _ * 18 * ( 2001 ) 1753 j.d .
barrow and h. kodama , _ int . j. mod .
_ * d * * 18 * ( 2002 ) h. kodama , ` gr - qc/0109064 ` s. hervik , _ class .
quantum grav . _
* 17 * ( 2000 ) 2765 .
s. hervik , _ class .
quantum grav . _
* 18 * ( 2001 ) 175 .
a. krasiski , _ inhomogeneous cosmological models _ , cambridge university press ( 1997 ) .
a. zecca , _ il nuovo cimento _ * 106b * ( 1990 ) 413 .
b. miller , _ astroph .
j. _ * 208 * ( 1976 ) 275 .
s.w . hawking and r. penrose , _ proc .
_ * a314 * ( 1970 ) 529 .
bonnor , _ phys .
lett . _ * 112a * ( 1985 ) 26 .
s. hayward , _ private communications _ ( 2000 ) .