We examine the thermodynamic state and cooling of the low-$z$ Circum-Galactic Medium (CGM) in five FIRE-2 galaxy formation simulations of Milky Way-mass galaxies. We find that the CGM in these simulations is generally multiphase and dynamic, with a wide spectrum of largely nonlinear density perturbations sourced by the accretion of gas from the IGM and winds from both the central and satellite galaxies. We investigate the origin of the multiphase structure of the CGM with a particle tracking analysis and find that most of the low entropy gas has cooled from the hot halo as a result of thermal instability triggered by these perturbations. The ratio of cooling to free-fall timescales $t_{\rm cool}/t_{\rm ff}$ in the hot component of the CGM spans a wide range $\sim 1-100$ at a given radius, but exhibits approximately constant median values $\sim 5-20$ at all radii $0.1 R_{\rm vir} < r < R_{\rm vir}$. These are similar to the $\approx 10-20$ value typically adopted as the thermal instability threshold in "precipitation" models of the ICM. Consequently, a one-dimensional model based on the assumption of a constant $t_{\rm cool}/t_{\rm ff}$ and hydrostatic equilibrium approximately reproduces the simulation number density and entropy profiles, but only if it assumes the metallicity profiles taken directly from the simulations. We explicitly show that the $t_{\rm cool}/t_{\rm ff}$ value of a gas parcel in the hot component of the CGM does not predict its probability of cooling and subsequently accreting onto the central galaxy. This suggests that the value of $t_{\rm cool}/t_{\rm ff}$ is a poor predictor of thermal stability in gaseous halos in which large-amplitude density perturbations are prevalent.

We use the FIRE-2 cosmological simulations to study the formation of a virial temperature, quasi-static gas phase in the circumgalactic medium (CGM) at redshifts 0<z<5, and how the formation of this virialized phase affects the evolution of galactic discs. We demonstrate that when the halo mass crosses ~10^12 M_sun, the cooling time of shocked gas in the inner CGM (~0.1 R_vir, where R_vir is the virial radius) exceeds the local free-fall time. The inner CGM then experiences a transition from on average sub-virial temperatures (T<<T_vir), large pressure fluctuations and supersonic inflow/outflow velocities, to virial temperatures (T~T_vir), uniform pressures and subsonic velocities. This transition occurs when the outer CGM (~0.5 R_vir) is already subsonic and has a temperature ~T_vir, indicating that the longer cooling times at large radii allow the outer CGM to virialize at lower halo masses than the inner CGM. This outside-in CGM virialization scenario is in contrast with inside-out scenarios commonly envisioned based on more idealized simulations. We demonstrate that virialization of the inner CGM coincides with abrupt changes in the properties of the central galaxy and its stellar feedback: the galaxy settles into a stable rotating disc, star formation transitions from `bursty' to `steady,' and stellar-driven galaxy-scale outflows are suppressed. Our results hence suggest that CGM virialization is initially associated with the formation of rotation-dominated thin galactic discs, rather than with the quenching of star formation as often assumed.

Pressure balance plays a central role in models of the interstellar medium (ISM), but whether and how pressure balance is realized in a realistic multiphase ISM is not yet well understood. We address this question using a set of FIRE-2 cosmological zoom-in simulations of Milky Way-mass disk galaxies, in which a multiphase ISM is self-consistently shaped by gravity, cooling, and stellar feedback. We analyze how gravity determines the vertical pressure profile as well as how the total ISM pressure is partitioned between different phases and components (thermal, dispersion/turbulence, and bulk flows). We show that, on average and consistent with previous more idealized simulations, the total ISM pressure balances the weight of the overlying gas. Deviations from vertical pressure balance increase with increasing galactocentric radius and with decreasing averaging scale. The different phases are in rough total pressure equilibrium with one another, but with large deviations from thermal pressure equilibrium owing to kinetic support in the cold and warm phases, which dominate the total pressure near the midplane. Bulk flows (e.g., inflows and fountains) are important at a few disk scale heights, while thermal pressure from hot gas dominates at larger heights. Overall, the total midplane pressure is well-predicted by the weight of the disk gas, and we show that it also scales linearly with the star formation rate surface density (Sigma_SFR). These results support the notion that the Kennicutt-Schmidt relation arises because Sigma_SFR and the gas surface density (Sigma_g) are connected via the ISM midplane pressure.

We revisit the question of 'hot mode' versus 'cold mode' accretion onto galaxies using steady-state cooling flow solutions and idealized 3D hydrodynamic simulations. We demonstrate that for the hot accretion mode to exist, the cooling time is required to be longer than the free-fall time near the radius where the gas is rotationally-supported, R_circ, i.e. the existence of the hot mode depends on physical conditions at the galaxy scale rather than on physical conditions at the halo scale. When allowing for the depletion of the halo baryon fraction relative to the cosmic mean, the longer cooling times imply that a virialized gaseous halo may form in halo masses below the threshold of ~10^12 M_sun derived for baryon-complete halos. We show that for any halo mass there is a maximum accretion rate for which the gas is virialized throughout the halo and can accrete via the hot mode of Mdot_crit ~ 0.7(v_c/100 km/s)^5.4 (R_circ / 10 kpc) (Z/Z_sun)^-0.9 M_sun / yr, where Z and v_c are the metallicity and circular velocity measured at R_circ. For accretion rates >~Mdot_crit the volume-filling gas phase can in principle be `transonic' -- virialized in the outer halo but cool and free-falling near the galaxy. We compare Mdot_crit to the average star formation rate (SFR) in halos at 0<z<10 implied by the stellar-mass halo-mass relation. For a plausible metallicity evolution with redshift, we find that SFR <~ Mdot_crit at most masses and redshifts, suggesting that the SFR of galaxies could be primarily sustained by the hot mode in halo masses well below the classic threshold of ~10^12 M_sun.

Galactic-scale winds driven by active galactic nuclei (AGN) are often invoked to suppress star formation in galaxy evolution models, but the mechanisms driving these outflows are hotly debated. Two key AGN feedback models are (1) radiation pressure accelerating cool gas and (2) a hot outflowing wind entraining the ISM. Highly ionized emission-line diagnostics represent a powerful means of differentiating these scenarios because of their sensitivity to the expected compression of the ISM clouds by the hot wind. Here, we report the first spatially resolved UV emission spectroscopy of a prototypical (radio-quiet) quasar-driven superwind around the obscured quasar SDSSJ1356+1026 at z=0.123. We observe ratios of OVI/CIV, NV/CIV, and CIV/HeII that are remarkably similar for outflowing gas clouds <100 pc and ~10 kpc from the nucleus. Such similarity is expected for clouds with AGN radiation pressure dominated dynamics. Comparing the observed line emission to models of clouds in balance with radiation pressure and/or a hot wind, we rule out the presence of a dynamically important hot wind and constrain the ratio of hot gas pressure to radiation pressure to P_hot/P_rad<0.25 both at <100 pc and ~10 kpc from the nucleus. Moreover, the predictions of the radiation pressure confined cloud models that best fit observed UV line ratios are consistent with the observed diffuse X-ray spectrum. These results indicate that this AGN superwind is driven by radiation pressure or was driven by a hot wind that has since dissipated despite on-going AGN activity.

We analyze the different fates of the circumgalactic medium (CGM) in FIRE-2 cosmological simulations, focusing on the redshifts z=0.25 and z=2 representative of recent surveys. Our analysis includes 21 zoom-in simulations covering the halo mass range Mh(z=0) ~ 10^10 - 10^12 Msun. We analyze both where the gas ends up after first leaving the CGM (its "proximate" fate), as well as its location at z=0 (its "ultimate" fate). Of the CGM at z=2, about half is found in the ISM or stars of the central galaxy by z=0 in Mh(z=2) ~ 5e11 Msun halos, but most of the CGM in lower-mass halos ends up in the IGM. This is so even though most of the CGM in M_h(z=2) ~ 5e10 Msun halos first accretes onto the central galaxy before being ejected into the IGM. On the other hand, most of the CGM mass at z=0.25 remains in the CGM by z=0 at all halo masses analyzed. Of the CGM gas that subsequently accretes onto the central galaxy in the progenitors of Mh(z=0) ~10^12 Msun halos, most of it is cool (T~10^4 K) at z=2 but hot (~Tvir) at z=0.25, consistent with the expected transition from cold mode to hot mode accretion. Despite the transition in accretion mode, at both z=0.25 and z=2 >~80% of the cool gas in Mh >~ 10^11 Msun halos will accrete onto a galaxy. We find that the metallicity of CGM gas is typically a poor predictor of both its proximate and ultimate fates. This is because there is in general little correlation between the origin of CGM gas and its fate owing to substantial mixing while in the CGM.

We investigate the impact of cosmic rays (CRs) on the circumgalactic medium (CGM) in FIRE-2 simulations, for ultra-faint dwarf through Milky Way (MW)-mass halos hosting star-forming (SF) galaxies. Our CR treatment includes injection by supernovae, anisotropic streaming and diffusion along magnetic field lines, collisional and streaming losses, with constant parallel diffusivity $\kappa\sim3\times10^{29}\,\mathrm{cm^2\ s^{-1}}$ chosen to match $\gamma$-ray observations. With this, CRs become more important at larger halo masses and lower redshifts, and dominate the pressure in the CGM in MW-mass halos at $z\lesssim 1-2$. The gas in these ``CR-dominated'' halos differs significantly from runs without CRs: the gas is primarily cool (a few $\sim10^{4}\,$K), and the cool phase is volume-filling and has a thermal pressure below that needed for virial or local thermal pressure balance. Ionization of the ``low'' and ``mid'' ions in this diffuse cool gas is dominated by photo-ionization, with O VI columns $\gtrsim 10^{14.5}\,\mathrm{cm^{-2}}$ at distances $\gtrsim 150\,\mathrm{kpc}$. CR and thermal gas pressure are locally anti-correlated, maintaining total pressure balance, and the CGM gas density profile is determined by the balance of CR pressure gradients and gravity. Neglecting CRs, the same halos are primarily warm/hot ($T\gtrsim 10^{5}\,$K) with thermal pressure balancing gravity, collisional ionization dominates, O VI columns are lower and Ne VIII higher, and the cool phase is confined to dense filaments in local thermal pressure equilibrium with the hot phase.

In several models of galaxy formation feedback occurs in cycles or mainly at high redshift. At times and in regions where feedback heating is ineffective, hot gas in the galaxy halo is expected to form a cooling flow, where the gas advects inward on a cooling timescale. Cooling flow solutions can thus be used as a benchmark for observations and simulations to constrain the timing and extent of feedback heating. Using analytic calculations and idealized 3D hydrodynamic simulations, we show that for a given halo mass and cooling function, steady-state cooling flows form a single-parameter family of solutions, while initially hydrostatic gaseous halos converge on one of these solutions within a cooling time. The solution is thus fully determined once either the mass inflow rate $\dot{M}$ or the total halo gas mass are known. In the Milky Way (MW) halo, a cooling flow with $\dot{M}$ equal to the star formation rate predicts a ratio of the cooling time to the free-fall time of ~10, similar to some feedback-regulated models. This solution also correctly predicts observed OVII and OVIII absorption columns, and the gas density profile implied by OVII and OVIII emission. These results suggest ongoing heating by feedback may be negligible in the inner MW halo. Extending similar solutions out to the cooling radius however underpredicts observed OVI columns around the MW and around other low-redshift star-forming galaxies. This can be reconciled with the successes of the cooling flow model with either a mechanism which preferentially heats the OVI-bearing outer halo, or alternatively if OVI traces cool photoionized gas beyond the accretion shock. We also demonstrate that the entropy profiles of some of the most relaxed clusters are reasonably well described by a cooling flow solution.

We use a particle tracking analysis to study the origins of the circumgalactic medium (CGM), separating it into (1) accretion from the intergalactic medium (IGM), (2) wind from the central galaxy, and (3) gas ejected from other galaxies. Our sample consists of 21 FIRE-2 simulations, spanning the halo mass range log(Mh/Msun) ~ 10-12 , and we focus on z=0.25 and z=2. Owing to strong stellar feedback, only ~L* halos retain a baryon mass >~50% of their cosmic budget. Metals are more efficiently retained by halos, with a retention fraction >~50%. Across all masses and redshifts analyzed >~60% of the CGM mass originates as IGM accretion (some of which is associated with infalling halos). Overall, the second most important contribution is wind from the central galaxy, though gas ejected or stripped from satellites can contribute a comparable mass in ~L* halos. Gas can persist in the CGM for billions of years, resulting in well-mixed halo gas. Sight lines through the CGM are therefore likely to intersect gas of multiple origins. For low-redshift ~L* halos, cool gas (T<10^4.7 K) is distributed on average preferentially along the galaxy plane, however with strong halo-to-halo variability. The metallicity of IGM accretion is systematically lower than the metallicity of winds (typically by >~1 dex), although CGM and IGM metallicities depend significantly on the treatment of subgrid metal diffusion. Our results highlight the multiple physical mechanisms that contribute to the CGM and will inform observational efforts to develop a cohesive picture.

The past decade has seen an explosion of discoveries and new insights into the diffuse gas within galaxies, galaxy clusters, and the filaments composing the Cosmic Web. A new decade will bring fresh opportunities to further this progress towards developing a comprehensive view of the composition, thermal state, and physical processes of diffuse gas in the Universe. Ultraviolet (UV) spectroscopy, probing diffuse 10^4-10^6 K gas at high spectral resolution, is uniquely poised to (1) witness environmental galaxy quenching processes in action, such as strangulation and tidal- and ram-pressure stripping, (2) directly account for the baryon content of galaxy clusters in the cold-warm (T<10^6 K) gas, (3) determine the phase structure and kinematics of gas participating in the equilibrium-regulating exchange of energy at the cores of galaxy clusters, and (4) map cold streams and filaments of the Cosmic Web that feed galaxies and clusters. With a substantial UV undertaking beyond the Hubble Space Telescope, all of the above would be achievable over the entire epoch of galaxy cluster formation. Such capabilities, coupled with already-planned advancements at other wavelengths, will transform extragalactic astronomy by revealing the dominant formation and growth mechanisms of gaseous halos over the mass spectrum, settling the debate between early- and late-time metal enrichment scenarios, and revealing how the ecosystems in which galaxies reside ultimately facilitate their demise.

Galaxies evolve under the influence of gas flows between their interstellar medium and their surrounding gaseous halos known as the circumgalactic medium (CGM). The CGM is a major reservoir of galactic baryons and metals, and plays a key role in the long cycles of accretion, feedback, and recycling of gas that drive star formation. In order to fully understand the physical processes at work within galaxies, it is therefore essential to have a firm understanding of the composition, structure, kinematics, thermodynamics, and evolution of the CGM. In this white paper we outline connections between the CGM and galactic star formation histories, internal kinematics, chemical evolution, quenching, satellite evolution, dark matter halo occupation, and the reionization of the larger-scale intergalactic medium in light of the advances that will be made on these topics in the 2020s. We argue that, in the next decade, fundamental progress on all of these major issues depends critically on improved empirical characterization and theoretical understanding of the CGM. In particular, we discuss how future advances in spatially-resolved CGM observations at high spectral resolution, broader characterization of the CGM across galaxy mass and redshift, and expected breakthroughs in cosmological hydrodynamic simulations will help resolve these major problems in galaxy evolution.

The observed spatial and kinematic overlap between soft X-ray emission and the Narrow Line Region (NLR) in obscured Active Galactic Nuclei (AGN) yields compelling evidence that relatively low-density gas co-exists with higher density gas on scales as large as 100s of pc. This is commonly interpreted as evidence for a constant gas pressure multiphase medium, likely produced by thermal instability. Alternatively, Radiation Pressure Compression (RPC) also leads to a density distribution, since a gas pressure (and hence density) gradient must arise within each cloud to counteract the incident ionising radiation pressure. RPC leads to a well-defined ionization distribution, and a Differential Emission Measure (DEM) distribution with a universal slope of $\sim-0.9$, weakly dependent on the gas properties and the illuminating radiation field. In contrast, a multiphase medium does not predict the form of the DEM. The observed DEMs of obscured AGN with XMM-Newton RGS spectra (the CHRESOS sample) are in striking agreement with the predicted RPC DEM, providing a clear signature that RPC is the dominant mechanism for the observed range of densities in the X-ray NLR. In contrast with the constant gas pressure multiphase medium, RPC further predicts an increasing gas pressure with decreasing ionization, which can be tested with future X-ray missions using density diagnostics.

We study the components of cool and warm/hot gas in the circumgalactic medium (CGM) of simulated galaxies and address the relative production of OVI by photoionization versus collisional ionization, as a function of halo mass, redshift, and distance from the galaxy halo center. This is done utilizing two different suites of zoom-in hydro-cosmological simulations, VELA (6 halos; $z>1$) and NIHAO (18 halos; to $z=0$), which provide a broad theoretical basis because they use different codes and physical recipes for star formation and feedback. In all halos studied in this work, we find that collisional ionization by thermal electrons dominates at high redshift, while photoionization of cool or warm gas by the metagalactic radiation takes over near $z\sim2$. In halos of $\sim 10^{12}M_{\odot}$ and above, collisions become important again at $z<0.5$, while photoionization remains significant down to $z=0$ for less massive halos. In halos with $M_{\textrm v}>3\times10^{11}~M_{\odot}$, at $z\sim 0$ most of the photoionized OVI is in a warm, not cool, gas phase ($T\lesssim 3\times 10^5$~K). We also find that collisions are dominant in the central regions of halos, while photoionization is more significant at the outskirts, around $R_{\textrm v}$, even in massive halos. This too may be explained by the presence of warm gas or, in lower mass halos, by cool gas inflows.

Large OVI columns are observed around star-forming, low-redshift ~L* galaxies, with a dependence on impact parameter indicating that most O^5+ particles reside beyond half the halo virial radius (>~100 kpc). In order to constrain the nature of the gas traced by OVI, we analyze additional observables of the outer halo, namely HI to OVI column ratios of 1-10, an absence of low-ion absorption, a mean differential extinction of E(B-V)~10^-3, and a linear relation between OVI column and velocity width. We contrast these observations with two physical scenarios: (1) OVI traces high-pressure (~30 cm^-3 K) collisionally-ionized gas cooling from a virially-shocked phase, and (2) OVI traces low-pressure (<~1 cm^-3 K) gas beyond the accretion shock, where the gas is in ionization and thermal equilibrium with the UV background. We demonstrate that the high-pressure scenario requires multiple gas phases to explain the observations, and a large deposition of energy at >~100 kpc to offset the energy radiated by the cooling gas. In contrast, the low-pressure scenario can explain all considered observations with a single gas phase in thermal equilibrium, provided that the baryon overdensity is comparable to the dark-matter overdensity, and that the gas is enriched to >~Z_sun/3 with an ISM-like dust-to-metal ratio. The low-pressure scenario implies that OVI traces a cool flow with mass flow rate of ~5 M_sun yr^-1, comparable to the star formation rate of the central galaxies. The OVI line widths are consistent with the velocity shear expected within this flow. The low-pressure scenario predicts a bimodality in absorption line ratios at ~100 kpc, due to the pressure jump across the accretion shock.

We present a new method to quantify the value of the escape fraction of ionizing photons, and the existence of ultra-faint galaxies clustered around brighter objects during the epoch of cosmic reionization, using the diffuse Ly$\alpha$, continuum and H$\alpha$ emission observed around galaxies at $z\sim6$. We model the surface brightness profiles of the diffuse halos considering the fluorescent emission powered by ionizing photons escaping from the central galaxies, and the nebular emission from satellite star-forming sources, by extending the formalisms developed in Mas-Ribas & Dijkstra (2016) and Mas-Ribas et al. (2017). The comparison between our predicted profiles and Ly$\alpha$ observations at $z=5.7$ and $z=6.6$ favors a low ionizing escape fraction, $f_{\rm esc}^{\rm ion}\sim5\%$, for galaxies in the range $-19\gtrsim M_{\rm UV} \gtrsim -21.5$. However, uncertainties and possible systematics in the observations do not allow for firm conclusions. We predict H$\alpha$ and rest-frame visible continuum observations with JWST, and show that JWST will be able to detect extended (a few tens of kpc) fluorescent H$\alpha$ emission powered by ionizing photons escaping from a bright, $L\gtrsim 5L^*$, galaxy. Such observations can differentiate fluorescent emission from nebular emission by satellite sources. We discuss how observations and stacking of several objects may provide unique constraints on the escape fraction for faint galaxies and/or the abundance of ultra-faint radiation sources.

The typical optical-UV continuum slopes observed in many type 1 AGN are redder than expected from thin accretion disk models. A possible resolution to this conundrum is that many AGN are reddened by dust along the line of sight. To explore this possibility, we stack 5000 SDSS AGN with luminosity L~10^45erg/s and redshift z~0.4 in bins of optical continuum slope alpha_opt and width of the broad H$\beta$ emission line. We measure the EW of the NaID absorption feature in each stacked spectrum. We find a linear relation between alpha_opt and EW(NaID), such that EW(NaID) increases as alpha_opt becomes redder. In the bin with the smallest H$\beta$ width, objects with the bluest slopes that are similar to accretion disk predictions are found to have EW(NaID)=0, supporting the line-of-sight dust hypothesis. This conclusion is also supported by the dependence of the $H\alpha/H\beta$ line ratio on alpha_opt. The implied relationship between alpha_opt and dust reddening is given by E(B-V)~0.2(-0.1-alpha_opt), and the implied reddening of a typical type 1 AGN with alpha_opt=-0.5 is E(B-V)~0.08mag. Photoionization calculations show that the dusty gas responsible for reddening is too ionized to produce the observed features. Therefore, we argue that the sodium absorption arises in regions of the host ISM which are shielded from the AGN radiation, and the correlation with alpha_opt arises since ISM columns along shielded and non-shielded sightlines are correlated. This scenario is supported by the similarity of the E(B-V)-NaID relation implied by our results with the relation in the Milky-Way found by previous studies.

We develop a new method to constrain the physical conditions in the cool (~10^4 K) circumgalactic medium (CGM) from measurements of ionic column densities, by assuming that the cool CGM spans a large range of gas densities and that small high-density clouds are hierarchically embedded in large low-density clouds. The new method combines the information available from different sightlines during the photoionization modeling, thus yielding tighter constraints on CGM properties compared to traditional methods which model each sightline individually. Applying this new technique to the COS-Halos survey of low-redshift ~L* galaxies, we find that we can reproduce all observed ion columns in all 44 galaxies in the sample, from the low-ions to OVI, with a single universal density structure for the cool CGM. The gas densities span the range 50 < \rho/\rho_mean < 5x10^5 (\rho_mean is the cosmic mean), while the physical size of individual clouds scales as ~\rho^-1, from ~35 kpc of the low density OVI clouds to ~6 pc of the highest density low-ion clouds. The deduced cloud sizes are too small for this density structure to be driven by self-gravity, thus its physical origin is unclear. The implied cool CGM mass within the virial radius is 1.3x10^10 M_sun (~1% of the halo mass), distributed rather uniformly over the four decades in density. The mean cool gas density profile scales as R^-1.0, where R is the distance from the galaxy center. We construct a 3D model of the cool CGM based on our results, which we argue provides a benchmark for the CGM structure in hydrodynamic simulations. Our results can be tested by measuring the coherence scales of different ions.

Quasar feedback models often predict an expanding hot gas bubble which drives a galaxy-scale outflow. In many circumstances this hot gas radiates inefficiently and is therefore difficult to observe directly. We present an indirect method to detect the presence of a hot bubble using hydrostatic photoionization calculations of the cold (~10^4 K) line-emitting gas. We compare our calculations with observations of the broad line region, the inner face of the torus, the narrow line region (NLR), and the extended NLR, and thus constrain the hot gas pressure at distances 0.1 pc -- 10 kpc from the center. We find that emission line ratios observed in the average quasar spectrum are consistent with radiation-pressure-dominated models on all scales. On scales <40 pc a dynamically significant hot gas pressure is ruled out, while on larger scales the hot gas pressure cannot exceed six times the local radiation pressure. In individual quasars, ~25% of quasars exhibit NLR ratios that are inconsistent with radiation-pressure-dominated models, though in these objects the hot gas pressure is also unlikely to exceed the radiation pressure by an order of magnitude or more. The derived upper limits on the hot gas pressure imply that the instantaneous gas pressure force acting on galaxy-scale outflows falls short of the time-averaged force needed to explain the large momentum fluxes \dot{p} >> L_AGN/c inferred for galaxy-scale outflows. This apparent discrepancy can be reconciled if optical quasars previously experienced a buried, fully-obscured phase during which the hot gas bubble was more effectively confined and during which galactic wind acceleration occurred.

We study the physical and kinematic properties of the narrow line region (NLR) of the nearest obscured quasar MRK 477 (z=0.037), using optical and near-infrared spectroscopy. We explore a diversity of aspects that provide a more complete understanding of the nature of this object, example of a type 2 quasar in the nearby Universe, as well as a starburst-AGN hybrid system [abridged].

The ionising continuum from active galactic nuclei (AGN) is fundamental for interpreting their broad emission lines and understanding their impact on the surrounding gas. Furthermore, it provides hints on how matter accretes onto supermassive black holes. Using HST's Wide Field Camera 3 we have constructed the first stacked ultraviolet (rest-frame wavelengths 600-2500\AA) spectrum of 53 luminous quasars at z=2.4, with a state-of-the-art correction for the intervening Lyman forest and Lyman continuum absorption. The continuum slope ($f_\nu \propto \nu^{\alpha_\nu}$) of the full sample shows a break at ~912\AA\ with spectral index $\alpha_\nu=-0.61\pm0.01$ at $\lambda>912$\AA\ and a softening at shorter wavelengths ($\alpha_\nu=-1.70 \pm 0.61$ at $\lambda\leq 912$\AA). Our analysis proves that a proper intergalactic medium absorption correction is required to establish the intrinsic continuum emission of quasars. We interpret our average ultraviolet spectrum in the context of photoionisation, accretion disk models, and quasar contribution to the ultraviolet background. We find that observed broad line ratios are consistent with those predicted assuming an ionising slope of $\alpha_\mathrm{ion}=$-2.0, similar to the observed ionising spectrum in the same wavelength range. The continuum break and softening are consistent with accretion disk plus X-ray corona models when black hole spin is taken into account. Our spectral energy distribution yields a 30% increase to previous estimates of the specific quasar emissivity, such that quasars may contribute significantly to the total specific Lyman limit emissivity estimated from the Ly$\alpha$ forest at z<3.2.

The broad line region (BLR) of luminous active galactic nuclei (AGN) is a prominent observational signature of the accretion flow around supermassive black holes, which can be used to measure their masses (M_BH) over cosmic history. Due to the <100 {\mu}as angular size of the BLR, current direct constraints on BLR kinematics are limited to those provided by reverberation mapping studies, which are most efficiently carried out on low-luminosity L and low-redshift z AGN. We analyze the possibility to measure the BLR size and study its kinematic structure using spectroastrometry, whereby one measures the spatial position centroid of emission line photons as a function of velocity. We calculate the expected spectroastrometric signal of a rotation-dominated BLR for various assumptions about the ratio of random to rotational motions, and the radial distribution of the BLR gas. We show that for hyper-luminous quasars at z < 2.5, the size of the low-ionization BLR can already be constrained with existing telescopes and adaptive optics systems, thus providing a novel method to spatially resolve the kinematics of the accretion flow at 10^3 -- 10^4 gravitational radii, and measure M_BH at the high-L end of the AGN family. With a 30m-class telescope, BLR spectroastrometry should be routinely detectable for much fainter quasars out to z ~ 6, and for various emission lines. This will enable kinematic M_BH measurements as a function of luminosity and redshift, providing a compelling science case for next generation telescopes.

The winds of ionized gas driven by Active Galactic Nuclei (AGN) can be studied through absorption lines in their X-ray spectra. A recurring feature of these outflows is their broad ionization distribution, including essentially all ionization levels (e.g., Fe^0+ to Fe^25+). This characteristic feature can be quantified with the absorption measure distribution (AMD), defined as the distribution of column density with ionization parameter |dN / dlog xi|. Observed AMDs extend over 0.1 < xi < 10^4 (cgs), and are remarkably similar in different objects. Power-law fits (|dN /dlog xi| ~ N_1 xi^a) yield N_1 = 3x10^{21} cm^-2 +- 0.4 dex and a = 0 -- 0.4. What is the source of this broad ionization distribution, and what sets the small range of observed $N_1$ and $a$? A common interpretation is a multiphase outflow, with a wide range of gas densities in a uniform gas pressure medium. However, the incident radiation pressure leads to a gas pressure gradient in the photoionized gas, and therefore to a broad range of ionization states within a single slab. We show that this compression of the gas by the radiation pressure leads to an AMD with |dN / dlog xi| = 8x10^21 xi^0.03 cm^-2, remarkably similar to that observed. The calculated values of $N_1$ and $a$ depend weakly on the gas metallicity, the ionizing spectral slope, the distance from the nucleus, the ambient density, and the total absorber column. Thus, radiation pressure compression (RPC) of the photoionized gas provides a natural explanation for the observed AMD. RPC predicts that the gas pressure increases with decreasing ionization, which can be used to test the validity of RPC in ionized AGN outflows.

A fraction of quasars present broad absorption lines, produced by outflowing gas with typical velocities of 3000 - 10,000 km/s. If the outflowing gas fills a significant fraction of the volume where it resides, then it will be highly ionized by the quasar due to its low density, and will not produce the observed UV absorption. The suggestion that the outflow is shielded from the ionizing radiation was excluded by recent observations. The remaining solution is a dense outflow with a filling factor $f<10^{-3}$. What produces such a small $f$? Here we point out that radiation pressure confinement (RPC) inevitably leads to gas compression and the formation of dense thin gas sheets/filaments, with a large gradient in density and ionization along the line of sight. The total column of ionized dustless gas is a few times $10^{22}$ cm$^{-2}$, consistent with the observed X-ray absorption and detectable P V absorption. The predicted maximal columns of various ions show a small dependence on the system parameters, and can be used to test the validity of RPC as a solution for the overionization problem. The ionization structure of the outflow implies that if the outflow is radiatively driven, then broad absorption line quasars should have $L/L_{\rm Edd} \gtrsim 0.1$.

We report on balanced SIS receivers covering the astronomical important 180-720 GHz submillimeter atmospheric window. To facilitate remote observations and automated spectral line surveys, fully synthesized local oscillators are employed. High-current-density Nb-AlN-Nb superconducting-insulating-superconducting (SIS) tunnel junctions are used as the mixing element. The measured double-sideband (DSB) 230 GHz receiver noise temperature, uncorrected for optics loss, ranges from 50K at 185 GHz, 33K at 246 GHz, to 51K at 280 GHz. In this frequency range the mixer has a DSB conversion gain of 0 +- 1.5 dB. The measured 460 GHz double-sideband receiver noise temperature, uncorrected for optics loss, is 32K at 400 GHz, 34K at 460 GHz, and 61K at 520 GHz. Similar to the 230 GHz balanced mixer, the DSB mixer conversion gain is 1 +- 1 dB. To help optimize performance, the mixer IF circuits and bias injection are entirely planar by design. Dual-frequency observation, by means of separating the incoming circular polarized electric field into two orthogonal components, is another important mode of operation offered by the new facility instrumentation. Instrumental stability is excellent supporting the LO noise cancellation properties of the balanced mixer configuration. In the spring of 2012 the dual-frequency 230/460 SIS receiver was successfully installed at Caltech Submillimeter Observatory (CSO), Mauna Kea, HI.

Active galactic nuclei (AGN) are characterized by similar broad emission lines properties at all luminosities ($10^{39}$-$10^{47}$ erg s$^{-1}$). What produces this similarity over a vast range of $10^8$ in luminosity? Photoionization is inevitably associated with momentum transfer to the photoionized gas. Yet, most of the photoionized gas in the Broad Line Region (BLR) follows Keplerian orbits, which suggests that the BLR originates from gas with a large enough column for gravity to dominate. The photoionized surface layer of the gas must develop a pressure gradient due to the incident radiation force. We present solutions for the structure of such a hydrostatic photoionized gas layer in the BLR. The gas is stratified, with a low-density highly-ionized surface layer, a density rise inwards, and a uniform-density cooler inner region, where the gas pressure reaches the incident radiation pressure. This radiation pressure confinement (RPC) of the photoionized layer leads to a universal ionization parameter $U\sim 0.1$ in the inner photoionized layer, independent of luminosity and distance. Thus, RPC appears to explain the universality of the BLR properties in AGN. We present predictions for the BLR emission per unit covering factor, as a function of distance from the ionizing source, for a range of ionizing continuum slopes and gas metallicity. The predicted mean strength of most lines (excluding H$\beta$), and their different average-emission radii, are consistent with the available observations.

We analyze the hydrostatic effect of AGN radiation pressure on optically thick gas in the host galaxy. We show that in luminous AGN, the radiation pressure likely confines the ionized layer of the illuminated gas. Radiation pressure confinement (RPC) has two main implications. First, the gas density near the ionization front is 7x10^4 L_{i,45} r_{50}^{-2} cm^{-3}, where L_{i,45} is the ionizing luminosity in units of 10^45 erg/s and r_{50} is the distance of the gas from the nucleus in units of 50 pc. Second, as shown by Dopita et al., the solution of the ionization structure within each slab is unique, independent of the ambient pressure. We show that the RPC density vs. distance relation is observed over a dynamical range of ~10^4 in distance, from sub-pc to kpc from the nucleus, and a range of ~10^8 in gas density, from 10^3 to 10^11 cm^{-3}. This relation implies that the radiative force of luminous AGN can compress giant molecular clouds in the host galaxy, and possibly affect the star formation rate. The unique ionization structure in RPC includes a highly ionized X-ray emitting surface, an intermediate layer which emits coronal lines, and a lower ionization inner layer which emits optical lines. This structure can explain the observed overlap of the extended X-ray and optical narrow line emission in nearby AGN. We further support RPC by comparing the predicted ratios of the narrow lines strength and narrow line widths with available observations. We suggest a new method, based on the narrow line widths, to estimate the black hole mass of low luminosity AGN.

We performed a 3D radiative MHD simulation of the solar atmosphere. This simulation shows a jet-like feature that shows similarities to the type II spicules observed for the first time with Hinode. Rapid Blueshifted Events (RBEs) on the solar disk are associated with these spicules. Observational results suggest they may contribute significantly in supplying the corona with hot plasma. We perform a detailed comparison of the properties of the simulated jet with those of type II spicules (observed with Hinode) and RBEs (with ground-based instruments). We analyze variety of synthetic emission and absorption lines from the simulations including chromospheric Ca II and Ha to TR and coronal temperatures (10E4 to several 10E6K). We compare their synthetic intensities, line profiles, Doppler shifts, line widths and asymmetries with observations from Hinode/SOT and EIS, SOHO/SUMER, SST and SDO/AIA. Many properties of the synthetic observables resemble the observations, and we describe in detail the physical processes that lead to these observables. Detailed analysis of the synthetic observables provides insight into how observations should be analyzed to derive information about physical variables in such a dynamic event. For example, we find that LOS superposition in the optically thin atmosphere requires the combination of Doppler shifts and spectral line asymmetry to determine the velocity in the jet. Other properties differ from the observations, especially in the chromospheric lines. The mass density of the part of the spicule with a chromospheric temperature is too low to produce significant opacity in chromospheric lines. These and other discrepancies are described in detail, and we discuss which mechanisms and physical processes may need to be included in the MHD simulations to mimic the thermodynamic processes of the chromosphere and corona, in particular to reproduce type II spicules.

We present the optical narrow line ratios in an SDSS based sample of 3,175 broad Ha selected type 1 AGN, and explore their positions in the BPT diagrams as a function of the AGN and the host properties. We find the following: 1. The luminosities of all measured narrow lines (Ha, Hb, [OIII], [NII], [SII], [OI]) show a Baldwin relation relative to the broad Ha luminosity L_bHa, with slopes in the range of 0.53-0.72. 2. About 20% of the type 1 AGN reside within the `Composite' and `SF' regions of the BPT diagrams. These objects also show excess narrow Ha and UV luminosities, for their L_bHa, consistent with contribution from star formation which dominates the narrow lines emission, as expected from their positions in the BPT diagrams. 3. The type 1 which reside within the AGN region in the BPT diagrams, are offset to lower [SII]/Ha and [NII]/Ha luminosity ratios, compared to type 2 AGN. This offset is a selection effect, related to the lower AGN/host luminosity selection of the type 2 AGN selected from the SDSS galaxy sample. 4. The [NII]/Ha and [NII]/[SII] ratios in type 1 AGN increase with the host mass, as expected if the mass-metallicity relation of quiescent galaxies holds for the AGN narrow line region. 5. The broad lines optical FeII is higher for a higher [NII]/Ha, at a fixed L_Bol and Eddington ratio L/L_Edd. This suggests that the broad line region metallicity is also related to the host mass. 6. The fraction of AGN which are LINERs increases sharply with decreasing L/L_Edd. This fraction is the same for type 1 and type 2 AGN. 7. The BPT position is unaffected by the amount of dust extinction of the optical-UV continuum, which suggests the extincting dust resides on scales larger than the NLR.

We explore the relative strength of the narrow emission lines in an SDSS based sample of broad H-alpha selected AGN, defined in paper I. We find a decrease in the narrow to broad H-alpha luminosity (L_bHa) ratio with increasing L_bHa, such that both L([OIII] lambda5007) and L(narrow H-alpha) scale as L_bHa^0.7 for 10^40 < L_bHa < 10^45 ergs s^-1. Following our earlier result that L_bHa \propto L_bol, this trend indicates that the relative narrow line luminosity decreases with increasing L_bol. We derive L_bol / 10^43 ergs s^-1 = 4000 (L([OIII]) / 10^43 ergs s^-1)^1.39. This implies that the bolometric correction factor, L_bol / L([OIII]), decreases from 3,000 at L_bol = 10^46.1 ergs s^-1 to 300 at L_bol = 10^42.5 ergs s^-1. At low luminosity, the narrow component dominates the observed H-alpha profile, and most type 1 AGN appear as intermediate type AGN. Partial obscuration or extinction cannot explain the dominance of intermediate type AGN at low luminosity, and the most likely mechanism is a decrease in the narrow line region covering factor with increasing L_bol. Deviations from the above trend occur in objects with L / L_Edd <~ 10^-2.6, probably due to the transition to LINERs with suppressed [OIII] emission, and in objects with M_BH > 10^8.5 M_Sun, probably due to the dominance of radio loud AGN, and associated enhanced [OIII] emission.

We present the emission properties of a sample of 3,579 type 1 AGN, selected based on the detection of broad H-alpha emission. The sample covers the range of black hole mass 10^6<M_BH/M_Sun<10^9.5 and luminosity in Eddington units 10^-3 < L/L_Edd < 1. Our main results are: 1. The distribution of the H-alpha FWHM values is independent of luminosity. 2. The observed mean optical-UV SED is well matched by a fixed shape SED of luminous quasars, which scales linearly with broad H-alpha luminosity, and a host galaxy contribution. 3. The host galaxy r-band (fibre) luminosity function follows well the luminosity function of inactive non-emission line galaxies (NEG), consistent with a fixed fraction of ~3% of NEG hosting an AGN, regardless of the host luminosity. 4. The optical-UV SED of the more luminous AGN shows a small dispersion, consistent with dust reddening of a blue SED, as expected for thermal thin accretion disc emission. 5. There is a rather tight relation of nuL_nu(2 keV) and broad H-alpha luminosity, which provides a useful probe for unobscured (true) type 2 AGN.

We analyze the emission properties of a new sample of 3,579 type 1 AGN, selected from the SDSS DR7 based on the detection of broad H-alpha emission. The sample extends over a broad H-alpha luminosity L_bHa of 10^40 - 10^44 erg s^-1 and a broad H-alpha FWHM of 1,000 - 25,000 km s^-1, which covers the range of black hole mass 10^6<M_BH/M_Sun<10^9.5 and luminosity in Eddington units 10^-3 < L/L_Edd < 1. We combine ROSAT, GALEX and 2MASS observations to form the SED from 2.2 mic to 2 keV. We find the following: 1. The distribution of the H-alpha FWHM values is independent of luminosity. 2. The observed mean optical-UV SED is well matched by a fixed shape SED of luminous quasars, which scales linearly with L_bHa, and a host galaxy contribution. 3. The host galaxy r-band (fibre) luminosity function follows well the luminosity function of inactive non-emission line galaxies (NEG), consistent with a fixed fraction of ~3% of NEG hosting an AGN, regardless of the host luminosity. 4. The hosts of lower luminosity AGN have a mean z band luminosity and u-z colour which are identical to NEG with the same redshift distribution. With increasing L_bHa the AGN hosts become bluer and less luminous than NEG. The implied increasing star formation rate with L_bHa is consistent with the relation for SDSS type 2 AGN of similar bolometric luminosity. 5. The optical-UV SED of the more luminous AGN shows a small dispersion, consistent with dust reddening of a blue SED, as expected for thermal thin accretion disc emission. 6. There is a rather tight relation of nuL_nu(2 keV) and L_bHa, which provides a useful probe for unobscured (true) type 2 AGN. 7. The primary parameter which drives the X-ray to UV emission ratio is the luminosity, rather than M_BH or L/L_Edd.

We report on an initial analysis of Herschel/HIFI observations of hydrogen chloride (HCl), hydrogen peroxide (H_2O_2), and molecular oxygen (O_2) in the Martian atmosphere performed on 13 and 16 April 2010 (L_s ~ 77{\deg}). We derived a constant volume mixing ratio of 1400 +/- 120 ppm for O_2 and determined upper limits of 200 ppt for HCl and 2 ppb for H_2O_2. Radiative transfer model calculations indicate that the vertical profile of O_2 may not be constant. Photochemical models determine the lowest values of H_2O_2 to be around L_s ~ 75{\deg} but overestimate the volume mixing ratio compared to our measurements.

In the framework of the Water in Star-forming regions with Herschel (WISH) key program, maps in water lines of several outflows from young stars are being obtained, to study the water production in shocks and its role in the outflow cooling. This paper reports the first results of this program, presenting a PACS map of the o-H2O 179 um transition obtained toward the young outflow L1157. The 179 um map is compared with those of other important shock tracers, and with previous single-pointing ISO, SWAS, and Odin water observations of the same source that allow us to constrain the water abundance and total cooling. Strong H2O peaks are localized on both shocked emission knots and the central source position. The H2O 179 um emission is spatially correlated with emission from H2 rotational lines, excited in shocks leading to a significant enhancement of the water abundance. Water emission peaks along the outflow also correlate with peaks of other shock-produced molecular species, such as SiO and NH3. A strong H2O peak is also observed at the location of the proto-star, where none of the other molecules have significant emission. The absolute 179 um intensity and its intensity ratio to the H2O 557 GHz line previously observed with Odin/SWAS indicate that the water emission originates in warm compact clumps, spatially unresolved by PACS, having a H2O abundance of the order of 10^-4. This testifies that the clumps have been heated for a time long enough to allow the conversion of almost all the available gas-phase oxygen into water. The total water cooling is ~10^-1 Lo, about 40% of the cooling due to H2 and 23% of the total energy released in shocks along the L1157 outflow.

We present results of a 3D model of optical to gamma-ray emission from the slot gap accelerator of a rotation-powered pulsar. Primary electrons accelerating to high-altitudes in the unscreened electric field of the slot gap reach radiation-reaction limited Lorentz factors of 2 x 10^7, while electron-positron pairs from lower-altitude cascades flow along field lines interior to the slot gap. The curvature, synchrotron and inverse Compton radiation of both primary electrons and pairs produce a broad spectrum of emission from infra-red to GeV energies. Both primaries and pairs undergo cyclotron resonant absorption of radio photons, allowing them to maintain significant pitch angles. Synchrotron radiation from pairs with a power-law energy spectrum with Lorentz factors 10^2 - 10^5, dominate the spectrum up to 10 MeV. Synchrotron and curvature radiation of primaries dominates from 10 MeV up to a few GeV. We examine the energy-dependent pulse profiles and phase-resolved spectra for parameters of the Crab pulsar as a function of magnetic inclination and viewing angle, comparing to broad-band data. In most cases, the pulse profiles are dominated by caustics on trailing field lines. We also explore the relation of the high-energy and the radio profiles, as well as the possibility of caustic formation in the radio cone emission. We find that the Crab pulsar profiles and spectrum can be reasonably well reproduced by a model with viewing angle 45 degrees and inclination angle 100 or 80 degrees. This model predicts that the slot gap emission below 200 MeV will exhibit correlations in time and phase with the radio emission.

GRLite and GRTensorJ are first and second generation graphical user interfaces to the computer algebra system GRTensorII. Current development centers on GRTensorJ, which provides fully customizable symbolic procedures that reduce many complex calculations to "elementary functions". Although still in development, GRTensorJ, which is now available (free of charge) over the internet, is sufficiently advanced to be of interest to researchers in general relativity and related fields.