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.

The G2 object has recently passed its pericenter passage in our Galactic Center. While the $Br_\gamma$ emission shows clear signs of tidal interaction, the change in the observed luminosity is only of about a factor of 2, in contention with all previous predictions. We present high resolution simulations performed with the moving mesh code, RICH, together with simple analytical arguments that reproduce the observed $Br_\gamma$ emission. In our model, G2 is a gas cloud that undergoes tidal disruption in a dilute ambient medium. We find that during pericenter passage, the efficient cooling of the cloud results in a vertical collapse, compressing the cloud by a factor of $\sim5000$. By properly taking into account the ionization state of the gas, we find that the cloud is UV starved and are able to reproduce the observed $Br_\gamma$ luminosity. For densities larger than $\approx500\;\mathrm{cm}^{-3}$ at pericenter, the cloud fragments, due to cooling instabilities and the emitted radiation is inconsistent with observations. For lower densities, the cloud survives the pericenter passage intact and its emitted radiation matches the observed lightcurve. From the duration of $Br_\gamma$ emission which contains both redshifted and blueshifted components, we show that the cloud is not spherical but rather elongated with a size ratio of 4 at year 2001. The simulated cloud's elongation grows as it travels towards pericenter and is consistent with observations, due to viewing angles. The simulation is also consistent with having a spherical shape at apocenter.

I explore the thermal evolution and ionization states in gas cooling from an initially hot state in the presence of external photoionizing radiation. I compute the equilibrium and nonequilibrium cooling efficiencies, heating rates, and ion fractions for low-density gas cooling while exposed to the ionizing metagalactic background radiation at various redshifts (z=0 - 3), for a range of temperatures (1e8-1e4 K), densities (1e-7 - 1e3 cm^-3) and metallicities (1e-3 - 2 times solar). The results indicate the existence of a threshold ionization parameter, above which the cooling efficiencies are very close to those in photoionization equilibrium (so that departures from equilibrium may be neglected), and below which the cooling efficiencies resemble those in collisional time-dependent gas cooling with no external radiation (and are thus independent of density)

With the advent of new wide-field, high-cadence optical transient surveys, our understanding of the diversity of core-collapse supernovae has grown tremendously in the last decade. However, the pre-supernova evolution of massive stars, that sets the physical backdrop to these violent events, is theoretically not well understood and difficult to probe observationally. Here we report the discovery of the supernova iPTF13dqy = SN 2013fs, a mere ~3 hr after explosion. Our rapid follow-up observations, which include multiwavelength photometry and extremely early (beginning at ~6 hr post-explosion) spectra, map the distribution of material in the immediate environment (<~ 10^15 cm) of the exploding star and establish that it was surrounded by circumstellar material (CSM) that was ejected during the final ~1 yr prior to explosion at a high rate, around 10^-3 solar masses per year. The complete disappearance of flash-ionised emission lines within the first several days requires that the dense CSM be confined to within <~ 10^15 cm, consistent with radio non-detections at 70--100 days. The observations indicate that iPTF13dqy was a regular Type II SN; thus, the finding that the probable red supergiant (RSG) progenitor of this common explosion ejected material at a highly elevated rate just prior to its demise suggests that pre-supernova instabilities may be common among exploding massive stars.

It is widely recognized that nucleosynthetic output of the first, Population III supernovae was a catalyst defining the character of subsequent stellar generations. Most of the work on the earliest enrichment was carried out assuming that the first stars were extremely massive and that the associated supernovae were unusually energetic, enough to completely unbind the baryons in the host cosmic minihalo and disperse the synthesized metals into the intergalactic medium. Recent work, however, suggests that the first stars may in fact have been somewhat less massive, with a characteristic mass scale of a few tens of solar masses. We present a cosmological simulation following the transport of the metals synthesized in a Population III supernova assuming that it had an energy of 1e51 ergs, compatible with standard Type II supernovae. A young supernova remnant is inserted in the first star's relic HII region in the free expansion phase and is followed for 40 Myr employing adaptive mesh refinement and Lagrangian tracer particle techniques. The supernova remnant remains partially trapped within the minihalo and the thin snowplow shell develops pronounced instability and fingering. Roughly half of the ejecta turn around and fall back toward the center of the halo, with 1% of the ejecta reaching the center in 30 kyr and 10% in 10 Myr. The average metallicity of the combined returning ejecta and the pristine filaments feeding into the halo center from the cosmic web is 0.001 - 0.01 Z_sun, but the two remain unmixed until accreting onto the central hydrostatic core that is unresolved at the end of the simulation. We conclude that if Population III stars had less extreme masses, they promptly enriched the host minihalos with metals and triggered Population II star formation.

Based on high-resolution, spatially resolved data of 10 z~2 star-forming galaxies from the SINS/zC-SINF survey and LUCI data for 12 additional galaxies, we probe the excitation properties of high-z galaxies and the impact of active galactic nuclei (AGN), shocks and photoionization. We explore how these spatially-resolved line ratios can inform our interpretation of integrated emission line ratios obtained at high redshift. Many of our galaxies fall in the `composite' region of the z~0 [NII]/Ha versus [OIII]/Hb diagnostic (BPT) diagram, between star-forming galaxies and those with AGN. Based on our resolved measurements, we find that some of these galaxies likely host an AGN, while others appear to be affected by the presence of shocks possibly caused by an outflow or from enhanced ionization parameter as compared with HII regions in normal local star-forming galaxies. We find that the Mass-Excitation (MEx) diagnostic, which separates purely star-forming and AGN hosting local galaxies in the [OIII]/Hb versus stellar mass plane, does not properly separate z~2 galaxies classified according to the BPT diagram. However, if we shift the galaxies based on the offset between the local and z~2 mass-metallicity relation (i.e. to the mass they would have at z~0 with the same metallicity), we find better agreement between the MEx and BPT diagnostics. Finally, we find that metallicity calibrations based on [NII]/Ha are more biased by shocks and AGN at high-z than the [OIII]/Hb/[NII]/Ha calibration.

We discuss recent improvements in the calculation of the radiative cooling in both collisionally and photo ionized plasmas. We are extending the spectral simulation code Cloudy so that as much as possible of the underlying atomic data is taken from external databases, some created by others, some developed by the Cloudy team. This paper focuses on recent changes in the treatment of many stages of ionization of iron, and discusses its extensions to other elements. The H-like and He-like ions are treated in the iso-electronic approach described previously. Fe II is a special case treated with a large model atom. Here we focus on Fe III through Fe XXIV, ions which are important contributors to the radiative cooling of hot, 1e5 to 1e7 K, plasmas and for X-ray spectroscopy. We use the Chianti atomic database to greatly expand the number of transitions in the cooling function. Chianti only includes lines that have atomic data computed by sophisticated methods. This limits the line list to lower excitation, longer wavelength, transitions. We had previously included lines from the Opacity Project database, which tends to include higher energy, shorter wavelength, transitions. These were combined with various forms of the g-bar approximation, a highly approximate method of estimating collision rates. For several iron ions the two databases are almost entirely complementary. We adopt a hybrid approach in which we use Chianti where possible, supplemented by lines from the Opacity Project for shorter wavelength transitions. The total cooling including the lightest thirty elements differs significantly from some previous calculations.

There is a consensus that Type-Ia supernovae (SNe Ia) arise from the thermonuclear explosion of white dwarf stars that accrete matter from a binary companion. However, direct observation of SN Ia progenitors is lacking, and the precise nature of the binary companion remains uncertain. A temporal series of high-resolution optical spectra of the SN Ia PTF 11kx reveals a complex circumstellar environment that provides an unprecedentedly detailed view of the progenitor system. Multiple shells of circumsteller are detected and the SN ejecta are seen to interact with circumstellar material (CSM) starting 59 days after the explosion. These features are best described by a symbiotic nova progenitor, similar to RS Ophiuchi.

Cold gas streaming along the dark-matter filaments of the cosmic web is predicted to be the major provider of resources for disc buildup, violent disk instability and star formation in massive galaxies in the early universe. We study to what extent these cold streams are traceable in the extended circum-galactic environment of galaxies via Ly alpha absorption and selected low ionisation metal absorption lines. We model the expected absorption signatures using high resolution zoom-in AMR cosmological simulations. In the postprocessing, we distinguish between self-shielded gas and unshielded gas. In the self-shielded gas, which is optically thick to Lyman continuum radiation, we assume pure collisional ionisation for species with an ionisation potential greater than 13.6 eV. In the optically thin, unshielded gas these species are also photoionised by the meta-galactic radiation. In addition to absorption of radiation from background quasars, we compute the absorption line profiles of radiation emitted by the galaxy at the centre of the same halo. We predict the strength of the absorption signal for individual galaxies without stacking. We find that the Ly alpha absorption profiles produced by the streams are consistent with observations of absorption and emission Ly alpha profiles in high redshift galaxies. Due to the low metallicities in the streams, and their low covering factors, the metal absorption features are weak and difficult to detect.

(Abridged). The optical light curve of some SNe may be powered by the outward diffusion of the energy deposited by the explosion shock in optically thick circumstellar matter (CSM). Recently, it was shown that the radiation-mediated and -dominated shock in an optically thick wind must transform into a collisionless shock and can produce hard X-rays. The X-rays are expected to peak at late times, relative to maximum visible light. Here we report on a search, using Swift and Chandra, for X-ray emission from 28 SNe that belong to classes whose progenitors are suspected to be embedded in dense CSM (IIn/Ibn/SLSN-I). Two SNe in our sample have X-ray properties that are roughly consistent with the expectation for X-rays from a collisionless shock in optically thick CSM. Therefore, we suggest that their optical light curves are powered by shock breakout in CSM. We show that two other events were too X-ray bright during the SN maximum optical light to be explained by the shock breakout model. We conclude that the light curves of some, but not all, type-IIn/Ibn SNe are powered by shock breakout in CSM. For the rest of the SNe in our sample, including all the SLSN-I events, our X-ray limits are not deep enough and were typically obtained at too early times to conclude about their nature. We argue that the optical light curves of SNe, for which the X-ray emission peaks at late times, are likely powered by the diffusion of shock energy from a dense CSM. We comment about the possibility to detect some of these events in radio.

We present ion-by-ion cooling efficiencies for low-density gas. We use Cloudy (ver. 08.00) to estimate the cooling efficiencies for each ion of the first 30 elements (H-Zn) individually. We present results for gas temperatures between 1e4 and 1e8K, assuming low densities and optically thin conditions. When nonequilibrium ionization plays a significant role the ionization states deviate from those that obtain in collisional ionization equilibrium (CIE), and the local cooling efficiency at any given temperature depends on specific non-equilibrium ion fractions. The results presented here allow for an efficient estimate of the total cooling efficiency for any ionic composition. We also list the elemental cooling efficiencies assuming CIE conditions. These can be used to construct CIE cooling efficiencies for non-solar abundance ratios, or to estimate the cooling due to elements not explicitly included in any nonequilibrium computation. All the computational results are listed in convenient online tables.

I present computations of the integrated column densities produced in the post-shock cooling layers and in the radiative precursors of partially-cooled fast shocks as a function of the shock age. The results are applicable to the shock-heated warm/hot intergalactic medium (WHIM) which is expected to be a major baryonic reservoir, and contain a large fraction of the so-called "missing baryons". My computations indicate that readily observable amounts of intermediate and high ions, such as CIV, NV, and OVI are created in the precursors of young shocks, for which the shocked gas remains hot and difficult to observe. I suggest that such precursors may provide a way to identify and estimate the "missing" baryonic mass associated with the shocks. The absorption-line signatures predicted here may be used to construct ion-ratio diagrams, which will serve as diagnostics for the photoionized gas in the precursors. In my numerical models, the time-evolution of the shock structure, self-radiation, and associated metal-ion column densities are computed by a series of quasi-static models, each appropriate for a different shock age. The shock code used in this work calculates the nonequilibrium ionization and cooling, follows the radiative transfer of the shock self-radiation through the post-shock cooling layers, takes into account the resulting photoionization and heating rates, follows the dynamics of the cooling gas, and self-consistently computes the photoionization states in the precursor gas. I present a complete set of the age-dependent post-shock and precursor columns for all ionization states of the elements H, He, C, N, O, Ne, Mg, Si, S, and Fe, as functions of the shock velocity, gas metallicity, and magnetic field. I present my numerical results in convenient online tables.

We present numerical computations of the equilibrium configurations of tidally-locked homogeneous binaries, rotating in circular orbits. Unlike the classical Roche approximations, we self-consistently account for the tidal and rotational deformations of both components, and relax the assumptions of ellipsoidal configurations and Keplerian rotation. We find numerical solutions for mass ratios q between 1e-3 and 1, starting at a small angular velocity for which tidal and rotational deformations are small, and following a sequence of increasing angular velocities. Each series terminates at an appropriate ``Roche limit'', above which no equilibrium solution can be found. Even though the Roche limit is crossed before the ``Roche lobe'' is filled, any further increase in the angular velocity will result in mass-loss. For close, comparable-mass binaries, we find that local deviations from ellipsoidal forms may be as large as 10-20%, and departures from Keplerian rotation are significant. We compute the light curves that arise from our equilibrium configurations, assuming their distance is >>1 AU (e.g. in the Kuiper Belt). We consider both backscatter (proportional to the projected area) and diffuse (Lambert) reflections. Backscatter reflection always yields two minima of equal depths. Diffuse reflection, which is sensitive to the surface curvature, generally gives rise to unequal minima. We find detectable intensity differences of up to 10% between our light curves and those arising from the Roche approximations. Finally, we apply our models to Kuiper Belt binary 2001 QG298, and find a nearly edge-on binary with a mass ratio q = 0.93 ^{+0.07}_{-0.03}, angular velocity Omega^2/G rho = 0.333+/-0.001 (statistical errors only), and pure diffuse reflection. For the observed period of 2001 QG298, these parameters imply a bulk density, rho = 0.72 +/- 0.04 g cm^-3.

We present computations of the ionization structure and metal-absorption properties of thermally conductive interface layers that surround evaporating warm spherical clouds, embedded in a hot medium. We rely on the analytical formalism of Dalton & Balbus to calculate the temperature profile in the evaporating gas, and explicitly solve the time-dependent ionization equations for H, He, C, N, O, Si, and S in the interface. We include photoionization by an external field. We estimate how departures from equilibrium ionization affect the resonance-line cooling efficiencies in the evaporating gas, and determine the conditions for which radiative losses may be neglected in the solution for the evaporation dynamics and temperature profile. Our results indicate that non-equilibrium cooling significantly increases the value of the saturation parameter at which radiative losses begin to affect the flow dynamics. As applications we calculate the ion fractions and projected column densities arising in the evaporating layers surrounding dwarf-galaxy-scale objects that are also photoionized by metagalactic radiation. We compare our results to the UV metal-absorption column densities observed in local highly-ionized metal-absorbers, located in the Galactic corona or intergalactic medium. Conductive interfaces significantly enhance the formation of high-ions such as C^3+, N^4+, and O^5+ relative to purely photoionized clouds, especially for clouds embedded in a high-pressure corona. However, the enhanced columns are still too low to account for the O VI columns (~1e14 cm^-2) observed in the local high-velocity absorbers. We find that O VI columns larger than ~1e13 cm^-2 cannot be produced in evaporating clouds. Our results do support the conclusion of Savage & Lehner, that absorption due to evaporating O VI likely occurs in the local interstellar medium, with characteristic columns of ~1e13 cm^-2.

Recent observations have revealed that some Type Ia supernovae exhibit narrow, time-variable Na I D absorption features. The origin of the absorbing material is controversial, but it may suggest the presence of circumstellar gas in the progenitor system prior to the explosion, with significant implications for the nature of the supernova progenitors. We present the third detection of such variable absorption, based on six epochs of high-resolution spectroscopy of the Type Ia supernova SN 2007le from Keck and the HET. The data span ~3 months, from 5 days before maximum light to 90 days after maximum. We find that one component of the Na D absorption lines strengthened significantly with time, indicating a total column density increase of ~2.5 x 10^12 cm^-2. The changes are most prominent after maximum light rather than at earlier times when the UV flux from the SN peaks. As with SN 2006X, we detect no change in the Ca II H&K lines over the same time period, rendering line-of-sight effects improbable and suggesting a circumstellar origin for the absorbing material. Unlike the previous two SNe exhibiting variable absorption, SN 2007le is not highly reddened (E_B-V = 0.27 mag), also pointing toward circumstellar rather than interstellar absorption. Photoionization models show that the data are consistent with a dense (10^7 cm^-3) cloud or clouds of gas located ~0.1 pc from the explosion. These results broadly support the single-degenerate scenario previously proposed to explain the variable absorption, with mass loss from a nondegenerate companion star responsible for providing the circumstellar gas. We also present tentative evidence for narrow Halpha emission associated with the SN, which will require followup observations at late times to confirm. [abridged]

In this paper we present computations of the integrated metal-ion column densities produced in the post-shock cooling layers behind fast, radiative shock-waves. For this purpose, we have constructed a new shock code that calculates the non-equilibrium ionization and cooling; follows the radiative transfer of the shock self-radiation through the post-shock cooling layers; takes into account the resulting photoionization and heating rates; follows the dynamics of the cooling gas; and self-consistently computes the initial photoionization state of the precursor gas. We discuss the shock structure and emitted radiation, and study the dependence on the shock velocity, magnetic field, and gas metallicity. We present a complete set of integrated post-shock and precursor metal-ion column densities of all ionization stages of the elements H, He, C, N, O, Ne, Mg, Si, S, and Fe, for shocks with velocities of 600 and ~2000 km/s, corresponding to initial post-shock temperatures of 5e6 and 5e7 K, cooling down to 1000 K. We consider shocks in which the magnetic field is negligible (B=0) so that the cooling occurs at approximately constant pressure ("isobaric"), and shocks in which the magnetic pressure dominates everywhere such that the cooling occurs at constant density (isochoric). We present results for gas metallicities Z ranging from 1e-3 to twice the solar abundance of heavy elements, and we study how the observational signatures of fast radiative shocks depend on Z. We present our numerical results in convenient online figures and tables, available at http://wise-obs.tau.ac.il/~orlyg/shocks/

We establish a simple model for the distribution of cold gas around L* galaxies using a large set of observational constraints on the properties of strong MgII absorber systems. Our analysis suggests that the halos of L* galaxies are filled with cool gaseous clouds having sizes of order 1kpc and densities of ~10^{-2} cm^{-3}. We then investigate the physical effects of cloud irradiation by a quasar and study the resulting spectral signatures. We show that quasar activity gives rise to (i) extended narrow-line emission on ~100kpc scales and (ii) an anisotropy in the properties of the absorbing gas arising from the geometry of the quasar radiation field. Provided that quasars reside in halos several times more massive than those of L* galaxies, our model predictions appear to be in agreement with observations of narrow emission-line nebulae around quasars and the recent detections of ~100kpc cold gaseous envelopes around those objects, suggesting a common origin for these phenomena. We discuss the implications of our results for understanding absorption systems, probing quasar environments at high redshifts, and testing the quasar unification scheme.

We present new computations of the equilibrium and non-equilibrium cooling efficiencies and ionization states for low-density radiatively cooling gas containing cosmic abundances of the elements H, He, C, N, O, Ne, Mg, Si, S, and Fe. We present results for gas temperatures between 1e4 and 1e8 K, assuming dust-free and optically thin conditions, and no external radiation. For non-equilibrium cooling we solve the coupled time-dependent ionization and energy loss equations for a radiating gas cooling from an initially hot, >5e6K equilibrium state, down to 1e4K. We present results for heavy element compositions ranging from 1e-3 to 2 times the elemental abundances in the Sun. We consider gas cooling at constant density (isochoric) and at constant pressure (isobaric). We calculate the critical column densities and temperatures at which radiatively cooling clouds make the dynamical transition from isobaric to isochoric evolution. We construct ion ratio diagnostics for the temperature and metallicity in radiatively cooling gas. We provide numerical estimates for the maximal cloud column densities for which the gas remains optically thin to the cooling radiation. We present our computational results in convenient on-line figures and tables (see http://wise-obs.tau.ac.il/~orlyg/cooling/).

(Abridged) We present computations of the metals photoionization structures of pressure supported gas clouds in gravitationally dominant DM minihalos in the extended Galactic halo or LG environment. We consider low-metallicity (0.1-0.3 solar) clouds that are photoionized by the metagalactic radiation field, and which are pressure-confined by a hot external medium. We study the properties of LCDM Burkert (and NFW) minihalos with circular velocities 12-30 km/s. We present results for the volume densities and projected column densities of low-ions, CII, NI, OI, SiII and SII, produced in the neutral cores, and high-ions, CIV, NV, OVI, SiIII, SiIV, and SIII produced in the ionized envelopes. We examine the possible relationships between the compact HI high-velocity clouds (CHVCs), dwarf galaxies, and high-velocity ionized CIV absorbers, within the context of photoionized minihalo models. In pressure-confined (P/k ~50 K/cm^3) mini-halo models for the CHVCs the photoionization states are much lower than in the CIV absorbers. However, for lower bounding pressures (<~1 K/cm^3) in more massive dwarf galaxy scale halos, the photoionization states in the outer envelopes resemble those observed in the CIV HVCs. An important exception is OVI for which our photoionization models underpredict the relative abundances by about an order-of-magnitude, suggesting that additional processes are at play as concluded by Sembach et al. For the cloud size-scales expected in median LCDM Burkert halos, a metallicity of 0.3 solar is required to produce absorption-line strengths comparable to those in the CIV HVCs. We argue that fully ionized and starless "dark galaxies" could be detectable in the local universe as UV metal line absorbers with ionization states similar to the CIV absorbers.

We present a radial-velocity study of the triple-lined system Gliese 644 and derive spectroscopic elements for the inner and outer orbits with periods of 2.9655 and 627 days. We also utilize old visual data, as well as modern speckle and adaptive optics observations, to derive a new astrometric solution for the outer orbit. These two orbits together allow us to derive masses for each of the three components in the system: M_A = 0.410 +/- 0.028 (6.9%), M_Ba = 0.336 +/- 0.016 (4.7%), and $M_Bb = 0.304 +/- 0.014 (4.7%) M_solar. We suggest that the relative inclination of the two orbits is very small. Our individual masses and spectroscopic light ratios for the three M stars in the Gliese 644 system provide three points for the mass-luminosity relation near the bottom of the Main Sequence, where the relation is poorly determined. These three points agree well with theoretical models for solar metallicity and an age of 5 Gyr. Our radial velocities for Gliese 643 and vB 8, two common-proper-motion companions of Gliese 644, support the interpretation that all five M stars are moving together in a physically bound group. We discuss possible scenarios for the formation and evolution of this configuration, such as the formation of all five stars in a sequence of fragmentation events leading directly to the hierarchical configuration now observed, versus formation in a small N cluster with subsequent dynamical evolution into the present hierarchical configuration.