Recent studies have shown that the passage of a massive satellite through the disk of a spiral galaxy can induce vertical wobbles in the disk and produce features such as in-plane rings and phase-space spirals. Here we analyze a high-resolution N-body simulation of a live stellar disk perturbed by the recent passage of a massive dwarf galaxy that induces such perturbations. We study the implications of the phase-space structures for the estimate of the matter density through traditional Jeans modelling. The dwarf satellite excites rapid time-variations in the potential, leading to a significant bias of the local matter surface density determined through such a method. In particular, while the Jeans modelling gives reasonable estimates in the most overdense regions of the disk, we show that it tends to overestimate the dynamical surface density in underdense regions. In these regions, the phase-space spiral is indeed more marked in the surface of section of height vs. vertical velocity. This prediction can be verified with future Gaia data releases. Our finding is highly relevant for future attempts at determining the dynamical surface density of the outer Milky Way disk as a function of radius. The outer disk of the Milky Way is indeed heavily perturbed, and \textit{Gaia} DR2 data have clearly shown that such phase-space perturbations are even present locally. While our results show that traditional Jeans modelling should give reliable results in overdense regions of the disk, the important biases in underdense regions call for the development of non-equilibrium methods to estimate the dynamical matter density locally and in the outer disk.

We reconstruct posterior distributions for the position (sky area and distance) of a simulated set of binary neutron-star gravitational-waves signals observed with Advanced LIGO and Advanced Virgo. We use a Dirichlet Process Gaussian-mixture model, a fully Bayesian non-parametric method that can be used to estimate probability density functions with a flexible set of assumptions. The ability to reliably reconstruct the source position is important for multimessenger astronomy, as recently demonstrated with GW170817. We show that for detector networks comparable to the early operation of Advanced LIGO and Advanced Virgo, typical localization volumes are $\sim10^4$--$10^5~\mathrm{Mpc^3}$ corresponding to $\sim10^2$--$10^3$ potential host galaxies. The localization volume is a strong function of the network signal-to-noise ratio, scaling roughly $\propto \varrho_{net}^{-6}$. Fractional localizations improve with the addition of further detectors to the network. Our Dirichlet Process Gaussian-mixture model can be adopted for localizing events detected during future gravitational-wave observing runs, and used to facilitate prompt multimessenger follow-up.

We study the role of cold gas in quenching star formation in the green valley by analysing ALMA $^{12}$CO (1-0) observations of three galaxies with resolved optical spectroscopy from the MaNGA survey. We present resolution-matched maps of the star formation rate and molecular gas mass. These data are used to calculate the star formation efficiency (SFE) and gas fraction ($f_{\rm~gas}$) for these galaxies separately in the central `bulge' regions and outer disks. We find that, for the two galaxies whose global specific star formation rate (sSFR) deviates most from the star formation main sequence, the gas fraction in the bulges is significantly lower than that in their disks, supporting an `inside-out' model of galaxy quenching. For the two galaxies where SFE can be reliably determined in the central regions, the bulges and disks share similar SFEs. This suggests that a decline in $f_{\rm~gas}$ is the main driver of lowered sSFR in bulges compared to disks in green valley galaxies. Within the disks, there exist common correlations between the sSFR and SFE and between sSFR and $f_{\rm~gas}$ on kpc scales -- the local SFE or $f_{\rm~gas}$ in the disks declines with local sSFR. Our results support a picture in which the sSFR in bulges is primarily controlled by $f_{\rm~gas}$, whereas both SFE and $f_{\rm~gas}$ play a role in lowering the sSFR in disks. A larger sample is required to confirm if the trend established in this work is representative of green valley as a whole.

The Astropy Project (http://astropy.org) is, in its own words, "a community effort to develop a single core package for Astronomy in Python and foster interoperability between Python astronomy packages." For five years this project has been managed, written, and operated as a grassroots, self-organized, almost entirely volunteer effort while the software is used by the majority of the astronomical community. Despite this, the project has always been and remains to this day effectively unfunded. Further, contributors receive little or no formal recognition for creating and supporting what is now critical software. This paper explores the problem in detail, outlines possible solutions to correct this, and presents a few suggestions on how to address the sustainability of general purpose astronomical software.

The modern merger hypothesis offers a method of forming a new elliptical galaxy through merging two equal-mass, gas-rich disk galaxies fuelling a nuclear starburst followed by efficient quenching and dynamical stabilization. A key prediction of this scenario is a central concentration of young stars during the brief phase of morphological transformation from highly-disturbed remnant to new elliptical galaxy. To test this aspect of the merger hypothesis, we use integral field spectroscopy to track the stellar Balmer absorption and 4000\AA\ break strength indices as a function of galactic radius for 12 massive (${\rm M_{*}}\ge10^{10}{\rm M_{\odot}}$), nearby (${\rm z}\le0.03$), visually-selected plausible new ellipticals with blue-cloud optical colours and varying degrees of morphological peculiarities. We find that these index values and their radial dependence correlate with specific morphological features such that the most disturbed galaxies have the smallest 4000\AA\ break strengths and the largest Balmer absorption values. Overall, two-thirds of our sample are inconsistent with the predictions of the modern merger hypothesis. Of these eight, half exhibit signatures consistent with recent minor merger interactions. The other half have star formation histories similar to local, quiescent early-type galaxies. Of the remaining four galaxies, three have the strong morphological disturbances and star-forming optical colours consistent with being remnants of recent, gas-rich major mergers, but exhibit a weak, central burst consistent with forming $\sim5\%$ of their stars. The final galaxy possesses spectroscopic signatures of a strong, centrally-concentrated starburst and quiescent core optical colours indicative of recent quenching (i.e., a post-starburst signature) as prescribed by the modern merger hypothesis.

We use the Sloan Digital Sky Survey to investigate the properties of massive elliptical galaxies in the local Universe (z\leq0.08) that have unusually blue optical colors. Through careful inspection, we distinguish elliptical from non-elliptical morphologies among a large sample of similarly blue galaxies with high central light concentrations (c_r\geq2.6). These blue ellipticals comprise 3.7 per cent of all c_r\geq2.6 galaxies with stellar masses between 10^10 and 10^11 h^{-2} {\rm M}_{\sun}. Using published fiber spectra diagnostics, we identify a unique subset of 172 non-star-forming ellipticals with distinctly blue urz colors and young (< 3 Gyr) light-weighted stellar ages. These recently quenched ellipticals (RQEs) have a number density of 2.7-4.7\times 10^{-5}\,h^3\,{\rm Mpc}^{-3} and sufficient numbers above 2.5\times10^{10} h^{-2} {\rm M}_{\sun} to account for more than half of the expected quiescent growth at late cosmic time assuming this phase lasts 0.5 Gyr. RQEs have properties that are consistent with a recent merger origin (i.e., they are strong `first-generation' elliptical candidates), yet few involved a starburst strong enough to produce an E+A signature. The preferred environment of RQEs (90 per cent reside at the centers of < 3\times 10^{12}\,h^{-1}{\rm M}_{\sun} groups) agrees well with the `small group scale' predicted for maximally efficient spiral merging onto their halo center and rules out satellite-specific quenching processes. The high incidence of Seyfert and LINER activity in RQEs and their plausible descendents may heat the atmospheres of small host halos sufficiently to maintain quenching.

This work investigates the effect of the aperture size on derived galaxy properties for which we have spatially-resolved optical spectra. We focus on some indicators of star formation activity and dust attenuation for spiral galaxies that have been widely used in previous work on galaxy evolution. We have used 104 spiral galaxies from the CALIFA survey for which 2D spectroscopy with complete spatial coverage is available. From the 3D cubes we have derived growth curves of the most conspicuous Balmer emission lines (Halpha, Hbeta) for circular apertures of different radii centered at the galaxy's nucleus after removing the underlying stellar continuum. We find that the Halpha flux (f(Halpha)) growth curve follows a well defined sequence with aperture radius showing low dispersion around the median value. From this analysis, we derive aperture corrections for galaxies in different magnitude and redshift intervals. Once stellar absorption is properly accounted for, the f(Halpha)/f(Hbeta) ratio growth curve shows a smooth decline, pointing towards the absence of differential dust attenuation as a function of radius. Aperture corrections as a function of the radius are provided in the interval [0.3,2.5]R_50. Finally, the Halpha equivalent width (EW(Halpha)) growth curve increases with the size of the aperture and shows a very large dispersion for small apertures. This large dispersion prevents the use of reliable aperture corrections for this quantity. In addition, this result suggests that separating star-forming and quiescent galaxies based on observed EW(Halpha) through small apertures is likely to result in low EW(Halpha) star-forming galaxies begin classified as quiescent.

We present the first public data release of the CALIFA survey. It consists of science-grade optical datacubes for the first 100 of eventually 600 nearby (0.005<z<0.03) galaxies, obtained with the integral-field spectrograph PMAS/PPak mounted on the 3.5m telescope at the Calar Alto observatory. The galaxies in DR1 already cover a wide range of properties in color-magnitude space, morphological type, stellar mass, and gas ionization conditions. This offers the potential to tackle a variety of open questions in galaxy evolution using spatially resolved spectroscopy. Two different spectral setups are available for each galaxy, (i) a low-resolution V500 setup covering the nominal wavelength range 3745-7500A with a spectral resolution of 6.0A (FWHM), and (ii) a medium-resolution V1200 setup covering the nominal wavelength range 3650-4840A with a spectral resolution of 2.3A (FWHM). We present the characteristics and data structure of the CALIFA datasets that should be taken into account for scientific exploitation of the data, in particular the effects of vignetting, bad pixels and spatially correlated noise. The data quality test for all 100 galaxies showed that we reach a median limiting continuum sensitivity of 1.0x10^-18erg/s/cm^2/A/arcsec^2 at 5635A and 2.2x10^-18erg/s/cm^2/A/arcsec^2 at 4500A for the V500 and V1200 setup respectively, which corresponds to limiting r and g band surface brightnesses of 23.6mag/arcsec^2 and 23.4mag/arcsec^2, or an unresolved emission-line flux detection limit of roughly 1x10^-17erg/s/cm^2/arcsec^2 and 0.6x10^-17erg/s/cm^2/arcsec^2, respectively. The median spatial resolution is 3.7", and the absolute spectrophotometric calibration is better than 15% (1sigma). We also describe the available interfaces and tools that allow easy access to this first public CALIFA data.

We present here the Calar Alto Legacy Integral Field Area (CALIFA) survey, which has been designed to provide a first step in this direction.We summarize the survey goals and design, including sample selection and observational strategy.We also showcase the data taken during the first observing runs (June/July 2010) and outline the reduction pipeline, quality control schemes and general characteristics of the reduced data. This survey is obtaining spatially resolved spectroscopic information of a diameter selected sample of $\sim600$ galaxies in the Local Universe (0.005< z <0.03). CALIFA has been designed to allow the building of two-dimensional maps of the following quantities: (a) stellar populations: ages and metallicities; (b) ionized gas: distribution, excitation mechanism and chemical abundances; and (c) kinematic properties: both from stellar and ionized gas components. CALIFA uses the PPAK Integral Field Unit (IFU), with a hexagonal field-of-view of $\sim1.3\sq\arcmin'$, with a 100% covering factor by adopting a three-pointing dithering scheme. The optical wavelength range is covered from 3700 to 7000 {\AA}, using two overlapping setups (V500 and V1200), with different resolutions: R\sim850 and R\sim1650, respectively. CALIFA is a legacy survey, intended for the community. The reduced data will be released, once the quality has been guaranteed. The analyzed data fulfill the expectations of the original observing proposal, on the basis of a set of quality checks and exploratory analysis. We conclude from this first look at the data that CALIFA will be an important resource for archaeological studies of galaxies in the Local Universe.

We present the results of indirect searches for Weakly Interacting Massive Particles (WIMPs) with 1679.6 live days of data from the Super-Kamiokande detector using neutrino-induced upward through-going muons. The search is performed by looking for an excess of high energy muon neutrinos from WIMP annihilations in the Sun, the core of the Earth, and the Galactic Center, as compared to the number expected from the atmospheric neutrino background. No statistically significant excess was seen. We calculate flux limits in various angular cones around each of the above celestial objects. We obtain conservative model-independent upper limits on WIMP-nucleon cross-section as a function of WIMP mass and compare these results with the corresponding results from direct dark matter detection experiments.

The Milagro telescope monitors the northern sky for 100 GeV to 100 TeV transient emission through continuous very high energy wide-field observations. The large effective area and ~100 GeV energy threshold of Milagro allow it to detect very high energy (VHE) gamma-ray burst emission with much higher sensitivity than previous instruments and a fluence sensitivity at VHE energies comparable to that of dedicated gamma-ray burst satellites at keV to MeV energies. Even in the absence of a positive detection, VHE observations can place important constraints on gamma-ray burst (GRB) progenitor and emission models. We present limits on the VHE flux of 40 s -- 3 h duration transients nearby to earth, as well as sensitivity distributions which have been corrected for gamma-ray absorption by extragalactic background light and cosmological effects. The sensitivity distributions suggest that the typical intrinsic VHE fluence of GRBs is similar or weaker than the keV -- MeV emission, and we demonstrate how these sensitivity distributions may be used to place observational constraints on the absolute VHE luminosity of gamma-ray bursts for any GRB emission and progenitor model.

In this paper, after a discussion of general properties of statistical tests, we present the construction of the most powerful hypothesis test for determining the existence of a new phenomenon in counting-type experiments where the observed Poisson process is subject to a Poisson distributed background with unknown mean.

In this paper we discuss several methods of significance calculation and point out the limits of their applicability. We then introduce a self consistent scheme for source detection and discuss some of its properties. The method allows incorporating background anisotropies by vetoing existing small scale regions on the sky and compensating for known large scale anisotropies. By giving an example using Milagro gamma ray observatory we demonstrate how the method can be employed to relax the detector stability assumption. Two practical implementations of the method are discussed. The method is universal and can be used with any large field-of-view detector, where the object of investigation, steady or transient, point or extended, traverses its field of view.

Evidence of TeV emission from GRB970417a has been previously reported using data from the Milagrito detector. Constraints on the TeV fluence and the energy spectrum are now derived using additional data from a scaler system that recorded the rate of signals from the Milagrito photomultipliers. This analysis shows that if emission from GRB970417a has been observed, it must contain photons with energies above 650 GeV. Some consequences of this observation are discussed.

Milagrito, a detector sensitive to very high energy gamma rays, monitored the northern sky from February 1997 through May 1998. With a large field of view and a high duty cycle, this instrument was well suited to perform a search for TeV gamma-ray bursts (GRBs). We report on a search made for TeV counterparts to GRBs observed by BATSE. BATSE detected 54 GRBs within the field of view of Milagrito during this period. An excess of events coincident in time and space with one of these bursts, GRB 970417a, was observed by Milagrito. The excess has a chance probability of $2.8 \times 10^{-5}$ of being a fluctuation of the background. The probability for observing an excess at least this large from any of the 54 bursts is $1.5 \times 10^{-3}$. No significant correlations were detected from the other bursts.

Milagrito, a detector sensitive to gamma-rays at TeV energies, monitored the northern sky during the period February 1997 through May 1998. With a large field of view and high duty cycle, this instrument was used to perform a search for TeV counterparts to gamma-ray bursts. Within the Milagrito field of view 54 gamma-ray bursts at keV energies were observed by the Burst And Transient Satellite Experiment (BATSE) aboard the Compton Gamma-Ray Observatory. This paper describes the results of a preliminary analysis to search for TeV emission correlated with BATSE detected bursts. Milagrito detected an excess of events coincident both spatially and temporally with GRB 970417a, with chance probability $2.8 \times 10^{-5}$ within the BATSE error radius. No other significant correlations were detected. Since 54 bursts were examined the chance probability of observing an excess with this significance in any of these bursts is $1.5 \times 10^{-3}$. The statistical aspects and physical implications of this result are discussed.

Previous works have shown that water Cherenkov detectors have superior sensitivity to those of scintillation counters as applied to detecting extensive air showers (EAS). This is in large part due to their much higher sensitivity to EAS photons which are more than five times more numerous than EAS electrons. Large area water Cherenkov detectors can be constructed relatively cheaply and operated reliably. A sparse detector array has been designed which uses these types of detectors to substantially increase the area over which the Milagro Gamma Ray Observatory collects EAS information. Improvements to the Milagro detector's performance characteristics and sensitivity derived from this array and preliminary results from a prototype array currently installed near the Milagro detector will be presented.

The Milagro detector is an air shower array which uses the water Cherenkov technique and is capable of continuously monitoring the sky at energies near 1 TeV. The detector consists of 20000 metric tons of pure water instrumented with 723 photo-multiplier tubes (PMTs). The PMTs are arranged in a two-layer structure on a lattice of 3 m spacing covering 5000 $m^2$ area. The direction of the shower is determined from the relative timing of the PMT signals, necessitating a common time reference and amplitude slewing corrections to improve the time resolution. The calibration system to provide these consists of a pulsed laser driving 30 diffusing light sources deployed in the pond to allow cross-calibration of the PMTs. The system is capable of calibrating times and the pulse-heights from the PMTs using the time-over-threshold technique. The absolute energy scale is provided using single muons passing through the detector. The description of the calibration system of the Milagro detector and its prototype Milagrito will be presented.

Milagrito, a prototype for the Milagro detector, operated for 15 months in 1997-8 and collected 8.9 billion events. It was the first extensive air shower (EAS) array sensitive to showers intiated by primaries with energy below 1 TeV. The shadows of the sun and moon observed with cosmic rays can be used to study systematic pointing shifts and measure the angular resolution of EAS arrays. Below a few TeV, the paths of cosmic rays coming toward the earth are bent by the helio- and geo-magnetic fields. This is expected to distort and displace the shadows of the sun and the moon. The moon shadow, offset from the nominal (undeflected) position, has been observed with high statistical significance in Milagrito. This can be used to establish energy calibrations, as well as to search for the anti-matter content of the VHE cosmic ray flux. The shadow of the sun has also been observed with high significance.

Solar Energetic Particles from the 6 November 1997 solar flare/CME(coronal mass ejection) with energies exceeding 10 GeV have been detected by Milagrito, a prototype of the Milagro Gamma Ray Observatory. While particle acceleration beyond 1 GeV at the Sun is well established, few data exist for protons or ions beyond 10 GeV. The Milagro observatory, a ground based water Cherenkov detector designed for observing very high energy gamma ray sources, can also be used to study the Sun. Milagrito, which operated for approximately one year in 1997/98, was sensitive to solar proton and neutron fluxes above ~5- 10 GeV. Milagrito operated in a scaler mode, which was primarily sensitive to muons, low energy photons, and electrons, and the detector operated in a mode sensitive to showers and high zenith angle muons. In its scaler mode, Milagrito registered a rate increase coincident with the 6 November 1997 ground level event observed by Climax and other neutron monitors. A preliminary analysis suggests the presence of >10 GeV particles.

The Milagrito water Cherenkov telescope operated for over a year. The most probable gamma-ray energy was ~1 TeV and the trigger rate was as high as 400 Hz. We have developed an efficient technique for searching the entire sky for short duration bursts of TeV photons. Such bursts may result from "traditional" gamma-ray bursts that were not in the field-of-view of any other instruments, the evaporation of primordial black holes, or some as yet undiscovered phenomenon. We have begun to search the Milagrito data set for bursts of duration 10 seconds. Here we will present the technique and the expected results. Final results will be presented at the conference.

Observing gamma ray bursts (GRBs) in the TeV energy range can be extremely valuable in providing insight to GRB radiation mechanisms and in constraining source distances. The Milagrito detector was an air shower array which used the water Cherenkov technique to search for TeV sources. Data from this detector was analyzed to look for a TeV component of GRBs coincident with low energy -rays detected by the BATSE instrument on the Compton Gamma Ray Observatory. A sample of 54 BATSE GRBs which were in the field of view of the Milagrito detector during its lifetime (February 1997 to May 1998) was used.

The Milagro water Cherenkov detector near Los Alamos, New Mexico, has been operated as a sky monitor at energies of a few TeV between February 1997 and April 1998. Serving as a test run for the full Milagro detector, Milagrito has taken data during the strong and long-lasting 1997 flare of Mrk 501. We present results from the analysis of Mrk 501 and compare the excess and background rates with expectations from the detector simulations.

The Milagro Gamma Ray Observatory is the world's first large-area water Cherenkov detector capable of continuously monitoring the sky at TeV energies. Located in northern New Mexico, Milagro will perform an all sky survey of the Northern Hemisphere at energies between ~250 GeV and 50 TeV. With a high duty cycle, large detector area (~5000 square meters), and a wide field-of-view (~1 sr), Milagro is uniquely capable of searching for transient and DC sources of high-energy gamma-ray emission. Milagro has been operating since February, 1999. The current status of the Milagro Observatory and initial results will be discussed.

The Milagrito water Cherenkov detector in the Jemez Mountains near Los Alamos, New Mexico took data from February 1997 to April 1998. Milagrito served as a prototype for the larger Milagro detector, which has just begun operations. Milagrito was the first large-aperture gamma-ray detector with sensitivity to gamma rays below 1 TeV. We report here on a search for steady emission from point sources over most of the northern sky using data from Milagrito.