We have previously focused on studying the electron-capture isotopes within the dynamic spiral-arms model and empirically derived the energy dependence of the electron attachment rate using the observation of $^{49}$Ti/$^{49}$V and $^{51}$V/$^{51}$Cr ratios in cosmic rays (Benyamin et al. 2017). We have also shown how this relation recovers the energy dependence seen in the lab measurements (Letaw et al. 1985). In this work we use this relation to construct the $^{44}$Ca/$^{44}$Ti ratio and place a lower limit on the amount of $^{44}$Ti that is required to be nucleosynthesized at the source. The results also imply that the acceleration process of the radioisotopes cannot be much longer than a century time scale (or else the required nucleosynthesized amount has to be correspondingly larger). We also provide a similar lower limit on the source $^{60}$Fe by comparing to the recently observed $^{60}$Fe/$^{56}$Fe (Binns et al. 2016).

We develop a fully three dimensional numerical code describing the diffusion of cosmic rays in the Milky Way. It includes the nuclear spallation chain up to Oxygen, and allows the study of various cosmic ray properties, such as the CR age, grammage traversed, and the ratio between secondary and primary particles. This code enables us to explore a model in which a large fraction of the cosmic ray acceleration takes place in the vicinity of galactic spiral arms and that these spiral arms are dynamic. We show that the effect of having dynamic spiral arms is to limit the age of cosmic rays at low energies. This is because at low energies the time since the last spiral arm passage governs the Cosmic Ray (CR) age, and not diffusion. Using the model, the observed spectral dependence of the secondary to primary ratio is recovered without requiring any further assumptions such as a galactic wind, re-acceleration or various assumptions on the diffusivity. In particular, we obtain a secondary to primary ratio which increases with energy below about 1 GeV.

Electron capture (EC) isotopes are known to provide constraints on the low energy behavior of cosmic rays (CRs), such as re-acceleration. Here we study the EC isotopes within the framework of the dynamic spiral-arms CR propagation model in which most of the CR sources reside in the galactic spiral arms. The model was previously used to explain the B/C and sub-Fe/Fe ratios \citep{BoverC,Iron}. We show that the known inconsistency between the $^{49}$Ti/$^{49}$V and $^{51}$V/$^{51}$Cr ratios remains also in the spiral-arms model. On the other hand, unlike the general wisdom in which the isotope ratios depend primarily on reacceleration, we find here that the ratio also depends on the halo size ($z_\mathrm{h}$) and in spiral-arms models also on the time since the last spiral arm passage ($\tau_\mathrm{arm}$). Namely, EC isotopes can in principle provide interesting constraints on the diffusion geometry. However, with the present uncertainties in the lab measurements of both the electron attachment rate and the fragmentation cross-sections, no meaningful constraint can be placed.

Every supernova hitherto observed has been considered to be the terminal explosion of a star. Moreover, all supernovae with absorption lines in their spectra show those lines decreasing in velocity over time, as the ejecta expand and thin, revealing slower moving material that was previously hidden. In addition, every supernova that exhibits the absorption lines of hydrogen has one main light-curve peak, or a plateau in luminosity, lasting approximately 100 days before declining. Here we report observations of iPTF14hls, an event that has spectra identical to a hydrogen-rich core-collapse supernova, but characteristics that differ extensively from those of known supernovae. The light curve has at least five peaks and remains bright for more than 600 days; the absorption lines show little to no decrease in velocity; and the radius of the line-forming region is more than an order of magnitude bigger than the radius of the photosphere derived from the continuum emission. These characteristics are consistent with a shell of several tens of solar masses ejected by the star at supernova-level energies a few hundred days before a terminal explosion. Another possible eruption was recorded at the same position in 1954. Multiple energetic pre-supernova eruptions are expected to occur in stars of 95-130 solar masses, which experience the pulsational pair instability. That model, however, does not account for the continued presence of hydrogen, or the energetics observed here. Another mechanism for the violent ejection of mass in massive stars may be required.

Most of the diffuse Galactic GeV gamma-ray emission is produced via collisions of cosmic ray (CR) protons with ISM protons. As such the observed spectra of the gamma-rays and the CRs should be strongly linked. Recent observations of Fermi-LAT exhibit a hardening of the gamma-ray spectrum at around a hundred GeV, between the Sagittarius and Carina tangents, and a further hardening at a few degrees above and below the Galactic plane. However, standard CR propagation models that assume a time independent source distribution and a location independent diffusion cannot give rise to a spatially dependent CR (and hence gamma-ray) spectral slopes. Here we consider a dynamic spiral arm model in which the distribution of CR sources is concentrated in the (dynamic) spiral arms, and we study the effects of this model on the $\pi^0$-decay produced gamma-ray spectra. Within this model, near the Galactic arms the observed gamma-ray spectral slope is not trivially related to the CR injection spectrum and energy dependence of the diffusion coefficient. We find unique signatures that agree with the Fermi-LAT observations. This model also provides a physical explanation for the difference between the local CR spectral slope and the CR slope inferred from the average gamma-ray spectrum.

The detection by Rest et al. (2012) of light echoes from $\eta$ Carinae has provided important new observational constraints on the nature of its 1840's era giant eruption. Spectra of the echoes suggest a relatively cool spectral temperature of about 5500K, lower than the lower limit of about 7000K suggested in the optically thick wind outflow analysis of Davidson (1987). This has lead to a debate about the viability of this steady wind model relative to alternative, explosive scenarios. Here we present an updated analysis of the wind outflow model using newer low-temperature opacity tabulations and accounting for the stronger mass loss implied by the $>$10 Msun mass now inferred for the Homunculus. A major conclusion is that, because of the sharp drop in opacity due to free electron recombination for $T<$6500K, a low temperature of about 5000K is compatible with, and indeed expected from, a wind with the extreme mass loss inferred for the eruption. Within a spherical gray model in radiative equilibrium, we derive spectral energy distributions for various assumptions for the opacity variation of the wind, providing a basis for comparisons with observed light echo spectra. The scaling results here are also potentially relevant for other highly optically thick outflows, including those from classical novae, giant eruptions of LBVs and SN Type IIn precursors. A broader issue therefore remains whether the complex, variable features observed from such eruptions are better understood in terms of a steady or explosive paradigm, or perhaps a balance of these idealizations.

Using a recent geochemical reconstruction of the Phanerozoic climate which exhibits a 32 Ma oscillation with a phase and the secondary modulation expected from the vertical the motion of the solar system perpendicular to the galactic plane (shaviv et al. 2014), we show that a kinematically cold strongly interacting disk dark matter (dDM) component is necessarily present in the disk. It has a local density $\rho_\mathrm{dDM} = 0.11 \pm 0.03$ M$_{\odot}/$pc$^3$. It is also consistent with the observed constraints on the total gravitating mass and the baryonic components, and it is the natural value borne from the Toomre stability criterion. It also has surface density $\Sigma_\mathrm{dDM} = 15 \pm 5$ M$_{\odot}/$pc$^2$ and a vertical velocity dispersion of $\sigma_{W} = 8.0 \pm 4.5$ km/s. A dense ("dinosaur killing") thin disk is ruled out. The "normal" halo dark matter (hDM) component should then have a local density $\rho_\mathrm{hDM} \lesssim 0.01$ M$_{\odot}/$pc$^3$. If the dDM component follows the baryons, its average density parameter is $\Omega_\mathrm{dDM} = 1.5 \pm 0.5\%$ and it comprises about 1/8 to 1/4 of Milky Way (MW) mass within the solar circle.

We analyze the motion of stars in the direction perpendicular to the galactic plane following a spiral arm passage. We show that the fast change in the vertical galactic potential causes a thermalized distribution to develop a distinctive "ringing"-like non-thermal signature. We use A and F-stars from the extended Hipparocos catalogue to show that a spiral arm passage (or passages), with an amplitude (or randomly combined amplitudes) of at least {\delta}{\rho} / {\rho} > ~0.15 must have taken place in the past (with {\rho} being the total background density). Presently, the local stellar distribution within ~100 pc of the plane, appears (at the 2.5{\sigma} level) to be contracting towards it.

The progenitors of some supernovae (SNe) exhibit outbursts with super-Eddington luminosities prior to their final explosions. This behavior is common among Type IIn SNe, but the driving mechanisms of these precursors are not yet well understood. SNHunt 275 was announced as a possible new SN during May 2015. Here we report on pre-explosion observations of the location of this event by the Palomar Transient Factory (PTF) and report the detection of a precursor about 500 days prior to the 2015 May activity (PTF 13efv). The observed velocities in the 2015 transient and its 2013 precursor absorption spectra are low (1000-2000 km/s), so it is not clear yet if the recent activity indeed marks the final disruption of the progenitor. Regardless of the nature of this event, we use the PTF photometric and spectral observations, as well as Swift-UVOT observations, to constrain the efficiency of the radiated energy relative to the total kinetic energy of the precursor. We find that, using an order-of-magnitude estimate and under the assumption of spherical symmetry, the ratio of the radiated energy to the kinetic energy is in the range of 0.04 to 3400.

The Boron to Carbon (B/C) and sub-Fe/Fe ratios provides an important clue on Cosmic Ray (CR) propagation within the Galaxy. These ratios estimate the grammage that the CR traverse as they propagate from their sources to Earth. Attempts to explain these ratios within the standard CR propagation models require ad hoc modifications and even with those these models necessitate inconsistent grammages to explain both ratios. As an alternative, physically motivated model, we have proposed that CR originate preferably within the galactic spiral arms. CR propagation from dynamic spiral arms has important imprints on various secondary to primary ratios, such as the B/C ratio and the positron fraction. We use our spiral arm diffusion model with the spallation network extended up to Nickel to calculate the sub-Fe/Fe ratio. We show that without any additional parameters the spiral arm model consistently explains both ratios with the same grammage, providing further evidence in favor of this model.

Stars with a different vertical motion relative to the galactic disk have a different average acceleration. According to Modified Newtonian Dynamics (MOND) theories they should therefore have a different average orbital velocity while revolving around the Milky Way. We show that this property can be used to constrain MOND theories by studying stars in the local neighborhood. With the {\sc Hipparcos} dataset we can only place marginal constraints. However, the forthcoming {\sc GAIA} catalogue with its significantly fainter cutoff should allow placing a stringent constraint.

The progenitor stars of several Type IIb supernovae (SNe) show indications for extended hydrogen envelopes. These envelopes might be the outcome of luminous energetic pre-explosion events, so-called precursor eruptions. We use the Palomar Transient Factory (PTF) pre-explosion observations of a sample of 27 nearby Type IIb SNe to look for such precursors during the final years prior to the SN explosion. No precursors are found when combining the observations in 15-day bins, and we calculate the absolute-magnitude-dependent upper limit on the precursor rate. At the 90% confidence level, Type IIb SNe have on average $<0.86$ precursors as bright as absolute $R$-band magnitude $-14$ in the final 3.5 years before the explosion and $<0.56$ events over the final year. In contrast, precursors among SNe IIn have a $\gtrsim 5$ times higher rate. The kinetic energy required to unbind a low-mass stellar envelope is comparable to the radiated energy of a few-weeks-long precursor which would be detectable for the closest SNe in our sample. Therefore, mass ejections, if they are common in such SNe, are radiatively inefficient or have durations longer than months. Indeed, when using 60-day bins a faint precursor candidate is detected prior to SN 2012cs ($\sim2$% false-alarm probability). We also report the detection of the progenitor of SN 2011dh which does not show detectable variability over the final two years before the explosion. The suggested progenitor of SN 2012P is still present, and hence is likely a compact star cluster, or an unrelated object.

We describe a new code to simulate the stellar evolution of a close interacting binary system. It is then used to calculate the evolution of a classical nova system composed of a 1.25 Msun Main-Sequence (MS) star and a 1.0 Msun white dwarf (WD) companion. The system begins as a well separated non-interacting binary system. Initially, the two stars evolve independently of each other. However, Roche lobe overflow begins as the MS star expands on its way to become a Red Giant. We follow the mass accreted onto the WD and the ensuing nuclear runaways for several thousand flashes. The main finding is that the Roche-Lobe mass transfer rate is modulated by oscillations in the MS star, with a period that is somewhat shorter than the thermal time scale of the star. This periodically modulates the rate of thermonuclear flashes on the WD, between once every 12000 yrs, such that the WD can cool, to once every 300 yrs, such that it cannot. The system is further complicated by the secular drift in the secondary modulation. Such secondary modulation could explain systems like T Pyxidis. Last, we find that the overall process of mass gain by the WD has an efficiency of roughly 9%, thus requiring a donor with an initial mass of larger than about 5 Msun MS for an initial 1 Msun WD, if the WD is to reach the Chandrasekhar mass.

There is a growing number of supernovae (SNe), mainly of Type IIn, which present an outburst prior to their presumably final explosion. These precursors may affect the SN display, and are likely related to some poorly charted phenomena in the final stages of stellar evolution. Here we present a sample of 16 SNe IIn for which we have Palomar Transient Factory (PTF) observations obtained prior to the SN explosion. By coadding these images taken prior to the explosion in time bins, we search for precursor events. We find five Type IIn SNe that likely have at least one possible precursor event, three of which are reported here for the first time. For each SN we calculate the control time. Based on this analysis we find that precursor events among SNe IIn are common: at the one-sided 99% confidence level, more than 50% of SNe IIn have at least one pre-explosion outburst that is brighter than absolute magnitude -14, taking place up to 1/3 yr prior to the SN explosion. The average rate of such precursor events during the year prior to the SN explosion is likely larger than one per year, and fainter precursors are possibly even more common. We also find possible correlations between the integrated luminosity of the precursor, and the SN total radiated energy, peak luminosity, and rise time. These correlations are expected if the precursors are mass-ejection events, and the early-time light curve of these SNe is powered by interaction of the SN shock and ejecta with optically thick circumstellar material.

The binary pulsar J0737-3039 is the only known system having two observable pulsars, thus offering a unique laboratory to test general relativity and explore pulsar physics. Based on the low eccentricity and the position within the galactic plane, Piran & Shaviv (2004, 2005) argued that pulsar B had a non-standard formation scenario with little or no mass ejection. They have also predicted that the system would have a very slow proper motion. Pulsar timing measurements (Kramer et al. 2006; Deller et al. 2009) confirmed this prediction. The recent observations of the alignment between the spin of pulsar A and the binary orbit is also in agreement with this scenario. Detailed simulations of the formation process of pulsar B enable us to show that its progenitor, just before the collapse, was a massive O-Ne-Mg white dwarf surrounded by a tenuous, 0.1-0.16 M_sun, envelope. This envelope was ejected when the white dwarf collapsed to form a neutron star. Pulsar B was born as a slow rotator (spin period ~ 1 s) and a kick received when the pulsar formed changed its spin direction to the current one. This realization sheds light on the angular momentum evolution of the progenitor star, a process which is strongly affected by interaction with the binary companion. The slow proper motion of the system also implies that the system must have undergone a phase of mass transfer in which Star A shed a significant fraction of its mass onto B.

Turbulent processes in the convective envelopes of the sun and stars have been shown to be a source of internal acoustic excitations. In single stars, acoustic waves having frequencies below a certain cutoff frequency propagate nearly adiabatically and are effectively trapped below the photosphere where they are internally reflected. This reflection essentially occurs where the local wavelength becomes comparable to the pressure scale height. In close binary stars, the sound speed is a constant on equipotentials, while the pressure scale height, which depends on the local effective gravity, varies on equipotentials and may be much greater near the inner Lagrangian point (L_1). As a result, waves reaching the vicinity of L_1 may propagate unimpeded into low density regions, where they tend to dissipate quickly due to non-linear and radiative effects. We study the three dimensional propagation and enhanced damping of such waves inside a set of close binary stellar models using a WKB approximation of the acoustic field. We find that these waves can have much higher damping rates in close binaries, compared to their non-binary counterparts. We also find that the relative distribution of acoustic energy density at the visible surface of close binaries develops a ring-like feature at specific acoustic frequencies and binary separations.

Various lines of evidence suggest that very massive stars experience extreme mass-loss episodes shortly before they explode as a supernova. Interestingly, several models predict such pre-explosion outbursts. Establishing a causal connection between these mass-loss episodes and the final supernova explosion will provide a novel way to study pre-supernova massive-star evolution. Here we report on observations of a remarkable mass-loss event detected 40 days prior to the explosion of the Type IIn supernova SN 2010mc (PTF 10tel). Our photometric and spectroscopic data suggest that this event is a result of an energetic outburst, radiating at least 6x10^47 erg of energy, and releasing about 0.01 Solar mass at typical velocities of 2000 km/s. We show that the temporal proximity of the mass-loss outburst and the supernova explosion implies a causal connection between them. Moreover, we find that the outburst luminosity and velocity are consistent with the predictions of the wave-driven pulsation model and disfavor alternative suggestions.

We study the dynamics of continuum driven winds from rotating stars, and develop an approximate analytical model. We then discuss the evolution of stellar angular momentum, and show that just above the Eddington limit, the winds are sufficiently concentrated towards the poles to spin up the star. A twin-lobe structure of the ejected nebula is seen to be a generic consequence of critical rotation. We find that if the pressure in such stars is sufficiently dominated by radiation, an equatorial ejection of mass will occur during eruptions. These results are then applied to {\eta}-Carinae. We show that if it began its life with a high enough angular momentum, the present day wind could have driven the star towards critical rotation, if it is the dominant mode of mass loss. We find that the shape and size of the Homunculus nebula, as given by our model, agree with recent observations. Moreover, the contraction expected due to the sudden increase in luminosity at the onset of the Great Eruption explains the equatorial "skirt" as well.

We study C/O white dwarfs with masses of 1.0 to 1.4 Msun accreting solar-composition material at very high accretion rates. We address the secular changes in the WDs, and in particular, the question whether accretion and the thermonuclear runaways result is net accretion or erosion. The present calculation is unique in that it follows a large number of cycles, thus revealing the secular evolution of the WD system. We find that counter to previous studies, accretion does not give rise to steady state burning. Instead, it produces cyclic thermonuclear runaways of two types. During most of the evolution, many small cycles of hydrogen ignition and burning build a helium layer over the surface of the white dwarf. This He layer gradually thickens and progressively becomes more degenerate. Once a sufficient amount of He has accumulated, several very large helium burning flashes take place and expel the accreted envelope, leaving no net mass accumulation. The results imply that such a system will not undergo an accretion induced collapse, nor will it lead to a SN Type Ia, unless a major new physical process is found.

Supermassive stars (SMS) are massive hydrogen objects, slowly radiating their gravitational binding energy. Such hypothetical primordial objects may have been the seed of the massive black holes (BHs) observed at the centre of galaxies. Under the standard picture, these objects can be approximately described as n=3 polytropes, and they are expected to shine extremely close to their Eddington luminosity. Once however, one considers the porosity induced by instabilities near the Eddington limit, which give rise to super-Eddington states, the standard picture should be modified. We study the structure, evolution and mass loss of these objects. We find the following. First, the evolution of SMSs is hastened due to their increased energy release. They accelerate continuum driven winds. If there is no rotational stabilization, these winds are insufficient to "evaporate" the objects, such that they can collapse to form a supermassive BHs, however, they do prevent SMSs from emitting a copious amount of ionizing radiation. If the SMSs are rotationally stabilized, the winds "evaporate" the objects until a normal sub-Eddington star remains, having a mass of a few 100Msun.

We generalize the problem of the semi-gray model to cases in which a non-negligible fraction of the stellar radiation falls on the long-wavelength range, and/or that the planetary long-wavelength emission penetrates into the transparent short wavelength domain of the absorption. Second, applying the most general assumptions and independently of any particular properties of an absorber, we show that the greenhouse effect saturates and any Earth-like planet has a maximal temperature which depends on the type of and distance to its main-sequence star, its albedo and the primary atmospheric components which determine the cutoff frequency below which the atmosphere is optically thick. For example, a hypothetical convection-less planet similar to Venus, that is optically thin in the visible, could have at most a surface temperature of 1200-1300K irrespective of the nature of the greenhouse gas. We show that two primary mechanisms are responsible for the saturation of the runaway greenhouse effect, depending on the value of the wavelength above which the atmosphere becomes optically thick. Unless this wavelength is small and resides in the optical region, saturation is achieved by radiating the thermal flux of the planet through the short wavelength tail of the thermal distribution. This has the observational implication, the radiation from such a planet should be skewed towards the NIR. Otherwise, saturation takes place by radiating through windows in the FIR.

We consider a scenario where supermassive black holes form through direct accumulation of gas at the centre of proto-galaxies. In the first stage, the accumulated gas forms a super-massive star whose core collapses when the nuclear fuel is exhausted, forming a black hole of $M_{\rm BH} \approx 100 M_{\sun}$. As the black hole starts accreting, it inflates the surrounding dense gas into an almost hydrostatic self-gravitating envelope, with at least 10-100 times the mass of the hole. We find that these "quasistars" suffer extremely high rates of mass loss through winds from their envelopes, in analogy to very massive stars such as eta-Carinae. Only for envelope masses greater than 2.8 \times 10^{5} (M_{\rm BH}/100 M_{\sun})^{9/11} is the envelope evaporation time-scale longer than the accretion time-scale of the black hole. This relation thus constitutes a "threshold growth line" above which quasistars can grow their internal black holes. Accretion rates can be 10 to 100 times the Eddington rate. The quasistars born in this "growth region" with 10^7-10^8 M_{\sun} can grow black holes with masses between 10^4$ to $10^5 M_{\sun}, before crossing the threshold growth line and dispersing their envelopes in less than $10^4$ yr. This scenario therefore predicts that massive black hole seeds can be found only in dark matter halos with total masses larger than about 10^9 M_{\sun}, which can provide sufficiently high accretion rates to form such massive quasistars.

We construct Super-Eddington Slim Disks models around both stellar and super-massive black holes by allowing the formation of a porous layer with a reduced effective opacity. We show that at high accretion rates, the inner part of the disks become radiation pressure dominated. However, unlike the standard scenario in which the disks become thick, super-Eddington disks remain slim. In addition, they accelerate a significant wind with a "thick disk" geometry. We show that above about 1.5 times the standard critical mass accretion rate (needed to release the Eddington luminosity), the net luminosity released is above Eddington. At above about 5 times the standard critical rate, the central BH accretes more than the Eddington accretion rate. Above about 20 m-dot_crit, the disk remains slim but the wind becomes spherical, and the present model breaks down.

Recent measurements of the positron/electron ratio in the cosmic ray (CR) flux exhibits an apparent anomaly, whereby this ratio increases between 10 and 100 GeV. In contrast, this ratio should decrease according to the standard scenario, in which CR positrons are secondaries formed by hadronic interactions between the primary CR protons and the interstellar medium (ISM). The positron excess is therefore interpreted as evidence for either an annihilation/decay of weakly interacting massive particles, or for a direct astrophysical source of pairs. The common feature of all proposed models is that they invoke new physics or new astrophysical sources. However, this line of argumentation relies implicitly on the assumption of a relatively homogeneous CR source distribution. Inhomogeneity of CR sources on a scale of order a kpc, can naturally explain this anomaly. If the nearest major CR source is about a kpc away, then low energy electrons (~GeV) can easily reach us. At higher energies (> 10 GeV), the source electrons cool via synchrotron and inverse-Compton before reaching the solar vicinity. Pairs formed in the local vicinity through the proton/ISM interactions can reach the solar system also at high energies, thus increasing the positron/electron ratio. A natural origin of source inhomogeneity is the strong concentration of supernovae to the galactic spiral arms. Assuming supernova remnants (SNRs) as the sole primary source of CRs, and taking into account their concentration near the galactic spiral arms, we consistently predict the observed positron fraction between 1 and 100 GeV, while abiding to different constraints such as the observed electron spectrum and the CRs cosmogenic age. An ATIC like spectral excess at ~600 GeV can be explained, in this picture, as the contribution of a few known nearby SNRs.

Recent measurements of the positron/electron ratio in the cosmic ray (CR) flux exhibits an apparent anomaly, whereby this ratio increases between 10 and 100 GeV. We show that inhomogeneity of CR sources on a scale of order a kpc, can naturally explain this anomaly. If the nearest major CR source is about a kpc away, then low energy electrons ($\sim 1$ GeV) can easily reach us. At higher energies ($\gtrsim 10$ GeV), the source electrons cool via synchrotron and inverse-Compton before reaching Earth. Pairs formed in the local vicinity through the proton/ISM interactions can reach Earth also at high energies, thus increasing the positron/electron ratio. A natural origin of source inhomogeneity is the strong concentration of supernovae in the galactic spiral arms. Assuming supernova remnants (SNRs) as the sole primary source of CRs, and taking into account their concentration near the galactic spiral arms, we consistently recover the observed positron fraction between 1 and 100 GeV. ATIC's electron excess at $\sim 600$ GeV is explained, in this picture, as the contribution of a few known nearby SNRs. The apparent coincident similarity between the cooling time of electrons at 10 GeV (where the positron/electron ratio upturn), $\sim 10$ Myr, and the CRs protons cosmogenic age at the same energy is predicted by this model.

Stars can produce steady-state winds through radiative driving as long as the mechanical luminosity of the wind does not exceed the radiative luminosity at its base. This upper bound on the mass loss rate is known as the photon-tiring limit. Once above this limit, the radiation field is unable to lift all the material out of the gravitational potential of the star, such that only part of it can escape and reach infinity. The rest stalls and falls back toward the stellar surface, making a steady-state wind impossible. Photon-tiring is not an issue for line-driven winds since they cannot achieve sufficiently high mass loss rates. It can however become important if the star exceeds the Eddington limit and continuum interaction becomes the dominant driving mechanism. This paper investigates the time-dependent behaviour of stellar winds that exceed the photon-tiring limit, using 1-D numerical simulations of a porosity moderated, continuum-driven stellar wind. We find that the regions close to the star show a hierarchical pattern of high density shells moving back and forth, unable to escape the gravitational potential of the star. At larger distances, the flow eventually becomes uniformly outward, though still quite variable. Typically, these winds have a very high density but a terminal flow speed well below the escape speed at the stellar surface. Since most of the radiative luminosity of the star is used to drive the stellar wind, such stars would appear much dimmer than expected from the super-Eddington energy generation at their core. The visible luminosity typically constitutes less then half of the total energy flow and can become as low as ten percent or less for those stars that exceed the photon-tiring limit by a large margin.

We present the results of numerical simulations of continuum-driven winds of stars that exceed the Eddington limit and compare these against predictions from earlier analytical solutions. Our models are based on the assumption that the stellar atmosphere consists of clumped matter, where the individual clumps have a much larger optical thickness than the matter between the clumps. This `porosity' of the stellar atmosphere reduces the coupling between radiation and matter, since photons tend to escape through the more tenuous gas between the clumps. This allows a star that formally exceeds the Eddington limit to remain stable, yet produce a steady outflow from the region where the clumps become optically thin. We have made a parameter study of wind models for a variety of input conditions in order to explore the properties of continuum-driven winds. The results show that the numerical simulations reproduce quite closely the analytical scalings. The mass loss rates produced in our models are much larger than can be achieved by line driving. This makes continuum driving a good mechanism to explain the large mass loss and flow speeds of giant outbursts, as observed in eta Carinae and other luminous blue variable (LBV) stars. Continuum driving may also be important in population III stars, since line driving becomes ineffective at low metalicities. We also explore the effect of photon tiring and the limits it places on the wind parameters.

Continuum driving is an effective method to drive a strong stellar wind. It is governed by two limits: the Eddington limit and the photon-tiring limit. A star must exceed the effective Eddington limit for continuum driving to overcome the stellar gravity. The photon-tiring limit places an upper limit on the mass loss rate that can be driven to infinity, given the energy available in the radiation field of the star. Since continuum driving does not require the presence of metals in the stellar atmosphere it is particularly suited to removing mass from low- and zero-metallicity stars and can play a crucial part in their evolution. Using a porosity length formalism we compute numerical simulations of super-Eddington, continuum driven winds to explore their behaviour for stars both below and above the photon-tiring limit. We find that below the photon tiring limit, continuum driving can produce a large, steady mass loss rate at velocities on the order of the escape velocity. If the star exceeds the photon-tiring limit, a steady solution is no longer possible. While the effective mass loss rate is still very large, the wind velocity is much smaller

Shortly after the discovery of the binary system PSR J0737-3039 we predicted that it has a small proper motion and suggested that this implies that the progenitor mass of the younger pulsar (B) must be around 1.45M_sun. This is in contradiction to standard evolutionary scenarios that suggest a He-star progenitor with M >2.1-2.3M_sun and in sharp contrast to a velocity of ~100 km/s predicted by these models. It requires a new previously unseen type of gravitational collapse. Our (generally ignored) prediction of a low peculiar motion was confirmed later by pulsar timing observations that found v_\bot < 30 km/s. Recently Willems et al. put forward an elaborate simulation suggesting that in spite of these observations, the most likely natal kick is large (50-170 km/s) and the corresponding most likely progenitor mass is large as well. We refute this claim and show that Willems et al. implicitly assume the result that they obtain.

It is generally accepted that neutron stars form in core collapse events that are accompanied by a supernovae (types II or Ib or Ic). Typical progenitors are, therefore, larger than $\sim 2.1 \Ms$. We suggest \cite{PS04,PS05} that the binary pulsar J0737-3039 provides evidence for a new formation channel: collapse of a light progenitor. This binary pulsar J0737-3039 has several remarkable features including among others: a very tight orbit with a Keplerian velocity of ~600km/sec, a low eccentricity, and a location ~50pc from the Galactic plane implying that the system has, at high likelihood, a small (compared to Keplerian) center of mass velocity. A significant mass loss during the formation of the second pulsar would have lead either to an eccentric orbit or to a large center of mass velocity or to both. Therefore, we can set a strong upper limit on the progenitor's mass. A progenitor more massive than 1.9$\Ms$ is ruled out (at 97% confidence). The kinematically favored option is of a progenitor mass around 1.45$\Ms$. Recent evidence for a rather low velocity proper motion supports this prediction and decreases the likelihood of the standard (high mass progenitor) scenario. Lack of variations in the pulses' profiles (which indicate no significant geodetic precession) provides further support for the no kick and low progenitor mass formation scenario.

The Milky Way spiral arm dynamics is studied using the birth place of open clusters. We separately study the nearby spiral arms, and find evidence for multiple spiral sets. In particular, the Sagittarius-Carina arm appears to be a superposition of two sets. The first has a pattern speed of Omega_P,Carina,1=16.5^{+1.2}_{-1.4sys}+/-1.1_{stat} (km/s)/kpc, while the second has Omega_P,Carina,2=29.8^{+0.6}_{-1.3sys}+/-1.3_{stat} (km/s)/kpc. The Perseus arm located at larger galactic radii exhibits only one clear set with Omega_P,Perseus=20.0^{+1.7}_{-1.2sys}+/-1.6_{stat} (km/s)/kpc. The Perseus and "slower" Carina arms are most probably part of an m=4 set. The Orion "armlet" appears to have a pattern speed of Omega_P,Orion=28.9 ^{+1.3}_{-1.2sys}+/-0.8_{stat} (km/s)/kpc. In other words, both the "faster" Carina arm and the Orion arm are nearly corotating with the solar system.

In the presence of strong magnetic fields, the vacuum becomes a birefringent medium. We show that this QED effect couples the direction of the polarization of photons leaving the NS surface, to the direction of the magnetic field along the ray's path. We analyze the consequences that this effect has on aligning the polarization vectors to generate large net polarizations, while considering thermal radiation originating from a thermal hydrogen atmosphere. Counter to previous predictions, we show that the thermal radiation should be highly polarized even in the optical. When detected, this polarization will be the first demonstration of vacuum birefringence. It could be used as a tool to prove the high magnetic field nature of AXPs and it could also be used to constrain physical NS parameters, such as $R/M$, to which the net polarization is sensitive.

Evolutionary scenarios suggest that the progenitor of the new binary pulsar J0737-3039B \cite{ref1,ref2} was a He-star with $M > 2.1-2.3~\Ms$ \cite{ref3,ref4}. We show that this case implies that the binary must have a large ($>120$ km/s) center of mass velocity. However, the location, $\sim 50$ pc from the Galactic plane, suggests that the system has, at high likelihood, a significantly smaller center of mass velocity and a progenitor more massive than 2.1~$\Ms$ is ruled out (at 97% c.l.). A progenitor mass around 1.45~$\Ms$, involving a new previously unseen gravitational collapse, is kinematically favored. The low mass progenitor is consistent with the recent scintillations based velocity measurement of 66$\pm 15$ km/s \cite{ref12new} (and which also rules out the high mass solution at 99% c.l.) and inconsistent with the higher earlier estimates of 141$\pm 8.5$ km/s \cite{ref11new}. Direct proper motion measurements, that should be available within a year or so, should better help to distinguish between the two scenarios.

We examine radiatively driven mass loss from stars near and above the Eddington limit (Ledd). We begin by reviewing the instabilities that are expected to form extensive structure near Ledd. We investigate how this "porosity" can reduce the effective coupling between the matter and radiation. Introducing a new "porosity-length'' formalism, we derive a simple scaling for the reduced effective opacity, and use this to derive an associated scaling for the porosity-moderated, continuum-driven mass loss rate from stars that formally exceed Ledd. For a simple super-Eddington model with a single porosity length that is assumed to be on the order of the gravitational scale height, the overall mass loss is similar to that derived in previous porosity work. This is much higher than is typical of line-driven winds, but is still only a few percent of the photon tiring limit--for which the luminosity becomes insufficient to carry the flow out of the gravitational potential. To obtain still stronger mass loss that approaches observationally inferred values near this limit, we introduce a power-law-porosity model in which the associated structure has a broad range of scales. We show that the mass loss rate can be enhanced by a factor that increases with the Eddington parameter Gamma, such that for moderately large Gamma (> 3-4), mass loss rates could approach the photon tiring limit. Together with the ability to drive quite fast outflow speeds, the derived mass loss could explain the large inferred mass loss and flow speeds of giant outbursts in eta Carinae and other LBV stars.

The observed orbital parameters of the recently discovered binary pulsar J0737-3039 are used to constrain the progenitor system. In particular, the small observed eccentricity and the small inferred peculiar velocity imply that during the formation of the second NS in the system, a very small mass Delta M <~ 0.15 Msun, was ejected. A progenitor more massive than 1.4 Msun is unlikely and a progenitor more massive than 2.3 Msun is practically ruled out. We therefore argue that the companion removed most of the progenitor envelope prior to its collapse. The second collapse may have been a SN of Type Ib, or more probably, of Type Ic. Alternatively, it could have been a bare collapse, an option which is kinematically favored. Future constraints on the peculiar motion could solidify this constraint on the mass loss. A future determination of the pulsars' spin vectors and proper motion could allow us the complete determination of the conditions during the second collapse.

Standard solar models predict a solar luminosity that gradually increased by about 30% over the past 4.5 billion years. Under the faint sun, Earth should have been frozen solid for most of its existence. Yet, running water is observed to have been present since very early in Earth's history. This enigma is known as the faint sun paradox. We show here that it can be partially resolved once we consider the cooling effect that cosmic rays are suspected to have on the global climate and that the younger sun must have had a stronger solar wind, such that it was more effective at stopping cosmic rays from reaching Earth. The paradox can then be completely resolved with the further contribution of modest greenhouse gas warming. When we add the cosmic ray flux modulation by a variable star formation rate in the Milky Way, we recover the long term glacial activity on Earth. As to the future, we find that the average global temperature will increase by typically 10K in the coming 2 Gyr.

We revisit the evaporation process of gas from dwarf galaxies after it has been photo-ionized by the UV flux from the first stars and AGNs and heated to T~10^4K or 2x10^4K respectively. Earlier estimates, based on the balance between pressure and gravity, indicated that dark haloes of virial velocity lower than Vevap ~ 11-13 km/s have lost most of their gas in a dynamical time. We follow the continuous evaporation by a thermal wind during the period when the ionizing flux was effective. We find that the critical virial velocity for significant evaporation is significantly higher. For example, if the ionization starts at z-ion=10 and is maintained until z=2, a mass loss of one e-fold occurs in haloes of Vevap ~ 25 (or 35 km/s) for T ~ 10^4K (or 2x10^4K). Haloes of Vevap ~ 21 km/s (or 29 km/s) lose one e-fold within the first Hubble time at z=10. Any dwarf galaxies with virial velocities smaller than Vevap must have formed their stars from a small fraction of their gas before z-ion, and then lost the rest of the gas by photo-evaporation. This may explain the gas-poor, low surface brightness dwarf spheroidal galaxies. By z<1, most of the IGM gas was evaporated at least once form dwarf galaxies, thus providing a lower bound to its metallicity.

In the presence of strong magnetic fields, the vacuum becomes a birefringent medium. We show that this QED effect decouples the polarization modes of photons leaving the NS surface. Both the total intensity and the intensity in each of the two modes is preserved along a ray's path through the neutron-star magnetosphere. We analyze the consequences that this effect has on aligning the observed polarization vectors across the image of the stellar surface to generate large net polarizations. Counter to previous predictions, we show that the thermal radiation of NSs should be highly polarized even in the optical. When detected, this polarization will be the first demonstration of vacuum birefringence. It could be used as a tool to prove the high magnetic field nature of AXPs and it could also be used to constrain physical NS parameters, such as $R/M$, to which the net polarization is sensitive.

The exact ionization state of Be7 in the solar core is crucial for the precise prediction of the solar B8 neutrino flux. We therefore examine the effect of pressure ionization on the ionization state of Be7 and all elements with 12 >= Z >= 4. We show that under the conditions prevailing in the solar core, one has to consider the effect of the nearest neighbor on the electronic structure of a given ion. To this goal, we first solve the Schroedinger and then the Kohn-Sham equations for an ion immersed in a dense plasma under conditions for which the mean interparticle distance is smaller than the Debye radius. The question of which boundary conditions should be imposed on the wave function is discussed, examined and found to be crucial. Contrary to previous estimates showing that Beryllium is partially ionized, we find that it is fully ionized. As a consequence, the predicted rate of the Be7 + e- reaction is reduced by 20-30%, depending on the exact solar model. Since Be7 is a trace element, its total production is controlled by the unchanged He4+He3 reaction rate, and its destruction is determined by the rate of electron capture. As the latter rate decreases when the Beryllium is fully ionized (relative to the case of partially ionized Be), the estimate for the abundance of Be7 increases and with it the B8 neutrino flux. The increase in phi_nu(B8) is by about 20-30%. The neutrino flux due to Be7 electron capture remains effectively unchanged because the change in the rate is compensated for by a change in the abundance. Hence the prediction for the ratio of phi_nu(B8) / phi_nu(Be7) changes as well.

The short term variability of the Galactic cosmic ray flux (CRF) reaching Earth has been previously associated with variations in the global low altitude cloud cover. This CRF variability arises from changes in the solar wind strength. However, cosmic ray variability also arises intrinsically from variable activity of and motion through the Milky Way. Thus, if indeed the CRF climate connection is real, the increased CRF witnessed while crossing the spiral arms could be responsible for a larger global cloud cover and a reduced temperature, thereby facilitating the occurrences of ice ages. This picture has been recently shown to be supported by various data (Shaviv, 2001). In particular, the variable CRF recorded in Iron meteorites appears to vary synchronously with the appearance ice ages. Here we expand upon the original treatment with a more thorough analysis and more supporting evidence. In particular, we discuss the cosmic ray diffusion model which considers the motion of the Galactic spiral arms. We also elaborate on the structure and dynamics of the Milky Way's spiral arms. In particular, we bring forth new argumentation using HI observations which imply that the galactic spiral arm pattern speed appears to be that which fits the glaciation period and the cosmic-ray flux record extracted from Iron meteorites. In addition, we show that apparent peaks in the star formation rate history, as deduced by several authors, coincides with particularly icy epochs, while the long period of 1 to 2 Gyr before present, during which no glaciations are known to have occurred, coincides with a significant paucity in the past star formation rate.

Several of the inconsistencies plaguing the field of novae are resolved once we consider novae to be steady state super-Eddington objects. In particular, we show that the super-Eddington shell burning state is a natural consequence of the equations of stellar structure, and that the predicted mass loss in the super-Eddington state agrees with nova observations. We also find that the transition phase of novae can be naturally explained as "stagnating" winds.

We construct a Galactic cosmic ray (CR) diffusion model while considering that CR sources reside predominantly in the Galactic spiral arms. We find that the CR flux (CRF) reaching the solar system should periodically increase each crossing of a Galactic spiral arm. We search for this signal in the CR exposure age record of Iron meteorites and confirm this prediction. We then check the hypothesis that climate, and in particular the temperature, is affected by the CRF to the extent that glaciations can be induced or completely hindered by possible climatic variations. We find that although the geological evidence for the occurrence of IAEs in the past Eon is not unequivocal, it appears to have a nontrivial correlation with the spiral arm crossings--agreeing in period and phase. Thus, a better timing study of glaciations could either confirm this result as an explanation to the occurrence of IAEs or refute a CRF climatic connection.

The thermal emission of strongly magnetized neutron-star atmospheres is thought to be highly polarized. However, because of the different orientations of the magnetic field over the surface of the neutron star (NS), it is commonly assumed that the net observed polarization will be significantly reduced as the polarization from different regions will cancel each other. We show that the birefringence of the magnetized QED vacuum decouples the polarization modes in the magnetosphere; therefore, the direction of the polarization follows the direction of the magnetic field up to a large distance from the stellar surface. At this distance, the rays that leave the surface and are destined for our detectors pass through only a small solid angle; consequently, the polarization direction of the emission originating in different regions will tend to align together. The net observed polarization of the thermal radiation of NSs should therefore be very large. Measurement of this polarization will be the first direct evidence of the birefringence of the magnetized vacuum due to QED and a direct probe of behavior of the vacuum at magnetic fields of order of and above the critical QED field of 4.4 x 10 13 G. The large observable polarization will also help us learn more about the atmospheric properties of NSs.

We present a model for steady state winds of systems with super-Eddington luminosities. These radiatively driven winds are expected to be optically thick and clumpy as they arise from an instability driven porous atmosphere. The model is then applied to derive the mass loss observed in bright classical novae. The main results are: 1) A general relation between the mass loss rate and the total luminosity in super-Eddington systems. 2) A quantitative agreement between the observed luminosity evolution which is used to predict both the mass loss and temperature evolution, and their observations. 3) An agreement between the predicted average integrated mass loss of novae as a function of WD mass and its observations. 4) A natural explanation for the `transition phase' of novae. 5) Agreement with eta Carinae which was used to double check the theory. The prediction for the mass shed in the star's great eruption agrees with observations to within the measurement error.

Atmospheres having a significant radiative support are shown to be intrinsically unstable at luminosities above a critical fraction Gamma_crit ~ 0.5-0.85 of the Eddington limit, with the exact value depending on the boundary conditions. Two different types of absolute radiation-hydrodynamic instabilities of acoustic waves are found to take place even in the electron scattering dominated limit. Both instabilities grow over dynamical time scales and both operate on non radial modes. One is stationary and arises only after the effects of the boundary conditions are taken into account, while the second is a propagating wave and is insensitive to the boundary conditions. Although a significant wind can be generated by these instabilities even below the classical Eddington luminosity limit, quasi-stable configurations can exist beyond the Eddington limit due to the generally reduced effective opacity. The study is done using a rigorous numerical linear analysis of a gray plane parallel atmosphere under the Eddington approximation. We also present more simplified analytical explanations.

We present a model for the steady state winds of super-Eddington systems. These radiatively driven winds are expected to be optically thick and clumpy as they arise from a porous atmosphere. The model is then used to derive the mass loss of bright classical novae. The long duration super-Eddington outflows that are clearly observed in at least two cases (Nova LMC 1988 #1 & Nova FH Ser) are naturally explained. Moreover, the predicted mass loss and temperature evolution agree nicely with the observations, as do additional features. eta-Car is then used to double check the theory which predicts the observed mass shed in the great eruption.

We analyze the wind generated by the great 20 year long super-Eddington outburst of eta-Carinae. We show that using classical stellar atmospheres and winds theory, it is impossible to construct a consistent wind model in which a sufficiently small amount of mass, like the one observed, is shed. One expects the super-Eddington luminosity to drive a thick wind with a mass loss rate substantially higher than the observed one. The easiest way to resolve the inconsistency is if we alleviate the implicit notion that atmospheres are homogeneous. An inhomogeneous atmosphere, or "porous", allows more radiation to escape while exerting a smaller average force. Consequently, such an atmosphere yields a considerably lower mass loss rate for the same total luminosity. Moreover, all the applications of the Eddington Luminosity as a strict luminosity limit should be revised, or at least reanalyzed carefully.

Extremely strong magnetic fields change the vacuum index of refraction. Although this polarization dependent effect is small for typical neutron stars, it is large enough to decouple the polarization states of photons traveling within the field. The photon states evolve adiabatically and follow the changing magnetic field direction. The combination of a rotating magnetosphere and a frequency dependent state decoupling predicts polarization phase lags between different wave bands, if the emission process takes place well within the light cylinder. This QED effect may allow observations to distinguish between different pulsar emission mechanisms and to reconstruct the structure of the magnetosphere.

Extremely strong magnetic fields change the vacuum index of refraction. This induces a lensing effect that is not unlike the lensing phenomenon in strong gravitational fields. The main difference between the two is the polarization dependency of the magnetic lensing, a behaviour that induces a handful of interesting effects. The main prediction is that the thermal emission of neutron stars with extremely strong magnetic fields is polarized - up to a few percent for the largest fields known. This potentially allows a direct method for measuring their magnetic fields.

We review and analyze the available information for nuclear fusion cross sections that are most important for solar energy generation and solar neutrino production. We provide best values for the low-energy cross-section factors and, wherever possible, estimates of the uncertainties. We also describe the most important experiments and calculations that are required in order to improve our knowledge of solar fusion rates.

The effects that large scale fluctuations had on small scale isothermal modes at the epoch of recombination are analysed. We find that: (a) Albeit the fact that primordial fluctuations were at this epoch still well in the linear regime, a significant nonlinear radiation hydrodynamic interaction could have taken place. (b) Short wavelength isothermal fluctuations are unstable. Their growth rate is exponential in the amplitude of the large scale fluctuations and is therefore very sensitive to the initial conditions. (c) The observed CMBR fluctuations are of order the limit above which the effect should be significant. Thus, according to their exact value, the effect may be negligible or lead to structure formation out of isothermal fluctuations within the period of recombination. (d) If the cosmological parameters are within the prescribed regime, the effect should be detectable through induced deviations in the Planck spectrum. (e) The sensitivity of the effect to the initial conditions provides a tool to set limits on various cosmological parameters with emphasis on the type and amplitude of the primordial fluctuation spectrum. (f) Under proper conditions, the effect may be responsible for the formation of sub-globular cluster sized objects at particularly high red shifts. (g) Under certain circumstances, it can also affect horizon sized large scale structure.

High energy cosmic ray jets from nearby mergers or accretion induced collapse (AIC) of neutron stars (NS) that hit the atmosphere can produce lethal fluxes of atmospheric muons at ground level, underground and underwater, destroy the ozone layer and radioactivate the environment. They could have caused most of the massive life extinctions on planet Earth in the past 600 My. Biological mutations due to ionizing radiations could have caused the fast appearance of new species after the massive extinctions. An early warning of future extinctions due to NS mergers may be obtained by identifying, mapping and timing all the nearby binary neutron stars systems. A warning of an approaching cosmic ray burst from a nearby NS merger/AIC may be provided by a very intense gamma ray burst which preceeds it.

Cosmic ray bursts (CRBs) from mergers or accretion induced collapse of neutron stars that hit an Earth-like planet closer than $\sim 1 kpc$ from the explosion produce lethal fluxes of atmospheric muons at ground level, underground and underwater. These CRBs also destroy the ozone layer and radioactivate the environment. The mean rate of such life devastating CRBs is one in 100 million years (Myr), consistent with the observed 5 ``great'' extinctions in the past 600 Myr. Unlike the previously suggested extraterrestrial extinction mechanisms the CRBs explain massive life extinction on the ground, underground and underwater and the higher survival levels of radiation resistant species and of terrain sheltered species. More distant mergers can cause smaller extinctions. Biological mutations caused by ionizing radiation produced by the CRB may explain a fast appearance of new species after mass extinctions. The CRB extinction predicts detectable enrichment of rock layers which formed during the extinction periods with cosmogenically produced radioactive nucleides such as $^{129}$I, $^{146}$Sm, $^{205}$Pb with and $^{244}$Pu. Tracks of high energy particles in rock layers on Earth and on the moon may also contain records of intense CRBs. An early warning of future extinctions due to neutron star mergers can be obtained by identifying, mapping and timing all the nearby binary neutron stars systems. A final warning of an approaching CRB from a nearby neutron stars merger will be provided by a gamma ray burst a few days before the arrival of the CRB.

Photoexcitation and ionization of partially ionized heavy atoms in highly relativistic flows by interstellar photons, followed by their reemission in radiative recombination and decay, boost star-light into beamed $\gamma$ rays along the flow direction. Repeated excitation/decay of highly relativistic baryonic ejecta from merger or accretion induced collapse of neutron stars in dense stellar regions (DSRs), like galactic cores, globular clusters and super star-clusters, can convert enough kinetic energy in such events in distant galaxies into cosmological gamma ray bursts (GRBs). The model predicts remarkably well all the main observed temporal and spectral properties of GRBs. Its success strongly suggests that GRBs are $\gamma$ ray tomography pictures of DSRs in galaxies at cosmological distances with unprecedented resolution: A time resolution of $dt\sim 1~ms$ in a GRB can resolve stars at a Hubble distance which are separated by only $D\sim 10^{10}cm$. This is equivalent to the resolving power of an optical telescope with a diameter larger than one Astronomical Unit!

We propose that repeated photoexcitation/ionization of high Z atoms of highly relativistic flows by star light in dense stellar regions followed by emission of decay/recombination photons, which are beamed and boosted to gamma ray energies in the observer frame, produce gamma ray bursts (GRBs). We show that this overlooked mechanism, which is able to convert efficiently baryonic kinetic energy release in merger or accretion induced collapse of neutron stars into cosmological GRBs, reproduces remarkably well all the main observed properties of GRBs.

We study a cosmological scenario for gamma ray bursts (GRBs) where relativistic flows interact with dense radiation fields. It is shown that this scenario is plausible in very dense stellar regions which are known to exist in collapsed cores of globular clusters or dense nuclei of galaxies. It yields a correct quantitative description of the temporal behavior of GRBs. Several other properties of GRBs are easily explained.

Blazar emission of gamma rays and cosmic ray production of gamma rays in gas-rich clusters have been proposed recently as alternative sources of the high energy extragalactic diffuse gamma ray background radiation. We show that these sources also produce a very different high energy extragalactic diffuse neutrino background radiation. This neutrino background may be detected by the new generation of large neutrino telescopes under construction and may be used to trace the origin of the extragalactic gamma radiation.

Most cosmological models for gamma-ray bursts invoke the production of a ``fireball'' in a compact region, as indicated by the short time variability of the observed GRBs. The high density of $e^+e^-$ pairs in such fireballs inevitably makes them opaque to gamma-rays and requires the gamma-ray emission to take place only after the fireball has expanded and became optically thin to gamma's. Many of the ``traditional'' scenarios explain the temporal variability of GRBs as being formed by growing instabilities in the expanding fireball. Here we explore this possiblity and examine its validity.

We show that cosmic rays in external galaxies, groups and clusters rich in gas, with an average flux similar to that observed in the Milky Way, could have produced the observed extragalactic diffuse gamma radiation.

Striking similarities exist between high energy gamma ray emission from active galactic nuclei (AGN) and gamma ray bursts (GRBs). They suggest that GRBs are generated by inverse Compton scattering from highly relativistic electrons in transient jets. Such jets may be produced along the axis of an accretion disk formed around stellar black holes (BH) or neutron stars (NS) in BH-NS and NS-NS mergers and in accretion induced collapse of magnetized white dwarfs (WD) or neutron stars in close binary systems. Such events may produce the cosmological GRBs. Transient jets formed by single old magnetized neutron stars in an extended Galactic halo may produce a local population of GRBs. Here we show that jet production of GRBs by inverse Compton scattering can explain quite simply the striking correlations that exist between various temporal features of GRBs, their duration histogram, the power spectrum of their complex multipeak light curves, their power-law high energy spectra and other features of GRBs. Some additional predictions are made including the expected polarization of gamma-rays in the bursts.