We derive the rates of capture, Ndot, of main sequence turn off stars by the central massive black hole in a sample of galaxies from Magorrian et al. 1998. The disruption rates are smaller than previously believed with Ndot ~ 10^-4 - 10^-7 per galaxy. A correlation between Ndot and black hole mass, M, is exploited to estimate the rate of tidal disruptions in the local universe. Assuming that all or most galaxies have massive black holes in their nuclei, this rate should be dominated by sub-Lstar galaxies. The rate of tidal disruptions could be high enough to be detected in supernova (or similar) monitoring campaigns---we estimate the rate of tidal disruptions to be 0.01 - 0.1 times the supernova rate. We have also estimated the rates of disruption of red giants, which may be significant (Ndot ~> 10^-4 y^-1 per galaxy) for M ~> 10^8 Msun, but are likely to be harder to observe---only of order 10^-4 times the supernova rate in the local universe. In calculating capture rates, we advise caution when applying scaling formulae by other authors, which are not applicable in the physical regime spanned by the galaxies considered here.

Optically-thick envelopes may form following the tidal disruption of a star by a massive black hole. Such envelopes would reprocess hard radiation from accretion close to the black hole into the UV and optical bands producing AGN-luminosity flares with duration ~1 year. We show that due to relativistic effects, the envelopes are convective. If convection is efficient, then the structure of the envelopes is similar to that described in previous work; however, the photospheric radius is shown to be very sensitive to the luminosity at the envelope base, suggesting that either the envelope collapses or the envelope expands to a maximum radius at which point a wind may set in. For an envelope without winds, we find a maximum photospheric radius of ~10^16 cm (i.e. minimum effective temperature ~6,000 K). The evolution of the envelopes is described based on simple energy arguments.

We have determined the location of the line-opacity modified Eddington limit for stars in the LMC using the most recent atmosphere models combined with a precise mapping to the HR Diagram through up-to-date stellar evolution calculations. While we find, in agreement with previous studies, that the shape of the modified Eddington limit qualitatively corresponds to the Humphreys-Davidson (HD) limit defined by the most luminous supergiants, the modified limit is actually a full magnitude higher than the upper luminosity limit observed for LMC stars. The observed limit is consistent with atmosphere models in which the maximum value of the ratio of the radiation force outwards to the gravitational force inwards, Y_max, is 0.9, i.e., the photospheres of stars at the observed luminosity limit are bound. As massive stars evolve, they move to higher, and therefore less stable values of Y_max, so mass loss, either sporadic or continuous, may halt their natural redward evolution as they approach the Y_max = 0.9 limit. We assess the metallicity dependence of this limit. If the limit does determine the most luminous stars, and the value of Y_max corresponding to the luminosity limit in the LMC is universal, then the brightest supergiants the SMC should be only marginally brighter (0.3 mag) than those of the LMC, in agreement with observations. Moreover, the brightest supergiants in M31 should be 0.75~mag fainter than those in the LMC.

Geometrically thick disks may form after tidal disruption events, and rapid accretion may lead to short flares followed by long-term, lower-level emission. Using a novel accretion disk code which relies primarily on global conservation laws and the assumption that viscosity is everywhere positive, a broad range of physically allowed evolutionary sequences of thick disks is investigated. The main result is that accretion in the thick disk phase can consume only a fraction of the initial disk material before the disk cools and becomes thin. This fraction is ~0.5-0.9 for disruptions around 10^6 to 10^7 M_\odot black holes and is sensitive to the mean angular momentum of the disk. The residual material will accrete in some form of thin disk over a longer period of time. The initial thick disk phase may reduce the dimming timescale of the disk by a factor of ~2 from estimates based on thin disks alone. Assuming an 0.5 M_\odot initial thick disk, even if the thin disks become advection dominated, the black hole mass to light ratio can rise above M_\odot/L_\odot = 1 in no less than 20 (0.1/\alpha) to 2000 (0.1/\alpha) years following a tidal disruption event, depending on the mass of the black hole and the initial conditions of the encounter. The long-term emission will be most prevalent around lower mass, 10^6 M_\odot black holes. If the tidal disruption rates in these galactic nuclei are ~10^{-4} yr^{-1}, then about 10% of the nuclei should exhibit the long-term UV/optical emission at a level of ~10^{38} ergs~s^{-1}.

Tidal disruption flares are differentiated into two classes -- those which are sub-Eddington and those which radiate near the Eddington limit. Flares from black holes above ~2 x 10^7 M_\odot will generally not radiate above the Eddington limit. For a Schwarzschild black hole, the maximum bolometric luminosity of a tidal disruption is ~L_Edd(5 x 10^7 M_\odot), substantially below the Eddington luminosities of the most massive disrupting black holes (~2 x 10^8 M_\odot). Bolometric corrections to the spectra of the brightest flares are found to be large (~7.5 mag). Nevertheless, the brightest flares are likely to have absolute magnitudes in excess of -19 in V and -21 in U (in the absence of reddening). Because the spectra are so blue, K-corrections may actually brighten the flares in optical bands. If such flares are as frequent as believed, they may soon be detected in low or high redshift supernovae searches. The He II ionizing radiation produced in the flares may dominate that which is produced by all other sources in the centers of quiescent galaxies, creating a steady state, highly ionized, fossil nebula with an extent of ~1 kpc which may be observable in recombination lines.

We show that the disruption of a star by a 10^6 solar mass black hole in a galactic nucleus could under favorable circumstances produce an optically-thick envelope that radiates with a thermal spectrum at the Eddington limit, 10^{44} erg/s, for tens of years. The low apparent temperature of this envelope, 10^4 K, would be easily detectable in optical surveys. If most galaxies harbor a massive black hole at their center, then the Sloan Digital Sky Survey might find hundreds of galaxies with nuclear activity of this type. Because the envelope is driven to shine near the Eddington limit, a measurement of the source redshift and total luminosity could yield the black hole mass.

The two-point correlation function, $\xi$, of Lyman-alpha forest is found to be large, $\xi = 1.8^{+1.6}_{-1.2}$, > 90% confidence level, on the scale of 250-500 km/s for a sample of absorbers (0 < z < 1.3) assembled from HST Key Project Observations. This correlation function is stronger than at high redshift (z > 1.7) where $\xi \approx 0.2$ for velocities > 250 km/s.

By comparing neutrino fluxes and central temperatures calculated from 1000 detailed numerical solar models, we derive improved scaling laws which show how each of the neutrino fluxes depends upon the central temperature (flux $\propto T^m$); we also estimate uncertainties for the temperature exponents. With the aid of a one-zone model of the sun, we derive expressions for the temperature exponents of the neutrino fluxes. For the most important neutrino fluxes, the exponents calculated with the one-zone model agree to within 20\% or better with the exponents extracted from the detailed numerical models. The one-zone model provides a physical understanding of the temperature dependence of the neutrino fluxes. For the $pp$ neutrino flux, the one-zone model explains the (initially-surprising) dependence of the flux upon a negative power of the temperature and suggests a new functional dependence. This new function makes explicit the strong anti-correlation between the $^7$Be and $pp$ neutrino fluxes. The one-zone model also predicts successfully the average linear relations between neutrino fluxes, but cannot predict the appreciable scatter in a $\Delta \phi_i/\phi_i$ versus $\Delta \phi_j/\phi_j$ diagram.

We examine the width of the gamma-ray burst luminosity function through the distribution of GRB peak fluxes as detected by the Pioneer Venus Orbiter (PVO) and the Burst and Transient Source Experiment (BATSE). The strength of the analysis is greatly enhanced by using a merged catalog of peak fluxes from both instruments with good cross-calibration of their sensitivities. The range of peak fluxes is increased by approximately a factor of 20 relative to the BATSE catalog. Thus, more sensitive investigations of the $\log N-\log P$ distribution are possible. We place constraints on the width of the luminosity function of gamma-ray bursts brighter than the BATSE completeness limit by comparing the intensity distribution in the merged catalog with those produced by a variety of spatial density and luminosity functions. For the models examined, $90\%$ of the {\em detectable\/} bursts have peak luminosities within a range of 10, indicating that the peak luminosities of gamma-ray bursts span a markedly less wide range of values than many other of their measurable properties. We also discuss for which slopes of a power-law luminosity function the observed width is at the upper end of the constrained range. This is important in determining the power-law slopes for which luminosity-duration correlations could be important.

Bursts from soft gamma repeaters have been shown to be super-Eddington by a factor of 1000 and have been persuasively associated with compact objects. Here, a model of super-Eddington radiation transfer on the surface of a strongly magnetic ($\geq 10^{13}$ gauss) neutron star is studied and related to the observational constraints on soft gamma repeaters. In strong magnetic fields, the cross-section to electron scattering is strongly suppressed in one polarization state, so super-Eddington fluxes can be radiated while the plasma remains in hydrostatic equilibrium. The model offers a somewhat natural explanation for the observation of similarity between spectra from bursts of varying intensity. The radiation produced in the model is found to be linearly polarized to about 1 part in 1000 in a direction determined by the local magnetic field, and the large intensity variations between bursts are understood as a change in the radiating area on the source. Therefore, the polarization may vary as a function of burst intensity, since the complex structure of the magnetic field may be more apparent for larger radiating areas. It is shown that for radiation transfer calculations in this limit of super-strong magnetic fields it is sufficient to solve the radiation transfer equations for the low opacity state rather than the coupled equations for both. With this approximation, standard stellar atmosphere techniques are utilized to calculate the model energy spectrum.

We examine the width of the gamma-ray burst (GRB) luminosity function through the distribution of GRB peak count rates, C$_{\rm{peak}}$, as detected by BATSE (\cite{batse:93}). In the context of galactic corona spatial distribution models, we attempt to place constraints on the characteristic width of the luminosity function by comparing the observed intensity distribution with those produced by a range of density and luminosity functions. We find that the intrinsic width of the luminosity function cannot be very well restricted. However, the distribution of intrinsic luminosities of {\it detected bursts} can be limited: we find that most observed bursts have luminosities that are in a range of one to two decades, but a significant population of undetected less luminous bursts cannot be excluded. These findings demonstrate that the assumption that GRB are standard candles is sufficient but not necessary to explain the observed intensity distribution. We show that the main reason for the relatively poor constraints is the fact that the bright-end part of the GRB flux distribution is not yet sampled by BATSE, and better sampling in the future may lead to significantly stronger constraints on the width of the luminosity function.

Femtolensing is a gravitational lensing effect in which the magnification is a function not only of the positions and sizes of the source and lens, but also of the wavelength of light. Femtolensing is the only known effect of $(10^{-13}-10^{-16} M_{\sun})$ dark-matter objects and may possibly be detectable in cosmological gamma-ray burst spectra. We present a new and efficient algorithm for femtolensing calculations in general potentials. The physical-optics results presented here differ at low frequencies from the semi-classical approximation, in which the flux is attributed to a finite number of mutually coherent images. At higher frequencies, our results agree well with the semi-classical predictions. Applying our method to a point-mass lens with external shear, we find complex events that have structure at both large and small spectral resolution. In this way, we show that femtolensing may be observable for lenses up to $10^{-11}$ solar masses, much larger than previously believed. Additionally, we discuss the possibility of a search for femtolensing of white dwarfs in the LMC at optical wavelengths.

Soft Gamma Repeaters (SGRs) are a class of rare, high-energy galactic transients that have episodes of short (~0.1 sec), soft (~30 keV), intense (~100 Crab), gamma-ray bursts. We report an analysis of the x-ray emission from 95 SGR1806-20 events observed by the International Cometary Explorer. The spectral shape remains remarkably constant for bursts that differ in intensity by a range of 50. Below 15 keV the number spectrum falls off rapidly such that we can estimate the total intensity of the events. Assuming that SGR1806-20 is associated with the supernova remnant G10.0-0.3 (Kulkarni and Frail, Murakami \etal), the brightest events had a total luminosity of ~1.8 x 10^42 erg sec^-1, a factor of 2 x 10^4 above the Eddington limit. A third of the emission was above 30 keV. There are at least three processes that are consistent with the spectral rollover below 15 keV. (1)The rollover is consistent with some forms of self absorption. Typical thermal temperatures are ~20 keV and require an emitting surface with a radius between 10 and 50 km. The lack of spectral variability implies that only the size of the emitting surface varies between events. If the process is thermal synchrotron the required magnetic field might be too small to confine the plasma against the super Eddington flux. (2)The low energy rollover could be due to photoelectric absorption by ~10^24 Hydrogen atoms cm^-2 of neutral material with a cosmic abundance assuming a continuum similar to TB with T= ~22 keV. (3) Emission in the two lowest harmonics from a 1.3 x 10^12 Gauss field would appear as Doppler broadened lines and fall off rapidly below 15 keV.

Interactions of grazing incidence, ultra high energy cosmic rays with the earth's atmosphere may provide a new method of studying energetic cosmic rays with gamma-ray satellites. It is found that these cosmic ray interactions may produce gamma-rays on millisecond time scales which may be detectable by satellites. An extremely low gamma-ray background for transient gamma-ray events and a large area of interaction, the earth's surface, make the scheme plausible. The effective cross section of detection of interactions for cosmic rays above 10^20 eV is found to be more than two orders of magnitude higher than earth based detection techniques. This method may eventually offer an efficient way of probing this region of the cosmic ray energy spectrum where events are scarce. In this paper, a conceptual model is presented for the production of short bursts of gamma-rays based on these grazing incidence encounters with the earth's atmosphere.