The Late Devonian was a protracted period of low speciation resulting in biodiversity decline, culminating in extinction events near the Devonian-Carboniferous boundary. Recent evidence indicates that the final extinction event may have coincided with a dramatic drop in stratospheric ozone, possibly due to a global temperature rise. Here we study an alternative possible cause for the postulated ozone drop: a nearby supernova explosion that could inflict damage by accelerating cosmic rays that can deliver ionizing radiation for up to $\sim 100$ kyr. We therefore propose that end-Devonian extinction was triggered by one or more supernova explosions at $\sim 20 \ \rm pc$, somewhat beyond the ``kill distance'' that would have precipitated a full mass extinction. Nearby supernovae are likely due to core-collapses of massive stars in clusters in the thin Galactic disk in which the Sun resides. Detecting any of the long-lived radioisotopes \sm146, \u235 or \pu244 in one or more end-Devonian extinction strata would confirm a supernova origin, point to the core-collapse explosion of a massive star, and probe supernova nucleosythesis. Other possible tests of the supernova hypothesis are discussed.

Deschamps and Mottez (hereafter DM) argue that the Gauss-Matuyama terrestrial magnetic field reversal may have left a vanishing main dipole moment to the field for a time of order 10,000 years. They say this may have allowed an enhanced cosmic ray flux, boosting the effect we proposed in Melott and Thomas (2019). We point out that the bulk of the cosmic ray flux from a nearby supernova should be too energetic, up to a million times more energetic than the limits of deflection by the terrestrial magnetic field. In fact, only those highly energetic ones will directly reach the troposphere, relevant for cloud-to-ground lightning. From Cosmic Explosions to Terrestrial Fires?: A Discussion. F. Deschamps and F. Mottez. J. Geology 128, online ahead of print. (2020) From Cosmic Explosions to Terrestrial Fires?: A Reply A.L. Melott and B.C. Thomas. J. Geology 128, online ahead of print. (2020) From cosmic explosions to terrestrial fires? (A.L. Melott and B.C. Thomas) Journal of Geology, 127, 475-481 10.1086/703418 (2019) [arXiv:1903.01501]

Motivated by the occurrence of a moderately nearby supernova near the beginning of the Pleistocene, possibly as part of a long-term series beginning in the Miocene, we investigate whether nitrate rainout resulting from the atmospheric ionization of enhanced cosmic ray flux could have, through its fertilizer effect, initiated carbon dioxide drawdown. Such a drawdown could possibly reduce the greenhouse effect and induce the climate change that led to the Pleistocene glaciations. We estimate that the nitrogen flux enhancement onto the surface from an event at 50 pc would be of order 10%, probably too small for dramatic changes. We estimate deposition of iron (another potential fertilizer) and find it is also too small to be significant. There are also competing effects of opposite sign, including muon irradiation and reduction in photosynthetic yield caused by UV increase from stratospheric ozone layer depletion, leading to an ambiguous result. However, if the atmospheric ionization induces a large increase in the frequency of lightning, as argued elsewhere, the amount of nitrate synthesis should be much larger, dominate over the other effects, and induce the climate change. More work needs to be done to clarify effects on lightning frequency.

The recent detection of 60Fe in the cosmic rays provides conclusive evidence that there is a recently synthesized component (few MY) in the GCRs (Binns et al. 2016). In addition, these nuclei must have been synthesized and accelerated in supernovae near the solar system, probably in the Sco-Cen OB association subgroups, which are about 100 pc distant from the Sun. Recent theoretical work on the production of r-process nuclei appears to indicate that it is difficult for SNe to produce the solar system abundances relative to iron of r-process elements with high atomic number (Z), including the actinides (Th, U, Np, Pu, and Cm). Instead, it is believed by many that the heaviest r-process nuclei, or perhaps even all r-process nuclei, are produced in binary neutron star mergers. Since we now know that there is at least a component of the GCRs that has been recently synthesized and accelerated, models of r-process production by SNe and BNSM can be tested by measuring the relative abundances of these ultra-heavy r-process nuclei, and especially the actinides, since they are radioactive and provide clocks that give the time interval from nucleosynthesis to detection at Earth. Since BNSM are believed to be much less frequent in our galaxy than SNe (roughly 1000 times less frequent, the ratios of the actinides, each with their own half-life, will enable a clear determination of whether the heaviest r-process nuclei are synthesized in SNe or in BNSM. In addition, the r-process nuclei for the charge range from 34 to 82 can be used to constrain models of r-process production in BNSM and SNe. Thus, GCRs become a multi-messenger component in the study of BNSM and SNe.

There is now solid experimental evidence of at least one supernova explosion within 100 pc of Earth within the last few million years, from measurements of the short-lived isotope 60Fe in widespread deep-ocean samples, as well as in the lunar regolith and cosmic rays. This is the first established example of a specific dated astrophysical event outside the Solar System having a measurable impact on the Earth, offering new probes of stellar evolution, nuclear astrophysics, the astrophysics of the solar neighborhood, cosmic-ray sources and acceleration, multi-messenger astronomy, and astrobiology. Interdisciplinary connections reach broadly to include heliophysics, geology, and evolutionary biology. Objectives for the future include pinning down the nature and location of the established near-Earth supernova explosions, seeking evidence for others, and searching for other short-lived isotopes such as 26Al and 244Pu. The unique information provided by geological and lunar detections of radioactive 60Fe to assess nearby supernova explosions make now a compelling time for the astronomy community to advocate for supporting multi-disciplinary, cross-cutting research programs.

Multiple lines of evidence point to one or more moderately nearby supernovae, with the strongest signal ~2.6 Ma. We build on previous work to argue for the likelihood of cosmic ray ionization of the atmosphere and electron cascades leading to more frequent lightning, and therefore an increase in nitrate deposition and in wildfires. The potential exists for a large increase in the pre-human nitrate flux onto the surface, which has previously been argued to lead to CO2 drawdown and cooling of the climate. Evidence for increased wildfires exists in an increase in soot and carbon deposits over the relevant period. The wildfires would have contributed to the transition from forest to savanna in northeast Africa, long argued to have been a factor in the evolution of hominin bipedalism.

Considerable data and analysis support the detection of one or more supernovae at a distance of about 50 pc, ~2.6 million years ago. This is possibly related to the extinction event around that time and is a member of a series of explosions which formed the Local Bubble in the interstellar medium. We build on previous work, and propagate the muon flux from supernova-initiated cosmic rays from the surface to the depths of the ocean. We find that the radiation dose from the muons will exceed the total present surface dose from all sources at depths up to a kilometer and will persist for at least the lifetime of marine megafauna. It is reasonable to hypothesize that this increase in radiation load may have contributed to a newly documented marine megafaunal extinction at that time.

Recent data indicate one or more moderately nearby supernovae in the early Pleistocene, with additional events likely in the Miocene. This has motivated more detailed computations, using new information about the nature of supernovae and the distances of these events to describe in more detail the sorts of effects that are indicated at the Earth. This short communication/review is designed to describe some of these effects so that they may possibly be related to changes in the biota around these times.

Recent 60Fe results have suggested that the estimated distances of supernovae in the last few million years should be reduced from 100 pc to 50 pc. Two events or series of events are suggested, one about 2.7 million years to 1.7 million years ago, and another may at 6.5 to 8.7 million years ago. We ask what effects such supernovae are expected to have on the terrestrial atmosphere and biota. Assuming that the Local Bubble was formed before the event being considered, and that the supernova and the Earth were both inside a weak, disordered magnetic field at that time, TeV-PeV cosmic rays at Earth will increase by a factor of a few hundred. Tropospheric ionization will increase proportionately, and the overall muon radiation load on terrestrial organisms will increase by a factor of 150. All return to pre-burst levels within 10kyr. In the case of an ordered magnetic field, effects depend strongly on the field orientation. The upper bound in this case is with a largely coherent field aligned along the line of sight to the supernova, in which case TeV-PeV cosmic ray flux increases are 10^4; in the case of a transverse field they are below current levels. We suggest a substantial increase in the extended effects of supernovae on Earth and in the lethal distance estimate; more work is needed.This paper is an explicit followup to Thomas et al. (2016). We also here provide more detail on the computational procedures used in both works.

Frolop and Scott (2016) claim significant 1-1 correspondence between anomalies in the cosmic microwave background (CMB) and the digits of pi, which they call 'Pi in the Sky'. They have without attribution republished the famous work of Joe Hill (Hill 1911), who first proposed this idea, then repudiated it.

Major discrepancies have been noted for some time between fossil ages and molecular divergence dates for a variety of taxa. Recently, systematic trends within avian clades have been uncovered. The trends show that the disparity is much larger for mitochondrial DNA than for nuclear DNA, also that it is larger for crown fossil dates than stem fossil dates. It has been argued that this pattern is largely inconsistent with incompleteness of the fossil record as the principal driver of the disparity. A case is presented that given the expected mutations from a fluctuating background of astrophysical radiation from such sources as supernovae, the rate of molecular clocks is variable and should increase back to a few Ma, before returning to the long-term average rate. This is a possible explanation for the disparity. One test of this hypothesis is to look for an acceleration of molecular clocks 2 to 2.5 Ma due to one or more moderately nearby supernovae known to have happened at that time. Another is to look for reduced disparity in benthic organisms of the deep ocean. In addition, due to the importance of highly penetrating muon irradiation, the disparity should be magnified for megafauna.

The record of reversals of the geomagnetic field has played an integral role in the development of plate tectonic theory. Statistical analyses of the reversal record are aimed at detailing patterns and linking those patterns to core-mantle processes. The geomagnetic polarity timescale is a dynamic record and new paleomagnetic and geochronologic data provide additional detail. In this paper, we examine the periodicity revealed in the reversal record back to 375 million years ago (Ma) using Fourier analysis. Four significant peaks were found in the reversal power spectra within the 16-40-million-year range (Myr). Plotting the function constructed from the sum of the frequencies of the proximal peaks yield a transient 26 Myr periodicity, suggesting chaotic motion with a periodic attractor. The possible 16 Myr periodicity, a previously recognized result, may be correlated with pulsation of mantle plumes and perhaps; more tentatively, with core-mantle dynamics originating near the large low shear velocity layers in the Pacific and Africa. Planetary magnetic fields shield against charged particles which can give rise to radiation at the surface and ionize the atmosphere, which is a loss mechanism particularly relevant to M stars. Understanding the origin and development of planetary magnetic fields can shed light on the habitable zone.

Duderstadt et al. [2016b] comment that the Melott et al. [2016] study of nitrate formation by solar proton events (SPEs) and comparison with the ice core archive is "fundamentally flawed," because it does not include pre-existing HNO3 in the stratosphere. We show that they exaggerate both the enhancement predicted by our findings and pre-industrial HNO3 levels in their model, and fail to prove this assertion. Our feasibility study matched expected SPE nitrate production with ground truth measurements. It is not clear that their approach is more realistic and absence of a detailed mechanism does not disprove our results. Models can be no better than the information they are provided and in this case there continue to be significant unknowns and uncertainties, especially in the role of polar stratospheric clouds (PSCs) and possible interactions with cosmic rays that constitute lower boundary conditions. Duderstadt et al. [2014; 2016a] used incomplete, poorly-constrained and incorrect initial and boundary conditions, and they continue to advocate on the basis of uncertain results. Meanwhile, Smart et al. [2014] identified a series of ice core nitrate spikes that have since been confirmed in 10Be by McCracken and Beer [2015]. Melott et al. [2016] computationally reproduced the ionization profile of the only major balloon measurements to date. We show that our calculated nitrate enhancement is consistent with measured results, given current levels of uncertainty, and that extreme SPEs can potentially produce occasional nitrate spikes with hundreds of percent increases. Instead of repeating old arguments to dismiss nitrates as proxies of SPEs, it is past time for a dedicated, fine-resolution, multi-parameter, replicate ice core field campaign to resolve this debate.

Recent results have strongly confirmed that multiple supernovae happened at distances ~100 pc consisting of two main events: one at 1.7 to 3.2 million years ago, and the other at 6.5 to 8.7 million years ago. These events are said to be responsible for excavating the Local Bubble in the interstellar medium and depositing 60Fe on Earth and the Moon. Other events are indicated by effects in the local cosmic ray (CR) spectrum. Given this updated and refined picture, we ask whether such supernovae are expected to have had substantial effects on the terrestrial atmosphere and biota. In a first cut at the most probable cases, combining photon and cosmic ray effects, we find that a supernova at 100 pc can have only a small effect on terrestrial organisms from visible light and that chemical changes such as ozone depletion are weak. However, tropospheric ionization right down to the ground due to the penetration of $\geq$TeV cosmic rays will increase by nearly an order of magnitude for thousands of years, and irradiation by muons on the ground and in the upper ocean will increase 20-fold, which will approximately triple the overall radiation load on terrestrial organisms. Such irradiation has been linked to possible changes in climate and increased cancer and mutation rates. This may be related to a minor mass extinction around the Pliocene-Pleistocene boundary, and further research on the effects is needed.

Solar energetic particles ionize the atmosphere, leading to production of nitrogen oxides. It has been suggested that some such events are visible as layers of nitrate in ice cores, yielding archives of energetic, high fluence solar proton events (SPEs). There has been controversy, due to slowness of transport for these species down from the upper stratosphere; past numerical simulations based on an analytic calculation have shown very little ionization below the mid stratosphere. These simulations suffer from deficiencies: they consider only soft SPEs and narrow energy ranges; spectral fits are poorly chosen; with few exceptions secondary particles in air showers are ignored. Using improved simulations that follow development of the proton-induced air shower, we find consistency with recent experiments showing substantial excess ionization down to 5 km. We compute nitrate available from the 23 February 1956 SPE, which had a high fluence, hard spectrum, and well-resolved associated nitrate peak in a Greenland ice core. For the first time, we find this event can account for ice core data with timely (~ 2 months) transport downward between 46 km and the surface, thus indicating an archive of high fluence, hard spectrum SPE covering the last several millennia. We discuss interpretations of this result, as well as the lack of a clearly-defined nitrate spike associated with the soft-spectrum 3-4 August 1972 SPE. We suggest that hard-spectrum SPEs, especially in the 6 months of polar winter, are detectable in ice cores, and that more work needs to be done to investigate this.

This work by Feng and papers which published its conclusions do not cite nor do they deal with objections by the author published in 2013-2014. There are many fundamental problems. We will summarize here the principal problems, as published by Melott and Bambach which render irrelevant most of the work presented by Feng.

It is known that a gamma ray burst (GRB) originating near the Earth could be devastating to life. The mechanism of ozone depletion and subsequent increased UVB exposure is the primary risk, but models also show increased nitrification culminating in nitric acid rainout. These effects are also expected after nearby supernovae and extreme solar proton events. In this work we considered specifically whether the increased nitric acid rainout from such events is a threat to modern terrestrial ecosystems. We also considered its potential benefit to early terrestrial Paleozoic ecosystems. We used established critical loads for nitrogen deposition in ecoregions of Europe and the US and compared them with previously predicted values of nitric acid rainout from a typical GRB within our galaxy. The predicted rainout was found to be too low to harm modern ecosystems, however, it is large compared with probable nitrate flux onto land prior to the invasion of plants. We suggest that this flux may have contributed nutrients to this invasion if, as hypothesized, the end-Ordovician extinction event were initiated by a GRB or other ionizing radiation event.

Cosmic rays are known to cause biological effects directly and through ionizing radiation produced by their secondaries. These effects have been detected in airline crews and other specific cases where members of the population are exposed to above average secondary fluxes. Recent work has found a correlation between solar particle events and congenital malformations. In this work we use the results of computational simulations to approximate the ionizing radiation from such events as well as longer term increases in cosmic ray flux. We find that the amounts of ionizing radiation produced by these events are insufficient to produce congenital malformations under the current paradigm regarding muon ionizing radiation. We believe that further work is needed to determine the correct ionizing radiation contribution of cosmogenic muons. We suggest that more extensive measurements of muon radiation effects may show a larger contribution to ionizing radiation dose than currently assumed.

Ice cores are archives of climate change and possibly large solar proton events (SPEs). Wolff et al. (2012) used a single event, a nitrate peak in the GISP2-H core, which McCracken et al. (2001a) time associated with the poorly quantified 1859 Carrington event, to discredit SPE-produced, impulsive nitrate deposition in polar ice. This is not the ideal test case. We critique the Wolff et al. analysis and demonstrate that the data they used cannot detect impulsive nitrate events because of resolution limitations. We suggest re-examination of the top of the Greenland ice sheet at key intervals over the last two millennia with attention to fine resolution and replicate sampling of multiple species. This will allow further insight into polar depositional processes on a sub-seasonal scale, including atmospheric sources, transport mechanisms to the ice sheet, post-depositional interactions, and a potential SPE association.

Firestone (2014) asserted evidence for numerous (23) nearby (d<300 pc) supernovae within the Middle and Late Pleistocene. If true, this would have strong implications for the irradiation of the Earth; at this rate, mass extinction level events due to supernovae would be more frequent than 100 Myr. However, there are numerous errors in the application of past research. The paper overestimates likely nitrate and 14C production from moderately nearby supernovae by about four orders of magnitude. Moreover, the results are based on wrongly selected (obsolete) nitrate and 14C datasets. The use of correct and up-to-date datasets does not confirm the claimed results. The claims in the paper are invalidated.

"The investigation into the possible effects of cosmic rays on living organisms will also offer great interest." - Victor F. Hess, Nobel Lecture, December 12, 1936 High-energy radiation bursts are commonplace in our Universe. From nearby solar flares to distant gamma ray bursts, a variety of physical processes accelerate charged particles to a wide range of energies, which subsequently reach the Earth. Such particles contribute to a number of physical processes occurring in the Earth system. A large fraction of the energy of charged particles gets deposited in the atmosphere, ionizing the atmosphere, causing changes in its chemistry and affecting the global electric circuit. Remaining secondary particles contribute to the background dose of cosmic rays on the surface and parts of the subsurface region. Life has evolved over the past ~ 3 billion years in presence of this background radiation, which itself has varied considerably during the period. As demonstrated by the Miller-Urey experiment, lightning plays a very important role in the formation of complex organic molecules, which are the building blocks of more complex structures forming life. There is growing evidence of increase in the lightning rate with increasing flux of charged particles. Is there a connection between enhanced rate of cosmic rays and the origin of life? Cosmic ray secondaries are also known to damage DNA and cause mutations, leading to cancer and other diseases. It is now possible to compute radiation doses from secondary particles, in particular muons and neutrons. Have the variations in cosmic ray flux affected the evolution of life on earth? We describe the mechanisms of cosmic rays affecting terrestrial life and review the potential implications of the variation of high-energy astrophysical radiation on the history of life on earth.

Gamma ray bursts (GBRs) are a class of highly energetic explosions emitting radiation in a very short timescale of a few seconds and with a very narrow opening angle. Although, all GRBs observed so far are extragalactic in origin, there is a high probability of a GRB of galactic origin beaming towards the Earth in the past ~ 0.5 Gyr. Such an intense burst of gamma rays would ionize the atmosphere and deplete the ozone layer. With depleted ozone, there will be an increased flux of solar UVB on the Earth's surface with harmful biological effects. In addition to the atmospheric damage, secondary particles produced by gamma ray-induced showers will reach the surface. Amongst all secondary particles, muons dominate the ground-level secondary particle flux (99% of the total number of particles) and are potentially of biological significance. Using the Monte Carlo simulation code CORSIKA, we modeled the air showers produced by gamma ray primaries up to 100 GeV. We found that the number of muons produced by hypothetical galactic GRBs significantly increases the total muon flux. However, since the muon production efficiency is extremely low for photon energy below 100 GeV, and because GRBs radiate strongly for only a very short time, we find that the biological radiation dose from secondary muons is negligible. The main mechanism of biological damage from GRBs is through solar UVB irradiation from the loss of ozone in the upper atmosphere.

Unfortunately, Liu et al. contains a number of errors and omissions which compromise its conclusions. These have to do with the amount of 14C which is necessary to deposit in the atmosphere in order to see a perturbation like that in 774 AD, and the ability of a comet to do so.

There may be structural principles pertaining to the general behavior of systems that lead to similarities in a variety of different contexts. Classic examples include the descriptive power of fractals, the importance of surface area to volume constraints, the universality of entropy in systems, and mathematical rules of growth and form. Documenting such overarching principles may represent a rejoinder to the Neodarwinian synthesis that emphasizes adaptation and competition. Instead, these principles could indicate the importance of constraint and structure on form and evolution. Here we document a potential example of a phenomenon suggesting congruent behavior of very different systems. We focus on the notion that universally there has been a tendency for more volatile entities to disappear from systems such that the net volatility in these systems tends to decline. We specifically focus on origination and extinction rates in the marine animal fossil record, the performance of stocks in the stock market, and the characters of stars and stellar systems. We consider the evidence that each is experiencing declining volatility, and also consider the broader significance of this.

Analysis of two independent data sets with increased taxonomic resolution (genera rather than families) using the revised 2012 time scale reveals that an extinction periodicity first detected by Raup and Sepkoski (1984) for only the post-Paleozoic actually runs through the entire Phanerozoic. Although there is not a local peak of extinction every 27 million years, an excess of the fraction of genus extinction by interval follows a 27 million year timing interval and differs from a random distribution at the p ~ 0.02 level. A 27 million year periodicity in the spectrum of interval lengths no longer appears, removing the question of a possible artifact arising from it. Using a method originally developed in Bambach (2006) we identify 19 intervals of marked extinction intensity, including mass extinctions, spanning the last 470 million years (and with another six present in the Cambrian) and find that 10 of the 19 lie within 3 Myr of the maxima in the spacing of the 27 Myr periodicity, which differs from a random distribution at the p = 0.004 level. These 19 intervals of marked extinction intensity also preferentially occur during decreasing diversity phases of a well-known 62 Myr periodicity in diversity (16 of 19, p = 0.002). Both periodicities appear to enhance the likelihood of increased severity of extinction, but the cause of neither periodicity is known. Variations in the strength of the many suggested causes of extinction coupled to the variation in combined effect of the two different periodicities as they shift in and out of phase is surely one of the reasons that definitive comparative study of the causes of major extinction events is so elusive.

We explore the idea that detectable excursions in 26Al may arise from direct deposition by any bolide, and excursions in 14C and 10Be abundances in the atmosphere may result from long-period comet impacts. This is very different from the usual processes of production by cosmic rays within Earths atmosphere. Long-period comets experience greatly increased cosmic ray flux beyond the protection of the suns magnetic field. We report the computed amount of 14C, 10Be, and 26Al present on long-period comets as a function of comet mass. We find that the amount of nuclide mass on large long-period comets entering the Earths atmosphere may be sufficient for creating anomalies in the records of 14C and 10Be from past impacts. In particular, the estimated mass of the proposed Younger Dryas comet is consistent with its having deposited sufficient isotopes to account for recorded 14C and 10Be increases at that time. The 26Al/10Be ratio is much larger in extraterrestrial objects than in the atmosphere, and so, we note that measuring this ratio in ice cores is a suitable definitive test for the Younger Dryas impact hypothesis, even if the hypothetical bolide is not a long-period comet and/or did not contribute to the 14C and 10Be increases.

A major revision of the geological timescale was published in 2012. We re-examine our past finding of a 27 Myr periodicity in marine extinction rates by re-assigning dates to the extinction data used previously. We find that the spectral power at this period is somewhat increased, and persists at a narrow bandwidth, which supports our previous contention that the Nemesis hypothesis is untenable as an explanation for the periodicity that was first noted by Raup and Sepkoski in the 1980s. We enumerate a number of problems in a recent study comparing extinction rates with time series models.

Under current conditions, the cosmic ray spectrum incident on the Earth is dominated by particles with energies < 1 GeV. Astrophysical sources including high energy solar flares, supernovae and gamma ray bursts produce high energy cosmic rays (HECRs) with drastically higher energies. The Earth is likely episodically exposed to a greatly increased HECR flux from such events, some of which lasting thousands to millions of years. The air showers produced by HECRs ionize the atmosphere and produce harmful secondary particles such as muons and neutrons. Neutrons currently contribute a significant radiation dose at commercial passenger airplane altitude. With higher cosmic ray energies, these effects will be propagated to ground level. This work shows the results of Monte Carlo simulations quantifying the neutron flux due to high energy cosmic rays at various primary energies and altitudes. We provide here lookup tables that can be used to determine neutron fluxes from primaries with total energies 1 GeV - 1 PeV. By convolution, one can compute the neutron flux for any arbitrary CR spectrum. Our results demonstrate that deducing the nature of primaries from ground level neutron enhancements would be very difficult.

We examine possible sources of a substantial increase in tree ring 14C measurements for the years AD 774-775. Contrary to claims regarding a coronal mass ejection (CME), the required CME energy is not several orders of magnitude greater than known solar events. We consider solar proton events (SPEs) with three different fluences and two different spectra. The data may be explained by an event with fluence about one order of magnitude beyond the October 1989 SPE. Two hard spectrum cases considered here result in moderate ozone depletion, so no mass extinction is implied, though we do predict increases in erythema and damage to plants from enhanced solar UV. We are able to rule out an event with a very soft spectrum that causes severe ozone depletion and subsequent biological impacts. Nitrate enhancements are consistent with their apparent absence in ice core data. The modern technological implications of such an event may be extreme, and considering recent confirmation of superflares on solar-type stars, this issue merits attention.

Atmospheric 14C production is a potential window into the energy of solar proton and other cosmic ray events. It was previously concluded that results from AD 774-775 are orders of magnitude greater than known solar events. We find that the coronal mass ejection energy based on 14C production is much smaller than claimed, but still substantially larger than the maximum historical Carrington Event of 1859. Such an event would cause great damage to modern technology, and in view of recent confirmation of superflares on solar-type stars, this issue merits attention.

One of the longstanding debates in the history of paleontology focuses on the issue of whether or not there have been long term cycles (operating over tens of millions of years) in biodiversity and extinction. Here we consider the history of this debate by connecting the skein from Grabau up to 2008. We focus on the evidence for periodicity that has emerged thus far, and conclude that there is indeed some evidence that periodicity may be real, though of course more work is needed. We also comment on possible causal mechanisms, focusing especially on the motion of our solar system in the Galaxy. Moreover, we consider the reasons why some scientists have opposed periodicity over the years. Finally, we consider the significance of this for our understanding of evolution and the history of life.

We find that the marine 87Sr/86Sr record shows a significant periodicity of 59.3 \pm 3 Myr. The 87Sr/86Sr record is 171{\deg} \pm 12{\deg}out of phase with a 62 (\pm 3) Myr periodicity previously reported in the record of marine-animal diversity. These periodicities are close to 58 (\pm 4) Myr cycles found for the number of gap-bounded sedimentary carbonate packages of North America We propose that these periodicities reflect the operation of a periodic "pulse of the Earth" in large-scale, Earth processes. These may be linked to mantle or plate-tectonic events, possibly uplift, which affects Earth's climate and oceans, and so the geochemistry, sedimentation, and biodiversity of the marine realm. Alternately, they may be linked to oscillation of the solar system normal to the plane of the galaxy.

A ~ 62 My periodicity in fossil biodiversity has been observed in independent studies of paleontology databases over ~0.5Gy. The period and phase of this biodiversity cycle coincides with the oscillation of our solar system normal to the galactic disk with an amplitude ~70 parsecs and a period ~64 My. Our Galaxy is falling toward the Virgo cluster, forming a galactic shock at the north end of our galaxy due to this motion, capable of accelerating particles and exposing our galaxy's northern side to a higher flux of cosmic rays. These high-energy particles strike the Earth's atmosphere initiating extensive air showers, ionizing the atmosphere by producing charged secondary particles. Secondary particles such as muons produced as a result of nuclear interactions are able to reach the ground and enhance the biological radiation dose. Using a Monte Carlo simulation package CORSIKA, we compute the biological dose resulting from enhanced muon exposure from cosmic rays and discuss their implications for terrestrial biodiversity variations.

On geological timescales, the Earth is likely to be exposed to an increased flux of high energy cosmic rays (HECRs) from astrophysical sources such as nearby supernovae, gamma ray bursts or by galactic shocks. Typical cosmic ray energies may be much higher than the ~ 1 GeV flux which normally dominates. These high-energy particles strike the Earth's atmosphere initiating an extensive air shower. As the air shower propagates deeper, it ionizes the atmosphere by producing charged secondary particles. Secondary particles such as muons and thermal neutrons produced as a result of nuclear interactions are able to reach the ground, enhancing the radiation dose. Muons contribute 85% to the radiation dose from cosmic rays. This enhanced dose could be potentially harmful to the biosphere. This mechanism has been discussed extensively in literature but has never been quantified. Here, we have developed a lookup table that can be used to quantify this effect by modeling terrestrial muon flux from any arbitrary cosmic ray spectra with 10 GeV - 1 PeV primaries. This will enable us to compute the radiation dose on terrestrial planetary surfaces from a number of astrophysical sources.

It is shown that the historical summary of the growth in size of N-body simulations as measured by particle number in this review is missing some key milestones. Size matters, because particle number with appropriate force smoothing is a key method to suppress unwanted discreteness, so that the initial conditions and equations of motion are appropriate to growth by gravitational instability in a Poisson-Vlasov system appropriate to a Universe with dark matter. Published strong constraints on what can be done are not included in the review.

Cosmic radiation backgrounds are a constraint on life, and their distribution will affect the Galactic Habitable Zone. Life on Earth has developed in the context of these backgrounds, and characterizing event rates will elaborate the important influences. This in turn can be a base for comparison with other potential life-bearing planets. In this review we estimate the intensities and rates of occurrence of many kinds of strong radiation bursts by astrophysical entities ranging from gamma-ray bursts at cosmological distances to the Sun itself. Many of these present potential hazards to the biosphere: on timescales long compared with human history, the probability of an event intense enough to disrupt life on the land surface or in the oceans becomes large. We enumerate the known sources of radiation and characterize their intensities at the Earth and rates or upper limits on these quantities. When possible, we estimate a "lethal interval", our best estimate of how often a major extinction-level event is probable given the current state of knowledge; we base these estimates on computed or expected depletion of stratospheric ozone. In general, moderate level events are dominated by the Sun, but the far more severe infrequent events are probably dominated by gamma-ray bursts and supernovae. We note for the first time that so-called "short-hard" gamma-ray bursts are a substantial threat, comparable in magnitude to supernovae and greater than that of the higher-luminosity long bursts considered in most past work. Given their precursors, short bursts may come with little or no warning.

We investigate evolutionary dynamics related to periodicity fossil biodiversity. Coherent periodic fluctuation in origination/extinction of marine genera that survive <45 million years is the source of an observed ~62 million year periodicity analyzed in Paper I. We also show that the evolutionary dynamics of "long-lived" genera (those that survive >45 million years) do not participate in the periodic fluctuation in diversity and differ from those of "short-lived" genera. The difference between the evolutionary dynamics of these 2 genera classes indicates that the periodic pattern is not an artifact of variation in quality of the geologic record. The interplay of these two previously undifferentiated systems, together with the secular increase in abundance of "long-lived" genera, is probably the source of heretofore unexplained differences in evolutionary dynamics between the Paleozoic and post-Paleozoic as reported by others. Testing for cycles similar to the 62 Myr cycle in fossil biodiversity superimposed on the long-term trends of the Phanerozoic as described in Paper I, we find a significant (but weaker) signal in sedimentary rock packages, particularly carbonates, which suggests a connection. The presence of a periodic pattern in evolutionary dynamics of the vulnerable "short-lived" component of marine fauna demonstrates that a long-term periodic fluctuation in environmental conditions capable of affecting evolution in the marine realm characterizes our planet. Coincidence in timing is more consistent with a common cause than sampling bias. A previously identified set of mass extinctions preferentially occur during the declining phase of the 62 Myr periodicity, supporting the idea that the periodicity relates to variation in biotically important stresses. Further work should focus on finding links to physical phenomena that might reveal the causal system or systems.

We use Fourier analysis and related techniques to investigate the question of periodicities in fossil biodiversity. These techniques are able to identify cycles superimposed on the long-term trends of the Phanerozoic. We review prior results and analyze data previously reduced and published. Joint time-series analysis of various reductions of the Sepkoski Data, Paleobiology Database, and Fossil Record 2 indicate the same periodicity in biodiversity of marine animals at 62 Myr. We have not found this periodicity in the terrestrial fossil record. We have found that the signal strength decreases with time because of the accumulation of apparently "resistant" long-lived genera. The existence of a 62-Myr periodicity despite very different treatment of systematic error, particularly sampling-strength biases, in all three major databases strongly argues for its reality in the fossil record.

The hypothesis of a companion object (Nemesis) orbiting the Sun was motivated by the claim of a terrestrial extinction periodicity, thought to be mediated by comet showers. The orbit of a distant companion to the Sun is expected to be perturbed by the Galactic tidal field and encounters with passing stars, which will induce variation in the period. We examine the evidence for the previously proposed periodicity, using two modern, greatly improved paleontological datasets of fossil biodiversity. We find that there is a narrow peak at 27 My in the cross-spectrum of extinction intensity time series between these independent datasets. This periodicity extends over a time period nearly twice that for which it was originally noted. An excess of extinction events are associated with this periodicity at 99% confidence. In this sense we confirm the originally noted feature in the time series for extinction. However, we find that it displays extremely regular timing for about 0.5 Gy. The regularity of the timing compared with earlier calculations of orbital perturbation would seem to exclude the Nemesis hypothesis as a causal factor.

The hypothesis of a companion object (Nemesis) orbiting the Sun was motivated by the claim of a terrestrial extinction periodicity, thought to be mediated by comet showers. The orbit of a distant companion to the Sun is expected to be perturbed by the Galactic tidal field and encounters with passing stars, which will induce variation in the period. We examine the evidence for the previously proposed periodicity, using two modern, greatly improved paleontological datasets of fossil biodiversity. We find that there is a narrow peak at 27 My in the cross-spectrum of extinction intensity time series between these independent datasets. This periodicity extends over a time period nearly twice that for which it was originally noted. An excess of extinction events are associated with this periodicity at 99% confidence. In this sense we confirm the originally noted feature in the time series for extinction. However, we find that it displays extremely regular timing for about 0.5 Gy. The regularity of the timing compared with earlier calculations of orbital perturbation would seem to exclude the Nemesis hypothesis as a causal factor.

Variation in high energy cosmic rays (HECRs) has been proposed to explain a 62 My periodicity in terrestrial fossil biodiversity. It has been suggested that the infall of our galaxy toward the Virgo cluster could generate an extragalactic shock, accelerating charged particles and exposing the earth to a flux of high energy cosmic rays (HECRs). The oscillation of the Sun perpendicular to the galactic plane could induce 62 My periodicity in the HECR flux on the Earth, with a magnitude much higher than the Galactic cosmic ray change we see in a solar cycle. This mechanism could potentially explain the observed 62 My periodicity in terrestrial biodiversity over the past 500 My. In addition to direct effects on life from secondaries, HECRs induced air showers ionize the atmosphere leading to changes in atmospheric chemistry and microphysical processes that can lead to cloud formation including low altitude cloud cover. An increase in ionization changes the global electric circuit which could enhance the formation of cloud condensation nuclei (CCN) through microphysical processes such as electroscavenging and ion mediated nucleation, leading to an increase in the cloud cover. This could increase the albedo and reduce the solar flux reaching the ground, reducing the global temperature. Using an existing model, we have calculated the enhancement in atmospheric ionization at low altitudes resulting from exposure to HECRs. We use a conservative model to estimate the change in low altitude cloud cover from this increased ionization.

We find agreement between models of atmospheric chemistry changes from ionization for the 1908 Tunguska airburst event and nitrate enhancement in GISP2H and GISP2 ice cores, plus an unexplained ammonium spike. We then consider a candidate cometary impact at the Younger Dryas onset (YD). The large estimated NO_x production and O_3 depletion are beyond accurate extrapolation, but the ice core peak is much lower, possibly because of insufficient sampling resolution. Ammonium and nitrate spikes have been attributed to biomass burning at YD onset in both GRIP and GISP2 ice cores. A similar result is well-resolved in Tunguska ice core data, but that forest fire was far too small to account for this. Direct input of ammonia from a comet into the atmosphere is adequate for YD ice core data, but not Tunguska data. An analog of the Haber process with hydrogen contributed by cometary or surface water, atmospheric nitrogen, high pressures, and possibly catalytic iron from a comet could in principle produce ammonia, accounting for the peaks in both data sets.

A variety of events such as gamma-ray bursts and supernovae may expose the Earth to an increased flux of high-energy cosmic rays, with potentially important effects on the biosphere. Existing atmospheric chemistry software does not have the capability of incorporating the effects of substantial cosmic ray flux above 10 GeV . An atmospheric code, the NASA-Goddard Space Flight Center two-dimensional (latitude, altitude) time-dependent atmospheric model (NGSFC), is used to study atmospheric chemistry changes. Using CORSIKA, we have created tables that can be used to compute high energy cosmic ray (10 GeV - 1 PeV) induced atmospheric ionization and also, with the use of the NGSFC code, can be used to simulate the resulting atmospheric chemistry changes. We discuss the tables, their uses, weaknesses, and strengths.

It has been suggested that galactic shock asymmetry induced by our galaxy's infall toward the Virgo Cluster may be a source of periodicity in cosmic ray exposure as the solar system oscillates perpendicular to the galactic plane. Here we investigate a mechanism by which cosmic rays might affect terrestrial biodiversity, ionization and dissociation in the atmosphere, resulting in depletion of ozone and a resulting increase in the dangerous solar UVB flux on the ground, with an improved ionization background computation averaged over a massive ensemble (about 7 x 10^5) shower simulations. We study minimal and full exposure to the postulated extragalactic background. The atmospheric effects are greater than with our earlier, simplified ionization model. At the lower end of the range effects are too small to be of serious consequence. At the upper end of the range, ~6 % global average loss of ozone column density exceeds that currently experienced due to effects such as accumulated chlorofluorocarbons. The intensity is less than a nearby supernova or galactic gamma-ray burst, but the duration would be about 10^6 times longer. Present UVB enhancement from current ozone depletion ~3% is a documented stress on the biosphere, but a depletion of the magnitude found at the upper end of our range would double the global average UVB flux. For estimates at the upper end of the range of the cosmic ray variability over geologic time, the mechanism of atmospheric ozone depletion may provide a major biological stress, which could easily bring about major loss of biodiversity. Future high energy astrophysical observations will resolve the question of whether such depletion is likely.

A 62 Myr periodicity is superimposed on other longer-term trends in fossil biodiversity. This cycle can be discerned in marine data based on the Sepkoski compendium, the Paleobiology Database, and the Fossil Record 2. The signal also exists in changes in sea level/sediment, but is much weaker than in biodiversity itself. A significant excess of 19 previously identified Phanerozoic mass extinctions occur on the declining phase of the 62 Myr cycle. appearance of the signal in sampling-standardized biodiversity data, it is likely not to be a sampling artifact, but either a consequence of sea-level changes or an additional effect of some common cause for them both. In either case, it is intriguing why both changes would have a regular pattern.

We re-examine past suggestions of a close link between terrestrial climate change and the Sun's transit of spiral arms in its path through the Milky Way galaxy. These links produced concrete fits, deriving the unknown spiral pattern speed from terrestrial climate correlations. We test these fits against new data on spiral structure based on CO data that does not make simplifying assumptions about symmetry and circular rotation. If we compare the times of these transits to changes in the climate of Earth, not only do the claimed correlations disappear, but also we find that they cannot be resurrected for any reasonable pattern speed.

Based on the intensity and rates of various kinds of intense ionizing radiation events such as supernovae and gamma-ray bursts, it is likely that the Earth has been subjected to one or extinction level events during the Phanerozoic. These induce changes in atmospheric chemistry so that the level of Solar ultraviolet-B radiation reaching the surface and near-surface waters may be doubled for up to a decade. This UVB level is known from experiment to be more than enough to kill off many kinds of organisms, particularly phytoplankton. It could easily induce a crash of the photosynthetic-based food chain in the oceans. Regularities in the latitudinal distribution of damage are apparent in simulations of the atmospheric changes. We previously proposed that the late Ordovician extinction is a plausible candidate for a contribution from an ionizing radiation event, based on environmental selectivity in trilobites. To test a null hypothesis based on this proposal, we confront latitudinal differential extinction rates predicted from the simulations with data from a published analysis of latitudinal gradients in the Ordovician extinction. The pattern of UVB damage always shows a strong maximum at some latitude, with substantially lower intensity to the north and south of this maximum. We find that the pattern of damage predicted from our simulations is consistent with the data assuming a burst approximately over the South Pole, and no further north than -75 degrees. We predict that any land mass (such as parts of north China, Laurentia, and New Guinea) which then lay north of the equator should be a refuge from UVB effects, and show a different pattern of extinction in the first strike of the end-Ordovician extinction, if induced by such a radiation event.

Time series analysis of fossil biodiversity of marine invertebrates in the Paleobiology Database (PBDB) shows a significant periodicity at approximately 63 My, in agreement with previous analyses based on the Sepkoski database. I discuss how this result did not appear in a previous analysis of the PBDB. The existence of the 63 My periodicity, despite very different treatment of systematic error in both PBDB and Sepkoski databases strongly argues for consideration of its reality in the fossil record. Cross-spectral analysis of the two datasets finds that a 62 My periodicity coincides in phase by 1.6 My, equivalent to better than the errors in either measurement. Consequently, the two data sets not only contain the same strong periodicity, but its peaks and valleys closely correspond in time. Two other spectral peaks appear in the PBDB analysis, but appear to be artifacts associated with detrending and with the increased interval length. Sampling-standardization procedures implemented by the PBDB collaboration suggest that the signal is not an artifact of sampling bias. Further work should focus on finding the cause of the 62 My periodicity.

Essential equivalence of the conclusions of Romeo et al. (arXiv:0804.0294v1) with earlier work is pointed out. A possible general explanation for these conclusions is suggested.

Data on the prevalence of bone cancer in dinosaurs is available from past radiological examination of preserved bones. We statistically test this data for consistency with rates extrapolated from information on bone cancer in modern vertebrates, and find that there is no evidence of a different rate. Thus, this test provides no support for a possible role of ionizing radiation in the K-T extinction event.

Recently Supernova 2006gy was noted as the most luminous ever recorded, with a total radiated energy of ~10^44 Joules. It was proposed that the progenitor may have been a massive evolved star similar to eta Carinae, which resides in our own galaxy at a distance of about 2.3 kpc. eta Carinae appears ready to detonate. Although it is too distant to pose a serious threat as a normal supernova, and given its rotation axis is unlikely to produce a Gamma-Ray Burst oriented toward the Earth, eta Carinae is about 30,000 times nearer than 2006gy, and we re-evaluate it as a potential superluminous supernova. We find that given the large ratio of emission in the optical to the X-ray, atmospheric effects are negligible. Ionization of the atmosphere and concomitant ozone depletion are unlikely to be important. Any cosmic ray effects should be spread out over ~10^4 y, and similarly unlikely to produce any serious perturbation to the biosphere. We also discuss a new possible effect of supernovae, endocrine disruption induced by blue light near the peak of the optical spectrum. This is a possibility for nearby supernovae at distances too large to be considered "dangerous" for other reasons. However, due to reddening and extinction by the interstellar medium, eta Carinae is unlikely to trigger such effects to any significant degree.

Wang and White (2007) have discussed some problems with N-body simulation methods. These problems are a special case of a more general problem which has been largely unacknowledged for approximately 25 years, and affects results of all dark matter simulations on small scales (the definition of 'small' varying with time). Extensive numerical experiments with multiple types of N-body codes have demonstrated that spurious fluctuations due to particle discreteness grow rapidly even in the presence of substantial small-scale power from the intended model spectrum, and modify the results on scales smaller than the mean comoving interparticle separation. This implies that the spatial resolution of such simulations is typically limited not by the force softening length, often referred to as the 'resolution', and not by the particle density in halos. Instead it is approximately N^-1/3, where N is the mean particle density, and of course depends on and improves very slowly with increased number of particles. This calls into question many results on smaller scales over more than two decades.

Medvedev and Melott (2007) have suggested that periodicity in fossil biodiversity may be induced by cosmic rays which vary as the Solar System oscillates normal to the galactic disk. We re-examine the evidence for a 62 million year (Myr) periodicity in biodiversity throughout the Phanerozoic history of animal life reported by Rohde & Mueller (2005), as well as related questions of periodicity in origination and extinction. We find that the signal is robust against variations in methods of analysis, and is based on fluctuations in the Paleozoic and a substantial part of the Mesozoic. Examination of origination and extinction is somewhat ambiguous, with results depending upon procedure. Origination and extinction intensity as defined by RM may be affected by an artifact at 27 Myr in the duration of stratigraphic intervals. Nevertheless, when a procedure free of this artifact is implemented, the 27 Myr periodicity appears in origination, suggesting that the artifact may ultimately be based on a signal in the data. A 62 Myr feature appears in extinction, when this same procedure is used. We conclude that evidence for a periodicity at 62 Myr is robust, and evidence for periodicity at approximately 27 Myr is also present, albeit more ambiguous.

Recent work has revealed a 62 (+/-) 3-million-year cycle in the fossil diversity in the past 542 My, however no plausible mechanism has been found. We propose that the cycle may be caused by modulation of cosmic ray (CR) flux by the Solar system vertical oscillation (64 My period) in the galaxy, the galactic north-south anisotropy of CR production in the galactic halo/wind/termination shock (due to the galactic motion toward the Virgo cluster), and the shielding by galactic magnetic fields. We revisit the mechanism of CR propagation and show that CR flux can vary by a factor of about 4.6 and reach a maximum at north-most displacement of the Sun. The very high statistical significance of (i) the phase agreement between Solar north-ward excursions and the diversity minima and (ii) the correlation of the magnitude of diversity drops with CR amplitudes through all cycles provide solid support for our model. Various observational predictions which can be used to confirm or falsify our hypothesis are presented.

We have modeled atmospheric effects, especially ozone depletion, due to a solar proton event which probably accompanied the extreme magnetic storm of 1-2 September 1859. We use an inferred proton fluence for this event as estimated from nitrate levels in Greenland ice cores. We present results showing production of odd nitrogen compounds and their impact on ozone. We also compute rainout of nitrate in our model and compare to values from ice core data.

Gamma-Ray Bursts (GRBs) directed at Earth from within a few kpc may have damaged the biosphere, primarily though changes in atmospheric chemistry which admit greatly increased Solar UV. However, GRBs are highly variable in spectrum and duration. Recent observations indicate that short (~0.1 s) burst GRBs, which have harder spectra, may be sufficiently abundant at low redshift that they may offer an additional significant effect. A much longer timescale is associated with shock breakout luminosity observed in the soft X-ray (~10^3 s) and UV (~10^5 s) emission, and radioactive decay gamma-ray line radiation emitted during the light curve phase of supernovae (~10^7 s). Here we generalize our atmospheric computations to include a broad range of peak photon energies and investigate the effect of burst duration while holding total fluence and other parameters constant. The results can be used to estimate the probable impact of various kinds of ionizing events (such as short GRBs, X-ray flashes, supernovae) upon the terrestrial atmosphere. We find that the ultimate intensity of atmospheric effects varies only slightly with burst duration from 10^-1 s to 10^8 s. Therefore, the effect of many astrophysical events causing atmospheric ionization can be approximated without including time development. Detailed modeling requires specification of the season and latitude of the event. Harder photon spectra produce greater atmospheric effects for spectra with peaks up to about 20 MeV, because of greater penetration into the stratosphere.

We have examined the Fourier power spectrum as well as the Hurst exponent of extinction, origination, and total biodiversity in the marine fossil record, using a recently improved geologic timescale. We find all of them strongly inconsistent with past claims of self-similarity as well as inconsistent with random walk behavior. Instead, they are dominated by low-frequency power, with approximate f^-2 power over one decade in frequency. The spectrum turns over at about 10^8 y, lending plausibility to connections with galactic dynamics. Even in the background of this low-frequency dominance, a previously noted 62 My biodiversity cycle stands out with better than 99% confidence above the noise level, accounting for about 35% of the total variance in the fossil biodiversity record.

We describe results of modeling the effects on Earth-like planets of long-duration gamma-ray bursts (GRBs) within a few kiloparsecs. A primary effect is generation of nitrogen oxide compounds which deplete ozone. Ozone depletion leads to an increase in solar UVB radiation at the surface, enhancing DNA damage, particularly in marine microorganisms such as phytoplankton. In addition, we expect increased atmospheric opacity due to buildup of nitrogen dioxide produced by the burst and enhanced precipitation of nitric acid. We review here previous work on this subject and discuss recent developments, including further discussion of our estimates of the rates of impacting GRBs and the possible role of short-duration bursts.

Stanek et al. (astro-ph/0604113) have noted that the four low-redshift long-duration gamma-ray bursts (LSB) observed to date all occurred in faint, metal-poor galaxies. Given this selection, they argue that it is improbable that there has been a substantial population of Milky Way galaxy bursts sufficiently recently to affect life on Earth. This argument ignores the heterogeneity of stellar populations in the Milky Way, with evidence for continuing mergers with low-metallicity dwarf galaxies; observational analysis that points to LSBs being hosted by such galaxies undergoing interaction; and the existence of a likely recent GRB remnant in our galaxy.

Gamma-Ray Bursts (GRBs) are likely to have made a number of significant impacts on the Earth during the last billion years. We have used a two-dimensional atmospheric model to investigate the effects on the Earth's atmosphere of GRBs delivering a range of fluences, at various latitudes, at the equinoxes and solstices, and at different times of day. We have estimated DNA damage levels caused by increased solar UVB radiation, reduction in solar visible light due to $\mathrm{NO_2}$ opacity, and deposition of nitrates through rainout of $\mathrm{HNO_3}$. For the ``typical'' nearest burst in the last billion years, we find globally averaged ozone depletion up to 38%. Localized depletion reaches as much as 74%. Significant global depletion (at least 10%) persists up to about 7 years after the burst. Our results depend strongly on time of year and latitude over which the burst occurs. We find DNA damage of up to 16 times the normal annual global average, well above lethal levels for simple life forms such as phytoplankton. The greatest damage occurs at low to mid latitudes. We find reductions in visible sunlight of a few percent, primarily in the polar regions. Nitrate deposition similar to or slightly greater than that currently caused by lightning is also observed, lasting several years. We discuss how these results support the hypothesis that the Late Ordovician mass extinction may have been initiated by a GRB.

It is likely that one or more gamma-ray bursts within our galaxy have strongly irradiated the Earth in the last Gy. This produces significant atmospheric ionization and dissociation, resulting in ozone depletion and DNA-damaging ultraviolet solar flux reaching the surface for up to a decade. Here we show the first detailed computation of two other significant effects. Visible opacity of NO2 is sufficient to reduce solar energy at the surface up to a few percent, with the greatest effect at the poles, which may be sufficient to initiate glaciation. Rainout of dilute nitric acid is could have been important for a burst nearer than our conservative nearest burst. These results support the hypothesis that the characteristics of the late Ordovician mass extinction are consistent with GRB initiation.

Based on cosmological rates, it is probable that at least once in the last Gy the Earth has been irradiated by a gamma-ray burst in our Galaxy from within 2 kpc. Using a two-dimensional atmospheric model we have performed the first computation of the effects upon the Earth's atmosphere of one such impulsive event. A ten second burst delivering 100 kJ/m^2 to the Earth penetrates to the stratosphere and results in globally averaged ozone depletion of 35%, with depletion reaching 55% at some latitudes. Significant global depletion persists for over 5 years after the burst. This depletion would have dramatic implications for life since a 50% decrease in ozone column density results in approximately three times the normal UVB flux. Widespread extinctions are likely, based on extrapolation from UVB sensitivity of modern organisms. Additional effects include a shot of nitrate fertilizer and NO2 opacity in the visible providing a cooling perturbation to the climate over a similar timescale. These results lend support to the hypothesis that a GRB may have initiated the late Ordovician mass extinction (Melott et al. 2004).

We compute the ozone distribution for a model terrestrial paleoatmosphere in which the present oxygen abundance is largely replaced by carbon dioxide, which we argue is a reasonable working assumption. In principle, the presence of carbon dioxide might supplement the ozone shield as compared with models based on nitrogen without high carbon dioxide abundance so that early life need not have been as UV-resistant as often assumed. An extrasolar planet with a high-CO2 atmosphere might contain enough O3 to be a source of false positive biomarkers. We find that the globally averaged O3 column density can be the same, or nearly four times higher (depending upon the O2 partial pressure) when CO2 is used in place of N2 as the replacement component for lowered O2 in a 1-atm terrestrial planet with solar radiation. The effect is important for making quantitative deductions from future data, but does not invalidate the use of O3 as a biomarker for free oxygen. These results make prospects for detection of extrasolar planetary O3 absorption somewhat better than before.

We have studied morphological evolution in clusters simulated in the adiabatic limit and with radiative cooling. Cluster morphology in the redshift range, $0 < z < 0.5$, is quantified by multiplicity and ellipticity. In terms of ellipticity, our result indicates slow evolution in cluster shapes compared to those observed in the X-ray and optical wavelengths. The result is consistent with Floor, Melott & Motl (2003). In terms of multiplicity, however, the result indicate relatively stronger evolution (compared to ellipticity but still weaker than observation) in the structure of simulated clusters suggesting that for comparative studies of simulation and observation, sub-structure measures are more sensitive than the shape measures. We highlight a few possibilities responsible for the discrepancy in the shape evolution of simulated and real clusters.

Gamma-ray bursts (hereafter GRB) produce a flux of radiation detectable across the observable Universe, and at least some of them are associated with galaxies. A GRB within our own Ggalaxy could do considerable damage to the Earth's biosphere; rate estimates suggest that a dangerously near GRB should occur on average two or more times per billion years. At least five times in the history of life, the Earth experienced mass extinctions that eliminated a large percentage of the biota. Many possible causes have been documented, and GRB may also have contributed. The late Ordovician mass extinction approximately 440 million years ago may be at least partly the result of a GRB. A special feature of GRB in terms of terrestrial effects is a nearly impulsive energy input of order 10 s. Due to expected severe depletion of the ozone layer, intense solar ultraviolet radiation would result from a nearby GRB, and some of the patterns of extinction and survivorship at this time may be attributable to elevated levels of UV radiation reaching the Earth. In addition a GRB could trigger the global cooling which occurs at the end of the Ordovician period that follows an interval of relatively warm climate. Intense rapid cooling and glaciation at that time, previously identified as the probable cause of this mass extinction, may have resulted from a GRB.

The rate of galaxy cluster eccentricity evolution is useful in understanding large scale structure. Rapid evolution for $z < $ 0.13 has been found in two different observed cluster samples. We present an analysis of projections of 41 clusters produced in hydrodynamic simulations augmented with radiative cooling and 43 clusters from adiabatic simulations. This new, larger set of simulated clusters strengthens the claims of previous eccentricity studies. We find very slow evolution in simulated clusters, significantly different from the reported rates of observational eccentricity evolution. We estimate the rate of change of eccentricity with redshift and compare the rates between simulated and observed clusters. We also use a variable aperture radius to compute the eccentricity, r$_{200}$. This method is much more robust than the fixed aperture radius used in previous studies. Apparently radiative cooling does not change cluster morphology on scales large enough to alter eccentricity. The discrepancy between simulated and observed cluster eccentricity remains. Observational bias or incomplete physics in simulations must be present to produce halos that evolve so differently.

We report on a new study of the velocity distribution in N-body simulations. We investigate the center-of-mass and internal kinetic energies of coarsening cells as a function of time, cell size and cell mass. By using self-similar cosmological models, we are able to derive theoretical predictions for comparison and to assess the influence of finite-size and resolution effects. The most interesting result is the discovery of a polytropic-like relationship between the average velocity dispersion (internal kinetic energy) and the mass density in an intermediate range of densities, K ~ rho^{2-eta}. The exponent eta measures the deviations from the virial prediction, eta_virial=0. For self-similar models, eta depends only on the spectral index of the initial power spectrum. We also study CDM models and confirm a previous result that the same polytropic-like dependence exists (astro-ph/0103313), with a time and coarsening length dependent eta. The dependence K(rho) is an important input for a recently proposed theoretical model of cosmological structure formation which improves over the standard dust model (pressureless fluid) by regularizing the density singularities.

We describe a new method of overcoming problems inherent in peculiar velocity surveys by using data compression as a filter with which to separate large-scale, linear flows from small-scale noise that biases the results systematically. We demonstrate the effectiveness of our method using realistic catalogs of galaxy velocities drawn from N--body simulations. Our tests show that a likelihood analysis of simulated catalogs that uses all of the information contained in the peculiar velocities results in a bias in the estimation of the power spectrum shape parameter $\Gamma$ and amplitude $\beta$, and that our method of analysis effectively removes this bias. We expect that this new method will cause peculiar velocity surveys to re--emerge as a useful tool to determine cosmological parameters.

We have used N-body simulations to develop two independent methods to quantify redshift distortions known as the Bull's Eye effect (large scale infall plus small scale virial motion). This effect depends upon the mass density, $\Omega_0$, so measuring it can in principle give an estimate of this important cosmological parameter. We are able to measure the effect and distinguish between its strength for high and low values of $\Omega_0$. Unlike other techniques which utilize redshift distortions, one of our methods is relatively insensitive to bias. In one approach, we use path lengths between contour crossings of the density field. The other is based upon percolation. We have found both methods to be successful in quantifying the effect and distinguishing between values of $\Omega_0$. However, only the path lengths method exhibits low sensitivity to bias.

The large-scale structure of the Universe is thought to evolve by a process of gravitational amplification from low-amplitude Gaussian noise generated in the early Universe. The later, non-linear stages of gravitation-induced clustering produce phase correlations with well-defined statistical properties. In particular, the distribution of phase differences $D$ between neighboring Fourier modes provides useful insights into the clustering phenomenon. Here we develop an approximate theory for the probability distribution $D$ and test it using a large battery of numerical simulations. We find a remarkable universal form for the distribution which is well described by theoretical arguments.

Analyses of peculiar velocity surveys face several challenges, including low signal--to--noise in individual velocity measurements and the presence of small--scale, nonlinear flows. This is the second in a series of papers in which we describe a new method of overcoming these problems by using data compression as a filter with which to separate large--scale, linear flows from small--scale noise that can bias results. We demonstrate the effectiveness of our method using realistic catalogs of galaxy velocities drawn from N--body simulations. Our tests show that a likelihood analysis of simulated catalogs that uses all of the information contained in the peculiar velocities results in a bias in the estimation of the power spectrum shape parameter $\Gamma$ and amplitude $\beta$, and that our method of analysis effectively removes this bias. We expect that this new method will cause peculiar velocity surveys to re--emerge as a useful tool to determine cosmological parameters.

Strong cluster eccentricity evolution for $z \le 0.13$ has appeared in a variety of observational data sets. We examine the evolution of eccentricity in simulated galaxy clusters using a variety of simulation methodologies, amplitude normalizations, and background cosmologies. We do not find find such evolution for $z < 0.1$ in any of our simulation ensembles. We suggest a systematic error in the form of a redshift-dependent selection effect in cluster catalogs or missing physics in cluster simulations important enough to modify the cluster morphology.

We have calculated lower limits for the global effective thermal conductivity in a sample of 165 Abell clusters. We assumed that cluster X-ray luminosity is compensated by a conductive heat flux which we then compare with an upper limit to the temperature gradient inferred from the cluster temperature and radius. This gives a lower limit to the thermal conductivity and therefore to the relative suppression from the Spitzer conductivity, kappa_{Sp}. Not a single cluster requires super-Spitzer values of thermal conduction to balance the observed X-ray luminosity. The suppression coefficient f=kappa/kappa_{Sp} is clustered in a range 10^{-2} < f < 0.4. A weak dependence of $f$ versus $z$ is observed over $ 0 < z <0.41$. Possible biases and/or selection effects are discussed.

We examine the orientations of rich galaxy cluster X-ray isophotes with respect to their rich nearest neighbors using existing samples of Abell cluster position angles measured from {\it Einstein} and {\it ROSAT} observations. We study a merged subset of these samples using updated and improved positions and redshifts for Abell/ACO clusters. We find high confidence for alignment, which increases as nearest neighbor distance is restricted. We conclude that there is a strong alignment signal in all this data, consistent with gravitational instability acting on Gaussian perturbations.

We present a new method for the analysis of peculiar velocity surveys which removes contributions to velocities from small scale, nonlinear velocity modes while retaining information about large scale motions. Our method utilizes Karhunen--Lo\`eve methods of data compression to construct a set of moments out of the velocities which are minimally sensitive to small scale power. The set of moments are then used in a likelihood analysis. We develop criteria for the selection of moments, as well as a statistic to quantify the overall sensitivity of a set of moments to small scale power. Although we discuss our method in the context of peculiar velocity surveys, it may also prove useful in other situations where data filtering is required.

We argue that there has been substantial evolution in the ellipticity of rich galaxy clusters between 0 < z < 0.1, and suggest that this is additional evidence for a low matter-density Universe.

We have examined whether the major axes of rich galaxy clusters tend to point toward their nearest neighboring cluster. We have used the data of Ulmer, McMillan, and Kowalski, who used position angles based on X-ray morphology. We also studied a subset of this sample with updated positions and distances from the MX Northern Abell Cluster Survey (for rich clusters ($R \geq 1$) with well known redshifts). A Kolmogorov-Smirnov (KS) test showed no significant signal for nonrandom angles on any scale $\leq 100h^{-1}$Mpc. However, refining the null hypothesis with the Wilcoxon rank-sum test, we found a high confidence signal for alignment. Confidence levels increase to a high of 99.997% as only near neighbors which are very close are considered. We conclude there is a strong alignment signal in the data, consistent with gravitational instability acting on Gaussian perturbations.

The sizes and shapes of voids in a galaxy survey depend not only on the physics of structure formation, but also on the sampling density of the survey and on the algorithm used to define voids. Using an N-body simulation with a CDM power spectrum, we study the properties of voids in samples with different number densities of galaxies, both in redshift space and in real space. When voids are defined as regions totally empty of galaxies, their characteristic volume is strongly dependent on sampling density; when they are defined as regions whose density is 0.2 times the mean galaxy density, the dependence is less strong. We compare two void-finding algorithms, one in which voids are nonoverlapping spheres, and one, based on the algorithm of Aikio and Mahonen, which does not predefine the shape of a void. Regardless of the algorithm chosen, the characteristic void size is larger in redshift space than in real space, and is larger for low sampling densities than for high sampling densities. We define an elongation statistic Q which measures the tendency of voids to be stretched or squashed along the line of sight. Using this statistic, we find that at sufficiently high sampling densities (comparable to the number densities of galaxies brighter than L_*), large voids tend to be slightly elongated along the line of sight in redshift space.

We examine the relationship between the total X-ray and optical luminosities of groups and clusters of galaxies taken from various samples in the literature. The clusters and groups were drawn from four different catalogs: (1) the Abell/ACO catalog, (2) the Edinburgh-Durham Cluster Catalog, (3) the RASS Bright Cluster Sample, and (4) galaxy groups selected from the CfA redshift survey. These catalogs represent a significant cross-section of cluster selection tchniques as well as a wide range in mass scale and can be considered statistically independent. We have calculated new cluster X-ray and optical luminosities if not available in the literature. Based on 126 systems, our analysis shows that the total optical luminosity of a cluster is directly proportional to the total X-ray luminosity over a wide mass range and across all four cluster samples studied herein. We also show that total cluster optical luminosity is a good indicator of cluster virial mass, but richness is not. Our results suggest that (1) the selection method of galaxy clusters may not be crucial, providing wider latitude in assembling catalogs, (2) the luminosity per baryon does not vary systematically between systems ranging from groups to very rich clusters. We propose that future optical catalogs of clusters use the total optical luminosity of a cluster, instead of galaxy richness, as luminosity appears to be a better measure of cluster mass and can be directly related to X-ray catalogs of clusters. This will also facilitate the first direct comparison between the growing number of X-ray and optical catalogs.

Galaxy clusters approximate a planar (FP) distribution in a three-dimensional parameter space which can be characterized by optical luminosity, half-light radius, and X-ray luminosity. Using a high-quality catalog of cluster redshifts, we find the nearest neighbor cluster for those common to an FP study and the cluster catalog. Examining scatter about the FP, we find 99.2% confidence that it is dependent on nearest neighbor distance. Our study of X-Ray clusters finds that those with high central gas densities are systematically closer to neighbor clusters. If we combine results here with those of Fritsch and Buchert, we find an explanation for some of our previous conclusions: Clusters in close proximity to other clusters are more likely to have massive cooling flows because they are more relaxed and have higher central gas densities.

We compare the statistical properties of structures normal and transverse to the line of sight which appear in observational data from redshift surveys. We present a statistic which can quantify this effect in a conceptually different way from standard analyses of distortions of the power-spectrum or correlation function. From tests with N-body experiments, we argue that this statistic represents a new, more direct and potentially powerful diagnostic of the cosmological density parameter Omega.

We apply Minkowski functionals and various derived measures to decipher the morphological properties of large-scale structure seen in simulations of gravitational evolution. Minkowski functionals of isodensity contours serve as tools to test global properties of the density field. Furthermore, we identify coherent objects at various threshold levels and calculate their partial Minkowski functionals. We propose a set of two derived dimensionless quantities, planarity and filamentarity, which reduce the morphological information in a simple and intuitive way. Several simulations of the gravitational evolution of initial power-law spectra provide a framework for systematic tests of our method.

It is now well established that galaxies are biased tracers of the distribution of matter, although it is still not known what form this bias takes. In local bias models the propensity for a galaxy to form at a point depends only on the overall density of matter at that point. Hierarchical scaling arguments allow one to build a fully-specified model of the underlying distribution of matter and to explore the effects of local bias in the regime of strong clustering. Using a generating-function method developed by Bernardeau & Schaeffer (1992), we show that hierarchical models lead one directly to the conclusion that a local bias does not alter the shape of the galaxy correlation function relative to the matter correlation function on large scales. This provides an elegant extension of a result first obtained by Coles (1993) for Gaussian underlying fields and confirms the conclusions of Scherrer & Weinberg (1998) obtained using a different approach. We also argue that particularly dense regions in a hierarchical density field display a form of bias that is different from that obtained by selecting such peaks in Gaussian fields: they are themselves hierarchically distributed with scaling parameters $S_p=p^{(p-2)}$. This kind of bias is also factorizable, thus in principle furnishing a simple test of this class of models.

With the availability of large-scale redshift survey data, it is now becoming possible to explore correlations between large-scale structure and the properties and morphologies of galaxy clusters. We investigate the spatial distributions of a 98% complete, volume-limited sample of nearby (z<0.1) Abell clusters with well-determined redshifts and find that cooling flow clusters with high mass-accretion rates have nearest neighbors which are much closer than those of other clusters in the sample (at the 99.8% confidence level), and reside in more crowded environments out to 30/h Mpc. Several possible explanations of this effect are discussed.

We develop a diagrammatic technique to represent the multi-point cumulative probability density function (CPDF) of mass fluctuations in terms of the statistical properties of individual collapsed objects and relate this to other statistical descriptors such as cumulants, cumulant correlators and factorial moments. We use this approach to establish key scaling relations describing various measurable statistical quantities if clustering follows a simple general scaling ansatz, as expected in hierarchical models. We test these detailed predictions against high-resolution numerical simulations. We show that, when appropriate variables are used, the count probability distribution function (CPDF) and void probability distribution function (VPF) shows clear scaling properties in the non-linear regime. Generalising the results to the two-point count probability distribution function (2CPDF), and the bivariate void probability function (2VPF) we find good match with numerical simulations. We explore the behaviour of the bias associated with collapsed objects in limit of large separations, finding that it depends only on the intrinsic scaling parameter associated with collapsed objects and that the bias for two different objects can be expressed as a product of the individual biases of the objects. Finally we use hierarchal ansatz to develop a novel technique for correcting finite-volume effects in the estimation of multi-point statistical quantities from observational data (abridged).

Within Abell galaxy clusters containing wide-angle tailed radio sources, there is evidence of a ``prevailing wind'' which directs the WAT jets. We study the alignment of WAT jets and nearby clusters to test the idea that this wind may be a fossil of drainage along large-scale supercluster axes. We also test this idea with a study of the alignment of WAT jets and supercluster axes. Statistical test neighbours indicate no alignment of WAT jets towards nearest clusters, but do indicate approximately 98% confidence in alignment with the long axis of the supercluster in which the cluster lies. We find a preferred scale for such superclusters of order 25 Mpc $h^{-1}$.

The highly nonlinear regime of gravitational clustering is characterized by the presence of scale-invariant form of many-body correlation functions. Useful insights can be obtained by investigating the consequences of a generic scaling {\em ansatz}. Extending earlier studies by Bernardeau & Schaeffer (1992) we calculate the detailed consequences of such scaling. We generalise the two-point cumulant correlators to a hierarchy of multi-point cumulant correlators (MCC). We introduce the concept of reduced cumulant correlators (RCC) and their related generating functions. Using analytical and diagrammatic method we show that every new vertex of the tree representation of higher-order correlations has its own reduced cumulant correlator. In the limit of large separations, MCCs of arbitrary order can be expressed in terms of RCCs of the same and lower order. We relate the generating functions of RCCs with the statistics of {\em collapsed} objects. We develop the hierarchy for the correlation functions of overdense regions. In this vein, we compute the lower-order cumulants and cumulant correlators for overdense regions. Our study shows how they vary as a function of the initial power spectrum of primordial density fluctuations (abridged).

The cumulant correlators, $C_{pq}$, are statistical quantities that generalise the better-known $S_p$ parameters; the former are obtained from the two-point probability distribution function of the density fluctuations while the latter describe only the one-point distribution. If galaxy clustering develops from Gaussian initial fluctuations and a small-angle approximation is adopted, standard perturbative methods suggest a particular hierarchical relationship of the $C_{pq}$ for projected clustering data, such as the APM survey. We establish the usefulness of the two-point cumulants for describing hierarchical clustering by comparing such calculations against available measurements from projected catalogs, finding very good agreement. We extend the idea of cumulant correlators to multi-point {\em generalised} cumulant correlators (related to the higher-order correlation functions). We extend previous studies in the highly nonlinear regime to express the generalised cumulant correlators in terms of the underlying ``tree amplitudes'' of hierarchical scaling models. Such considerations lead to a technique for determining these hierarchical amplitudes, to arbitrary order, from galaxy catalogs and numerical simulations. Knowledge of these amplitudes yields important clues about the phenomenology of gravitational clustering. For instance, we show that three-point cumulant correlator can be used to separate the tree amplitudes up to sixth order. We also combine the particular hierarchical {\em ansatz} of Bernardeau & Schaeffer (1992) with extended and hyperextended perturbation theory to infer values of the tree amplitudes in the highly nonlinear regime.

We compare the statistical properties of structures normal and transverse to the line of sight which appear in theoretical N-body simulations of structure formation, and seem also to be present in observational data from redshift surveys. We present a statistic which can quantify this effect in a conceptually different way from standard analyses of distortions of the power-spectrum or correlation function. From tests with $N$--body experiments, we argue that this statistic represents a new and potentially powerful diagnostic of the cosmological density parameter, $\Omega_0$.

Shape dependence of higher order correlations introduces complication in direct determination of these quantities. For this reason theoretical and observational progress has been restricted in calculating one point distribution functions and related moments. Methods based on factorial moments of two point count probability distribution (cumulant correlators) were recently shown to be efficient in subtracting discreteness effects and extracting useful information from galaxy catalogs. We use these cumulant correlators $Q_{NM}$ to study clustering in scale free simulations both in two and three dimensions. Using this method in the highly nonlinear regime we were able to separate hierarchal amplitudes $Q_3, R_a and R_b$ associated with different tree graphs contributing to third and fourth order correlation functions. They were found to increase with power on large scales. Results based on factorial moments of one point cumulants show very good agreement with that on cumulant correlators. Comparison of simulation results with perturbation theory and extended perturbation theory were found to be in reasonable agreement. Comparisons were also made against predictions from hierarchal models for higher order correlations. We argue that finite volume corrections are very important for computation of cumulant correlators.

We report on a series of tests of agreement between three types of N-body simulations: PM, P$^3$M, and Tree codes. We find good agreement in both the individual and the statistical properties only on scales larger than the mean interparticle separation. As a result, we question most numerical results at and below below galaxy scales, either concerning primordial dark matter or baryonic matter coupled to it by gravitation. Mass resolution (specifically, sufficient to have the mean interparticle separation comparable to the smoothing scales) appears be necessary for results to converge.

We explore some of the effects that discreteness and two-body scattering may have on N-body simulations with ``realistic'' cosmological initial conditions. We use an identical subset of particles from the initial conditions for a $128^3$ Particle-Mesh (PM) calculation as the initial conditions for a variety P$^3$M and Tree code runs. We investigate the effect of mass resolution (the mean interparticle separation) since most ``high resolution'' codes only have high resolution in gravitational force. The phase-insensitive two--point statistics, such as the power spectrum (autocorrelation) are somewhat affected by these variations, but phase-sensitive statistics show greater differences. Results converge at the mean interparticle separation scale of the lowest mass-resolution code. As more particles are added, but the force resolution is held constant, the P$^3$M and the Tree runs agree more and more strongly with each other and with the PM run which had the same initial conditions. This shows high particle density is necessary for correct time evolution, since many different results cannot all be correct. However, they do not so converge to a PM run which continued the fluctuations to small scales. Our results show that ignoring them is a major source of error on comoving scales of the missing wavelengths. This can be resolved by putting in a high particle density. Since the codes never agree well on scales below the mean comoving interparticle separation, we find little justification for quantitative predictions on this scale. Some measures vary by 50%, but others can be off by a factor of three or more. Our results suggest possible problems with the density of galaxy halos, formation of early generation objects such as QSO absorber clouds, etc.

The bispectrum, the three-point function of density fluctuations in Fourier space, is the lowest order statistic that carries information about the spatial coherence of large-scale structures. For Gaussian initial conditions, when the density fluctuation amplitude is small, tree-level (leading order) perturbation theory predicts a characteristic dependence of the bispectrum on the shape of the triangle formed by the three wave vectors. This configuration dependence provides a signature of gravitational instability, and departures from it in galaxy catalogs can be interpreted as due to bias, that is, non-gravitational effects. On the other hand, N-body simulations indicate that the reduced three-point function becomes relatively shape-independent in the strongly non-linear regime. In order to understand this non-linear transition and assess the domain of reliability of shape-dependence as a probe of bias, we calculate the one-loop (next-to-leading order) corrections to the bispectrum in perturbation theory. We compare these results with measurements in numerical simulations with scale-free and CDM initial power spectra. We find that the one-loop corrections account very well for the departures from the tree-level results measured in numerical simulations on weakly non-linear scales. In this regime, the reduced bispectrum qualitatively retains its tree-level shape, but the amplitude can change significantly. At smaller scales, the reduced bispectrum in the simulations starts to flatten, an effect which can be partially understood from the one-loop results. In the strong clustering regime, where perturbation theory breaks down entirely, the simulation results confirm that the reduced bispectrum has almost no dependence on triangle shape, in rough agreement with the hierarchical ansatz.

Using an ensemble of high resolution 2D numerical simulations, we explore the scaling properties of cosmological density fluctuations in the non-linear regime. We study the scaling behaviour of the usual $N$--point volume-averaged correlations, and also examine the scaling of the entire probability density function (pdf) of the fluctuations. We focus on two important issues: (i) whether the scaling behaviour of 2D clustering is consistent with what one infer from radial collapse arguments; and (ii) whether there is any evidence from these high-resolution simulations that a regime of stable clustering is ever entered. We find that the answers are (i) yes and (ii) no. We further find that the behaviour of the highly non-linear regime in these simulations suggests the existence of a regime where the correlation function is independent of the initial power spectrum.

Perturbation Theory (PT) applied to a cosmological density field with Gaussian initial fluctuations suggests a specific hierarchy for the correlation functions when the variance is small. In particular quantitative predictions have been made for the moments and the shape of the one-point probability distribution function (PDF) of the top-hat smoothed density. In this paper we perform a series of systematic checks of these predictions against N-body computations both in 2D and 3D with a wide range of featureless power spectra. In agreement with previous studies, we found that the reconstructed PDF-s work remarkably well down to very low probabilities, even when the variance approaches unity. Our results for 2D reproduce the features for the 3D dynamics. In particular we found that the PT predictions are more accurate for spectra with less power on small scales. The nonlinear regime has been explored with various tools, PDF-s, moments and Void Probability Function (VPF). These studies have been done with unprecedented dynamical range, especially for the 2D case, allowing in particular more robust determinations of the asymptotic behaviour of the VPF. We have also introduced a new method to determine the moments based on the factorial moments. Results using this method and taking into account the finite volume effects are presented.

The Lagrangian perturbation theory on Friedmann-Lemaitre cosmologies is compared with numerical simulations (tree-, adaptive P$^3$M- and PM codes). In previous work we have probed the large-scale performance of the Lagrangian perturbation solutions up to the third order by studying their cross- correlations with N-body simulations for various power spectra (Buchert etal 1994, Melott etal 1995, Weiss etal 1996). Thereby, spatial optimization techniques were applied by (high-frequency-)filtering of the initial power spectra. In this work the novel method of temporal optimization [Shifted-Time- Approximation (STA) and Frozen-Time-Approximation (FTA)] is investigated and used. The method is designed to compensate the native property of Lagrangian perturbation solutions to delay the collapse of structures. The method can be treated analytically. Applying the STA and FTA prescriptions a significant improvement of the performance of Lagrangian perturbation schemes (as measured by cross-correlation statistics) is observed. Using this tool we investigate a local study of special clustering models of dark matter as candidates for typical elements of the large-scale structure in the Universe, and so also focus on the performance of the perturbation solutions on smaller scales at high-spatial resolution. The models analyzed were presented in (Buchert etal 1996) and allow studying typical features of the clustering process in the weakly non-linear regime. The spatial and temporal limits of applicability of the solutions at second and third order are determined and compared with the first-order solution, which is equivalent to the ``Zel'dovich approximation'' (Zel'dovich 1970, 1973) for the type of initial data analyzed.

Two-body scattering and other discreteness effects are unimportant in cosmological gravitational clustering in most scenarios, since the dark matter has a small particle mass. The collective field should determine evolution: Two-body scattering in simulations violates the Poisson-Vlasov equations. We test this in PM, P$^3$M, Tree, and NGPM codes, noting that a collisionless code will preserve the one-dimensional character of plane wave collapse. We find collisionality vanishing as the softening parameter approaches the mean interparticle separation. Solutions for the problem are suggested, involving greater computer power, PM-based nested grid codes, and a more conservative approach to resolution claims.

We investigate a distortion in redshift-space which causes galaxies to appear to lie in walls concentric about the observer, forming a rough bull's-eye pattern. We simulate what an observer would see in a thin slice of redshift-space, including a magnitude limit and constant slice angle. The result is an enhanced ring of galaxies encircling the observer at a distance roughly corresponding to the peak of the selection function. This ring is an artificial enhancement of weak features in real-space. This may explain visually prominent features such as the ``Great Wall'' and periodicity found in deep narrow fields.

We test tree-level perturbation theory for Gaussian initial conditions with power spectra $P(k)\propto k^n$ by comparing the probability distribution function (PDF) for the density predicted by the Local Lagrangian Approximation (LLA) with the results of numerical gravitational clustering simulations. Our results indicate that our approximation correctly reproduces the evolved density PDF for $-3 \leq n \leq-1$ power spectra up to the weakly nonlinear regime, while it shows marginal agreement for power indices n=0 and +1 in the linear regime and poor agreement beyond this point. This suggests that tree-level perturbation theory (as realized in the Local Lagrangian Approximation) can accurately predict the density distribution function for $-3 \leq n \leq -1$ but fails for $n \ge 0$.

In this study we use simulations of 128$^3$ particles to study the ellipticity and orientation of clusters of galaxies in N-body simulations of differing power-law initial spectra (P(k) \propto k^n ,n = +1, 0, -1, -2$), and density parameters ($\Omega_0 = 0.2$ to 1.0). Furthermore, unlike most theoretical studies we mimic most observers by removing all particles which lie at distances greater than 2 1/h Mpc from the cluster center of mass. We computed the axial ratio and the principal axes using the inertia tensor of each cluster. The mean ellipticity of clusters increases strongly with increasing $n$. We also find that clusters tend to become more spherical at smaller radii. We compared the orientation of a cluster to the orientation of neighboring clusters as a function of distance (correlation). In addition, we considered whether a cluster's major axis tends to lie along the line connecting it to a neighboring cluster, as a function of distance (alignment). Both alignments and correlations were computed in three dimensions and in projection to mimic observational surveys. Our results show that significant alignments exist for all spectra at small separations ($D < 15 h^{-1}$ Mpc) but drops off at larger distance in a strongly $n-$dependent way.

We present results of tests of the evolution of small ``fluid elements'' in cosmological N--body simulations, to examine the validity of their treatment as particles. We find that even very small elements typically collapse along one axis while expanding along another, often to twice or more their initial comoving diameter. This represents a possible problem for high--resolution uses of such simulations.

Using two-dimensional numerical simulations of gravitational clustering, with initial power spectra of a power-law form with index n, we compare the properties of voids in real space to their properties in redshift space. Both the void probability function (VPF) and the underdense probability function (UPF) are enhanced in redshift space. The enhancement is greatest in the limit n=-2. The VPF and UPF treat voids statistically and assume that two-dimensional voids are circular; we present an algorithm which detects individual voids. Voids found by this algorithm are ellipses whose enclosed density of galaxies falls below a threshold density. When voids are identified using this algorithm, the mean void size and the maximum void size both increase in going from real space to redshift space. The increase is greatest in the limit n=-2. In redshift space, the principal axes of the largest voids in the n=2 and n=0 simulations show a statistically significant tendency (at the 95% confidence level) to be distributed anisotropically, relative to the line of sight from the origin to the void center.

The potential is a constant to linear order in cosmological gravitational clustering. In this Letter we present results of testing the conjecture, proposed by Pauls and Melott (1995), that the effect of nonlinear evolution on the potential can be better described by smoothing it on the scale of nonlinearity. We show two-dimensional simulations consisting of an N-body code in which particles are accelerated not by their mutual attraction, but by the gradient of the initial potential smoothed on the current scale of nonlinearity. This approximation produces results considerably improved over using a constant potential to move particles, and it is generally better than any approximation we have tested, lending further support to the smoothing description of the evolved potential.

Since the early seventies, there was a conjecture that the first collapse of a selfgravitating dust--like medium (appropriate approximation for nonbaryonic dark matter) results in the formation of a ``pancake" object, that is a thin surface. The conjecture has been based on the Zel'dovich approximate solution of the nonlinear gravitational instability of a generic smooth density perturbation. Recent works cast doubt on the Zel'dovich conjecture, suggesting that the first collapse might be point--like or filament--like rather than pancake--like. Our $N$--body simulations show first pancake collapse. We can reject with 97\% confidence the Bayesian prior that the other kinds of collapse are more or equally probable.

Approximations to the exact solutions for gravitational instability in the expanding Universe are extremely useful for understanding the evolution of large--scale structure. We report on a series of tests of Newtonian Lagrangian perturbation schemes using N--body simulations for various power--spectra with scale--independent indices in the range $-3$ to $+1$. The models have been evolved deeply into the non--linear regime of structure formation in order to probe the dynamical and statistical performance of the Lagrangian perturbation schemes (whose first--order solution contains as a subset the celebrated ``Zel'dovich--approximation'', hereafter ZA). These tests reveal properties of the approximations at stages beyond the obvious validity of perturbation theory. Recently, another series of tests of different analytical and semi--numerical approximations for large--scale structure was conducted with the result that ZA displays the best dynamical performance in comparison with the N--body simulations, if the initial data were smoothed before evolving the model, i.e., a truncated form of ZA (TZA). We show in this contribution that the excellent performance of TZA can be further improved by going to second order in the Lagrangian perturbation approach. The truncated second--order Lagrangian scheme provides a useful improvement over TZA especially for negative power indices, which suggests it will be very useful for modelling standard scenarios such as ``Cold--'', ``Hot--'' and ``Mixed--Dark--Matter''.

To explain the rich structure of voids, clusters, sheets, and filaments apparent in the Universe, we present evidence for the convergence of the two classic approaches to gravitational clustering, the ``pancake'' and ``hierarchical'' pictures. We compare these two models by looking at agreement between individual structures -- the ``pancakes'' which are characteristic of the Zel'dovich Approximation (ZA) and also appear in hierarchical N-body simulations. We find that we can predict the orientation and position of N-body simulation objects rather well, with decreasing accuracy for increasing large-$k$ (small scale) power in the initial conditions. We examined an N-body simulation with initial power spectrum $P(k) \propto k^3$, and found that a modified version of ZA based on the smoothed initial potential worked well in this extreme hierarchical case, implying that even here very low-amplitude long waves dominate over local clumps (although we can see the beginning of the breakdown expected for $k^4$). In this case the correlation length of the initial potential is extremely small initially, but grows considerably as the simulation evolves. We show that the nonlinear gravitational potential strongly resembles the smoothed initial potential. This explains why ZA with smoothed initial conditions reproduces large-scale structure so well, and probably why our Universe has a coherent large-scale structure.

We study the development of gravitational instability in the strongly non-linear regime. For this purpose we use a number of statistical indicators such as filamentary statistics, spectrum of overdense/underdense regions and the void probability function, each of which probes a particular aspect of gravitational clustering. We use these statistical indicators to discriminate between different approximations to gravitational instability which we test against N-body simulations. The approximations which we test are, the truncated Zel'dovich approximation (TZ), the adhesion model (AM), and the frozen flow (FF) and linear potential (LP) approximations. Of these we find that FF and LP break down relatively early, soon after the non-linear length scale exceeds $R_*$ -- the mean distance between peaks of the gravitational potential. The reason for this break down is easy to understand, particles in FF are constrained to follow the streamlines of the initial velocity field. Shell crossing is absent in this case and structure gradually freezes as particles begin to collect near minima of the gravitational potential. In LP particles follow the lines of force of the primordial potential, oscillating about its minima at late times when the non-linear length scale $k_{\rm NL}^{-1}\simeq R_*$. Unlike FF and LP the adhesion model (and to some extent TZ) continues to give accurate results even at late times when $k_{\rm NL}^{-1} \ge R_*$. This is because both AM and TZ use

In recent years there has been a developing realization that the interesting large--scale structure of voids, ``pancakes", and filaments in the Universe is a consequence of the efficacy of an approximation scheme for cosmological gravitation clustering proposal by Zel'dovich in 1970. However, this scheme was only supposed to apply to smoothed initial conditions, and its efficacy for hierarchical clustering has not been justified. We show that this is explained by the fact that the gravitationally evolved potential from N--body simulations closely resembles the smoothed potential of the initial conditions. The resulting ``hierarchical pancaking" picture combines features of the former Soviet and Western theoretical pictures for galaxy and large--scale structure formation.

We report on a series of tests of Newtonian Lagrangian perturbation schemes using N--body simulations for various power--spectra with scale--independent indices in the range $-3$ to $+1$. The models have been evolved deeply into the non--linear regime of structure formation in order to probe the dynamical and statistical performance of the Lagrangian perturbation schemes (whose first--order solution contains as a subset the celebrated ``Zel'dovich--approximation'' (hereafter ZA). These tests reveal properties of the approximations at stages beyond the obvious validity of perturbation theory. Recently, another series of tests of different analytical and semi--numerical approximations for large--scale structure was conducted with the result that ZA displays the best dynamical performance in comparison with the N--body simulations, if the initial data were smoothed before evolving the model, i.e., a truncated form of ZA (TZA). We show in this Letter that the excellent performance of TZA can be further improved by going to second order in the Lagrangian perturbation approach. The truncated second--order Lagrangian scheme provides a useful improvement over TZA especially for negative power indices, which suggests it will be very useful for modelling standard scenarios such as ``Cold--'', ``Hot--'' and ``Mixed--Dark--Matter''.

We examine the effect of a threshold bias on the power spectrum and the bispectrum in an ensemble of numerical simulations (Gaussian initial perturbations with power law spectra P(k) \sim k^n, n=+1, 0, -1, -2) and compare our results with theoretical predictions. Our simulations are evolved sufficiently that on the scale where we apply the threshold the rms fluctuation has developed significantly into the nonlinear regime. Thus, predictions based on perturbation theory do not necessarily apply. Nevertheless, we find our results for the power spectrum, biased power simply amplified by a numerical factor, follow predicted trends, far beyond the regime where perturbation theory is expected to be valid. We find that the biased bispectrum continues to follow the so-called hierarchical form, with reduced three-point amplitude Q \approx 1 in the strongly nonlinear regime, independent of initial spectrum. In the quasi-linear perturbative regime the three-point amplitude depends on configuration shape, a behavior that is found to give useful information about the amount of bias without information about the unbiased matter distribution.

We present results showing an improvement of the accuracy of perturbation theory as applied to cosmological structure formation for a useful range of scales. The Lagrangian theory of gravitational instability of Friedmann--Lema\^\i tre cosmogonies investigated and solved up to the third order in the series of papers by Buchert (1989, 1992, 1993), Buchert \& Ehlers (1993), Buchert (1994), Ehlers \& Buchert (1994), is compared with numerical simulations. In this paper we study the dynamics of hierarchical models as a second step. In the first step (Buchert, Melott and Wei{\ss} 1994) we analyzed the performance of the Lagrangian schemes for pancake models, i.e., models which initially have a truncated power spectrum. We here explore whether the results found for pancake models carry over to hierarchical models which are evolved deeply into the non--linear regime............We find that for spectra with negative power--index the second--order scheme performs considerably better than TZA in terms of statistics which probe the dynamics, and slightly better in terms of low--order statistics like the power--spectrum. In cases with much small--scale power the gain from the higher--order schemes is small, but still measurable. However, in contrast to the results found for pancake models, where the higher--order schemes get worse than TZA at late non--linear stages and on small scales, we here find that the second--order model is as robust as TZA, retaining the improvement at later stages and on smaller scales. In view of these

I report on controlled comparison of gravitational approximation schemes linear/lognormal/adhesion/frozen-flow/Zel'dovich(ZA) and ZA's second--order generalization. In the last two cases we also created new versions of the approximation by truncation, i.e., by finding an optimum smoothing window (see text) for the initial conditions. The Zel'dovich approximation, with optimized initial smoothing, worked extremely well. Its second-order generalization was slightly better. The success of our best-choice was a result of the treatment of the phases of nonlinear Fourier components. The adhesion approximation produced the most accurate nonlinear power spectrum and density distribution, but its phase errors suggest mass condensations were moved somewhat incorrectly. Due to its better reproduction of the mass density distribution function and power spectrum, adhesion might be preferred for some uses. We recommend either n-body simulations or our modified versions of ZA, depending on the purpose. Modified ZA can rapidly generate large numbers of realizations of model universes with good accuracy down to galaxy group (or smaller) mass scales.

We have recently learned that the Zeldovich approximation can be successfully used for a far wider range of gravitational instability scenarios than formerly proposed; we study here how to extend this range. In previous work we studied the accuracy of several analytic approximations to gravitational clustering in the mildly nonlinear regime. We found that the ``truncated Zel'dovich approximation" (TZA) was better than any other (except in one case the ordinary Zeldovich approximation) over a wide range from linear to mildly nonlinear ($\sigma \sim 3$) regimes. TZA sets Fourier amplitudes equal to zero for {\it all} wavenumbers greater than $k_{n\ell}$, where $k_{n\ell}$ marks the transition to the nonlinear regime. Here, we study crosscorrelation of generalized TZA with a group of $n$--body simulations for three shapes of window function: sharp $k$--truncation (as in CMS), tophat in coordinate space, or a Gaussian. We also study the crosscorrelation as a function of initial scale within each window type. We find $k$--truncation, which was so much better than other things tried in CMS, is the {it worst} of these three window shapes. We find that a Gaussian window $e^{-k^2/2k_G^2}$ applied to the intial Fourier amplitudes is the best choice. It produces a greatly improved crosscorrelation all cases we studied. The optimum choice of $k_G$ for the Gaussian window is (spectrum-- dependent) 1--1.5 times $k_{n\ell}$, with $k_{n\ell}$ defined by (3). Although all three windows produce similar power spectra and density distribution functions after application of the Zeldovich approximation, phase agreement with the $n$--body simulation is better for the Gaussian window. We ascribe Gaussian window success to its superior treatment of phase evolution.

Davidsen et al. (1991) have argued that the failure to detect uv photons from the dark matter DM) in cluster A665 excludes the decaying neutrino hypothesis. Sciama et al. (1993) argued that because of high central concentration the DM in that cluster must be baryonic. We study the DM profile in clusters of galaxies simulated using the Harrison--Zel'dovich spectrum of density fluctuations, and an amplitude previously derived from numerical simulations (Melott 1984b; Anninos et al. 1991) and in agreement with microwave background fluctuations (Smoot et al. 1992). We find that with this amplitude normalization cluster neutrino DM densities are comparable to observed cluster DM values. We conclude that given this normalization, the cluster DM should be at least largely composed of neutrinos. The constraint of Davidsen et al. can be somewhat weakened by the presence of baryonic DM; but it cannot be eliminated given our assumptions.

We present in this paper a quantitative method for defining void size in large-scale structure based on percolation threshold density. Beginning with two-dimensional gravitational clustering simulations smoothed to the threshold of nonlinearity, we perform percolation analysis to determine the large scale structure. The resulting objective definition of voids has a natural scaling property, is topologically interesting, and can be applied immediately to redshift surveys.

We quantitatively compare a particle implementation of the adhesion approximation to fully non--linear, numerical nbody simulations. Our primary tool, cross--correlation of nbody simulations with the adhesion approximation, indicates good agreement, better than that found by the same test performed with the Zel'dovich approximation (hereafter ZA). However, the cross--correlation is not as good as that of the truncated Zel'dovich approximation (TZA), obtained by applying the Zel'dovich approximation after smoothing the initial density field with a Gaussian filter. We confirm that the adhesion approximation produces an excessively filamentary distribution. Relative to the nbody results, we also find that: (a) the power spectrum obtained from the adhesion approximation is more accurate than that from ZA or TZA, (b) the error in the phase angle of Fourier components is worse than that from TZA, and (c) the mass distribution function is more accurate than that from ZA or TZA. It appears that adhesion performs well statistically, but that TZA is more accurate dynamically, in the sense of moving mass to the right place. Subject Heading: Galaxies, formation, clustering--large--scale structure of the Universe

We present results showing an improvement of the accuracy of perturbation theory as applied to cosmological structure formation for a useful range of quasilinear scales. The Lagrangian theory of gravitational instability of an Einstein-de Sitter dust cosmogony investigated and solved up to the third order in the series of papers by Buchert (1989, 1992, 1993a), Buchert \& Ehlers (1993), Buchert (1993b), Ehlers \& Buchert (1993), is compared with numerical simulations. In this paper we study the dynamics of pancake models as a first step. In previous work (Coles \etal 1993, Melott \etal 1993, Melott 1993) the accuracy of several analytical approximations for the modeling of large-scale structure in the mildly non-linear regime was analyzed in the same way, allowing for direct comparison of the accuracy of various approximations. In particular, the ``Zel'dovich approximation'' (Zel'dovich 1970, 1973, hereafter ZA) as a subclass of the first-order Lagrangian perturbation solutions was found to provide an excellent approximation to the density field in the mildly non-linear regime (i.e. up to a linear r.m.s. density contrast of $\sigma \approx 2$). The performance of ZA in hierarchical clustering models can be greatly improved by truncating the initial power spectrum (smoothing the initial data). We here explore whether this approximation can be further improved with higher-order corrections in the displacement mapping from homogeneity. We study a single pancake model (truncated power-spectrum with power-index $n=-1$) using cross-correlation statistics employed in

We investigate the accuracy of the frozen--flow approximation (FFA), recently proposed by Matarrese \etal (1992), for following the nonlinear evolution of cosmological density fluctuations under gravitational instability. We compare a number of statistics between results of the FFA and nbody simulations, including those used by Melott, Pellman \& Shandarin (1993) to test the Zel'dovich approximation. The FFA performs reasonably well in a statistical sense, e.g. in reproducing the counts--in--cell distribution, at small scales, but it does poorly in the crosscorrelation with nbody which means it is generally not moving mass to the right place, especially in models with high small--scale power.

We calculate reduced moments $\overline \xi_q$ of the matter density fluctuations, up to order $q=5$, from counts in cells produced by Particle--Mesh numerical simulations with scale--free Gaussian initial conditions. We use power--law spectra $P(k) \propto k^n$ with indices $n=-3,~-2,~-1,~0,~1$. Due to the supposed absence of characteristic times or scales in our models, all quantities are expected to depend on a single scaling variable. For each model, the moments at all times can be expressed in terms of the variance $\overline \xi_2$, alone. We look for agreement with the hierarchical scaling ansatz, according to which $\overline \xi_q \propto \overline \xi_2^{~q-1}$. For $n\leq -2$ models we find strong deviations from the hierarchy, which are mostly due to the presence of boundary problems in the simulations. A small, residual signal of deviation from the hierarchical scaling is however also found in $n \geq -1$ models. For the first time, due to our large dynamic range and careful checks of scaling and shot--noise effects, we are able to detect evolution away from the perturbation theory result.