Utilizing more than 100 long gamma-ray bursts (LGRBs) in the Swift-Ryan-2012 sample that include the observed redshifts and jet angles, Le & Mehta performed a timely study of the rate-density of LGRBs with an assumed broken power law GRB spectrum and obtained a GRB-burst-rate functional form that gives acceptable fits to the preSwift and Swift redshift and jet angle distributions. The results indicated an excess of LGRBs at redshift below $z \sim 2$ in the Swift sample. In this work, we are investigating if the excess is caused by the cosmological Hubble constant $H_0$, the $\gamma$-ray energy released $E_\gamma$, the low- and high-energy indices ($\alpha, \beta$) of the Band function, the minimum and maximum jet angles $\theta_{\rm j,min}$ and $\theta_{\rm j, max}$, or that the excess is due to a biased in the Swift-Ryan-2012 sample. Our analyses indicate that none of the above physical parameters resolved the excess problem, but suggesting that the Swift-Ryan-2012~sample is biased with possible afterglow selection effect. The following model physical parameter values provide the best fit to the Swift-Ryan-2012 and preSwift~samples: the Hubble constant $H_0 = 72 \, {\rm km s^{-1} Mpc^{-1}}$, the energy released $E_\gamma \sim 4.47 \times 10^{51}$ erg, the energy indices $\alpha \sim 0.9$ and $\beta \sim -2.13$, the jet angles of $\theta_{\rm j,max} \sim 0.8$ rad, and $\theta_{\rm j, min} \sim 0.065$ and $\sim 0.04$ rad for preSwift and Swift, respectively, $s \sim -1.55$ the jet angle power-law index, and a GRB formation rate that is similar to the Hopkins & Beacom observed star formation history and as extended by Li. Using the Swift Gamma-Ray Burst Host Galaxy Legacy Survey (SHOALS) Swift-Perley LGRB sample and applying the same physical parameter values as above, however, our model provides consistent results with this data set and indicating no excess of LGRBs at any redshift.

The determination of the neutrino mass is one of the major challenges in astroparticle physics today. Direct neutrino mass experiments, based solely on the kinematics of beta-decay, provide a largely model-independent probe to the neutrino mass scale. The Karlsruhe Tritium Neutrino (KATRIN) experiment is designed to directly measure the effective electron antineutrino mass with a sensitivity of 0.2 eV 90% CL. In this work we report on the first operation of KATRIN with tritium which took place in 2018. During this commissioning phase of the tritium circulation system, excellent agreement of the theoretical prediction with the recorded spectra was found and stable conditions over a time period of 13 days could be established. These results are an essential prerequisite for the subsequent neutrino mass measurements with KATRIN in 2019.

We report on the neutrino mass measurement result from the first four-week science run of the Karlsruhe Tritium Neutrino experiment KATRIN in spring 2019. Beta-decay electrons from a high-purity gaseous molecular tritium source are energy analyzed by a high-resolution MAC-E filter. A fit of the integrated electron spectrum over a narrow interval around the kinematic endpoint at 18.57 keV gives an effective neutrino mass square value of $(-1.0^{+0.9}_{-1.1})$ eV$^2$. From this we derive an upper limit of 1.1 eV (90$\%$ confidence level) on the absolute mass scale of neutrinos. This value coincides with the KATRIN sensitivity. It improves upon previous mass limits from kinematic measurements by almost a factor of two and provides model-independent input to cosmological studies of structure formation.

The Space Technology 7 Disturbance Reduction System (ST7-DRS) is a NASA technology demonstration payload that operated from January 2016 through July of 2017 on the European Space Agency's LISA Pathfinder spacecraft. The joint goal of the NASA and ESA missions was to validate key technologies for a future space-based gravitational wave observatory targeting the source-rich milliHertz band. The two primary components of ST7-DRS are a micropropulsion system based on colloidal micro-Newton thrusters (CMNTs) and a control system that simultaneously controls the attitude and position of the spacecraft and the two free-flying test masses (TMs). This paper presents our main experimental results and summarizes the overall the performance of the CMNTs and control laws. We find that the CMNT performance to be consistent with pre-flight predictions, with a measured system thrust noise on the order of $100\,\textrm{nN}/\sqrt{\textrm{Hz}}$ in the $1\,\textrm{mHz}\leq f \leq 30\,\textrm{mHz}$ band. The control system maintained the TM-spacecraft separation with an RMS error of less than 2$\,$nm and a noise spectral density of less than $3\,\textrm{nm}/\sqrt{\textrm{Hz}}$ in the same band. Thruster calibration measurements yield thrust values consistent with the performance model and ground-based thrust-stand measurements, to within a few percent. We also report a differential acceleration noise between the two test masses with a spectral density of roughly $3\,\textrm{fm}/\textrm{s}^2/\sqrt{\textrm{Hz}}$ in the $1\,\textrm{mHz}\leq f \leq 30\,\textrm{mHz}$ band, slightly less than twice as large as the best performance reported with the baseline LISA Pathfinder configuration and below the current requirements for the Laser Interferometer Space Antenna (LISA) mission.

A High-performance Atmospheric Radiation Package (HARP) is developed for studying multiple-scattering planetary atmospheres. HARP is an open-source program written in C++ that utilizes high-level data structure and parallel-computing algorithms. It is generic in three aspects. First, the construction of the model atmospheric profile is generic. The program can either take in an atmospheric profile or construct an adiabatic thermal and compositional profile, taking into account the clouds and latent heat release due to condensation. Second, the calculation of opacity is generic, based on line-by-line molecular transitions and tabulated continuum data, along with a table of correlated-k opacity provided as an option to speed up the calculation of energy fluxes. Third, the selection of the solver for the radiative transfer equation is generic. The solver is not hardwired in the program. Instead, based on the purpose, a variety of radiative transfer solvers can be chosen to couple with the atmosphere model and the opacity model. We use the program to investigate the radiative heating and cooling rates of all four giant planets in the Solar System. Our Jupiter's result is consistent with previous publications. Saturn has a nearly perfect balance between the heating rate and cooling rate. Uranus has the least radiative fluxes because of the lack of CH4 and its photochemical products. Both Uranus and Neptune suffer from a severe energy deficit in their stratospheres. Possible ways to resolve this issue are discussed. Finally, we recalculate the radiative time constants of all four giant planet atmospheres and find that the traditional values from (Conrath BJ, Gierasch PJ, Leroy SS. Temperature and Circulation in the Stratosphere of the Outer Planets. Icar. 1990;83:255-81) are significantly overestimated.

Using our theory for the production of relativistic outflows, we estimate the jet launching radius and the inferred mass accretion rate for 52 low-power radio-loud AGNs based on the observed jet powers. Our analysis indicates that (1) a significant fraction of the accreted energy is required to convert the accreted mass to relativistic energy particles for the production of the jets near the event horizon, (2) the jets launching radius moves radially toward the horizon as the mass accretion rate or jets power increases, and (3) no jet/outflow formation is possible beyond 44 gravitational radii.

Le & Dermer developed a gamma-ray burst (GRB) model to fit the redshift and the jet opening angle distributions measured with pre-Swift and Swift missions and showed that GRBs do not follow the star formation rate. Their fitted results were obtained without the opening angle distribution from Swift with an incomplete Swift sample, and that the calculated jet opening angle distribution was obtained by assuming a flat $\nu F_\nu$ spectrum. In this paper, we revisit the work done by Le & Dermer with an assumed broken power law GRB spectrum. Utilizing more than 100 GRBs in the Swift sample that include both the observed estimated redshifts and jet opening angles, we obtain a GRB burst rate functional form that gives acceptable fits to the pre-Swift and Swift redshift and jet opening angle distributions with an indication that an excess of GRBs exists at low redshift below $z \approx 2$. The mean redshifts and jet opening angles for pre-Swift~(Swift) are $\langle z \rangle \sim 1.3$ ($1.7$) and $\langle \theta_{\rm j} \rangle \sim 7^o$ ($11^o$), respectively. Assuming a GRB rate density (SFR9), similar to the Hopkins & Beacom star formation history and as extended by Li, the fraction of high-redshift GRBs is estimated to be below 10% and 5% at $z \geq 4$ and $z\geq 5$, respectively, and below 10% at $z \leq 1$.

Depending on the values of the energy and angular momentum per unit mass in the gas supplied at large radii, inviscid advection-dominated accretion flows can display velocity profiles with either preshock deceleration or preshock acceleration. Nakayama has shown that these two types of flow configurations are expected to have different stability properties. By employing the Chevalier \& Imamura linearization method and the Nakayama instability boundary conditions, we discover that there are regions of parameter space where disks/shocks with outflows can be stable or unstable. In regions of instability, we find that preshock deceleration is always unstable to the zeroth mode with zero frequency of oscillation, but is always stable to the fundamental mode and overtones. Furthermore, we also find that preshock acceleration is always unstable to the zeroth mode and that the fundamental mode and overtones become increasingly less stable as the shock location moves away from the horizon when the disk half-height expands above $\sim 12$ gravitational radii at the shock radius. In regions of stability, we demonstrate the zeroth mode to be stable for the velocity profiles that exhibit preshock acceleration and deceleration. Moreover, for models that are linearly unstable, our model suggests the possible existence of quasi-periodic oscillations (QPOs) with ratios 2:3 and 3:5. These ratios are believed to occur in stellar and supermassive black hole candidates, for example, in GRS 1915+105 and Sgr A*, respectively. We expect that similar QPO ratios also exist in regions of stable shocks.

We present evidence for ultraviolet/optical microlensing in the gravitationally lensed quasar Q0957+561. We combine new measurements from our optical monitoring campaign at the United States Naval Observatory, Flagstaff (USNO) with measurements from the literature and find that the time-delay-corrected r-band flux ratio m_A - m_B has increased by ~0.1 magnitudes over a period of five years beginning in the fall of 2005. We apply our Monte Carlo microlensing analysis procedure to the composite light curves, obtaining a measurement of the optical accretion disk size, log {(r_s/cm)[cos(i)/0.5]^{1/2}} = 16.2^{+0.5}_{-0.6}, that is consistent with the quasar accretion disk size - black hole mass relation.

We develop and discuss global accretion solutions for viscous ADAF disks containing centrifugally supported isothermal shock waves. The fact that such shocks can exist at all in ADAF disks is a new result. Interestingly, we find that isothermal shocks can form even when the level of viscous dissipation is relatively high. In order to better understand this phenomenon, we explore all possible combinations of the fundamental flow parameters, such as specific energy, specific angular momentum, and viscosity, to obtain the complete family of global solutions. This procedure allows us to identify the region of the parameter space where isothermal shocks can exist in viscous ADAF disks. The allowed region is maximized in the inviscid case, and it shrinks as the level of viscous dissipation increases. Adopting the canonical value gamma=1.5 for the ratio of specific heats, we find that the shock region disappears completely when the Shakura-Sunyaev viscosity parameter alpha exceeds the critical value ~0.27. This establishes for the first time that steady ADAF disks containing shocks can exist even for relatively high levels of viscous dissipation. If an isothermal shock is present in the disk, it would have important implications for the acceleration of energetic particles that can escape to power the relativistic jets commonly observed around underfed, radio-loud black holes. In two specific applications, we confirm that the kinetic luminosity lost from the disk at the isothermal shock location is sufficient to power the observed relativistic outflows in M87 and Sgr A*.

In this Letter, we present a new self-consistent theory for the production of the relativistic outflows observed from radio-loud black hole candidates and active galaxies as a result of particle acceleration in hot, viscous accretion disks containing standing, centrifugally-supported isothermal shocks. This is the first work to obtain the structure of such disks for a relatively large value of the Shakura-Sunyaev viscosity parameter ($\alpha=0.1$), and to consider the implications of the shock for the acceleration of relativistic particles in viscous disks. In our approach, the hydrodynamics and the particle acceleration are coupled and the solutions are obtained self-consistently based on a rigorous mathematical method. We find that particle acceleration in the vicinity of the shock can provide enough energy to power the observed relativistic jet in M87.

The distribution function describing the acceleration of relativistic particles in an advection-dominated accretion disk is analyzed using a transport formalism that includes first-order Fermi acceleration, advection, spatial diffusion, and the escape of particles through the upper and lower surfaces of the disk. When a centrifugally-supported shock is present in the disk, the concentrated particle acceleration occurring in the vicinity of the shock channels a significant fraction of the binding energy of the accreting gas into a population of relativistic particles. These high-energy particles diffuse vertically through the disk and escape, carrying away both energy and entropy and allowing the remaining gas to accrete. The dynamical structure of the disk/shock system is computed self-consistently using a model previously developed by the authors that successfully accounts for the production of the observed relativistic outflows (jets) in M87 and \SgrA. This ensures that the rate at which energy is carried away from the disk by the escaping relativistic particles is equal to the drop in the radial energy flux at the shock location, as required for energy conservation. We investigate the influence of advection, diffusion, and acceleration on the particle distribution by computing the nonthermal Green's function, which displays a relatively flat power-law tail at high energies. We also obtain the energy distribution for the particles escaping from the disk, and we conclude by discussing the spectrum of the observable secondary radiation produced by the escaping particles.

Results of a phenomenological model to estimate the GRB detection rate by the Fermi Gamma ray Space Telescope are reported. This estimate is based on the BATSE 4B GRB fluence distribution, the mean ratio of fluences measured at 100 MeV - 5 GeV with EGRET and at 20 keV - 2 MeV with BATSE, and the mean EGRET GRB spectrum for the 5 EGRET spark-chamber GRBs. For a 10% fluence ratio and a number spectral index alpha_1 = -2 at 100 MeV - 5 GeV energies, we estimate a rate of ~ 20 and 4 GRBs per yr in the Fermi Large Area Telescope field of view with at least 5 photons with energy E > 100 MeV and E > 1 GeV, respectively. We also estimate ~ 1.5 GRBs per yr in the Fermi FoV where at least 1 photon with energy E > 10 GeV is detected. For these parameters, we estimate = 1 - 2 GRBs per year detected with the Fermi telescope with more than 100 gamma rays with E >~ 100$ MeV. Comparison predictions for alpha_1 = -2.2, different fluence ratios, and the AGILE gamma-ray satellite are made. Searches for different classes of GRBs using a diagram plotting 100 MeV - 10 GeV fluence vs. 20 keV - 20 MeV fluence is considered as a way to search for separate classes of GRBs and, specifically, spectral differences between the short-hard and long duration GRB classes, and for hard components in GRBs.

Relativistic black-hole jet sources are leading candidates for high energy (>> TeV) neutrino production. The relations defining (a) efficient photopion losses of cosmic-ray protons on target photons and (b) gamma-gamma opacity of gamma rays through that same target photon field imply clear multiwavelength predictions for when and at what energies blazars and GRBs should be most neutrino bright and gamma-ray dim. The use of multiwavelength observations to test the standard relativistic jet model for these source is illustrated.

We summarize our model for long-duration gamma ray bursts (GRBs) that fits the redshift (z) distributions measured with Swift and missions before Swift, and the pre-Swift GRB jet opening-angle distribution inferred from achromatic breaks in the optical light curves. We find that the comoving rate density of GRB sources exhibits positive evolution to $z \gtrsim 3$ -- 5, whereas the star formation rate inferred from measurements of the blue and UV luminosity density peaks at $z\sim 1$ -- 3. The mean intrinsic beaming factor of GRBs is found to be $\approx 34$ -- 42, and we predict that the mean GRB optical jet opening half-angle measured with Swift is $\approx 10^\circ$. We estimate the number of GRBs per year that GLAST is expected to observe based on ratios of BATSE and EGRET GRB fluences.

A simple physical model for long-duration gamma ray bursts (GRBs) is used to fit the redshift (z) and the jet opening-angle distributions measured with earlier GRB missions and with Swift. The effect of different sensitivities for GRB triggering is sufficient to explain the difference in the z distributions of the pre-Swift and Swift samples, with mean redshifts of <z> ~1.5 and <z>~2.7, respectively. Assuming that the emission properties of GRBs do not change with time, we find that the data can only be fitted if the comoving rate-density of GRB sources exhibits positive evolution to z >~ 3-5. The mean intrinsic beaming factor of GRBs is found to range from ~34-42, with the Swift average opening half-angle <\theta_j> ~10 degree, compared to the pre-Swift average of <\theta_j> ~7 degree. Within the uniform jet model, the GRB luminosity function is proportional to L^{-3.25}_*, as inferred from our best fit to the opening angle distribution. Because of the unlikely detection of several GRBs with z <~ 0.25, our analysis indicates that low redshift GRBs represent a different population of GRBs than those detected at higher redshifts. Neglecting possible metallicity effects on GRB host galaxies, we find that ~1 GRB occurs every 600,000 yrs in a local L_* spiral galaxy like the Milky Way. The fraction of high-redshift GRBs is estimated at 8-12% and 2.5-6% at z >= 5 and z >= 7, respectively, assuming continued positive evolution of the GRB rate density to high redshifts.

Stochastic acceleration of charged particles due to interactions with magnetohydrodynamic (MHD) plasma waves is the dominant process leading to the formation of the high-energy electron and ion distributions in a variety of astrophysical systems. Collisions with the waves influence both the energization and the spatial transport of the particles, and therefore it is important to treat these two aspects of the problem in a self-consistent manner. We solve the representative Fokker-Planck equation to obtain a new, closed-form solution for the time-dependent Green's function describing the acceleration and escape of relativistic ions interacting with Alfven or fast-mode waves characterized by momentum diffusion coefficient $D(p)\propto p^q$ and mean particle escape timescale $t_esc(p) \propto p^{q-2}$, where $p$ is the particle momentum and $q$ is the power-law index of the MHD wave spectrum. In particular, we obtain solutions for the momentum distribution of the ions in the plasma and also for the momentum distribution of the escaping particles, which may form an energetic outflow. The general features of the solutions are illustrated via examples based on either a Kolmogorov or Kraichnan wave spectrum. The new expressions complement the results obtained by Park and Petrosian, who presented exact solutions for the hard-sphere scattering case ($q=2$) in addition to other scenarios in which the escape timescale has a power-law dependence on the momentum. Our results have direct relevance for models of high-energy radiation and cosmic-ray production in astrophysical environments such as $\gamma$-ray bursts, active galaxies, and magnetized coronae around black holes.

Relativistic outflows (jets) of matter are commonly observed from systems containing black holes. The strongest outflows occur in the radio-loud systems, in which the accretion disk is likely to have an advection-dominated structure. In these systems, it is clear that the binding energy of the accreting gas is emitted primarily in the form of particles rather than radiation. However, no comprehensive model for the disk structure and the associated outflows has yet been produced. In particular, none of the existing models establishes a direct physical connection between the presence of the outflows and the action of a microphysical acceleration mechanism operating in the disk. In this paper we explore the possibility that the relativistic protons powering the jet are accelerated at a standing, centrifugally-supported shock in the underlying accretion disk via the first-order Fermi mechanism. The theoretical analysis employed here parallels the early studies of cosmic-ray acceleration in supernova shock waves, and the particle acceleration and disk structure are treated in a coupled, self-consistent manner based on a rigorous mathematical approach. We find that first-order Fermi acceleration at standing shocks in advection-dominated disks proves to be a very efficient means for accelerating the jet particles. Using physical parameters appropriate for M87 and SgrA*, we verify that the jet kinetic luminosities computed using our model agree with estimates based on observations of the sources.

In this Letter, we suggest that the relativistic protons powering the outflows emanating from radio-loud systems containing black holes are accelerated at standing, centrifugally-supported shocks in hot, advection-dominated accretion disks. Such disks are ideal sites for first-order Fermi acceleration at shocks because the gas is tenuous, and consequently the mean free path for particle-particle collisions generally exceeds the thickness of the disk. The accelerated particles are therefore able to avoid thermalization, and as a result a small fraction of them achieve very high energies and escape from the disk. In our approach the hydrodynamics and the particle acceleration are coupled and the solutions are obtained self-consistently based on a rigorous mathematical treatment. The theoretical analysis of the particle transport parallels the early studies of cosmic-ray acceleration in supernova shock waves. We find that particle acceleration in the vicinity of the shock can extract enough energy to power a relativistic jet. Using physical parameters appropriate for M87 and Sgr A*, we confirm that the jet kinetic luminosities predicted by the theory agree with the observational estimates.

The structure of the inner region of an advection-dominated accretion disk around a nonrotating black hole is explored by applying asymptotic analysis in the region just outside the event horizon. We assume that the viscous transport is described by the standard Shakura-Sunyaev prescription throughout the disk, including the inner region close to the horizon. One of our goals is to explore the self-consistency of this assumption by analyzing the causality of the viscous transport near the black hole. The effects of general relativity are incorporated in an approximate manner by utilizing a pseudo-Newtonian gravitational potential. Analysis of the conservation equations yields unique asymptotic forms for the behaviors of the radial inflow velocity, the density, the sound speed, and the angular velocity. The specific behaviors are determined by three quantities; namely, the accreted specific energy, the accreted specific angular momentum, and the accreted specific entropy. The additional requirement of passage through a sonic point further constrains the problem, leaving only two free parameters. Our detailed results confirm that the Shakura-Sunyaev viscosity yields a well-behaved flow structure in the inner region that satisfies the causality constraint. We also show that the velocity distribution predicted by our pseudo-Newtonian model agrees with general relativity in the vicinity of the horizon. The asymptotic expressions we derive provide convenient boundary conditions for the development of global models via numerical integration of the conservation equations.

We present a spectral analysis of a 30d, near-continuous observation of the Seyfert 1 galaxy NGC 7469 with RXTE. Daily integrations show strong spectral changes during the observation. Ourmain result is that we find the X-ray spectral index to be correlated with the UV flux. Furthermore, the broad-band X-ray photon flux is also correlated with the UV continuum. These correlations point towards a model in which the X-rays originate via thermal Comptonization of UV seed photons. Furthermore, the UV is also correlated with the extrapolation of the X-ray power law into the soft X-ray/EUV region. Our data analysis therefore re-opens the possibility the the UV photons and their variability arise from reprocessing, as long as the primary source of heating is photo-electric absorption in the reprocessor, rather than Compton downscattering, A coherent picture of the X-ray/UV variability can therefore be constructed whereby absorption and reprocessing of EUV/soft X-rays in a standard accretion disk produce a variable seed photon distribution which are in turn up-scattered into the X-ray band. We also find a significant correlation between the 2-10 keV flux and the 6.4 keV iron K-line suggesting that at least some portion of the line originates within ~1 lt day of the X-ray continuum source. Neither the power law photon index nor the Compton reflection component are correlated with the 2-10 keV flux. The latter is not correlated with the iron K line flux either. We do find an apparent correlation between the X-ray spectral index and the strength of the Compton reflection component. In an Appendix we show, however, that this can be produced by a combination of statistical and systematic errors. We conclude the apparent variations in the Compton reflection component may be an artifact of these effects.