The All-sky Medium Energy Gamma-ray Observatory (AMEGO) is a probe class mission concept that will provide essential contributions to multimessenger astrophysics in the late 2020s and beyond. AMEGO combines high sensitivity in the 200 keV to 10 GeV energy range with a wide field of view, good spectral resolution, and polarization sensitivity. Therefore, AMEGO is key in the study of multimessenger astrophysical objects that have unique signatures in the gamma-ray regime, such as neutron star mergers, supernovae, and flaring active galactic nuclei. The order-of-magnitude improvement compared to previous MeV missions also enables discoveries of a wide range of phenomena whose energy output peaks in the relatively unexplored medium-energy gamma-ray band.

We present a likelihood profile stacking technique based on the {\it Fermi}-Large Area Telescope (LAT) data to explore the $\gamma$-ray characteristics of {\it Fermi}-LAT undetected astrophysical populations. The pipeline is applied to a sample of $\gamma$-ray unresolved extreme blazars, i.e., sources with the highest synchrotron peak frequencies ($\nu_{\rm Syn}^{\rm peak}\geqslant 10^{17}$ Hz), and we report a cumulative $\gamma$-ray detection with more than 32$\sigma$ confidence for 2 degrees of freedom. Comparing the generated stacked $\gamma$-ray spectrum with the sensitivity limits of the upcoming Cherenkov Telescope Array (CTA), we find that the {\it Fermi}-LAT undetected population of such extreme blazars, on average, may remain well below the CTA detection threshold due to their faintness and extragalactic background light (EBL) absorption. However, $\gamma$-ray detected blazars belonging to the same class are promising candidates for CTA observations. The EBL corrected stacked spectra of these sources do not show any softening up to 1 TeV. This finding suggests the inverse Compton peak of extreme blazars to lie above 1 TeV, thus indicating a hard intrinsic TeV spectrum. Our analysis also predicts that at 100 GeV, at least $\sim$10\% of the diffuse extragalactic $\gamma$-ray background originates from the $\gamma$-ray undetected extreme blazars. These results highlight the effectiveness of the developed stacking technique to explore the uncharted territory of $\gamma$-ray undetected astrophysical objects.

One of the predictions of the $\Lambda$CDM cosmological framework is the hierarchical formation of structure, giving rise to dark matter (DM) halos and subhalos. When the latter are massive enough they retain gas (i.e., baryons) and become visible. This is the case of the dwarf satellite galaxies in the Milky Way (MW). Below a certain mass, halos may not accumulate significant amounts of baryons and remain completely dark. However, if DM particles are Weakly Interacting Massive Particles (WIMPs), we expect them to annihilate in subhalos, producing gamma rays which can be detected with the Fermi satellite. Using the three most recent point-source Fermi Large Area Telescope (LAT) catalogs (3FGL, 2FHL and 3FHL), we search for DM subhalo candidates among the unidentified sources, i.e., sources with no firm association to a known astrophysical object. We apply several selection criteria based on the expected properties of the DM-induced emission from subhalos, which allow us to significantly reduce the list of potential candidates. Then, by characterizing the minimum detection flux of the instrument and comparing our sample to predictions from the Via Lactea II (VL-II) N-body cosmological simulation, we place conservative and robust constraints on the $\langle\sigma v\rangle-m_{DM}$ parameter space. For annihilation via the $\tau^+\tau^-$ channel, we put an upper limit of $4\times 10^{-26}~(5\times 10^{-25})~cm^3~s^{-1}$ for a mass of 10 (100) GeV. A critical improvement over previous treatments is the repopulation we made to include low-mass subhalos below the VL-II mass resolution. With more advanced subhalo candidate filtering the sensitivity reach of our method can potentially improve these constraints by a factor 3 (2) for $\tau^+\tau^-$ ($b \bar{b}$) channel.

MeV blazars are the most luminous persistent sources in the Universe and emit most of their energy in the MeV band. These objects display very large jet powers and accretion luminosities and are known to host black holes with a mass often exceeding $10^9 M_{\odot}$. An MeV survey, performed by a new generation MeV telescope which will bridge the entire energy and sensitivity gap between the current generation of hard X-ray and gamma-ray instruments, will detect $>$1000 MeV blazars up to a redshift of $z=5-6$. Here we show that this would allow us: 1) to probe the formation and growth mechanisms of supermassive black holes at high redshifts, 2) to pinpoint the location of the emission region in powerful blazars, 3) to determine how accretion and black hole spin interplay to power the jet.

We test different physically motivated models for the spectral shape of the $\gamma$-ray emission in a sample of 128 blazars with known redshifts detected by the Fermi Large Area Telescope (LAT) at energies above 50 GeV. The first nine years of LAT data in the energy range from 300 MeV to 2 TeV are analyzed in order to extend the spectral energy coverage of the 2FHL blazars in our sample. We compare these spectral data to four leptonic models for the production of $\gamma$-rays through Compton scattering by a population of electrons with different spectral shapes. In the first three models we consider Compton scattering in the Thomson regime with different acceleration mechanisms for the electrons. In the fourth model we consider Compton scattering by a pure power law distribution of electrons with spectral curvature due to scattering in the Klein-Nishina regime. The majority of blazar $\gamma$-ray spectra are preferentially fit with either a power law with exponential cut-off in the Thomson regime or a power law electron distribution with Compton scattering in the Klein-Nishina regime, while a log-parabola with a low-energy power-law and broken power-law spectral shape in the Thomson regime appears systematically disfavoured, which is likely a consequence of the restriction to pure Thomson scattering which we imposed on those models. This finding may be an indication that the $\gamma$-ray emission from FSRQs in the 2FHL catalog is dominated by Compton scattering of radiation from the dusty torus, while in the case of BL Lac objects, it is dominated by synchrotron self-Compton radiation.

The light emitted by all galaxies over the history of the Universe produces the extragalactic background light (EBL) at ultraviolet, optical, and infrared wavelengths. The EBL is a source of opacity for $\gamma$ rays via photon-photon interactions, leaving an imprint in the spectra of distant $\gamma$-ray sources. We measure this attenuation using {739} active galaxies and one gamma-ray burst detected by the {\it Fermi} Large Area Telescope. This allows us to reconstruct the evolution of the EBL and determine the star-formation history of the Universe over 90\% of cosmic time. Our star-formation history is consistent with independent measurements from galaxy surveys, peaking at redshift $z\sim2$. Upper limits of the EBL at the epoch of re-ionization suggest a turnover in the abundance of faint galaxies at $z\sim 6$.

The HAWC (High Altitude Water Cherenkov) collaboration recently published their 2HWC catalog, listing 39 very high energy (VHE; >100~GeV) gamma-ray sources based on 507 days of observation. Among these, there are nineteen sources that are not associated with previously known TeV sources. We have studied fourteen of these sources without known counterparts with VERITAS and Fermi-LAT. VERITAS detected weak gamma-ray emission in the 1~TeV-30~TeV band in the region of DA 495, a pulsar wind nebula coinciding with 2HWC J1953+294, confirming the discovery of the source by HAWC. We did not find any counterpart for the selected fourteen new HAWC sources from our analysis of Fermi-LAT data for energies higher than 10 GeV. During the search, we detected GeV gamma-ray emission coincident with a known TeV pulsar wind nebula, SNR G54.1+0.3 (VER J1930+188), and a 2HWC source, 2HWC J1930+188. The fluxes for isolated, steady sources in the 2HWC catalog are generally in good agreement with those measured by imaging atmospheric Cherenkov telescopes. However, the VERITAS fluxes for SNR G54.1+0.3, DA 495, and TeV J2032+4130 are lower than those measured by HAWC and several new HAWC sources are not detected by VERITAS. This is likely due to a change in spectral shape, source extension, or the influence of diffuse emission in the source region.

e-ASTROGAM (enhanced ASTROGAM) is a breakthrough Observatory space mission, with a detector composed by a Silicon tracker, a calorimeter, and an anticoincidence system, dedicated to the study of the non-thermal Universe in the photon energy range from 0.3 MeV to 3 GeV - the lower energy limit can be pushed to energies as low as 150 keV for the tracker, and to 30 keV for calorimetric detection. The mission is based on an advanced space-proven detector technology, with unprecedented sensitivity, angular and energy resolution, combined with polarimetric capability. Thanks to its performance in the MeV-GeV domain, substantially improving its predecessors, e-ASTROGAM will open a new window on the non-thermal Universe, making pioneering observations of the most powerful Galactic and extragalactic sources, elucidating the nature of their relativistic outflows and their effects on the surroundings. With a line sensitivity in the MeV energy range one to two orders of magnitude better than previous generation instruments, e-ASTROGAM will determine the origin of key isotopes fundamental for the understanding of supernova explosion and the chemical evolution of our Galaxy. The mission will provide unique data of significant interest to a broad astronomical community, complementary to powerful observatories such as LIGO-Virgo-GEO600-KAGRA, SKA, ALMA, E-ELT, TMT, LSST, JWST, Athena, CTA, IceCube, KM3NeT, and LISA.

BL Lacertae (BL Lac) objects are the prominent members of the third {\it Fermi} Large Area Telescope catalog of $\gamma$-ray sources. Half of the BL Lac population ($\sim$ 300) lack redshift measurements, which is due to the absence of lines in their optical spectrum, thereby making it difficult to utilize spectroscopic methods. Our photometric drop-out technique can be used to establish the redshift for a fraction of these sources. This work employed 6 filters mounted on the $Swift$-UVOT and 4 optical filters on two telescopes, the 0.65 m SARA-CTIO in Chile and 1.0 m SARA-ORM in the Canary Islands, Spain. A sample of 15 sources was extracted from the $Swift$ archival data for which 6 filter UVOT observations were conducted. By complementing the {\it Swift} observations with the SARA ones, we were able to discover two high redshift sources: 3FGL J1155.4-3417 and 3FGL J1156.7-2250 at $z=1.83^{+0.10}_{-0.13}$ and $z=1.73^{+0.11}_{-0.19}$ , respectively, resulting from the dropouts in the powerlaw template fits to these data. The discoveries add to the important (26 total) sample of high-redshift BL Lacs. While the sample of high-z BL Lacs is still rather small, these objects do not seem to fit well within known schemes of the blazar population and represent the best probes of the extragalactic background light.

The extragalactic background light (EBL), from ultraviolet to infrared wavelengths, is predominantly due to emission from stars, accreting black holes and reprocessed light due to Galactic dust. The EBL can be studied through the imprint it leaves, via $\gamma$-$\gamma$ absorption of high-energy photons, in the spectra of distant $\gamma$-ray sources. The EBL has been probed through the search for the attenuation it produces in the spectra of BL Lacertae (BL Lac) objects and individual $\gamma$-ray bursts (GRBs). GRBs have significant advantages over blazars for the study of the EBL especially at high redshifts. Here we analyze a combined sample of twenty-two GRBs, detected by the Fermi Large Area Telescope between 65 MeV and 500 GeV. We report a marginal detection (at the ~2.8 $\sigma$ level) of the EBL attenuation in the stacked spectra of the source sample. This measurement represents a first constraint of the EBL at an effective redshift of ~1.8. We combine our results with prior EBL constraints and conclude that Fermi-LAT is instrumental to constrain the UV component of the EBL. We discuss the implications on existing empirical models of EBL evolution.

We analyze the spectra of 300,000 luminous red galaxies (LRGs) with stellar masses $M_* \gtrsim 10^{11} M_{\odot}$ from the SDSS-III Baryon Oscillation Spectroscopic Survey (BOSS). By studying their star-formation histories, we find two main evolutionary paths converging into the same quiescent galaxy population at $z\sim0.55$. Fast-growing LRGs assemble $80\%$ of their stellar mass very early on ($z\sim5$), whereas slow-growing LRGs reach the same evolutionary state at $z\sim1.5$. Further investigation reveals that their clustering properties on scales of $\sim$1-30 Mpc are, at a high level of significance, also different. Fast-growing LRGs are found to be more strongly clustered and reside in overall denser large-scale structure environments than slow-growing systems, for a given stellar-mass threshold. Our results imply a dependence of clustering on stellar-mass assembly history (naturally connected to the mass-formation history of the corresponding halos) for a homogeneous population of similar mass and color, which constitutes a strong observational evidence of galaxy assembly bias.

We describe and report the status of a neutrino-triggered program in IceCube that generates real-time alerts for gamma-ray follow-up observations by atmospheric-Cherenkov telescopes (MAGIC and VERITAS). While IceCube is capable of monitoring the whole sky continuously, high-energy gamma-ray telescopes have restricted fields of view and in general are unlikely to be observing a potential neutrino-flaring source at the time such neutrinos are recorded. The use of neutrino-triggered alerts thus aims at increasing the availability of simultaneous multi-messenger data during potential neutrino flaring activity, which can increase the discovery potential and constrain the phenomenological interpretation of the high-energy emission of selected source classes (e.g. blazars). The requirements of a fast and stable online analysis of potential neutrino signals and its operation are presented, along with first results of the program operating between 14 March 2012 and 31 December 2015.

Determining redshifts for BL Lacertae (BL Lac) objects using the traditional spectroscopic method is challenging due to the absence of strong emission lines in their optical spectra. We employ the photometric dropout technique to determine redshifts for this class of blazars using the combined 13 broad-band filters from Swift-UVOT and the multi-channel imager GROND at the MPG 2.2 m telescope at ESO's La Silla Observatory. The wavelength range covered by these 13 filters extends from far ultraviolet to the near-Infrared. We report results on 40 new Fermi detected BL Lacs with the photometric redshifts determinations for 5 sources, with 3FGL J1918.2-4110 being the most distant in our sample at z=2.16. Reliable upper limits are provided for 20 sources in this sample. Using the highest energy photons for these Fermi-LAT sources, we evaluate the consistency with the Gamma-ray horizon due to the extragalactic background light.

LOFT-P is a concept for a NASA Astrophysics Probe-Class (<$1B) X-ray timing mission, based on the LOFT concept originally proposed to ESAs M3 and M4 calls. LOFT-P requires very large collecting area (>6 m^2, >10x RXTE), high time resolution, good spectral resolution, broad-band spectral coverage (2-30 keV), highly flexible scheduling, and an ability to detect and respond promptly to time-critical targets of opportunity. It addresses science questions such as: What is the equation of state of ultra dense matter? What are the effects of strong gravity on matter spiraling into black holes? It would be optimized for sub-millisecond timing to study phenomena at the natural timescales of neutron star surfaces and black hole event horizons and to measure mass and spin of black holes. These measurements are synergistic to imaging and high-resolution spectroscopy instruments, addressing much smaller distance scales than are possible without very long baseline X-ray interferometry, and using complementary techniques to address the geometry and dynamics of emission regions. A sky monitor (2-50 keV) acts as a trigger for pointed observations, providing high duty cycle, high time resolution monitoring of the X-ray sky with ~20 times the sensitivity of the RXTE All-Sky Monitor, enabling multi-wavelength and multi-messenger studies. A probe-class mission concept would employ lightweight collimator technology and large-area solid-state detectors, technologies which have been recently greatly advanced during the ESA M3 study. Given the large community interested in LOFT (>800 supporters, the scientific productivity of this mission is expected to be very high, similar to or greater than RXTE (~2000 refereed publications). We describe the results of a study, recently completed by the MSFC Advanced Concepts Office, that demonstrates that LOFT-P is feasible within a NASA probe-class mission budget.

The Geminga pulsar, one of the brighest gamma-ray sources, is a promising candidate for emission of very-high-energy (VHE > 100 GeV) pulsed gamma rays. Also, detection of a large nebula have been claimed by water Cherenkov instruments. We performed deep observations of Geminga with the MAGIC telescopes, yielding 63 hours of good-quality data, and searched for emission from the pulsar and pulsar wind nebula. We did not find any significant detection, and derived 95% confidence level upper limits. The resulting upper limits of 5.3 x 10^{-13} TeV cm^{-2} s^{-1} for the Geminga pulsar and 3.5 x 10^{-12} TeV cm^{-2} s^{-1} for the surrounding nebula at 50 GeV are the most constraining ones obtained so far at VHE. To complement the VHE observations, we also analyzed 5 years of Fermi-LAT data from Geminga, finding that the sub-exponential cut-off is preferred over the exponential cut-off that has been typically used in the literature. We also find that, above 10 GeV, the gamma-ray spectra from Geminga can be described with a power law with index softer than 5. The extrapolation of the power-law Fermi-LAT pulsed spectra to VHE goes well below the MAGIC upper limits, indicating that the detection of pulsed emission from Geminga with the current generation of Cherenkov telescopes is very difficult.

We report for the first time a gamma-ray and multi-wavelength nearly-periodic oscillation in an active galactic nucleus. Using the Fermi Large Area Telescope (LAT) we have discovered an apparent quasi-periodicity in the gamma-ray flux (E >100 MeV) from the GeV/TeV BL Lac object PG 1553+113. The marginal significance of the 2.18 +/-0.08 year-period gamma-ray cycle is strengthened by correlated oscillations observed in radio and optical fluxes, through data collected in the OVRO, Tuorla, KAIT, and CSS monitoring programs and Swift UVOT. The optical cycle appearing in ~10 years of data has a similar period, while the 15 GHz oscillation is less regular than seen in the other bands. Further long-term multi-wavelength monitoring of this blazar may discriminate among the possible explanations for this quasi-periodicity.

This paper is the second in a series describing the southern Galactic Disk Survey (GDS) performed at the Universit\"atssternwarte Bochum near Cerro Armazones in Chile. Haas et al. (2012, Paper I) presented the survey design and a the characteristics of the observations and data. They identified ~2200 variable stars in an area of 50 square degrees with more than 50 observations in 2011. Here we present the first complete version of the GDS covering all 268 fields with 1323 square degrees along the Galactic disk including revised data from Paper I. The individual fields were observed up to 272 times and comprise a maximum time span between September 2010 and May 2015. We detect a total of 64151 variable sources, which are presented in a catalog including some of their properties and their light curves. A comparison with the International Variable Star Index (VSX) and All Sky Automated Survey (ASAS) indicates that 56794 of these sources are previously unknown variables. Furthermore, we present UBVr'i'z' photometry for all sources within the GDS, resulting in a new multi-color catalog of nearly 16x10^6 sources detected in at least one filter. Both the GDS and the near-infrared VISTA Variables in the Via Lactea survey (VVV) complement each other in the overlap area of about 300 square degrees enabling future comparison studies.

We estimate the galaxy density along lines of sight to hard extragalactic gamma-ray sources by correlating source positions on the sky with a void catalog based on the Sloan Digital Sky Survey (SDSS). Extragalactic gamma-ray sources that are detected at very high energy (VHE; E>100 GeV) or have been highlighted as VHE-emitting candidates in the Fermi Large Area Telescope hard source catalog (together referred to as "VHE-like" sources) are distributed along underdense lines of sight at the 2.4 sigma level. There is also a less suggestive correlation for the Fermi hard source population (1.7 sigma). A correlation between 10-500 GeV flux and underdense fraction along the line of sight for VHE-like and Fermi hard sources is found at 2.4 sigma and 2.6 sigma, respectively. The preference for underdense sight lines is not displayed by gamma-ray emitting galaxies within the second Fermi catalog, containing sources detected above 100 MeV, or the SDSS DR7 quasar catalog. We investigate whether this marginal correlation might be a result of lower extragalactic background light (EBL) photon density within the underdense regions and find that, even in the most extreme case of a entirely underdense sight line, the EBL photon density is only 2% less than the nominal EBL density. Translating this into gamma-ray attenuation along the line of sight for a highly attenuated source with opacity tau(E,z) ~5, we estimate that the attentuation of gamma-rays decreases no more than 10%. This decrease, although non-neglible, is unable to account for the apparent hard source correlation with underdense lines of sight.

The discovery of rapidly variable Very High Energy (VHE; E > 100 GeV) gamma-ray emission from 4C +21.35 (PKS 1222+216) by MAGIC on 2010 June 17, triggered by the high activity detected by the Fermi Large Area Telescope (LAT) in high energy (HE; E > 100 MeV) gamma-rays, poses intriguing questions on the location of the gamma-ray emitting region in this flat spectrum radio quasar. We present multifrequency data of 4C +21.35 collected from centimeter to VHE during 2010 to investigate the properties of this source and discuss a possible emission model. The first hint of detection at VHE was observed by MAGIC on 2010 May 3, soon after a gamma-ray flare detected by Fermi-LAT that peaked on April 29. The same emission mechanism may therefore be responsible for both the HE and VHE emission during the 2010 flaring episodes. Two optical peaks were detected on 2010 April 20 and June 30, close in time but not simultaneous with the two gamma-ray peaks, while no clear connection was observed between the X-ray an gamma-ray emission. An increasing flux density was observed in radio and mm bands from the beginning of 2009, in accordance with the increasing gamma-ray activity observed by Fermi-LAT, and peaking on 2011 January 27 in the mm regime (230 GHz). We model the spectral energy distributions (SEDs) of 4C +21.35 for the two periods of the VHE detection and a quiescent state, using a one-zone model with the emission coming from a very compact region outside the broad line region. The three SEDs can be fit with a combination of synchrotron self-Compton and external Compton emission of seed photons from a dust torus, changing only the electron distribution parameters between the epochs. The fit of the optical/UV part of the spectrum for 2010 April 29 seems to favor an inner disk radius of <6 gravitational radii, as one would expect from a prograde-rotating Kerr black hole.

We present results from near-infrared spectroscopy of 26 emission-line galaxies at z ~ 2 obtained with the FIRE spectrometer on the Magellan Baade telescope. The sample was selected from the WISP survey, which uses the near-infrared grism of the Hubble Space Telescope Wide Field Camera 3 to detect emission-line galaxies over 0.3 < z < 2.3. Our FIRE follow-up spectroscopy (R~5000) over 1.0-2.5 micron permits detailed measurements of physical properties of the z~2 emission-line galaxies. Dust-corrected star formation rates for the sample range from ~5-100 M_sun yr-1. We derive a median metallicity for the sample of ~0.45 Z_sun, and the estimated stellar masses range from ~10^8.5 - 10^9.5 M_sun. The average ionization parameters measured for the sample are typically much higher than what is found for local star-forming galaxies. We derive composite spectra from the FIRE sample, from which we infer typical nebular electron densities of ~100-400 cm^-3. Based on the location of the galaxies and composite spectra on BPT diagrams, we do not find evidence for significant AGN activity in the sample. Most of the galaxies as well as the composites are offset in the BPT diagram toward higher [O III]/H-beta at a given [N II]/H-alpha, in agreement with other observations of z > 1 star-forming galaxies, but composite spectra derived from the sample do not show an appreciable offset from the local star-forming sequence on the [O III]/H-beta versus [S II]/H-alpha diagram. We infer a high nitrogen-to-oxygen abundance ratio from the composite spectrum, which may contribute to the offset of the high-redshift galaxies from the local star-forming sequence in the [O III]/H-beta versus [N II]/H-alpha diagram. We speculate that the elevated nitrogen abundance could result from substantial numbers of Wolf-Rayet stars in starbursting galaxies at z~2. (Abridged)

We identified the z~2 Lyman break galaxies using deep HST ultraviolet (F275W/F336W) imaging of Abell 1689. Because of the imaging depth and the large magnification provided by the cluster, we detect galaxies 100x fainter (-19.5< M_1500 <-13) than previous surveys at this redshift. We are able to calculate the intrinsic sensitivity of the observations as a function of source plane position, allowing determinations of effective volume as a function of luminosity. We fit the faint-end slope of the luminosity function to be alpha = -1.74 +/-0.08, consistent with the values obtained for 2.5 < z < 6. There is no turnover in the luminosity function down to MUV = -13. The trend of increasingly redder UV spectral slopes with luminosity at higher redshifts is observed in our sample, but with redder slopes at all luminosities and average reddening of < E(B - V) >= 0.15. We assume the stars in these galaxies are metal poor (0.2Z_sun) compared to their brighter counterparts (Z_sun), resulting in bluer assumed intrinsic UV slopes and larger derived dust extinction. The total UV luminosity density at z ~ 2 is 4.31x10^26 erg/s/Hz/Mpc^3, more than 70% of which is emitted by galaxies in the luminosity range of our sample. We determine the star formation rate density at z ~ 2 (assuming constant dust extinction correction of 4.2 over all luminosities and a Kroupa IMF) of 0.148 M/yr/Mpc^3, significantly higher than previous determinations because of the additional population of fainter galaxies and the larger dust correction factors.[abridged]

We present near-infrared emission line counts and luminosity functions from the HST WFC3 Infrared Spectroscopic Parallels (WISP) program for 29 fields (0.037 deg^2) observed using both the G102 and G141 grisms. Altogether we identify 1048 emission line galaxies with observed equivalent widths greater than 40 Angstroms, 467 of which have multiple detected emission lines. The WISP survey is sensitive to fainter flux levels (3-5x10^-17 ergs/s/cm^2) than the future space near-infrared grism missions aimed at baryonic acoustic oscillation cosmology (1-4x10^-16 ergs/s/cm^2), allowing us to probe the fainter emission line galaxies that the shallower future surveys may miss. Cumulative number counts of 0.7<z<1.5 galaxies reach 10,000 deg^-2 above an H-alpha flux of 2x10^-16 ergs/s/cm^2. H-alpha-emitting galaxies with comparable [OIII] flux are roughly 5 times less common than galaxies with just H-alpha emission at those flux levels. Galaxies with low H-alpha/[OIII] ratios are very rare at the brighter fluxes that future near-infrared grism surveys will probe; our survey finds no galaxies with H-alpha/[OIII] < 0.95 that have H-alpha flux greater than 3x10^-16 ergs/s/cm^2. Our H-alpha luminosity function contains a comparable number density of faint line emitters to that found by the NICMOS near-infrared grism surveys, but significantly fewer (factors of 3-4 less) high luminosity emitters. We also find that our high redshift (z=0.9-1.5) counts are in agreement with the high redshift (z=1.47) narrow band H-alpha survey of HiZELS (Sobral et al. 2013), while our lower redshift luminosity function (z=0.3-0.9) falls slightly below their z=0.84 result. The evolution in both the H-alpha luminosity function from z=0.3--1.5 and the [OIII] luminosity function from z=0.7-2.3 is almost entirely in the L* parameter, which steadily increases with redshift over those ranges.

We present the first robust measurement of the high redshift mass-metallicity (MZ) relation at 10^{8}< M/M_{\sun} < 10^{10}, obtained by stacking spectra of 83 emission-line galaxies with secure redshifts between 1.3 < z < 2.3. For these redshifts, infrared grism spectroscopy with the Hubble Space Telescope Wide Field Camera 3 is sensitive to the R23 metallicity diagnostic: ([OII]3726,3729 + [OIII] 4959,5007)/H\beta. Using spectra stacked in four mass quartiles, we find a MZ relation that declines significantly with decreasing mass, extending from 12+log(O/H) = 8.8 at M=10^{9.8} M_{\sun} to 12+log(O/H)= 8.2 at M=10^{8.2} M_{\sun}. After correcting for systematic offsets between metallicity indicators, we compare our MZ relation to measurements from the stacked spectra of galaxies with M>10^{9.5} M_{\sun} and z~2.3. Within the statistical uncertainties, our MZ relation agrees with the z~2.3 result, particularly since our somewhat higher metallicities (by around 0.1 dex) are qualitatively consistent with the lower mean redshift z=1.76 of our sample. For the masses probed by our data, the MZ relation shows a steep slope which is suggestive of feedback from energy-driven winds, and a cosmological downsizing evolution where high mass galaxies reach the local MZ relation at earlier times. In addition, we show that our sample falls on an extrapolation of the star-forming main sequence (the SFR-M_{*} relation) at this redshift. This result indicates that grism emission-line selected samples do not have preferentially high SFRs. Finally, we report no evidence for evolution of the mass-metallicity-SFR plane; our stack-averaged measurements show excellent agreement with the local relation.

We combine Hubble Space Telescope (HST) G102 & G141 NIR grism spectroscopy with HST/WFC3-UVIS, HST/WFC3-IR and Spitzer/IRAC[3.6\mu m] photometry to assemble a sample of massive (log(M_star/M_sun) ~ 11) and quenched galaxies at z~1.5. Our sample of 41 galaxies is the largest with G102+G141 NIR spectroscopy for quenched sources at these redshifts. In contrast to the local Universe, z~1.5 quenched galaxies in the high-mass range have a wide range of stellar population properties. We find their SEDs are well fitted with exponentially decreasing SFHs, and short star-formation time-scales (\tau<100Myr). Quenched galaxies also show a wide distribution in ages, between 1-4Gyr. In the (u-r)_0-versus-mass space quenched galaxies have a large spread in rest-frame color at a given mass. Most quenched galaxies populate the z~1.5 red-sequence (RS), but an important fraction of them (32%) have substantially bluer colors. Although with a large spread, we find that the quenched galaxies ON the RS have older median ages (3.1Gyr) than the quenched galaxies OFF the RS (1.5Gyr). We also show that a rejuvenated SED cannot reproduce the observed stacked spectra of (the bluer) quenched galaxies OFF the RS. We derive the upper limit on the fraction of massive galaxies ON the RS at z~1.5 to be <43%. We speculate that the young quenched galaxies OFF the RS are in a transition phase between vigorous star formation at z>2 and the z~1.5 RS. According to their estimated ages, the time required for quenched galaxies OFF the RS to join their counterparts ON the z~1.5 RS is of the order of ~1Gyr.

The shutdown of star formation in galaxies is generally termed `quenching'. Although quenching may occur through a variety of processes, the exact mechanism(s) that is in fact responsible for quenching is still in question. This paper addresses quenching by searching for traces of possible quenching processes through their effects on galaxy structural parameters such as surface stellar mass density and Sersic index (n). We analyze the rest-frame U-B color correlations versus these structural parameters using a sample of galaxies in the redshift range 0.5< z<0.8 from the DEEP2/AEGIS survey. We find that Sersic index (n) has the smallest overlap region among all tested parameters and resembles a step-function with a threshold value of n=2.3. There exists, however, a significant population of outliers with blue colors yet high n values that seem to contradict this behavior. We hypothesize that their Sersic values may be distorted by bursts of star formation, AGNs, and/or poor fits, leading us to consider central surface stellar mass density as an alternative to Sersic index. Not only does it correct the outliers, it also forms a tight relationship with color, suggesting that the innermost structure of galaxies is most physically linked with quenching. Furthermore, at z~0.65, the majority of the blue cloud galaxies cannot simply fade onto the red sequence since their GIM2D bulge masses are only half as large on average as the bulge masses of similar red sequence galaxies, thus demonstrating that stellar mass must absolutely increase at the centers of galaxies as they quench. We discuss a two-stage model for quenching in which galaxy star formation rates are controlled by their dark halos while they are still in the blue cloud and a second quenching process sets in later, associated with the central stellar mass build-up.

Context: The blazar Markarian 421 is one of the brightest TeV gamma-ray sources of the northern sky. From December 2007 until June 2008 it was intensively observed in the very high energy (VHE, E > 100 GeV) band by the single-dish Major Atmospheric Gamma-ray Imaging Cherenkov telescope (MAGIC-I). Aims: We aimed to measure the physical parameters of the emitting region of the blazar jet during active states. Methods: We performed a dense monitoring of the source in VHE with MAGIC-I, and also collected complementary data in soft X-rays and optical-UV bands; then, we modeled the spectral energy distributions (SED) derived from simultaneous multi-wavelength data within the synchrotron self--compton (SSC) framework. Results: The source showed intense and prolonged gamma-ray activity during the whole period, with integral fluxes (E > 200 GeV) seldom below the level of the Crab Nebula, and up to 3.6 times this value. Eight datasets of simultaneous optical-UV (KVA, Swift/UVOT), soft X-ray (Swift/XRT) and MAGIC-I VHE data were obtained during different outburst phases. The data constrain the physical parameters of the jet, once the spectral energy distributions obtained are interpreted within the framework of a single-zone SSC leptonic model. Conclusions: The main outcome of the study is that within the homogeneous model high Doppler factors (40 <= delta <= 80) are needed to reproduce the observed SED; but this model cannot explain the observed short time-scale variability, while it can be argued that inhomogeneous models could allow for less extreme Doppler factors, more intense magnetic fields and shorter electron cooling times compatible with hour or sub-hour scale variability.

Star-forming dwarfs are studied to elucidate the physical underpinnings of their fundamental plane. It is confirmed that residuals in the Tully-Fisher relation are correlated with surface brightness, but that even after accommodating the surface brightness dependence through the dwarf fundamental plane, residuals in absolute magnitude are far larger than expected from observational errors. Rather, a more fundamental plane is identified which connects the potential to HI line width and surface brightness. Residuals correlate with the axis ratio in a way which can be accommodated by recognizing the galaxies to be oblate spheroids viewed at varying angles. Correction of surface brightnesses to face-on leads to a correlation among the potential, line width, and surface brightness for which residuals are entirely attributable to observational uncertainties. The mean mass-to-light ratio of the diffuse component of the galaxies is constrained to be 0.88 +/- 0.20 in Ks. Blue compact dwarfs lie in the same plane as dwarf irregulars. The dependence of the potential on line width is less strong than expected for virialized systems, but this may be because surface brightness is acting as a proxy for variations in the mass-to-light ratio from galaxy to galaxy. Altogether, the observations suggest that gas motions are predominantly disordered and isotropic, that they are a consequence of gravity, not turbulence, and that the mass and scale of dark matter haloes scale with the amount and distribution of luminous matter. The tight relationship between the potential and observables offers the promise of determining distances to unresolved star-forming dwarfs to an accuracy comparable to that provided by the Tully-Fisher relation for spirals.

We combine a semi-analytic model of galaxy formation with simple analytic recipes describing the absorption and re-emission of starlight by dust in the interstellar medium of galaxies. We use the resulting models to predict galaxy counts and luminosity functions from the far-ultraviolet to the sub-mm, from redshift five to the present, and compare with an extensive compilation of observations. We find that in order to reproduce the rest-UV and optical luminosity functions at high redshift, we must assume an evolving normalization in the dust-to-metal ratio, implying that galaxies of a given bolometric luminosity (or metal column density) must be less extinguished than their local counterparts. In our best-fit model, we find remarkably good agreement with observations from rest-frame 1500 Angstroms to 250 microns. At longer wavelengths, most dramatically in the sub-mm, our models underpredict the number of bright galaxies by a large factor. The models reproduce the observed total IR luminosity function fairly well. We show the results of varying several ingredients of the models, including various aspects of the dust attenuation recipe, the dust emission templates, and the cosmology. We use our models to predict the integrated Extragalactic Background Light (EBL), and compare with an observationally-motivated EBL model and with other available observational constraints.

Discoveries of very high energy (VHE) photons from distant blazars suggest that, after correction by extragalactic background light (EBL) absorption, there is a flatness or even a turn-up in their spectra at the highest energies that cannot be easily explained by the standard framework. Here, it is shown that a possible solution to this problem is achieved by assuming the existence of axion-like particles (ALPs) with masses ~1 neV. The ALP scenario is tested making use of observations of the highest redshift blazars known in the VHE energy regime, namely 3C 279, 3C 66A, PKS 1222+216 and PG 1553+113. In all cases, better fits to the observed spectra are found when including ALPs rather than considering EBL only. Interestingly, quite similar critical energies for photon/ALP conversions are also derived, independently of the source considered.

Data from (non-) attenuation of gamma rays from active galactic nuclei (AGN) and gamma ray bursts (GRBs) give upper limits on the extragalactic background light (EBL) from the UV to the mid-IR that are only a little above the lower limits from observed galaxies. These upper limits now rule out some EBL models and purported observations, with improved data likely to provide even stronger constraints. We present EBL calculations both based on multiwavelength observations of thousands of galaxies and also based on semi-analytic models, and show that they are consistent with these lower limits from observed galaxies and with the gamma-ray upper limit constraints. Such comparisons "close the loop" on cosmological galaxy formation models, since they account for all the light, including that from galaxies too faint to see. We compare our results with those of other recent works, and discuss the implications of these new EBL calculations for gamma ray attenuation. Catching a few GRBs with groundbased atmospheric Cherenkov Telescope (ACT) arrays or water Cherenkov detectors could provide important new constraints on the high-redshift star formation history of the universe.

The extragalactic background light (EBL) is of fundamental importance both for understanding the entire process of galaxy evolution and for gamma-ray astronomy. However, the overall spectrum of the EBL between 0.1 and 1000 microns has never been determined directly neither from observed luminosity functions (LFs), over a wide redshift range, nor from any multiwavelength observation of galaxy spectral energy distributions (SEDs). The evolving, overall spectrum of the EBL is derived here utilizing a novel method based on observations only. The changing fractions of quiescent galaxies, star-forming galaxies, starburst galaxies and active galactic nucleus (AGN) galaxies from redshift 0.2 to 1 are estimated, and two alternative extrapolations of SED types to higher redshifts are considered. This allows calculation of the evolving EBL. The EBL uncertainties in our modelling based directly on the data are quantified, and their consequences for attenuation of very-high-energy gamma-rays due to pair production on the EBL are discussed. It is concluded that the EBL seems well constrained from the UV to the mid-IR at an intensity level roughly matching galaxy count data. Independent efforts from IR and gamma-ray astronomy are needed in order to reduce the uncertainties in the far-IR.

Context. 3C 279, the first quasar discovered to emit VHE gamma-rays by the MAGIC telescope in 2006, was reobserved by MAGIC in January 2007 during a major optical flare and from December 2008 to April 2009 following an alert from the Fermi space telescope on an exceptionally high gamma -ray state. Aims. The January 2007 observations resulted in a detection on January 16 with significance 5.2 sigma, corresponding to a F(> 150 GeV) (3.8 \pm 0.8) \cdot 10^-11 ph cm^-2 s^-1 while the overall data sample does not show significant signal. The December 2008 - April 2009 observations did not detect the source. We study the multiwavelength behavior of the source at the epochs of MAGIC observations, collecting quasi-simultaneous data at optical and X-ray frequencies and for 2009 also gamma-ray data from Fermi. Methods. We study the light curves and spectral energy distribution of the source. The spectral energy distributions of three observing epochs (including the February 2006, which has been previously published in Albert et al. 2008a) are modeled with one-zone inverse Compton models and the emission on January 16, 2007 also with two zone model and with a lepto-hadronic model. Results. We find that the VHE gamma-ray emission detected in 2006 and 2007 challenges standard one-zone model, based on relativistic electrons in a jet scattering broad line region photons, while the other studied models fit the observed spectral energy distribution more satisfactorily.

The extragalactic background light (EBL) is of fundamental importance both for understanding the entire process of galaxy evolution and for gamma-ray astronomy, but the overall spectrum of the EBL between 0.1-1000 microns has never been determined directly from galaxy spectral energy distribution (SED) observations over a wide redshift range. The evolving, overall spectrum of the EBL is derived here utilizing a novel method based on observations only. This is achieved from the observed evolution of the rest-frame K-band galaxy luminosity function up to redshift 4 (Cirasuolo et al. 2010), combined with a determination of galaxy SED-type fractions. These are based on fitting SWIRE templates to a multiwavelength sample of about 6000 galaxies in the redshift range from 0.2 to 1 from the All-wavelength Extended Groth Strip International Survey (AEGIS). The changing fractions of quiescent galaxies, star-forming galaxies, starburst galaxies and AGN galaxies in that redshift range are estimated, and two alternative extrapolations of SED-types to higher redshifts are considered. This allows calculation of the evolution of the luminosity densities from the UV to the IR, the evolving star formation rate density of the universe, the evolving contribution to the bolometric EBL from the different galaxy populations including AGN galaxies and the buildup of the EBL. Our EBL calculations are compared with those from a semi-analytic model, from another observationally-based model and observational data. The EBL uncertainties in our modeling based directly on the data are quantified, and their consequences for attenuation of very high energy gamma-rays due to pair production on the EBL are discussed. It is concluded that the EBL is well constrained from the UV to the mid-IR, but independent efforts from infrared and gamma-ray astronomy are needed in order to reduce the uncertainties in the far-IR.

Axion Like Particles (ALPs), postulated to solve the strong-CP problem, are predicted to couple with photons in the presence of magnetic fields, which may lead to a significant change in the observed spectra of gamma-ray sources such as AGNs. Here we simultaneously consider in the same framework both the photon/axion mixing that takes place in the gamma-ray source and that one expected to occur in the intergalactic magnetic fields. We show that photon/axion mixing could explain recent puzzles regarding the observed spectra of distant gamma-ray sources as well as the recently published lower limit to the EBL intensity. We finally summarize the different signatures expected and discuss the best strategy to search for ALPs with the Fermi satellite and current Cherenkov telescopes like CANGAROO, HESS, MAGIC and VERITAS.

Axion Like Particles (ALPs) are predicted to couple with photons in the presence of magnetic fields. This effect may lead to a significant change in the observed spectra of gamma-ray sources such as AGNs. Here we carry out a detailed study that for the first time simultaneously considers in the same framework both the photon/axion mixing that takes place in the gamma-ray source and that one expected to occur in the intergalactic magnetic fields. An efficient photon/axion mixing in the source always means an attenuation in the photon flux, whereas the mixing in the intergalactic medium may result in a decrement and/or enhancement of the photon flux, depending on the distance of the source and the energy considered. Interestingly, we find that decreasing the value of the intergalactic magnetic field strength, which decreases the probability for photon/axion mixing, could result in an increase of the expected photon flux at Earth if the source is far enough. We also find a 30% attenuation in the intensity spectrum of distant sources, which occurs at an energy that only depends on the properties of the ALPs and the intensity of the intergalactic magnetic field, and thus independent of the AGN source being observed. Moreover, we show that this mechanism can easily explain recent puzzles in the spectra of distant gamma-ray sources... [ABRIDGED] The consequences that come from this work are testable with the current generation of gamma-ray instruments, namely Fermi (formerly known as GLAST) and imaging atmospheric Cherenkov telescopes like CANGAROO, HESS, MAGIC and VERITAS.

This Conference is the fifth of a series of Workshops on High Energy Gamma- ray Experiments, following the Conferences held in Perugia 2003, Bari 2004, Cividale del Friuli 2005, Elba Island 2006. This year the focus was on the use of gamma-ray to study the Dark Matter component of the Universe, the origin and propagation of Cosmic Rays, Extra Large Spatial Dimensions and Tests of Lorentz Invariance.

GAW, acronym for Gamma Air Watch, is a Research and Development experiment in the TeV range, whose main goal is to explore the feasibility of large field of view Imaging Atmospheric Cherenkov Telescopes. GAW is an array of three relatively small telescopes (2.13 m diameter) which differs from the existing and presently planned projects in two main features: the adoption of a refractive optics system as light collector and the use of single photoelectron counting as detector working mode. The optics system allows to achieve a large field of view (24x24 squared degrees) suitable for surveys of large sky regions. The single photoelectron counting mode in comparison with the charge integration mode improves the sensitivity by permitting also the reconstruction of events with a small number of collected Cherenkov photons. GAW, which is a collaboration effort of Research Institutes in Italy, Portugal and Spain, will be erected in the Calar Alto Observatory (Sierra de Los Filabres - Andalucia, Spain), at 2150 m a.s.l.). The first telescope will be settled within Autumn 2007. This paper shows the main characteristics of the experiment and its expected performance.

There is evidence of a scale-invariant matter distribution up to scales over 10 Megaparsecs. We review scaling (fractal or multifractal) models of large scale structure and their observational evidence. We conclude that the dynamics of cosmological structure formation seems to be driven to a multifractal attractor. This supports previous studies, which we review, of structure formation by means of the renormalization group within a hydrodynamic formulation.

We present a series of cosmological N-body simulations which make use of the hydrodynamic approach to the evolution of structures (Dominguez 2000). This approach addresses explicitly the existence of a finite spatial resolution and the dynamical effect of subresolution degrees of freedom. We adapt this method to cosmological simulations of the standard LCDM structure formation scenario and study the effects induced at redshift z=0 by this novel approach on the large-scale clustering patterns as well as (individual) dark matter halos. Comparing these simulations to usual N-body simulations, we find that (i) the new (hydrodynamic) model entails a proliferation of low--mass halos, and (ii) dark matter halos have a higher degree of rotational support. These results agree with the theoretical expectation about the qualitative behaviour of the "correction terms" introduced by the hydrodynamic approach: these terms act as a drain of inflow kinetic energy and a source of vorticity by the small-scale tidal torques and shear stresses.

The notion of `adhesion' has been advanced for the phenomenon of stabilization of large-scale structure emerging from gravitational instability of a cold medium. Recently, the physical origin of adhesion has been identified: a systematic derivation of the equations of motion for the density and the velocity fields leads naturally to the key equation of the `adhesion approximation' - however, under a set of strongly simplifying assumptions. In this work, we provide an evaluation of the current status of adhesive gravitational clustering and a clear explanation of the assumptions involved. Furthermore, we propose systematic generalizations with the aim to relax some of the simplifying assumptions. We start from the general Newtonian evolution equations for self-gravitating particles on an expanding Friedmann background and recover the popular `dust model' (pressureless fluid), which breaks down after the formation of density singularities; then we investigate, in a unified framework, two other models which, under the restrictions referred to above, lead to the `adhesion approximation'. We apply the Eulerian and Lagrangian perturbative expansions to these new models and, finally, we discuss some non-perturbative results that may serve as starting points for workable approximations of non-linear structure formation in the multi-stream regime. In particular, we propose a new approximation that includes, in limiting cases, the standard `adhesion model' and the Eulerian as well as Lagrangian first-order approximations.

Recently D. Vollick [Phys. Rev. D68, 063510 (2003)] has shown that the inclusion of the 1/R curvature terms in the gravitational action and the use of the Palatini formalism offer an alternative explanation for cosmological acceleration. In this work we show not only that this model of Vollick does not have a good Newtonian limit, but also that any f(R) theory with a pole of order n in R=0 and its second derivative respect to R evaluated at Ro is not zero, where Ro is the scalar curvature of background, does not have a good Newtonian limit.

We report on a new study of the velocity distribution in N-body simulations. We investigate the center-of-mass and internal kinetic energies of coarsening cells as a function of time, cell size and cell mass. By using self-similar cosmological models, we are able to derive theoretical predictions for comparison and to assess the influence of finite-size and resolution effects. The most interesting result is the discovery of a polytropic-like relationship between the average velocity dispersion (internal kinetic energy) and the mass density in an intermediate range of densities, K ~ rho^{2-eta}. The exponent eta measures the deviations from the virial prediction, eta_virial=0. For self-similar models, eta depends only on the spectral index of the initial power spectrum. We also study CDM models and confirm a previous result that the same polytropic-like dependence exists (astro-ph/0103313), with a time and coarsening length dependent eta. The dependence K(rho) is an important input for a recently proposed theoretical model of cosmological structure formation which improves over the standard dust model (pressureless fluid) by regularizing the density singularities.

This work reports on a study of the spatially coarse-grained velocity dispersion in cosmological N-body simulations (OCDM and Lambda-CDM models) as a function of time (redshifts z=0-4) and of the coarsening length (0.6-20 Mpc/h). The main result is the discovery of a polytropic relationship I_1 ~ rho^{2-eta} between the velocity-dispersion kinetic energy density of the coarsening cells, I_1, and their mass density, rho. The exponent eta, dependent on time and coarsening scale, is a compact measure of the deviations from the naive virial prediction eta_virial=0. This relationship supports the ``polytropic assumption'' which has been employed in theoretical models for the growth of cosmological structure by gravitational instability.

In astro-ph/0109199, the initial conditions (IC's) of cosmological N-body simulations by the Virgo Consortium are analyzed and it is concluded that the density fluctuations are rather different from the desired ones. We have repeated the analysis of the IC's using our own code and the code provided by the authors of astro-ph/0109199, obtaining results that disprove the criticisms.

We present the study of ten random realizations of a density field characterized by a cosmological power spectrum P(k) at redshift z=50. The reliability of such initial conditions for n-body simulations are tested with respect to their correlation properties. The power spectrum P(k), and the mass variance sigmaM(r) do not show detectable deviations from the desired behavior in the intermediate range of scales between the mean interparticle distance and the simulation volume. The estimator for xi(r) is too noisy to detect any reliable signal at the initial redshift z=50. The particle distributions are then evolved forward until z=0. This allows us to explore the cosmic variance stemming from the random nature of the initial conditions. With cosmic variance we mean the fact that a simulation represents a single realization of the stochastic initial conditions whereas the real Universe contains many realizations of regions of the size of the box; this problem affects most importantly the scales at about the fundamental mode. We study morphological descriptors of the matter distribution such as the genus, as well as the internal properties of the largest object(s) forming in the box. We find that the scatter is at least comparable to the scatter in the fundamental mode.

This work reports on the application of the Eulerian perturbation theory to a recently proposed model of cosmological structure formation by gravitational instability (astro-ph/0009414). Its physical meaning is discussed in detail and put in perspective of previous works. The model incorporates in a systematic fashion corrections to the popular dust model due to multistreaming and, more generally, the small-scale, virialized degrees of freedom. It features a time-dependent length scale L(t) estimated to be L/r0 ~ 0.1 (r0(t) is the nonlinear scale, at which the density variance =1). The model provides a new angle on the dust model and allows to overcome some of its limitations. Thus, the scale L(t) works as a physically meaningful short-distance cutoff for the divergences appearing in the perturbation expansion of the dust model when there is too much initial power on small scales. The model also incorporates the generation of vorticity by tidal forces; according to the perturbational result, the filtered vorticity for standard CDM initial conditions should be significant today only at scales below ~ 1 Mpc/h.

The effect of the large-scale cosmological expansion on small systems is studied in the light of modern cosmological models of large-scale structure. We identify certain assumptions of earlier works which render them unrealistic regarding these cosmological models. The question is reanalyzed by dropping these assumptions to conclude that a given small system can experience either an expansion or a contraction of cosmological origin.

The Minkowski functionals are a mathematical tool to quantify morphological features of patterns. Some applications to the matter distribution in galaxy catalogues and N-body simulations are reviewed, with an emphasis on the effects of cosmic variance. The conclusions are that (i) the observed large-scale morphology is sensitive to cosmic variance on scales much larger than the nonlinear length (approx. 8 Mpc/h), and (ii) the large-scale morphology predicted by simulations is thus affected by finite-size effects, but nonetheless a Lambda-CDM model is favored.

A hydrodynamic formulation of the evolution of large-scale structure in the Universe is presented. It relies on the spatially coarse-grained description of the dynamical evolution of a many-body gravitating system. Because of the assumed irrelevance of short-range (``collisional'') interactions, the way to tackle the hydrodynamic equations is essentially different from the usual case. The main assumption is that the influence of the small scales over the large-scale evolution is weak: this idea is implemented in the form of a large-scale expansion for the coarse-grained equations. This expansion builds a framework in which to derive in a controlled manner the popular ``dust'' model (as the lowest-order term) and the ``adhesion'' model (as the first-order correction). It provides a clear physical interpretation of the assumptions involved in these models and also the possibility to improve over them.

There is an ongoing debate in cosmology about the value of the length scale at which ``homogeneity'' in the matter distribution is reached or even if such a scale exists. In the wake of this debate, we intend in this letter to clarify the meaning of the statement transition to homogeneity and of the concept of correlation length. We show that there are two scales each associated to the fractal and to the homogeneous regimes of the matter distribution, respectively. The distinction between both scales has deep consequences, for there can be a regime which despite having small fluctuations around the average density exhibits large clusters of galaxies.

We propose a phenomenological generalization of the models of large-scale structure formation in the Universe by gravitational instability in two ways: we include pressure forces to model multi-streaming, and noise to model fluctuations due to neglected short-scale physical processes. We show that pressure gives rise to a viscous-like force of the same character as that one introduced in the ``adhesion model'', while noise leads to a roughening of the density field yielding a scaling behavior of its correlations.

We show how the interplay of non-linear dynamics, self-gravity, and fluctuations leads to self-affine behavior of matter density correlations quite generically, i.e., with a power-law exponent whose value does not depend in a very direct way on the specific details of the dynamics. This we do by means of the Renormalization Group, a powerful analytic tool for extracting asymptotic behavior of many-body systems

The generally held view that a model of large-scale structure, formed by collisionless matter in the Universe, can be based on the matter model ``dust'' fails in the presence of multi-stream flow, i.e., velocity dispersion. We argue that models for large-scale structure should rather be constructed for a flow which describes the average motion of a multi-stream system. We present a clearcut reasoning how to approach the problem and derive an evolution equation for the mean peculiar-velocity relative to background solutions of Friedmann-Lema\^\i tre type. We consider restrictions of the nonlinear problem and show that the effect of velocity dispersion gives rise to an effective viscosity of non-dissipative gravitational origin. We discuss subcases which arise naturally from this approach: the ``sticky particle model'' and the ``adhesion approximation''. We also construct a novel approximation that features adhesive action in the multi-stream regime while conserving momentum, which was considered a drawback of the standard approximation based on Burger's equation. We finally argue that the assumptions made to obtain these models should be relaxed and we discuss how this can be achieved.

As is well known, structure formation in the Universe at times after decoupling can be described by hydrodynamic equations. These are shown here to be equivalent to a generalization of the stochastic Kardar--Parisi--Zhang equation with time-- dependent viscosity in epochs of dissipation. As a consequence of the Dynamical Critical Scaling induced by noise and fluctuations, these equations describe the fractal behavior (with a scale dependent fractal dimension) observed at the smaller scales for the galaxy--to--galaxy correlation function and $also$ the Harrison--Zel'dovich spectrum at decoupling. By a Renormalization Group calculation of the two--point correlation function between galaxies in the presence of (i) the expansion of the Universe and (ii) non--equilibrium, we can account, from first principles, for the main features of the observed shape of the power spectrum.