Context. Variability caused by stellar activity represents a challenge to the discovery and characterization of terrestrial exoplanets and complicates the interpretation of atmospheric planetary signals. Aims. We aim to use a detailed modeling tool to reproduce the effect of active regions on radial velocity measurements, which aids the identification of the key parameters that have an impact on the induced variability. Methods. We analyzed the effect of stellar activity on radial velocities as a function of wavelength by simulating the impact of the properties of spots, shifts induced by convective motions, and rotation. We focused our modeling effort on the active star YZ CMi (GJ 285), which was photometrically and spectroscopically monitored with CARMENES and the Telescopi Joan Or\'o. Results. We demonstrate that radial velocity curves at different wavelengths yield determinations of key properties of active regions, including spot filling factor, temperature contrast, and location, thus solving the degeneracy between them. Most notably, our model is also sensitive to convective motions. Results indicate a reduced convective shift for M dwarfs when compared to solar-type stars (in agreement with theoretical extrapolations) and points to a small global convective redshift instead of blueshift. Conclusions. Using a novel approach based on simultaneous chromatic radial velocities and light curves, we can set strong constraints on stellar activity, including an elusive parameter such as the net convective motion effect.

We present $NuSTAR$ observations of four active galactic nuclei (AGN) located within 15 Mpc. These AGN, namely ESO 121-G6, NGC 660, NGC 3486 and NGC 5195, have observed X-ray luminosities of $L_{\rm 2-10\ keV, obs} \lesssim$ 10$^{39}$ erg s$^{-1}$, classifying them as low luminosity AGN (LLAGN). We perform broadband X-ray spectral analysis for the AGN by combining our $NuSTAR$ data with $Chandra$ or $XMM-Newton$ observations to directly measure their column densities ($N_{\rm H}$) and infer their intrinsic power. We complement our X-ray data with archival and new high angular resolution mid-infrared (mid-IR) data for all objects, except NGC 5195. Based on our X-ray spectral analysis, we found that both ESO 121-G6 and NGC 660 are heavily obscured ($N_{\rm H}$ > 10$^{23}$ cm$^{-2}$; $L_{\rm 2-10\ keV,\ int} \sim$ 10$^{41}$ erg s$^{-1}$), and NGC 660 may be Compton-thick. We also note that the X-ray flux and spectral slope for ESO 121-G6 have significantly changed over the last decade, indicating significant changes in the obscuration and potentially accretion rate. On the other hand, NGC 3486 and NGC 5195 appear to be unobscured and just mildly obscured, respectively, with $L_{\rm 2-10\ keV,\ int} <$ 10$^{39}$ erg s$^{-1}$; i.e., genuine LLAGN. Both of the heavily obscured AGN have $L_{\rm bol} >$ 10$^{41}$ erg s$^{-1}$ and $\lambda_{\rm Edd} \gtrsim$ 10$^{-3}$, and are detected in high angular resolution mid-IR imaging, indicating the presence of obscuring dust on nuclear scale. NGC 3486 however, is undetected in high-resolution mid-IR imaging, and the current data do not provide stringent constraints on the presence or absence of obscuring nuclear dust in the AGN.

On 2019 August 14, the LIGO and Virgo interferometers detected a high-significance event labelled S190814bv. Preliminary analysis of the GW data suggests that the event was likely due to the merger of a compact binary system formed by a BH and a NS. ElectromagNetic counterparts of GRAvitational wave sources at the VEry Large Telescope (ENGRAVE) collaboration members carried out an intensive multi-epoch, multi-instrument observational campaign to identify the possible optical/near infrared counterpart of the event. In addition, the ATLAS, GOTO, GRAWITA-VST, Pan-STARRS and VINROUGE projects also carried out a search on this event. Our observations allow us to place limits on the presence of any counterpart and discuss the implications for the kilonova (KN) possibly generated by this NS-BH merger, and for the strategy of future searches. Altogether, our observations allow us to exclude a KN with large ejecta mass $M\gtrsim 0.1\,\mathrm{M_\odot}$ to a high ($>90\%$) confidence, and we can exclude much smaller masses in a subsample of our observations. This disfavours the tidal disruption of the neutron star during the merger. Despite the sensitive instruments involved in the campaign, given the distance of S190814bv we could not reach sufficiently deep limits to constrain a KN comparable in luminosity to AT 2017gfo on a large fraction of the localisation probability. This suggests that future (likely common) events at a few hundreds Mpc will be detected only by large facilities with both high sensitivity and large field of view. Galaxy-targeted observations can reach the needed depth over a relevant portion of the localisation probability with a smaller investment of resources, but the number of galaxies to be targeted in order to get a fairly complete coverage is large, even in the case of a localisation as good as that of this event.

One distinctive feature of low-luminosity active galactic nuclei (LLAGN) is the relatively weak reflection features they may display in the X-ray spectrum, which can result from the disappearance of the torus with decreasing accretion rates. Some material, however, must surround the active nucleus, i.e., the accretion flow itself and, possibly, a flattened-out or thinned torus. In this work, we study whether reflection is indeed absent or undetectable due to its intrinsically weak features together with the low statistics inherent to LLAGN. Here we focus on NGC 3718 ($L/L_{\rm Edd}\sim10^{-5}$) combining observations from XMM--Newton and the deepest to date NuSTAR (0.5--79 keV) spectrum of a LLAGN, to constrain potential reflectors, and analyze how the fitted coronal parameters depend on the reflection model. We test models representing both an accretion disc (Relxill) and a torus-like (MYTorus and Borus) neutral reflector. From a statistical point of view, reflection is not required, but its inclusion allows to place strong constraints on the geometry and physical features of the surroundings: both neutral reflectors (torus) tested should be Compton thin ($N_H<10^{23.2}$cm$^{-2}$) and preferentially cover a large fraction of the sky. If the reflected light instead arises from an ionized reflector, a highly ionized case is preferred. These models produce an intrinsic power-law spectral index in the range [1.81--1.87], where the torus models result in steeper slopes. The cut-off energy of the power-law emission also changes with the inclusion of reflection models, resulting in constrained values for the disc reflectors and unconstrained values for torus reflectors.

The He I infrared (IR) triplet at 10830 \AA is known as an activity indicator in solar-type stars and has become a primary diagnostic in exoplanetary transmission spectroscopy. He I lines are a tracer of the stellar extreme-ultraviolet irradiation from the transition region and corona. We study the variability of the He I IR triplet lines in a spectral time series of 319 M~dwarf stars that was obtained with the CARMENES high-resolution optical and near-infrared spectrograph at Calar Alto. We detect He I IR line variability in 18% of our sample stars, all of which show H$\alpha$ in emission. Therefore, we find detectable He I IR variability in 78% of the sub-sample of stars with H$\alpha$ emission. Detectable variability is strongly concentrated in the latest spectral sub-types, where the He I IR lines during quiescence are typically weak. The fraction of stars with detectable He I IR variation remains lower than 10% for stars earlier than M3.0 V, while it exceeds 30% for the later spectral sub-types. Flares are accompanied by particularly pronounced line variations, including strongly broadened lines with red and blue asymmetries. However, we also find evidence for enhanced He I IR absorption, which is potentially associated with increased high-energy irradiation levels at flare onset. Generally, He I IR and H$\alpha$ line variations tend to be correlated, with H$\alpha$ being the most sensitive indicator in terms of pseudo-equivalent width variation. This makes the He I IR triplet a favourable target for planetary transmission spectroscopy.

We present a characterization of the physical properties of a sample of 35 securely-detected, dusty galaxies in the deep ALMA 1.2-mm image obtained as part of the ALMA Spectroscopic Survey in the {\it Hubble} Ultra Deep Field (ASPECS) Large Program. This sample is complemented by 26 additional sources identified via an optical/infrared source positional prior. Using their well-characterized spectral energy distributions, we derive median stellar masses and star formation rates (SFR) of $4.8\times10^{10}~M_\odot$ and 30 $M_\odot$ yr$^{-1}$, and interquartile ranges of $(2.4-11.7)\times10^{10}~M_\odot$ and $20-50~M_\odot$ yr$^{-1}$, respectively. We derive a median spectroscopic redshift of 1.8 with an interquartile range $1.1-2.6$, significantly lower than sub-millimeter galaxies detected in shallower, wide-field surveys. We find that 59\%$\pm$13\%, 6\%$\pm$4\%, and 34\%$\pm$9\% of our sources are within, above and below $\pm0.4$ dex from the SFR- stellar mass relation or main-sequence (MS), respectively. The ASPECS galaxies closely follow the SFR-molecular gas mass relation and other previously established scaling relations, confirming a factor of five increase of the gas-to-stellar mass ratio from $z=0.5$ to $z=2.5$, and a mild evolution of the gas depletion timescales with a typical value of 0.7 Gyr at $z=1-3$. ASPECS galaxies located significantly below the MS, a poorly exploited parameter space, have low gas-to-stellar mass ratios of $\sim0.1-0.2$ and long depletion timescales $>1$ Gyr. Galaxies along the MS dominate the cosmic density of molecular gas at all redshifts. Systems above the main sequence have an increasing contribution to the total gas reservoirs from $z<1$ to $z=2.5$, while the opposite is found for galaxies below the MS.

Hot Dust-Obscured Galaxies (Hot DOGs) are among the most luminous galaxies in the Universe. Powered by highly obscured, possibly Compton-thick, active galactic nuclei (AGNs), Hot DOGs are characterized by SEDs that are very red in the mid-IR yet dominated by the host galaxy stellar emission in the UV and optical. An earlier study identified a sub-sample of Hot DOGs with significantly enhanced UV emission. One target, W0204-0506, was studied in detail and, based on Chandra observations, it was concluded that the enhanced emission was most likely due to either extreme unobscured star-formation (${\rm SFR}>1000~M_{\odot}~\rm yr^{-1}$) or to light from the highly obscured AGN scattered by gas or dust into our line of sight. Here, we present a follow-up study of W0204-0506 as well as two more Hot DOGs with excess UV emission. For the two new objects we obtained Chandra/ACIS-S observations, and for all three targets we obtained HST/WFC3 F555W and F160W imaging. We conclude that the excess UV emission is primarily dominated by light from the central highly obscured, hyper-luminous AGN that has been scattered into our line of sight. We cannot rule out, however, that star-formation may significantly contribute to the UV excess of W0204-0506.

The high-resolution, dual channel, visible and near-infrared spectrograph CARMENES offers exciting opportunities for stellar and exoplanetary research on M dwarfs. In this work we address the challenge of reaching the highest radial velocity precision possible with a complex, actively cooled, cryogenic instrument, such as the near-infrared channel. We describe the performance of the instrument and the work flow used to derive precise Doppler measurements from the spectra. The capability of both CARMENES channels to detect small exoplanets is demonstrated with the example of the nearby M5.0 V star CD Cet (GJ 1057), around which we announce a super-Earth ($4.0\pm0.4\,M_\oplus$) companion on a 2.29 d orbit.

The He I infrared (IR) line at a vacuum wavelength of 10833 A is a diagnostic for the investigation of atmospheres of stars and planets orbiting them. For the first time, we study the behavior of the He I IR line in a set of chromospheric models for M-dwarf stars, whose much denser chromospheres may favor collisions for the level population over photoionization and recombination, which are believed to be dominant in solar-type stars. For this purpose, we use published PHOENIX models for stars of spectral types M2 V and M3 V and also compute new series of models with different levels of activity following an ansatz developed for the case of the Sun. We perform a detailed analysis of the behavior of the He I IR line within these models. We evaluate the line in relation to other chromospheric lines and also the influence of the extreme ultraviolet (EUV) radiation field. The analysis of the He I IR line strengths as a function of the respective EUV radiation field strengths suggests that the mechanism of photoionization and recombination is necessary to form the line for inactive models, while collisions start to play a role in our most active models. Moreover, the published model set, which is optimized in the ranges of the Na I D2, H$\alpha$, and the bluest Ca II IR triplet line, gives an adequate prediction of the He I IR line for most stars of the stellar sample. Because especially the most inactive stars with weak He I IR lines are fit worst by our models, it seems that our assumption of a 100% filling factor of a single inactive component no longer holds for these stars.

Enhanced emission from the dense gas tracer HCN (relative to HCO$^+ $) has been proposed as a signature of active galactic nuclei (AGN). In a previous single-dish millimeter line survey we identified galaxies with HCN/HCO$ ^+ $ (1-0) intensity ratios consistent with those of many AGN but whose mid-infrared spectral diagnostics are consistent with little to no ( $\lesssim15\% $) contribution of an AGN to the bolometric luminosity. To search for putative heavily obscured AGN, we present and analyze \nustar hard X-ray (3-79 keV) observations of four such galaxies from the Great Observatories All-sky LIRG Survey. We find no X-ray evidence for AGN in three of the systems and place strong upper limits on the energetic contribution of any heavily obscured ($N_{\rm H}>10^{24}$ cm$^{-2}$) AGN to their bolometric luminosity. The X-ray flux upper limits are presently an order of magnitude below what XDR-driven chemistry model predict are necessary to drive HCN enhancements. In a fourth system we find a hard X-ray excess consistent with the presence of an AGN, but contributing only $\sim3\%$ of the bolometric luminosity. It is also unclear if the AGN is spatially associated with the HCN enhancement. We further explore the relationship between HCN/HCO$^+$ (for several $\mathrm{J}_\mathrm{upper}$ levels) and $L_\mathrm{AGN}/L_\mathrm{IR}$ for a larger sample of systems in the literature. We find no evidence for correlations between the line ratios and the AGN fraction derived from X-rays, indicating that HCN/HCO$^+$ intensity ratios are not driven by the energetic dominance of AGN, nor are they reliable indicators of whether SMBH accretion is ongoing.

We report on radial velocity time series for two M0.0V stars, GJ338B and GJ338A, using the CARMENES spectrograph, complemented by ground-telescope photometry from Las Cumbres and Sierra Nevada observatories. We aim to explore the presence of small planets in tight orbits using the spectroscopic radial velocity technique. We obtained 159 and 70 radial velocity measurements of GJ338B and A, respectively, with the CARMENES visible channel. We also compiled additional relative radial velocity measurements from the literature and a collection of astrometric data that cover 200 a of observations to solve for the binary orbit. We found dynamical masses of 0.64$\pm$0.07M$_\odot$ for GJ338B and 0.69$\pm$0.07M$_\odot$ for GJ338A. The CARMENES radial velocity periodograms show significant peaks at 16.61$\pm$0.04 d (GJ338B) and 16.3$^{+3.5}_{-1.3}$ d (GJ338A), which have counterparts at the same frequencies in CARMENES activity indicators and photometric light curves. We attribute these to stellar rotation. GJ338B shows two additional, significant signals at 8.27$\pm$0.01 and 24.45$\pm$0.02 d, with no obvious counterparts in the stellar activity indices. The former is likely the first harmonic of the star's rotation, while we ascribe the latter to the existence of a super-Earth planet with a minimum mass of 10.27$^{+1.47}_{-1.38}$$M_{\oplus}$ orbiting GJ338B. GJ338B b lies inside the inner boundary of the habitable zone around its parent star. It is one of the least massive planets ever found around any member of stellar binaries. The masses, spectral types, brightnesses, and even the rotational periods are very similar for both stars, which are likely coeval and formed from the same molecular cloud, yet they differ in the architecture of their planetary systems.

For years, the standard procedure to measure radial velocities (RVs) of spectral observations consisted in cross-correlating the spectra with a binary mask, that is, a simple stellar template that contains information on the position and strength of stellar absorption lines. The cross-correlation function (CCF) profiles also provide several indicators of stellar activity. We present a methodology to first build weighted binary masks and, second, to compute the CCF of spectral observations with these masks from which we derive radial velocities and activity indicators. These methods are implemented in a python code that is publicly available. To build the masks, we selected a large number of sharp absorption lines based on the profile of the minima present in high signal-to-noise ratio (S/N) spectrum templates built from observations of reference stars. We computed the CCFs of observed spectra and derived RVs and the following three standard activity indicators: full-width-at-half-maximum as well as contrast and bisector inverse slope.We applied our methodology to CARMENES high-resolution spectra and obtain RV and activity indicator time series of more than 300 M dwarf stars observed for the main CARMENES survey. Compared with the standard CARMENES template matching pipeline, in general we obtain more precise RVs in the cases where the template used in the standard pipeline did not have enough S/N. We also show the behaviour of the three activity indicators for the active star YZ CMi and estimate the absolute RV of the M dwarfs analysed using the CCF RVs.

HD 209458b is an exoplanet with an upper atmosphere undergoing blow-off escape that has mainly been studied using measurements of the Ly-alpha absorption. Recently, high-resolution measurements of absorption in the He I triplet line at 10830 angstroms of several exoplanets (including HD 209458b) have been reported, creating a new opportunity to probe escaping atmospheres. We aim to better understand the atmospheric regions of HD 209458b from where the escape originates. We developed a 1D hydrodynamic model with spherical symmetry for the HD 209458 b thermosphere coupled with a non-local thermodynamic model for the population of the He triplet state. In addition, we performed high-resolution radiative transfer calculations of synthetic spectra for the He triplet lines and compared them with the measured absorption spectrum in order to retrieve information about the atmospheric parameters. We find that the measured spectrum constrains the [H]/[H$^{+}$] transition altitude occurring in the range of 1.2 to 1.9Rp. H is almost fully ionised at altitudes above 2.9Rp. We also find that the X-ray and EUV absorption takes place at effective radii from 1.16 to 1.30Rp, and that the He triplet peak density occurs at altitudes from 1.04 to 1.60Rp. Additionally, the averaged mmw is confined to the 0.61-0.73 g/mole interval, and the thermospheric H/He ratio should be larger than 90/10, and most likely approximately 98/2. We also provide a one-to-one relationship between mass-loss rate and temperature. Based on the energy-limited escape approach and assuming heating efficiencies of 0.1-0.2, we find a mass-loss rate in the range of (0.42-1.00)$\times 10^{11}$ g/s and a corresponding temperature range of 7125 to 8125K. The analysis of the measured He triplet absorption spectrum significantly constrains the thermospheric structure of HD 209458b and advances our knowledge of its escaping atmosphere.

We present the discovery and characterisation of two transiting planets observed by the Transiting Exoplanet Survey Satellite (TESS) orbiting the nearby (d ~ 22 pc), bright (J ~ 9 mag) M3.5 dwarf LP 729-54 (TOI-732). We confirm both planets and their association with LP 729-54 via ground-based photometry and determine their masses using precise radial velocities measured with the CARMENES spectrograph. Precise stellar parameters determined from CARMENES high resolution spectra confirm that LP 729-54 is a mid-M dwarf with an effective temperature of T_eff = 3360 +\- 51 K, a surface gravity of log(g) = 4.81 +/- 0.04 (cgs), and an iron abundance of [Fe/H] = 0.09 +/- 0.16 dex, with an inferred mass of M_star = 0.379 +/- 0.016 M_sun and a radius of R_star = 0.382 +/- 0.012 R_sun. The ultra-short period planet LP 729-54 b (P_b = 0.77 d) with a radius of 1.35^{+0.06}_{-0.06} R_earth, a mass of 2.34^{+0.24}_{-0.23} M_earth, and a bulk density of 5.24^{+0.94}_{-0.81} g cm^{-3} joins the population of Earth-size planets with rocky, terrestrial composition. The outer planet, LP 729-54 c, with an orbital period of 12.25 d, radius of 2.42^{+0.10}_{-0.10} R_earth, mass of 6.29^{+0.63}_{-0.61} M_earth, and mean density of 2.45^{+0.44}_{-0.37} g cm^{-3} belongs to the population of dense sub-Neptunes. With the two planets located on opposite sides of the radius gap, this planetary system is an excellent target for testing planetary formation, evolution and atmospheric models. In particular, LP 729-54 c is an ideal object for atmospheric studies with the James Webb Space Telescope.

Using the deepest 1.2 mm continuum map to date in the Hubble Ultra Deep Field obtained as part of the ALMA Spectroscopic Survey (ASPECS) large program, we measure the cosmic density of dust and implied gas (H$_{2}+$H I) mass in galaxies as a function of look-back time. We do so by stacking the contribution from all $H$-band selected galaxies above a given stellar mass in distinct redshift bins, $\rho_{\rm dust}(M_\ast>M,z)$ and $\rho_{\rm gas}(M_\ast>M,z)$. At all redshifts, $\rho_{\rm dust}(M_\ast>M,z)$ and $\rho_{\rm gas}(M_\ast>M,z)$ grow rapidly as $M$ decreases down to $10^{10}\,M_\odot$, but this growth slows down towards lower stellar masses. This flattening implies that at our stellar mass-completeness limits ($10^8\,M_\odot$ and $10^{8.9}\,M_\odot$ at $z\sim0.4$ and $z\sim3$), both quantities converge towards the total cosmic dust and gas mass densities in galaxies. The cosmic dust and gas mass densities increase at early cosmic time, peak around $z\sim2$, and decrease by a factor $\sim4$ and 7, compared to the density of dust and molecular gas in the local universe, respectively. The contribution of quiescent galaxies -- i.e., with little on-going star-formation -- to the cosmic dust and gas mass densities is minor ($\lesssim10\%$). The redshift evolution of the cosmic gas mass density resembles that of the star-formation rate density, as previously found by CO-based measurements. This confirms that galaxies have relatively constant star-formation efficiencies (within a factor $\sim2$) across cosmic time. Our results also imply that by $z\sim0$, a large fraction ($\sim90\%$) of dust formed in galaxies across cosmic time has been destroyed or ejected to the intergalactic medium.

We present a new semi-empirical model for the dust continuum number counts of galaxies at 1.1 millimeter and 850 \micron. Our approach couples an observationally motivated model for the stellar mass and SFR distribution of galaxies with empirical scaling relations to predict the dust continuum flux density of these galaxies. Without a need to tweak the IMF, the model reproduces the currently available observations of the 1.1 millimeter and 850 \micron number counts, including the observed flattening in the 1.1 millimeter number counts below 0.3 mJy \citep{Gonzalez2019numbercounts} and the number counts in discrete bins of different galaxy properties. Predictions of our work include : (1) the galaxies that dominate the number counts at flux densities below 1 mJy (3 mJy) at 1.1 millimeter (850 $\mu$m) have redshifts between $z=1$ and $z=2$, stellar masses of $\sim 5\times10^{10}~\rm{M}_\odot$, and dust masses of $\sim 10^{8}~\rm{M}_\odot$; (2) the flattening in the observed 1.1 millimeter number counts corresponds to the knee of the 1.1 millimeter luminosity function. A similar flattening is predicted for the number counts at 850 $\mu$m; (3) the model reproduces the redshift distribution of current 1.1 millimeter detections; (4) to efficiently detect large numbers of galaxies through their dust continuum, future surveys should scan large areas once reaching a 1.1 millimeter flux density of 0.1 mJy rather than integrating to fainter fluxes. Our modeling framework also suggests that the amount of information on galaxy physics that can be extracted from the 1.1 millimeter and 850 $\mu$m number counts is almost exhausted.

We present the results from the 1.2 mm continuum image obtained as part of the ALMA Spectroscopic Survey in the Hubble Ultra Deep Field (ASPECS). The 1.2 mm continuum image has a size of 2.9 (4.2) arcmin$^2$ within a primary beam response of 50% (10%) and a rms value of $9.3\thinspace{\rm\mu Jy\thinspace beam^{-1}}$. We detect 35 sources at high significance (Fidelity $\geq0.5$), 32 of these have well characterized near-infrared HST counterparts. We estimate the 1.2 mm number counts to flux levels of $<30\thinspace{\rm\mu Jy}$ in two different ways: we first use the detected sources to constrain the number counts and find a significant flattening of the counts below $S_\nu \sim 0.1$ mJy. In a second approach, we constrain the number counts by using a probability of deflection statistics (P(D)) analysis. For this latter approach, we describe new methods to accurately measure the noise in interferometric imaging (employing jack-knifing in the cube and in the visibility plane). This independent measurement confirms the flattening of the number counts. Our analysis of the differential number counts shows that we are detecting $\sim$93% ($\sim$100% if we include the lower fidelity detections) of the total continuum dust emission associated to galaxies in the HUDF. The ancillary data allows us to study the dependence of the 1.2 mm number counts on redshift ($z=0-4$), galaxy dust mass (${\rm M}_{\rm dust}=10^{7}-10^{9}{\rm M}_{\odot}$), stellar mass (${\rm M}_{*}=10^{9}-10^{12}{\rm M}_{\odot}$), and star-formation rate (${\rm SFR}=1-1000\thinspace{\rm M}_{\odot}\thinspace{\rm yr^{-1}}$). In an accompanying paper we show that the number counts are crucial to constrain galaxy evolution models and the understanding of star-forming galaxies at high redshift.

The nearby ultra-compact multiplanetary system YZ Ceti consists of at least three planets. The orbital period of each planet is the subject of discussion in the literature due to strong aliasing in the radial velocity data. The stellar activity of this M dwarf also hampers significantly the derivation of the planetary parameters. With an additional 229 radial velocity measurements obtained since the discovery publication, we reanalyze the YZ Ceti system and resolve the alias issues. We use model comparison in the framework of Bayesian statistics and periodogram simulations based on a method by Dawson and Fabrycky to resolve the aliases. We discuss additional signals in the RV data, and derive the planetary parameters by simultaneously modeling the stellar activity with a Gaussian process regression model. To constrain the planetary parameters further we apply a stability analysis on our ensemble of Keplerian fits. We resolve the aliases: the three planets orbit the star with periods of $2.02$ d, $3.06$ d, and $4.66$ d. We also investigate an effect of the stellar rotational signal on the derivation of the planetary parameters, in particular the eccentricity of the innermost planet. Using photometry we determine the stellar rotational period to be close to $68$ d. From the absence of a transit event with TESS, we derive an upper limit of the inclination of $i_\mathrm{max} = 87.43$ deg. YZ Ceti is a prime example of a system where strong aliasing hindered the determination of the orbital periods of exoplanets. Additionally, stellar activity influences the derivation of planetary parameters and modeling them correctly is important for the reliable estimation of the orbital parameters in this specific compact system. Stability considerations then allow additional constraints to be placed on the planetary parameters.

Context. GJ 1148 is an M-dwarf star hosting a planetary system composed of two Saturn-mass planets in eccentric orbits with periods of 41.38 and 532.02 days. Aims. We reanalyze the orbital configuration and dynamics of the GJ 1148 multi-planetary system based on new precise radial velocity (RV) measurements taken with CARMENES. Methods. We combined new and archival precise Doppler measurements from CARMENES with those available from HIRES for GJ 1148 and modeled these data with a self-consistent dynamical model. We studied the orbital dynamics of the system using the secular theory and direct N-body integrations. The prospects of potentially habitable moons around GJ 1148 b were examined. Results. The refined dynamical analyses show that the GJ 1148 system is long-term stable in a large phase-space of orbital parameters with an orbital configuration suggesting apsidal alignment, but not in any particular high-order mean-motion resonant commensurability. GJ 1148 b orbits inside the optimistic habitable zone (HZ). We find only a narrow stability region around the planet where exomoons can exist. However, in this stable region exomoons exhibit quick orbital decay due to tidal interaction with the planet. Conclusions. The GJ 1148 planetary system is a very rare M-dwarf planetary system consisting of a pair of gas giants, the inner of which resides in the HZ. We conclude that habitable exomoons around GJ 1148 b are very unlikely to exist.

The CDF-S is the deepest X-ray image available and will remain so for the near future. We provide a spectroscopic (64.5%; 64% with spectral classifications) and photometric redshift catalog for the full 7 Ms sample, but much of our analysis focuses on the central (off-axis angles <5.7') region, which contains a large, faint ALMA sample of 75 >4.5-sigma 850 micron sources. We measure the 850 micron fluxes at the X-ray positions using the ALMA images, where available, or an ultradeep SCUBA-2 map. We find that the full X-ray sample produces ~10% of the 850 micron extragalactic background light. We separate the submillimeter detected X-ray sources into star-forming galaxies and AGNs using a star formation rate (SFR) versus X-ray luminosity calibration for high SFR galaxies. We confirm this separation using the X-ray photon indices. We measure the X-ray fluxes at the accurate positions of the 75 ALMA sources and detect 70% at >3-sigma in either the 0.5-2 or 2-7 keV bands. However, many of these may produce both their X-ray and submillimeter emission by star formation. Indeed, we find that only 20% of the ALMA sources have intermediate X-ray luminosities (rest-frame 8-28 keV luminosities of 10^42.5-10^44 erg/s), and none has a high X-ray luminosity (>10^44 erg/s). Conversely, after combining the CDF-S with the CDF-N, we find extreme star formation (SFR>300 solar masses per yr) in some intermediate X-ray luminosity sources but not in any high X-ray luminosity sources. We argue that the quenching of star formation in the most luminous AGNs may be a consequence of the clearing of gas in these sources.

We search for high-redshift (z>4.5) X-ray AGNs in the deep central (off-axis angle <5.7') region of the 7 Ms Chandra Deep Field-South X-ray image. We compile an initial candidate sample from direct X-ray detections. We then probe more deeply in the X-ray data by using pre-selected samples with high spatial resolution NIR/MIR (HST 1.6 micron and Spitzer 4.5 micron) and submillimeter (ALMA 850 micron) observations. The combination of the NIR/MIR and submillimeter pre-selections allows us to find X-ray sources with a wide range of dust properties and spectral energy distributions (SEDs). We use the SEDs from the optical to the submillimeter to determine if previous photometric redshifts were plausible. Only five possible z>5 X-ray AGNs are found, all of which might also lie at lower redshifts. If they do lie at high redshifts, then two are Compton-thick AGNs, and three are ALMA 850 micron sources. We find that (i) the number density of X-ray AGNs is dropping rapidly at high redshifts, (ii) the detected AGNs do not contribute significantly to the photoionization at z>5, and (iii) the measured X-ray light density over z=5-10 implies a very low black hole accretion density with very little growth in the black hole mass density in this redshift range.

We present the results of the latest NuSTAR monitoring campaign of the Compton-thick Seyfert 2 galaxy NGC 1068, composed of four $\sim$50 ks observations performed between July 2017 and February 2018 to search for flux and spectral variability on timescales from 1 to 6 months. We detect one unveiling and one eclipsing event with timescales less than 27 and 91 days, respectively, ascribed to Compton-thick material with $N_H=(1.8\pm0.8)\times10^{24}$ cm$^{-2}$ and $N_H\geq(2.4\pm0.5)\times10^{24}$ cm$^{-2}$ moving across our line of sight. This gas is likely located in the innermost part of the torus or even further inward, thus providing further evidence of the clumpy structure of the circumnuclear matter in this source. Taking advantage of simultaneous Swift-XRT observations, we also detected a new flaring ULX, at a distance $d\sim$30" (i.e. $\sim$2 kpc) from the nuclear region of NGC 1068, with a peak X-ray intrinsic luminosity of $(3.0\pm0.4)\times10^{40}$ erg s$^{-1}$ in the 2-10 keV band.

Variability has proven to be a powerful tool to detect active galactic nuclei (AGN) in multi-epoch surveys. The new-generation facilities expected to become operational in the next few years will mark a new era in time-domain astronomy and their wide-field multi-epoch campaigns will favor extensive variability studies. We present our analysis of AGN variability in the second half of the VST survey of the Wide Chandra Deep Field South (W-CDF-S), performed in the r band and covering a 2 sq. deg. area. The analysis complements a previous work, in which the first half of the area was investigated. We provide a reliable catalog of variable AGN candidates, which will be critical targets in future variability studies. We selected a sample of optically variable sources and made use of infrared data from the Spitzer mission to validate their nature by means of color-based diagnostics. We obtain a sample of 782 AGN candidates among which 12 are classified as supernovae, 54 as stars, and 232 as AGN. We estimate a contamination $\lesssim 20\%$ and a completeness $\sim 38\%$ with respect to mid-infrared selected samples.

We have conducted 22 GHz radio imaging at 1" resolution of 100 low-redshift AGN selected at 14-195 keV by the Swift-BAT. We find a radio core detection fraction of 96%, much higher than lower-frequency radio surveys. Of the 96 radio-detected AGN, 55 have compact morphologies, 30 have morphologies consistent with nuclear star formation, and 11 have sub-kpc to kpc-scale jets. We find that the total radio power does not distinguish between nuclear star formation and jets as the origin of the radio emission. For 87 objects, we use optical spectroscopy to test whether AGN physical parameters are distinct between radio morphological types. We find that X-ray luminosities tend to be higher if the 22 GHz morphology is jet-like, but find no significant difference in other physical parameters. We find that the relationship between the X-ray and core radio luminosities is consistent with the $L_R/L_X \sim 10^{-5}$ of coronally active stars. We further find that the canonical fundamental planes of black hole activity systematically over-predict our radio luminosities, particularly for objects with star formation morphologies.

We present the highest resolution --- 15 pc (0.03'') --- ALMA $^{12}$CO(2-1) line emission and 1.3mm continuum maps, tracers of the molecular gas and dust, respectively, in the nearby merging galaxy system NGC 6240, that hosts two supermassive black holes growing simultaneously. These observations provide an excellent spatial match to existing Hubble optical and near-infrared observations of this system. A significant molecular gas mass, $\sim$9$\times$10$^9$M$_\odot$, is located in between the two nuclei, forming a clumpy stream kinematically dominated by turbulence, rather than a smooth rotating disk as previously assumed from lower resolution data. Evidence for rotation is seen in the gas surrounding the southern nucleus, but not in the northern one. Dynamical shells can be seen, likely associated with nuclear supernovae remnants. We further detect the presence of significant high velocity outflows, some of them reaching velocities $>$500 km/s, affecting a significant fraction, $\sim$11\% of the molecular gas in the nuclear region. Inside the spheres of influence of the northern and southern supermassive black holes we find molecular masses of 7.4$\times$10$^8$M$_\odot$ and 3.3$\times$10$^9$M$_\odot$, respectively. We are thus directly imaging the reservoir of gas that can accrete onto each supermassive black hole. These new ALMA maps highlight the critical need for high resolution observations of molecular gas in order to understand the feeding of supermassive black holes and its connection to galaxy evolution in the context of a major galaxy merger.

The Hubble Frontier Fields offer an exceptionally deep window into the high-redshift universe, covering a substantially larger area than the Hubble Ultra-Deep field at low magnification and probing 1--2 mags deeper in exceptional high-magnification regions. We aim to leverage Atacama Large Millimetre Array (ALMA) band 6 ($\approx$263\,GHz) mosaics in the central portions of five Frontier Fields to characterize the infrared (IR) properties of $1582$ ultraviolet (UV)-selected Lyman-Break Galaxies (LBGs) at redshifts of $z {\sim}$2--8. We investigated individual and stacked fluxes and IR excess (IRX) values of the LBG sample as functions of stellar mass ($\mathrm{M}_{\bigstar}$), redshift, UV luminosity and slope $\beta$, and lensing magnification. Two (2) LBG candidates were individually detected above a significance of $4.1{-}\sigma$, while stacked samples of the remaining LBG candidates yielded no significant detections. We investigated our detections and upper limits in the context of the IRX-$\mathrm{M}_{\bigstar}$ and IRX-$\beta$ relations, probing at least one dex lower in stellar mass than past studies have done. Our upper limits exclude substantial portions of parameter space and they are sufficiently deep in a handful of cases to create mild tension with the typically assumed attenuation and consensus relations. We observe a clear and smooth trend between $\mathrm{M}_{\bigstar}$ and $\beta$, which extends to low masses and blue (low) $\beta$ values, consistent with expectations from previous works.

We present VLT/XSHOOTER rest-frame UV-optical spectra of 10 Hot Dust-Obscured Galaxies (Hot DOGs) at $z\sim2$ to investigate AGN diagnostics and to assess the presence and effect of ionized gas outflows. Most Hot DOGs in this sample are narrow-line dominated AGN (type 1.8 or higher), and have higher Balmer decrements than typical type 2 quasars. Almost all (8/9) sources show evidence for ionized gas outflows in the form of broad and blueshifted [O III] profiles, and some sources have such profiles in H$\alpha$ (5/7) or [O II] (3/6). Combined with the literature, these results support additional sources of obscuration beyond the simple torus invoked by AGN unification models. Outflow rates derived from the broad [O III] line ($\rm \gtrsim10^{3}\,M_{\odot}\,yr^{-1}$) are greater than the black hole accretion and star formation rates, with feedback efficiencies ($\sim0.1-1\%$) consistent with negative feedback to the host galaxy's star formation in merger-driven quasar activity scenarios. We find the broad emission lines in luminous, obscured quasars are often better explained by outflows within the narrow line region, and caution that black hole mass estimates for such sources in the literature may have substantial uncertainty. Regardless, we find lower bounds on the Eddington ratio for Hot DOGs near unity.

Two bright X-ray transients were reported from the Chandra Deep Field South archival data, namely CDF-S XT1 and XT2. Whereas the nature of the former is not identified, the latter was suggested as an excellent candidate for a rapidly spinning magnetar born from a binary neutron star (BNS) merger. Here we propose a unified model to interpret both transients within the framework of the BNS merger magnetar model. According to our picture, CDF-S XT2 is observed from the "free zone" where the magnetar spindown powered X-ray emission escapes freely, whereas CDF-S XT1 originates from the "trapped zone" where the X-ray emission is initially blocked by the dynamical ejecta and becomes transparent after the ejecta is pushed to a distance where Thomson optical depth drops below unity. We fit the magnetar model to the light curves of both transients and derived consistent parameters for the two events, with magnetic field, initial spin period and X-ray emission efficiency being ($B_p=10^{16}\,G$, $P=1.2\,\rm ms$, $\eta = 0.001$) and ($B_p=10^{15.8}\,G$, $P=4.4\, \rm ms$, $\eta = 0.001$) for XT1 and XT2, respectively. The "isotropic equivalent" ejecta mass of XT1 is $M_{\rm ej} \sim 10^{-3}$ $M_{\odot}$, while it is not constrained for XT2. Our results suggest that more extreme magnetar parameters are required to have XT1 detected from the trapped zone. The model parameters for both events are generally consistent with those derived from SGRB X-ray plateau observations. The host galaxy properties of both transients are also consistent with those of SGRBs. The event rate densities of both XT1 and XT2 are consistent with that of BNS mergers.

The properties of the molecular gas can shed light on the physical conditions of quasar host galaxies and the effect of feedback from accreting supermassive black holes. We present a new CO(2-1) survey of 23 z<0.1 Palomar-Green quasars conducted with the Atacama Large Millimeter/submillimeter Array. CO emission was successfully detected in 91% (21/23) of the objects, from which we derive CO luminosities, molecular gas masses, and velocity line widths. Together with CO(1-0) measurements in the literature for 32 quasars (detection rate 53%), there are 15 quasars with both CO(1-0) and CO(2-1) measurements and in total 40 sources with CO measurements. We find that the line ratio R_21 = L'_CO(2-1)/L'_CO(1-0) is subthermal, broadly consistent with nearby galaxies and other quasars previously studied. No clear correlation is found between R_21 and the intensity of the interstellar radiation field or the luminosity of the active nucleus. As with the general galaxy population, quasar host galaxies exhibit a strong, tight, linear L_IR-L_CO relation, with a normalization consistent with that of starburst systems. We investigate the molecular-to-total gas mass fraction with the aid of total gas masses inferred from dust masses previously derived from infrared observations. Although the scatter is considerable, the current data do not suggest that the CO-to-H_2 conversion factor of quasar host galaxies significantly differs from that of normal star-forming galaxies.

We present a detailed study of ionized outflows in a large sample of ~650 hard X-ray detected AGN. Using optical spectroscopy from the BAT AGN Spectroscopic Survey (BASS) we are able to reveal the faint wings of the [OIII] emission lines associated with outflows covering, for the first time, an unexplored range of low AGN bolometric luminosity at low redshift (z~0.05). We test if and how the incidence and velocity of ionized outflow is related to AGN physical parameters: black hole mass, gas column density, Eddington Ratio, [OIII], X-ray, and bolometric luminosities. We find a higher occurrence of ionized outflows in type 1.9 (55%) and type 1 AGN (46%) with respect to type 2 AGN (24%). While outflows in type 2 AGN are evenly balanced between blue and red velocity offsets with respect to the [OIII] narrow component, they are almost exclusively blueshifted in type 1 and type 1.9 AGN. We observe a significant dependence between the outflow occurrence and accretion rate, which becomes relevant at high Eddington ratios (> -1.7). We interpret such behaviour in the framework of covering factor-Eddington ratio dependence. We don't find strong trends of the outflow maximum velocity with AGN physical parameters, as an increase with bolometric luminosity can be only identified when including samples of AGN at high luminosity and high redshift taken from literature.

Ultra-hot Jupiters are emerging as a new class of exoplanets. Studying their chemical compositions and temperature structures will improve the understanding of their mass loss rate as well as their formation and evolution. We present the detection of ionized calcium in the two hottest giant exoplanets - KELT-9b and WASP-33b. By utilizing transit datasets from CARMENES and HARPS-N observations, we achieved high confidence level detections of Ca II using the cross-correlation method. We further obtain the transmission spectra around the individual lines of the Ca II H&K doublet and the near-infrared triplet, and measure their line profiles. The Ca II H&K lines have an average line depth of 2.02 $\pm$ 0.17 % (effective radius of 1.56 Rp) for WASP-33b and an average line depth of 0.78 $\pm$ 0.04 % (effective radius of 1.47 Rp) for KELT-9b, which indicates that the absorptions are from very high upper atmosphere layers close to the planetary Roche lobes. The observed Ca II lines are significantly deeper than the predicted values from the hydrostatic models. Such a discrepancy is probably a result of hydrodynamic outflow that transports a significant amount of Ca II into the upper atmosphere. The prominent Ca II detection with the lack of significant Ca I detection implies that calcium is mostly ionized in the upper atmospheres of the two planets.

The He I infrared (IR) triplet at 10830 AA is an important activity indicator for the Sun and in solar-type stars, however, it has rarely been studied in relation to M dwarfs to date. In this study, we use the time-averaged spectra of 319 single stars with spectral types ranging from M0.0 V to M9.0 V obtained with the CARMENES high resolution optical and near-infrared spectrograph at Calar Alto to study the properties of the He I IR triplet lines. In quiescence, we find the triplet in absorption with a decrease of the measured pseudo equivalent width (pEW) towards later sub-types. For stars later than M5.0 V, the He I triplet becomes undetectable in our study. This dependence on effective temperature may be related to a change in chromospheric conditions along the M dwarf sequence. When an emission in the triplet is observed, we attribute it to flaring. The absence of emission during quiescence is consistent with line formation by photo-ionisation and recombination, while flare emission may be caused by collisions within dense material. The He I triplet tends to increase in depth according to increasing activity levels, ultimately becoming filled in; however, we do not find a correlation between the pEW(He IR) and X-ray properties. This behaviour may be attributed to the absence of very inactive stars ($L_{\rm X}/L_{\rm bol}$ < -5.5) in our sample or to the complex behaviour with regard to increasing depth and filling in.

Dynamical black hole mass measurements in some gas-rich galaxy mergers indicate that they are overmassive relative to their host galaxy properties. Overmassive black holes in these systems present a conflict with the standard progression of galaxy merger - quasar evolution; an alternative explanation is that a nuclear concentration of molecular gas driven inward by the merger is affecting these dynamical black hole mass estimates. We test for the presence of such gas near the two black holes in NGC 6240 using long-baseline ALMA Band 6 observations (beam size 0"06 $\times$ 0"03 or 30 pc$\times$15 pc). We find (4.2-9.8) $\times10^{7}$ M$_{\odot}$ and (1.2-7.7) $\times10^{8}$ M$_{\odot}$ of molecular gas within the resolution limit of the original black hole mass measurements for the north and south black holes, respectively. In the south nucleus, this measurement implies that 6-89% of the original black hole mass measurement actually comes from molecular gas, resolving the tension in the original black hole scaling relations. For the north, only 5% to 11% is coming from molecular gas, suggesting the north black hole is actually overmassive. Our analysis provides the first measurement of significant molecular gas masses contaminating dynamical black hole mass measurements. These high central molecular gas densities further present a challenge to theoretical black hole accretion prescriptions, which often assume accretion proceeds rapidly through the central 10 pc.

We present X-ray spectra spanning 18 years of evolution for SN\,1996cr, one of the five nearest SNe detected in the modern era. {\it Chandra} HETG exposures in 2000, 2004, and 2009 allow us to resolve spectrally the velocity profiles of Ne, Mg, Si, S, and Fe emission lines and monitor their evolution as tracers of the ejecta-circumstellar medium (CSM) interaction. To explain the diversity of X-ray line profiles, we explore several possible geometrical models. Based on the highest signal-to-noise 2009 epoch, we find that a polar geometry with two distinct opening angle configurations and internal obscuration can successfully reproduce all of the observed line profiles. The best fit model consists of two plasma components: (1) a mildly absorbed (2$\times$10$^{21}$\,cm$^{-2}$), cooler ($\approx$2\,keV) with high Ne, Mg, Si, and S abundances associated with a wide polar interaction region (half-opening angle $\approx$58$^{\circ}$); (2) a moderately absorbed (2$\times$10$^{22}$\,cm$^{-2}$), hotter (>~20\,keV) plasma with high Fe abundances and strong internal obscuration associated with a narrow polar interaction region (half-opening angle $\approx$20$^{\circ}$). We extend this model to seven further epochs with lower signal-to-noise ratio and/or lower spectral-resolution between 2000-2018, yielding several interesting trends in absorption, flux, geometry and expansion velocity. We argue that the hotter and colder components are associated with reflected and forward shocks, respectively, at least at later epochs. We discuss the physical implications of our results and plausible explosion scenarios to understand the X-ray data of SN\,1996cr.

Statistical analyses from exoplanet surveys around low-mass stars indicate that super-Earth and Neptune-mass planets are more frequent than gas giants around such stars, in agreement with core accretion theory of planet formation. Using precise radial velocities derived from visual and near-infrared spectra, we report the discovery of a giant planet with a minimum mass of 0.46 Jupiter masses in an eccentric 204-day orbit around the very low-mass star GJ 3512. Dynamical models show that the high eccentricity of the orbit is most likely explained from planet-planet interactions. The reported planetary system challenges current formation theories and puts stringent constraints on the accretion and migration rates of planet formation and evolution models, indicating that disc instability may be more efficient in forming planets than previously thought.

Context. Teegarden's Star is the brightest and one of the nearest ultra-cool dwarfs in the solar neighbourhood. For its late spectral type (M7.0V), the star shows relatively little activity and is a prime target for near-infrared radial velocity surveys such as CARMENES. Aims. As part of the CARMENES search for exoplanets around M dwarfs, we obtained more than 200 radial-velocity measurements of Teegarden's Star and analysed them for planetary signals. Methods. We find periodic variability in the radial velocities of Teegarden's Star. We also studied photometric measurements to rule out stellar brightness variations mimicking planetary signals. Results. We find evidence for two planet candidates, each with $1.1M_\oplus$ minimum mass, orbiting at periods of 4.91 and 11.4 d, respectively. No evidence for planetary transits could be found in archival and follow-up photometry. Small photometric variability is suggestive of slow rotation and old age. Conclusions. The two planets are among the lowest-mass planets discovered so far, and they are the first Earth-mass planets around an ultra-cool dwarf for which the masses have been determined using radial velocities.

X-ray emission from QSOs has been used to assess SMBH accretion properties up to $z$~6. However, at $z>6$ only ~15 QSOs are covered by sensitive X-ray observations, preventing a statistically significant investigation of the X-ray properties of QSOs in the first Gyr of the Universe. We present new Chandra observations of 10 $z>6$ QSOs, selected to have virial black-hole mass estimates from Mg II line spectroscopy. Adding archival X-ray data for an additional 15 $z>6$ QSOs, we investigate the X-ray properties of the QSO population in the first Gyr of the Universe, focusing in particular on the $L_{UV}-L_{X}$ relation, which is traced by the $\alpha_{ox}$ parameter, and the shape of their X-ray spectra. We performed photometric analyses to derive estimates of the X-ray luminosities, and thus the $\alpha_{ox}$ values and bolometric corrections ($K_{bol}=L_{bol}/L_{X}$). We compared the resulting $\alpha_{ox}$ and $K_{bol}$ distributions with the results found for QSO samples at lower redshift. Finally, we performed a basic X-ray spectral analysis of the brightest $z>6$ QSOs to derive their individual photon indices, and joint spectral analysis of the whole sample to estimate the average photon index. We confirm a lack of significant evolution of $\alpha_{ox}$ with redshift, extending the results from previous works up to $z>6$, and the trend of an increasing bolometric correction with increasing luminosity found for QSOs at lower redshifts. The average power-law photon index of our sample ($\Gamma=2.20_{-0.34}^{+0.39}$ and $\Gamma=2.13_{-0.13}^{+0.13}$ for sources with $<30$ and $>30$ net counts, respectively) is slightly steeper than, but still consistent with, typical QSOs at $z=1-6$. All these results point toward a lack of substantial evolution of the inner accretion-disk/hot-corona structure in QSOs from low redshift to $z>6$. Our data hint at generally high Eddington ratios at $z>6$.

Aims: We aim at detecting H$_2$O in the atmosphere of the hot Jupiter HD 209458 b and perform a multi-band study in the near infrared with CARMENES. Methods: The H$_2$O absorption lines from the planet's atmosphere are Doppler-shifted due to the large change in its radial velocity during transit. This shift is of the order of tens of km s$^{-1}$, whilst the Earth's telluric and the stellar lines can be considered quasi-static. We took advantage of this to remove the telluric and stellar lines using SYSREM, a principal component analysis algorithm. The residual spectra contain the signal from thousands of planetary molecular lines well below the noise level. We retrieve this information by cross-correlating the spectra with models of the atmospheric absorption. Results: We find evidence of H$_2$O in HD 209458 b with a signal-to-noise ratio (S/N) of 6.4. The signal is blueshifted by --5.2 $^{+2.6}_{-1.3}$ km s$^{-1}$, which, despite the error bars, is a firm indication of day-to-night winds at the terminator of this hot Jupiter. Additionally, we performed a multi-band study for the detection of H$_2$O individually from the three NIR bands covered by CARMENES. We detect H$_2$O from its 1.0 $\mu$m band with a S/N of 5.8, and also find hints from the 1.15 $\mu$m band, with a low S/N of 2.8. No clear planetary signal is found from the 1.4 $\mu$m band. Conclusions: Our significant signal from the 1.0 $\mu$m band in HD 209458 b represents the first detection of H$_2$O from this band, the bluest one to date. The unfavorable observational conditions might be the reason for the inconclusive detection from the stronger 1.15 and 1.4 $\mu$m bands. H$_2$O is detected from the 1.0 $\mu$m band in HD 209458 b, but hardly in HD 189733 b, which supports a stronger aerosol extinction in the latter.

Context: Recently, the He I triplet at 10830 \r{A} has been rediscovered as an excellent probe of the extended and possibly evaporating atmospheres of close-in transiting planets. This has already resulted in detections of this triplet in the atmospheres of a handful of planets, both from space and from the ground. However, while a strong signal is expected for the hot Jupiter HD 209458 b, only upper limits have been obtained so far. Aims: Our goal is to measure the helium excess absorption from HD 209458 b and assess the extended atmosphere of the planet and possible evaporation. Methods: We obtained new high-resolution spectral transit time-series of HD 209458 b using CARMENES at the 3.5 m Calar Alto telescope, targeting the He I triplet at 10830 \r{A} at a spectral resolving power of 80 400. The observed spectra were corrected for stellar absorption lines using out of transit data, for telluric absorption using the molecfit software, and for the sky emission lines using simultaneous sky measurements through a second fibre. Results: We detect He I absorption at a level of 0.91 $\pm$ 0.10 % (9 $\sigma$) at mid-transit. The absorption follows the radial velocity change of the planet during transit, unambiguously identifying the planet as the source of the absorption. The core of the absorption exhibits a net blueshift of 1.8 $\pm$ 1.3 km s$^{-1}$. Possible low-level excess absorption is seen further blueward from the main absorption near the centre of the transit, which could be caused by an extended tail. However, this needs to be confirmed. Conclusions: Our results further support a close relationship between the strength of planetary absorption in the helium triplet lines and the level of ionising, stellar X-ray and extreme-UV irradiation.

We report the detection of a transiting Earth-size planet around GJ 357, a nearby M2.5V star, using data from the Transiting Exoplanet Survey Satellite (TESS). GJ 357 b (TOI-562.01) is a transiting, hot, Earth-sized planet (Teq=525+-11 K) with a radius of Rb=1.217+-0.084 Re and an orbital period of Pb=3.93 d. Precise stellar radial velocities from CARMENES and PFS, as well as archival data from HIRES, UVES, and HARPS also display a 3.93-day periodicity, confirming the planetary nature and leading to a planetary mass of Mb=1.84+-0.31 Me. In addition to the radial velocity signal for GJ 357 b, more periodicities are present in the data indicating the presence of two further planets in the system: GJ 357 c, with a minimum mass of Mc=3.40+-0.46 Me in a 9.12 d orbit, and GJ 357 d, with a minimum mass of Md=6.1+-1.0 Me in a 55.7 d orbit inside the habitable zone. The host is relatively inactive and exhibits a photometric rotation period of Prot=78+-2 d. GJ 357 b is to date the second closest transiting planet to the Sun, making it a prime target for further investigations such as transmission spectroscopy. Therefore, GJ 357 b represents one of the best terrestrial planets suitable for atmospheric characterization with the upcoming JWST and ground-based ELTs.

We present a unification model for a clumpy obscurer in active galactic nuclei (AGN) and investigate the properties of the resulting X-ray spectrum. Our model is constructed to reproduce the column density distribution of the AGN population and cloud eclipse events in terms of their angular sizes and frequency. We developed and release a generalised Monte Carlo X-ray radiative transfer code, XARS, to compute X-ray spectra of obscurer models. The geometry results in strong Compton scattering, causing soft photons to escape also along Compton-thick sight lines. This makes our model spectra very similar to the Brightman & Nandra TORUS model. However, only if we introduce an additional Compton-thick reflector near the corona, we achieve good fits to NuSTAR spectra. This additional component in our model can be interpreted as part of the dust-free broad-line region, an inner wall or rim, or a warped disk. It cannot be attributed to a simple disk because the reflector must simultaneously block the line of sight to the corona and reflect its radiation. We release our model as an Xspec table model and present corresponding CLUMPY infrared spectra, paving the way for self-consistent multi-wavelength analyses.

We present the results obtained with VLT/MUSE on the faint-end of the Lyman-alpha luminosity function (LF) based on deep observations of four lensing clusters. The precise aim of the present study is to further constrain the abundance of Lyman-alpha emitters (LAEs) by taking advantage of the magnification provided by lensing clusters. We blindly selected a sample of 156 LAEs, with redshifts between $2.9 \le z \le 6.7$ and magnification-corrected luminosities in the range $ 39 \lesssim \log L_{Ly_{\alpha}}$ [erg s$^{-1}$] $\lesssim 43$. The price to pay to benefit from magnification is a reduction of the effective volume of the survey, together with a more complex analysis procedure. To properly take into account the individual differences in detection conditions (including lensing configurations, spatial and spectral morphologies) when computing the LF, a new method based on the 1/Vmax approach was implemented. The LAE LF has been obtained in four different redshift bins with constraints down to $\log L_{Ly_{\alpha}} = 40.5$. From our data only, no significant evolution of LF mean slope can be found. When performing a Schechter analysis including data from the literature to complete the present sample a steep faint-end slope was measured varying from $\alpha = -1.69^{+0.08}_{-0.08}$ to $\alpha = -1.87^{+0.12}_{-0.12}$ between the lowest and the highest redshift bins. The contribution of the LAE population to the star formation rate density at $z \sim 6$ is $\lesssim 50$% depending on the luminosity limit considered, which is of the same order as the Lyman-break galaxy (LBG) contribution. The evolution of the LAE contribution with redshift depends on the assumed escape fraction of Lyman-alpha photons, and appears to slightly increase with increasing redshift when this fraction is conservatively set to one. (abridged)

We analyze the interstellar medium properties of a sample of sixteen bright CO line emitting galaxies identified in the ALMA Spectroscopic Survey in the Hubble Ultra Deep Field (ASPECS) Large Program. This CO$-$selected galaxy sample is complemented by a couple of additional CO line emitters in the UDF that are identified based on their MUSE optical spectroscopic redshifts. The ASPECS CO$-$selected galaxies cover a larger range of star-formation rates and stellar masses compared to literature CO emitting galaxies at $z>1$ for which scaling relations have been established previously. Most of ASPECS CO-selected galaxies follow these established relations in terms of gas depletion timescales and gas fractions as a function of redshift, as well as the star-formation rate-stellar mass relation (`galaxy main sequence'). However, we find that $\sim30\%$ of the galaxies (5 out of 16) are offset from the galaxy main sequence at their respective redshift, with $\sim12\%$ (2 out of 16) falling below this relationship. Some CO-rich galaxies exhibit low star-formation rates, and yet show substantial molecular gas reservoirs, yielding long gas depletion timescales. Capitalizing on the well-defined cosmic volume probed by our observations, we measure the contribution of galaxies above, below, and on the galaxy main sequence to the total cosmic molecular gas density at different lookback times. We conclude that main sequence galaxies are the largest contributor to the molecular gas density at any redshift probed by our observations (z$\sim$1$-$3). The respective contribution by starburst galaxies above the main sequence decreases from z$\sim$2.5 to z$\sim$1, whereas we find tentative evidence for an increased contribution to the cosmic molecular gas density from the passive galaxies below the main sequence.

We present precise photospheric parameters of 282 M dwarfs determined from fitting the most recent version of PHOENIX models to high-resolution CARMENES spectra in the visible (0.52 - 0.96 $\mu$m) and near-infrared wavelength range (0.96 - 1.71 $\mu$m). With its aim to search for habitable planets around M dwarfs, several planets of different masses have been detected. The characterization of the target sample is important for the ability to derive and constrain the physical properties of any planetary systems that are detected. As a continuation of previous work in this context, we derived the fundamental stellar parameters effective temperature, surface gravity, and metallicity of the CARMENES M-dwarf targets from PHOENIX model fits using a $\chi^2$ method. We calculated updated PHOENIX stellar atmosphere models that include a new equation of state to especially account for spectral features of low-temperature stellar atmospheres as well as new atomic and molecular line lists. We show the importance of selecting magnetically insensitive lines for fitting to avoid effects of stellar activity in the line profiles. For the first time, we directly compare stellar parameters derived from multiwavelength range spectra, simultaneously observed for the same star. In comparison with literature values we show that fundamental parameters derived from visible spectra and visible and near-infrared spectra combined are in better agreement than those derived from the same spectra in the near-infrared alone.

While theoretical arguments predict that most of the early growth of supermassive black holes (SMBHs) happened during heavily obscured phases of accretion, current methods used for selecting $z>6$ quasars (QSOs) are strongly biased against obscured QSOs, thus considerably limiting our understanding of accreting SMBHs during the first Gyr of the Universe from an observational point of view. We report the $Chandra$ discovery of the first heavily obscured QSO candidate in the early universe, hosted by a close ($\approx5$ kpc) galaxy pair at $z=6.515$. One of the members is an optically classified type 1 QSO, PSO167-13. The companion galaxy was first detected as a [C II] emitter by ALMA. An X-ray source is significantly ($P=0.9996$) detected by $Chandra$ in the 2-5 keV band, with $<1.14$ net counts in the 0.5-2 keV band, although the current positional uncertainty does not allow a conclusive association with either PSO167-13 or its companion galaxy. From X-ray photometry and hardness-ratio arguments, we estimated an obscuring column density of $N_H>2\times10^{24}\,\mathrm{cm^{-2}}$ and $N_H>6\times10^{23}\,\mathrm{cm^{-2}}$ at $68\%$ and $90\%$ confidence levels, respectively. Thus, regardless of which of the two galaxies is associated with the X-ray emission, this source is the first heavily obscured QSO candidate at $z>6$.

Recent works have discovered two fast ($\approx 10$ ks) extragalactic X-ray transients in the Chandra Deep Field-South (CDF-S XT1 and XT2). These findings suggest that a large population of similar extragalactic transients might exist in archival X-ray observations. We develop a method that can effectively detect such transients in a single Chandra exposure, and systematically apply it to Chandra surveys of CDF-S, CDF-N, DEEP2, UDS, COSMOS, and E-CDF-S, totaling 19~Ms of exposure. We find 13 transient candidates, including CDF-S XT1 and XT2. With the aid of available excellent multiwavelength observations, we identify the physical nature of all these candidates. Aside from CDF-S XT1 and XT2, the other 11 sources are all stellar objects, and all of them have $z$-band magnitudes brighter than 20. We estimate an event rate of $59^{+77}_{-38}\ \rm{evt\ yr^{-1} deg^{-2}}$ for CDF-S XT-like transients with 0.5-7 keV peak fluxes $\log F_{\rm peak} \gtrsim -12.6$ (erg cm$^{-2}$ s$^{-1}$). This event rate translates to $\approx 15^{+20}_{-10}$ transients existing among Chandra archival observations at Galactic latitudes $|b|>20^{\circ}$, which can be probed in future work. Future missions such as Athena and the Einstein Probe with large grasps (effective area $\times$ field of view) are needed to discover a large sample ($\sim$ thousands) of fast extragalactic X-ray transients.

The scatter in the relationship between the strength of [CII] 158$\mu$m emission and the star formation rate at high-redshift has been the source of much recent interest. Although the relationship is well-established locally, several intensely star-forming galaxies have been found whose [CII] 158$\mu$m emission is either weak, absent or spatially offset from the young stars. Here we present new ALMA data for the two most distant, gravitationally-lensed and spectroscopically-confirmed galaxies, A2744\_YD4 at $z=$8.38 and MACS1149\_JD1 at $z=$9.11, both of which reveal intense [OIII] 88$\mu$m emission. In both cases we provide stringent upper limits on the presence of [CII] 158$\mu$m with respect to [OIII] 88$\mu$m. We review possible explanations for this apparent redshift-dependent [CII] deficit in the context of our recent hydrodynamical simulations. Our results highlight the importance of using several emission line diagnostics with ALMA to investigate the nature of the interstellar medium in early galaxies.

The analysis of the variability of active galactic nuclei (AGNs) at different wavelengths and the study of possible correlations among different spectral windows are nowadays a major field of inquiry. Optical variability has been largely used to identify AGNs in multivisit surveys. The strength of a selection based on optical variability lies in the chance to analyze data from surveys of large sky areas by ground-based telescopes. However the effectiveness of optical variability selection, with respect to other multiwavelength techniques, has been poorly studied down to the depth expected from next generation surveys. Here we present the results of our r-band analysis of a sample of 299 optically variable AGN candidates in the VST survey of the COSMOS field, counting 54 visits spread over three observing seasons spanning > 3 yr. This dataset is > 3 times larger in size than the one presented in our previous analysis (De Cicco et al. 2015), and the observing baseline is ~8 times longer. We push towards deeper magnitudes (r(AB) ~23.5 mag) compared to past studies; we make wide use of ancillary multiwavelength catalogs in order to confirm the nature of our AGN candidates, and constrain the accuracy of the method based on spectroscopic and photometric diagnostics. We also perform tests aimed at assessing the relevance of dense sampling in view of future wide-field surveys. We demonstrate that the method allows the selection of high-purity (> 86%) samples. We take advantage of the longer observing baseline to achieve great improvement in the completeness of our sample with respect to X-ray and spectroscopically confirmed samples of AGNs (59%, vs. ~15% in our previous work), as well as in the completeness of unobscured and obscured AGNs. The effectiveness of the method confirms the importance to develop future, more refined techniques for the automated analysis of larger datasets.

Although M dwarfs are known for high levels of stellar activity, they are ideal targets for the search of low-mass exoplanets with the radial velocity (RV) method. We report the discovery of a planetary-mass companion around LSPM J2116+0234 (M3.0 V) and confirm the existence of a planet orbiting GJ 686 (BD+18 3421; M1.0 V). The discovery of the planet around LSPM J2116+0234 is based on CARMENES RV observations in the visual and near-infrared channels. We confirm the planet orbiting around GJ 686 by analyzing the RV data spanning over two decades of observations from CARMENES VIS, HARPS-N, HARPS, and HIRES. We find planetary signals at 14.44 and 15.53 d in the RV data for LSPM J2116+0234 and GJ 686, respectively. Additionally, the RV, photometric time series, and various spectroscopic indicators show hints of variations of 42 d for LSPM J2116+0234 and 37 d for GJ 686, which we attribute to the stellar rotation periods. The orbital parameters of the planets are modeled with Keplerian fits together with correlated noise from the stellar activity. A mini-Neptune with a minimum mass of 11.8 Me orbits LSPM J2116+0234 producing an RV semi-amplitude of 6.19 m/s, while a super-Earth of mass 6.6 Me orbits GJ 686 and produces an RV semi-amplitude of 3.0 m/s. Both LSPM J2116+0234 and GJ 686 have planetary companions populating the regime of exoplanets with masses lower than 15 Me and orbital periods <20 d.

M dwarfs are known to generate the strongest magnetic fields among main-sequence stars with convective envelopes, but the link between the magnetic fields and underlying dynamo mechanisms, rotation, and activity still lacks a consistent picture. In this work we measure magnetic fields from the high-resolution near-infrared spectra taken with the CARMENES radial-velocity planet survey in a sample of 29 active M dwarfs and compare our results against stellar parameters. We use the state-of-the-art radiative transfer code to measure total magnetic flux densities from the Zeeman broadening of spectral lines and filling factors. We detect strong kG magnetic fields in all our targets. In 16 stars the magnetic fields were measured for the first time. Our measurements are consistent with the magnetic field saturation in stars with rotation periods P<4d. The analysis of the magnetic filling factors reveal two different patterns of either very smooth distribution or a more patchy one, which can be connected to the dynamo state of the stars and/or stellar mass. Our measurements extend the list of M dwarfs with strong surface magnetic fields. They also allow us to better constrain the interplay between the magnetic energy, stellar rotation, and underlying dynamo action. The high spectral resolution and observations at near-infrared wavelengths are the beneficial capabilities of the CARMENES instrument that allow us to address important questions about the stellar magnetism.

Neutron star-neutron star mergers are known to be associated with short gamma-ray bursts. If the neutron star equation of state is sufficiently stiff, at least some of such mergers will leave behind a supramassive or even a stable neutron star that spins rapidly with a strong magnetic field (i.e., a magnetar). Such a magnetar signature may have been observed as the X-ray plateau following a good fraction (up to 50%) of short gamma-ray bursts, and it has been expected that one may observe short gamma-ray burst-less X-ray transients powered by double neutron star mergers. A fast X-ray transient (CDF-S XT1) was recently found to be associated with a faint host galaxy whose redshift is unknown. Its X-ray and host-galaxy properties allow several possibleexplanations including a short gamma-ray burst seen off axis, a low-luminosity gamma-ray burst at high redshift, or a tidal disruption event involving an intermediate mass black hole and a white dwarf. Here we report a second X-ray transient, CDF-S XT2, that is associated with a galaxy at redshift z = 0.738. The light curve is fully consistent with being powered by a millisecond magnetar. More intriguingly, CDF-S XT2 lies in the outskirts of its star-forming host galaxy with a moderate offset from the galaxy center, as short bursts often do. The estimated event rate density of similar X-ray transients, when corrected to the local value, is consistent with the double neutron star merger rate density inferred from the detection of GW170817.

We used data from the QUEST-La Silla Active Galactic Nuclei (AGN) variability survey to construct light curves for 208,583 sources over $\sim 70$ deg$^2$, with a a limiting magnitude $r \sim 21$. Each light curve has at least 40 epochs and a length of $\geq 200$ days. We implemented a Random Forest algorithm to classify our objects as either AGN or non-AGN according to their variability features and optical colors, excluding morphology cuts. We tested three classifiers, one that only includes variability features (RF1), one that includes variability features and also $r-i$ and $i-z$ colors (RF2), and one that includes variability features and also $g-r$, $r-i$, and $i-z$ colors (RF3). We obtained a sample of high probability candidates (hp-AGN) for each classifier, with 5,941 candidates for RF1, 5,252 candidates for RF2, and 4,482 candidates for RF3. We divided each sample according to their $g-r$ colors, defining blue ($g-r\leq 0.6$) and red sub-samples ($g-r>0.6$). We find that most of the candidates known from the literature belong to the blue sub-samples, which is not necessarily surprising given that, unlike for many literature studies, we do not cut our sample to point-like objects. This means that we can select AGN that have a significant contribution from redshifted starlight in their host galaxies. In order to test the efficiency of our technique we performed spectroscopic follow-up, confirming the AGN nature of 44 among 54 observed sources (81.5\% of efficiency). From the campaign we concluded that RF2 provides the purest sample of AGN candidates.

We present an analysis of a sample of 69 local obscured Swift/Burst Alert Telescope active galactic nuclei (AGNs) with X-ray spectra from NuSTAR and infrared (IR) spectral energy distributions from Herschel and WISE. We combine this X-ray and IR phenomenological modeling and find a significant correlation between reflected hard X-ray emission and IR AGN emission, with suggestive indications that this correlation may be stronger than the one between intrinsic hard X-ray and IR emissions. This relation between the IR and reflected X-ray emission suggests that both are the result of the processing of intrinsic emission from the corona and accretion disk by the same structure. We explore the resulting implications on the underlying distribution of covering fraction for all AGNs, by generating mock observables for the reflection parameter and IR luminosity ratio using empirical relations found for the covering fraction with each quantity. We find that the observed distributions of the reflection parameter and IR-to-X-ray ratio are reproduced with broad distributions centered around covering fractions of at least ~40%-50%, whereas narrower distributions match our observations only when centered around covering fractions of ~70%-80%. Our results are consistent with both independent estimates of the covering fractions of individual objects and the typical covering fraction obtained on the basis of obscured fractions for samples of AGNs. These results suggest that the level of reprocessing in AGNs, including X-ray reflection, is related in a relatively straightforward way to the geometry of the obscuring material.

We determine the radii and masses of 293 nearby, bright M dwarfs of the CARMENES survey. This is the first time that such a large and homogeneous high-resolution (R>80 000) spectroscopic survey has been used to derive these fundamental stellar parameters. We derived the radii using Stefan-Boltzmann's law. We obtained the required effective temperatures $T_{\rm eff}$ from a spectral analysis and we obtained the required luminosities L from integrated broadband photometry together with the Gaia DR2 parallaxes. The mass was then determined using a mass-radius relation that we derived from eclipsing binaries known in the literature. We compared this method with three other methods: (1) We calculated the mass from the radius and the surface gravity log g, which was obtained from the same spectral analysis as $T_{\rm eff}$. (2) We used a widely used infrared mass-magnitude relation. (3) We used a Bayesian approach to infer stellar parameters from the comparison of the absolute magnitudes and colors of our targets with evolutionary models. Between spectral types M0V and M7V our radii cover the range $0.1\,R_{\normalsize\odot}<R<0.6\,R_{\normalsize\odot}$ with an error of 2-3% and our masses cover $0.09\,{\mathcal M}_{\normalsize\odot}<{\mathcal M}<0.6\,{\mathcal M}_{\normalsize\odot}$ with an error of 3-5%. We find good agreement between the masses determined with these different methods for most of our targets. Only the masses of very young objects show discrepancies. This can be well explained with the assumptions that we used for our methods.

We present an updated strong-lensing analysis of the massive cluster Abell 370 (A370), continuing the work first presented in Lagattuta et al. (2017). In this new analysis, we take advantage of the deeper imaging data from the Hubble Space Telescope (HST) Frontier Fields program, as well as a large spectroscopic mosaic obtained with the Multi-Unit Spectroscopic Explorer (MUSE). Thanks to the extended coverage of this mosaic, we probe the full 3D distribution of galaxies in the field, giving us a unique picture of the extended structure of the cluster and its surroundings. Our final catalog contains 584 redshifts, representing the largest spectroscopic catalog of A370 to date. Constructing the model, we measure a total mass distribution that is quantitatively similar to our previous work -- though to ensure a low rms error in the model fit, we invoke a significantly large external shear term. Using the redshift catalog, we search for other bound groups of galaxies, which may give rise to a more physical interpretation of this shear. We identify three structures in narrow redshift ranges along the line of sight, highlighting possible infalling substructures into the main cluster halo. We also discover additional substructure candidates in low-resolution imaging at larger projected radii. More spectroscopic coverage of these regions (pushing close to the A370 virial radius) and more extended, high-resolution imaging will be required to investigate this possibility, further advancing the analysis of these interesting developments.

Small planets around low-mass stars often show orbital periods in a range that corresponds to the temperate zones of their host stars which are therefore of prime interest for planet searches. Surface phenomena such as spots and faculae create periodic signals in radial velocities and in observational activity tracers in the same range, so they can mimic or hide true planetary signals. We aim to detect Doppler signals corresponding to planetary companions, determine their most probable orbital configurations, and understand the stellar activity and its impact on different datasets. We analyze 22 years of data of the M1.5V-type star Gl49 (BD+61 195) including HARPS-N and CARMENES spectrographs, complemented by APT2 and SNO photometry. Activity indices are calculated from the observed spectra, and all datasets are analyzed with periodograms and noise models. We investigate how the variation of stellar activity imprints on our datasets. We further test the origin of the signals and investigate phase shifts between the different sets. To search for the best-fit model we maximize the likelihood function in a Markov Chain Monte Carlo approach. As a result of this study, we are able to detect the super-Earth Gl49b with a minimum mass of 5.6 Ms. It orbits its host star with a period of 13.85d at a semi-major axis of 0.090 au and we calculate an equilibrium temperature of 350 K and a transit probability of 2.0%. The contribution from the spot-dominated host star to the different datasets is complex, and includes signals from the stellar rotation at 18.86d, evolutionary time-scales of activity phenomena at 40-80d, and a long-term variation of at least four years.

We discuss the nature and physical properties of gas-mass selected galaxies in the ALMA spectroscopic survey (ASPECS) of the Hubble Ultra Deep Field (HUDF). We capitalize on the deep optical integral-field spectroscopy from the MUSE HUDF Survey and multi-wavelength data to uniquely associate all 16 line-emitters, detected in the ALMA data without preselection, with rotational transitions of carbon monoxide (CO). We identify ten as CO(2-1) at $1 < z < 2$, five as CO(3-2) at $2 < z < 3$ and one as CO(4-3) at $z = 3.6$. Using the MUSE data as a prior, we identify two additional CO(2-1)-emitters, increasing the total sample size to 18. We infer metallicities consistent with (super-)solar for the CO-detected galaxies at $z \le 1.5$, motivating our choice of a Galactic conversion factor between CO luminosity and molecular gas mass for these galaxies. Using deep Chandra imaging of the HUDF, we determine an X-ray AGN fraction of 20% and 60% among the CO-emitters at $z \sim 1.4$ and $z \sim 2.6$, respectively. Being a CO-flux limited survey, ASPECS-LP detects molecular gas in galaxies on, above and below the main sequence (MS) at $z \sim 1.4$. For stellar masses $\ge 10^{10} (10^{10.5})$ M$_{\odot}$, we detect about 40% (50%) of all galaxies in the HUDF at $1 < z < 2$ ($2 < z < 3$). The combination of ALMA and MUSE integral-field spectroscopy thus enables an unprecedented view on MS galaxies during the peak of galaxy formation.

We use the results from the ALMA large program ASPECS, the spectroscopic survey in the Hubble Ultra Deep Field (HUDF), to constrain CO luminosity functions of galaxies and the resulting redshift evolution of $\rho$(H$_2$). The broad frequency range covered enables us to identify CO emission lines of different rotational transitions in the HUDF at $z>1$. We find strong evidence that the CO luminosity function evolves with redshift, with the knee of the CO luminosity function decreasing in luminosity by an order of magnitude from $\sim$2 to the local universe. Based on Schechter fits, we estimate that our observations recover the majority (up to $\sim$90%, depending on the assumptions on the faint end) of the total cosmic CO luminosity at $z$=1.0-3.1. After correcting for CO excitation, and adopting a Galactic CO-to-H$_2$ conversion factor, we constrain the evolution of the cosmic molecular gas density $\rho$(H$_2$): this cosmic gas density peaks at $z\sim1.5$ and drops by factor of $6.5_{-1.4}^{+1.8}$ to the value measured locally. The observed evolution in $\rho$(H$_2$) therefore closely matches the evolution of the cosmic star formation rate density $\rho_{\rm SFR}$. We verify the robustness of our result with respect to assumptions on source inclusion and/or CO excitation. As the cosmic star formation history can be expressed as the product of the star formation efficiency and the cosmic density of molecular gas, the similar evolution of $\rho$(H$_2$) and $\rho_{\rm SFR}$ leaves only little room for a significant evolution of the average star formation efficiency in galaxies since $z\sim 3$ (85% of cosmic history).

The ALMA SPECtroscopic Survey in the {\it Hubble} Ultra Deep Field is an ALMA large program that obtained a frequency scan in the 3\,mm band to detect emission lines from the molecular gas in distant galaxies. We here present our search strategy for emission lines and continuum sources in the HUDF. We compare several line search algorithms used in the literature, and critically account for the line-widths of the emission line candidates when assessing significance. We identify sixteen emission lines at high fidelity in our search. Comparing these sources to multi-wavelength data we find that all sources have optical/infrared counterparts. Our search also recovers candidates that have lower significance that can be used statistically to derive, e.g. the CO luminosity function. We apply the same detection algorithm to obtain a sample of six 3 mm continuum sources. All of these are also detected in the 1.2 mm continuum with optical/near-infrared counterparts. We use the continuum sources to compute 3 mm number counts in the sub-mJy regime, and find them to be higher by an order of magnitude than expected for synchrotron-dominated sources. However, the number counts are consistent with those derived at shorter wavelengths (0.85--1.3\,mm) once extrapolating to 3\,mm with a dust emissivity index of $\beta=1.5$, dust temperature of 35\,K and an average redshift of $z=2.5$. These results represent the best constraints to date on the faint end of the 3 mm number counts.

One of the last missing pieces in the puzzle of galaxy formation and evolution through cosmic history is a detailed picture of the role of the cold gas supply in the star-formation process. Cold gas is the fuel for star formation, and thus regulates the buildup of stellar mass, both through the amount of material present through a galaxy's gas mass fraction, and through the efficiency at which it is converted to stars. Over the last decade, important progress has been made in understanding the relative importance of these two factors along with the role of feedback, and the first measurements of the volume density of cold gas out to redshift 4, (the "cold gas history of the Universe") has been obtained. To match the precision of measurements of the star formation and black-hole accretion histories over the coming decades, a two orders of magnitude improvement in molecular line survey speeds is required compared to what is possible with current facilities. Possible pathways towards such large gains include significant upgrades to current facilities like ALMA by 2030 (and beyond), and eventually the construction of a new generation of radio-to-millimeter wavelength facilities, such as the next generation Very Large Array (ngVLA) concept.

Chromospheric modeling of observed differences in stellar activity lines is imperative to fully understand the upper atmospheres of late-type stars. We present one-dimensional parametrized chromosphere models computed with the atmosphere code PHOENIX using an underlying photosphere of 3500 K. The aim of this work is to model chromospheric lines of a sample of 50 M2-3 dwarfs observed in the framework of the CARMENES, the Calar Alto high-Resolution search for M dwarfs with Exo-earths with Near-infrared and optical Echelle Spectrographs, exoplanet survey. The spectral comparison between observed data and models is performed in the chromospheric lines of Na I D2, H$\alpha$, and the bluest Ca II infrared triplet line to obtain best-fit models for each star in the sample. We find that for inactive stars a single model with a VAL C-like temperature structure is sufficient to describe simultaneously all three lines adequately. Active stars are rather modeled by a combination of an inactive and an active model, also giving the filling factors of inactive and active regions. Moreover, the fitting of linear combinations on variable stars yields relationships between filling factors and activity states, indicating that more active phases are coupled to a larger portion of active regions on the surface of the star.

The CARMENES survey is searching for Earth-like planets orbiting M dwarfs using the radial velocity method. Studying the stellar activity of the target stars is important to avoid false planet detections and to improve our understanding of the atmospheres of late-type stars. In this work we present measurements of activity indicators at visible and near-infrared wavelengths for 331 M dwarfs observed with CARMENES. Our aim is to identify the activity indicators that are most sensitive and easiest to measure, and the correlations among these indicators. We also wish to characterise their variability. Using a spectral subtraction technique, we measured pseudo-equivalent widths of the He I D3, H$\alpha$, He I $\lambda$10833 {\AA}, and Pa$\beta$ lines, the Na I D doublet, and the Ca II infrared triplet, which have a chromospheric component in active M dwarfs. In addition, we measured an index of the strength of two TiO and two VO bands, which are formed in the photosphere. We also searched for periodicities in these activity indicators for all sample stars using generalised Lomb-Scargle periodograms. We find that the most slowly rotating stars of each spectral subtype have the strongest H$\alpha$ absorption. H$\alpha$ is correlated most strongly with He I D3, whereas Na I D and the Ca II infrared triplet are also correlated with H$\alpha$. He I $\lambda$10833 {\AA} and Pa$\beta$ show no clear correlations with the other indicators. The TiO bands show an activity effect that does not appear in the VO bands. We find that the relative variations of H$\alpha$ and He I D3 are smaller for stars with higher activity levels, while this anti-correlation is weaker for Na I D and the Ca II infrared triplet, and is absent for He I $\lambda$10833 {\AA} and Pa$\beta$. Periodic variation with the rotation period most commonly appears in the TiO bands, H$\alpha$, and in the Ca II infrared triplet.

We use spectra from CARMENES, the Calar Alto high-Resolution search for M dwarfs with Exo-earths with Near-infrared and optical Echelle Spectrographs, to search for periods in chromospheric indices in 16 M0 to M2 dwarfs. We measure spectral indices in the H$\alpha$, the Ca II infrared triplet (IRT), and the Na I D lines to study which of these indices are best-suited to find rotation periods in these stars. Moreover, we test a number of different period-search algorithms, namely the string length method, the phase dispersion minimisation, the generalized Lomb-Scargle periodogram, and the Gaussian process regression with quasi-periodic kernel. We find periods in four stars using H$\alpha$ and in five stars using the Ca II IRT, two of which have not been found before. Our results show that both H$\alpha$ and the Ca II IRT lines are well suited for period searches, with the Ca II IRT index performing slightly better than H$\alpha$. Unfortunately, the Na I D lines are strongly affected by telluric airglow, and we could not find any rotation period using this index. Further, different definitions of the line indices have no major impact on the results. Comparing the different search methods, the string length method and the phase dispersion minimisation perform worst, while Gaussian process models produce the smallest numbers of false positives and non-detections.

We investigate the abundance of Super-Massive Black Hole (SMBH) seeds in primordial galaxy halos. We explore the assumption that dark matter halos outgrowing a critical halo mass M_c have some probability p of having spawned a SMBH seed. Current observations of local, intermediate-mass galaxies constrain these parameters: For $M_c=10^{11}M_\odot$, all halos must be seeded, but when adopting smaller M_c masses the seeding can be much less efficient. The constraints also put lower limits on the number density of black holes in the local and high-redshift Universe. Reproducing z~6 quasar space densities depends on their typical halo mass, which can be constrained by counting nearby Lyman Break Galaxies and Lyman Alpha Emitters. For both observables, our simulations demonstrate that single-field predictions are too diverse to make definitive statements, in agreement with mixed claims in the literature. If quasars are not limited to the most massive host halos, they may represent a tiny fraction (~10^-5) of the SMBH population. Finally, we produce a wide range of predictions for gravitational events from SMBH mergers. We define a new diagnostic diagram for LISA to measure both SMBH space density and the typical delay between halo merger and black hole merger. While previous works have explored specific scenarios, our results hold independent of the seed mechanism, seed mass, obscuration, fueling methods and duty cycle.

We present the discovery and characterization of a new transiting planet from Campaign 17 of the Kepler extended mission K2. HD 119130 b is a warm sub-Neptune on a 17-d orbit around a bright (V = 9.9 mag) solar-like G3 V star with a mass and radius of $M_\star = 1.00\pm0.03\,\mathrm{M_\odot}$ and $R_\star = 1.09\pm0.03\,\mathrm{R_\odot}$, respectively. We model simultaneously the K2 photometry and CARMENES spectroscopic data and derive a radius of $R_\mathrm{p} = 2.63_{-0.10}^{+0.12}\,\mathrm{R_\oplus}$ and mass of $M_\mathrm{p} = 24.5_{-4.4}^{+4.4}\,\mathrm{M_\oplus}$, yielding a mean density of $\rho_\mathrm{p} = 7.4_{-1.5}^{+1.6}\,\mathrm{g\,cm^{-3}}$, which makes it one of the densest sub-Neptune planets known to date. We also detect a linear trend in radial velocities of HD 119130 ($\dot{\gamma}_{\rm RV}= -0.40^{+0.07}_{-0.07}\,\mathrm{m\,s^{-1}\,d^{-1}}$) that suggests a long-period companion with a minimum mass on the order of $33\,\mathrm{M_\oplus}$. If confirmed, it would support a formation scenario of HD 119130 b by migration caused by Kozai-Lidov oscillations.

Hot gas giant exoplanets can lose part of their atmosphere due to strong stellar irradiation, affecting their physical and chemical evolution. Studies of atmospheric escape from exoplanets have mostly relied on space-based observations of the hydrogen Lyman-{\alpha} line in the far ultraviolet which is strongly affected by interstellar absorption. Using ground-based high-resolution spectroscopy we detect excess absorption in the helium triplet at 1083 nm during the transit of the Saturn-mass exoplanet WASP-69b, at a signal-to-noise ratio of 18. We measure line blue shifts of several km/s and post transit absorption, which we interpret as the escape of part of the atmosphere trailing behind the planet in comet-like form. [Additional notes by authors: Furthermore, we provide upper limits for helium signals in the atmospheres of the exoplanets HD 209458b, KELT-9b, and GJ 436b. We investigate the host stars of all planets with detected helium signals and those of the three planets we derive upper limits for. In each case we calculate the X-ray and extreme ultraviolet flux received by these planets. We find that helium is detected in the atmospheres of planets (orbiting the more active stars and) receiving the larger amount of irradiation from their host stars.]

We present three transit observations of HD 189733 b obtained with the high-resolution spectrograph CARMENES at Calar Alto. A strong absorption signal is detected in the near-infrared He I triplet at 10830 \AA{} in all three transits. During mid-transit, the mean absorption level is $0.88\pm0.04$ % measured in a $\pm$10 km s$^{-1}$ range at a net blueshift of $-3.5\pm0.4$ km s$^{-1}$ (10829.84--10830.57 \AA{}). The absorption signal exhibits radial velocities of $+6.5\pm3.1$ km s$^{-1}$ and $-12.6\pm1.0$ km s$^{-1}$ during ingress and egress, respectively; measured in the planetary rest frame. We show that stellar activity related pseudo-signals interfere with the planetary atmospheric absorption signal. They could contribute as much as 80% of the observed signal and might also affect the radial velocity signature, but pseudo-signals are very unlikely to explain the entire signal. The observed line ratio between the two unresolved and the third line of the He I triplet is $2.8\pm0.2$, which strongly deviates from the value expected for an optically thin atmospheres. When interpreted in terms of absorption in the planetary atmosphere, this favors a compact helium atmosphere with an extent of only 0.2 planetary radii and a substantial column density on the order of $4\times 10^{12}$ cm$^{-2}$. The observed radial velocities can be understood either in terms of atmospheric circulation with equatorial superrotation or as a sign of an asymmetric atmospheric component of evaporating material. We detect no clear signature of ongoing evaporation, like pre- or post-transit absorption, which could indicate material beyond the planetary Roche lobe, or radial velocities in excess of the escape velocity. These findings do not contradict planetary evaporation, but only show that the detected helium absorption in HD 189733 b does not trace the atmospheric layers that show pronounced escape signatures.

Aims: We explore the capabilities of CARMENES for characterizing hot-Jupiter atmospheres by targeting multiple water bands, in particular, those at 1.15 and 1.4 $\mu$m. Hubble Space Telescope observations suggest that this wavelength region is relevant for distinguishing between hazy/cloudy and clear atmospheres. Methods: We observed one transit of the hot Jupiter HD 189733 b with CARMENES. Telluric and stellar absorption lines were removed using Sysrem, which performs a principal component analysis including proper error propagation. The residual spectra were analysed for water absorption with cross-correlation techniques using synthetic atmospheric absorption models. Results: We report a cross-correlation peak at a signal-to-noise ratio (SNR) of 6.6, revealing the presence of water in the transmission spectrum of HD 189733 b. The absorption signal appeared slightly blueshifted at -3.9 $\pm$ 1.3 kms$^{-1}$. We measured the individual cross-correlation signals of the water bands at 1.15 and 1.4 $\mu$m, finding cross-correlation peaks at SNRs of 4.9 and 4.4, respectively. The 1.4 $\mu$m feature is consistent with that observed with the Hubble Space Telescope. Conclusions: The water bands studied in this work have been mainly observed in a handful of planets from space. The ability of also detecting them individually from the ground at higher spectral resolution can provide insightful information to constrain the properties of exoplanet atmospheres. Although the current multiband detections can not yet constrain atmospheric haze models for HD 189733 b, future observations at higher signal-to-noise ratio could provide an alternative way to achieve this aim.

This white paper specifies the footprints, cadence requirements, and total-depth requirements needed to allow the most-successful AGN studies in the four currently selected LSST Deep-Drilling Fields (DDFs): ELAIS-S1, XMM-LSS, CDF-S, and COSMOS. The information provided on cadence and total-depth requirements will also likely be applicable to enabling effective AGN science in any additional DDFs that are chosen.

We quantify the luminosity contribution of active galactic nuclei (AGN) to the 12 $\mu$m, mid-infrared (MIR; 5-38 $\mu$m), and the total IR (5-1000 $\mu$m) emission in the local AGN detected in the all-sky 70-month Swift/Burst Alert Telescope (BAT) ultra hard X-ray survey. We decompose the IR spectral energy distributions (SEDs) of 587 objects into AGN and starburst components using AGN torus and star-forming galaxy templates. This enables us to recover the AGN torus emission also for low-luminosity end, down to $\log (L_{14-150}/{\rm erg}~{\rm s}^{-1}) \simeq 41$, which typically have significant host galaxy contamination. We find that the luminosity contribution of the AGN to the 12 $\mu$m, the MIR, and the total IR band is an increasing function of the 14-150 keV luminosity. We also find that for the most extreme cases, the IR pure-AGN emission from the torus can extend up to 90 $\mu$m. The obtained total IR AGN luminosity through the IR SED decomposition enables us to estimate the fraction of the sky obscured by dust, i.e., the dust covering factor. We demonstrate that the median of the dust covering factor is always smaller than that of the X-ray obscuration fraction above the AGN bolometric luminosity of $\log (L_{\rm bol}/{\rm erg}~{\rm s}^{-1}) \simeq 42.5$. Considering that X-ray obscuration fraction is equivalent to the covering factor coming from both the dust and gas, it indicates that an additional neutral gas component, along with the dusty torus, is responsible for the absorption of X-ray emission.

We present Atacama Large Millimeter/submillimeter Array (ALMA) resolved observations of molecular gas in galaxies up to $z=0.35$ to characterise the role of global galactic dynamics on the global interstellar medium (ISM) properties. These observations consist of a sub-sample of 39 galaxies taken from the Valpara\'iso ALMA Line Emission Survey (VALES). From the CO($J=1-0)$ emission line, we quantify the kinematic parameters by modelling the velocity fields. We find that the IR luminosity increases with the rotational to dispersion velocity ratio ($V_{\rm rot}/\sigma_v$, corrected for inclination). We find a dependence between $V_{\rm rot}/\sigma_v$ and the [CII]/IR ratio, suggesting that the so-called `[CII] deficit' is related to the dynamical state of the galaxies. We find that global pressure support is needed to reconcile the dynamical mass estimates with the stellar masses in our systems with low $V_{\rm rot}/\sigma_v$ values. The star formation rate (SFR) is weakly correlated with the molecular gas fraction ($f_{\rm H_2}$) in our sample, suggesting that the release of gravitational energy from cold gas may not be the main energy source of the turbulent motions seen in the VALES galaxies. By defining a proxy of the `star formation efficiency' parameter as the SFR divided by the CO luminosity (SFE$'\equiv$ SFR/L$'_{\rm CO}$), we find a constant SFE$'$ per crossing time ($t_{\rm cross}$). We suggest that $t_{\rm cross}$ may be the controlling timescale in which the star formation occurs in dusty $z\sim0.03-0.35$ galaxies.

We report a serendipitous hard X-ray (3-24 keV) detection of the blazar LBQS 1319+0039. The X-ray spectrum is consistent with powerlaw emission with a photon index of $\Gamma = 1.78$. The observed emission and published redshift imply an observed 2-10 keV luminosity of $3.4\times10^{45}$ erg s$^{-1}$. This is consistent with extrapolations from a previous soft X-ray detection, thus there is no evidence for strong variability.

[abridged] Characterizing the number counts of faint, dusty star-forming galaxies is currently a challenge even for deep, high-resolution observations in the FIR-to-mm regime. They are predicted to account for approximately half of the total extragalactic background light at those wavelengths. Searching for dusty star-forming galaxies behind massive galaxy clusters benefits from strong lensing, enhancing their measured emission while increasing spatial resolution. Derived number counts depend, however, on mass reconstruction models that properly constrain these clusters. We estimate the 1.1 mm number counts along the line of sight of three galaxy clusters, i.e. Abell 2744, MACSJ0416.1-2403 and MACSJ1149.5+2223, which are part of the ALMA Frontier Fields Survey. We perform detailed simulations to correct these counts for lensing effects. We use several publicly available lensing models for the galaxy clusters to derive the intrinsic flux densities of our sources. We perform Monte Carlo simulations of the number counts for a detailed treatment of the uncertainties in the magnifications and adopted source redshifts. We find an overall agreement among the number counts derived for the different lens models, despite their systematic variations regarding source magnifications and effective areas. Our number counts span ~2.5 dex in demagnified flux density, from several mJy down to tens of uJy. Our number counts are consistent with recent estimates from deep ALMA observations at a 3$\sigma$ level. Below $\approx$ 0.1 mJy, however, our cumulative counts are lower by $\approx$ 1 dex, suggesting a flattening in the number counts. In our deepest ALMA mosaic, we estimate number counts for intrinsic flux densities $\approx$ 4 times fainter than the rms level. This highlights the potential of probing the sub-10 uJy population in larger samples of galaxy cluster fields with deeper ALMA observations.

The bulk of the X-ray emission in Active Galactic Nuclei (AGN) is produced very close to the accreting supermassive black hole (SMBH), in a corona of hot electrons which up scatters optical and ultraviolet photons from the accretion flow. The cutoff energy ($E_{\rm C}$) of the primary X-ray continuum emission carries important information on the physical characteristics of the X-ray emitting plasma, but little is currently known about its potential relation with the properties of accreting SMBHs. Using the largest broad-band (0.3-150 keV) X-ray spectroscopic study available to date, we investigate how the corona is related to the AGN luminosity, black hole mass and Eddington ratio ($\lambda_{\rm Edd}$). Assuming a slab corona the median values of the temperature and optical depth of the Comptonizing plasma are $kT_{\rm e}=105 \pm 18$ keV and $\tau=0.25\pm0.06$, respectively. When we properly account for the large number of $E_{\rm C}$ lower limits, we find a statistically significant dependence of the cutoff energy on the Eddington ratio. In particular, objects with $ \lambda_{\rm Edd}>0.1$ have a significantly lower median cutoff energy ($E_{\rm C}=160\pm41$ keV) than those with $\lambda_{\rm Edd}\leq 0.1$ ($E_{\rm C}=370\pm51$ keV). This is consistent with the idea that radiatively compact coronae are also cooler, because they tend to avoid the region in the temperature-compactness parameter space where runaway pair production would dominate. We show that this behaviour could also straightforwardly explain the suggested positive correlation between the photon index ($\Gamma$) and the Eddington ratio, being able to reproduce the observed slope of the $\Gamma-\lambda_{\rm Edd}$ trend.

Long gamma-ray bursts are produced during the dramatic deaths of massive stars with very short lifetimes, meaning that they explode close to the birth place of their progenitors. During a short period they become the most luminous objects observable in the Universe, being perfect beacons to study high-redshift star-forming regions. To use the afterglow of GRB 161023A at a redshift $z=2.710$ as a background source to study the environment of the explosion and the intervening systems along its line-of-sight. r the first time, we complement UV/Optical/NIR spectroscopy with millimetre spectroscopy using ALMA, which allows us to probe the molecular content of the host galaxy. The X-shooter spectrum shows a plethora of absorption features including fine-structure and metastable transitions of Fe, Ni, Si, C and O. We present photometry ranging from 43 s to over 500 days after the burst. We infer a host-galaxy metallicity of [Zn/H] $=-1.11\pm0.07$, which corrected for dust depletion results in [X/H] $=-0.94\pm0.08$. We do not detect molecular features in the ALMA data, but we derive limits on the molecular content of $log(N_{CO}/cm^{-2})<15.7$ and $log(N_{HCO+}/cm^{-2})<13.2$, which are consistent with those that we obtain from the optical spectra, $log(N_{H_2}/cm^{-2})<15.2$ and $log(N_{CO}/cm^{-2})<14.5$. Within the host galaxy we detect three velocity systems through UV/Optical/NIR absorption spectroscopy, all with levels that were excited by the GRB afterglow. We determine the distance from these systems to the GRB to be in the range between 0.7 and 1.0 kpc. The sight-line to GRB 161023A shows 9 independent intervening systems, most of them with multiple components. (Abridged)

Despite their activity, low-mass stars are of particular importance for the search of exoplanets by the means of Doppler spectroscopy, as planets with lower masses become detectable. We report on the discovery of a planetary companion around HD 180617, a bright J = 5.58 mag, low-mass M = 0.45 M_{sun} star of spectral type M2.5 V. The star, located at a distance of 5.9 pc, is the primary of the high proper motion binary system containing vB 10, a star with one of the lowest masses known in most of the twentieth century. Our analysis is based on new radial velocity (RV) measurements made at red-optical wavelengths provided by the high-precision spectrograph CARMENES, which was designed to carry out a survey for Earth-like planets around M dwarfs. The available CARMENES data are augmented by archival Doppler measurements from HIRES and HARPS. Altogether, the RVs span more than 16 years. The modeling of the RV variations, with a semi-amplitude of K = 2.85-0.25/+0.16m/s yields a Neptune-like planet with a minimum mass of 12.2-1.4/+1.0 M_{Earth} on a 105.90-0.10/+0.09d circumprimary orbit, which is partly located in the host star's habitable zone. The analysis of time series of common activity indicators does not show any dependence on the detected RV signal. The discovery of HD 180617 b not only adds information to a currently hardly filled region of the mass-period diagram of exoplanets around M dwarfs, but the investigated system becomes the third known binary consisting of M dwarfs and hosting an exoplanet in an S-type configuration. Its proximity makes it an attractive candidate for future studies.

It is well known that environment affects galaxy evolution, which is broadly related to supermassive black hole (SMBH) growth. We investigate whether SMBH evolution also depends on host-galaxy local (sub-Mpc) and global ($\approx 1-10$ Mpc) environment. We construct the surface-density field (local environment) and cosmic web (global environment) in the COSMOS field at $z=0.3-3.0$. The environments in COSMOS range from the field to clusters ($M_\mathrm{halo} \lesssim 10^{14}\ M_\odot$), covering the environments where ${\approx 99\%}$ of galaxies in the Universe reside. We measure sample-averaged SMBH accretion rate ($\overline{\mathrm{BHAR}}$) from X-ray observations, and study its dependence on overdensity and cosmic-web environment at different redshifts while controlling for galaxy stellar mass ($M_\star$). Our results show that $\overline{\mathrm{BHAR}}$ does not significantly depend on overdensity or cosmic-web environment once $M_\star$ is controlled, indicating that environment-related physical mechanisms (e.g. tidal interaction and ram-pressure stripping) might not significantly affect SMBH growth. We find that $\overline{\mathrm{BHAR}}$ is strongly related to host-galaxy $M_\star$, regardless of environment.

Context: Redshifts are fundamental for our understanding of extragalactic X-ray sources. Ambiguous counterpart associations, expensive optical spectroscopy and/or multimission multiwavelength coverage to resolve degeneracies make estimation often difficult in practice. Aims: We attempt to constrain redshifts of obscured Active Galactic Nuclei (AGN) using only low-resolution X-ray spectra. Methods: Our XZ method fits AGN X-ray spectra with a moderately complex spectral model incorporating a corona, torus obscurer and warm mirror. Using the Bayesian X-ray Astronomy (BXA) package, we constrain redshift, column density, photon index and luminosity simultaneously. The redshift information primarily comes from absorption edges in Compton-thin AGN, and from the Fe K$\alpha$ fluorescent line in heavily obscured AGN. A new generic background fitting method allows us to extract more information from limited numbers of source counts. Results: We derive redshift constraints for 74/321 hard-band detected sources in the Chandra deep field South. Comparing with spectroscopic redshifts, we find an outlier fraction of 8%, indicating that our model assumptions are valid. For three Chandra deep fields, we release our XZ redshift estimates. Conclusions: The independent XZ estimate is easy to apply and effective for a large fraction of obscured AGN in todays deep surveys without the need for any additional data. Comparing to different redshift estimation methods, XZ can resolve degeneracies in photometric redshifts, help to detect potential association problems and confirm uncertain single-line spectroscopic redshifts. With high spectral resolution and large collecting area, this technique will be highly effective for Athena/WFI observations.

Thorium hollow cathode lamps (HCLs) are used as frequency calibrators for many high resolution astronomical spectrographs, some of which aim for Doppler precision at the 1 m/s level. We aim to determine the most suitable combination of elements (Th or U, Ar or Ne) for wavelength calibration of astronomical spectrographs, to characterize differences between similar HCLs, and to provide a new U line-list. We record high resolution spectra of different HCLs using a Fourier transform spectrograph: (i) U-Ne, U-Ar, Th-Ne, and Th-Ar lamps in the spectral range from 500 to 1000 nm and U-Ne and U-Ar from 1000 to 1700 nm; (ii) we systematically compare the number of emission lines and the line intensity ratio for a set of 12 U-Ne HCLs; and (iii) we record a master spectrum of U-Ne to create a new U line-list. Uranium lamps show more lines suitable for calibration than Th lamps from 500 to 1000 nm. The filling gas of the lamps significantly affects their performance because Ar and Ne lines contaminate different spectral regions. We find differences (up to 88 %) in the line intensity of U lines in different lamps from the same batch. We find 8239 isolated lines between 500 and 1700 nm that we attribute to U, 3379 of which were not contained in earlier line-lists. The U line-list is available at the http://www.astro.physik.uni-goettingen.de/research/U_atlas . We suggest using a combination of U-Ne and U-Ar lamps to wavelength-calibrate astronomical spectrographs up to 1000 nm. From 1000 to 1700 nm, U-Ne shows better properties. The differences in line strength between different HCLs underline the importance of characterizing HCLs in the laboratory. The new 3379 U lines can significantly improve the radial velocity precision of astronomical spectrographs.

We analyze the >4-sigma sources in the most sensitive 100 arcmin^2 area (rms <0.56 mJy) of a SCUBA-2 850 micron survey of the GOODS-S and present the 75 band 7 ALMA sources (>4.5-sigma) obtained from high-resolution interferometric follow-up observations. The SCUBA-2---and hence ALMA---samples should be complete to 2.25 mJy. Of the 53 SCUBA-2 sources in this complete sample, only five have no ALMA detections, while 13% (68% confidence range 7-19%) have multiple ALMA counterparts. Color-based high-redshift dusty galaxy selection techniques find at most 55% of the total ALMA sample. In addition to using literature spectroscopic and optical/NIR photometric redshifts, we estimate FIR photometric redshifts based on an Arp 220 template. We identify seven z>4 candidates. We see the expected decline with redshift of the 4.5 micron and 24 micron to 850 micron flux ratios, confirming these as good diagnostics of z>4 candidates. We visually classify 52 ALMA sources, finding 44% (68% confidence range 35-53%) to be apparent mergers. We calculate rest-frame 2-8 keV and 8-28 keV luminosities using the 7 Ms Chandra X-ray image. Nearly all of the ALMA sources detected at 0.5-2 keV are consistent with a known X-ray luminosity to 850 micron flux relation for star-forming galaxies, while most of those detected at 2-7 keV are moderate luminosity AGNs that lie just above the 2-7 keV detection threshold. The latter largely have substantial obscurations of log N_H = 23-24 cm^-2, but two of the high-redshift candidates may even be Compton thick.

We present a measurement of the two-point autocorrelation function of photometrically-selected, high-$z$ quasars over $\sim$ 100 deg$^2$ on the Sloan Digitial Sky Survey Stripe 82 field. Selection is performed using three machine-learning algorithms, trained on known high-$z$ quasar colors, in a six-dimensional, optical/mid-infrared color space. Optical data from the Sloan Digitial Sky Survey is combined with overlapping deep mid-infrared data from the \emph{Spitzer} IRAC Equatorial Survey and the \emph{Spitzer}-HETDEX Exploratory Large-area survey. The selected quasar sample consists of 1378 objects and contains both spectroscopically-confirmed quasars and photometrically-selected quasar candidates. These objects span a redshift range of $2.9 \leq z \leq 5.1$ and are generally fainter than $i=20.2$; a regime which has lacked sufficient number density to perform autocorrelation function measurements of photometrically-classified quasars. We compute the angular correlation function of these data, marginally detecting quasar clustering. We fit a single power-law with an index of $\delta = 1.39 \pm 0.618$ and amplitude of $\theta_0 = 0.71 \pm 0.546$ arcmin. A dark-matter model is fit to the angular correlation function to estimate the linear bias. At the average redshift of our survey ($\langle z \rangle = 3.38$) the bias is $b = 6.78 \pm 1.79$. Using this bias, we calculate a characteristic dark-matter halo mass of 1.70--9.83$\times 10^{12}h^{-1} M_{\odot}$. Our bias estimate suggests that quasar feedback intermittently shuts down the accretion of gas onto the central super-massive black hole at early times. If confirmed, these results hint at a level of luminosity dependence in the clustering of quasars at high-$z$.

A fundamental quest of modern astronomy is to locate the earliest galaxies and study how they influenced the intergalactic medium a few hundred million years after the Big Bang. The abundance of star-forming galaxies is known to decline from redshifts of about 6 to 10, but a key question is the extent of star formation at even earlier times, corresponding to the period when the first galaxies might have emerged. Here we present spectroscopic observations of MACS1149-JD1, a gravitationally lensed galaxy observed when the Universe was less than four per cent of its present age. We detect an emission line of doubly ionized oxygen at a redshift of $9.1096\pm0.0006$, with an uncertainty of one standard deviation. This precisely determined redshift indicates that the red rest-frame optical colour arises from a dominant stellar component that formed about 250 million years after the Big Bang, corresponding to a redshift of about 15. Our results indicate the it may be possible to detect such early episodes of star formation in similar galaxies with future telescopes.

We present an X-ray point-source catalog from the XMM-Large Scale Structure survey region (XMM-LSS), one of the XMM-Spitzer Extragalactic Representative Volume Survey (XMM-SERVS) fields. We target the XMM-LSS region with $1.3$ Ms of new XMM-Newton AO-15 observations, transforming the archival X-ray coverage in this region into a 5.3 deg$^2$ contiguous field with uniform X-ray coverage totaling $2.7$ Ms of flare-filtered exposure, with a $46$ ks median PN exposure time. We provide an X-ray catalog of 5242 sources detected in the soft (0.5-2 keV), hard (2-10 keV), and/or full (0.5-10 keV) bands with a 1% expected spurious fraction determined from simulations. A total of 2381 new X-ray sources are detected compared to previous source catalogs in the same area. Our survey has flux limits of $1.7\times10^{-15}$, $1.3\times10^{-14}$, and $6.5\times10^{-15}$ erg cm$^{-2}$ s$^{-1}$ over 90% of its area in the soft, hard, and full bands, respectively, which is comparable to those of the XMM-COSMOS survey. We identify multiwavelength counterpart candidates for 99.9% of the X-ray sources, of which 93% are considered as reliable based on their matching likelihood ratios. The reliabilities of these high-likelihood-ratio counterparts are further confirmed to be $\approx 97$% reliable based on deep Chandra coverage over $\approx 5$% of the XMM-LSS region. Results of multiwavelength identifications are also included in the source catalog, along with basic optical-to-infrared photometry and spectroscopic redshifts from publicly available surveys. We compute photometric redshifts for X-ray sources in 4.5 deg$^2$ of our field where forced-aperture multi-band photometry is available; $>70$% of the X-ray sources in this subfield have either spectroscopic or high-quality photometric redshifts.

Previous simulations predicted the activity-induced radial-velocity (RV) variations of M dwarfs to range from $\sim1$ cm/s to $\sim1$ km/s, depending on various stellar and activity parameters. We investigate the observed relations between RVs, stellar activity, and stellar parameters of M dwarfs by analyzing CARMENES high-resolution visual-channel spectra ($0.5$$-$$1$$\mu$m), which were taken within the CARMENES RV planet survey during its first $20$ months of operation. During this time, $287$ of the CARMENES-sample stars were observed at least five times. From each spectrum we derived a relative RV and a measure of chromospheric H$\alpha$ emission. In addition, we estimated the chromatic index (CRX) of each spectrum, which is a measure of the RV wavelength dependence. Despite having a median number of only $11$ measurements per star, we show that the RV variations of the stars with RV scatter of $>10$ m/s and a projected rotation velocity $v \sin{i}>2$ km/s are caused mainly by activity. We name these stars `active RV-loud stars' and find their occurrence to increase with spectral type: from $\sim3\%$ for early-type M dwarfs (M$0.0$$-$$2.5$V) through $\sim30\%$ for mid-type M dwarfs (M$3.0$$-$$5.5$V) to $>50\%$ for late-type M dwarfs (M$6.0$$-$$9.0$V). Their RV-scatter amplitude is found to be correlated mainly with $v \sin{i}$. For about half of the stars, we also find a linear RV$-$CRX anticorrelation, which indicates that their activity-induced RV scatter is lower at longer wavelengths. For most of them we can exclude a linear correlation between RV and H$\alpha$ emission. Our results are in agreement with simulated activity-induced RV variations in M dwarfs. The RV variations of most active RV-loud M dwarfs are likely to be caused by dark spots on their surfaces, which move in and out of view as the stars rotate.

We analyse a sample of 8 highly magnified galaxies at redshift 0.6<z<1.5 observed with MUSE, exploring the resolved properties of these galaxies at sub-kiloparsec scales. Combining multi-band HST photometry and MUSE spectra, we derive the stellar mass, global star formation rates, extinction and metallicity from multiple nebular lines, concluding that our sample is representative of z~1 star-forming galaxies. We derive the 2D kinematics of these galaxies from the [OII] emission and model it with a new method that accounts for lensing effects and fits multiple images simultaneously. We use these models to calculate the 2D beam-smearing correction and derive intrinsic velocity dispersion maps. We find them to be fairly homogeneous, with relatively constant velocity dispersions between 15 - 80 km/s and Gini coefficient of <0.3. We do not find any evidence for higher (or lower) velocity dispersions at the positions of bright star-forming clumps. We derive resolved maps of dust attenuation and attenuation-corrected star formation rates from emission lines for two objects in the sample. We use this information to study the relation between resolved star formation rate and velocity dispersion. We find that these quantities are not correlated, and the high-velocity dispersions found for relatively low star-forming densities seems to indicate that, at sub-kiloparsec scales, turbulence in high-z discs is mainly dominated by gravitational instability rather than stellar feedback.

The CARMENES radial velocity (RV) survey is observing 324 M dwarfs to search for any orbiting planets. In this paper, we present the survey sample by publishing one CARMENES spectrum for each M dwarf. These spectra cover the wavelength range 520--1710nm at a resolution of at least $R > 80,000$, and we measure its RV, H$\alpha$ emission, and projected rotation velocity. We present an atlas of high-resolution M-dwarf spectra and compare the spectra to atmospheric models. To quantify the RV precision that can be achieved in low-mass stars over the CARMENES wavelength range, we analyze our empirical information on the RV precision from more than 6500 observations. We compare our high-resolution M-dwarf spectra to atmospheric models where we determine the spectroscopic RV information content, $Q$, and signal-to-noise ratio. We find that for all M-type dwarfs, the highest RV precision can be reached in the wavelength range 700--900nm. Observations at longer wavelengths are equally precise only at the very latest spectral types (M8 and M9). We demonstrate that in this spectroscopic range, the large amount of absorption features compensates for the intrinsic faintness of an M7 star. To reach an RV precision of 1ms$^{-1}$ in very low mass M dwarfs at longer wavelengths likely requires the use of a 10m class telescope. For spectral types M6 and earlier, the combination of a red visual and a near-infrared spectrograph is ideal to search for low-mass planets and to distinguish between planets and stellar variability. At a 4m class telescope, an instrument like CARMENES has the potential to push the RV precision well below the typical jitter level of 3-4ms$^{-1}$.

CARMENES is a spectrograph for radial velocity surveys of M dwarfs with the aim of detecting Earth-mass planets orbiting in the habitable zones of their host stars. To ensure an optimal use of the CARMENES Guaranteed Time Observations, in this paper we investigate the correlation of activity and rotation for approximately 2200 M dwarfs, ranging in spectral type from M0.0 V to M9.0 V. We present new high-resolution spectroscopic observations with FEROS, CAFE, and HRS of approximately 500 M dwarfs. For each new observation, we determined its radial velocity and measured its Halpha activity index and its rotation velocity. Additionally, we have multiple observations of many stars to investigate if there are any radial velocity variations due to multiplicity. The results of our survey confirm that early-M dwarfs are Halpha inactive with low rotational velocities and that late-M dwarfs are Halpha active with very high rotational velocities. The results of this high-resolution analysis comprise the most extensive catalogue of rotation and activity in M dwarfs currently available.

Context: The main goal of the CARMENES survey is to find Earth-mass planets around nearby M-dwarf stars. Seven M-dwarfs included in the CARMENES sample had been observed before with HIRES and HARPS and either were reported to have one short period planetary companion (GJ15A, GJ176, GJ436, GJ536 and GJ1148) or are multiple planetary systems (GJ581 and GJ876). Aims: We aim to report new precise optical radial velocity measurements for these planet hosts and test the overall capabilities of CARMENES. Methods: We combined our CARMENES precise Doppler measurements with those available from HIRES and HARPS and derived new orbital parameters for the systems. Bona-fide single planet systems are fitted with a Keplerian model. The multiple planet systems were analyzed using a self-consistent dynamical model and their best fit orbits were tested for long-term stability. Results: We confirm or provide supportive arguments for planets around all the investigated stars except for GJ15A, for which we find that the post-discovery HIRES data and our CARMENES data do not show a signal at 11.4 days. Although we cannot confirm the super-Earth planet GJ15Ab, we show evidence for a possible long-period ($P_{\rm c}$ = 7025$_{-629}^{+972}$ d) Saturn-mass ($m_{\rm c} \sin i$ = 51.8$_{-5.8}^{+5.5}M_\oplus$) planet around GJ15A. In addition, based on our CARMENES and HIRES data we discover a second planet around GJ1148, for which we estimate a period $P_{\rm c}$ = 532.6$_{-2.5}^{+4.1}$ d, eccentricity $e_{\rm c}$ = 0.34$_{-0.06}^{+0.05}$ and minimum mass $m_{\rm c} \sin i$ = 68.1$_{-2.2}^{+4.9}M_\oplus$. Conclusions: The CARMENES optical radial velocities have similar precision and overall scatter when compared to the Doppler measurements conducted with HARPS and HIRES. We conclude that CARMENES is an instrument that is up to the challenge of discovering rocky planets around low-mass stars.

We present optical and near-IR Integral Field Unit (IFU) and ALMA band 6 observations of the nearby dual Active Galactic Nuclei (AGN) Mrk 463. At a distance of 210 Mpc, and a nuclear separation of $\sim$4 kpc, Mrk 463 is an excellent laboratory to study the gas dynamics, star formation processes and supermassive black hole (SMBH) accretion in a late-stage gas-rich major galaxy merger. The IFU observations reveal a complex morphology, including tidal tails, star-forming clumps, and emission line regions. The optical data, which map the full extent of the merger, show evidence for a biconical outflow and material outflowing at $>$600 km s$^{-1}$, both associated with the Mrk 463E nucleus, together with large scale gradients likely related to the ongoing galaxy merger. We further find an emission line region $\sim$11 kpc south of Mrk 463E that is consistent with being photoionized by an AGN. Compared to the current AGN luminosity, the energy budget of the cloud implies a luminosity drop in Mrk 463E by a factor 3-20 over the last 40,000 years. The ALMA observations of $^{12}$CO(2-1) and adjacent 1mm continuum reveal the presence of $\sim$10$^{9}$M$_\odot$ in molecular gas in the system. The molecular gas shows velocity gradients of $\sim$800 km/s and $\sim$400 km/s around the Mrk 463E and 463W nuclei, respectively. We conclude that in this system the infall of $\sim$100s $M_\odot$/yr of molecular gas is in rough balance with the removal of ionized gas by a biconical outflow being fueled by a relatively small, $<$0.01% of accretion onto each SMBH.

We discuss the spectral analysis of a sample of 63 Active Galactic Nuclei (AGN) detected above a limiting flux of $S(8-24$ keV$)=7\times10^{-14}$ erg/s/cm$^2$ in the multi-tiered NuSTAR Extragalactic Survey program. The sources span a redshift range z=0-2.1 (median $\langle$z$\rangle=$0.58). The spectral analysis is performed over the broad 0.5-24 keV energy range, combining NuSTAR with Chandra and/or XMM-Newton data and employing empirical and physically motivated models. This constitutes the largest sample of AGN selected at $ >10$ keV to be homogeneously spectrally analyzed at these flux levels. We study the distribution of spectral parameters such as photon index, column density ($N_{\rm H}$), reflection parameter ($R$) and 10-40 keV luminosity ($L_{X}$). Heavily obscured ($\log [N_{\rm H}/cm^{-2}]\ge23$) and Compton Thick (CT; $\log [N_{\rm H}/cm^{-2}]\ge24$) AGN constitute $\sim$25% (15-17 sources) and $\sim$2-3% ( 1-2 sources) of the sample, respectively. The observed $N_{\rm H}$ distribution fairly agrees with predictions of Cosmic X-ray Background population synthesis models (CXBPSM). We estimate the intrinsic fraction of AGN as a function of $N_{H}$, accounting for the bias against obscured AGN in a flux-selected sample. The fraction of CT AGN relative to $\log [N_{\rm H}/cm^{-2}]=20-24$ AGN is poorly constrainted, formally in the range 2-56% (90% upper limit of 66%). We derived a fraction ($f_{abs}$) of obscured AGN ($\log [N_{\rm H}/cm^{-2}]=22-24$) as a function of $L_{X}$ in agreement with CXBPSM and previous z < 1 X-ray determinations. Furthermore $f_{abs}$ at z=0.1-0.5 and $\log (L_{X}/erg$ $s^{-1})\approx43.6-44.3$ agrees with observational measurements/trends obtained over larger redshift intervals. We report a significant anti-correlation of $R$ with $L_{X}$ (confirmed by our companion paper on stacked spectra) with considerable scatter around the median $R$ values.

We present 5 years of optical and infrared data of the black hole candidate MAXI J1659-152 covering its 2010 outburst, decay and quiescence. Combining optical data taken during the outburst decay, we obtain an orbital period of 2.414 $\pm$ 0.005 h, in perfect agreement with the value previously measured from X-ray dips. In addition, we detect a clear H$\alpha$ excess in MAXI J1659-152 with data taken during the outburst decay. We also detect a single hump modulation most likely produced by irradiation. Assuming that the maximum occurs at orbital phase 0.5, we constrain the phase of the X-ray dips to be ~ 0.65. We also detect the quiescent optical counterpart at r' = 24.20 $\pm$ 0.08, I = 23.32 $\pm$ 0.02 and H = 20.7 $\pm$ 0.1. These magnitudes provide colour indices implying an M2-M5 donor star assuming 60% contribution from a disc component in the r'-band.

The most heavily-obscured, luminous quasars might represent a specific phase of the evolution of actively accreting supermassive black holes and their host galaxies, possibly related to mergers. We investigated a sample of the most luminous quasars at $z\approx 1-3$ in the GOODS fields, selected in the mid-infrared band through detailed spectral energy distribution (SED) decomposition. The vast majority of these quasars (~80%) are obscured in the X-ray band and ~30% of them to such an extent, that they are undetected in some of the deepest (2 and 4 Ms) Chandra X-ray data. Although no clear relation is found between the star-formation rate of the host galaxies and the X-ray obscuration, we find a higher incidence of heavily-obscured quasars in disturbed/merging galaxies compared to the unobscured ones, thus possibly representing an earlier stage of evolution, after which the system is relaxing and becoming unobscured.

We report on the first star discovered to host a planet detected by radial velocity (RV) observations obtained within the CARMENES survey for exoplanets around M dwarfs. HD 147379 ($V = 8.9$ mag, $M = 0.58 \pm 0.08$ M$_{\odot}$), a bright M0.0V star at a distance of 10.7 pc, is found to undergo periodic RV variations with a semi-amplitude of $K = 5.1\pm0.4$ m s$^{-1}$ and a period of $P = 86.54\pm0.06$ d. The RV signal is found in our CARMENES data, which were taken between 2016 and 2017, and is supported by HIRES/Keck observations that were obtained since 2000. The RV variations are interpreted as resulting from a planet of minimum mass $m_{\rm p}\sin{i} = 25 \pm 2$ M$_{\oplus}$, 1.5 times the mass of Neptune, with an orbital semi-major axis $a = 0.32$ au and low eccentricity ($e < 0.13$). HD 147379b is orbiting inside the temperate zone around the star, where water could exist in liquid form. The RV time-series and various spectroscopic indicators show additional hints of variations at an approximate period of 21.1d (and its first harmonic), which we attribute to the rotation period of the star.

We present one year of optical and near-infrared photometry and spectroscopy of the Type IIP SNe 2013K and 2013am. Both objects are affected by significant extinction, due to their location in dusty regions of their respective host galaxies, ESO 009-10 and NGC 3623 (M65). From the photospheric to nebular phases, these objects display spectra congruent with those of underluminous Type IIP SNe (i.e. the archetypal SNe 1997D or 2005cs), showing low photospheric velocities (~2 X 10**3 km/s at 50d) together with features arising from Ba II which are particularly prominent in faint SNe IIP. The peak V-band magnitudes of SN 2013K (-15.6 mag) and SN 2013am (-16.2 mag) are fainter than standard-luminosity Type IIP SNe. The ejected Nickel masses are 0.012+-0.010 and 0.015+-0.006 Msol for SN 2013K and SN 2013am, respectively. The physical properties of the progenitors at the time of explosion are derived through hydrodynamical modelling. Fitting the bolometric curves, the expansion velocity and the temperature evolution, we infer total ejected masses of 12 and 11.5 Msol, pre-SN radii of ~460 and ~360 Rsol, and explosion energies of 0.34 foe and 0.40 foe for SN 2013K and SN 2013am. Late time spectra are used to estimate the progenitor masses from the strength of nebular emission lines, which turn out to be consistent with red supergiant progenitors of ~15 Msol. For both SNe, a low-energy explosion of a moderate-mass red supergiant star is therefore the favoured scenario.

Hard X-ray ($\geq 10$ keV) observations of Active Galactic Nuclei (AGN) can shed light on some of the most obscured episodes of accretion onto supermassive black holes. The 70-month Swift/BAT all-sky survey, which probes the 14-195 keV energy range, has currently detected 838 AGN. We report here on the broad-band X-ray (0.3-150 keV) characteristics of these AGN, obtained by combining XMM-Newton, Swift/XRT, ASCA, Chandra, and Suzaku observations in the soft X-ray band ($\leq 10$ keV) with 70-month averaged Swift/BAT data. The non-blazar AGN of our sample are almost equally divided into unobscured ($N_{\rm H}< 10^{22}\rm cm^{-2}$) and obscured ($N_{\rm H}\geq 10^{22}\rm cm^{-2}$) AGN, and their Swift/BAT continuum is systematically steeper than the 0.3-10 keV emission, which suggests that the presence of a high-energy cutoff is almost ubiquitous. We discuss the main X-ray spectral parameters obtained, such as the photon index, the reflection parameter, the energy of the cutoff, neutral and ionized absorbers, and the soft excess for both obscured and unobscured AGN.

The SUperluminous Supernova Host galaxIES (SUSHIES) survey aims to provide strong new constraints on the progenitors of superluminous supernovae (SLSNe) by understanding the relationship to their host galaxies. We present the photometric properties of 53 H-poor and 16 H-rich SLSN host galaxies out to $z\sim4$. We model their spectral energy distributions to derive physical properties, which we compare with other galaxy populations. At low redshift, H-poor SLSNe are preferentially found in very blue, low-mass galaxies with high average specific star-formation rates. As redshift increases, the host population follows the general evolution of star-forming galaxies towards more luminous galaxies. After accounting for secular evolution, we find evidence for differential evolution in galaxy mass, but not in the $B$-band and the far UV luminosity ($3\sigma$ confidence). Most remarkable is the scarcity of hosts with stellar masses above $10^{10}~M_\odot$ for both classes of SLSNe. In the case of H-poor SLSNe, we attribute this to a stifled production efficiency above $\sim0.4$ solar metallicity. However, we argue that, in addition to low metallicity, a short-lived stellar population is also required to regulate the SLSN production. H-rich SLSNe are found in a very diverse population of star-forming galaxies. Still, the scarcity of massive hosts suggests a stifled production efficiency above $\sim0.8$ solar metallicity. The large dispersion of the H-rich SLSNe host properties is in stark contrast to those of gamma-ray burst, regular core-collapse SN, and H-poor SLSNe host galaxies. We propose that multiple progenitor channels give rise to this sub-class.

Binary neutron-star mergers (BNSMs) are among the most readily detectable gravitational-wave (GW) sources with LIGO. They are also thought to produce short $\gamma$-ray bursts (SGRBs), and kilonovae that are powered by r-process nuclei. Detecting these phenomena simultaneously would provide an unprecedented view of the physics during and after the merger of two compact objects. Such a Rosetta Stone event was detected by LIGO/Virgo on 17 August 2017 at a distance of $\sim 44$ Mpc. We monitored the position of the BNSM with ALMA at 338.5 GHz and GMRT at 1.4 GHz, from 1.4 to 44 days after the merger. Our observations rule out any afterglow more luminous than $3\times 10^{26}~{\rm erg\,s}^{-1}\,{\rm Hz}^{-1}$ in these bands, probing $>$2--4 dex fainter than previous SGRB limits. We match these limits, in conjunction with public data announcing the appearance of X-ray and radio emission in the weeks after the GW event, to templates of off-axis afterglows. Our broadband modeling suggests that GW170817 was accompanied by a SGRB and that the GRB jet, powered by $E_{\rm AG,\,iso}\sim10^{50}$~erg, had a half-opening angle of $\sim20^\circ$, and was misaligned by $\sim41^\circ$ from our line of sight. The data are also consistent with a more collimated jet: $E_{\rm AG,\,iso}\sim10^{51}$~erg, $\theta_{1/2,\,\rm jet}\sim5^\circ$, $\theta_{\rm obs}\sim17^\circ$. This is the most conclusive detection of an off-axis GRB afterglow and the first associated with a BNSM-GW event to date. Assuming a uniform top-hat jet, we use the viewing angle estimates to infer the initial bulk Lorentz factor and true energy release of the burst.

Context: The CARMENES survey is a high-precision radial velocity (RV) programme that aims to detect Earth-like planets orbiting low-mass stars. Aims: We develop least-squares fitting algorithms to derive the RVs and additional spectral diagnostics implemented in the SpEctrum Radial Velocity Analyser (SERVAL), a publicly available python code. Methods: We measured the RVs using high signal-to-noise templates created by coadding all available spectra of each star.We define the chromatic index as the RV gradient as a function of wavelength with the RVs measured in the echelle orders. Additionally, we computed the differential line width by correlating the fit residuals with the second derivative of the template to track variations in the stellar line width. Results: Using HARPS data, our SERVAL code achieves a RV precision at the level of 1m/s. Applying the chromatic index to CARMENES data of the active star YZ CMi, we identify apparent RV variations induced by stellar activity. The differential line width is found to be an alternative indicator to the commonly used full width half maximum. Conclusions: We find that at the red optical wavelengths (700--900 nm) obtained by the visual channel of CARMENES, the chromatic index is an excellent tool to investigate stellar active regions and to identify and perhaps even correct for activity-induced RV variations.

Gravitational waves were discovered with the detection of binary black hole mergers and they should also be detectable from lower mass neutron star mergers. These are predicted to eject material rich in heavy radioactive isotopes that can power an electromagnetic signal called a kilonova. The gravitational wave source GW170817 arose from a binary neutron star merger in the nearby Universe with a relatively well confined sky position and distance estimate. Here we report observations and physical modelling of a rapidly fading electromagnetic transient in the galaxy NGC4993, which is spatially coincident with GW170817 and a weak short gamma-ray burst. The transient has physical parameters broadly matching the theoretical predictions of blue kilonovae from neutron star mergers. The emitted electromagnetic radiation can be explained with an ejected mass of 0.04 +/- 0.01 Msol, with an opacity of kappa <= 0.5 cm2/gm at a velocity of 0.2 +/- 0.1c. The power source is constrained to have a power law slope of beta = -1.2 +/- 0.3, consistent with radioactive powering from r-process nuclides. We identify line features in the spectra that are consistent with light r-process elements (90 < A < 140). As it fades, the transient rapidly becomes red, and emission may have contribution by a higher opacity, lanthanide-rich ejecta component. This indicates that neutron star mergers produce gravitational waves, radioactively powered kilonovae, and are a nucleosynthetic source of the r-process elements.

We systematically analyze X-ray variability of active galactic nuclei (AGNs) in the 7~Ms \textit{Chandra} Deep Field-South survey. On the longest timescale ($\approx~17$ years), we find only weak (if any) dependence of X-ray variability amplitudes on energy bands or obscuration. We use four different power spectral density (PSD) models to fit the anti-correlation between normalized excess variance ($\sigma^2_{\rm nxv}$) and luminosity, and obtain a best-fit power law index $\beta=1.16^{+0.05}_{-0.05}$ for the low-frequency part of AGN PSD. We also divide the whole light curves into 4 epochs in order to inspect the dependence of $\sigma^2_{\rm nxv}$ on these timescales, finding an overall increasing trend. The analysis of these shorter light curves also infers a $\beta$ of $\sim 1.3$ that is consistent with the above-derived $\beta$, which is larger than the frequently-assumed value of $\beta=1$. We then investigate the evolution of $\sigma^2_{\rm nxv}$. No definitive conclusion is reached due to limited source statistics but, if present, the observed trend goes in the direction of decreasing AGN variability at fixed luminosity toward large redshifts. We also search for transient events and find 6 notable candidate events with our considered criteria. Two of them may be a new type of fast transient events, one of which is reported here for the first time. We therefore estimate a rate of fast outbursts $\langle\dot{N}\rangle = 1.0^{+1.1}_{-0.7}\times 10^{-3}~\rm galaxy^{-1}~yr^{-1}$ and a tidal disruption event~(TDE) rate $\langle\dot{N}_{\rm TDE}\rangle=8.6^{+8.5}_{-4.9}\times 10^{-5}~\rm galaxy^{-1}~yr^{-1}$ assuming the other four long outbursts to be TDEs.

We present a study of the average X-ray spectral properties of the sources detected by the NuSTAR extragalactic survey, comprising observations of the E-CDFS, EGS and COSMOS fields. The sample includes 182 NuSTAR sources (64 detected at 8-24 keV), with 3-24 keV fluxes ranging between $f_{\rm 3-24 keV}\approx10^{-14}$ and $6\times10^{-13}$ erg/cm$^2$/s ($f_{\rm 8-24 keV}\approx3\times10^{-14}-3\times10^{-13}$ erg/cm$^2$/s) and redshifts of $z=0.04-3.21$. We produce composite spectra from the Chandra+NuSTAR data ($E\approx2-40$ keV, rest frame) for all the sources with redshift identifications (95%) and investigate the intrinsic, average spectra of the sources, divided into broad-line (BL) and narrow-line (NL) AGN, and also in different bins of X-ray column density and luminosity. The average power-law photon index for the whole sample is $\Gamma=1.65_{-0.03}^{+0.03}$, flatter than $\Gamma\approx1.8$ typically found for AGN. While the spectral slope of BL and X-ray unabsorbed AGN is consistent with typical values ($\Gamma=1.79_{-0.01}^{+0.01}$), a significant flattening is seen in NL AGN and heavily-absorbed sources ($\Gamma=1.60_{-0.05}^{+0.08}$ and $\Gamma=1.38_{-0.12}^{+0.12}$, respectively), likely due to the effect of absorption and to the contribution from Compton reflection to the high-energy flux (E>10 keV). We find that the typical reflection fraction in our spectra is $R\approx0.5$ (for $\Gamma=1.8$), with a tentative indication of an increase of the reflection strength with column density. While there is no significant evidence for a dependence of the photon index with X-ray luminosity in our sample, we find that $R$ decreases with luminosity, with relatively high levels of reflection ($R\approx1.2$) for $L_{\rm 10-40 keV}<10^{44}$ erg/s and $R\approx0.3$ for $L_{\rm 10-40 keV}>10^{44}$ erg/s AGN, assuming $\Gamma=1.8$.

Most sub-mm emission line studies of galaxies to date have targeted sources with known redshifts where the frequencies of the lines are well constrained. Recent blind line scans circumvent the spectroscopic redshift requirement, which could represent a selection bias. Our aim is to detect emission lines present in continuum oriented observations. The detection of such lines provides spectroscopic redshift and yields properties of the galaxies. We perform a search for emission lines in the ALMA observations of five Frontier Fields clusters and assess the reliability of our detection by associating line candidates with detected galaxies in deep near-infrared imaging. We find 26 significant emission lines candidates, with observed line fluxes between 0.2-4.6 Jy km s$^{-1}$ and velocity dispersions (FWHM) of 25-600 km s$^{-1}$. Nine of these candidates lie nearby to near-infrared sources, boosting their reliability; in six cases the observed line frequency and strength are consistent with expectations given the photometric redshift and properties of the galaxy counterparts. We present redshift identifications, magnifications and molecular gas estimates for the galaxies with identified lines. We show that two of these candidates likely originate from starburst galaxies, one of which is a jellyfish galaxy, while another two are consistent with being main sequence galaxies based in their depletion times. This work highlights the degree to which serendipitous emission lines can be discovered in large mosaic continuum observations when deep ancillary data are available. The low number of high-significance line detections, however, confirms that such surveys are not as optimal as blind line scans. We stress that Monte Carlo simulations should be used to assess the line detections significances, since using the negative noise suffers from stochasticity and incurs larger uncertainties.

The large majority of the accreting supermassive black holes in the Universe are obscured by large columns of gas and dust. The location and evolution of this obscuring material have been the subject of intense research in the past decades, and are still highly debated. A decrease in the covering factor of the circumnuclear material with increasing accretion rates has been found by studies carried out across the electromagnetic spectrum. The origin of this trend has been suggested to be driven either by the increase in the inner radius of the obscuring material with incident luminosity due to the sublimation of dust; by the gravitational potential of the black hole; by radiative feedback; or by the interplay between outflows and inflows. However, the lack of a large, unbiased and complete sample of accreting black holes, with reliable information on gas column density, luminosity and mass, has left the main physical mechanism regulating obscuration unclear. Using a systematic multi-wavelength survey of hard X-ray-selected black holes, here we show that radiation pressure on dusty gas is indeed the main physical mechanism regulating the distribution of the circumnuclear material. Our results imply that the bulk of the obscuring dust and gas in these objects is located within the sphere of influence of the black hole (i.e., a few to tens of parsecs), and that it can be swept away even at low radiative output rates. The main physical driver of the differences between obscured and unobscured accreting black holes is therefore their mass-normalized accretion rate.

Core collapse supernova (CCSN) rates suffer from large uncertainties as many CCSNe exploding in regions of bright background emission and significant dust extinction remain unobserved. Such a shortfall is particularly prominent in luminous infrared galaxies (LIRGs), which have high star formation (and thus CCSN) rates and host bright and crowded nuclear regions, where large extinctions and reduced search detection efficiency likely lead to a significant fraction of CCSNe remaining undiscovered. We present the first results of project SUNBIRD (Supernovae UNmasked By Infra-Red Detection), where we aim to uncover CCSNe that otherwise would remain hidden in the complex nuclear regions of LIRGs, and in this way improve the constraints on the fraction that is missed by optical seeing-limited surveys. We observe in the near-infrared 2.15 {\mu}m $K_s$-band, which is less affected by dust extinction compared to the optical, using the multi-conjugate adaptive optics imager GeMS/GSAOI on Gemini South, allowing us to achieve a spatial resolution that lets us probe close in to the nuclear regions. During our pilot program and subsequent first full year we have discovered three CCSNe and one candidate with projected nuclear offsets as small as 200 pc. When compared to the total sample of LIRG CCSNe discovered in the near-IR and optical, we show that our method is singularly effective in uncovering CCSNe in nuclear regions and we conclude that the majority of CCSNe exploding in LIRGs are not detected as a result of dust obscuration and poor spatial resolution.

The low distance, $z=0.0591$, of GRB 100316D and its association with SN 2010bh represent two important motivations for studying this host galaxy and the GRB's immediate environment with the Integral-Field Spectrographs like VLT/MUSE. Its large field-of-view allows us to create 2D maps of gas metallicity, ionization level, and the star-formation rate distribution maps, as well as to investigate the presence of possible host companions. The host is a late-type dwarf irregular galaxy with multiple star-forming regions and an extended central region with signatures of on-going shock interactions. The GRB site is characterized by the lowest metallicity, the highest star-formation rate and the youngest ($\sim$ 20-30 Myr) stellar population in the galaxy, which suggest a GRB progenitor stellar population with masses up to 20 -- 40 $M_{\odot}$. We note that the GRB site has an offset of $\sim$660pc from the most luminous SF region in the host. The observed SF activity in this galaxy may have been triggered by a relatively recent gravitational encounter between the host and a small undetected ($L_{H\alpha} \leq 10^{36}$ erg/s) companion.

We have discovered heavy obscuration in the dual active galactic nucleus (AGN) in the Swift/Burst Alert Telescope (BAT) source SWIFT J2028.5+2543 using Nuclear Spectroscopic Telescope Array (NuSTAR). While an early XMM-Newton study suggested the emission was mainly from NGC 6921, the superior spatial resolution of NuSTAR above 10 keV resolves the Swift/BAT emission into two sources associated with the nearby galaxies MCG +04-48-002 and NGC 6921 (z = 0.014) with a projected separation of 25.3 kpc (91"). NuSTAR's sensitivity above 10 keV finds both are heavily obscured to Compton-thick levels (N H=(1-2)x10^24 cm-2) and contribute equally to the BAT detection ({L}10-50 {keV}{{int}}= 6x10^42 erg s-1). The observed luminosity of both sources is severely diminished in the 2-10 keV band, illustrating the importance of >10 keV surveys like those with NuSTAR and Swift/BAT. Compared to archival X-ray data, MCG +04-48-002 shows significant variability (>3) between observations. Despite being bright X-ray AGNs, they are difficult to detect using optical emission-line diagnostics because MCG +04-48-002 is identified as a starburst/composite because of the high rates of star formation from a luminous infrared galaxy while NGC 6921 is only classified as a LINER using line detection limits. SWIFT J2028.5+2543 is the first dual AGN resolved above 10 keV and is the second most heavily obscured dual AGN discovered to date in the X-rays other than NGC 6240.

We identify sources with extremely hard X-ray spectra (i.e., with photon indices of Gamma<0.6 in the 13 sq. deg. NuSTAR serendipitous survey, to search for the most highly obscured AGNs detected at >10 keV. Eight extreme NuSTAR sources are identified, and we use the NuSTAR data in combination with lower energy X-ray observations (from Chandra, Swift XRT, and XMM-Newton) to characterize the broad-band (0.5-24 keV) X-ray spectra. We find that all of the extreme sources are highly obscured AGNs, including three robust Compton-thick (CT; N_H > 1.5e24 cm^-2) AGNs at low redshift (z<0.1), and a likely-CT AGN at higher redshift (z=0.16). Most of the extreme sources would not have been identified as highly obscured based on the low energy (<10 keV) X-ray coverage alone. The multiwavelength properties (e.g., optical spectra and X-ray/MIR luminosity ratios) provide further support for the eight sources being significantly obscured. Correcting for absorption, the intrinsic rest-frame 10-40 keV luminosities of the extreme sources cover a broad range, from ~ 5 x 10^42 to 10^45 erg s^-1. The estimated number counts of CT AGNs in the NuSTAR serendipitous survey are in broad agreement with model expectations based on previous X-ray surveys, except for the lowest redshifts (z<0.07) where we measure a high CT fraction of f_CT^obs = 30 (+16 -12) %. For the small sample of CT AGNs, we find a high fraction of galaxy major mergers (50 +/- 33%) compared to control samples of "normal" AGNs.

We study the X-ray variability properties of distant AGNs in the Chandra Deep Field-South region over 17 years, up to $z\sim 4$, and compare them with those predicted by models based on local samples. We use the results of Monte Carlo simulations to account for the biases introduced by the discontinuous sampling and the low-count regime. We confirm that variability is an ubiquitous property of AGNs, with no clear dependence on the density of the environment. The variability properties of high-z AGNs, over different temporal timescales, are most consistent with a Power Spectral Density (PSD) described by a broken (or bending) power-law, similar to nearby AGNs. We confirm the presence of an anti-correlation between luminosity and variability, resulting from the dependence of variability on BH mass and accretion rate. We explore different models, finding that our acceptable solutions predict that BH mass influences the value of the PSD break frequency, while the Eddington ratio $\lambda_{Edd}$ affects the PSD break frequency and, possibly, the PSD amplitude as well. We derive the evolution of the average $\lambda_{Edd}$ as a function of redshift, finding results in agreement with measurements based on different estimators. The large statistical uncertainties make our results consistent with a constant Eddington ratio, although one of our models suggest a possible increase of $\lambda_{Edd}$ with lookback time up to $z\sim 2-3$. We conclude that variability is a viable mean to trace the accretion history of supermassive BHs, whose usefulness will increase with future, wide-field/large effective area X-ray missions.

[abridged] The Hubble and Spitzer Space Telescope surveys of the Frontier Fields (FF) provide extremely deep images around six massive, strong-lensing clusters of galaxies. The ALMA FF survey aims to cover the same fields at 1.1mm, with maps reaching (unlensed) sensitivities of $<$70$\mu$Jy, in order to explore the properties of background dusty star-forming galaxies. We report on the multi-wavelength photometric analysis of all 12 significantly detected ($>$5$\sigma$) sources in the first three FF clusters observed by ALMA, based on data from Hubble and Spitzer, the VLT and Herschel. We measure the total photometry in all available bands and determine the photometric redshifts and the physical properties of the counterparts via SED-fitting. In particular, we carefully estimate the FIR photometry using 1.1mm priors to limit the misidentification of blended FIR counterparts, which strongly affect some flux estimates in previous FIR catalogs. We identify robust near-infrared (NIR) counterparts for all 11 sources with K$_s$ detection, the majority of which are quite red, with eight having $F814W-K_s\gtrsim 4$ and five having $F160W-[4.5]\gtrsim3$. From the FIR point of view, all our objects have $z_{phot}$$\sim$1--3, whereas based on the optical SED one object prefers a high-$z$ solution ($z\geq\ $7). Five objects among our sample have spectroscopic redshifts from the GLASS survey for which we can reproduce their SEDs with existing templates. This verification confirms the validity of our photometric redshift methodology. The mean redshift of our sample is $z_{phot}$=1.99$\pm$0.27. All 1.1mm selected objects are massive (10.0$<\log[M_{\star}(M_{\odot})]<$ 11.5), with high star formation rates ($<\log[SFR(M_{\odot}/yr)]> \approx$1.6) and high dust contents (8.1 $<\log[M_{dust} (M_{\odot})]<$8.8), consistent with previous ALMA surveys.

We present $NuSTAR$ X-ray observations of the active galactic nucleus (AGN) in NGC 7674. The source shows a flat X-ray spectrum, suggesting that it is obscured by Compton-thick gas columns. Based upon long-term flux dimming, previous work suggested the alternate possibility that the source is a recently switched-off AGN with the observed X-rays being the lagged echo from the torus. Our high-quality data show the source to be reflection-dominated in hard X-rays, but with a relatively weak neutral Fe K$\alpha$ emission line (equivalent width [EW] of $\approx$ 0.4 keV) and a strong Fe XXVI ionised line (EW $\approx$ 0.2 keV). We construct an updated long-term X-ray light curve of NGC 7674 and find that the observed 2-10 keV flux has remained constant for the past $\approx$ 20 years, following a high flux state probed by $Ginga$. Light travel time arguments constrain the minimum radius of the reflector to be $\sim$ 3.2 pc under the switched-off AGN scenario, $\approx$ 30 times larger than the expected dust sublimation radius, rendering this possibility unlikely. A patchy Compton-thick AGN (CTAGN) solution is plausible, requiring a minimum line-of-sight column density ($N_{\rm H}$) of 3 $\times$ 10$^{24}$ cm$^{-2}$ at present, and yields an intrinsic 2-10 keV luminosity of (3-5) $\times$ 10$^{43}$ erg s$^{-1}$. Realistic uncertainties span the range of $\approx$ (1-13) $\times$ 10$^{43}$ erg s$^{-1}$. The source has one of the weakest fluorescence lines amongst {\em bona fide} CTAGN, and is potentially a local analogue of bolometrically luminous systems showing complex neutral and ionised Fe emission. It exemplifies the difficulty of identification and proper characterisation of distant CTAGN based on the strength of the neutral Fe K$\alpha$ line.

We present a detailed spectral analysis of the brightest Active Galactic Nuclei (AGN) identified in the 7Ms Chandra Deep Field South (CDF-S) survey over a time span of 16 years. Using a model of an intrinsically absorbed power-law plus reflection, with possible soft excess and narrow Fe K$\alpha$ line, we perform a systematic X-ray spectral analysis, both on the total 7Ms exposure and in four different periods with lengths of 2-21 months. With this approach, we not only present the power-law slopes, column densities $N_H$, observed fluxes, and absorption-corrected 2-10~keV luminosities $L_X$ for our sample of AGNs, but also identify significant spectral variabilities among them on time scales of years. We find that the $N_H$ variabilities can be ascribed to two different types of mechanisms, either flux-driven or flux-independent. We also find that the correlation between the narrow Fe line EW and $N_H$ can be well explained by the continuum suppression with increasing $N_H$. Accounting for the sample incompleteness and bias, we measure the intrinsic distribution of $N_H$ for the CDF-S AGN population and present re-selected subsamples which are complete with respect to $N_H$. The $N_H$-complete subsamples enable us to decouple the dependences of $N_H$ on $L_X$ and on redshift. Combining our data with that from C-COSMOS, we confirm the anti-correlation between the average $N_H$ and $L_X$ of AGN, and find a significant increase of the AGN obscured fraction with redshift at any luminosity. The obscured fraction can be described as $f_{obscured}\thickapprox 0.42\ (1+z)^{0.60}$.

We present X-ray bolometric correction factors, $\kappa_{Bol}$ ($\equiv L_{Bol}/L_X$), for Compton-thick (CT) active galactic nuclei (AGN) with the aim of testing AGN torus models, probing orientation effects, and estimating the bolometric output of the most obscured AGN. We adopt bolometric luminosities, $L_{Bol}$, from literature infrared (IR) torus modeling and compile published intrinsic 2--10 keV X-ray luminosities, $L_{X}$, from X-ray torus modeling of NuSTAR data. Our sample consists of 10 local CT AGN where both of these estimates are available. We test for systematic differences in $\kappa_{Bol}$ values produced when using two widely used IR torus models and two widely used X-ray torus models, finding consistency within the uncertainties. We find that the mean $\kappa_{Bol}$ of our sample in the range $L_{Bol}\approx10^{42}-10^{45}$ erg/s is log$_{10}\kappa_{Bol}=1.44\pm0.12$ with an intrinsic scatter of $\sim0.2$ dex, and that our derived $\kappa_{Bol}$ values are consistent with previously established relationships between $\kappa_{Bol}$ and $L_{Bol}$ and $\kappa_{Bol}$ and Eddington ratio. We investigate if $\kappa_{Bol}$ is dependent on $N_H$ by comparing our results on CT AGN to published results on less-obscured AGN, finding no significant dependence. Since many of our sample are megamaser AGN, known to be viewed edge-on, and furthermore under the assumptions of AGN unification whereby unobscured AGN are viewed face-on, our result implies that the X-ray emitting corona is not strongly anisotropic. Finally, we present $\kappa_{Bol}$ values for CT AGN identified in X-ray surveys as a function of their observed $L_X$, where an estimate of their intrinsic $L_{X}$ is not available, and redshift, useful for estimating the bolometric output of the most obscured AGN across cosmic time.

We investigate the dependence of black-hole accretion rate (BHAR) on host-galaxy star formation rate (SFR) and stellar mass ($M_*$) in the CANDELS/GOODS-South field in the redshift range of $0.5\leq z < 2.0$. Our sample consists of $\approx 18000$ galaxies, allowing us to probe galaxies with $0.1 \lesssim \mathrm{SFR} \lesssim 100\ M_\odot\ \mathrm{yr}^{-1}$ and/or $10^8 \lesssim M_* \lesssim 10^{11}\ M_{\odot}$. We use sample-mean BHAR to approximate long-term average BHAR. Our sample-mean BHARs are derived from the $Chandra$ Deep Field-South 7 Ms observations, while the SFRs and $M_*$ have been estimated by the CANDELS team through SED fitting. The average BHAR is correlated positively with both SFR and $M_*$, and the BHAR-SFR and BHAR-$M_*$ relations can both be described acceptably by linear models with a slope of unity. However, BHAR appears to be correlated more strongly with $M_*$ than SFR. This result indicates that $M_*$ is the primary host-galaxy property related to black-hole growth, and the apparent BHAR-SFR relation is largely a secondary effect due to the star-forming main sequence. Among our sources, massive galaxies ($M_* \gtrsim 10^{10} M_{\odot}$) have significantly higher BHAR/SFR ratios than less-massive galaxies, indicating the former have higher black-hole fueling efficiency and/or higher SMBH occupation fraction than the latter. Our results can naturally explain the observed proportionality between $M_{\rm BH}$ and $M_*$ for local giant ellipticals, and suggest their $M_{\rm BH}/M_*$ is higher than that of local star-forming galaxies. Among local star-forming galaxies, massive systems might have higher $M_{\rm BH}/M_*$ compared to dwarfs.

We present ALMA and ATCA observations of the luminous blue variable \rmc. The radio maps show for the first time the core of the nebula and evidence that the nebula is strongly asymmetric with a Z-pattern shape. Hints of this morphology are also visible in the archival \emph{HST} $\rm H\alpha$ image, which overall resembles the radio emission. The emission mechanism in the outer nebula is optically thin free-free in the radio. At high frequencies, a component of point-source emission appears at the position of the star, up to the ALMA frequencies. The rising flux density distribution ($S_{\nu}\sim \nu^{0.78\pm0.05}$) of this object suggests thermal emission from the ionized stellar wind and indicates a departure from spherical symmetry with $n_{e}(r)\propto r^{-2}$. We examine different scenarios to explain this excess of thermal emission from the wind and show that this can arise from a bipolar outflow, supporting the suggestion by other authors that the stellar wind of \rmc is aspherical. We fit the data with two collimated ionized wind models and we find that the mass-loss rate can be a factor of two or more smaller than in the spherical case. We also fit the photometry obtained by IR space telescopes and deduce that the mid- to far-IR emission must arise from extended, cool ($\sim80\,\rm K$) dust within the outer ionized nebula. Finally we discuss two possible scenarios for the nebular morphology: the canonical single star expanding shell geometry, and a precessing jet model assuming presence of a companion star.

Mergers of galaxies are thought to cause significant gas inflows to the inner parsecs, which can activate rapid accretion onto supermassive black holes (SMBHs), giving rise to Active Galactic Nuclei (AGN). During a significant fraction of this process, SMBHs are predicted to be enshrouded by gas and dust. Studying 52 galactic nuclei in infrared-selected local Luminous and Ultra-luminous infrared galaxies in different merger stages in the hard X-ray band, where radiation is less affected by absorption, we find that the amount of material around SMBHs increases during the last phases of the merger. We find that the fraction of Compton-thick (CT, $N_{\rm\,H}\geq 10^{24}\rm\,cm^{-2}$) AGN in late merger galaxies is higher ($f_{\rm\,CT}=65^{+12}_{-13}\%$) than in local hard X-ray selected AGN ($f_{\rm\,CT}=27\pm 4\%$), and that obscuration reaches its maximum when the nuclei of the two merging galaxies are at a projected distance of $D_{12}\simeq0.4-10.8$ kiloparsecs ($f_{\rm\,CT}=77_{-17}^{+13}\%$). We also find that all AGN of our sample in late merger galaxies have $N_{\rm\,H}> 10^{23}\rm\,cm^{-2}$, which implies that the obscuring material covers $95^{+4}_{-8}\%$ of the X-ray source. These observations show that the material is most effectively funnelled from the galactic scale to the inner tens of parsecs during the late stages of galaxy mergers, and that the close environment of SMBHs in advanced mergers is richer in gas and dust with respect to that of SMBHs in isolated galaxies, and cannot be explained by the classical AGN unification model in which the torus is responsible for the obscuration.

We report on the detailed analysis of a gravitationally-lensed Y-band dropout, A2744_YD4, selected from deep Hubble Space Telescope imaging in the Frontier Field cluster Abell 2744. Band 7 observations with the Atacama Large Millimeter Array (ALMA) indicate the proximate detection of a significant 1mm continuum flux suggesting the presence of dust for a star-forming galaxy with a photometric redshift of $z\simeq8$. Deep X-SHOOTER spectra confirms the high redshift identity of A2744_YD4 via the detection of Lyman $\alpha$ emission at a redshift $z$=8.38. The association with the ALMA detection is confirmed by the presence of [OIII] 88$\mu$m emission at the same redshift. Although both emission features are only significant at the 4 $\sigma$ level, we argue their joint detection and the positional coincidence with a high redshift dropout in the HST images confirms the physical association. Analysis of the available photometric data and the modest gravitational magnification ($\mu\simeq2$) indicates A2744_YD4 has a stellar mass of $\sim$ 2$\times$10$^9$ M$_{\odot}$, a star formation rate of $\sim20$ M$_{\odot}$/yr and a dust mass of $\sim$6$\times$10$^{6}$ M$_{\odot}$. We discuss the implications of the formation of such a dust mass only $\simeq$200 Myr after the onset of cosmic reionisation.

We present a catalog of hard X-ray sources in a square-degree region surveyed by NuSTAR in the direction of the Norma spiral arm. This survey has a total exposure time of 1.7 Ms, and typical and maximum exposure depths of 50 ks and 1 Ms, respectively. In the area of deepest coverage, sensitivity limits of $5\times10^{-14}$ and $4\times10^{-14}$ erg s$^{-1}$ cm$^{-2}$ in the 3-10 and 10-20 keV bands, respectively, are reached. Twenty-eight sources are firmly detected and ten are detected with low significance; eight of the 38 sources are expected to be active galactic nuclei. The three brightest sources were previously identified as a low-mass X-ray binary, high-mass X-ray binary, and pulsar wind nebula. Based on their X-ray properties and multi-wavelength counterparts, we identify the likely nature of the other sources as two colliding wind binaries, three pulsar wind nebulae, a black hole binary, and a plurality of cataclysmic variables (CVs). The CV candidates in the Norma region have plasma temperatures of $\approx$10-20 keV, consistent with the Galactic Ridge X-ray emission spectrum but lower than temperatures of CVs near the Galactic Center. This temperature difference may indicate that the Norma region has a lower fraction of intermediate polars relative to other types of CVs compared to the Galactic Center. The NuSTAR log$N$-log$S$ distribution in the 10-20 keV band is consistent with the distribution measured by Chandra at 2-10 keV if the average source spectrum is assumed to be a thermal model with $kT\approx15$~keV, as observed for the CV candidates.

We present a sample of 10 low-mass active galactic nuclei (AGNs) selected from the 40-month NuSTAR serendipitous survey. The sample is selected to have robust NuSTAR detections at $3 - 24$~keV, to be at $z < 0.3$, and to have optical r-band magnitudes at least 0.5~mag fainter than an $L_\star$ galaxy at its redshift. The median values of absolute magnitude, stellar mass and 2--10 X-ray luminosity of our sample are $\langle M_r\rangle = -20.03$, $\langle M_\star\rangle = 4.6\times10^{9}M_\odot$, and $\langle L_{2-10\mathrm{keV}}\rangle = 3.1\times10^{42}$ erg s$^{-1}$, respectively. Five objects have detectable broad H$\alpha$ emission in their optical spectra, indicating black-hole masses of $(1.1-10.4)\times 10^6 M_\odot$. We find that $30^{+17}_{-10}\%$ of the galaxies in our sample do not show AGN-like optical narrow emission lines, and one of the ten galaxies in our sample, J115851+4243.2, shows evidence for heavy X-ray absorption. This result implies that a non-negligible fraction of low-mass galaxies might harbor accreting massive black holes that are missed by optical spectroscopic surveys and $<10$ keV X-ray surveys. The mid-IR colors of our sample also indicate these optically normal low-mass AGNs cannot be efficiently identified with typical AGN selection criteria based on WISE colors. While the hard ($>10$ keV) X-ray selected low-mass AGN sample size is still limited, our results show that sensitive NuSTAR observations are capable of probing faint hard X-ray emission originating from the nuclei of low-mass galaxies out to moderate redshift ($z<0.3$), thus providing a critical step in understanding AGN demographics in low-mass galaxies.

We report on the detection of a remarkable new fast high-energy transient found in the Chandra Deep Field-South, robustly associated with a faint ($m_{\rm R}=27.5$ mag, $z_{\rm ph}$$\sim$2.2) host in the CANDELS survey. The X-ray event is comprised of 115$^{+12}_{-11}$ net 0.3-7.0 keV counts, with a light curve characterised by a $\approx$100 s rise time, a peak 0.3-10 keV flux of $\approx$5$\times$10$^{-12}$ erg s$^{-1}$ cm$^{-2}$, and a power-law decay time slope of $-1.53\pm0.27$. The average spectral slope is $\Gamma=1.43^{+0.23}_{-0.13}$, with no clear spectral variations. The \hbox{X-ray} and multi-wavelength properties effectively rule out the vast majority of previously observed high-energy transients. A few theoretical possibilities remain: an "orphan" X-ray afterglow from an off-axis short-duration Gamma-ray Burst (GRB) with weak optical emission; a low-luminosity GRB at high redshift with no prompt emission below $\sim$20 keV rest-frame; or a highly beamed Tidal Disruption Event (TDE) involving an intermediate-mass black hole and a white dwarf with little variability. However, none of the above scenarios can completely explain all observed properties. Although large uncertainties exist, the implied rate of such events is comparable to those of orphan and low-luminosity GRBs as well as rare TDEs, implying the discovery of an untapped regime for a known transient class, or a new type of variable phenomena whose nature remains to be determined.

We search for high-redshift dropout galaxies behind the Hubble Frontier Fields (HFF) galaxy cluster MACS J1149.5+2223, a powerful cosmic lens that has revealed a number of unique objects in its field. Using the deep images from the Hubble and Spitzer space telescopes, we find 11 galaxies at z>7 in the MACS J1149.5+2223 cluster field, and 11 in its parallel field. The high-redshift nature of the bright z~9.6 galaxy MACS1149-JD, previously reported by Zheng et al., is further supported by non-detection in the extremely deep optical images from the HFF campaign. With the new photometry, the best photometric redshift solution for MACS1149-JD reduces slightly to z=9.44 +/- 0.12. The young galaxy has an estimated stellar mass of (7 +/- 2)X10E8 Msun, and was formed at z=13.2 +1.9-1.6 when the universe was ~300 Myr old. Data available for the first four HFF clusters have already enabled us to find faint galaxies to an intrinsic magnitude of M(UV) ~ -15.5, approximately a factor of ten deeper than the parallel fields.

IC883 is a luminous infrared galaxy (LIRG) classified as a starburst-active galactic nucleus (AGN) composite. In a previous study we detected a low-luminosity AGN (LLAGN) radio candidate. Here we report on our radio follow-up at three frequencies which provides direct and unequivocal evidence of the AGN activity in IC883. Our analysis of archival X-ray data, together with the detection of a transient radio source with luminosity typical of bright supernovae, give further evidence of the ongoing star formation activity, which dominates the energetics of the system. At sub-parsec scales, the radio nucleus has a core-jet morphology with the jet being a newly ejected component showing a subluminal proper motion of 0.6c-1c. The AGN contributes less than two per cent of the total IR luminosity of the system. The corresponding Eddington factor is ~1E-3, suggesting this is a low-accretion rate engine, as often found in LLAGNs. However, its high bolometric luminosity (~10E44erg/s) agrees better with a normal AGN. This apparent discrepancy may just be an indication of the transition nature of the nucleus from a system dominated by star-formation, to an AGN-dominated system. The nucleus has a strongly inverted spectrum and a turnover at ~4.4GHz, thus qualifying as a candidate for the least luminous (L_5.0GHz ~ 6.3E28erg/s/Hz) and one of the youngest (~3000yr) gigahertz-peaked spectrum (GPS) sources. If the GPS origin for the IC883 nucleus is confirmed, then advanced mergers in the LIRG category are potentially key environments to unveil the evolution of GPS sources into more powerful radio galaxies.

NGC 1448 is one of the nearest luminous galaxies ($L_{8-1000\mu m} >$ 10$^{9} L_{\odot}$) to ours ($z$ $=$ 0.00390), and yet the active galactic nucleus (AGN) it hosts was only recently discovered, in 2009. In this paper, we present an analysis of the nuclear source across three wavebands: mid-infrared (MIR) continuum, optical, and X-rays. We observed the source with the Nuclear Spectroscopic Telescope Array (NuSTAR), and combined this data with archival Chandra data to perform broadband X-ray spectral fitting ($\approx$0.5-40 keV) of the AGN for the first time. Our X-ray spectral analysis reveals that the AGN is buried under a Compton-thick (CT) column of obscuring gas along our line-of-sight, with a column density of $N_{\rm H}$(los) $\gtrsim$ 2.5 $\times$ 10$^{24}$ cm$^{-2}$. The best-fitting torus models measured an intrinsic 2-10 keV luminosity of $L_{2-10\rm{,int}}$ $=$ (3.5-7.6) $\times$ 10$^{40}$ erg s$^{-1}$, making NGC 1448 one of the lowest luminosity CTAGNs known. In addition to the NuSTAR observation, we also performed optical spectroscopy for the nucleus in this edge-on galaxy using the European Southern Observatory New Technology Telescope. We re-classify the optical nuclear spectrum as a Seyfert on the basis of the Baldwin-Philips-Terlevich diagnostic diagrams, thus identifying the AGN at optical wavelengths for the first time. We also present high spatial resolution MIR observations of NGC 1448 with Gemini/T-ReCS, in which a compact nucleus is clearly detected. The absorption-corrected 2-10 keV luminosity measured from our X-ray spectral analysis agrees with that predicted from the optical [OIII]$\lambda$5007\AA\ emission line and the MIR 12$\mu$m continuum, further supporting the CT nature of the AGN.

We present the first full catalog and science results for the NuSTAR serendipitous survey. The catalog incorporates data taken during the first 40 months of NuSTAR operation, which provide ~20Ms of effective exposure time over 331 fields, with an areal coverage of 13 sq deg, and 497 sources detected in total over the 3-24 keV energy range. There are 276 sources with spectroscopic redshifts and classifications, largely resulting from our extensive campaign of ground-based spectroscopic followup. We characterize the overall sample in terms of the X-ray, optical, and infrared source properties. The sample is primarily comprised of active galactic nuclei (AGNs), detected over a large range in redshift from z = 0.002 - 3.4 (median of <z> = 0.56), but also includes 16 spectroscopically confirmed Galactic sources. There is a large range in X-ray flux, from log( f_3-24keV / erg s^-1 cm^-2 ) ~ -14 to -11, and in rest-frame 10-40 keV luminosity, from log( L_10-40keV / erg s^-1 ) ~ 39 to 46, with a median of 44.1. Approximately 79% of the NuSTAR sources have lower energy (<10 keV) X-ray counterparts from XMM-Newton, Chandra, and Swift/XRT. The mid-infrared (MIR) analysis, using WISE all-sky survey data, shows that MIR AGN color selections miss a large fraction of the NuSTAR-selected AGN population, from ~15% at the highest luminosities (Lx > 10^44 erg s^-1) to ~80% at the lowest luminosities (Lx < 10^43 erg s^-1). Our optical spectroscopic analysis finds that the observed fraction of optically obscured AGNs (i.e., the Type 2 fraction) is F_Type2 = 53(+14-15)%, for a well-defined subset of the 8-24 keV selected sample. This is higher, albeit at a low significance level, than the Type 2 fraction measured for redshift- and luminosity-matched AGNs selected by <10 keV X-ray missions.

We present X-ray source catalogs for the $\approx7$ Ms exposure of the Chandra Deep Field-South (CDF-S), which covers a total area of 484.2 arcmin$^2$. Utilizing WAVDETECT for initial source detection and ACIS Extract for photometric extraction and significance assessment, we create a main source catalog containing 1008 sources that are detected in up to three X-ray bands: 0.5-7.0 keV, 0.5-2.0 keV, and 2-7 keV. A supplementary source catalog is also provided including 47 lower-significance sources that have bright ($K_s\le23$) near-infrared counterparts. We identify multiwavelength counterparts for 992 (98.4%) of the main-catalog sources, and we collect redshifts for 986 of these sources, including 653 spectroscopic redshifts and 333 photometric redshifts. Based on the X-ray and multiwavelength properties, we identify 711 active galactic nuclei (AGNs) from the main-catalog sources. Compared to the previous $\approx4$ Ms CDF-S catalogs, 291 of the main-catalog sources are new detections. We have achieved unprecedented X-ray sensitivity with average flux limits over the central $\approx1$ arcmin$^2$ region of $\approx1.9\times10^{-17}$, $6.4\times10^{-18}$, and $2.7\times10^{-17}$ erg cm$^{-2}$ s$^{-1}$ in the three X-ray bands, respectively. We provide cumulative number-count measurements observing, for the first time, that normal galaxies start to dominate the X-ray source population at the faintest 0.5-2.0 keV flux levels. The highest X-ray source density reaches $\approx50\,500$ deg$^{-2}$, and $47\%\pm4\%$ of these sources are AGNs ($\approx23\,900$ deg$^{-2}$).

Hot, Dust-Obscured Galaxies (Hot DOGs), selected from the WISE all sky infrared survey, host some of the most powerful Active Galactic Nuclei (AGN) known, and might represent an important stage in the evolution of galaxies. Most known Hot DOGs are at $z> 1.5$, due in part to a strong bias against identifying them at lower redshift related to the selection criteria. We present a new selection method that identifies 153 Hot DOG candidates at $z\sim 1$, where they are significantly brighter and easier to study. We validate this approach by measuring a redshift $z=1.009$, and an SED similar to higher redshift Hot DOGs for one of these objects, WISE J1036+0449 ($L_{\rm\,Bol}\simeq 8\times 10^{46}\rm\,erg\,s^{-1}$), using data from Keck/LRIS and NIRSPEC, SDSS, and CSO. We find evidence of a broadened component in MgII, which, if due to the gravitational potential of the supermassive black hole, would imply a black hole mass of $M_{\rm\,BH}\simeq 2 \times 10^8 M_{\odot}$, and an Eddington ratio of $\lambda_{\rm\,Edd}\simeq 2.7$. WISE J1036+0449 is the first Hot DOG detected by NuSTAR, and the observations show that the source is heavily obscured, with a column density of $N_{\rm\,H}\simeq(2-15)\times10^{23}\rm\,cm^{-2}$. The source has an intrinsic 2-10 keV luminosity of $\sim 6\times 10^{44}\rm\,erg\,s^{-1}$, a value significantly lower than that expected from the mid-infrared/X-ray correlation. We also find that the other Hot DOGs observed by X-ray facilities show a similar deficiency of X-ray flux. We discuss the origin of the X-ray weakness and the absorption properties of Hot DOGs. Hot DOGs at $z\lesssim1$ could be excellent laboratories to probe the characteristics of the accretion flow and of the X-ray emitting plasma at extreme values of the Eddington ratio.

The "torus" obscurer of Active Galactic Nuclei (AGN) is poorly understood in terms of its density, substructure and physical mechanisms. Large X-ray surveys provide model boundary constraints, for both Compton-thin and Compton-thick levels of obscuration, as obscured fractions are mean covering factors $f_{\text{cov}}$. However, a major remaining uncertainty is host galaxy obscuration. In Paper I we discovered a relation of $N_H \propto M_{\star}^{1/3}$ for the obscuration of galaxy-scale gas. Here we apply this observational relation to the AGN population, and find that galaxy-scale gas is responsible for a luminosity-independent fraction of Compton-thin AGN, but does not produce Compton-thick columns. With the host galaxy obscuration understood, we present a model of the remaining, nuclear obscurer which is consistent with a range of observations. Our radiation-lifted torus model consists of a Compton-thick component ($f_{\text{cov}}\sim35\%$) and a Compton-thin component ($f_{\text{cov}}\sim40\%$), which depends on both black hole mass and luminosity. This provides a useful summary of observational constraints for torus modellers who attempt to reproduce this behaviour. It can also be employed as a sub-grid recipe in cosmological simulations which do not resolve the torus. We also investigate host-galaxy X-ray obscuration inside cosmological, hydro-dynamic simulations (EAGLE, Illustris). The obscuration from ray-traced galaxy gas can agree with observations, but is highly sensitive to the chosen feedback assumptions.

An important constraint for galaxy evolution models is how much gas resides in galaxies, in particular at the peak of star formation z=1-3. We attempt a novel approach by letting long-duration Gamma Ray Bursts (LGRBs) x-ray their host galaxies and deliver column densities to us. This requires a good understanding of the obscurer and biases introduced by incomplete follow-up observations. We analyse the X-ray afterglow of all 844 Swift LGRBs to date for their column density $N_H$. To derive the population properties we propagate all uncertainties in a consistent Bayesian methodology. The $N_H$ distribution covers the $10^{20-23}\mathrm{cm}^{-2}$ range and shows no evolutionary effect. Higher obscurations, e.g. Compton-thick columns, could have been detected but are not observed. The $N_H$ distribution is consistent with sources randomly populating a ellipsoidal gas cloud of major axis $N_H^\text{major}=10^{23}\mathrm{cm}^{-2}$ with 0.22 dex intrinsic scatter between objects. The unbiased SHOALS survey of afterglows and hosts allows us to constrain the relation between Spitzer-derived stellar masses and X-ray derived column densities $N_H$. We find a well-constrained powerlaw relation of $N_H=10^{21.7}\mathrm{cm}^{-2}\times\left(M_{\star}/10^{9.5}M_{\odot}\right)^{1/3}$, with 0.5 dex intrinsic scatter between objects. The Milky Way and the Magellanic clouds also follow this relation. From the geometry of the obscurer, its stellar mass dependence and comparison with local galaxies we conclude that LGRBs are primarily obscured by galaxy-scale gas. Ray tracing of simulated Illustris galaxies reveals a relation of the same normalisation, but a steeper stellar-mass dependence and mild redshift evolution. Our new approach provides valuable insight into the gas residing in high-redshift galaxies.

We analyse high-quality NuSTAR observations of the local (z = 0.011) Seyfert 2 active galactic nucleus (AGN) IC 3639, in conjunction with archival Suzaku and Chandra data. This provides the first broadband X-ray spectral analysis of the source, spanning nearly two decades in energy (0.5-30 keV). Previous X-ray observations of the source below 10 keV indicated strong reflection/obscuration on the basis of a pronounced iron fluorescence line at 6.4 keV. The hard X-ray energy coverage of NuSTAR, together with self-consistent toroidal reprocessing models, enables direct broadband constraints on the obscuring column density of the source. We find the source to be heavily Compton-thick (CTK) with an obscuring column in excess of $3.6\times10^{24}$ cm$^{-2}$, unconstrained at the upper end. We further find an intrinsic 2-10 keV luminosity of $\textrm{log}_{10}(L_{\textrm{2-10 keV}} \textrm{[erg s}^{-1}]) = 43.4^{+0.6}_{-1.1}$ to 90% confidence, almost 400 times the observed flux, and consistent with various multi-wavelength diagnostics. Such a high intrinsic to observed flux ratio in addition to an Fe-K$\alpha$ fluorescence line equivalent width exceeding 2 keV is extreme amongst known bona fide CTK AGN, which we suggest are both due to the high level of obscuration present around IC 3639. Our study demonstrates that broadband spectroscopic modelling with NuSTAR enables large corrections for obscuration to be carried out robustly, and emphasises the need for improved modelling of AGN tori showing intense iron fluorescence.

Dusty star-forming galaxies are among the most prodigious systems at high redshift (z>1), characterized by high star formation rates and huge dust reservoirs. The bright end of this population has been well characterized in recent years, but considerable uncertainties remain for fainter dusty star-forming galaxies, which are responsible for the bulk of star formation at high redshift and thus play a key role in galaxy growth and evolution. In this first paper of our series, we describe our methods for finding high redshift faint dusty galaxies using millimeter observations with ALMA. We obtained ALMA 1.1 mm mosaic images for three strong-lensing galaxy clusters from the Frontier Fields survey. The 2'x2' mosaics overlap with the deep HST WFC3/IR footprints and encompass the high magnification regions of each cluster. The combination of extremely high ALMA sensitivity and the magnification power of these clusters allows us to systematically probe the sub-mJy population of dusty star-forming galaxies over a large surveyed area. We present a description of the reduction and analysis of the ALMA continuum observations for the galaxy clusters Abell 2744 (z=0.308), MACSJ0416.1-2403 (z=0.396) and MACSJ1149.5+2223 (z=0.543), for which we reach observed rms sensitivities of 55, 59 and 71 $\mu$Jy beam$^{-1}$ respectively. We detect 12 dusty star-forming galaxies at S/N>5.0 across the three clusters, all of them presenting coincidence with NIR detected counterparts in the HST images. None of the sources fall close to the lensing caustics, thus they are not strongly lensed. The observed 1.1 mm flux densities for the total sample of galaxies range from 0.41 to 2.82 mJy, with observed effective radii spanning <0.05" to 0.37"$\pm$0.21". The lensing-corrected sizes of the detected sources appear to be in the same range as those measured in brighter samples, albeit with possibly larger dispersion.

Context. Recent NuSTAR observations revealed a somewhat unexpected increase in the X-ray flux of the nucleus of NGC 1068. We expect the infrared emission of the dusty torus to react on the intrinsic changes of the accretion disk. Aims. We aim to investigate the origin of the X-ray variation by investigating the response of the mid-infrared environment. Methods. We obtained single-aperture and interferometric mid-infrared measurements and directly compared the measurements observed before and immediately after the X-ray variations. The average correlated and single-aperture fluxes as well as the differential phases were directly compared to detect a possible change in the structure of the nuclear emission on scales of $\sim$ 2 pc. Results. The flux densities and differential phases of the observations before and during the X-ray variation show no significant change over a period of ten years. Possible minor variations in the infrared emission are $\lesssim$ 8 %. Conclusions. Our results suggest that the mid-infrared environment of NGC 1068 has remained unchanged for a decade. The recent transient change in the X-rays did not cause a significant variation in the infrared emission. This independent study supports previous conclusions that stated that the X-ray variation detected by NuSTAR observations is due to X-ray emission piercing through a patchy section of the dusty region.

$Context:$ Intermediate-Mass Black Holes (IMBHs) are thought to be the seeds of early Supermassive Black Holes (SMBHs). While $\gtrsim$100 IMBH and small SMBH candidates have been identified in recent years, few have been robustly confirmed to date, leaving their number density in considerable doubt. Placing firmer constraints both on the methods used to identify and confirm IMBHs/SMBHs, as well as characterizing the range of host environments that IMBHs/SMBHs likely inhabit is therefore of considerable interest and importance. Additionally, finding significant numbers of IMBHs in metal-poor systems would be particularly intriguing, since such systems may represent local analogs of primordial galaxies, and therefore could provide clues of early accretion processes. $Aims:$ Here we study in detail several candidate Active Galactic Nuclei (AGN) found in metal-poor hosts. $Methods:$ We utilize new X-ray and optical observations to characterize these metal-poor AGN candidates and compare them against known AGN luminosity relations and well-characterized IMBH/SMBH samples. $Results:$ Despite having clear broad optical emission lines that are long-lived ($\gtrsim$10--13\,yr), these candidate AGN appear to lack associated strong X-ray and hard UV emission, lying at least 1--2 dex off the known AGN correlations. If they are IMBHs/SMBHs, our constraints imply that they either are not actively accreting, their accretion disks are fully obscured along our line-of-sight, or their accretion disks are not producing characteristic high energy emission. Alternatively, if they are not AGN, then their luminous broad emission lines imply production by extreme stellar processes. The latter would have profound implications on the applicability of broad lines for mass estimates of massive black holes.

We study the molecular gas properties of high-$z$ galaxies observed in the ALMA Spectroscopic Survey (ASPECS) that targets a $\sim1$ arcmin$^2$ region in the Hubble Ultra Deep Field (UDF), a blind survey of CO emission (tracing molecular gas) in the 3mm and 1mm bands. Of a total of 1302 galaxies in the field, 56 have spectroscopic redshifts and correspondingly well-defined physical properties. Among these, 11 have infrared luminosities $L_{\rm{}IR}>10^{11}$ L$_\odot$, i.e. a detection in CO emission was expected. Out these, 7 are detected at various significance in CO, and 4 are undetected in CO emission. In the CO-detected sources, we find CO excitation conditions that are lower than typically found in starburst/SMG/QSO environments. We use the CO luminosities (including limits for non-detections) to derive molecular gas masses. We discuss our findings in context of previous molecular gas observations at high redshift (star-formation law, gas depletion times, gas fractions): The CO-detected galaxies in the UDF tend to reside on the low-$L_{\rm{}IR}$ envelope of the scatter in the $L_{\rm{}IR}-L'_{\rm{}CO}$ relation, but exceptions exist. For the CO-detected sources, we find an average depletion time of $\sim$ 1 Gyr, with significant scatter. The average molecular-to-stellar mass ratio ($M_{\rm{}H2}$/$M_*$) is consistent with earlier measurements of main sequence galaxies at these redshifts, and again shows large variations among sources. In some cases, we also measure dust continuum emission. On average, the dust-based estimates of the molecular gas are a factor $\sim$2-5$\times$ smaller than those based on CO. Accounting for detections as well as non-detections, we find large diversity in the molecular gas properties of the high-redshift galaxies covered by ASPECS.

In this paper we use ASPECS, the ALMA Spectroscopic Survey in the {\em Hubble} Ultra Deep Field (UDF) in band 3 and band 6, to place blind constraints on the CO luminosity function and the evolution of the cosmic molecular gas density as a function of redshift up to $z\sim 4.5$. This study is based on galaxies that have been solely selected through their CO emission and not through any other property. In all of the redshift bins the ASPECS measurements reach the predicted `knee' of the CO luminosity function (around $5\times10^{9}$ K km/s pc$^2$). We find clear evidence of an evolution in the CO luminosity function with respect to $z\sim 0$, with more CO luminous galaxies present at $z\sim 2$. The observed galaxies at $z\sim 2$ also appear more gas-rich than predicted by recent semi-analytical models. The comoving cosmic molecular gas density within galaxies as a function of redshift shows a factor 3-10 drop from $z \sim 2$ to $z \sim 0$ (with significant error bars), and possibly a decline at $z>3$. This trend is similar to the observed evolution of the cosmic star formation rate density. The latter therefore appears to be at least partly driven by the increased availability of molecular gas reservoirs at the peak of cosmic star formation ($z\sim2$).

We present the 3-8 keV and 8-24 keV number counts of active galactic nuclei (AGN) identified in the NuSTAR extragalactic surveys. NuSTAR has now resolved 33-39% of the X-ray background in the 8-24 keV band, directly identifying AGN with obscuring columns up to approximately 1e25 /cm2. In the softer 3-8 keV band the number counts are in general agreement with those measured by XMM-Newton and Chandra over the flux range 5e-15 < S(3 - 8 keV)/(erg/cm2/s) < 1e-12 probed by NuSTAR. In the hard 8-24 keV band NuSTAR probes fluxes over the range 2e-14 < S(8-24 keV)/(erg/cm2/s) < 1e-12, a factor of approximately 100 fainter than previous measurements. The 8-24 keV number counts match predictions from AGN population synthesis models, directly confirming the existence of a population of obscured and/or hard X-ray sources inferred from the shape of the integrated cosmic X-ray background. The measured NuSTAR counts lie significantly above simple extrapolation with a Euclidian slope to low flux of the Swift/BAT 15-55 keV number counts measured at higher fluxes S(15-55 keV) > 1e-11 erg/cm2/s, reflecting the evolution of the AGN population between the Swift/BAT local (z<0.1) sample and NuSTAR's z~1 sample. CXB synthesis models, which account for AGN evolution, lie above the Swift/BAT measurements, suggesting that they do not fully capture the evolution of obscured AGN at low redshifts.

We present an analysis of a deep (1$\sigma$=13 $\mu$Jy) cosmological 1.2-mm continuum map based on ASPECS, the ALMA Spectroscopic Survey in the Hubble Ultra Deep Field. In the 1 arcmin$^2$ covered by ASPECS we detect nine sources at $>3.5\sigma$ significance at 1.2-mm. Our ALMA--selected sample has a median redshift of $z=1.6\pm0.4$, with only one galaxy detected at z$>$2 within the survey area. This value is significantly lower than that found in millimeter samples selected at a higher flux density cut-off and similar frequencies. Most galaxies have specific star formation rates similar to that of main sequence galaxies at the same epoch, and we find median values of stellar mass and star formation rates of $4.0\times10^{10}\ M_\odot$ and $\sim40~M_\odot$ yr$^{-1}$, respectively. Using the dust emission as a tracer for the ISM mass, we derive depletion times that are typically longer than 300 Myr, and we find molecular gas fractions ranging from $\sim$0.1 to 1.0. As noted by previous studies, these values are lower than using CO--based ISM estimates by a factor $\sim$2. The 1\,mm number counts (corrected for fidelity and completeness) are in agreement with previous studies that were typically restricted to brighter sources. With our individual detections only, we recover $55\pm4\%$ of the extragalactic background light (EBL) at 1.2 mm measured by the Planck satellite, and we recover $80\pm7\%$ of this EBL if we include the bright end of the number counts and additional detections from stacking. The stacked contribution is dominated by galaxies at $z\sim1-2$, with stellar masses of (1-3)$\times$10$^{10}$ M$_\odot$. For the first time, we are able to characterize the population of galaxies that dominate the EBL at 1.2 mm.

On March 20, 2015, we obtained 159 spectra of the Sun as a star with the solar telescope and the FTS at the Institut f\"ur Astrophysik G\"ottingen, 76 spectra were taken during partial solar eclipse. We obtained RVs using $I_2$ as wavelength reference and determined the RM curve with a peak-to-peak amplitude of almost 1.4 km s$^{-1}$ at typical RV precision better than 1 m s$^{-1}$. We modeled disk-integrated solar RVs using surface velocities, limb darkening, and information about convective blueshift from 3D magneto-hydrodynamic simulations. We confirm that convective blueshift is crucial to understand solar RVs during eclipse. Our best model reproduced the observations to within a relative precision of 10% with residuals less than 30 m s$^{-1}$. We cross-checked parameterizations of velocity fields using a Dopplergram from the Solar Dynamics Observatory and conclude that disk-integration of the Dopplergram does not provide correct information about convective blueshift necessary for m s$^{-1}$ RV work. As main limitation for modeling RVs during eclipses, we identified limited knowledge about convective blueshift and line shape as functions of solar limb angle. We suspect that our model line profiles are too shallow at limb angles larger than $\mu = 0.6$ resulting in incorrect weighting of the velocities across the solar disk. Alternative explanations cannot be excluded like suppression of convection in magnetic areas and undiscovered systematics during eclipse observations. Accurate observations of solar line profiles across the solar disk are suggested. We publish our RVs taken during solar eclipse as a benchmark curve for codes calculating the RM effect and for models of solar surface velocities and line profiles.

Aims. We study the long-term variability of the well-known Seyfert 2 galaxy Mrk 1210 (a.k.a. UGC 4203, or the Phoenix galaxy). Methods. The source was observed by many X-ray facilities in the last 20 years. Here we present a NuSTAR observation and put the results in context of previously published observations. Results. NuSTAR observed Mrk 1210 in 2012 for 15.4 ks. The source showed Compton-thin obscuration similar to that observed by Chandra, Suzaku, BeppoSAX and XMM-Newton over the past two decades, but different from the first observation by ASCA in 1995, in which the active nucleus was caught in a low flux state - or obscured by Compton-thick matter, with a reflection-dominated spectrum. Thanks to the high-quality hard X-ray spectrum obtained with NuSTAR and exploiting the long-term spectral coverage spanning 16.9 years, we can precisely disentangle the transmission and reflection components and put constraints on both the intrinsic long-term variability and hidden nucleus scenarios. In the former case, the distance between the reflector and the source must be at least ~ 2 pc, while in the latter one the eclipsing cloud may be identified with a water maser-emitting clump.

We perform long-term ($\approx 15$ yr, observed-frame) X-ray variability analyses of the 68 brightest radio-quiet active galactic nuclei (AGNs) in the 6 Ms $Chandra$ Deep Field-South (CDF-S) survey; the majority are in the redshift range of $0.6-3.1$, providing access to penetrating rest-frame X-rays up to $\approx 10-30$ keV. Twenty-four of the 68 sources are optical spectral type I AGNs, and the rest (44) are type II AGNs. The time scales probed in this work are among the longest for X-ray variability studies of distant AGNs. Photometric analyses reveal widespread photon-flux variability: $90\%$ of AGNs are variable above a 95% confidence level, including many X-ray obscured AGNs and several optically classified type II quasars. We characterize the intrinsic X-ray luminosity ($L_{\rm{X}}$) and absorption ($N_{\rm{H}}$) variability via spectral fitting. Most (74%) sources show $L_{\rm{X}}$ variability; the variability amplitudes are generally smaller for quasars. A Compton-thick candidate AGN shows variability of its high-energy X-ray flux, indicating the size of reflecting material to be $\lesssim 0.3$ pc. $L_{\rm{X}}$ variability is also detected in a broad absorption line (BAL) quasar. The $N_{\rm{H}}$ variability amplitude for our sample appears to rise as time separation increases. About 16% of sources show $N_{\rm{H}}$ variability. One source transitions from an X-ray unobscured to obscured state while its optical classification remains type I; this behavior indicates the X-ray eclipsing material is not large enough to obscure the whole broad-line region.

We present a search for [CII] line and dust continuum emission from optical dropout galaxies at $z>6$ using ASPECS, our ALMA Spectroscopic Survey in the Hubble Ultra-Deep Field (UDF). Our observations, which cover the frequency range $212-272$ GHz, encompass approximately the range $6<z<8$ for [CII] line emission and reach a limiting luminosity of L$_{\rm [CII]}\sim$(1.6-2.5)$\times$10$^{8}$ L$_{\odot}$. We identify fourteen [CII] line emitting candidates in this redshift range with significances $>$4.5 $\sigma$, two of which correspond to blind detections with no optical counterparts. At this significance level, our statistical analysis shows that about 60\% of our candidates are expected to be spurious. For one of our blindly selected [CII] line candidates, we tentatively detect the CO(6-5) line in our parallel 3-mm line scan. None of the line candidates are individually detected in the 1.2 mm continuum. A stack of all [CII] candidates results in a tentative detection with $S_{1.2mm}=14\pm5\mu$Jy. This implies a dust-obscured star formation rate (SFR) of $(3\pm1)$ M$_\odot$ yr$^{-1}$. We find that the two highest--SFR objects have candidate [CII] lines with luminosities that are consistent with the low-redshift $L_{\rm [CII]}$ vs. SFR relation. The other candidates have significantly higher [CII] luminosities than expected from their UV--based SFR. At the current sensitivity it is unclear whether the majority of these sources are intrinsically bright [CII] emitters, or spurious sources. If only one of our line candidates was real (a scenario greatly favored by our statistical analysis), we find a source density for [CII] emitters at $6<z<8$ that is significantly higher than predicted by current models and some extrapolations from galaxies in the local universe.

We exploit the 7 Ms \textit{Chandra} observations in the \chandra\,Deep Field-South (\mbox{CDF-S}), the deepest X-ray survey to date, coupled with CANDELS/GOODS-S data, to measure the total X-ray emission arising from 2076 galaxies at $3.5\leq z < 6.5$. This aim is achieved by stacking the \textit{Chandra} data at the positions of optically selected galaxies, reaching effective exposure times of $\geq10^9\mathrm{s}$. We detect significant ($>3.7\sigma$) X-ray emission from massive galaxies at $z\approx4$. We also report the detection of massive galaxies at $z\approx5$ at a $99.7\%$ confidence level ($2.7\sigma$), the highest significance ever obtained for X-ray emission from galaxies at such high redshifts. No significant signal is detected from galaxies at even higher redshifts. The stacking results place constraints on the BHAD associated with the known high-redshift galaxy samples, as well as on the SFRD at high redshift, assuming a range of prescriptions for X-ray emission due to X- ray binaries. We find that the X-ray emission from our sample is likely dominated by processes related to star formation. Our results show that low-rate mass accretion onto SMBHs in individually X-ray-undetected galaxies is negligible, compared with the BHAD measured for samples of X-ray detected AGN, for cosmic SMBH mass assembly at high redshift. We also place, for the first time, constraints on the faint-end of the AGN X-ray luminosity function ($\mathrm{logL_X\sim42}$) at $z>4$, with evidence for fairly flat slopes. The implications of all of these findings are discussed in the context of the evolution of the AGN population at high redshift.

We make use of deep 1.2mm-continuum observations (12.7microJy/beam RMS) of a 1 arcmin^2 region in the Hubble Ultra Deep Field to probe dust-enshrouded star formation from 330 Lyman-break galaxies spanning the redshift range z=2-10 (to ~2-3 Msol/yr at 1sigma over the entire range). Given the depth and area of ASPECS, we would expect to tentatively detect 35 galaxies extrapolating the Meurer z~0 IRX-beta relation to z>~2 (assuming T_d~35 K). However, only 6 tentative detections are found at z>~2 in ASPECS, with just three at >3sigma. Subdividing z=2-10 galaxies according to stellar mass, UV luminosity, and UV-continuum slope and stacking the results, we only find a significant detection in the most massive (>10^9.75 Msol) subsample, with an infrared excess (IRX=L_{IR}/L_{UV}) consistent with previous z~2 results. However, the infrared excess we measure from our large selection of sub-L* (<10^9.75 Msol) galaxies is 0.11(-0.42)(+0.32) and 0.14(-0.14)(+0.15) at z=2-3 and z=4-10, respectively, lying below even an SMC IRX-beta relation (95% confidence). These results demonstrate the relevance of stellar mass for predicting the IR luminosity of z>~2 galaxies. We furthermore find that the evolution of the IRX-stellar mass relationship depends on the evolution of the dust temperature. If the dust temperature increases monotonically with redshift (as (1+z)^0.32) such that T_d~44-50 K at z>=4, current results are suggestive of little evolution in this relationship to z~6. We use these results to revisit recent estimates of the z>~3 SFR density. One less obvious implication is in interpreting the high Halpha EWs seen in z~5 galaxies: our results imply that star-forming galaxies produce Lyman-continuum photons at twice the efficiency (per unit UV luminosity) as implied in conventional models. Star-forming galaxies can then reionize the Universe, even if the escape fraction is <10%.

We present the rationale for and the observational description of ASPECS: The ALMA SPECtroscopic Survey in the Hubble Ultra-Deep Field (UDF), the cosmological deep field that has the deepest multi-wavelength data available. Our overarching goal is to obtain an unbiased census of molecular gas and dust continuum emission in high-redshift (z$>$0.5) galaxies. The $\sim$1$'$ region covered within the UDF was chosen to overlap with the deepest available imaging from HST. Our ALMA observations consist of full frequency scans in band 3 (84-115 GHz) and band 6 (212-272 GHz) at approximately uniform line sensitivity ($L'_{\rm CO}\sim$2$\times$10$^{9}$ K km/s pc$^2$), and continuum noise levels of 3.8 $\mu$Jy beam$^{-1}$ and 12.7 $\mu$Jy beam$^{-1}$, respectively. The molecular surveys cover the different rotational transitions of the CO molecule, leading to essentially full redshift coverage. The [CII] emission line is also covered at redshifts $6.0<z<8.0$. We present a customized algorithm to identify line candidates in the molecular line scans, and quantify our ability to recover artificial sources from our data. Based on whether multiple CO lines are detected, and whether optical spectroscopic redshifts as well as optical counterparts exist, we constrain the most likely line identification. We report 10 (11) CO line candidates in the 3mm (1mm) band, and our statistical analysis shows that $<$4 of these (in each band) are likely spurious. Less than 1/3 of the total CO flux in the low-J CO line candidates are from sources that are not associated with an optical/NIR counterpart. We also present continuum maps of both the band 3 and band 6 observations. The data presented here form the basis of a number of dedicated studies that are presented in subsequent papers.

SN 2005kd is among the most luminous supernovae (SNe) to be discovered at X-ray wavelengths. We have re-analysed all good angular resolution (better than $20"$ FWHM PSF) archival X-ray data for SN 2005kd. The data reveal an X-ray light curve that decreases as t$^{-1.62 \pm 0.06}$. Our modelling of the data suggests that the early evolution is dominated by emission from the forward shock in a high-density medium. Emission from the radiative reverse shock is absorbed by the cold dense shell formed behind the reverse shock. Our results suggest a progenitor with a mass-loss rate towards the end of its evolution of $\ge$ 4.3 $\times$ 10$^{-4} M_{\odot} \,{\rm yr}^{-1}$, for a wind velocity of 10 km s$^{-1}$, at 4.0 $\times$ 10$^{16}$ cm. This mass-loss rate is too high for most known stars, except perhaps hypergiant stars. A higher wind velocity would lead to a correspondingly higher mass-loss rate. A Luminous Blue Variable star undergoing a giant eruption could potentially fulfill this requirement, but would need a high mass-loss rate lasting for several hundred years, and need to explain the plateau observed in the optical light curve. The latter could perhaps be due to the ejecta expanding in the dense circumstellar material at relatively small radii. These observations are consistent with the fact that Type IIn SNe appear to expand into high density and high mass-loss rate environments, and also suggest rapid variability in the wind mass-loss parameters within at least the last 5000 years of stellar evolution prior to core-collapse.

Very massive stars in the final phases of their lives often show unpredictable outbursts that can mimic supernovae, so-called, "SN impostors", but the distinction is not always straigthforward. Here we present observations of a luminous blue variable (LBV) in NGC 2770 in outburst over more than 20 years that experienced a possible terminal explosion as type IIn SN in 2015, named SN 2015bh. This possible SN or "main event" was preceded by a precursor peaking $\sim$ 40 days before maximum. The total energy release of the main event is $\sim$1.8$\times$10$^{49}$ erg, which can be modeled by a $<$ 0.5 M$_\odot$ shell plunging into a dense CSM. All emission lines show a single narrow P-Cygni profile during the LBV phase and a double P-Cygni profile post maximum suggesting an association of this second component with the possible SN. Since 1994 the star has been redder than during a typical S-Dor like outburst. SN 2015bh lies within a spiral arm of NGC 2770 next to a number of small star-forming regions with a metallicity of $\sim$ 0.5 solar and a stellar population age of 7-10 Myr. SN 2015bh shares many similarities with SN 2009ip, which, together with other examples may form a new class of objects that exhibit outbursts a few decades prior to "hyper-eruption" or final core-collapse. If the star survives this event it is undoubtedly altered, and we suggest that these "zombie stars" may evolve from an LBV to a Wolf Rayet star over a very short timescale of only a few years. The final fate of these types of massive stars can only be determined with observations years after the possible SN.

We study the geometry of the AGN obscurer in IRAS 09104+4109, an IR-luminous, radio-intermediate FR-I source at $z=0.442$, using infrared data from Spitzer and Herschel, X-ray data from Nustar, Swift, Suzaku, and Chandra, and an optical spectrum from Palomar. The infrared data imply a total rest-frame 1-1000$\mu$m luminosity of $5.5\times10^{46}$erg s$^{-1}$ and require both an AGN torus and starburst model. The AGN torus has an anisotropy-corrected IR luminosity of $4.9\times10^{46}$erg s$^{-1}$, and a viewing angle and half opening angle both of approximately $36$ degrees from pole-on. The starburst has a star formation rate of $(110\pm34)$M$_{\odot}$ yr$^{-1}$ and an age of $<50$Myr. These results are consistent with two epochs of luminous activity in 09104: one approximately $150$Myr ago, and one ongoing. The X-ray data suggest a photon index of $\Gamma \simeq 1.8$ and a line-of-sight column of $N_{\rm H} \simeq 5\times10^{23}$cm$^{-2}$. This argues against a reflection-dominated hard X-ray spectrum, which would have implied a much higher $N_{\rm H}$ and luminosity. The X-ray and infrared data are consistent with a bolometric AGN luminosity of $L_{\rm bol}\sim(0.5-2.5)\times10^{47}$erg s$^{-1}$. The X-ray and infrared data are further consistent with coaligned AGN obscurers in which the line of sight "skims" the torus. This is also consistent with the optical spectra, which show both coronal iron lines and broad lines in polarized but not direct light. Combining constraints from the X-ray, optical, and infrared data suggests that the AGN obscurer is within a vertical height of $20$pc, and a radius of $125$pc, of the nucleus.

We present near-infrared integral-field spectroscopic observations targeting H$\alpha$ in eight sub-millimeter galaxies (SMGs) at $z$=1.3-2.5 using VLT/SINFONI, obtaining significant detections for six of them. The star formation rates derived from the H$\alpha$ emission are $\sim$100 M$_\odot$yr$^{-1}$, which account for only $\sim$ 20-30\% of the infrared-derived values, thus suggesting that these systems are very dusty. Two of these systems present [NII]/H$\alpha$ ratios indicative of the presence of an Active Galactic Nucleus (AGN). We mapped the spatial distribution and kinematics of the star forming regions in these galaxies on kpc-scales. In general, the H$\alpha$ morphologies tend to be highly irregular and/or clumpy, showing spatial extents of $\sim$3-11~kpc. We find evidence for significant spatial offsets, of $\sim$0.1-0.4$"$ or 1.2-3.4 kpc, between the H$\alpha$ and the continuum emission in three of the sources. Performing a kinemetry analysis we conclude that the majority of the sample is not consistent with disk-like rotation-dominated kinematics. Instead, they tend to show irregular and/or clumpy and turbulent velocity and velocity dispersion fields. This can be interpreted as evidence for scenario in which these extreme star formation episodes are triggered by galaxy-galaxy interactions and major mergers. In contrast to recent results for SMGs, these sources appear to follow the same relations between gas and star forming rate densities as less luminous and/or normal star forming galaxies.

We present constraints on the molecular outflows in a sample of five Hyper-Luminous Infrared Galaxies using Herschel observations of the OH doublet at 119 {\mu}m. We have detected the OH doublet in three cases: one purely in emission and two purely in absorption. The observed emission profile has a significant blueshifted wing suggesting the possibility of tracing an outflow. Out of the two absorption profiles, one seems to be consistent with the systemic velocity while the other clearly indicates the presence of a molecular outflow whose maximum velocity is about ~1500 km/s. Our analysis shows that this system is in general agreement with previous results on Ultra-luminous Infrared Galaxies and QSOs, whose outflow velocities do not seem to correlate with stellar masses or starburst luminosities (star formation rates). Instead the galaxy outflow likely arises from an embedded AGN.

We present the first survey results of hard X-ray point sources in the Galactic Center (GC) region by NuSTAR. We have discovered 70 hard (3-79 keV) X-ray point sources in a 0.6 deg^2 region around Sgr A* with a total exposure of 1.7 Ms, and 7 sources in the Sgr B2 field with 300 ks. We identify clear Chandra counterparts for 58 NuSTAR sources and assign candidate counterparts for the remaining 19. The NuSTAR survey reaches X-ray luminosities of ~4 x and ~8 x 10^32 erg s^-1 at the GC (8 kpc) in the 3-10 and 10-40 keV bands, respectively. The source list includes three persistent luminous X-ray binaries and the likely run-away pulsar called the Cannonball. New source-detection significance maps reveal a cluster of hard (>10 keV) X-ray sources near the Sgr A diffuse complex with no clear soft X-ray counterparts. The severe extinction observed in the Chandra spectra indicates that all the NuSTAR sources are in the central bulge or are of extragalactic origin. Spectral analysis of relatively bright NuSTAR sources suggests that magnetic cataclysmic variables constitute a large fraction (>40-60%). Both spectral analysis and logN-logS distributions of the NuSTAR sources indicate that the X-ray spectra of the NuSTAR sources should have kT > 20 keV on average for a single temperature thermal plasma model or an average photon index of Gamma = 1.5 - 2 for a power-law model. These findings suggest that the GC X-ray source population may contain a larger fraction of X-ray binaries with high plasma temperatures than the field population.

In this paper we report the results of an X-ray monitoring campaign on the heavily obscured Seyfert galaxy Markarian 3 carried out between the fall of 2014 and the spring of 2015 with NuSTAR, Suzaku and XMM-Newton. The hard X-ray spectrum of Markarian 3 is variable on all the time scales probed by our campaign, down to a few days. The observed continuum variability is due to an intrinsically variable primary continuum seen in transmission through a large, but still Compton-thin column density (N_H~0.8-1.1$\times$10$^{24}$ cm$^{-2}$). If arranged in a spherical-toroidal geometry, the Compton scattering matter has an opening angle ~66 degrees and is seen at a grazing angle through its upper rim (inclination angle ~70 degrees). We report a possible occultation event during the 2014 campaign. If the torus is constituted by a system of clouds sharing the same column density, this event allows us to constrain their number (17$\pm$5) and individual column density, [~(4.9$\pm$1.5)$\times$10$^{22}$ cm$^{-2}$]. The comparison of IR and X-ray spectroscopic results with state-of-the art "torus" models suggests that at least two thirds of the X-ray obscuring gas volume might be located within the dust sublimation radius. We report also the discovery of an ionized absorber, characterised by variable resonant absorption lines due to He- and H-like iron. This discovery lends support to the idea that moderate column density absorbers could be due to clouds evaporated at the outer surface of the torus, possibly accelerated by the radiation pressure due to the central AGN emission leaking through the patchy absorber.

We exploit the widely-separated images of the lensed quasar SDSS J1029+2623 ($z_{em}$=2.197, $\theta =22^{\prime\prime}\!\!.5$) to observe its outflowing wind through two different sightlines. We present an analysis of three observations, including two with the Subaru telescope in 2010 February (Misawa et al. 2013) and 2014 April (Misawa et al. 2014), separated by 4 years, and one with the Very Large Telescope, separated from the second Subaru observation by $\sim$2 months. We detect 66 narrow absorption lines (NALs), of which 24 are classified as intrinsic NALs that are physically associated with the quasar based on partial coverage analysis. The velocities of intrinsic NALs appear to cluster around values of $v_{ej}$ $\sim$ 59,000, 43,000, and 29,000 km/s, which is reminiscent of filamentary structures obtained by numerical simulations. There are no common intrinsic NALs at the same redshift along the two sightlines, implying that the transverse size of the NAL absorbers should be smaller than the sightline distance between two lensed images. In addition to the NALs with large ejection velocities of $v_{ej}$ > 1,000 km/s, we also detect broader proximity absorption lines (PALs) at $z_{abs}$ $\sim$ $z_{em}$. The PALs are likely to arise in outflowing gas at a distance of r $\leq$ 620 pc from the central black hole with an electron density of $n_e$ $\geq$ 8.7$\times$10$^{3}$ cm$^{-3}$. These limits are based on the assumption that the variability of the lines is due to recombination. We discuss the implications of these results on the three-dimensional structure of the outflow.

We present a new metric that uses the spectral curvature (SC) above 10 keV to identify Compton-thick AGN in low-quality Swift BAT X-ray data. Using NuSTAR, we observe nine high SC-selected AGN. We find that high-sensitivity spectra show the majority are Compton-thick (78% or 7/9) and the remaining two are nearly Compton-thick (NH~5-8x10^23 cm^-2). We find the SC_bat and SC_nustar measurements are consistent, suggesting this technique can be applied to future telescopes. We tested the SC method on well-known Compton-thick AGN and find it is much more effective than broad band ratios (e.g. 100% using SC vs. 20% using 8-24/3-8 keV). Our results suggest that using the >10 keV emission may be the only way to identify this population since only two sources show Compton-thick levels of excess in the OIII to X-ray emission ratio (F_OIII/F_2-10 keV>1) and WISE colors do not identify most of them as AGN. Based on this small sample, we find that a higher fraction of these AGN are in the final merger stage than typical BAT AGN. Additionally, these nine obscured AGN have, on average, four times higher accretion rates than other BAT-detected AGN (Edd ratio=0.068 compared to 0.016). The robustness of SC at identifying Compton-thick AGN implies a higher fraction of nearby AGN may be Compton-thick (~22%) and the sum of black hole growth in Compton-thick AGN (Eddington ratio times population percentage), is nearly as large as mildly obscured and unobscured AGN.

We study in detail the evolution of the 2015 outburst of GS 1354-64 (BW Cir) at optical, UV and X-ray wavelengths using Faulkes Telescope South/LCOGT, SMARTS and Swift. The outburst was found to stay in the hard X-ray state, albeit being anomalously luminous with a peak luminosity of L$_{X} >$ 0.15 L$_{Edd}$, which could be the most luminous hard state observed in a black hole X-ray binary. We found that the optical/UV emission is tightly correlated with the X-ray emission, consistent with accretion disc irradiation and/or a jet producing the optical emission. The X-ray spectra can be fitted well with a Comptonisation model, and show softening towards the end of the outburst. In addition, we detect a QPO in the X-ray lightcurves with increasing centroid frequency during the peak and decay periods of the outburst. The long-term optical lightcurves during quiescence show a statistically significant, slow rise of the source brightness over the 7 years prior to the 2015 outburst. This behaviour as well as the outburst evolution at all wavelengths studied can be explained by the disc instability model with irradiation and disc evaporation/condensation.

We present measurements of the evolution of normal-galaxy X-ray emission from $z \approx$ 0-7 using local galaxies and galaxy samples in the 6 Ms Chandra Deep Field-South (CDF-S) survey. The majority of the CDF-S galaxies are observed at rest-frame energies above 2 keV, where the emission is expected to be dominated by X-ray binary (XRB) populations; however, hot gas is expected to provide small contributions to the observed- frame < 1 keV emission at $z < 1$. We show that a single scaling relation between X-ray luminosity ($L_{\rm X}$) and star-formation rate (SFR) is insufficient for characterizing the average X-ray emission at all redshifts. We establish that scaling relations involving not only SFR, but also stellar mass ($M_\star$) and redshift, provide significantly improved characterizations of the average X-ray emission from normal galaxy populations at $z \approx$ 0-7. We further provide the first empirical constraints on the redshift evolution of X-ray emission from both low-mass XRB (LMXB) and high-mass XRB (HMXB) populations and their scalings with $M_\star$ and SFR, respectively. We find $L_{\rm 2-10 keV}({\rm LMXB})/M_\star \propto (1 + z)^{2-3}$ and $L_{\rm 2-10 keV}({\rm HMXB})$/SFR $\propto (1 + z)$, and show that these relations are consistent with XRB population-synthesis model predictions, which attribute the increase in LMXB and HMXB scaling relations with redshift as being due to declining host galaxy stellar ages and metallicities, respectively. We discuss how emission from XRBs could provide an important source of heating to the intergalactic medium in the early Universe, exceeding that of active galactic nuclei.

PG1247+267 is one of the most luminous known quasars at $z\sim2$ and is a strongly super-Eddington accreting SMBH candidate. We obtained NuSTAR data of this intriguing source in December 2014 with the aim of studying its high-energy emission, leveraging the broad band covered by the new NuSTAR and the archival XMM-Newton data. Several measurements are in agreement with the super-Eddington scenario for PG1247+267: the soft power law ($\Gamma=2.3\pm0.1$); the weak ionized Fe emission line and a hint of the presence of outflowing ionized gas surrounding the SMBH. The presence of an extreme reflection component is instead at odds with the high accretion rate proposed for this quasar. This can be explained with three different scenarios; all of them are in good agreement with the existing data, but imply very different conclusions: i) a variable primary power law observed in a low state, superimposed on a reflection component echoing a past, higher flux state; ii) a power law continuum obscured by an ionized, Compton thick, partial covering absorber; and iii) a relativistic disk reflector in a lamp-post geometry, with low coronal height and high BH spin. The first model is able to explain the high reflection component in terms of variability. The second does not require any reflection to reproduce the hard emission, while a rather low high-energy cutoff of $\sim100$ keV is detected for the first time in such a high redshift source. The third model require a face-on geometry, which may affect the SMBH mass and Eddington ratio measurements. Deeper X-ray broad-band data are required in order to distinguish between these possibilities.

We report photometric and spectroscopic observations of the optical transient LSQ13zm. Historical data reveal the presence of an eruptive episode (that we label as `2013a') followed by a much brighter outburst (`2013b') three weeks later, that we argue to be the genuine supernova explosion. This sequence of events closely resemble those observed for SN2010mc and (in 2012) SN2009ip. The absolute magnitude reached by LSQ13zm during 2013a ($M_R=-14.87\pm0.25\,\rm{mag}$) is comparable with those of supernova impostors, while that of the 2013b event ($M_R=-18.46\pm0.21\,\rm{mag}$) is consistent with those of interacting supernovae. Our spectra reveal the presence of a dense and structured circumstellar medium, probably produced through numerous pre-supernova mass-loss events. In addition, we find evidence for high-velocity ejecta, with a fraction of gas expelled at more than 20000\kms. The spectra of LSQ13zm show remarkable similarity with those of well-studied core-collapse supernovae. From the analysis of the available photometric and spectroscopic data, we conclude that we first observed the last event of an eruptive sequence from a massive star, likely a Luminous Blue Variable, which a short time later exploded as a core-collapse supernova. The detailed analysis of archival images suggest that the host galaxy is a star-forming Blue Dwarf Compact Galaxy.

We describe the first data release from the Spitzer-IRAC Equatorial Survey (SpIES); a large-area survey of 115 deg^2 in the Equatorial SDSS Stripe 82 field using Spitzer during its 'warm' mission phase. SpIES was designed to probe sufficient volume to perform measurements of quasar clustering and the luminosity function at z > 3 to test various models for "feedback" from active galactic nuclei (AGN). Additionally, the wide range of available multi-wavelength, multi-epoch ancillary data enables SpIES to identify both high-redshift (z > 5) quasars as well as obscured quasars missed by optical surveys. SpIES achieves 5{\sigma} depths of 6.13 {\mu}Jy (21.93 AB magnitude) and 5.75 {\mu}Jy (22.0 AB magnitude) at 3.6 and 4.5 microns, respectively - depths significantly fainter than WISE. We show that the SpIES survey recovers a much larger fraction of spectroscopically-confirmed quasars (98%) in Stripe 82 than are recovered by WISE (55%). This depth is especially powerful at high-redshift (z > 3.5), where SpIES recovers 94% of confirmed quasars, whereas WISE only recovers 25%. Here we define the SpIES survey parameters and describe the image processing, source extraction, and catalog production methods used to analyze the SpIES data. In addition to this survey paper, we release 234 images created by the SpIES team and three detection catalogs: a 3.6 {\mu}m-only detection catalog containing 6.1 million sources, a 4.5 {\mu}m-only detection catalog containing 6.5 million sources, and a dual-band detection catalog containing 5.4 million sources.

In this paper we present the results of our search for and study of $z \gtrsim 6$ galaxy candidates behind the third Frontier Fields (FF) cluster, MACSJ0717.5+3745, and its parallel field, combining data from Hubble and Spitzer. We select 39 candidates using the Lyman Break technique, for which the clear non-detection in optical make the extreme mid-$z$ interlopers hypothesis unlikely. We also take benefit from $z \gtrsim 6$ samples selected using previous Frontier Fields datasets of Abell 2744 and MACS0416 to improve the constraints on the properties of very high-redshift objects. We compute the redshift and the physical properties, such emission lines properties, star formation rate, reddening, and stellar mass for all Frontier Fields objects from their spectral energy distribution using templates including nebular emission lines. We study the relationship between several physical properties and confirm the trend already observed in previous surveys for evolution of star formation rate with galaxy mass, and between the size and the UV luminosity of our candidates. The analysis of the evolution of the UV Luminosity Function with redshift seems more compatible with an evolution of density. Moreover, no robust $z\ge$8.5 object is selected behind the cluster field, and few $z$$\sim$9 candidates have been selected in the two previous datasets from this legacy survey, suggesting a strong evolution in the number density of galaxies between $z$$\sim$8 and 9. Thanks to the use of the lensing cluster, we study the evolution of the star formation rate density produced by galaxies with L$>$0.03L$^{\star}$, and confirm the strong decrease observed between $z$$\sim$8 and 9.

Heavily obscured accretion is believed to represent an important stage in the growth of supermassive black holes, and to play an important role in shaping the observed spectrum of the Cosmic X-ray Background (CXB). Hard X-ray (E$>$10 keV) selected samples are less affected by absorption than samples selected at lower energies, and are therefore one of the best ways to detect and identify Compton-thick (CT, $\log N_{\rm\,H}\geq 24$) Active Galactic Nuclei (AGN). In this letter we present the first results of the largest broad-band (0.3-150 keV) X-ray spectral study of hard X-ray selected AGN to date, focusing on the properties of heavily obscured sources. Our sample includes the 834 AGN (728 non-blazar, average redshift $z\simeq 0.055$) reported in the 70-months catalog of the all-sky hard X-ray Swift/BAT survey. We find 55 CT AGN, which represent $7.6^{+1.1}_{-2.1}\%$ of our non-blazar sample. Of these, 26 are reported as candidate CT AGN for the first time. We correct for selection bias and derive the intrinsic column density distribution of AGN in the local Universe in two different luminosity ranges. We find a significant decrease in the fraction of obscured Compton-thin AGN for increasing luminosity, from $46\pm3\%$ (for $\log L_{\rm\,14-195} = 40-43.7$) to $39\pm3\%$ (for $\log L_{\rm\,14-195} = 43.7-46$). A similar trend is also found for CT AGN. The intrinsic fraction of CT AGN with $\log N_{\rm\,H}=24-25$ normalised to unity in the $\log N_{\rm H} = 20-25$ range is $27\pm4\%$, and is consistent with the observed value obtained for AGN located within 20 Mpc.

We present improved point-source catalogs for the 2 Ms Chandra Deep Field-North (CDF-N) and the 250 ks Extended Chandra Deep Field-South (E-CDF-S), implementing a number of recent improvements in Chandra source-cataloging methodology. For the CDF-N/E-CDF-S, we provide a main catalog that contains 683/1003 X-ray sources detected with wavdetect at a false-positive probability threshold of $10^{-5}$ that also satisfy a binomial-probability source-selection criterion of $P<0.004$/$P<0.002$. Such an approach maximizes the number of reliable sources detected: a total of 196/275 main-catalog sources are new compared to the Alexander et al. (2003) CDF-N/Lehmer et al. (2005) E-CDF-S main catalogs. We also provide CDF-N/E-CDF-S supplementary catalogs that consist of 72/56 sources detected at the same wavdetect threshold and having $P$ of $0.004-0.1$/$0.002-0.1$ and $K_s\le22.9/K_s\le22.3$ mag counterparts. For all $\approx1800$ CDF-N and E-CDF-S sources, including the $\approx500$ newly detected ones (these being generally fainter and more obscured), we determine X-ray source positions utilizing centroid and matched-filter techniques; we also provide multiwavelength identifications, apparent magnitudes of counterparts, spectroscopic and/or photometric redshifts, basic source classifications, and estimates of observed AGN and galaxy source densities around respective field centers. Simulations show that both the CDF-N and E-CDF-S main catalogs are highly reliable and reasonably complete. Background and sensitivity analyses indicate that the on-axis mean flux limits reached represent a factor of $\approx1.5-2.0$ improvement over the previous CDF-N and E-CDF-S limits. We make our data products publicly available.

Aims. Study the connection between the masing disk and obscuring torus in Seyfert 2 galaxies. Methods. We present a uniform X-ray spectral analysis of the high energy properties of 14 nearby megamaser Active Galactic Nuclei observed by NuSTAR. We use a simple analytical model to localize the maser disk and understand its connection with the torus by combining NuSTAR spectral parameters with available physical quantities from VLBI mapping. Results. Most of the sources analyzed are heavily obscured, showing a column density in excess of $\sim 10^{23}$ cm$^{-2}$. In particular, $79\%$ are Compton-thick ($N_{\rm H} > 1.5 \times 10^{24}$ cm$^{-2}$). Using column densities measured by NuSTAR, with the assumption that the torus is the extension of the maser disk, and further assuming a reasonable density profile, the torus dimensions can be predicted. They are found to be consistent with mid-IR interferometry parsec-scale observations of Circinus and NGC 1068. In this picture, the maser disk is intimately connected to the inner part of the torus. It is probably made of a large number of molecular clouds connecting the torus and the outer part of the accretion disk, giving rise to a thin disk rotating in most cases in Keplerian or sub-Keplerian motion. This toy model explains the established close connection between water megamaser emission and nuclear obscuration as a geometric effect.

Hot Dust-Obscured Galaxies (Hot DOGs) are a population of hyper-luminous infrared galaxies identified by the WISE mission from their very red mid-IR colors, and characterized by hot dust temperatures ($T>60~\rm K$). Several studies have shown clear evidence that the IR emission in these objects is powered by a highly dust-obscured AGN that shows close to Compton-thick absorption at X-ray wavelengths. Thanks to the high AGN obscuration, the host galaxy is easily observable, and has UV/optical colors usually consistent with those of a normal galaxy. Here we discuss a sub-population of 8 Hot DOGs that show enhanced rest-frame UV/optical emission. We discuss three scenarios that might explain the excess UV emission: (i) unobscured light leaked from the AGN by reflection over the dust or by partial coverage of the accretion disk; (ii) a second unobscured AGN in the system; or (iii) a luminous young starburst. X-ray observations can help discriminate between these scenarios. We study in detail the blue excess Hot DOG WISE J020446.13-050640.8, which was serendipitously observed by Chandra/ACIS-I for 174.5 ks. The X-ray spectrum is consistent with a single, hyper-luminous, highly absorbed AGN, and is strongly inconsistent with the presence of a secondary unobscured AGN. Based on this, we argue that the excess blue emission in this object is most likely either due to reflection or a co-eval starburst. We favor the reflection scenario as the unobscured star-formation rate needed to power the UV/optical emission would be $\gtrsim 1000~\rm M_{\odot}~\rm yr^{-1}$. Deep polarimetry observations could confirm the reflection hypothesis.

We present the results of five NuSTAR observations of the type 2 active galactic nucleus (AGN) in IC 751, three of which were performed simultaneously with XMM-Newton or Swift/XRT. We find that the nuclear X-ray source underwent a clear transition from a Compton-thick ($N_{\rm\,H}\simeq 2\times 10^{24}\rm\,cm^{-2}$) to a Compton-thin ($N_{\rm\,H}\simeq 4\times 10^{23}\rm\,cm^{-2}$) state on timescales of $\lesssim 3$ months, which makes IC 751 the first changing-look AGN discovered by NuSTAR. Changes of the line-of-sight column density at a $\sim2\sigma$ level are also found on a time-scale of $\sim 48$ hours ($\Delta N_{\rm\,H}\sim 10^{23}\rm\,cm^{-2}$). From the lack of spectral variability on timescales of $\sim 100$ ks we infer that the varying absorber is located beyond the emission-weighted average radius of the broad-line region, and could therefore be related either to the external part of the broad-line region or a clumpy molecular torus. By adopting a physical torus X-ray spectral model, we are able to disentangle the column density of the non-varying absorber ($N_{\rm\,H}\sim 3.8\times 10^{23}\rm\,cm^{-2}$) from that of the varying clouds [$N_{\rm\,H}\sim(1-150)\times10^{22}\rm\,cm^{-2}$], and to constrain that of the material responsible for the reprocessed X-ray radiation ($N_{\rm\,H} \sim 6 \times 10^{24}\rm\,cm^{-2}$). We find evidence of significant intrinsic X-ray variability, with the flux varying by a factor of five on timescales of a few months in the 2-10 and 10-50 keV band.

We report the detection of a heavily obscured Active Galactic Nucleus (AGN) in the luminous infrared galaxy (LIRG) NGC 6286, identified in a 17.5 ks NuSTAR observation. The source is in an early merging stage, and was targeted as part of our ongoing NuSTAR campaign observing local luminous and ultra-luminous infrared galaxies in different merger stages. NGC 6286 is clearly detected above 10 keV and, by including the quasi-simultaneous Swift/XRT and archival XMM-Newton and Chandra data, we find that the source is heavily obscured [$N_{\rm\,H}\simeq (0.95-1.32)\times 10^{24}\rm\,cm^{-2}$], with a column density consistent with being Compton-thick [CT, $\log (N_{\rm\,H}/\rm cm^{-2})\geq 24$]. The AGN in NGC 6286 has a low absorption-corrected luminosity ($L_{2-10\rm\,keV}\sim 3-20\times 10^{41}\rm\,erg\,s^{-1}$) and contributes $\lesssim$1\% to the energetics of the system. Because of its low-luminosity, previous observations carried out in the soft X-ray band ($<10$ keV) and in the infrared did not notice the presence of a buried AGN. NGC 6286 has multi-wavelength characteristics typical of objects with the same infrared luminosity and in the same merger stage, which might imply that there is a significant population of obscured low-luminosity AGN in LIRGs that can only be detected by sensitive hard X-ray observations.

We assemble a sample of 24 hydrogen-poor super-luminous supernovae (SLSNe). Parameterizing the light curve shape through rise and decline timescales shows that the two are highly correlated. Magnetar-powered models can reproduce the correlation, with the diversity in rise and decline rates driven by the diffusion timescale. Circumstellar interaction models can exhibit a similar rise-decline relation, but only for a narrow range of densities, which may be problematic for these models. We find that SLSNe are approximately 3.5 magnitudes brighter and have light curves 3 times broader than SNe Ibc, but that the intrinsic shapes are similar. There are a number of SLSNe with particularly broad light curves, possibly indicating two progenitor channels, but statistical tests do not cleanly separate two populations. The general spectral evolution is also presented. Velocities measured from Fe II are similar for SLSNe and SNe Ibc, suggesting that diffusion time differences are dominated by mass or opacity. Flat velocity evolution in most SLSNe suggests a dense shell of ejecta. If opacities in SLSNe are similar to other SNe Ibc, the average ejected mass is higher by a factor 2-3. Assuming $\kappa=0.1\,$cm$^2\,$g$^{-1}$, we estimate a mean (median) SLSN ejecta mass of 10$\,$M$_\odot$ (6$\,$M$_\odot$), with a range of 3-30$\,$M$_\odot$. Doubling the assumed opacity brings the masses closer to normal SNe Ibc, but with a high-mass tail. The most probable mechanism for generating SLSNe seems to be the core-collapse of a very massive hydrogen-poor star, forming a millisecond magnetar.

We present the infrared (IR) and X-ray properties of a sample of 33 mid-IR luminous quasars ($\nu$L(6 micron)>6x10$^{44}$ erg/s) at redshift z~1-3, identified through detailed spectral energy distribution analyses of distant star-forming galaxies, using the deepest IR data from Spitzer and Herschel in the GOODS-Herschel fields. The aim is to constrain the fraction of obscured, and Compton-thick (CT, N$_H$>1.5x10$^{24}$ cm$^{-2}$) quasars at the peak era of nuclear and star-formation activities. Despite being very bright in the mid-IR band, ~30% of these quasars are not detected in the extremely deep 2 Ms and 4 Ms Chandra X-ray data available in these fields. X-ray spectral analysis of the detected sources reveals that the majority (~67%) are obscured by column densities N$_H$>10$^{22}$ cm$^{-2}$; this fraction reaches ~80% when including the X-ray undetected sources (9 out of 33), which are likely to be the most heavily-obscured, CT quasars. We constrain the fraction of CT quasars in our sample to be ~24-48%, and their space density to be $\Phi$=(6.7$\pm$2.2)x10$^{-6}$ Mpc$^{-3}$. From the investigation of the quasar host galaxies in terms of star-formation rates (SFRs) and morphological distortions, as a sign of galaxy mergers/interactions, we do not find any direct relation between SFRs and quasar luminosity or X-ray obscuration. On the other hand, there is tentative evidence that the most heavily-obscured quasars have, on average, more disturbed morphologies than the unobscured/moderately-obscured quasar hosts, which preferentially live in undisturbed systems. However, the fraction of quasars with disturbed morphology amongst the whole sample is ~40%, suggesting that galaxy mergers are not the main fuelling mechanism of quasars at z~2.

We present the first results from the KMOS AGN Survey at High redshift (KASHz), a VLT/KMOS integral-field spectroscopic survey of z>0.6 AGN. We present galaxy-integrated spectra of 89 X-ray AGN (Lx=10^42-10^45 erg/s), for which we observed [O III] (z=1.1-1.7) or Halpha emission (z=0.6-1.1). The targets have X-ray luminosities representative of the parent AGN population and we explore the emission-line luminosities as a function of X-ray luminosity. For the [O III] targets, ~50 per cent have ionised gas velocities indicative of gas that is dominated by outflows and/or highly turbulent material (i.e., overall line-widths >~600 km/s). The most luminous half (i.e., Lx>6x10^43 erg/s) have a >~2 times higher incidence of such velocities. On the basis of our results, we find no evidence that X-ray obscured AGN are more likely to host extreme kinematics than unobscured AGN. Our KASHz sample has a distribution of gas velocities that is consistent with a luminosity-matched sample of z<0.4 AGN. This implies little evolution in the prevalence of ionised outflows, for a fixed AGN luminosity, despite an order-of-magnitude decrease in average star-formation rates over this redshift range. Furthermore, we compare our Halpha targets to a redshift-matched sample of star-forming galaxies and despite a similar distribution of Halpha luminosities and likely star-formation rates, we find extreme ionised gas velocities are up to ~10x more prevalent in the AGN-host galaxies. Our results reveal a high prevalence of extreme ionised gas velocities in high-luminosity X-ray AGN and imply that the most powerful ionised outflows in high-redshift galaxies are driven by AGN activity.

To provide the census of the sources contributing to the X-ray background peak above 10 keV, NuSTAR is performing extragalactic surveys using a three-tier "wedding cake" approach. We present the NuSTAR survey of the COSMOS field, the medium sensitivity and medium area tier, covering 1.7 deg2 and overlapping with both Chandra and XMM-Newton data. This survey consists of 121 observations for a total exposure of ~3 Ms. To fully exploit these data, we developed a new detection strategy, carefully tested through extensive simulations. The survey sensitivity at 20% completeness is 5.9, 2.9 and 6.4 x 10^-14 erg/cm2/s in the 3-24 keV, 3-8 keV and 8-24 keV bands, respectively. By combining detections in 3 bands, we have a sample of 91 NuSTAR sources with 10^42 -10^45.5 erg/s luminosities and redshift z=0.04-2.5. Thirty two sources are detected in the 8-24 keV band with fluxes ~100 times fainter than sources detected by Swift-BAT. Of the 91 detections, all but four are associated with a Chandra and/or XMM-Newton point-like counterpart. One source is associated with an extended lower energy X-ray source. We present the X-ray (hardness ratio and luminosity) and optical-to-X-ray properties. The observed fraction of candidate Compton-thick AGN measured from the hardness ratio is between 13-20%. We discuss the spectral properties of NuSTAR J100259+0220.6 (ID 330) at z=0.044, with the highest hardness ratio in the entire sample. The measured column density exceeds 10^24 /cm2, implying the source is Compton-thick. This source was not previously recognized as such without the >10 keV data.

We present the first direct measurements of the rest-frame 10-40 keV X-ray luminosity function (XLF) of Active Galactic Nuclei (AGNs) based on a sample of 94 sources at 0.1 < z <3, selected at 8-24 keV energies from sources in the NuSTAR extragalactic survey program. Our results are consistent with the strong evolution of the AGN population seen in prior, lower-energy studies of the XLF. However, different models of the intrinsic distribution of absorption, which are used to correct for selection biases, give significantly different predictions for the total number of sources in our sample, leading to small, systematic differences in our binned estimates of the XLF. Adopting a model with a lower intrinsic fraction of Compton-thick sources and a larger population of sources with column densities N_H ~ 10^{23-24} /cm2 or a model with stronger Compton reflection component (with a relative normalization of R ~ 2 at all luminosities) can bring extrapolations of the XLF from 2-10 keV into agreement with our NuSTAR sample. Ultimately, X-ray spectral analysis of the NuSTAR sources is required to break this degeneracy between the distribution of absorbing column densities and the strength of the Compton reflection component and thus refine our measurements of the XLF. Furthermore, the models that successfully describe the high-redshift population seen by NuSTAR tend to over-predict previous, high-energy measurements of the local XLF, indicating that there is evolution of the AGN population that is not fully captured by the current models.

We present initial results and the source catalog from the NuSTAR survey of the Extended Chandra Deep Field South (hereafter, ECDFS) - currently the deepest contiguous component of the NuSTAR extragalactic survey program. The survey covers the full ~30 arcmin x 30 arcmin area of this field to a maximum depth of ~360 ks (~220 ks when corrected for vignetting at 3-24 keV), reaching sensitivity limits of ~1.3 x 10^-14 erg/cm2/s (3-8 keV), ~3.4 x 10^-14 erg/cm2/s (8-24 keV) and ~3.0 x 10^-14 erg/cm2/s (3-24 keV). Fifty four (54) sources are detected over the full field, although five of these are found to lie below our significance threshold once contaminating flux from neighboring (i.e., blended) sources is taken into account. Of the remaining 49 that are significant, 19 are detected in the 8-24 keV band. The 8-24 keV to 3-8 keV band ratios of the twelve sources that are detected in both bands span the range 0.39-1.7, corresponding to a photon index range of Gamma ~ 0.5-2.3, with a median photon index of 1.70 +/- 0.52. The redshifts of the 49 sources in our main sample span the range z = 0.21-2.7, and their rest-frame 10-40 keV luminosities (derived from the observed 8-24 keV fluxes) span the range L(10-40 keV) ~ (0.7-300) x 10^43 erg/s, sampling below the "knee" of the X-ray luminosity function out to z ~ 0.8-1. Finally, we identify one NuSTAR source that has neither a Chandra nor an XMM-Newton counterpart, but that shows evidence of nuclear activity at infrared wavelengths, and thus may represent a genuine, new X-ray source detected by NuSTAR in the ECDFS.

We present a NuSTAR and XMM-Newton monitoring campaign in 2014/2015 of the Compton-thick Seyfert 2 galaxy, NGC 1068. During the August 2014 observation, we detect with NuSTAR a flux excess above 20 keV ($32\pm6 \%$) with respect to the December 2012 observation and to a later observation performed in February 2015. We do not detect any spectral variation below 10 keV in the XMM-Newton data. The transient excess can be explained by a temporary decrease of the column density of the obscuring material along the line of sight (from N$_{\rm H}\simeq10^{25}$ cm$^{-2}$ to N$_{\rm H}=6.7\pm1.0\times10^{24}$ cm$^{-2}$), which allows us for the first time to unveil the direct nuclear radiation of the buried AGN in NGC 1068 and to infer an intrinsic 2-10 keV luminosity L$_{\rm X}=7^{+7}_{-4} \times 10^{43}$ erg s$^{-1}$.

NGC 7582 is a well-studied X-ray bright Seyfert 2 with moderately heavy ($N_{\text{H}}\sim10^{23}-10^{24}$~cm$^{-2}$), highly variable absorption and strong reflection spectral features. The spectral shape changed around the year 2000, dropping in observed flux and becoming much more highly absorbed. Two scenarios have been put forth to explain this spectral change: 1) the central X-ray source partially ``shut off'' around this time, decreasing in intrinsic luminosity, with a delayed decrease in reflection features due to the light-crossing time of the Compton-thick material or 2) the source became more heavily obscured, with only a portion of the power law continuum leaking through. NuSTAR observed NGC~7582 twice in 2012, two weeks apart, in order to quantify the reflection using high-quality data above 10 keV. We find that the most plausible scenario is that NGC 7582 has recently become more heavily absorbed by a patchy torus with a covering fraction of $\sim\,80-90\%$ and an equatorial column density of $\sim 3 \times10^{24}$ cm$^{-2}$. We find the need for an additional highly variable full-covering absorber with $N_{\text{H}}= 4-6 \times10^{23}$ cm$^{-2}$ in the line of sight, possibly associated with a hidden broad line region.

Using data from four deep fields (COSMOS, AEGIS, ECDFS, and CDFN), we study the correlation between the position of galaxies in the star formation rate (SFR) versus stellar mass plane and local environment at $z<1.1$. To accurately estimate the galaxy SFR, we use the deepest available Spitzer/MIPS 24 and Herschel/PACS datasets. We distinguish group environments ( $M_{halo}\sim$10$^{12.5-14.2}$$M_{\odot}$) based on the available deep X-ray data and lower halo mass environments based on the local galaxy density. We confirm that the Main Sequence (MS) of star forming galaxies is not a linear relation and there is a flattening towards higher stellar masses ( $M_*>10^{10.4-10.6}$ $M_{\odot}$), across all environments. At high redshift ( $0.5<z<1.1$ ), the MS varies little with environment. At low redshift ( $0.15<z<0.5$ ), group galaxies tend to deviate from the mean MS towards the region of quiescence with respect to isolated galaxies and less-dense environments. We find that the flattening of the MS toward low SFR is due to an increased fraction of bulge dominated galaxies at high masses. Instead, the deviation of group galaxies from the MS at low redshift is caused by a large fraction of red disk dominated galaxies which are not present in the lower density environments. Our results suggest that above a mass threshold ( $\sim10^{10.4}-10^{10.6}$$M_{\odot}$ ) stellar mass, morphology and environment act together in driving the evolution of the SF activity towards lower level. The presence of a dominating bulge and the associated quenching processes are already in place beyond $z\sim$1. The environmental effects appear, instead, at lower redshifts and have a long time-scale.

During the last ~50 years, the population of black hole candidates in X-ray binaries has increased considerably with 59 Galactic objects detected in transient low-mass X-ray binaries, plus a few in persistent systems (including ~5 extragalactic binaries). We collect near-infrared, optical and X-ray information spread over hundreds of references in order to study the population of black holes in X-ray transients as a whole. We present the most updated catalogue of black hole transients, which contains X-ray, optical and near-infrared observations together with their astrometric and dynamical properties. It provides new useful information in both statistical and observational parameters providing a thorough and complete overview of the black hole population in the Milky Way. Analysing the distances and spatial distribution of the observed systems, we estimate a total population of ~1300 Galactic black hole transients. This means that we have already discovered less than ~5% of the total Galactic distribution. The complete version of this catalogue will be continuously updated online and in the Virtual Observatory, including finding charts and data in other wavelengths.

We present VLT VIMOS, Keck DEIMOS and Keck LRIS multi-object spectra of 367 sources in the field of the z ~ 3.09 protocluster SSA22. Sources are spectroscopically classified via template matching, allowing new identifications for 206 extragalactic sources, including 36 z > 2 Lyman-break galaxies (LBGs) and Lyman-\alpha\ emitters (LAEs), 8 protocluster members, and 94 X-ray sources from the ~ 400 ks Chandra deep survey of SSA22. Additionally, in the area covered by our study, we have increased by ~ 4, 13, and 6 times the number of reliable redshifts of sources at 1.0 < z < 2.0, at z > 3.4, and with X-Ray emission, respectively. We compare our results with past spectroscopic surveys of SSA22 to investigate the completeness of the LBGs and the X-Ray properties of the new spectroscopically-classified sources in the SSA22 field.

We searched for z > 7 Lyman-break galaxies (LBGs) in the optical-to-mid-infrared Hubble Frontier Field and associated parallel field observations of the strong-lensing cluster MACS J0416-2403. We discovered 22 candidates, of which six lie at z > 9 and one lies at z > 10. Based on the Hubble and Spitzer photometry, all have secure photometric redshifts and a negligible probability of being at lower redshifts, according to their peak probability ratios, R. This substantial increase in the number of known high-redshift galaxies allows a solid determination of the luminosity function at z > 8. The number of high-z candidates in the parallel field is considerably higher than that in the Abell 2744 parallel field. Our candidates have median stellar masses of log(M_*) ~ 8.40^{+0.55}_{-0.31}~Msun, SFRs of ~ 1.6^{+0.5}_{-0.4} Msun yr^-1, and SFR-weighted ages of < 310^{+70}_{-140} Myr. Finally, we are able to put strong constraints on the z = 7,8,9 and 10 luminosity functions. One of the objects in the cluster field is a z ~ 10 candidate, with a magnification of mu ~ 20 +- 13. This object is likely the faintest z ~ 10 object known to date, allowing a first look into the extreme faint-end (L ~ 0.04L*) of the z ~ 10 luminosity function.

We present a broad-band (~0.3-70 keV) spectral and temporal analysis of NuSTAR observations of the luminous infrared galaxy NGC 6240, combined with archival Chandra, XMM-Newton and BeppoSAX data. NGC 6240 is a galaxy in a relatively early merger state with two distinct nuclei separated by ~1."5. Previous Chandra observations have resolved the two nuclei, showing that they are both active and obscured by Compton-thick material. Although they cannot be resolved by NuSTAR, thanks to the unprecedented quality of the NuSTAR data at energies >10 keV, we clearly detect, for the first time, both the primary and the reflection continuum components. The NuSTAR hard X-ray spectrum is dominated by the primary continuum piercing through an absorbing column density which is mildly optically thick to Compton scattering (tau ~ 1.2, N_H ~ 1.5 x 10^(24) cm^-2). We detect moderate hard X-ray (> 10 keV) flux variability up to 20% on short (15-20 ksec) timescales. The amplitude of the variability is maximum at ~30 keV and is likely to originate from the primary continuum of the southern nucleus. Nevertheless, the mean hard X-ray flux on longer timescales (years) is relatively constant. Moreover, the two nuclei remain Compton-thick, although we find evidence of variability of the material along the line of sight with column densities N_H <~ 2 x 10^(23) cm-2 over long (~3-15 years) timescales. The observed X-ray emission in the NuSTAR energy range is fully consistent with the sum of the best-fit models of the spatially resolved Chandra spectra of the two nuclei.

We present the first sub-arcminute images of the Galactic Center above 10 keV, obtained with NuSTAR. NuSTAR resolves the hard X-ray source IGR J17456-2901 into non-thermal X-ray filaments, molecular clouds, point sources and a previously unknown central component of hard X-ray emission (CHXE). NuSTAR detects four non-thermal X-ray filaments, extending the detection of their power-law spectra with $\Gamma\sim1.3$-$2.3$ up to ~50 keV. A morphological and spectral study of the filaments suggests that their origin may be heterogeneous, where previous studies suggested a common origin in young pulsar wind nebulae (PWNe). NuSTAR detects non-thermal X-ray continuum emission spatially correlated with the 6.4 keV Fe K$\alpha$ fluorescence line emission associated with two Sgr A molecular clouds: MC1 and the Bridge. Broad-band X-ray spectral analysis with a Monte-Carlo based X-ray reflection model self-consistently determined their intrinsic column density ($\sim10^{23}$ cm$^{-2}$), primary X-ray spectra (power-laws with $\Gamma\sim2$) and set a lower limit of the X-ray luminosity of Sgr A* flare illuminating the Sgr A clouds to $L_X \stackrel{>}{\sim} 10^{38}$ erg s$^{-1}$. Above ~20 keV, hard X-ray emission in the central 10 pc region around Sgr A* consists of the candidate PWN G359.95-0.04 and the CHXE, possibly resulting from an unresolved population of massive CVs with white dwarf masses $M_{\rm WD} \sim 0.9 M_{\odot}$. Spectral energy distribution analysis suggests that G359.95-0.04 is likely the hard X-ray counterpart of the ultra-high gamma-ray source HESS J1745-290, strongly favoring a leptonic origin of the GC TeV emission.

We present a NuSTAR, Chandra, and XMM--Newton survey of nine of the nearest ultraluminous infrared galaxies (ULIRGs). The unprecedented sensitivity of NuSTAR at energies above 10 keV enables spectral modeling with far better precision than was previously possible. Six of the nine sources observed were detected sufficiently well by NuSTAR to model in detail their broadband X-ray spectra, and recover the levels of obscuration and intrinsic X-ray luminosities. Only one source (IRAS 13120--5453) has a spectrum consistent with a Compton--thick AGN, but we cannot rule out that a second source (Arp 220) harbors an extremely highly obscured AGN as well. Variability in column density (reduction by a factor of a few compared to older observations) is seen in IRAS 05189--2524 and Mrk 273, altering the classification of these border-line sources from Compton-thick to Compton-thin. The ULIRGs in our sample have surprisingly low observed fluxes in high energy (>10 keV) X-rays, especially compared to their bolometric luminosities. They have lower ratios of unabsorbed 2--10 keV to bolometric luminosity, and unabsorbed 2--10 keV to mid-IR [O IV] line luminosity than do Seyfert 1 galaxies. We identify IRAS 08572+3915 as another candidate intrinsically X-ray weak source, similar to Mrk 231. We speculate that the X-ray weakness of IRAS 08572+3915 is related to its powerful outflow observed at other wavelengths.

In active galactic nuclei (AGN)-galaxy co-evolution models, AGN winds and outflows are often invoked to explain why super-massive black holes and galaxies stop growing efficiently at a certain phase of their lives. They are commonly referred to as the leading actors of feedback processes. Evidence of ultra-fast (v>0.05c) outflows in the innermost regions of AGN has been collected in the past decade by sensitive X-ray observations for sizable samples of AGN, mostly at low redshift. Here we present ultra-deep XMM-Newton and Chandra spectral data of an obscured (Nh~2x10^{23} cm^-2), intrinsically luminous (L2-10keV~4x10^{44} erg/s) quasar (named PID352) at z~1.6 (derived from the X-ray spectral analysis) in the Chandra Deep Field-South. The source is characterized by an iron emission and absorption line complex at observed energies of E~2-3 keV. While the emission line is interpreted as being due to neutral iron (consistent with the presence of cold absorption), the absorption feature is due to highly ionized iron transitions (FeXXV, FeXXVI) with an outflowing velocity of 0.14^{+0.02}_{-0.06}c, as derived from photoionization models. The mass outflow rate - ~2 Msun/yr - is similar to the source accretion rate, and the derived mechanical energy rate is ~9.5x10^{44} erg/s, corresponding to 9% of the source bolometric luminosity. PID352 represents one of the few cases where indications of X-ray outflowing gas have been observed at high redshift thus far. This wind is powerful enough to provide feedback on the host galaxy.

We present two NuSTAR observations of the local Seyfert 2 active galactic nucleus (AGN) and an ultraluminous X-ray source (ULX) candidate in NGC 5643. Together with archival data from Chandra, XMM-Newton and Swift-BAT, we perform a high-quality broadband spectral analysis of the AGN over two decades in energy ($\sim$0.5-100 keV). Previous X-ray observations suggested that the AGN is obscured by a Compton-thick (CT) column of obscuring gas along our line-of-sight. However, the lack of high-quality $\gtrsim$ 10 keV observations, together with the presence of a nearby X-ray luminous source, NGC 5643 X-1, had left significant uncertainties in the characterization of the nuclear spectrum. NuSTAR now enables the AGN and NGC 5643 X-1 to be separately resolved above 10 keV for the first time and allows a direct measurement of the absorbing column density toward the nucleus. The new data show that the nucleus is indeed obscured by a CT column of $N_{\rm{H}}$ $\gtrsim$ 5 $\times$ 10$^{24}$ cm$^{-2}$. The range of 2-10 keV absorption-corrected luminosity inferred from the best fitting models is $L_{2-10,\rm{int}} =$ (0.8-1.7) $\times$ 10$^{42}$ erg s$^{-1}$, consistent with that predicted from multiwavelength intrinsic luminosity indicators. We also study the NuSTAR data for NGC 5643 X-1, and show that it exhibits evidence for a spectral cut-off at energy, $E$ $\sim$ 10 keV, similar to that seen in other ULXs observed by NuSTAR. Along with the evidence for significant X-ray luminosity variations in the 3-8 keV band from 2003-2014, our results further strengthen the ULX classification of NGC 5643 X-1.

We present the characterization and initial results from the QUEST-La Silla AGN variability survey. This is an effort to obtain well sampled optical light curves in extragalactic fields with unique multi-wavelength observations. We present photometry obtained from 2010 to 2012 in the XMM-COSMOS field, which was observed over 150 nights using the QUEST camera on the ESO-Schmidt telescope. The survey uses a broadband filter, the $Q$-band, similar to the union of the $g$ and the $r$ filters, achieving an intrinsic photometric dispersion of $0.05$ mag, and a systematic error of $0.05$ mag in the zero-point. Since some detectors of the camera show significant non-linearity, we use a linear correlation to fit the zero-points as a function of the instrumental magnitudes, thus obtaining a good correction to the non-linear behavior of these detectors. We obtain good photometry to an equivalent limiting magnitude of $r\sim 20.5$. Studying the optical variability of X-ray detected sources in the XMM-COSMOS field, we find that the survey is $\sim75-80$% complete to magnitudes $r\sim20$, and $\sim67$% complete to a magnitude $r\sim21$. The determination and parameterization of the structure function (${SF}_{norm}(\tau) = A \tau^{\gamma}$) of the variable sources shows that most BL AGN are characterized by $A > 0.1$ and $\gamma > 0.025$. It is further shown that variable NL AGN and GAL sources occupying the same parameter space in $A$ and $\gamma$ are very likely to correspond to obscured or low luminosity AGN. Our samples are, however, small, and we expect to revisit these results using larger samples with longer light curves obtained as part of our ongoing survey.

Galactic Center (GC) molecular cloud Sgr B2 is the best manifestation of an X-ray reflection nebula (XRN) reprocessing a past giant outburst from the supermassive black hole Sgr A*. Alternatively, Sgr B2 could be illuminated by low-energy cosmic ray electrons (LECRe) or protons (LECRp). In 2013, NuSTAR for the first time resolved Sgr B2 hard X-ray emission on sub-arcminute scales. Two prominent features are detected above 10 keV - a newly emerging cloud G0.66-0.13 and the central 90" radius region containing two compact cores Sgr B2(M) and Sgr B2(N) surrounded by diffuse emission. It is inconclusive whether the remaining level of Sgr B2 emission is still decreasing or has reached a constant background level. A decreasing Fe K$\alpha$ emission can be best explained by XRN while a constant background emission can be best explained by LECRp. In the XRN scenario, the 3-79 keV Sgr B2 spectrum can well constrain the past Sgr A* outburst, resulting in an outburst spectrum with a peak luminosity of $L_{3-79\rm~keV} \sim 5\times10^{38} \rm~erg~s^{-1}$ derived from the maximum Compton-scattered continuum and the Fe K$\alpha$ emission consistently. The XRN scenario is preferred by the fast variability of G0.66-0.13, which could be a molecular clump located in the Sgr B2 envelope reflecting the same Sgr A* outburst. In the LECRp scenario, we derived the required CR ion power $dW/dt=(1-4)\times10^{39}\rm~erg~s^{-1}$ and the CR ionization rate $\zeta_{H}=(6-10)\times 10^{-15}\rm~H^{-1}~s^{-1}$. The Sgr B2 background level X-ray emission will be a powerful tool to constrain GC CR population.

We present an analysis of the 3-79 keV NuSTAR spectrum of the low-luminosity active galactic nucleus NGC 7213. In agreement with past observations, we find a lower limit to the high-energy cut-off of Ec > 140 keV, no evidence for a Compton-reflected continuum, and the presence of an iron Kalpha complex, possibly produced in the broad-line region. From the application of the MYTorus model, we find that the line-emitting material is consistent with the absence of a significant Compton reflection if arising from a Compton-thin torus of gas with a column density of 5.0(+2.0,-1.6) cm^-2. We report variability of the equivalent width of the iron lines on the time-scale of years using archival observations from XMM-Newton, Chandra and Suzaku. This analysis suggests a possible contribution from dusty gas. A fit with a Comptonization model indicates the presence of a hot corona with a temperature kTe > 40 keV and an optical depth tau < 1, assuming a spherical geometry.

Over the past decades hollow-cathode lamps have been calibration standards for spectroscopic measurements. Advancing to cm/s radial velocity precisions with the next generation of instruments requires more suitable calibration sources with more lines and less dynamic range problems. Fabry-Perot interferometers provide a regular and dense grid of lines and homogeneous amplitudes making them good candidates for next generation calibrators. We investigate the usefulness of Fabry-Perot etalons in wavelength calibration, present an algorithm to incorporate the etalon spectrum in the wavelength solution and examine potential problems. The quasi periodic pattern of Fabry-Perot lines is used along with a hollow-cathode lamp to anchor the numerous spectral features on an absolute scale. We test our method with the HARPS spectrograph and compare our wavelength solution to the one derived from a laser frequency comb. The combined hollow-cathode lamp/etalon calibration overcomes large distortion (50 m/s) in the wavelength solution of the HARPS data reduction software. Direct comparison to the laser frequency comb bears differences of only maximum 10 m/s. Combining hollow-cathode lamps with Fabry-Perot Interferometers can lead to substantial improvements in the wavelength calibration of echelle spectrographs. Etalons can provide economical alternatives to the laser frequency comb, especially for smaller projects.

The intrinsic column density (NH) distribution of quasars is poorly known. At the high obscuration end of the quasar population and for redshifts z<1, the X-ray spectra can only be reliably characterized using broad-band measurements which extend to energies above 10 keV. Using the hard X-ray observatory NuSTAR, along with archival Chandra and XMM-Newton data, we study the broad-band X-ray spectra of nine optically selected (from the SDSS), candidate Compton-thick (NH > 1.5e24 cm^-2) type 2 quasars (CTQSO2s); five new NuSTAR observations are reported herein, and four have been previously published. The candidate CTQSO2s lie at z<0.5, have observed [OIII] luminosities in the range 8.4 < log (L_[OIII]/L_solar) < 9.6, and show evidence for extreme, Compton-thick absorption when indirect absorption diagnostics are considered. Amongst the nine candidate CTQSO2s, five are detected by NuSTAR in the high energy (8-24 keV) band: two are weakly detected at the ~ 3 sigma confidence level and three are strongly detected with sufficient counts for spectral modeling (>~ 90 net source counts at 8-24 keV). For these NuSTAR-detected sources direct (i.e., X-ray spectral) constraints on the intrinsic AGN properties are feasible, and we measure column densities ~2.5-1600 times higher and intrinsic (unabsorbed) X-ray luminosities ~10-70 times higher than pre-NuSTAR constraints from Chandra and XMM-Newton. Assuming the NuSTAR-detected type 2 quasars are representative of other Compton-thick candidates, we make a correction to the NH distribution for optically selected type 2 quasars as measured by Chandra and XMM-Newton for 39 objects. With this approach, we predict a Compton-thick fraction of f_CT = 36^{+14}_{-12} %, although higher fractions (up to 76%) are possible if indirect absorption diagnostics are assumed to be reliable.

We present new NuSTAR and Chandra observations of NGC 3393, a galaxy reported to host the smallest separation dual AGN resolved in the X-rays. While past results suggested a 150 pc separation dual AGN, three times deeper Chandra imaging, combined with adaptive optics and radio imaging suggest a single, heavily obscured, radio-bright AGN. Using VLA and VLBA data, we find an AGN with a two-sided jet rather than a dual AGN and that the hard X-ray, UV, optical, NIR, and radio emission are all from a single point source with a radius <0.2". We find that the previously reported dual AGN is most likely a spurious detection resulting from the low number of X-ray counts (<160) at 6-7 keV and Gaussian smoothing of the data on scales much smaller than the PSF (0.25" vs. 0.80" FWHM). We show that statistical noise in a single Chandra PSF generates spurious dual peaks of the same separation (0.55$\pm$0.07" vs. 0.6") and flux ratio (39$\pm$9% vs. 32% of counts) as the purported dual AGN. With NuSTAR, we measure a Compton-thick source (NH=$2.2\pm0.4\times10^{24}$ cm$^{-2}$) with a large torus half-opening angle, {\theta}=79 which we postulate results from feedback from strong radio jets. This AGN shows a 2-10 keV intrinsic to observed flux ratio of 150. Using simulations, we find that even the deepest Chandra observations would severely underestimate the intrinsic luminosity of NGC 3393 above z>0.2, but would detect an unobscured AGN of this luminosity out to high redshift (z=5).

The Public European Southern Observatory Spectroscopic Survey of Transient Objects (PESSTO) began as a public spectroscopic survey in April 2012. We describe the data reduction strategy and data products which are publicly available through the ESO archive as the Spectroscopic Survey Data Release 1 (SSDR1). PESSTO uses the New Technology Telescope with EFOSC2 and SOFI to provide optical and NIR spectroscopy and imaging. We target supernovae and optical transients brighter than 20.5mag for classification. Science targets are then selected for follow-up based on the PESSTO science goal of extending knowledge of the extremes of the supernova population. The EFOSC2 spectra cover 3345-9995A (at resolutions of 13-18 Angs) and SOFI spectra cover 0.935-2.53 micron (resolutions 23-33 Angs) along with JHK imaging. This data release contains spectra from the first year (April 2012 - 2013), consisting of all 814 EFOSC2 spectra and 95 SOFI spectra (covering 298 distinct objects), in standard ESO Phase 3 format. We estimate the accuracy of the absolute flux calibrations for EFOSC2 to be typically 15%, and the relative flux calibration accuracy to be about 5%. The PESSTO standard NIR reduction process does not yet produce high accuracy absolute spectrophotometry but the SOFI JHK imaging will improve this. Future data releases will focus on improving the automated flux calibration of the data products.

We present X-ray timing and spectral analyses of simultaneous 150 ks Nuclear Spectroscopic Telescope Array (NuSTAR) and Suzaku X-ray observations of the Seyfert 1.5 galaxy NGC 4151. We disentangle the continuum emission, absorption, and reflection properties of the active galactic nucleus (AGN) by applying inner accretion disk reflection and absorption-dominated models. With a time-averaged spectral analysis, we find strong evidence for relativistic reflection from the inner accretion disk. We find that relativistic emission arises from a highly ionized inner accretion disk with a steep emissivity profile, which suggests an intense, compact illuminating source. We find a preliminary, near-maximal black hole spin a>0.9 accounting for statistical and systematic modeling errors. We find a relatively moderate reflection fraction with respect to predictions for the lamp post geometry, in which the illuminating corona is modeled as a point source. Through a time-resolved spectral analysis, we find that modest coronal and inner disk reflection flux variation drives the spectral variability during the observations. We discuss various physical scenarios for the inner disk reflection model, and we find that a compact corona is consistent with the observed features.

Ultra-deep observations of ECDF-S with Chandra and XMM-Newton enable a search for extended X-ray emission down to an unprecedented flux of $2\times10^{-16}$ ergs s$^{-1}$ cm$^{-2}$. We present the search for the extended emission on spatial scales of 32$^{\prime\prime}$ in both Chandra and XMM data, covering 0.3 square degrees and model the extended emission on scales of arcminutes. We present a catalog of 46 spectroscopically identified groups, reaching a redshift of 1.6. We show that the statistical properties of ECDF-S, such as logN-logS and X-ray luminosity function are broadly consistent with LCDM, with the exception that dn/dz/d$\Omega$ test reveals that a redshift range of $0.2<z<0.5$ in ECDF-S is sparsely populated. The lack of nearby structure, however, makes studies of high-redshift groups particularly easier both in X-rays and lensing, due to a lower level of clustered foreground. We present one and two point statistics of the galaxy groups as well as weak-lensing analysis to show that the detected low-luminosity systems are indeed low-mass systems. We verify the applicability of the scaling relations between the X-ray luminosity and the total mass of the group, derived for the COSMOS survey to lower masses and higher redshifts probed by ECDF-S by means of stacked weak lensing and clustering analysis, constraining any possible departures to be within 30% in mass. Abridged.

Superluminous supernovae (SLSNe) are very bright explosions that were only discovered recently and that show a preference for occurring in faint dwarf galaxies. Understanding why stellar evolution yields different types of stellar explosions in these environments is fundamental in order to both uncover the elusive progenitors of SLSNe and to study star formation in dwarf galaxies. In this paper, we present the first results of our project to study SUperluminous Supernova Host galaxIES, focusing on the sample for which we have obtained spectroscopy. We show that SLSNe-I and SLSNe-R (hydrogen-poor) often (~50% in our sample) occur in a class of galaxies that is known as Extreme Emission Line Galaxies (EELGs). The probability of this happening by chance is negligible and we therefore conclude that the extreme environmental conditions and the SLSN phenomenon are related. In contrast, SLSNe-II (hydrogen-rich) occur in more massive, more metal-rich galaxies with softer radiation fields. Therefore, if SLSNe-II constitute a uniform class, their progenitor systems are likely different from those of H-poor SLSNe. Gamma-ray bursts (GRBs) are, on average, not found in as extreme environments as H-poor SLSNe. We propose that H-poor SLSNe result from the very first stars exploding in a starburst, even earlier than GRBs. This might indicate a bottom-light initial mass function in these systems. SLSNe present a novel method of selecting candidate EELGs independent of their luminosity.

The covering factor of Compton-thick obscuring material associated with the torus in active galactic nuclei (AGN) is at present best understood through the fraction of sources exhibiting Compton-thick absorption along the line of sight ($N_{H}>1.5\times10^{24}$ cm$^{-2}$) in the X-ray band, which reveals the average covering factor. Determining this Compton-thick fraction is difficult however, due to the extreme obscuration. With its spectral coverage at hard X-rays ($>$10 keV), NuSTAR is sensitive to the AGN covering factor since Compton scattering of X-rays off optically thick material dominates at these energies. We present a spectral analysis of 10 AGN observed with NuSTAR where the obscuring medium is optically thick to Compton scattering, so called Compton-thick (CT) AGN. We use the torus models of Brightman & Nandra which predict the X-ray spectrum from reprocessing in a torus and include the torus opening angle as a free parameter and aim to determine the covering factor of the Compton-thick gas in these sources individually. Across the sample we find mild to heavy Compton-thick columns, with $N_{H}$ measured from $10^{24}-10^{26}$ cm$^{-2}$, and a wide range of covering factors, where individual measurements range from 0.2-0.9. We find that the covering factor, $f_{c}$, is a strongly decreasing function of the intrinsic 2-10 keV luminosity, $L_{X}$, where $f_{c}=(-0.41\pm0.13)$log$_{10}$($L_{X}$/erg s$^{-1}$)$+18.31\pm5.33$, across more than two orders of magnitude in $L_{X}$ (10$^{41.5}-10^{44}$ erg s$^{-1}$). The covering factors measured here agree well with the obscured fraction as a function of $L_{X}$ as determined by studies of local AGN with $L_{X}>10^{42.5}$ erg s$^{-1}$.

We present photometric and spectroscopic observations of the interacting transient SN 2009ip taken during the 2013 and 2014 observing seasons. We characterise the photometric evolution as a steady and smooth decline in all bands, with a decline rate that is slower than expected for a solely $^{56}$Co-powered supernova at late phases. No further outbursts or eruptions were seen over a two year period from 2012 December until 2014 December. SN 2009ip remains brighter than its historic minimum from pre-discovery images. Spectroscopically, SN 2009ip continues to be dominated by strong, narrow ($\lesssim$2000 km~s$^{-1}$) emission lines of H, He, Ca, and Fe. While we make tenuous detections of [Fe~{\sc ii}] $\lambda$7155 and [O~{\sc i}] $\lambda\lambda$6300,6364 lines at the end of 2013 June and the start of 2013 October respectively, we see no strong broad nebular emission lines that could point to a core-collapse origin. In general, the lines appear relatively symmetric, with the exception of our final spectrum in 2014 May, when we observe the appearance of a redshifted shoulder of emission at +550 km~s$^{-1}$. The lines are not blue-shifted, and we see no significant near- or mid-infrared excess. From the spectroscopic and photometric evolution of SN 2009ip until 820 days after the start of the 2012a event, we still see no conclusive evidence for core-collapse, although whether any such signs could be masked by ongoing interaction is unclear.

We present multi-frequency (1-8 GHz) VLA data, combined with VIMOS IFU data and HST imaging, of a z=0.085 radio-quiet type 2 quasar (with L(1.4GHz)~5e23 W/Hz and L(AGN)~2e45 erg/s). Due to the morphology of its emission-line region, the target (J1430+1339) has been referred to as the Teacup AGN in the literature. We identify "bubbles" of radio emission that are extended ~10-12 kpc to both the east and west of the nucleus. The edge of the brighter eastern bubble is co-spatial with an arc of luminous ionized gas. We also show that the Teacup AGN hosts a compact radio structure, located ~0.8 kpc from the core position, at the base of the eastern bubble. This radio structure is co-spatial with an ionized outflow with an observed velocity of v=-740 km/s. This is likely to correspond to a jet, or possibly a quasar wind, interacting with the interstellar medium at this position. The large-scale radio bubbles appear to be inflated by the central AGN, which indicates that the AGN can also interact with the gas on >~10 kpc scales. Our study highlights that even when a quasar is formally "radio-quiet" the radio emission can be extremely effective for observing the effects of AGN feedback.

Multi-wavelength analysis of the young massive cluster VVV CL077 is presented for the first time. Our Chandra survey of this region enabled the detection of three X-ray emitting stellar members of the cluster, as well as a possible diffuse X-ray component that extends a few arcseconds from the cluster core with an intrinsic flux of (9+/-3)x10^-14 erg cm^-2 s^-1 in the 0.5-10 keV band. Infrared spectra we obtained for two of these X-ray point sources show absorption lines typical of the atmospheres of massive O stars. The X-ray spectrum from the visible extent of VVV CL077 i.e., a 15"-radius around the cluster, can be modeled with an absorbed power law with nH = (6+/-4)x10^22 cm^-2 and gamma = 2+/-1. In addition, the X-ray core of VVV CL077 coincides with diffuse emission seen in the infrared band and with a local maximum in the radio continuum map. A possible association with a neighboring H II region would place VVV CL077 at a distance of around 11 kpc; on the far side of the Norma Arm. At this distance, the cluster is 0.8 pc wide with a mass density of (1-4)x10^3 Msol pc^-3.

(Abridged) The HST Frontier Fields project started at the end of 2013 with the aim of providing extremely deep images of 6 massive galaxy clusters. One of the main goals of this program is to push several telescopes to their limits in order to provide the best current view of the earliest stages of the Universe. We present a detailed analysis of $z$$\sim$8 objects behind the HFFs lensing cluster, MACS0416-2403, combining 0.3-1.6 $\mu$m imaging from HST, ground-based $K_s$ imaging from VLT HAWK-I, and 3.6 $\mu$m and 4.5 $\mu$m Spitzer Space Telescope. We apply the classical Lyman Break technique, combining non-detection criteria and color-selection. In order to avoid contamination by mid-$z$ interlopers, we require a strong break between optical and near infrared data. The luminosity function at $z$$\sim$8 is computed using a MC method taking advantage of the SED-fitting results. A piece of cautionary information is gleaned from new deep optical photometry of a previously identified $z$$\sim$8 galaxy in this cluster, which is now firmly detected as a mid-$z$ interloper with a strong $\approx1.5$ mag Balmer break (between F606W and F125W). Using the SED of this interloper, we estimate the contamination rate of our MACS0416$-$2403 sample, as well as for previous samples in Abell 2744 based on HFF data. Our selection recovers 4 robust objects with m$_{F160W}$ ranging from 26.0 to 27.9 AB and located in modest amplification regions ($\mu<$2.4). Two of the objects display a secondary break between the IRAC 3.6 $\mu$m and 4.5 $\mu$m bands which could be associated to the Balmer break or emission lines at $z$$\sim$8. The candidates generally have star formation rates around $\sim$10 M$_{\odot}$/yr and sizes ranging from 0.2 to 0.5 kpc, in good agreement with previous observations and expectations for objects in the early Universe.

We present NuSTAR observations of the bright Seyfert 2 galaxy NGC 2110 obtained in 2012, when the source was at the highest flux level ever observed, and in 2013, when the source was at a more typical flux level. We include archival observations from other X-ray satellites, namely XMM-Newton, Suzaku, BeppoSAX, Chandra and Swift. Simultaneous NuSTAR and Swift broad band spectra (in the 3-80 keV range) indicate a cutoff energy $E_{\rm c}>210$ keV, with no detectable contribution from Compton reflection. NGC 2110 is one of the very few sources where no evidence for distant Compton thick scattering is found and, by using temporal information collected over more than a decade, we investigate variations of the iron K$\alpha$ line on time scales of years. The Fe K$\alpha$ line is likely the sum of two components: one constant (originating from distant Compton-thick material) and the other one variable and linearly correlated with the source flux (possibly arising from Compton-thin material much closer to the black hole).

We report on observations of NGC1068 with NuSTAR, which provide the best constraints to date on its $>10$~keV spectral shape. We find no strong variability over the past two decades, consistent with its Compton-thick AGN classification. The combined NuSTAR, Chandra, XMM-Newton, and Swift-BAT spectral dataset offers new insights into the complex reflected emission. The critical combination of the high signal-to-noise NuSTAR data and a spatial decomposition with Chandra allow us to break several model degeneracies and greatly aid physical interpretation. When modeled as a monolithic (i.e., a single N_H) reflector, none of the common Compton-reflection models are able to match the neutral fluorescence lines and broad spectral shape of the Compton reflection. A multi-component reflector with three distinct column densities (e.g., N_H~1.5e23, 5e24, and 1e25 cm^{-2}) provides a more reasonable fit to the spectral lines and Compton hump, with near-solar Fe abundances. In this model, the higher N_H components provide the bulk of the Compton hump flux while the lower N_H component produces much of the line emission, effectively decoupling two key features of Compton reflection. We note that ~30% of the neutral Fe Kalpha line flux arises from >2" (~140 pc), implying that a significant fraction of the <10 keV reflected component arises from regions well outside of a parsec-scale torus. These results likely have ramifications for the interpretation of poorer signal-to-noise observations and/or more distant objects [Abridged].

We derive photometric redshifts (\zp) for sources in the entire ($\sim0.4$ deg$^2$) Hawaii-Hubble Deep Field-North (\hdfn) field with the EAzY code, based on point spread function-matched photometry of 15 broad bands from the ultraviolet (\bandu~band) to mid-infrared (IRAC 4.5 $\mu$m). Our catalog consists of a total of 131,678 sources. We evaluate the \zp~quality by comparing \zp~with spectroscopic redshifts (\zs) when available, and find a value of normalized median absolute deviation \sigm$=$0.029 and an outlier fraction of 5.5\% (outliers are defined as sources having $\rm |\zp - \zs|/(1+\zs) > 0.15$) for non-X-ray sources. More specifically, we obtain \sigm$=0.024$ with 2.7\% outliers for sources brighter than $R=23$~mag, \sigm$=0.035$ with 7.4\% outliers for sources fainter than $R=23$~mag, \sigm$=$0.026 with 3.9\% outliers for sources having $z<1$, and \sigm$=$0.034 with 9.0\% outliers for sources having $z>1$. Our \zp\ quality shows an overall improvement over an earlier \zp\ work that focused only on the central \hdfn\ area. We also classify each object as star or galaxy through template spectral energy distribution fitting and complementary morphological parametrization, resulting in 4959 stars and 126,719 galaxies. Furthermore, we match our catalog with the 2~Ms {\it Chandra} Deep Field-North main \xray~catalog. For the 462 matched non-stellar \xray~sources (281 having \zs), we improve their \zp~quality by adding three additional AGN templates, achieving \sigm$=0.035$ and an outlier fraction of 12.5\%. We make our catalog publicly available presenting both photometry and \zp, and provide guidance on how to make use of our catalog.

We present a catalog of 1415 X-ray sources identified in the Norma arm region Chandra survey (NARCS), which covers a 2 deg x 0.8 deg region in the direction of the Norma spiral arm to a depth of $\approx$20 ks. Of these sources, 1130 are point-like sources detected with $\geq3\sigma$ confidence in at least one of three energy bands (0.5-10, 0.5-2, and 2-10 keV), five have extended emission, and the remainder are detected at low significance. Since most sources have too few counts to permit individual classification, they are divided into five spectral groups defined by their quantile properties. We analyze stacked spectra of X-ray sources within each group, in conjunction with their fluxes, variability, and infrared counterparts, to identify the dominant populations in our survey. We find that $\sim$50% of our sources are foreground sources located within 1-2 kpc, which is consistent with expectations from previous surveys. Approximately 20% of sources are likely located in the proximity of the Scutum-Crux and near Norma arm, while 30% are more distant, in the proximity of the far Norma arm or beyond. We argue that a mixture of magnetic and nonmagnetic CVs dominates the Scutum-Crux and near Norma arms, while intermediate polars (IPs) and high-mass stars (isolated or in binaries) dominate the far Norma arm. We also present the cumulative number count distribution for sources in our survey that are detected in the hard energy band. A population of very hard sources in the vicinity of the far Norma arm and active galactic nuclei dominate the hard X-ray emission down to $f_X\approx10^{-14}$ erg cm$^{-2}$ s$^{-1}$, but the distribution curve flattens at fainter fluxes. We find good agreement between the observed distribution and predictions based on other surveys.

We study the evolution of the total star formation (SF) activity, total stellar mass and halo occupation distribution in massive halos by using one of the largest X-ray selected sample of galaxy groups with secure spectroscopic identification in the major blank field surveys (ECDFS, CDFN, COSMOS, AEGIS). We provide an accurate measurement of SFR for the bulk of the star-forming galaxies using very deep mid-infrared Spitzer MIPS and far-infrared Herschel PACS observations. For undetected IR sources, we provide a well-calibrated SFR from SED fitting. We observe a clear evolution in the level of SF activity in galaxy groups. The total SF activity in the high redshift groups (0.5<z<1.1) is higher with respect to the low redshift (0.15<z<0.5) sample at any mass by 0.8+/-0.12 dex. A milder difference (0.35+/-0.1 dex) is observed between the low redshift bin and the groups at z~0. We show that the level of SF activity is declining more rapidly in the more massive halos than in the more common lower mass halos. We do not observe any evolution in the halo occupation distribution and total stellar mass- halo mass relations in groups. The picture emerging from our findings suggests that the galaxy population in the most massive systems is evolving faster than galaxies in lower mass halos, consistently with a "halo downsizing" scenario.

We study the relationship between the structure and star-formation rate (SFR) of X-ray selected low and moderate luminosity active galactic nuclei (AGNs) in the two Chandra Deep Fields, using Hubble Space Telescope imaging from the Cosmic Assembly Near Infrared Extragalactic Legacy Survey (CANDELS) and deep far-infrared maps from the PEP+GOODS-Herschel survey. We derive detailed distributions of structural parameters and FIR luminosities from carefully constructed control samples of galaxies, which we then compare to those of the AGNs. At z~1, AGNs show slightly diskier light profiles than massive inactive (non-AGN) galaxies, as well as modestly higher levels of gross galaxy disturbance (as measured by visual signatures of interactions and clumpy structure). In contrast, at z~2, AGNs show similar levels of galaxy disturbance as inactive galaxies, but display a red central light enhancement, which may arise due to a more pronounced bulge in AGN hosts or due to extinguished nuclear light. We undertake a number of tests of these alternatives, but our results do not strongly favour one interpretation over the other. The mean SFR and its distribution among AGNs and inactive galaxies are similar at z>1.5. At z<1, however, clear and significant enhancements are seen in the SFRs of AGNs with bulge-dominated light profiles. These trends suggest an evolution in the relation between nuclear activity and host properties with redshift, towards a minor role for mergers and interactions at z>1.5.

The deflection angles of lensed sources increase with their distance behind a given lens. We utilize this geometric effect to corroborate the $z_{phot}\simeq9.8$ photometric redshift estimate of a faint near-IR dropout, triply-imaged by the massive galaxy cluster Abell 2744 in deep Hubble Frontier Fields images. The multiple images of this source follow the same symmetry as other nearby sets of multiple images which bracket the critical curves and have well defined redshifts (up to $z_{spec}\simeq3.6$), but with larger deflection angles, indicating that this source must lie at a higher redshift. Similarly, our different parametric and non-parametric lens models all require this object be at $z\gtrsim4$, with at least 95\% confidence, thoroughly excluding the possibility of lower-redshift interlopers. To study the properties of this source we correct the two brighter images for their magnifications, leading to a SFR of $\sim0.3 M_{\odot}$/yr, a stellar mass of $\sim4\times10^{7} M_{\odot}$, and an age of $\lesssim220$ Myr (95\% confidence). The intrinsic apparent magnitude is 29.9 AB (F160W), and the rest-frame UV ($\sim1500 \AA$) absolute magnitude is $M_{UV,AB}=-17.6$. This corresponds to $\sim0.1 L^{*}_{z=8}$ ($\sim0.2 L^{*}_{z=10}$, adopting $dM^{*}/dz\sim0.45$), making this candidate one of the least luminous galaxies discovered at $z\sim10$.

We present X-ray spectral analyses for three Seyfert 2 active galactic nuclei, NGC 424, NGC 1320, and IC 2560, observed by NuSTAR in the 3-79 keV band. The high quality hard X-ray spectra allow detailed modeling of the Compton reflection component for the first time in these sources. Using quasi-simultaneous NuSTAR and Swift/XRT data, as well as archival XMM-Newton data, we find that all three nuclei are obscured by Compton-thick material with column densities in excess of ~5 x $10^{24}$ cm$^{-2}$, and that their X-ray spectra above 3 keV are dominated by reflection of the intrinsic continuum on Compton-thick material. Due to the very high obscuration, absorbed intrinsic continuum components are not formally required by the data in any of the sources. We constrain the intrinsic photon indices and the column density of the reflecting medium through the shape of the reflection spectra. Using archival multi-wavelength data we recover the intrinsic X-ray luminosities consistent with the broadband spectral energy distributions. Our results are consistent with the reflecting medium being an edge-on clumpy torus with a relatively large global covering factor and overall reflection efficiency of the order of 1%. Given the unambiguous confirmation of the Compton-thick nature of the sources, we investigate whether similar sources are likely to be missed by commonly used selection criteria for Compton-thick AGN, and explore the possibility of finding their high-redshift counterparts.

We report NuSTAR observations of a sample of six X-ray weak broad absorption line (BAL) quasars. These targets, at z=0.148-1.223, are among the optically brightest and most luminous BAL quasars known at z<1.3. However, their rest-frame 2 keV luminosities are 14 to >330 times weaker than expected for typical quasars. Our results from a pilot NuSTAR study of two low-redshift BAL quasars, a Chandra stacking analysis of a sample of high-redshift BAL quasars, and a NuSTAR spectral analysis of the local BAL quasar Mrk 231 have already suggested the existence of intrinsically X-ray weak BAL quasars, i.e., quasars not emitting X-rays at the level expected from their optical/UV emission. The aim of the current program is to extend the search for such extraordinary objects. Three of the six new targets are weakly detected by NuSTAR with <45 counts in the 3-24 keV band, and the other three are not detected. The hard X-ray (8-24 keV) weakness observed by NuSTAR requires Compton-thick absorption if these objects have nominal underlying X-ray emission. However, a soft stacked effective photon index ({\Gamma}~1.8) for this sample disfavors Compton-thick absorption in general. The uniform hard X-ray weakness observed by NuSTAR for this and the pilot samples selected with <10 keV weakness also suggests that the X-ray weakness is intrinsic in at least some of the targets. We conclude that the NuSTAR observations have likely discovered a significant population (>33%) of intrinsically X-ray weak objects among the BAL quasars with significantly weak <10 keV emission. We suggest that intrinsically X-ray weak quasars might be preferentially observed as BAL quasars.

We report the discovery of 24 Lyman-break candidates at 7<z<10.5, in the Hubble Frontier Fields (HFF) imaging data of Abell 2744 (z=0.308), plus Spiter/IRAC data and archival ACS data. The sample includes a triple image system with a photometric redshift of z~7.4. This high redshift is geometrically confirmed by our lens model corresponding to deflection angles that are 12% larger than the lower-redshift systems used to calibrate the lens model at z=2.019. The majority of our high-redshift candidates are not expected to be multiply lensed given their locations in the image plane and the brightness of foreground galaxies, but are magnified by factors of ~1.3-15, so that we are seeing further down the luminosity function than comparable deep field imaging. It is apparent that the redshift distribution of these sources does not smoothly extend over the full redshift range accessible at z<12, but appears to break above z=9. Nine candidates are clustered within a small region of 20" across, representing a potentially unprecedented concentration. Given the poor statistics, however, we must await similar constraints from the additional HFF clusters to properly examine this trend. The physical properties of our candidates are examined using the range of lens models developed for the HFF program by various groups including our own, for a better estimate of underlying systematics. Our spectral-energy-distribution fits for the brightest objects suggest stellar masses of ~1E9 solar mass, star-formation rates of ~4 solar mass per yr, and a typical formation redshift of z<19.