Black holes formed in the early universe, prior to the formation of stars, can exist as dark matter and also contribute to the black hole merger events observed in gravitational waves. We set a new limit on the abundance of primordial black holes (PBHs) by considering interactions of PBHs with the interstellar medium, which result in the heating of gas. We examine generic heating mechanisms, including emission from the accretion disk, dynamical friction, and disk outflows. Using the data from the Leo T dwarf galaxy, we set a new cosmology-independent limit on the abundance of PBHs in the mass range $\mathcal{O}(1) M_{\odot}-10^7 M_{\odot}$.

In this addendum to the article JCAP 1912 (2019) no.12, 047 (arXiv:1909.13328) on the cosmological dependence of non-resonantly produced sterile neutrinos we discuss, using an analytic treatment, the parameter regions of large active-sterile neutrino mixing angles where sterile neutrinos can approach thermalization. We show that these additional considerations affect only large active-sterile neutrino mixing already rejected by different limits. Hence, the allowed sterile neutrino parameter regions are unaffected.

The detection of a sterile neutrino could constitute the first observation of a particle that could have been produced before Big-Bang Nucleosynthesis (BBN), and could provide information about the yet untested pre-BBN era. The cosmological evolution in this era could be drastically different than typically assumed in what constitutes the standard cosmology, as happens in a variety of motivated particle models. In this work we assess the sensitivity to different pre-BBN cosmologies in which entropy is conserved of 0.01 eV to 1 MeV mass sterile neutrinos produced in the early Universe via resonant active-sterile oscillations, which requires a large lepton asymmetry. We identify mass ranges where it is possible to have two populations of the same sterile neutrino, one with a colder and one with a hotter momentum spectra, which is in principle an observable effect. Furthermore, we show the regions in mass and mixing where fully resonant production (i.e. simultaneously coherent and adiabatic) can occur. We find that in several of the cosmologies we consider, including the standard one, for a lepton asymmetry larger than $\sim10^{-4}$ fully resonantly produced sterile neutrinos in the eV-mass range can evade all cosmological constraints.

We discuss how a laboratory detection of a sterile neutrino not only would constitute a fundamental discovery of a new particle, but could also provide an indication of the evolution of the Universe before Big-Bang Nucleosynthesis (BBN), a fundamental discovery in cosmology. These "visible" sterile neutrinos could be detected in experiments such as KATRIN/TRISTAN and HUNTER in the keV mass range and PTOLEMY, KATRIN and reactor neutrino experiments in the eV mass range. Standard assumptions are usually made to compute the relic abundance and momentum distribution of particles produced before the temperature of the Universe was 5 MeV, an epoch from which there are no observed remnants thus far. However, non-standard pre-BBN cosmologies based on other assumptions that are equally in agreement with all existing data can arise in some theoretical models. We revisit the production of 0.01 eV to 1 MeV sterile neutrinos via non-resonant active-sterile flavor oscillations in several pre-BBN cosmologies. We give general equations for models in which the expansion of the Universe is parametrized by its amplitude and temperature power and where entropy is conserved, which include kination and scalar tensor models as special cases.

A laboratory detection of a sterile neutrino could provide the first indication of the evolution of the Universe before Big-Bang Nucleosynthesis (BBN), an epoch yet untested. Such "visible" sterile neutrinos are observable in upcoming experiments such as KATRIN/TRISTAN and HUNTER in the keV mass range and PTOLEMY and others in the eV mass range. A set of standard assumptions is typically made about cosmology before the temperature of the Universe was 5 MeV. However, non-standard pre-BBN cosmologies based on alternative assumptions could arise in motivated theoretical models and are equally in agreement with all existing data. We revisit the production of sterile neutrinos of mass 0.01 eV to 1 MeV in two examples of such models: scalar-tensor and low reheating temperature pre-BBN cosmologies. In both of them, the putative 3.5 keV X-ray signal line corresponds to a sterile neutrino with a mixing large enough to be tested in upcoming laboratory experiments. Additionally, the cosmological/astrophysical upper limits on active-sterile neutrino mixings are significantly weaker than in the standard pre-BBN cosmology, which shows that these limits are not robust. For example, in the scalar-tensor case the potential signal regions implied by the LSND and MiniBooNE short-baseline as well as the DANSS and NEOS reactor experiments are entirely not bound by cosmological restrictions. Our work highlights that sterile neutrinos may constitute a sensitive probe of the pre-BBN epoch.

In this work, BV RI light curves of 55 Type II supernovae (SNe II) from the Lick Observatory Supernova Search program obtained with the Katzman Automatic Imaging Telescope and the 1 m Nickel telescope from 2006 to 2018 are presented. Additionally, more than 150 spectra gathered with the 3 m Shane telescope are published. We conduct an analyse of the peak absolute magnitudes, decline rates, and time durations of different phases of the light and colour curves. Typically, our light curves are sampled with a median cadence of 5.5 days for a total of 5093 photometric points. In average V-band plateau declines with a rate of 1.29 mag (100 days)-1, which is consistent with previously published samples. For each band, the plateau slope correlates with the plateau length and the absolute peak magnitude: SNe II with steeper decline have shorter plateau duration and are brighter. A time-evolution analysis of spectral lines in term of velocities and pseudoequivalent widths is also presented in this paper. Our spectroscopic sample ranges between 1 and 200 days post-explosion and has a median ejecta expansion velocity at 50 days post-explosion of 6500 km/s (Halpha line) and a standard dispersion of 2000 km/s. Nebular spectra are in good agreement with theoretical models using a progenitor star having a mass <16 Msol. All the data are available to the community and will help to understand SN II diversity better, and therefore to improve their utility as cosmological distance indicators.

We present BVRI and unfiltered light curves of 93 Type Ia supernovae (SNe Ia) from the Lick Observatory Supernova Search (LOSS) follow-up program conducted between 2005 and 2018. Our sample consists of 78 spectroscopically normal SNe Ia, with the remainder divided between distinct subclasses (three SN 1991bg-like, three SN 1991T-like, four SNe Iax, two peculiar, and three super-Chandrasekhar events), and has a median redshift of 0.0192. The SNe in our sample have a median coverage of 16 photometric epochs at a cadence of 5.4 days, and the median first observed epoch is ~4.6 days before maximum B-band light. We describe how the SNe in our sample are discovered, observed, and processed, and we compare the results from our newly developed automated photometry pipeline to those from the previous processing pipeline used by LOSS. After investigating potential biases, we derive a final systematic uncertainty of 0.03 mag in BVRI for our dataset. We perform an analysis of our light curves with particular focus on using template fitting to measure the parameters that are useful in standardising SNe Ia as distance indicators. All of the data are available to the community, and we encourage future studies to incorporate our light curves in their analyses.

The Gamma-Ray Integrated Detectors (GRID) is a space mission concept dedicated to monitoring the transient gamma-ray sky in the energy range from 10 keV to 2 MeV using scintillation detectors onboard CubeSats in low Earth orbits. The primary targets of GRID are the gamma-ray bursts (GRBs) in the local universe. The scientific goal of GRID is, in synergy with ground-based gravitational wave (GW) detectors such as LIGO and VIRGO, to accumulate a sample of GRBs associated with the merger of two compact stars and study jets and related physics of those objects. It also involves observing and studying other gamma-ray transients such as long GRBs, soft gamma-ray repeaters, terrestrial gamma-ray flashes, and solar flares. With multiple CubeSats in various orbits, GRID is unaffected by the Earth occultation and serves as a full-time and all-sky monitor. Assuming a horizon of 200 Mpc for ground-based GW detectors, we expect to see a few associated GW-GRB events per year. With about 10 CubeSats in operation, GRID is capable of localizing a faint GRB like 170817A with a 90% error radius of about 10 degrees, through triangulation and flux modulation. GRID is proposed and developed by students, with considerable contribution from undergraduate students, and will remain operated as a student project in the future. The current GRID collaboration involves more than 20 institutes and keeps growing. On August 29th, the first GRID detector onboard a CubeSat was launched into a Sun-synchronous orbit and is currently under test.

Two separated sequences of blue straggler stars (BSSs) have been revealed by Ferraro et al. (2009) in the color-magnitude diagram (CMD) of the Milky Way globular cluster M30. Their presence has been suggested to be related to the two BSS formation channels (namely, collisions and mass-transfer in close binaries) operating within the same stellar system. The blue sequence was indeed found to be well reproduced by collisional BSS models. In contrast, no specific models for mass transfer BSSs were available for an old stellar system like M30. Here we present binary evolution models, including case-B mass transfer and binary merging, specifically calculated for this cluster. We discuss in detail the evolutionary track of a $0.9+0.5 M_\odot$ binary, which spends approximately 4 Gyr in the BSS region of the CMD of a 13 Gyr old cluster. We also run Monte-Carlo simulations to study the distribution of mass transfer BSSs in the CMD and to compare it with the observational data. Our results show that: (1) the color and magnitude distribution of synthetic mass transfer BSSs defines a strip in the CMD that nicely matches the observed red BSS sequence, thus providing strong support to the mass transfer origin for these stars; (2) the CMD distribution of synthetic BSSs never attains the observed location of the blue BSS sequence, thus reinforcing the hypothesis that the latter formed through a different channel (likely collisions); (3) most ($\sim 60\%$) of the synthetic BSSs are produced by mass-transfer models, while the remaining $< 40\%$ requires the contribution from merger models.

We present the results of a LIGO search for gravitational waves (GWs) associated with GRB 051103, a short-duration hard-spectrum gamma-ray burst (GRB) whose electromagnetically determined sky position is coincident with the spiral galaxy M81, which is 3.6 Mpc from Earth. Possible progenitors for short-hard GRBs include compact object mergers and soft gamma repeater (SGR) giant flares. A merger progenitor would produce a characteristic GW signal that should be detectable at the distance of M81, while GW emission from an SGR is not expected to be detectable at that distance. We found no evidence of a GW signal associated with GRB 051103. Assuming weakly beamed gamma-ray emission with a jet semi-angle of 30 deg we exclude a binary neutron star merger in M81 as the progenitor with a confidence of 98%. Neutron star-black hole mergers are excluded with > 99% confidence. If the event occurred in M81 our findings support the the hypothesis that GRB 051103 was due to an SGR giant flare, making it the most distant extragalactic magnetar observed to date.

The gravitational-wave (GW) sky may include nearby pointlike sources as well as astrophysical and cosmological stochastic backgrounds. Since the relative strength and angular distribution of the many possible sources of GWs are not well constrained, searches for GW signals must be performed in a model-independent way. To that end we perform two directional searches for persistent GWs using data from the LIGO S5 science run: one optimized for pointlike sources and one for arbitrary extended sources. The latter result is the first of its kind. Finding no evidence to support the detection of GWs, we present 90% confidence level (CL) upper-limit maps of GW strain power with typical values between 2-20x10^-50 strain^2 Hz^-1 and 5-35x10^-49 strain^2 Hz^-1 sr^-1 for pointlike and extended sources respectively. The limits on pointlike sources constitute a factor of 30 improvement over the previous best limits. We also set 90% CL limits on the narrow-band root-mean-square GW strain from interesting targets including Sco X-1, SN1987A and the Galactic Center as low as ~7x10^-25 in the most sensitive frequency range near 160 Hz. These limits are the most constraining to date and constitute a factor of 5 improvement over the previous best limits.

In this paper, a grid of the binary evolution models are calculated for the study of blue straggler (BS) population in intermediate age ($\log Age$=7.85-8.95) star clusters. The BS formation via mass transfer and merging is studied systematically using our models. Both Case A and B close binary evolutionary tracks are calculated in a large range of parameters. The results show that BSs formed via Case B are generally bluer and even more luminous than those produced by Case A. Furthermore, the larger range in orbital separations of Case B models provide a probability of producing more BSs than Case A. Based on the grid of models, several Monte-Carlo simulations of BS populations in the clusters in the age range are carried out. The results show that BSs formed via different channels populate different areas in color magnitude diagram(CMD). The locations of BSs in CMD for a number of clusters are compared to our simulations as well. In order to investigate the influence of mass transfer efficiency in the models and simulations, a set of models are also calculated by implementing a constant mass transfer efficiency, $\beta$=0.5 during Roche lobe overflow (Case A binary evolution excluded). The result shows BSs can be formed via mass transfer at any given age in both cases. However, the distributions of the BS populations on CMD are different.

We present direct upper limits on continuous gravitational wave emission from the Vela pulsar using data from the Virgo detector's second science run. These upper limits have been obtained using three independent methods that assume the gravitational wave emission follows the radio timing. Two of the methods produce frequentist upper limits for an assumed known orientation of the star's spin axis and value of the wave polarization angle of, respectively, $1.9\ee{-24}$ and $2.2\ee{-24}$, with 95% confidence. The third method, under the same hypothesis, produces a Bayesian upper limit of $2.1\ee{-24}$, with 95% degree of belief. These limits are below the indirect {\it spin-down limit} of $3.3\ee{-24}$ for the Vela pulsar, defined by the energy loss rate inferred from observed decrease in Vela's spin frequency, and correspond to a limit on the star ellipticity of $\sim 10^{-3}$. Slightly less stringent results, but still well below the spin-down limit, are obtained assuming the star's spin axis inclination and the wave polarization angles are unknown.

Soft gamma repeaters (SGRs) and anomalous X-ray pulsars (AXPs) are thought to be magnetars: neutron stars powered by extreme magnetic fields. These rare objects are characterized by repeated and sometimes spectacular gamma-ray bursts. The burst mechanism might involve crustal fractures and excitation of non-radial modes which would emit gravitational waves (GWs). We present the results of a search for GW bursts from six galactic magnetars that is sensitive to neutron star f-modes, thought to be the most efficient GW emitting oscillatory modes in compact stars. One of them, SGR 0501+4516, is likely ~1 kpc from Earth, an order of magnitude closer than magnetars targeted in previous GW searches. A second, AXP 1E 1547.0-5408, gave a burst with an estimated isotropic energy >10^{44} erg which is comparable to the giant flares. We find no evidence of GWs associated with a sample of 1279 electromagnetic triggers from six magnetars occurring between November 2006 and June 2009, in GW data from the LIGO, Virgo, and GEO600 detectors. Our lowest model-dependent GW emission energy upper limits for band- and time-limited white noise bursts in the detector sensitive band, and for f-mode ringdowns (at 1090 Hz), are 3.0x10^{44} d_1^2 erg and 1.4x10^{47} d_1^2 erg respectively, where d_1 = d_{0501} / 1 kpc and d_{0501} is the distance to SGR 0501+4516. These limits on GW emission from f-modes are an order of magnitude lower than any previous, and approach the range of electromagnetic energies seen in SGR giant flares for the first time.

The physical mechanisms responsible for pulsar timing glitches are thought to excite quasi-normal mode oscillations in their parent neutron star that couple to gravitational wave emission. In August 2006, a timing glitch was observed in the radio emission of PSR B0833-45, the Vela pulsar. At the time of the glitch, the two co-located Hanford gravitational wave detectors of the Laser Interferometer Gravitational-wave observatory (LIGO) were operational and taking data as part of the fifth LIGO science run (S5). We present the first direct search for the gravitational wave emission associated with oscillations of the fundamental quadrupole mode excited by a pulsar timing glitch. No gravitational wave detection candidate was found. We place Bayesian 90% confidence upper limits of 6.3e-21 to 1.4e-20 on the peak intrinsic strain amplitude of gravitational wave ring-down signals, depending on which spherical harmonic mode is excited. The corresponding range of energy upper limits is 5.0e44 to 1.3e45 erg.

We present a search for periodic gravitational waves from the neutron star in the supernova remnant Cassiopeia A. The search coherently analyzes data in a 12-day interval taken from the fifth science run of the Laser Interferometer Gravitational-Wave Observatory. It searches gravitational wave frequencies from 100 to 300 Hz, and covers a wide range of first and second frequency derivatives appropriate for the age of the remnant and for different spin-down mechanisms. No gravitational wave signal was detected. Within the range of search frequencies, we set 95% confidence upper limits of 0.7--1.2e-24 on the intrinsic gravitational wave strain, 0.4--4e-4 on the equatorial ellipticity of the neutron star, and 0.005--0.14 on the amplitude of r-mode oscillations of the neutron star. These direct upper limits beat indirect limits derived from energy conservation and enter the range of theoretical predictions involving crystalline exotic matter or runaway r-modes. This is the first gravitational wave search to present upper limits on r-modes.

We present the results of a search for gravitational-wave bursts associated with 137 gamma-ray bursts (GRBs) that were detected by satellite-based gamma-ray experiments during the fifth LIGO science run and first Virgo science run. The data used in this analysis were collected from 2005 November 4 to 2007 October 1, and most of the GRB triggers were from the Swift satellite. The search uses a coherent network analysis method that takes into account the different locations and orientations of the interferometers at the three LIGO-Virgo sites. We find no evidence for gravitational-wave burst signals associated with this sample of GRBs. Using simulated short-duration (<1 s) waveforms, we set upper limits on the amplitude of gravitational waves associated with each GRB. We also place lower bounds on the distance to each GRB under the assumption of a fixed energy emission in gravitational waves, with typical limits of D ~ 15 Mpc (E_GW^iso / 0.01 M_o c^2)^1/2 for emission at frequencies around 150 Hz, where the LIGO-Virgo detector network has best sensitivity. We present astrophysical interpretations and implications of these results, and prospects for corresponding searches during future LIGO-Virgo runs.

We present an up-to-date, comprehensive summary of the rates for all types of compact binary coalescence sources detectable by the Initial and Advanced versions of the ground-based gravitational-wave detectors LIGO and Virgo. Astrophysical estimates for compact-binary coalescence rates depend on a number of assumptions and unknown model parameters, and are still uncertain. The most confident among these estimates are the rate predictions for coalescing binary neutron stars which are based on extrapolations from observed binary pulsars in our Galaxy. These yield a likely coalescence rate of 100 per Myr per Milky Way Equivalent Galaxy (MWEG), although the rate could plausibly range from 1 per Myr per MWEG to 1000 per Myr per MWEG. We convert coalescence rates into detection rates based on data from the LIGO S5 and Virgo VSR2 science runs and projected sensitivities for our Advanced detectors. Using the detector sensitivities derived from these data, we find a likely detection rate of 0.02 per year for Initial LIGO-Virgo interferometers, with a plausible range between 0.0002 and 0.2 per year. The likely binary neutron-star detection rate for the Advanced LIGO-Virgo network increases to 40 events per year, with a range between 0.4 and 400 per year.

Progenitor scenarios for short gamma-ray bursts (short GRBs) include coalescenses of two neutron stars or a neutron star and black hole, which would necessarily be accompanied by the emission of strong gravitational waves. We present a search for these known gravitational-wave signatures in temporal and directional coincidence with 22 GRBs that had sufficient gravitational-wave data available in multiple instruments during LIGO's fifth science run, S5, and Virgo's first science run, VSR1. We find no statistically significant gravitational-wave candidates within a [-5, +1) s window around the trigger time of any GRB. Using the Wilcoxon-Mann-Whitney U test, we find no evidence for an excess of weak gravitational-wave signals in our sample of GRBs. We exclude neutron star-black hole progenitors to a median 90% CL exclusion distance of 6.7 Mpc.

We present a search for gravitational waves from 116 known millisecond and young pulsars using data from the fifth science run of the LIGO detectors. For this search ephemerides overlapping the run period were obtained for all pulsars using radio and X-ray observations. We demonstrate an updated search method that allows for small uncertainties in the pulsar phase parameters to be included in the search. We report no signal detection from any of the targets and therefore interpret our results as upper limits on the gravitational wave signal strength. The most interesting limits are those for young pulsars. We present updated limits on gravitational radiation from the Crab pulsar, where the measured limit is now a factor of seven below the spin-down limit. This limits the power radiated via gravitational waves to be less than ~2% of the available spin-down power. For the X-ray pulsar J0537-6910 we reach the spin-down limit under the assumption that any gravitational wave signal from it stays phase locked to the X-ray pulses over timing glitches, and for pulsars J1913+1011 and J1952+3252 we are only a factor of a few above the spin-down limit. Of the recycled millisecond pulsars several of the measured upper limits are only about an order of magnitude above their spin-down limits. For these our best (lowest) upper limit on gravitational wave amplitude is 2.3x10^-26 for J1603-7202 and our best (lowest) limit on the inferred pulsar ellipticity is 7.0x10^-8 for J2124-3358.

We present direct upper limits on gravitational wave emission from the Crab pulsar using data from the first nine months of the fifth science run of the Laser Interferometer Gravitational-wave Observatory (LIGO). These limits are based on two searches. In the first we assume that the gravitational wave emission follows the observed radio timing, giving an upper limit on gravitational wave emission that beats indirect limits inferred from the spin-down and braking index of the pulsar and the energetics of the nebula. In the second we allow for a small mismatch between the gravitational and radio signal frequencies and interpret our results in the context of two possible gravitational wave emission mechanisms.

We present nine color CCD intermediate-band spectrophotometry of a two square degree field centered on the old open cluster M67, from 3890$\rm \AA$ to nearly 1$\mu$. These observations are taken as a part of the BATC (Beijing-Arizona-Taipei-Connecticut) Color Survey of the Sky, for both scientific and calibration reasons. With these data we show that the BATC survey can reach its goal of obtaining spectrophotometry to a zero point accuracy of 0.01 mag, and down to V = 21 with 0.3 mag random error. We fit the color-magnitude diagrams (CMDs) with Worthey's theoretical models. The net result is the excellent fit of the 4.0 Gyr, [Fe/H] = $-0.10$ model to our data, including a good fit to the main sequence (MS) turn-off. Our data are consistent with a toy model with 50\% of the stars in M67 being binaries and a random distribution of binary mass-ratios, although other models with different mass-ratio distributions cannot be ruled out. The spatial distribution and mass function (MF) of stars in M67 show marked effects of dynamical evolution and evaporation of stars from the cluster. Blue stragglers and binary stars are the most condensed within the cluster, with degree of condensation depending on mass.We find M67 to have an elongated shape, oriented at an angle of $15^{\circ}$ relative to the galactic plane. Within its tidal radius, the observed MF of M67 between 1.2 $\rm M_\odot$ and $\rm 0.8 M_\odot$ has a Salpeter slope $\rm \eta = -1.93 \pm 0.66$. For stars of mass below 0.8 $\rm M_\odot$, $\rm \eta \sim 0$. It is plausible that the leveling-off of the MF at lower masses is a result of evaporation of lower mass stars in this mass range at a rate of one every $\sim 10^7$ years. If so, it is plausible that the IMF of M67 has the canonical field value of $\rm \eta = -2.0$.