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.

We present a survey of extremely high-velocity outflows (EHVOs) in quasars, defined by speeds between 0.1c and 0.2c. This region of the parameter space has not been included in previous surveys, but it might present the biggest challenge for theoretical models and it might be a large contributor to feedback due to the outflows' potentially large kinetic power. Using the Sloan Digital Sky Survey, we find 40 quasar spectra with broad EHVO CIV absorption, 10 times more than the number of previously known cases. We characterize the EHVO absorption and find that in 26 cases, CIV is accompanied by NV and/or OVI absorption. We find that EHVO quasars lack HeII emission and have overall larger bolometric luminosities and black hole masses than those of their parent sample and BALQSOs, while we do not find significant differences in their Eddington ratios. We also report a trend toward larger black hole masses as the velocity of the outflowing gas increases in the BALQSOs in our sample. The overall larger Lbol and lack of HeII emission of EHVO quasars suggest that radiation is likely driving these outflows. We find a potential evolutionary effect as EHVO quasars seem to be more predominant at large redshifts. We estimate that the kinetic power of these outflows may be similar to or even larger than that of the outflows from BALQSOs as the velocity factor increases this parameter by 1-2.5 orders of magnitude. Further study of EHVO quasars will help improve our understanding of quasar physics.

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 report the formation of bound star clusters in a sample of high-resolution cosmological zoom-in simulations of z>5 galaxies from the FIRE project. We find that bound clusters preferentially form in high-pressure clouds with gas surface densities over 10^4 Msun pc^-2, where the cloud-scale star formation efficiency is near unity and young stars born in these regions are gravitationally bound at birth. These high-pressure clouds are compressed by feedback-driven winds and/or collisions of smaller clouds/gas streams in highly gas-rich, turbulent environments. The newly formed clusters follow a power-law mass function of dN/dM~M^-2. The cluster formation efficiency is similar across galaxies with stellar masses of ~10^7-10^10 Msun at z>5. The age spread of cluster stars is typically a few Myrs and increases with cluster mass. The metallicity dispersion of cluster members is ~0.08 dex in [Z/H] and does not depend on cluster mass significantly. Our findings support the scenario that present-day old globular clusters (GCs) were formed during relatively normal star formation in high-redshift galaxies. Simulations with a stricter/looser star formation model form a factor of a few more/fewer bound clusters per stellar mass formed, while the shape of the mass function is unchanged. Simulations with a lower local star formation efficiency form more stars in bound clusters. The simulated clusters are larger than observed GCs due to finite resolution. Our simulations are among the first cosmological simulations that form bound clusters self-consistently in a wide range of high-redshift galaxies.

Radiative feedback (RFB) from stars plays a key role in galaxies, but remains poorly-understood. We explore this using high-resolution, multi-frequency radiation-hydrodynamics (RHD) simulations from the Feedback In Realistic Environments (FIRE) project. We study ultra-faint dwarf through Milky Way mass scales, including H+He photo-ionization; photo-electric, Lyman Werner, Compton, and dust heating; and single+multiple scattering radiation pressure (RP). We compare distinct numerical algorithms: ray-based LEBRON (exact when optically-thin) and moments-based M1 (exact when optically-thick). The most important RFB channels on galaxy scales are photo-ionization heating and single-scattering RP: in all galaxies, most ionizing/far-UV luminosity (~1/2 of lifetime-integrated bolometric) is absorbed. In dwarfs, the most important effect is photo-ionization heating from the UV background suppressing accretion. In MW-mass galaxies, meta-galactic backgrounds have negligible effects; but local photo-ionization and single-scattering RP contribute to regulating the galactic star formation efficiency and lowering central densities. Without some RFB (or other 'rapid' FB), resolved GMCs convert too-efficiently into stars, making galaxies dominated by hyper-dense, bound star clusters. This makes star formation more violent and 'bursty' when SNe explode in these hyper-clustered objects: thus, including RFB 'smoothes' SFHs. These conclusions are robust to RHD methods, but M1 produces somewhat stronger effects. Like in previous FIRE simulations, IR multiple-scattering is rare (negligible in dwarfs, ~10% of RP in massive galaxies): absorption occurs primarily in 'normal' GMCs with A_v~1.

The correlation between galaxies' integrated stellar masses and star formation rates (the `star formation main sequence'; SFMS) is a well-established scaling relation. Recently, surveys have found a relationship between the star formation rate and stellar mass surface densities on kpc and sub-kpc scales (the `resolved SFMS'; rSFMS). In this work, we demonstrate that the rSFMS emerges naturally in FIRE-2 zoom-in simulations of Milky Way-mass galaxies. We make SFR and stellar mass maps of the simulated galaxies at a variety of spatial resolutions and star formation averaging time-scales and fit the rSFMS using multiple methods from the literature. While the absolute value of the SFMS slope depends on the fitting method, the slope is steeper for longer star formation time-scales and lower spatial resolutions regardless of the fitting method employed. We present a toy model that quantitatively captures the dependence of the simulated galaxies' rSFMS slope on spatial resolution and use it to illustrate how this dependence can be used to constrain the characteristic mass of star-forming clumps.

We present models of CO(1-0) emission from Milky Way-mass galaxies at redshift zero in the FIRE-2 cosmological zoom-in simulations. We calculate the molecular abundances by post-processing the simulations with an equilibrium chemistry solver while accounting for the effects of local sources, and determine the emergent CO(1-0) emission using a line radiative transfer code. We find that the results depend strongly on the shielding length assumed, which in our models sets the attenuation of the incident UV radiation field. Commonly used choices for the shielding length, such as the Jeans length, result in CO abundances that are too high at a given H$_2$ abundance. We find that a model with a distribution of shielding lengths, which has a median shielding length of $\sim 3$ pc in cold gas ($T < 300$ K) for both CO and H$_{2}$, is able to reproduce both the observed CO(1-0) luminosity and inferred CO-to-H$_{2}$ conversion factor at a given star formation rate compared with observations. We suggest that this short shielding length can be thought of as a subgrid model which controls the amount of radiation that penetrates giant molecular clouds.

Supermassive black holes (SMBHs) are thought to provide energy that prevents catastrophic cooling in the centers of massive galaxies and galaxy clusters. However, it remains unclear how this "feedback" process operates. We use high-resolution optical data to study the kinematics of multi-phase filamentary structures by measuring the velocity structure function (VSF) of the filaments over a wide range of scales in the centers of three nearby galaxy clusters: Perseus, Abell 2597 and Virgo. We find that the motions of the filaments are turbulent in all three clusters studied. There is a clear correlation between features of the VSFs and the sizes of bubbles inflated by SMBH driven jets. Our study demonstrates that SMBHs are the main driver of turbulent gas motions in the centers of galaxy clusters and suggests that this turbulence is an important channel for coupling feedback to the environment. Our measured amplitude of turbulence is in good agreement with Hitomi Doppler line broadening measurement and X-ray surface brightness fluctuation analysis, suggesting that the motion of the cold filaments is well-coupled to that of the hot gas. The smallest scales we probe are comparable to the mean free path in the intracluster medium (ICM). Our direct detection of turbulence on these scales provides the clearest evidence to date that isotropic viscosity is suppressed in the weakly-collisional, magnetized intracluster plasma.

We present the first measurement of the lifetimes of Giant Molecular Clouds (GMCs) in cosmological simulations at $z = 0$, using the Latte suite of FIRE-2 simulations of Milky Way-mass galaxies. We track GMCs with total gas mass $\gtrsim 10^5$ M$_\odot$ at high spatial ($\sim1$ pc), mass ($7100$ M$_{\odot}$), and temporal (1 Myr) resolution. Our simulated GMCs are consistent with the distribution of masses for massive GMCs in the Milky Way and nearby galaxies. We find GMC lifetimes of $5-7$ Myr, or 1-2 freefall times, on average, with less than 1\% of clouds living longer than 20 Myr. We find increasing GMC lifetimes with galactocentric radius, implying that environment affects the evolutionary cycle of GMCs. However, our GMC lifetimes show no systematic dependence on GMC mass or amount of star formation. These results are broadly consistent with inferences from the literature and provide an initial investigation into ultimately understanding the physical processes that govern GMC lifetimes in a cosmological setting.

We study the spatially resolved (sub-kpc) gas velocity dispersion ($\sigma$)--star formation rate (SFR) relation in the FIRE-2 (Feedback in Realistic Environments) cosmological simulations. We specifically focus on Milky Way mass disk galaxies at late times. In agreement with observations, we find a relatively flat relationship, with $\sigma \approx 15-30$ km/s in neutral gas across 3 dex in SFRs. We show that higher dense gas fractions (ratios of dense gas to neutral gas) and SFRs are correlated at constant $\sigma$. Similarly, lower gas fractions (ratios of gas to stellar mass) are correlated with higher $\sigma$ at constant SFR. The limits of the $\sigma$-$\Sigma_{\rm SFR}$ relation correspond to the onset of strong outflows. We see evidence of "on-off" cycles of star formation in the simulations, corresponding to feedback injection timescales of 10-100 Myr, where SFRs oscillate about equilibrium SFR predictions. Finally, SFRs and velocity dispersions in the simulations agree well with feedback-regulated and marginally stable gas disk (Toomre's $Q =1$) model predictions, and the data effectively rule out models assuming that gas turns into stars at (low) constant efficiency (i.e., ${\rm 1\%}$ per free-fall time). And although the simulation data do not entirely exclude gas accretion/gravitationally powered turbulence as a driver of $\sigma$, it appears to be strongly subdominant to stellar feedback in the simulated galaxy disks.

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.

The Cerro Chajnantor Atacama Telescope-prime (CCAT-prime) is a new 6-m, off-axis, low-emissivity, large field-of-view submillimeter telescope scheduled for first light in the last quarter of 2021. In summary, (a) CCAT-prime uniquely combines a large field-of-view (up to 8-deg), low emissivity telescope (< 2%) and excellent atmospheric transmission (5600-m site) to achieve unprecedented survey capability in the submillimeter. (b) Over five years, CCAT-prime first generation science will address the physics of star formation, galaxy evolution, and galaxy cluster formation; probe the re-ionization of the Universe; improve constraints on new particle species; and provide for improved removal of dust foregrounds to aid the search for primordial gravitational waves. (c) The Observatory is being built with non-federal funds (~ \$40M in private and international investments). Public funding is needed for instrumentation (~ \$8M) and operations (\$1-2M/yr). In return, the community will be able to participate in survey planning and gain access to curated data sets. (d) For second generation science, CCAT-prime will be uniquely positioned to contribute high-frequency capabilities to the next generation of CMB surveys in partnership with the CMB-S4 and/or the Simons Observatory projects or revolutionize wide-field, sub-millimetter line intensity mapping surveys.

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.

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.

We develop a framework for joint constraints on the CO luminosity function based on power spectra (PS) and voxel intensity distributions (VID), and apply this to simulations of COMAP, a CO intensity mapping experiment. This Bayesian framework is based on a Markov chain Monte Carlo (MCMC) sampler coupled to a Gaussian likelihood with a joint PS + VID covariance matrix computed from a large number of fiducial simulations, and re-calibrated with a small number of simulations per MCMC step. The simulations are based on dark matter halos from fast peak patch simulations combined with the $L_\text{CO}(M_\text{halo})$ model of Li et al. (2016). We find that the relative power to constrain the CO luminosity function depends on the luminosity range of interest. In particular, the VID is more sensitive at both small and large luminosities, while the PS is more sensitive at intermediate luminosities. The joint analysis is superior to using either observable separately. When averaging over CO luminosities ranging between $L_\text{CO} = 10^4-10^7L_\odot$, and over 10 cosmological realizations of COMAP Phase 2, the uncertainties (in dex) are larger by 58 % and 30 % for the PS and VID, respectively, when compared to the joint analysis (PS + VID). This method is generally applicable to any other random field, with a complicated likelihood, as long a fast simulation procedure is available.

Line-intensity mapping (LIM or IM) is an emerging field of observational work, with strong potential to fit into a larger effort to probe large-scale structure and small-scale astrophysical phenomena using multiple complementary tracers. Taking full advantage of such complementarity means, in part, undertaking line-intensity surveys with galaxy surveys in mind. We consider the potential for detection of a cross-correlation signal between COMAP and blind surveys based on photometric redshifts (as in COSMOS) or based on spectroscopic data (as with the HETDEX survey of Lyman-$\alpha$ emitters). We find that obtaining $\sigma_z/(1+z)\lesssim0.003$ accuracy in redshifts and $\gtrsim10^{-4}$ sources per Mpc$^3$ with spectroscopic redshift determination should enable a CO-galaxy cross spectrum detection significance at least twice that of the CO auto spectrum. Either a future targeted spectroscopic survey or a blind survey like HETDEX may be able to meet both of these requirements.

Recent analyses of the gamma-ray spectrum from the ultra-luminous infrared galaxy Arp 220 have revealed a discrepancy in the cosmic ray energy injection rates derived from the gamma-rays versus the radio emission. While the observed radio emission is consistent with the star formation rate inferred from infrared observations, a significantly higher cosmic ray population is necessary to accurately model the measured gamma-ray flux. To resolve this discrepancy between the radio and gamma-ray observations, we find that we must increase the cosmic ray energy injection rate and account for an infrared optical depth greater than unity. Raising the energy injection rate naturally raises the total gamma-ray flux but also raises the radio flux unless there is also an increase in the energy loss rate for cosmic ray leptons. A optically thick medium results in an increase in energy losses via inverse Compton for cosmic ray leptons and preserves agreement with submillimeter, millimeter, and infrared wavelength observations.

The Feedback In Realistic Environments (FIRE) project explores feedback in cosmological galaxy formation simulations. Previous FIRE simulations used an identical source code (FIRE-1) for consistency. Motivated by the development of more accurate numerics - including hydrodynamic solvers, gravitational softening, and supernova coupling algorithms - and exploration of new physics (e.g. magnetic fields), we introduce FIRE-2, an updated numerical implementation of FIRE physics for the GIZMO code. We run a suite of simulations and compare against FIRE-1: overall, FIRE-2 improvements do not qualitatively change galaxy-scale properties. We pursue an extensive study of numerics versus physics. Details of the star-formation algorithm, cooling physics, and chemistry have weak effects, provided that we include metal-line cooling and star formation occurs at higher-than-mean densities. We present new resolution criteria for high-resolution galaxy simulations. Most galaxy-scale properties are robust to numerics we test, provided: (1) Toomre masses are resolved; (2) feedback coupling ensures conservation, and (3) individual supernovae are time-resolved. Stellar masses and profiles are most robust to resolution, followed by metal abundances and morphologies, followed by properties of winds and circum-galactic media (CGM). Central (~kpc) mass concentrations in massive (L*) galaxies are sensitive to numerics (via trapping/recycling of winds in hot halos). Multiple feedback mechanisms play key roles: supernovae regulate stellar masses/winds; stellar mass-loss fuels late star formation; radiative feedback suppresses accretion onto dwarfs and instantaneous star formation in disks. We provide all initial conditions and numerical algorithms used.

We study the implementation of mechanical feedback from supernovae (SNe) and stellar mass loss in galaxy simulations, within the Feedback In Realistic Environments (FIRE) project. We present the FIRE-2 algorithm for coupling mechanical feedback, which can be applied to any hydrodynamics method (e.g. fixed-grid, moving-mesh, and mesh-less methods), and black hole as well as stellar feedback. This algorithm ensures manifest conservation of mass, energy, and momentum, and avoids imprinting 'preferred directions' on the ejecta. We show that it is critical to incorporate both momentum and thermal energy of mechanical ejecta in a self-consistent manner, accounting for SNe cooling radii when they are not resolved. Using idealized simulations of single SN explosions, we show that the FIRE-2 algorithm, independent of resolution, reproduces converged solutions in both energy and momentum. In contrast, common 'fully-thermal' (energy-dump) or 'fully-kinetic' (particle-kicking) schemes in the literature depend strongly on resolution: when applied at mass resolution >100 solar masses, they diverge by orders-of-magnitude from the converged solution. In galaxy-formation simulations, this divergence leads to orders-of-magnitude differences in galaxy properties, unless those models are adjusted in a resolution-dependent way. We show that all models that individually time-resolve SNe converge to the FIRE-2 solution at sufficiently high resolution. However, in both idealized single-SN simulations and cosmological galaxy-formation simulations, the FIRE-2 algorithm converges much faster than other sub-grid models without re-tuning parameters.

We investigate star forming scaling relations using Bayesian inference on a comprehensive data sample of low- (z<0.1) and high-redshift (1<z<5) star forming regions. This full data set spans a wide range of host galaxy stellar mass ($M_{*} \sim10^6-10^{11} M_{\odot}$) and clump star formation rates (SFR $ \sim10^{-5}-10^2 M_\odot yr^{-1}$). We fit the power-law relationship between the size (r$_{H\alpha}$) and luminosity (L$_{H\alpha}$) of the star forming clumps using the Bayesian statistical modeling tool Stan that makes use of Markov Chain Monte Carlo (MCMC) sampling techniques. Trends in the scaling relationship are explored for the full sample and subsets based on redshift and selection effects between samples. In our investigation we find no evidence of redshift evolution of the size-luminosity scaling relationship, nor a difference in slope between lensed and unlensed data. There is evidence of a break in the scaling relationship between high and low star formation rate surface density ($\Sigma_{SFR}$) clumps. The size-luminosity power law fit results are L$_{H\alpha}\sim$ r$_{H\alpha}^{2.8}$ and L$_{H\alpha}\sim$ r$_{H\alpha}^{1.7}$ for low and high $\Sigma_{SFR}$ clumps, respectively. We present a model where star forming clumps form at locations of gravitational instability and produce an ionized region represented by the Str\"{o}mgren radius. A radius smaller than the scale height of the disk results in a scaling relationship of $L \propto r^3$ (high $\Sigma_{SFR}$ clumps), and a scaling of $L \propto r^2$ (low $\Sigma_{SFR}$ clumps) if the radius is larger than the disk scale height.

Measurements of the star formation efficiency (SFE) of giant molecular clouds (GMCs) in the Milky Way generally show a large scatter, which could be intrinsic or observational. We use magnetohydrodynamic simulations of GMCs (including feedback) to forward-model the relationship between the true GMC SFE and observational proxies. We show that individual GMCs trace broad ranges of observed SFE throughout collapse, star formation, and disruption. Low measured SFEs (<<1%) are "real" but correspond to early stages, the true "per-freefall" SFE where most stars actually form can be much larger. Very high (>>10%) values are often artificially enhanced by rapid gas dispersal. Simulations including stellar feedback reproduce observed GMC-scale SFEs, but simulations without feedback produce 20x larger SFEs. Radiative feedback dominates among mechanisms simulated. An anticorrelation of SFE with cloud mass is shown to be an observational artifact. We also explore individual dense "clumps" within GMCs and show that (with feedback) their bulk properties agree well with observations. Predicted SFEs within the dense clumps are ~2x larger than observed, possibly indicating physics other than feedback from massive (main sequence) stars is needed to regulate their collapse.

Young massive star clusters spanning $\sim 10^4 - 10^8 M_\odot$ in mass have been observed to have similar surface brightness profiles. Recent hydrodynamical simulations of star cluster formation have also produced star clusters with this structure. We argue analytically that this type of mass distribution arises naturally in the relaxation from a hierarchically-clustered distribution of stars into a monolithic star cluster through hierarchical merging. We show that arbitrary initial profiles will tend to converge to a universal profile under hierarchical merging, owing to phase-space mixing obeying certain conservation constraints. We perform $N$-body simulations of a pairwise merger of model star clusters and find that mergers readily produce the shallow surface brightness profiles observed in young massive clusters. Finally, we simulate the relaxation of a hierarchically-clustered mass distribution constructed from an idealized fragmentation model. Assuming only power-law spatial and kinematic scaling relations, these numerical experiments are able to reproduce the surface density profiles of observed young massive star clusters. Thus we provide physical motivation for the structure of young massive clusters within the paradigm of hierarchical star formation. This has important implications for the structure of nascent globular clusters.

The Gemini Infrared Multi-Object Spectrograph (GIRMOS) is a powerful new instrument being built to facility-class standards for the Gemini telescope. It takes advantage of the latest developments in adaptive optics and integral field spectrographs. GIRMOS will carry out simultaneous high-angular-resolution, spatially-resolved infrared ($1-2.4$ $\mu$m) spectroscopy of four objects within a two-arcminute field-of-regard by taking advantage of multi-object adaptive optics. This capability does not currently exist anywhere in the world and therefore offers significant scientific gains over a very broad range of topics in astronomical research. For example, current programs for high redshift galaxies are pushing the limits of what is possible with infrared spectroscopy at $8-10$-meter class facilities by requiring up to several nights of observing time per target. Therefore, the observation of multiple objects simultaneously with adaptive optics is absolutely necessary to make effective use of telescope time and obtain statistically significant samples for high redshift science. With an expected commissioning date of 2023, GIRMOS's capabilities will also make it a key followup instrument for the James Webb Space Telescope when it is launched in 2021, as well as a true scientific and technical pathfinder for future Thirty Meter Telescope (TMT) multi-object spectroscopic instrumentation. In this paper, we will present an overview of this instrument's capabilities and overall architecture. We also highlight how this instrument lays the ground work for a future TMT early-light instrument.

We present the detailed science case, and brief descriptions of the telescope design, site, and first light instrument plans for a new ultra-wide field submillimeter observatory, CCAT-prime, that we are constructing at a 5600 m elevation site on Cerro Chajnantor in northern Chile. Our science goals are to study star and galaxy formation from the epoch of reionization to the present, investigate the growth of structure in the Universe, improve the precision of B-mode CMB measurements, and investigate the interstellar medium and star formation in the Galaxy and nearby galaxies through spectroscopic, polarimetric, and broadband surveys at wavelengths from 200 um to 2 mm. These goals are realized with our two first light instruments, a large field-of-view (FoV) bolometer-based imager called Prime-Cam (that has both camera and an imaging spectrometer modules), and a multi-beam submillimeter heterodyne spectrometer, CHAI. CCAT-prime will have very high surface accuracy and very low system emissivity, so that combined with its wide FoV at the unsurpassed CCAT site our telescope/instrumentation combination is ideally suited to pursue this science. The CCAT-prime telescope is being designed and built by Vertex Antennentechnik GmbH. We expect to achieve first light in the spring of 2021.

CCAT-prime will be a 6-meter aperture telescope operating from sub-mm to mm wavelengths, located at 5600 meters elevation on Cerro Chajnantor in the Atacama Desert in Chile. Its novel crossed-Dragone optical design will deliver a high throughput, wide field of view capable of illuminating much larger arrays of sub-mm and mm detectors than can existing telescopes. We present an overview of the motivation and design of Prime-Cam, a first-light instrument for CCAT-prime. Prime-Cam will house seven instrument modules in a 1.8 meter diameter cryostat, cooled by a dilution refrigerator. The optical elements will consist of silicon lenses, and the instrument modules can be individually optimized for particular science goals. The current design enables both broadband, dual-polarization measurements and narrow-band, Fabry-Perot spectroscopic imaging using multichroic transition-edge sensor (TES) bolometers operating between 190 and 450 GHz. It also includes broadband kinetic induction detectors (KIDs) operating at 860 GHz. This wide range of frequencies will allow excellent characterization and removal of galactic foregrounds, which will enable precision measurements of the sub-mm and mm sky. Prime-Cam will be used to constrain cosmology via the Sunyaev-Zeldovich effects, map the intensity of [CII] 158 $\mu$m emission from the Epoch of Reionization, measure Cosmic Microwave Background polarization and foregrounds, and characterize the star formation history over a wide range of redshifts. More information about CCAT-prime can be found at www.ccatobservatory.org.

One important result from recent large integral field spectrograph (IFS) surveys is that the intrinsic velocity dispersion of galaxies traced by star-forming gas increases with redshift. Massive, rotation-dominated discs are already in place at z~2, but they are dynamically hotter than spiral galaxies in the local Universe. Although several plausible mechanisms for this elevated velocity dispersion (e.g. star formation feedback, elevated gas supply, or more frequent galaxy interactions) have been proposed, the fundamental driver of the velocity dispersion enhancement at high redshift remains unclear. We investigate the origin of this kinematic evolution using a suite of cosmological simulations from the FIRE (Feedback In Realistic Environments) project. Although IFS surveys generally cover a wider range of stellar masses than in these simulations, the simulated galaxies show trends between intrinsic velocity dispersion, SFR, and redshift in agreement with observations. In both the observed and simulated galaxies, intrinsic velocity dispersion is positively correlated with SFR. Intrinsic velocity dispersion increases with redshift out to z~1 and then flattens beyond that. In the FIRE simulations, intrinsic velocity dispersion can vary significantly on timescales of <100 Myr. These variations closely mirror the time evolution of the SFR and gas inflow rate. By cross-correlating pairs of intrinsic velocity dispersion, gas inflow rate, and SFR, we show that increased gas inflow leads to subsequent enhanced star formation, and enhancements in intrinsic velocity dispersion tend to temporally coincide with increases in gas inflow rate and SFR.

We present an analysis of the global and spatially-resolved Kennicutt-Schmidt (KS) star formation relation in the FIRE (Feedback In Realistic Environments) suite of cosmological simulations, including halos with $z = 0$ masses ranging from $10^{10}$ -- $10^{13}$ M$_{\odot}$. We show that the KS relation emerges and is robustly maintained due to the effects of feedback on local scales regulating star-forming gas, independent of the particular small-scale star formation prescriptions employed. We demonstrate that the time-averaged KS relation is relatively independent of redshift and spatial averaging scale, and that the star formation rate surface density is weakly dependent on metallicity and inversely dependent on orbital dynamical time. At constant star formation rate surface density, the `Cold \& Dense' gas surface density (gas with $T < 300$~K and $n > 10$~cm$^{-3}$, used as a proxy for the molecular gas surface density) of the simulated galaxies is $\sim$0.5~dex less than observed at $\sim$kpc scales. This discrepancy may arise from underestimates of the local column density at the particle-scale for the purposes of shielding in the simulations. Finally, we show that on scales larger than individual giant molecular clouds, the primary condition that determines whether star formation occurs is whether a patch of the galactic disk is thermally Toomre-unstable (not whether it is self-shielding): once a patch can no longer be thermally stabilized against fragmentation, it collapses, becomes self-shielding, cools, and forms stars, regardless of epoch or environment.

A maximum stellar surface density $\Sigma_{max} \sim 3 \times 10^5\,{\rm M_\odot\,pc^{-2}}$ is observed across all classes of dense stellar systems (e.g. star clusters, galactic nuclei, etc.), spanning $\sim 8$ orders of magnitude in mass. It has been proposed that this characteristic scale is set by some dynamical feedback mechanism preventing collapse beyond a certain surface density. However, simple analytic models and detailed simulations of star formation moderated by feedback from massive stars argue that feedback becomes {\it less} efficient at higher surface densities (with the star formation efficiency increasing as $\sim \Sigma/\Sigma_{crit}$). We therefore propose an alternative model wherein stellar feedback becomes ineffective at moderating star formation above some $\Sigma_{crit}$, so the supply of star-forming gas is rapidly converted to stars before the system can contract to higher surface density. We show that such a model -- with $\Sigma_{crit}$ taken directly from the theory -- naturally predicts the observed $\Sigma_{max}$. $\Sigma_{max}\sim 100\Sigma_{crit}$ because the gas consumption time is longer than the global freefall time even when feedback is ineffective. Moreover the predicted $\Sigma_{max}$ is robust to spatial scale and metallicity, and is preserved even if multiple episodes of star formation/gas inflow occur. In this context, the observed $\Sigma_{max}$ directly tells us where feedback fails.

We present a suite of 3D multi-physics MHD simulations following star formation in isolated turbulent molecular gas disks ranging from 5 to 500 parsecs in radius. These simulations are designed to survey the range of surface densities between those typical of Milky Way GMCs ($\sim 10^2 M_\odot\,pc^{-2}}$) and extreme ULIRG environments ($\sim 10^2 M_\odot\,pc^{-2}}$) so as to map out the scaling of the cloud-scale star formation efficiency (SFE) between these two regimes. The simulations include prescriptions for supernova, stellar wind, and radiative feedback, which we find to be essential in determining both the instantaneous per-freefall ($\epsilon_{ff}$) and integrated ($\epsilon_{int}$) star formation efficiencies. In all simulations, the gas disks form stars until a critical stellar surface density has been reached and the remaining gas is blown out by stellar feedback. We find that surface density is a good predictor of $\epsilon_{int}$, as suggested by analytic force balance arguments from previous works. SFE eventually saturates to $\sim 1$ at high surface density. We also find a proportional relationship between $\epsilon_{ff}$ and $\epsilon_{int}$, implying that star formation is feedback-moderated even over very short time-scales in isolated clouds. These results have implications for star formation in galactic disks, the nature and fate of nuclear starbursts, and the formation of bound star clusters. The scaling of $\epsilon_{ff}$ with surface density is not consistent with the notion that $\epsilon_{ff}$ is always $\sim 1\%$ on the scale of GMCs, but our predictions recover the $\sim 1\%$ value for GMC parameters similar to those found in sprial galaxies, including our own.

Using a state-of-the-art cosmological simulation of merging proto-galaxies at high redshift from the FIRE project, with explicit treatments of star formation and stellar feedback in the interstellar medium, we investigate the formation of star clusters and examine one of the formation hypothesis of present-day metal-poor globular clusters. We find that frequent mergers in high-redshift proto-galaxies could provide a fertile environment to produce long-lasting bound star clusters. The violent merger event disturbs the gravitational potential and pushes a large gas mass of ~> 1e5-6 Msun collectively to high density, at which point it rapidly turns into stars before stellar feedback can stop star formation. The high dynamic range of the reported simulation is critical in realizing such dense star-forming clouds with a small dynamical timescale, t_ff <~ 3 Myr, shorter than most stellar feedback timescales. Our simulation then allows us to trace how clusters could become virialized and tightly-bound to survive for up to ~420 Myr till the end of the simulation. Because the cluster's tightly-bound core was formed in one short burst, and the nearby older stars originally grouped with the cluster tend to be preferentially removed, at the end of the simulation the cluster has a small age spread.

We present high angular resolution multi-wavelength data of the 3C 298 radio-loud quasar host galaxy (z=1.439) taken using the W.M. Keck Observatory OSIRIS integral field spectrograph with adaptive optics, Atacama Large Millimeter/submillimeter Array (ALMA), Hubble Space Telescope (HST) WFC3, and the Very Large Array (VLA). Extended emission is detected in the rest-frame optical nebular emission lines H$\beta$, [OIII], H$\alpha$, [NII], and [SII], as well as molecular lines CO (J=3-2) and (J=5-4). Along the path of 3C 298's relativistic jets we detect conical outflows of ionized gas with velocities up to 1700 km s$^{-1}$ and outflow rate of 450-1500 M$_\odot$yr$^{-1}$. Near the spatial center of the conical outflow, CO (J=3-2) emission shows a molecular gas disc with a total molecular mass ($\rm M_{H_{2}}$) of 6.6$\pm0.36\times10^{9}$M$_{\odot}$. On the molecular disc's blueshifted side we observe a molecular outflow with a rate of 2300 M$_\odot$yr$^{-1}$ and depletion time scale of 3 Myr. We detect no narrow H$\alpha$ emission in the outflow regions, suggesting a limit on star formation of 0.3 M$_\odot$yr$^{-1}$kpc$^{-2}$. Quasar driven winds are evacuating the molecular gas reservoir thereby directly impacting star formation in the host galaxy. The observed mass of the supermassive black hole is $10^{9.37-9.56}$M$_{\odot}$ and we determine a dynamical bulge mass of 1-1.7$\rm\times10^{10}\frac{R}{1.6 kpc}$ M$_{\odot}$. The bulge mass of 3C 298 resides 2-2.5 orders of magnitude below the expected value from the local M$\rm_{bulge}-M_{BH}$ relationship. A second galactic disc observed in nebular emission is offset from the quasar by 9 kpc suggesting the system is an intermediate stage merger. These results show that galactic scale negative feedback is occurring early in the merger phase of 3C 298, well before the coalescence of the galactic nuclei and assembly on the local relationship.

The VIS instrument on board the Euclid mission is a weak-lensing experiment that depends on very precise shape measurements of distant galaxies obtained by a large CCD array. Due to the harsh radiative environment outside the Earth's atmosphere, it is anticipated that the CCDs over the mission lifetime will be degraded to an extent that these measurements will only be possible through the correction of radiation damage effects. We have therefore created a Monte Carlo model that simulates the physical processes taking place when transferring signal through a radiation-damaged CCD. The software is based on Shockley-Read-Hall theory, and is made to mimic the physical properties in the CCD as closely as possible. The code runs on a single electrode level and takes three dimensional trap position, potential structure of the pixel, and multi-level clocking into account. A key element of the model is that it also takes device specific simulations of electron density as a direct input, thereby avoiding to make any analytical assumptions about the size and density of the charge cloud. This paper illustrates how test data and simulated data can be compared in order to further our understanding of the positions and properties of the individual radiation-induced traps.

The light curve of GW170817 is surprisingly blue and bright. Assuming that the event is a binary neutron star merger, we argue that blueness and brightness of the light curve is the result of ejecta that contains an substantial amount of thermal energy. To achieve this, the ejecta must be reheated at a substantial distance (1 to 2000 solar radii) from the merger to avoid losing the energy to adiabatic cooling. We show that this reheating can occur if the merger occurs in a hierarchical triple system where the outer star has evolved and filled its Roche lobe. The outer star feeds mass to the inner binary, forming a circumbinary disc, driving the inner binary to merge. Because the outer star fills its Roche lobe, a substantial fraction of the dynamical ejecta collides with the evolved star, reheating the ejecta in the process. We suggest that the process of mass transfer in heirarchical triples tends to form coplanar triple systems such as PSR J0337+1715, and may provide electromagnetic counterparts to binary black hole mergers.

Following the first two annual intensity mapping workshops at Stanford in March 2016 and Johns Hopkins in June 2017, we report on the recent advances in theory, instrumentation and observation that were presented in these meetings and some of the opportunities and challenges that were identified looking forward. With preliminary detections of CO, [CII], Lya and low-redshift 21cm, and a host of experiments set to go online in the next few years, the field is rapidly progressing on all fronts, with great anticipation for a flood of new exciting results. This current snapshot provides an efficient reference for experts in related fields and a useful resource for nonspecialists. We begin by introducing the concept of line-intensity mapping and then discuss the broad array of science goals that will be enabled, ranging from the history of star formation, reionization and galaxy evolution to measuring baryon acoustic oscillations at high redshift and constraining theories of dark matter, modified gravity and dark energy. After reviewing the first detections reported to date, we survey the experimental landscape, presenting the parameters and capabilities of relevant instruments such as COMAP, mmIMe, AIM-CO, CCAT-p, TIME, CONCERTO, CHIME, HIRAX, HERA, STARFIRE, MeerKAT/SKA and SPHEREx. Finally, we describe recent theoretical advances: different approaches to modeling line luminosity functions, several techniques to separate the desired signal from foregrounds, statistical methods to analyze the data, and frameworks to generate realistic intensity map simulations.

We present molecular gas mass estimates for a sample of 13 local galaxies whose kinematic and star forming properties closely resemble those observed in $z\approx 1.5$ main-sequence galaxies. Plateau de Bure observations of the CO[1-0] emission line and Herschel Space Observatory observations of the dust emission both suggest molecular gas mass fractions of ~20%. Moreover, dust emission modeling finds $T_{dust}<$30K, suggesting a cold dust distribution compared to their high infrared luminosity. The gas mass estimates argue that $z\sim$0.1 DYNAMO galaxies not only share similar kinematic properties with high-z disks, but they are also similarly rich in molecular material. Pairing the gas mass fractions with existing kinematics reveals a linear relationship between $f_{gas}$ and $\sigma$/$v_{c}$, consistent with predictions from stability theory of a self-gravitating disk. It thus follows that high gas velocity dispersions are a natural consequence of large gas fractions. We also find that the systems with lowest depletion times ($\sim$0.5 Gyr) have the highest ratios of $\sigma$/$v_{c}$ and more pronounced clumps, even at the same high molecular gas fraction.

In the context of the ESA M5 (medium mission) call we proposed a new satellite mission, Theia, based on relative astrometry and extreme precision to study the motion of very faint objects in the Universe. Theia is primarily designed to study the local dark matter properties, the existence of Earth-like exoplanets in our nearest star systems and the physics of compact objects. Furthermore, about 15 $\%$ of the mission time was dedicated to an open observatory for the wider community to propose complementary science cases. With its unique metrology system and "point and stare" strategy, Theia's precision would have reached the sub micro-arcsecond level. This is about 1000 times better than ESA/Gaia's accuracy for the brightest objects and represents a factor 10-30 improvement for the faintest stars (depending on the exact observational program). In the version submitted to ESA, we proposed an optical (350-1000nm) on-axis TMA telescope. Due to ESA Technology readiness level, the camera's focal plane would have been made of CCD detectors but we anticipated an upgrade with CMOS detectors. Photometric measurements would have been performed during slew time and stabilisation phases needed for reaching the required astrometric precision.

We use cosmological hydrodynamic simulations with stellar feedback from the FIRE project to study the physical nature of Lyman limit systems (LLSs) at z<1. At these low redshifts, LLSs are closely associated with dense gas structures surrounding galaxies, such as galactic winds, dwarf satellites, and cool inflows from the intergalactic medium. Our analysis is based on 14 zoom-in simulations covering the halo mass range M_h~10^9-10^13 Msun at z=0, which we convolve with the dark matter halo mass function to produce cosmological statistics. We find that the majority of cosmologically-selected LLSs are associated with halos in the mass range 10^10 < M_h < 10^12 Msun. The incidence and HI column density distribution of simulated absorbers with columns 10^16.2 < N_HI < 2x10^20 cm^-2 are consistent with observations. High-velocity outflows (with radial velocity exceeding the halo circular velocity by a factor >~2) tend to have higher metallicities ([X/H] ~ -0.5) while very low metallicity ([X/H] < -2) LLSs are typically associated with gas infalling from the intergalactic medium. However, most LLSs occupy an intermediate region in metallicity-radial velocity space, for which there is no clear trend between metallicity and radial kinematics. Metal-enriched inflows arise in the FIRE simulations as a result of galactic winds that fall back onto galaxies at low redshift. The overall simulated LLS metallicity distribution has a mean (standard deviation) [X/H] = -0.9 (0.4) and does not show significant evidence for bimodality, in contrast to recent observational studies but consistent with LLSs arising from halos with a broad range of masses and metallicities.

TRAPPIST-1 is a late M-dwarf orbited by seven Earth-sized planets with orbital period ratios near a chain of mean motion resonances. Due to uncertain system parameters, most orbital configurations drawn from the inferred posterior distribution are unstable on short timescales of $\sim$ 0.5 Myr, even when including the eccentricity damping effect of tides. By contrast, we show that most physically plausible resonant configurations generated through disk migration are stable even without tidal dissipation on timescales of at least 50 Myr ($10^{10}$ orbits), an increase of at least two orders of magnitude. This result, together with the remarkable chain of period ratios in the system, provide strong evidence for convergent migration naturally emplacing the system near an equilibrium configuration forced by the resonant chain. We provide an openly available database of physically plausible initial conditions for TRAPPIST-1 generated through parametrized planet-disk interactions, as well as bit-by-bit reproducible N-body integrations over $10^9-10^{10}$ orbits.

We present an analysis of the flow of metals through the circumgalactic medium (CGM) in the Feedback in Realistic Environments (FIRE) simulations of galaxy formation, ranging from isolated dwarfs to $L*$ galaxies. We find that nearly all metals produced in high-redshift galaxies are carried out in winds that reach $0.25 R_{\rm vir}$. When measured at $0.25 R_{\rm vir}$ the metallicity of outflows is slightly higher than the interstellar medium (ISM) metallicity. Many metals thus reside in the CGM. Cooling and recycling from this reservoir determine the metal budget in the ISM. The outflowing metal flux decreases by a factor of $\sim2-5$ between $0.25 R_{\rm vir}$ and $R_{\rm vir}$. Furthermore, outflow metallicity is typically lower at $R_{\rm vir}$ owing to dilution of the remaining outflow by metal-poor material swept up from the CGM. The inflow metallicity at $R_{\rm vir}$ is generally low, but outflow and inflow metallicities are similar in the inner halo. At low redshift, massive galaxies no longer generate outflows that reach the CGM, causing a divergence in CGM and ISM metallicity. Dwarf galaxies continue to generate outflows, although they preferentially retain metal ejecta. In all but the least massive galaxy considered, a majority of the metals are within the halo at $z=0$. We measure the fraction of metals in CGM, ISM and stars, and quantify the thermal state of CGM metals in each halo. The total amount of metals in the low-redshift CGM of two simulated $L*$ galaxies is consistent with estimates from the COS halos survey, while for the other two it appears to be lower.

Efficient thermalization of overlapping supernovae within star-forming galaxies may produce a supernova-heated fluid that drives galactic winds. For fiducial assumptions about the timescale for cloud shredding from high-resolution simulations (which neglect magnetic fields) we show that cool clouds with temperature from $T_{c}\sim 10^{2}-10^{4}$ K seen in emission and absorption in galactic winds cannot be accelerated to observed velocities by the ram pressure of a hot wind. Taking into account both the radial structure of the hot flow and gravity, we show that this conclusion holds over a wide range of galaxy, cloud, and hot wind properties. This finding calls into question the prevailing picture whereby the cool atomic gas seen in galactic winds is entrained and accelerated by the hot flow. Given these difficulties with ram pressure acceleration, we discuss alternative models for the origin of high velocity cool gas outflows. Another possibility is that magnetic fields in cool clouds are sufficiently important that they prolong the cloud's life. For $T_{c}=10^{3}$\,K and $10^{4}$\,K clouds, we show that if conductive evaporation can be neglected, the cloud shredding timescale must be $\sim15$ and 5 times longer, respectively, than the values from hydrodynamical simulations in order for cool cloud velocities to reach those seen in observations.

The Atacama Large Millimeter/submilimeter Array (ALMA) recently revealed a set of nearly concentric gaps in the protoplanetary disk surrounding the young star HL Tau. If these are carved by forming gas giants, this provides the first set of orbital initial conditions for planets as they emerge from their birth disks. Using N-body integrations, we have followed the evolution of the system for 5 Gyr to explore the possible outcomes. We find that HL Tau initial conditions scaled down to the size of typically observed exoplanet orbits naturally produce several populations in the observed exoplanet sample. First, for a plausible range of planetary masses, we can match the observed eccentricity distribution of dynamically excited radial velocity giant planets with eccentricities $>$ 0.2. Second, we roughly obtain the observed rate of hot Jupiters around FGK stars. Finally, we obtain a large efficiency of planetary ejections of $\approx 2$ per HL Tau-like system, but the small fraction of stars observed to host giant planets makes it hard to match the rate of free-floating planets inferred from microlensing observations. In view of upcoming GAIA results, we also provide predictions for the expected mutual inclination distribution, which is significantly broader than the absolute inclination distributions typically considered by previous studies.

We present results from IROCKS (Intermediate Redshift OSIRIS Chemo-Kinematic Survey) for sixteen z~1 and one z~1.4 star-forming galaxies. All galaxies were observed with OSIRIS with the laser guide star adaptive optics system at Keck Observatory. We use rest-frame nebular Ha emission lines to trace morphologies and kinematics of ionized gas in star-forming galaxies on sub-kiloparsec physical scales. We observe elevated velocity dispersions (sigma > 50 km/s) seen in z > 1.5 galaxies persist at z~1 in the integrated galaxies. Using an inclined disk model and the ratio of v/sigma, we find that 1/3 of the z~1 sample are disk candidates while the other 2/3 of the sample are dominated by merger-like and irregular sources. We find that including extra attenuation towards HII regions derived from stellar population synthesis modeling brings star formation rates (SFR) using Ha and stellar population fit into a better agreement. We explore properties of compact Ha sub-component, or "clump," at z~1 and find that they follow a similar size-luminosity relation as local HII regions but are scaled-up by an order of magnitude with higher luminosities and sizes. Comparing the z~1 clumps to other high-redshift clump studies, we determine that the clump SFR surface density evolves as a function of redshift. This may imply clump formation is directly related to the gas fraction in these systems and support disk fragmentation as their formation mechanism since gas fraction scales with redshift.

This study presents a catalog of 8107 molecular clouds that covers the entire Galactic plane and includes 98% of the $^{12}$CO emission observed within $b\pm5^\circ$. The catalog was produced using a hierarchical cluster identification method applied to the result of a Gaussian decomposition of the Dame et al. data. The total H$_2$ mass in the catalog is $1.2\times10^9 M_\odot$, in agreement with previous estimates. We find that 30% of the sight lines intersect only a single cloud, with another 25% intersecting only two clouds. The most probable cloud size is $R\sim30$ pc. We find that $M\propto R^{2.2\pm0.2}$, with no correlation between the cloud surface density, $\Sigma$, and $R$. In contrast with the general idea, we find a rather large range of values of $\Sigma$, from 2 to $300 M_\odot$ pc$^{-2}$, and a systematic decrease with increasing Galactic radius, $R_{\rm gal}$. The cloud velocity dispersion and the normalization $\sigma_0=\sigma_v/R^{1/2}$ both decrease systematically with $R_{\rm gal}$. When studied over the whole Galactic disk, there is a large dispersion in the line width-size relation, and a significantly better correlation between $\sigma_v$ and $\Sigma\,R$. The normalization of this correlation is constant to better than a factor of two for $R_{\rm gal}<20$ kpc. This relation is used to disentangle the ambiguity between near and far kinematic distances. We report a strong variation of the turbulent energy injection rate. In the outer Galaxy it may be maintained by accretion through the disk and/or onto the clouds, but neither source can drive the 100 times higher cloud-averaged injection rate in the inner Galaxy.

The requirement that planetary systems be dynamically stable is often used to vet new discoveries or set limits on unconstrained masses or orbital elements. This is typically carried out via computationally expensive N-body simulations. We show that characterizing the complicated and multi-dimensional stability boundary of tightly packed systems is amenable to machine learning methods. We find that training an XGBoost machine learning algorithm on physically motivated features yields an accurate classifier of stability in packed systems. On the stability timescale investigated ($10^7$ orbits), it is 3 orders of magnitude faster than direct N-body simulations. Optimized machine learning classifiers for dynamical stability may thus prove useful across the discipline, e.g., to characterize the exoplanet sample discovered by the upcoming Transiting Exoplanet Survey Satellite (TESS). This proof of concept motivates investing computational resources to train algorithms capable of predicting stability over longer timescales and over broader regions of phase space.

The morphology of massive star-forming galaxies at high redshift is often dominated by giant clumps of mass ~10^8-10^9 Msun and size ~100-1000 pc. Previous studies have proposed that giant clumps might have an important role in the evolution of their host galaxy, particularly in building the central bulge. However, this depends on whether clumps live long enough to migrate from their original location in the disc or whether they get disrupted by their own stellar feedback before reaching the centre of the galaxy. We use cosmological hydrodynamical simulations from the FIRE (Feedback in Realistic Environments) project that implement explicit treatments of stellar feedback and ISM physics to study the properties of these clumps. We follow the evolution of giant clumps in a massive (stellar mass ~10^10.8 Msun at z=1), discy, gas-rich galaxy from redshift z>2 to z=1. Even though the clumpy phase of this galaxy lasts over a gigayear, individual gas clumps are short-lived, with mean lifetime of massive clumps of ~20 Myr. During that time, they turn between 0.1% and 20% of their gas into stars before being disrupted, similar to local GMCs. Clumps with M>10^7 Msun account for ~20% of the total star formation in the galaxy during the clumpy phase, producing ~10^10 Msun of stars. We do not find evidence for net inward migration of clumps within the galaxy. The number of giant clumps and their mass decrease at lower redshifts, following the decrease in the overall gas fraction and star-formation rate.

Motivated by the nonlinear star formation efficiency found in recent numerical simulations by a number of workers, we perform high-resolution adaptive mesh refinement simulations of star formation in self-gravitating turbulently driven gas. As we follow the collapse of this gas, we find that the character of the flow changes at two radii, the disk radius $r_d$, and the radius $r_*$ where the enclosed gas mass exceeds the stellar mass. Accretion starts at large scales and works inwards. In line with recent analytical work, we find that the density evolves to a fixed attractor, $\rho(r,t ) \rightarrow \rho(r)$, for $r_d<r<r_*$; mass flows through this structure onto a sporadically gravitationally unstable disk, and from thence onto the star. In the bulk of the simulation box we find that the random motions $v_T \sim r^p$ with $p \sim 0.5$, in agreement with Larson's size-linewidth relation. In the vicinity of massive star forming regions we find $ p \sim 0.2-0.3$, as seen in observations. For $r<r_*$, $v_T$ increases inward, with $p=-1/2$. Finally, we find that the total stellar mass $M_*(t)\sim t^2$ in line with previous numerical and analytic work that suggests a nonlinear rate of star formation.

We use cosmological simulations from the FIRE (Feedback In Realistic Environments) project to study the baryon cycle and galaxy mass assembly for central galaxies in the halo mass range $M_{\rm halo} \sim 10^{10} - 10^{13} M_{\odot}$. By tracing cosmic inflows, galactic outflows, gas recycling, and merger histories, we quantify the contribution of physically distinct sources of material to galaxy growth. We show that in situ star formation fueled by fresh accretion dominates the early growth of galaxies of all masses, while the re-accretion of gas previously ejected in galactic winds often dominates the gas supply for a large portion of every galaxy's evolution. Externally processed material contributes increasingly to the growth of central galaxies at lower redshifts. This includes stars formed ex situ and gas delivered by mergers, as well as smooth intergalactic transfer of gas from other galaxies, an important but previously under-appreciated growth mode. By $z=0$, wind transfer, i.e. the exchange of gas between galaxies via winds, can dominate gas accretion onto $\sim L^{*}$ galaxies over fresh accretion and standard wind recycling. Galaxies of all masses re-accrete >50% of the gas ejected in winds and recurrent recycling is common. The total mass deposited in the intergalactic medium per unit stellar mass formed increases in lower mass galaxies. Re-accretion of wind ejecta occurs over a broad range of timescales, with median recycling times ($\sim 100-350$ Myr) shorter than previously found. Wind recycling typically occurs at the scale radius of the halo, independent of halo mass and redshift, suggesting a characteristic recycling zone around galaxies that scales with the size of the inner halo and the galaxy's stellar component.

Due to its proximity, the Milky Way nuclear star cluster provides us with a wealth of data not available in other galactic nuclei. In particular, with adaptive optics, we can observe the detailed properties of individual stars, which can offer insight into the origin and evolution of the cluster. We summarize work on the central parsec of the Galactic center based on imaging and spectroscopic observations at the Keck and Gemini telescopes. These observations include stellar positions in two dimension and the velocity in three dimensions. Spectroscopic observations also enable measurements of the physical properties of individual stars, such as the spectral type and in some cases the effective temperature, metallicity, and surface gravity. We present a review of our latest measurements of the density profiles and luminosity functions of the young and old stars in this region. These observations show a complex stellar population with a young (4-6 Myr) compact star cluster in the central 0.5 pc embedded in an older and much more massive nuclear star cluster. Surprisingly, the old late-type giants do not show a cusp profile as long been expected from theoretical work. The solution to the missing cusp problem may offer us insight into the dynamical evolution of the cluster. Finally, we also discuss recent work on the metallicity of stars in this region and how they might be used to trace their origin. The nuclear star cluster shows a large variation in metallicity ([M/Fe]). The majority of the stars have higher than solar metallicity, with about 6\% having [M/Fe] $< -0.5$. These observations indicate that the NSC was not built from the globular clusters that we see today. The formation of the nuclear star cluster is more likely from the inward migration of gas originating in the disk of the Milky Way.

We present ALMA CO(1-0) and CO(3-2) observations of the brightest cluster galaxy (BCG) in the 2A 0335+096 galaxy cluster (z = 0.0346). The total molecular gas mass of (1.13+/-0.15) x 10^9 M_sun is divided into two components: a nuclear region and a 7 kpc long dusty filament. The central molecular gas component accounts for (3.2+/-0.4) x 10^8 M_sun of the total supply of cold gas. Instead of forming a rotationally-supported ring or disk, it is composed of two distinct, blueshifted clumps south of the nucleus and a series of low-significance redshifted clumps extending toward a nearby companion galaxy. The velocity of the redshifted clouds increases with radius to a value consistent with the companion galaxy, suggesting that an interaction between these galaxies <20 Myr ago disrupted a pre-existing molecular gas reservoir within the BCG. Most of the molecular gas, (7.8+/-0.9) x 10^8 M_sun, is located in the filament. The CO emission is co-spatial with a 10^4 K emission-line nebula and soft X-rays from 0.5 keV gas, indicating that the molecular gas has cooled out of the intracluster medium over a period of 25-100 Myr. The filament trails an X-ray cavity, suggesting that the gas has cooled from low entropy gas that has been lifted out of the cluster core and become thermally unstable. We are unable to distinguish between inflow and outflow along the filament with the present data. Cloud velocities along the filament are consistent with gravitational free-fall near the plane of the sky, although their increasing blueshifts with radius are consistent with outflow.

Star formation on galactic scales is known to be a slow process, but whether it is slow on smaller scales is uncertain. We cross-correlate 5469 giant molecular clouds (GMCs) from a new all-sky catalog with 256 star forming complexes (SFCs) to build a sample of 191 SFC-GMC complexes---collections of multiple clouds each matched to 191 SFCs. The total mass in stars harbored by these clouds is inferred from WMAP free-free fluxes. We measure the GMC mass, the virial parameter, the star formation efficiency $\epsilon$ and the star formation rate per free-fall time $\epsilon_{\rm ff}$. Both $\epsilon$ and $\epsilon_{\rm ff}$ range over 3--4 orders of magnitude. We find that 68.3% of the clouds fall within $\sigma_{\log\epsilon}=0.79\pm0.22\,{\rm dex}$ and $\sigma_{\log\epsilon_{\rm ff}}=0.91\pm0.22\,{\rm dex}$ about the median. Compared to these observed scatters, a simple model with a time independent $\epsilon_{\rm ff}$ that depends on the host GMC properties predicts $\sigma_{\log\epsilon_{\rm ff}}=$0.12-0.24. Allowing for a time-variable $\epsilon_{\rm ff}$, we can recover the large dispersion in the rate of star formation. This strongly suggests that star formation in the Milky Way is a dynamic process on GMC scales. We also show that the surface star formation rate profile of the Milky Way correlates well with the molecular gas surface density profile.

Euclid-VIS is the large format visible imager for the ESA Euclid space mission in their Cosmic Vision program, scheduled for launch in 2020. Together with the near infrared imaging within the NISP instrument, it forms the basis of the weak lensing measurements of Euclid. VIS will image in a single r+i+z band from 550-900 nm over a field of view of ~0.5 deg2. By combining 4 exposures with a total of 2260 sec, VIS will reach to deeper than mAB=24.5 (10sigma) for sources with extent ~0.3 arcsec. The image sampling is 0.1 arcsec. VIS will provide deep imaging with a tightly controlled and stable point spread function (PSF) over a wide survey area of 15000 deg2 to measure the cosmic shear from nearly 1.5 billion galaxies to high levels of accuracy, from which the cosmological parameters will be measured. In addition, VIS will also provide a legacy dataset with an unprecedented combination of spatial resolution, depth and area covering most of the extra-Galactic sky. Here we will present the results of the study carried out by the Euclid Consortium during the period up to the Critical Design Review.

We use cosmological simulations from the Feedback In Realistic Environments (FIRE) project, which implement a comprehensive set of stellar feedback processes, to study ultra-violet (UV) metal line emission from the circum-galactic medium of high-redshift (z=2-4) galaxies. Our simulations cover the halo mass range Mh ~ 2x10^11 - 8.5x10^12 Msun at z=2, representative of Lyman break galaxies. Of the transitions we analyze, the low-ionization C III (977 A) and Si III (1207 A) emission lines are the most luminous, with C IV (1548 A) and Si IV (1394 A) also showing interesting spatially-extended structures. The more massive halos are on average more UV-luminous. The UV metal line emission from galactic halos in our simulations arises primarily from collisionally ionized gas and is strongly time variable, with peak-to-trough variations of up to ~2 dex. The peaks of UV metal line luminosity correspond closely to massive and energetic mass outflow events, which follow bursts of star formation and inject sufficient energy into galactic halos to power the metal line emission. The strong time variability implies that even some relatively low-mass halos may be detectable. Conversely, flux-limited samples will be biased toward halos whose central galaxy has recently experienced a strong burst of star formation. Spatially-extended UV metal line emission around high-redshift galaxies should be detectable by current and upcoming integral field spectrographs such as the Multi Unit Spectroscopic Explorer (MUSE) on the Very Large Telescope and Keck Cosmic Web Imager (KCWI).

Direct detection of black hole-neutron star pairs is anticipated with the advent of aLIGO. Electromagnetic counterparts may be crucial for a confident gravitational-wave detection as well as for extraction of astronomical information. Yet black hole-neutron star pairs are notoriously dark and so inaccessible to telescopes. Contrary to this expectation, a bright electromagnetic transient can occur in the final moments before merger as long as the neutron star is highly magnetized. The orbital motion of the neutron star magnet creates a Faraday flux and corresponding power available for luminosity. A spectrum of curvature radiation ramps up until the rapid injection of energy ignites a fireball, which would appear as an energetic blackbody peaking in the x ray to $\gamma$ rays for neutron star field strengths ranging from $10^{12}$G to $10^{16}$G respectively and a $10M_{\odot} $ black hole. The fireball event may last from a few milliseconds to a few seconds depending on the neutron star magnetic-field strength, and may be observable with the Fermi Gamma-Ray Burst Monitor with a rate up to a few per year for neutron star field strengths $\gtrsim 10^{14}$G. We also discuss a possible decaying post-merger event which could accompany this signal. As an electromagnetic counterpart to these otherwise dark pairs, the black-hole battery should be of great value to the development of multi-messenger astronomy in the era of aLIGO.

Observations reveal that quasar host halos at z~2 have large covering fractions of cool dense gas (>~60% for Lyman limit systems within a projected virial radius). Most simulations have so far have failed to explain these large observed covering fractions. We analyze a new set of 15 simulated massive halos with explicit stellar feedback from the FIRE project, covering the halo mass range M_h~2x10^12-10^13 Msun at z=2. This extends our previous analysis of the circum-galactic medium of high-redshift galaxies to more massive halos. AGN feedback is not included in these simulations. We find Lyman limit system covering fractions consistent with those observed around quasars. The large HI covering fractions arise from star formation-driven galactic winds, including winds from low-mass satellite galaxies that interact with cosmological filaments. We show that it is necessary to resolve these satellite galaxies and their winds to reproduce the large Lyman limit system covering fractions observed in quasar-mass halos. Our simulations predict that galaxies occupying dark matter halos of mass similar to quasars but without a luminous AGN should have Lyman limit system covering fractions comparable to quasars.

Empirical constraints on reionization require galactic ionizing photon escape fractions fesc>20%, but recent high-resolution radiation-hydrodynamic calculations have consistently found much lower values ~1-5%. While these models include strong stellar feedback and additional processes such as runaway stars, they almost exclusively consider stellar evolution models based on single (isolated) stars, despite the fact that most massive stars are in binaries. We re-visit these calculations, combining radiative transfer and high-resolution cosmological simulations from the Feedback in Realistic Environments (FIRE) project. For the first time, we use a stellar evolution model that includes a physically and observationally motivated treatment of binaries (the BPASS model). Binary mass transfer and mergers enhance the population of massive stars at late times (>3 Myr) after star formation, which in turn strongly enhances the late-time ionizing photon production (especially at low metallicities). These photons are produced after feedback from massive stars has carved escape channels in the ISM, and so efficiently leak out of galaxies. As a result, the time-averaged "effective" escape fraction (ratio of escaped ionizing photons to observed 1500 A photons) increases by factors 4-10, sufficient to explain reionization. While important uncertainties remain, we conclude that binary evolution may be critical for understanding the ionization of the Universe.

We present ALMA observations of the CO(1-0) and CO(3-2) line emission tracing filaments of cold molecular gas in the central galaxy of the cluster PKS0745-191. The total molecular gas mass of 4.6 +/- 0.3 x 10^9 solar masses, assuming a Galactic X_{CO} factor, is divided roughly equally between three filaments each extending radially 3-5 kpc from the galaxy centre. The emission peak is located in the SE filament roughly 1 arcsec (2 kpc) from the nucleus. The velocities of the molecular clouds in the filaments are low, lying within +/-100 km/s of the galaxy's systemic velocity. Their FWHMs are less than 150 km/s, which is significantly below the stellar velocity dispersion. Although the molecular mass of each filament is comparable to a rich spiral galaxy, such low velocities show that the filaments are transient and the clouds would disperse on <10^7 yr timescales unless supported, likely by the indirect effect of magnetic fields. The velocity structure is inconsistent with a merger origin or gravitational free-fall of cooling gas in this massive central galaxy. If the molecular clouds originated in gas cooling even a few kpc from their current locations their velocities would exceed those observed. Instead, the projection of the N and SE filaments underneath X-ray cavities suggests they formed in the updraft behind bubbles buoyantly rising through the cluster atmosphere. Direct uplift of the dense gas by the radio bubbles appears to require an implausibly high coupling efficiency. The filaments are coincident with low temperature X-ray gas, bright optical line emission and dust lanes indicating that the molecular gas could have formed from lifted warmer gas that cooled in situ.

We study the interaction of feedback from active galactic nuclei (AGN) and a multi-phase interstellar medium (ISM), in simulations including explicit stellar feedback, multi-phase cooling, accretion-disk winds, and Compton heating. We examine radii ~0.1-100 pc around a black hole (BH), where the accretion rate onto the BH is determined and where AGN-powered winds and radiation couple to the ISM. We conclude: (1) The BH accretion rate is determined by exchange of angular momentum between gas and stars in gravitational instabilities. This produces accretion rates ~0.03-1 Msun/yr, sufficient to power luminous AGN. (2) The gas disk in the galactic nucleus undergoes an initial burst of star formation followed by several Myrs where stellar feedback suppresses the star formation rate (SFR). (3) AGN winds injected at small radii with momentum fluxes ~L/c couple efficiently to the ISM and have dramatic effects on ISM properties within ~100 pc. AGN winds suppress the nuclear SFR by factors ~10-30 and BH accretion rate by factors ~3-30. They increase the outflow rate from the nucleus by factors ~10, consistent with observational evidence for galaxy-scale AGN-driven outflows. (4) With AGN feedback, the predicted column density distribution to the BH is consistent with observations. Absent AGN feedback, the BH is isotropically obscured and there are not enough optically-thin sightlines to explain Type-I AGN. A 'torus-like' geometry arises self-consistently as AGN feedback evacuates gas in polar regions.

We examine the stability of feedback-regulated star formation (SF) in galactic nuclei and contrast it to SF in extended discs. In galactic nuclei the dynamical time becomes shorter than the time over which feedback from young stars evolves. We argue analytically that the balance between stellar feedback and gravity is unstable in this regime. We study this using numerical simulations with pc-scale resolution and explicit stellar feedback taken from stellar evolution models. The nuclear gas mass, young stellar mass, and SFR within the central ~100 pc (the short-timescale regime) never reach steady-state, but instead go through dramatic, oscillatory cycles. Stars form until a critical surface density of young stars is present (such that feedback overwhelms gravity), at which point they begin to expel gas from the nucleus. Since the dynamical times are shorter than the stellar evolution times, the stars do not die as the gas is expelled, but continue to push, triggering a runaway quenching of star formation in the nucleus. However the expelled gas is largely not unbound from the galaxy, but goes into a galactic fountain which re-fills the nuclear region after the massive stars from the previous burst cycle have died off (~50 Myr timescale). On large scales (>1 kpc), the galaxy-scale gas content and SFR is more stable. We examine the consequences of this episodic nuclear star formation for the Kennicutt-Schmidt (KS) relation: while a tight KS relation exists on ~1 kpc scales in good agreement with observations, the scatter increases dramatically in smaller apertures centered on galactic nuclei.

The Wide Field InfraRed Survey Telescope-Astrophysics Focused Telescope Asset (WFIRST-AFTA) mission is a 2.4-m class space telescope that will be used across a swath of astrophysical research domains. JPL will provide a high-contrast imaging coronagraph instrument - one of two major astronomical instruments. In order to achieve the low noise performance required to detect planets under extremely low flux conditions, the electron multiplying charge-coupled device (EMCCD) has been baselined for both of the coronagraph's sensors - the imaging camera and integral field spectrograph. JPL has established an EMCCD test laboratory in order to advance EMCCD maturity to technology readiness level-6. This plan incorporates full sensor characterization, including read noise, dark current, and clock-induced charge. In addition, by considering the unique challenges of the WFIRST space environment, degradation to the sensor's charge transfer efficiency will be assessed, as a result of damage from high-energy particles such as protons, electrons, and cosmic rays. Science-grade CCD201-20 EMCCDs have been irradiated to a proton fluence that reflects the projected WFIRST orbit. Performance degradation due to radiation displacement damage is reported, which is the first such study for a CCD201-20 that replicates the WFIRST conditions. In addition, techniques intended to identify and mitigate radiation-induced electron trapping, such as trap pumping, custom clocking, and thermal cycling, are discussed.

We study the abundance of substructure in the matter density near galaxies using ALMA Science Verification observations of the strong lensing system SDP.81. We present a method to measure the abundance of subhalos around galaxies using interferometric observations of gravitational lenses. Using simulated ALMA observations, we explore the effects of various systematics, including antenna phase errors and source priors, and show how such errors may be measured or marginalized. We apply our formalism to ALMA observations of SDP.81. We find evidence for the presence of a $M=10^{8.96\pm 0.12} M_{\odot}$ subhalo near one of the images, with a significance of $6.9\sigma$ in a joint fit to data from bands 6 and 7; the effect of the subhalo is also detected in both bands individually. We also derive constraints on the abundance of dark matter subhalos down to $M\sim 2\times 10^7 M_{\odot}$, pushing down to the mass regime of the smallest detected satellites in the Local Group, where there are significant discrepancies between the observed population of luminous galaxies and predicted dark matter subhalos. We find hints of additional substructure, warranting further study using the full SDP.81 dataset (including, for example, the spectroscopic imaging of the lensed carbon monoxide emission). We compare the results of this search to the predictions of $\Lambda$CDM halos, and find that given current uncertainties in the host halo properties of SDP.81, our measurements of substructure are consistent with theoretical expectations. Observations of larger samples of gravitational lenses with ALMA should be able to improve the constraints on the abundance of galactic substructure.

We present an analysis of the galaxy-scale gaseous outflows from the FIRE (Feedback in Realistic Environments) simulations. This suite of hydrodynamic cosmological zoom simulations resolves formation of star-forming giant molecular clouds to $z=0$, and features an explicit stellar feedback model on small scales. Our simulations reveal that high redshift galaxies undergo bursts of star formation followed by powerful gusts of galactic outflows that eject much of the ISM and temporarily suppress star formation. At low redshift, however, sufficiently massive galaxies corresponding to L*-progenitors develop stable disks and switch into a continuous and quiescent mode of star formation that does not drive outflows far into the halo. Mass-loading factors for winds in L*-progenitors are $\eta \approx 10$ at high redshift, but decrease to $\eta \ll 1$ at low redshift. Although lower values of $\eta$ are expected as halos grow in mass over time, we show that the strong suppression of outflows with decreasing redshift cannot be explained by mass evolution alone. Circumgalactic outflow velocities are variable and broadly distributed, but typically range between one and three times the circular velocity of the halo. Much of the ejected material builds a reservoir of enriched gas within the circumgalactic medium, some of which could be later recycled to fuel further star formation. However, a fraction of the gas that leaves the virial radius through galactic winds is never regained, causing most halos with mass $M_h \le 10^{12} M_{\odot}$ to be deficient in baryons compared to the cosmic mean by $z=0$.

We present ultra-high resolution cosmological hydrodynamic simulations of $M_*\simeq10^{4-6.3}M_{\odot}$ dwarf galaxies that form within $M_{v}=10^{9.5-10}M_{\odot}$ dark matter halos. Our simulations rely on the FIRE implementation of star formation feedback and were run with high enough force and mass resolution to directly resolve stellar and dark matter structure on the ~200 pc scales of interest for classical and ultra-faint dwarfs in the Local Group. The resultant galaxies sit on the $M_*$ vs. $M_{v}$ relation required to match the Local Group stellar mass function. They have bursty star formation histories and also form with half-light radii and metallicities that broadly match those observed for local dwarfs at the same stellar mass. We demonstrate that it is possible to create a large (~1 kpc) dark matter core in a cosmological simulation of an $M_*\simeq10^{6.5}M_{\odot}$ dwarf galaxy that resides within an $M_{v}=10^{10}M_{\odot}$ halo -- precisely the scale of interest for resolving the Too Big to Fail problem. However, these large cores are not ubiquitous and appear to correlate closely with the star formation histories of the dwarfs: dark matter cores are largest in systems that form their stars late ($z\lesssim2$), after the early epoch of cusp building mergers has ended. Our $M_*\simeq10^4M_{\odot}$ dwarf retains a cuspy dark matter halo density profile that matches almost identically that of a dark-matter only run of the same system. Despite forming in a field environment, this very low mass dwarf has observable properties that match closely to those of ultra-faint satellite galaxies of the Milky Way, including a uniformly old stellar population (>10 Gyr). Though ancient, most of the stars in our ultra-faint form after reionization; the UV field acts mainly to suppress fresh gas accretion, not to boil away gas that is already present in the proto-dwarf.

Was PS1-10jh (Gezari et al. 2012), an optical/UV transient discovered by the Pan-STARRS Medium Deep Survey, the tidal disruption of a star by a massive black hole? We address two aspects of the problem: the composition of the putative disrupted object (using the spectroscopic data), and the energetics of the observed gas and radiation (using the photometric data). We perform photoionization calculations and compare with the observed lower limit of the line ratio L(He II 4686)/L(Halpha) > 5 to argue that this event was not the disruption of a solar-type star, and instead was likely the disruption of a helium core (as first proposed by Gezari et al. 2012). Disruption of such a dense object requires a relatively small central BH, M_BH <~ 2 x 10^5 M_sun. We use the photometric data to infer that PS1-10jh comprised an outflow of ~ 0.01 M_sun of gas, escaping from the BH at ~1000 km/s, and we propose that this outflow was driven primarily by radiation pressure trapped by Thomson and resonance line scattering. The large ratio of radiated energy to kinetic energy, E_rad/E_K ~ 10^4, together with the large value of E_rad ~ 2 x 10^(51) erg, suggests that the outflow was shocked at large radius (perhaps similar to super-luminous supernovae or the internal shock model for gamma-ray bursts). We describe puzzles in the physics of PS1-10jh, and discuss how this event may help us understand future tidal disruptions and super-Eddington accretion events as well.

We present a series of high-resolution (20-2000 Msun, 0.1-4 pc) cosmological zoom-in simulations at z~6 from the Feedback In Realistic Environment (FIRE) project. These simulations cover halo masses 10^9-10^11 Msun and rest-frame ultraviolet magnitude Muv = -9 to -19. These simulations include explicit models of the multi-phase ISM, star formation, and stellar feedback, which produce reasonable galaxy properties at z = 0-6. We post-process the snapshots with a radiative transfer code to evaluate the escape fraction (fesc) of hydrogen ionizing photons. We find that the instantaneous fesc has large time variability (0.01%-20%), while the time-averaged fesc over long time-scales generally remains ~5%, considerably lower than the estimate in many reionization models. We find no strong dependence of fesc on galaxy mass or redshift. In our simulations, the intrinsic ionizing photon budgets are dominated by stellar populations younger than 3 Myr, which tend to be buried in dense birth clouds. The escaping photons mostly come from populations between 3-10 Myr, whose birth clouds have been largely cleared by stellar feedback. However, these populations only contribute a small fraction of intrinsic ionizing photon budgets according to standard stellar population models. We show that fesc can be boosted to high values, if stellar populations older than 3 Myr produce more ionizing photons than standard stellar population models (as motivated by, e.g., models including binaries). By contrast, runaway stars with velocities suggested by observations can enhance fesc by only a small fraction. We show that "sub-grid" star formation models, which do not explicitly resolve star formation in dense clouds with n >> 1 cm^-3, will dramatically over-predict fesc.

We present new near-infrared spectroscopic observations of the outer edges of the young stellar cluster around the supermassive black hole at the Galactic center. The observations show a break in the surface-density profile of young stars at approximately 13 arcsec (0.52 pc). These observations spectroscopically confirm previous suggestions of a break based on photometry. Using Gemini North's Near-Infrared Integral Field Spectrometer (NIFS) we are able to detect and separate early- and late-type stars with a 75% completeness at Ks = 15.5. We sample a region with radii between 7" to 23" (0.28 pc to 0.92 pc) from Sgr A*, and present new spectral classifications of 144 stars brighter than Ks = 15.5, where 140 stars are late-type (> 1 Gyr) and only four stars are early-type (young, 4-6 Myr). A broken power-law fit of the early-type surface-density matches well with our data and previously published values. The projected surface-density of late-type stars is also measured and found to be consistent with previous results. We find that the observed early-type surface-density profile is inconsistent with the theory of the young stars originating from a tightly bound infalling cluster, as no significant trail of young stars is found at radii above 13". We also note that either a simple disk instability criterion or a cloud-cloud collision could explain the location of the outer edge, though we lack information to make conclusive remarks on either alternative. If this break in surface-density represents an edge to the young stellar cluster it would set an important scale for the most recent episode of star formation at the Galactic center.

Planets in the habitable zone of lower-mass stars are often assumed to be in a state of tidally synchronized rotation, which would considerably affect their putative habitability. Although thermal tides cause Venus to rotate retrogradely, simple scaling arguments tend to attribute this peculiarity to the massive Venusian atmosphere. Using a global climate model, we show that even a relatively thin atmosphere can drive terrestrial planets' rotation away from synchronicity. We derive a more realistic atmospheric tide model that predicts four asynchronous equilibrium spin states, two being stable, when the amplitude of the thermal tide exceeds a threshold that is met for habitable Earth-like planets with a 1-bar atmosphere around stars more massive than 0.5-0.7Msun. Thus, many recently discovered terrestrial planets could exhibit asynchronous spin-orbit rotation, even with a thin atmosphere.

We use high-resolution cosmological zoom-in simulations from the FIRE project to make predictions for the covering fractions of neutral hydrogen around galaxies at z=2-4. These simulations resolve the interstellar medium of galaxies and explicitly implement a comprehensive set of stellar feedback mechanisms. Our simulation sample consists of 16 main halos covering the mass range M_h~10^9-6x10^12 Msun at z=2, including 12 halos in the mass range M_h~10^11-10^12 Msun corresponding to Lyman break galaxies (LBGs). We process our simulations with a ray tracing method to compute the ionization state of the gas. Galactic winds increase the HI covering fractions in galaxy halos by direct ejection of cool gas from galaxies and through interactions with gas inflowing from the intergalactic medium. Our simulations predict HI covering fractions for Lyman limit systems (LLSs) consistent with measurements around z~2-2.5 LBGs; these covering fractions are a factor ~2 higher than our previous calculations without galactic winds. The fractions of HI absorbers arising in inflows and in outflows are on average ~50% but exhibit significant time variability, ranging from ~10% to ~90%. For our most massive halos, we find a factor ~3 deficit in the LLS covering fraction relative to what is measured around quasars at z~2, suggesting that the presence of a quasar may affect the properties of halo gas on ~100 kpc scales. The predicted covering fractions, which decrease with time, peak at M_h~10^11-10^12 Msun, near the peak of the star formation efficiency in dark matter halos. In our simulations, star formation and galactic outflows are highly time dependent; HI covering fractions are also time variable but less so because they represent averages over large areas.

It is typically assumed that radiation pressure driven winds are accelerated to an asymptotic velocity of V ~ v_esc, where v_esc is the escape velocity from the central source. We note that this is not the case for dusty shells and clouds. Instead, if the shell or cloud is initially optically-thick to the UV emission from the source of luminosity L, then there is a significant boost in V that reflects the integral of the momentum absorbed as it is accelerated. For shells reaching a generalized Eddington limit, we show that V ~ (4R_UV L/M_sh c)^1/2, in both point-mass and isothermal-sphere potentials, where R_UV is the radius where the shell becomes optically-thin to UV photons, and M_sh is the mass of the shell. The asymptotic velocity significantly exceeds v_esc for typical parameters, and can explain the ~1000-2000km/s outflows observed from rapidly star-forming galaxies and active galactic nuclei if the surrounding halo has low gas density. Similarly fast outflows from massive stars can be accelerated on few - 10^3 yr timescales. These results carry over to clouds that subtend only a small fraction of the solid angle from the source of radiation and that expand as a consequence of their internal sound speed. We further consider the dynamics of shells that sweep up a dense circumstellar or circumgalactic medium. We calculate the "momentum ratio" Mdot v/(L/c) in the shell limit and show that it can only significantly exceed ~2 if the effective optical depth of the shell to re-radiated FIR photons is much larger than unity. We discuss simple prescriptions for the properties of galactic outflows for use in large-scale cosmological simulations. We also briefly discuss applications to the dusty ejection episodes of massive stars, the disruption of giant molecular clouds, and AGN.

We present numerical evidence of dynamic star formation in which the accreted stellar mass grows superlinearly with time, roughly as $t^2$. We perform simulations of star formation in self-gravitating hydrodynamic and magneto-hydrodynamic turbulence that is continuously driven. By turning the self-gravity of the gas in the simulations on or off, we demonstrate that self-gravity is the dominant physical effect setting the mass accretion rate at early times before feedback effects take over, contrary to theories of turbulence-regulated star formation. We find that gravitational collapse steepens the density profile around stars, generating the power-law tail on what is otherwise a lognormal density probability distribution function. Furthermore, we find turbulent velocity profiles to flatten inside collapsing regions, altering the size-linewidth relation. This local flattening reflects enhancements of turbulent velocity on small scales, as verified by changes to the velocity power spectra. Our results indicate that gas self-gravity dynamically alters both density and velocity structures in clouds, giving rise to a time-varying star formation rate. We find that a substantial fraction of the gas that forms stars arrives via low density flows, as opposed to accreting through high density filaments.

Recent first detections of the cross-correlation of the thermal Sunyaev-Zel'dovich (tSZ) signal in Planck cosmic microwave background (CMB) temperature maps with gravitational lensing maps inferred from the Planck CMB data and the CFHTLenS galaxy survey provide new probes of the relationship between baryons and dark matter. Using cosmological hydrodynamics simulations, we show that these cross-correlation signals are dominated by contributions from hot gas in the intracluster medium (ICM), rather than diffuse, unbound gas located beyond the virial radius (the "missing baryons"). Thus, these cross-correlations offer a tool with which to study the ICM over a wide range of halo masses and redshifts. In particular, we show that the tSZ - CMB lensing cross-correlation is more sensitive to gas in lower-mass, higher-redshift halos and gas at larger cluster-centric radii than the tSZ - galaxy lensing cross-correlation. Combining these measurements with primary CMB data will constrain feedback models through their signatures in the ICM pressure profile. We forecast the ability of ongoing and future experiments to constrain such ICM parameters, including the mean amplitude of the pressure - mass relation, the redshift evolution of this amplitude, and the mean outer logarithmic slope of the pressure profile. The results are promising, with $\approx 5-20$% precision constraints achievable with upcoming experiments, even after marginalizing over cosmological parameters.

In this work we extend our dynamical study of Ultra Compact X-ray Binaries (UCXB) 4U 1820-30 from Prodan and Murray 2012 to three more UCXBs in globular clusters: 4U 1850-087, 4U 0513-40 and M15 X-2. These three UCXBs have orbital periods < 20 mins. Two of them, 4U 1850-087 and 4U 0513-40, have suspected luminosity variations of order of ~ 1yr. There is insufficient observational data to make any statements regarding the long periodicity in the light curve of M15 X-2 at this point. The properties of these three systems are quite similar to 4U 1820-30, which prompt us to model their dynamics in the same manner. As in the case of 4U 1820-30, we interpret the suspected long periods as the period of small oscillations around a stable fixed point in the Kozai resonance. We provide a lower limit on the tidal dissipation factor Q which is in agreement with results obtained for the case of 4U 1820-30.

We present a series of high-resolution cosmological simulations of galaxy formation to z=0, spanning halo masses ~10^8-10^13 M_sun, and stellar masses ~10^4-10^11. Our simulations include fully explicit treatment of both the multi-phase ISM (molecular through hot) and stellar feedback. The stellar feedback inputs (energy, momentum, mass, and metal fluxes) are taken directly from stellar population models. These sources of stellar feedback, with zero adjusted parameters, reproduce the observed relation between stellar and halo mass up to M_halo~10^12 M_sun (including dwarfs, satellites, MW-mass disks, and small groups). By extension, this leads to reasonable agreement with the stellar mass function for M_star<10^11 M_sun. We predict weak redshift evolution in the M_star-M_halo relation, consistent with current constraints to z>6. We find that the M_star-M_halo relation is insensitive to numerical details, but is sensitive to the feedback physics. Simulations with only supernova feedback fail to reproduce the observed stellar masses, particularly in dwarf and high-redshift galaxies: radiative feedback (photo-heating and radiation pressure) is necessary to disrupt GMCs and enable efficient coupling of later supernovae to the gas. Star formation rates agree well with the observed Kennicutt relation at all redshifts. The galaxy-averaged Kennicutt relation is very different from the numerically imposed law for converting gas into stars in the simulation, and is instead determined by self-regulation via stellar feedback. Feedback reduces star formation rates considerably and produces a reservoir of gas that leads to rising late-time star formation histories significantly different from the halo accretion history. Feedback also produces large short-timescale variability in galactic SFRs, especially in dwarfs. Many of these properties are not captured by common 'sub-grid' galactic wind models.

We present a model of star formation in self-gravitating turbulent gas. We treat the turbulent velocity $v_T$ as a dynamical variable, and assume that it is adiabatically heated by the collapse. The theory predicts the run of density, infall velocity, and turbulent velocity, and the rate of star formation in compact massive gas clouds. The turbulent pressure is dynamically important at all radii, a result of the adiabatic heating. The system evolves toward a coherent spatial structure with a fixed run of density, $\rho(r,t)\to\rho(r)$; mass flows through this structure onto the central star or star cluster. We define the sphere of influence of the accreted matter by $m_*=M_g(r_*)$, where $m_*$ is the stellar plus disk mass in the nascent star cluster and $M_g(r)$ is the gas mass inside radius $r$. The density is given by a broken power law with a slope $-1.5$ inside $r_*$ and $\sim -1.6$ to $-1.8$ outside $r_*$. Both $v_T$ and the infall velocity $|u_r|$ decrease with decreasing $r$ for $r>r_*$; $v_T(r)\sim r^p$, the size-linewidth relation, with $p\approx0.2-0.3$, explaining the observation that Larson's Law is altered in massive star forming regions. The infall velocity is generally smaller than the turbulent velocity at $r>r_*$. For $r<r_*$, the infall and turbulent velocities are again similar, and both increase with decreasing $r$ as $r^{-1/2}$, with a magnitude about half of the free-fall velocity. The accreted (stellar) mass grows super-linearly with time, $\dot M_*=\phi M_{\rm cl}(t/\tau_{ff})^2$, with $\phi$ a dimensionless number somewhat less than unity, $M_{\rm cl}$ the clump mass and $\tau_{ff}$ the free-fall time of the clump. We suggest that small values of p can be used as a tracer of convergent collapsing flows.

We describe the long-term evolution of compact systems of terrestrial planets, using a set of simulations that match the statistical properties of the observed exoplanet distribution. The evolution is driven by tidal dissipation in the planetary interiors, but the systems evolve as a whole due to secular gravitational interactions. We find that, for Earth-like dissipation levels, planetary orbits can be circularised out to periods of order 100 days, an order of magnitude larger than is possible for single planets. The resulting distribution of eccentricities is a qualitative match to that inferred from transit timing variations, with a minority of non-zero eccentricities maintained by particular secular configurations. The coupling of the tidal and secular processes enhance the inward migration of the innermost planets in these systems, and can drive them to short orbital periods. Resonant interactions of both the mean motion and secular variety are observed, although the interactions are not strong enough to drive systemic instability in most cases. However, we demonstrate that these systems can easily be driven unstable if coupled to giant planets on longer period orbits.

We report ALMA Early Science observations of the Abell 1835 brightest cluster galaxy (BCG) in the CO (3-2) and CO (1-0) emission lines. We detect $5\times 10^{10}~\rm M_\odot$ of molecular gas within 10 kpc of the BCG. Its ensemble velocity profile width of $\sim 130 ~\rm km~s^{-1}$ FWHM is too narrow for the molecular cloud sto be supported in the galaxy by dynamic pressure. The gas may instead be supported in a rotating, turbulent disk oriented nearly face-on. Roughly $10^{10}~\rm M_\odot$ of molecular gas is projected $3-10 ~\rm kpc$ to the north-west and to the east of the nucleus with line of sight velocities lying between $-250 ~\rm km~s^{-1}$ to $+480 ~\rm km~s^{-1}$ with respect to the systemic velocity. The high velocity gas may be either inflowing or outflowing. However, the absence of high velocity gas toward the nucleus that would be expected in a steady inflow, and its bipolar distribution on either side of the nucleus, are more naturally explained as outflow. Star formation and radiation from the AGN are both incapable of driving an outflow of this magnitude. If so, the molecular outflow may be associated a hot outflow on larger scales reported by Kirkpatrick and colleagues. The molecular gas flow rate of approximately $200~\rm M_\odot ~yr^{-1}$ is comparable to the star formation rate of $100-180~\rm M_\odot ~yr^{-1}$ in the central disk. How radio bubbles would lift dense molecular gas in their updrafts, how much gas will be lost to the BCG, and how much will return to fuel future star formation and AGN activity are poorly understood. Our results imply that radio-mechanical (radio mode) feedback not only heats hot atmospheres surrounding elliptical galaxies and BCGs, it is able to sweep higher density molecular gas away from their centers.

The star formation rate (SFR) in the Central Molecular Zone (CMZ, i.e. the central 500 pc) of the Milky Way is lower by a factor of >10 than expected for the substantial amount of dense gas it contains, which challenges current star formation theories. In this paper, we quantify which physical mechanisms could be responsible. On scales larger than the disc scale height, the low SFR is found to be consistent with episodic star formation due to secular instabilities or possibly variations of the gas inflow along the Galactic bar. The CMZ is marginally Toomre-stable when including gas and stars, but highly Toomre-stable when only accounting for the gas, indicating a low condensation rate of self-gravitating clouds. On small scales, we find that the SFR in the CMZ may be caused by an elevated critical density for star formation due to the high turbulent pressure. The existence of a universal density threshold for star formation is ruled out. The HI-H$_2$ phase transition of hydrogen, the tidal field, a possible underproduction of massive stars due to a bottom-heavy initial mass function, magnetic fields, and cosmic ray or radiation pressure feedback also cannot individually explain the low SFR. We propose a self-consistent cycle of star formation in the CMZ, in which the effects of several different processes combine to inhibit star formation. The rate-limiting factor is the slow evolution of the gas towards collapse - once star formation is initiated it proceeds at a normal rate. The ubiquity of star formation inhibitors suggests that a lowered central SFR should be a common phenomenon in other galaxies. We discuss the implications for galactic-scale star formation and supermassive black hole growth, and relate our results to the star formation conditions in other extreme environments.

We use our variable Eddington tensor (VET) radiation hydrodynamics code to perform two-dimensional simulations to study the impact of radiation forces on atmospheres composed of dust and gas. Our setup closely follows that of Krumholz & Thompson, assuming that dust and gas are well-coupled and that the radiation field is characterized by blackbodies with temperatures >~ 80 K, as might be found in ultraluminous infrared galaxies. In agreement with previous work, we find that Rayleigh-Taylor instabilities develop in radiation supported atmospheres, leading to inhomogeneities that limit momentum exchange between radiation and dusty gas, and eventually providing a near balance of the radiation and gravitational forces. However, the evolution of the velocity and spatial distributions of the gas differs significantly from previous work, which utilized a less accurate flux-limited diffusion (FLD) method. Our VET simulations show continuous net acceleration of the gas, with no steady-state reached by the end of the simulation. In contrast, FLD results show little net acceleration of the gas and settle in to a quasi-steady, turbulent state with low velocity dispersion. The discrepancies result primarily from the inability of FLD to properly model the variation of the radiation field around structures that are less than a few optical depths across. We conclude that radiation feedback remains a viable mechanism for driving high-Mach number turbulence. We discuss implications for observed systems and global numerical simulations of feedback, but more realistic setups are needed to make robust observational predictions and assess the prospect of launching outflows with radiation.

Galactic superwinds may be driven by very hot outflows generated by overlapping supernovae within the host galaxy. We use the Chevalier & Clegg (CC85) wind model and the observed correlation between X-ray luminosities of galaxies and their SFRs to constrain the mass loss rates (\dot{M}_hot) across a wide range of star formation rates (SFRs), from dwarf starbursts to ultra-luminous infrared galaxies. We show that for fixed thermalization efficiency and mass loading rate, the X-ray luminosity of the hot wind scales as L_X ~ SFR^2, significantly steeper than is observed for star-forming galaxies: L_X ~ SFR. Using this difference we constrain the mass-loading and thermalization efficiency of hot galactic winds. For reasonable values of the thermalization efficiency (<~ 1) and for SFR >~ 10 M_sun/yr we find that \dot{M}_hot/SFR <~ 1, significantly lower than required by integrated constraints on the efficiency of stellar feedback in galaxies, and potentially too low to explain observations of winds from rapidly star-forming galaxies. In addition, we highlight the fact that heavily mass-loaded winds cannot be described by the adiabatic CC85 model because they become strongly radiative.

We report the first detection of a correlation between gravitational lensing by large scale structure and the thermal Sunyaev-Zeldovich (tSZ) effect. Using the mass map from the Canada France Hawaii Telescope Lensing Survey (CFHTLenS) and a newly-constructed tSZ map from Planck, we measure a non-zero correlation between the two maps out to one degree angular separation on the sky, with an overall significance of 6 sigma. The tSZ maps are formed in a manner that removes primary cosmic microwave background fluctuations and minimizes residual contamination by galactic and extragalactic dust emission, and by CO line emission. We perform numerous tests to show that our measurement is immune to these residual contaminants. The resulting correlation function is consistent with the existence of a warm baryonic gas tracing the large scale structure with a bias b_gas. Given the shape of the lensing kernel, our signal sensitivity peaks at a redshift z~0.4, where half a degree separation on the sky corresponds to a physical scale of ~10 Mpc. The amplitude of the signal constrains the product (b_gas/1)(T_e / 0.1 keV)(n_e / 1 m^-3)=2.01\pm 0.31\pm 0.21, at redshift zero. Our study suggests that a substantial fraction of the "missing" baryons in the universe may reside in a low density warm plasma that traces dark matter.

We survey the present landscape in submillimetre astronomy for Canada and describe a plan for continued engagement in observational facilities to ~2020. Building on Canada's decadal Long Range Plan process, we emphasize that continued involvement in a large, single-dish facility is crucial given Canada's substantial investment in ALMA and numerous PI-led submillimetre experiments. In particular, we recommend: i) an extension of Canadian participation in the JCMT until at least the unique JCMT Legacy Survey program is able to realize the full scientific potential provided by the world-leading SCUBA-2 instrument; and ii) involvement of the entire Canadian community in CCAT, with a large enough share in the partnership for Canadian astronomers to participate at all levels of the facility. We further recommend continued participation in ALMA development, involvement in many focused PI-led submillimetre experiments, and partnership in SPICA.

White dwarf-white dwarf (WD-WD) mergers may lead to type Ia supernovae events. Thompson (2011) suggested that many such binaries are produced in hierarchical triple systems. The tertiary induces eccentricity oscillations in the inner binary via the Kozai-Lidov mechanism, driving the binary to high eccentricities, and significantly reducing the gravitational wave merger timescale (T_GW) over a broad range of parameter space. Here, we investigate the role of tidal forces in these systems. We show that tidal effects are important in the regime of moderately high initial relative inclination between the inner binary and the outer tertiary. For 85 < i_0 < 90 degrees (prograde) and 97 < i_0 < 102 degrees (retrograde), tides combine with GW radiation to dramatically decrease T_GW. In the regime of high inclinations between 91 < i_0 < 96 degrees, the inner binary likely suffers a direct collision, as in the work of Katz & Dong (2012) and tidal effects do not play an important role.

Coalescing black hole (BH) binaries forming in the dense core of globular clusters (GCs) are expected to be one the brightest sources of gravitational wave (GW) radiation for the next generation of ground-based laser interferometers. Favorable conditions for merger are initiated by the Kozai resonance in which the gravitational interaction with a third distant object, typically another BH, induces quasi-periodic variations of the inner BH binary eccentricity. In this paper we perform high precision N-body simulations of the long term evolution of hierarchical BH triples and investigate the conditions that lead to the merging of the BH binary and the way it might become an observable source of GW radiation. We find that the secular orbit average treatment, adopted in previous works, does not reliably describe the dynamics of these systems if the binary is orbited by the outer BH on a highly inclined orbit at a moderate distance. During the high eccentricity phase of a Kozai cycle the torque due to the outer BH can drive the binary to extremely large eccentricities in a fraction of the binary's orbital period. This occurs before relativistic terms become important to the evolution and allows the binary GW signal to reach large GW frequencies (>~10 Hz) at high eccentricities. We show that ~50 % of coalescing BH binaries driven by the Kozai mechanism in GCs will have eccentricities larger than 0.1, with 10 % of them being extremely eccentric, (1-e)<~10^-5, when they first chirp in the frequency band of ground based laser interferometers. This implies that a large fraction of such GW sources could be missed if conventional quasi-circular templates are used for analysis of GW detectors data. The efficient detection of all coalescing BH binaries in GCs will therefore require template banks of eccentric inspiral waveforms for matched-filtering and dedicated search strategies.

We report ALMA Early Science observations of the Abell 1835 brightest cluster galaxy (BCG) in the CO (3-2) and CO (1-0) emission lines. We detect 5E10 solar masses of molecular gas within 10 kpc of the BCG. Its velocity width of ~130 km/s FWHM is too narrow to be supported by dynamical pressure. The gas may instead be supported in a rotating, turbulent disk oriented nearly face-on. The disk is forming stars at a rate of 100-180 solar masses per year. Roughly 1E10 solar masses of molecular gas is projected 3-10 kpc to the north-west and to the east of the nucleus with line of sight velocities lying between -250 km/s to +480 km/s with respect to the systemic velocity. Although inflow cannot be ruled out, the rising velocity gradient with radius is consistent with a broad, bipolar outflow driven by radio jets or buoyantly rising X-ray cavities. The molecular outflow may be associated with an outflow of hot gas in Abell 1835 seen on larger scales. Molecular gas is flowing out of the BCG at a rate of approximately 200 solar masses per year, which is comparable to its star formation rate. How radio bubbles lift dense molecular gas in their updrafts, how much gas will be lost to the BCG, and how much will return to fuel future star formation and AGN activity are poorly understood. Our results imply that radio-mechanical (radio mode) feedback not only heats hot atmospheres surrounding elliptical galaxies and BCGs, it is able to sweep higher density molecular gas away from their centers.

We report ALMA Early Science CO(1-0) and CO(3-2) observations of the brightest cluster galaxy (BCG) in Abell 1664. The BCG contains 1.1x10^{10} solar masses of molecular gas divided roughly equally between two distinct velocity systems: one from -250 to +250 km/s centred on the BCG's systemic velocity and a high velocity system blueshifted by 570 km/s with respect to the systemic velocity. The BCG's systemic component shows a smooth velocity gradient across the BCG center with velocity proportional to radius suggestive of solid body rotation about the nucleus. However, the mass and velocity structure are highly asymmetric and there is little star formation coincident with a putative disk. It may be an inflow of gas that will settle into a disk over several 10^8 yr. The high velocity system consists of two gas clumps, each ~2 kpc across, located to the north and southeast of the nucleus. Each has a line of sight velocity spread of 250-300 km/s. The velocity of the gas in the high velocity system tends to increase towards the BCG center and could signify a massive high velocity flow onto the nucleus. However, the velocity gradient is not smooth and these structures are also coincident with low optical-UV surface brightness regions, which could indicate dust extinction associated with each clump. If so, the high velocity gas would be projected in front of the BCG and moving toward us along the line of sight in a massive outflow most likely driven by the AGN. A merger origin is unlikely but cannot be ruled out.

We use simulations with realistic models for stellar feedback to study galaxy mergers. These high resolution (1 pc) simulations follow formation and destruction of individual GMCs and star clusters. The final starburst is dominated by in situ star formation, fueled by gas which flows inwards due to global torques. The resulting high gas density results in rapid star formation. The gas is self gravitating, and forms massive (~10^10 M_sun) GMCs and subsequent super-starclusters (masses up to 10^8 M_sun). However, in contrast to some recent simulations, the bulk of new stars which eventually form the central bulge are not born in superclusters which then sink to the center of the galaxy, because feedback efficiently disperses GMCs after they turn several percent of their mass into stars. Most of the mass that reaches the nucleus does so in the form of gas. The Kennicutt-Schmidt law emerges naturally as a consequence of feedback balancing gravitational collapse, independent of the small-scale star formation microphysics. The same mechanisms that drive this relation in isolated galaxies, in particular radiation pressure from IR photons, extend over seven decades in SFR to regulate star formation in the most extreme starbursts (densities >10^4 M_sun/pc^2). Feedback also drives super-winds with large mass loss rates; but a significant fraction of the wind material falls back onto the disks at later times, leading to higher post-starburst SFRs in the presence of stellar feedback. Strong AGN feedback is required to explain sharp cutoffs in star formation rate. We compare the predicted relic structure, mass profile, morphology, and efficiency of disk survival to simulations which do not explicitly resolve GMCs or feedback. Global galaxy properties are similar, but sub-galactic properties and star formation rates can differ significantly.

We consider the effects of different star formation criteria on galactic scales, in high-resolution simulations with explicitly resolved GMCs and stellar feedback. We compare: (1) a self-gravity criterion (based on the local virial parameter and the assumption that self-gravitating gas collapses to high density in a free-fall time), (2) a fixed density threshold, (3) a molecular-gas law, (4) a temperature threshold, (5) a Jeans-instability requirement, (6) a criteria that cooling times be shorter than dynamical times, and (7) a convergent-flow criterion. We consider these both MW-like and high-density (starburst) galaxies. With feedback present, all models produce identical integrated star formation rates (SFRs), in agreement with the Kennicutt relation. Without feedback all produce orders-of-magnitude excessive SFRs. This is totally dependent on feedback and independent of the SF law. However, the spatial and density distribution of SF depend strongly on the SF criteria. Because cooling rates are generally fast and gas is turbulent, criteria (4)-(7) are 'weak' and spread SF uniformly over the disk (above densities n~0.01-0.1 cm^-3). A molecular criterion (3) localizes to higher densities, but still a wide range; for Z Z_solar, it is similar to a density threshold at n~1 cm^-3 (well below mean densities in the MW center or starbursts). Fixed density thresholds (2) can always select the highest densities, but must be adjusted for simulation resolution and galaxy properties; the same threshold that works in a MW-like simulation will select nearly all gas in a starburst. Binding criteria (1) tend to adaptively select the largest over-densities, independent of galaxy model or resolution, and automatically predict clustered SF. We argue that this SF model is most physically-motivated and presents significant numerical advantages in large-dynamic range simulations.

Rapid accretion of cold gas plays a crucial role in getting gas into galaxies. It has been suggested that this accretion proceeds along narrow streams that might also directly drive the turbulence in galactic gas, dynamical disturbances, and bulge formation. In cosmological simulations, however, it is impossible to isolate and hence disentangle the effect of accretion from internal instabilities and mergers. Moreover, in most cosmological simulations, the phase structure and turbulence in the ISM arising from stellar feedback are treated in a sub-grid manner, so that feedback cannot generate ISM turbulence. In this paper we therefore test the effects of cold streams in extremely high-resolution simulations of otherwise isolated galaxy disks using detailed models for star formation and feedback; we then include or exclude mock cold flows falling onto the galaxies with accretion rates, velocities and geometry set to maximize their effect on the disk. We find: (1) Turbulent velocity dispersions in gas disks are identical with or without the cold flow; the energy injected by the flow is dissipated where it meets the disk. (2) In runs without stellar feedback, the presence of a cold flow has essentially no effect on runaway local collapse, resulting in star formation rates (SFRs) that are far too large. (3) Disks in runs with feedback and cold flows have higher SFRs, but only insofar as they have more gas. (4) Because flows are extended relative to the disk, they do not trigger strong resonant responses and so induce weak morphological perturbation (bulge formation via instabilities is not accelerated). (5) However, flows can thicken the disk by direct contribution of out-of-plane streams. We conclude that while inflows are critical over cosmological timescales to determine the supply and angular momentum of gas disks, they have weak instantaneous dynamical effects on galaxies.

We study galaxy super-winds driven in major mergers, using pc-resolution simulations with detailed models for stellar feedback that can self-consistently follow the formation/destruction of GMCs and generation of winds. The models include molecular cooling, star formation at high densities in GMCs, and gas recycling and feedback from SNe (I&II), stellar winds, and radiation pressure. We study mergers of systems from SMC-like dwarfs and Milky Way analogues to z~2 starburst disks. Multi-phase super-winds are generated in all passages, with outflow rates up to ~1000 M_sun/yr. However, the wind mass-loading efficiency (outflow rate divided by SFR) is similar to that in isolated galaxy counterparts of each merger: it depends more on global galaxy properties (mass, size, escape velocity) than on the dynamical state of the merger. Winds tend to be bi- or uni-polar, but multiple 'events' build up complex morphologies with overlapping, differently-oriented bubbles/shells at a range of radii. The winds have complex velocity and phase structure, with material at a range of speeds up to ~1000 km/s, and a mix of molecular, ionized, and hot gas that depends on galaxy properties and different feedback mechanisms. These simulations resolve a problem in some 'sub-grid' models, where simple wind prescriptions can dramatically suppress merger-induced starbursts. But despite large mass-loading factors (>~10) in the winds, the peak SFRs are comparable to those in 'no wind' simulations. Wind acceleration does not act equally, so cold dense gas can still lose angular momentum and form stars, while blowing out gas that would not have participated in the starburst in the first place. Considerable wind material is not unbound, and falls back on the disk at later times post-merger, leading to higher post-starburst SFRs in the presence of stellar feedback. This may require AGN feedback to explain galaxy quenching.

We show that the mass fraction of GMC gas (n>100 cm^-3) in dense (n>>10^4 cm^-3) star-forming clumps, observable in dense molecular tracers (L_HCN/L_CO(1-0)), is a sensitive probe of the strength and mechanism(s) of stellar feedback. Using high-resolution galaxy-scale simulations with pc-scale resolution and explicit models for feedback from radiation pressure, photoionization heating, stellar winds, and supernovae (SNe), we make predictions for the dense molecular gas tracers as a function of GMC and galaxy properties and the efficiency of stellar feedback. In models with weak/no feedback, much of the mass in GMCs collapses into dense sub-units, predicting L_HCN/L_CO(1-0) ratios order-of-magnitude larger than observed. By contrast, models with feedback properties taken directly from stellar evolution calculations predict dense gas tracers in good agreement with observations. Changing the strength or timing of SNe tends to move systems along, rather than off, the L_HCN-L_CO relation (because SNe heat lower-density material, not the high-density gas). Changing the strength of radiation pressure (which acts efficiently in the highest density gas), however, has a much stronger effect on L_HCN than on L_CO. We predict that the fraction of dense gas (L_HCN/L_CO(1-0)) increases with increasing GMC surface density; this drives a trend in L_HCN/L_CO(1-0) with SFR and luminosity which has tentatively been observed. Our results make specific predictions for enhancements in the dense gas tracers in unusually dense environments such as ULIRGs and galactic nuclei (including the galactic center).

Mapping Mg II resonance emission scattered by galactic winds offers a means to determine the spatial extent and density of the warm outflow. Using Keck/LRIS spectroscopy, we have resolved scattered Mg II emission to the east of 32016857, a star-forming galaxy at z =0.9392 with an outflow. The Mg II emission from this galaxy exhibits a P-Cygni profile, extends further than both the continuum and [O II] emission along the eastern side of the slit, and has a constant Doppler shift along the slit which does not follow the velocity gradient of the nebular [O II] emission. Using the Sobolev approximation, we derive the density of Mg+ ions at a radius of 12 to 18 kpc in the outflow. We model the ionization correction and find that much of the outflowing Mg is in Mg++. We estimate that the total mass flux could be as large as 330 - 500 solar masses per year, with the largest uncertainties coming from the depletion of Mg onto grains and the clumpiness of the warm outflow. We show that confining the warm clouds with a hot wind reduces the estimated mass flux of the warm outflow and indicates amass-loading factor near unity in the warm phase alone. Based on the high blue luminosities that distinguish 32016857 and TKRS 4389, described by Rubin et al. 2011, from other galaxies with P-Cygni emission, we suggest that, as sensitivity to diffuse emission improves, scattering halos may prove to be a generic property of star-forming galaxies at intermediate redshifts.

The European Space Agency (ESA) Gaia satellite has 106 CCD image sensors which will suffer from increased charge transfer inefficiency (CTI) as a result of radiation damage. To aid the mitigation at low signal levels, the CCD design includes Supplementary Buried Channels (SBCs, otherwise known as `notches') within each CCD column. We present the largest published sample of Gaia CCD SBC Full Well Capacity (FWC) laboratory measurements and simulations based on 13 devices. We find that Gaia CCDs manufactured post-2004 have SBCs with FWCs in the upper half of each CCD that are systematically smaller by two orders of magnitude (<50 electrons) compared to those manufactured pre-2004 (thousands of electrons). Gaia's faint star (13 < G < 20 mag) astrometric performance predictions by Prod'homme et al. and Holl et al. use pre-2004 SBC FWCs as inputs to their simulations. However, all the CCDs already integrated onto the satellite for the 2013 launch are post-2004. SBC FWC measurements are not available for one of our five post-2004 CCDs but the fact it meets Gaia's image location requirements suggests it has SBC FWCs similar to pre-2004. It is too late to measure the SBC FWCs onboard the satellite and it is not possible to theoretically predict them. Gaia's faint star astrometric performance predictions depend on knowledge of the onboard SBC FWCs but as these are currently unavailable, it is not known how representative of the whole focal plane the current predictions are. Therefore, we suggest Gaia's initial in-orbit calibrations should include measurement of the onboard SBC FWCs. We present a potential method to do this. Faint star astrometric performance predictions based on onboard SBC FWCs at the start of the mission would allow satellite operating conditions or CTI software mitigation to be further optimised to improve the scientific return of Gaia.

We present a Monte Carlo model for the structure of low mass (total mass < 25 earth mass) planetary systems that form by the in situ gravitational assembly of planetary embryos into final planets. Our model includes distributions of mass, eccentricity, inclination and period spacing that are based on the simulation of a disk of 20 earth masses, forming planets around a solar mass star, and assuming a power law surface density distribution $\propto a^{-1.5}$. The output of the Monte Carlo model is then subjected to the selection effects that mimic the observations of a transiting planet search such as that performed by the Kepler satellite. The resulting comparison of the output to the properties of the observed sample yields an encouraging agreement in terms of the relative frequencies of multiple planet systems and the distribution of the mutual inclinations, when moderate tidal circularisation is taken into account. The broad features of the period distribution and radius distribution can also be matched within this framework, although the model underpredicts the distribution of small period ratios. This likely indicates that some dissipation is still required in the formation process. The most striking deviation between model and observations is in the ratio of single to multiple systems, in that there are roughly 50% more single planet candidates observed than are produced in any model population. This suggests that some systems must suffer additional attrition to reduce the number of planets or increase the range of inclinations.

We use numerical simulations of isolated galaxies to study the effects of stellar feedback on the formation and evolution of giant star-forming gas 'clumps' in high-redshift, gas-rich galaxies. Such galactic disks are unstable to the formation of bound gas-rich clumps whose properties initially depend only on global disk properties, not the microphysics of feedback. In simulations without stellar feedback, clumps turn an order-unity fraction of their mass into stars and sink to the center, forming a large bulge and kicking most of the stars out into a much more extended stellar envelope. By contrast, strong radiative stellar feedback disrupts even the most massive clumps after they turn ~10-20% of their mass into stars, in a timescale of ~10-100 Myr, ejecting some material into a super-wind and recycling the rest of the gas into the diffuse ISM. This suppresses the bulge formation rate by direct 'clump coalescence' by a factor of several. However, the galactic disks do undergo significant internal evolution in the absence of mergers: clumps form and disrupt continuously and torque gas to the galactic center. The resulting evolution is qualitatively similar to bar/spiral evolution in simulations with a more homogeneous ISM.

We investigate how strong lensing of dusty, star-forming galaxies by foreground galaxies can be used as a probe of dark matter halo substructure. We find that spatially resolved spectroscopy of lensed sources allows dramatic improvements to measurements of lens parameters. In particular we find that modeling of the full, three-dimensional (angular position and radial velocity) data can significantly facilitate substructure detection, increasing the sensitivity of observables to lower mass subhalos. We carry out simulations of lensed dusty sources observed by early ALMA (Cycle 1) and use a Fisher matrix analysis to study the parameter degeneracies and mass detection limits of this method. We find that, even with conservative assumptions, it is possible to detect galactic dark matter subhalos of ~ 10^8 M_{\odot} with high significance in most lensed DSFGs. Specifically, we find that in typical DSFG lenses, there is a ~ 55 % probability of detecting a substructure with M>10^8 M_{\odot} with more than 5 sigma detection significance in each lens, if the abundance of substructure is consistent with previous lensing results. The full ALMA array, with its significantly enhanced sensitivity and resolution, should improve these estimates considerably. Given the sample of ~100 lenses provided by surveys like the South Pole Telescope, our understanding of dark matter substructure in typical galaxy halos is poised to improve dramatically over the next few years.

The conversion of gas into stars is a fundamental process in astrophysics and cosmology. Stars are known to form from the gravitational collapse of dense clumps in interstellar molecular clouds, and it has been proposed that the resulting star formation rate is proportional to either the amount of mass above a threshold gas surface density, or the gas volume density. These star-formation prescriptions appear to hold in nearby molecular clouds in our Milky Way Galaxy's disk as well as in distant galaxies where the star formation rates are often much larger. The inner 500 pc of our Galaxy, the Central Molecular Zone (CMZ), contains the largest concentration of dense, high-surface density molecular gas in the Milky Way, providing an environment where the validity of star-formation prescriptions can be tested. Here we show that by several measures, the current star formation rate in the CMZ is an order-of-magnitude lower than the rates predicted by the currently accepted prescriptions. In particular, the region 1 deg < l < 3.5 deg, |b| < 0.5 deg contains ~10^7 Msun of dense molecular gas -- enough to form 1000 Orion-like clusters -- but the present-day star formation rate within this gas is only equivalent to that in Orion. In addition to density, another property of molecular clouds, such as the amplitude of turbulent motions, must be included in the star-formation prescription to predict the star formation rate in a given mass of molecular gas.

Euclid-VIS is a large format visible imager for the ESA Euclid space mission in their Cosmic Vision program, scheduled for launch in 2019. Together with the near infrared imaging within the NISP instrument it forms the basis of the weak lensing measurements of Euclid. VIS will image in a single r+i+z band from 550-900 nm over a field of view of ~0.5 deg2. By combining 4 exposures with a total of 2240 sec, VIS will reach to V=24.5 (10{\sigma}) for sources with extent ~0.3 arcsec. The image sampling is 0.1 arcsec. VIS will provide deep imaging with a tightly controlled and stable point spread function (PSF) over a wide survey area of 15000 deg2 to measure the cosmic shear from nearly 1.5 billion galaxies to high levels of accuracy, from which the cosmological parameters will be measured. In addition, VIS will also provide a legacy imaging dataset with an unprecedented combination of spatial resolution, depth and area covering most of the extra-Galactic sky. Here we will present the results of the study carried out by the Euclid Consortium during the Euclid Definition phase.

Feedback from massive stars is believed to play a critical role in driving galactic super-winds that enrich the IGM and shape the galaxy mass function and mass-metallicity relation. In previous papers, we introduced new numerical methods for implementing stellar feedback on sub-GMC through galactic scales in galaxy simulations. This includes radiation pressure (UV through IR), SNe (Type-I & II), stellar winds ('fast' O-star through 'slow' AGB winds), and HII photoionization. Here, we show that these feedback mechanisms drive galactic winds with outflow rates as high as ~10-20 times the galaxy SFR. The mass-loading efficiency (wind mass loss rate divided by SFR) scales inversely with circular velocity, consistent with momentum-conservation expectations. We study the contributions of each feedback mechanism to galactic winds in a range of galaxy models, from SMC-like dwarfs & MW-analogues to z~2 clumpy disks. In massive, gas-rich systems (local starbursts and high-z galaxies), radiation pressure dominates the wind generation. For MW-like spirals and dwarf galaxies the gas densities are much lower, and shock-heated gas from SNe and stellar winds dominates production of large-scale outflows. In all models, however, winds have a multi-phase structure that depends on interactions between multiple feedback mechanisms operating on different spatial & time scales: any single mechanism fails to reproduce the winds observed. We provide fitting functions for wind mass-loading and velocities as a function of galaxy properties, for use in cosmological simulations and semi-analytic models. These differ from typically-adopted formulae with explicit dependence on gas surface density that can be very important in both low-density dwarf galaxies and high-density gas-rich galaxies.

We present numerical methods for including stellar feedback in galaxy-scale simulations. We include heating by SNe (I & II), gas recycling and shock-heating from O-star & AGB winds, HII photoionization, and radiation pressure from stellar photons. The energetics and time-dependence are taken directly from stellar evolution models. We implement these in simulations with pc-scale resolution, modeling galaxies from SMC-like dwarfs and MW analogues to massive z~2 starburst disks. Absent feedback, gas cools and collapses without limit. With feedback, the ISM reaches a multi-phase steady state in which GMCs continuously form, disperse, and re-form. Our primary results include: (1) Star forming galaxies generically self-regulate at Toomre Q~1. Most of the volume is in diffuse hot gas with most of the mass in dense GMC complexes. The phase structure and gas mass at high densities are much more sensitive probes of stellar feedback physics than integrated quantities (Toomre Q or gas velocity dispersion). (2) Different feedback mechanisms act on different scales: radiation & HII pressure are critical to prevent runaway collapse of dense gas in GMCs. SNe and stellar winds dominate the dynamics of volume-filling hot gas; however this primarily vents out of the disk. (3) The galaxy-averaged SFR is determined by feedback. For given feedback efficiency, restricting star formation to molecular gas or modifying the cooling function has little effect; but changing feedback mechanisms directly translates to shifts off the Kennicutt-Schmidt relation. (4) Self-gravity leads to marginally-bound GMCs with an ~M^-2 mass function with a cutoff at the Jeans mass; they live a few dynamical times before being disrupted by stellar feedback and turn ~1-10% of their mass into stars (increasing from dwarfs through starburst galaxies). Low-mass GMCs are preferentially unbound.

IRAS 09104+4109 is a rare example of a dust enshrouded type 2 QSO in the centre of a cool-core galaxy cluster. Previous observations of this z=0.44 system showed that as well as powering the hyper-luminous infrared emission of the cluster-central galaxy, the QSO is associated with a double-lobed radio source. However, the steep radio spectral index and misalignment between the jets and ionised optical emission suggested that the orientation of the QSO had recently changed. We use a combination of new, multi-band Giant Metrewave Radio Telescope observations and archival radio data to confirm that the jets are no longer powered by the QSO, and estimate their age to be 120-160 Myr. This is in agreement with the ~70-200 Myr age previously estimated for star-formation in the galaxy. Previously unpublished Very Long Baseline Array data reveal a 200 pc scale double radio source in the galaxy core which is more closely aligned with the current QSO axis and may represent a more recent period of jet activity. These results suggest that the realignment of the QSO, the cessation of jet activity, and the onset of rapid star-formation may have been caused by a gas-rich galaxy merger. A Chandra X-ray observation confirms the presence of cavities associated with the radio jets, and we estimate the energy required to inflate them to be ~7.7x10^60 erg. The mechanical power of the jets is sufficient to balance radiative cooling in the cluster, provided they are efficiently coupled to the intra-cluster medium (ICM). We find no evidence of direct radiative heating and conclude that the QSO either lacks the radiative luminosity to heat the ICM, or that it requires longer than 100-200 Myr to significantly impact its environment. [Abridged]

We analyze Spitzer GLIMPSE, MSX, and WMAP images of the Milky Way to identify 8 micron and free-free sources in the Galaxy. Seventy-two of the eighty-eight WMAP sources have coverage in the GLIMPSE and MSX surveys suitable for identifying massive star forming complexes (SFC). We measure the ionizing luminosity functions of the SFCs and study their role in the turbulent motion of the Galaxy's molecular gas. We find a total Galactic free-free flux f_{\nu} = 46177.6 Jy; the 72 WMAP sources with full 8 micron coverage account for 34263.5 Jy (~75%), with both measurements made at \nu=94GHz (W band). We find a total of 280 SFCs, of which 168 have unique kinematic distances and free-free luminosities. We use a simple model for the radial distribution of star formation to estimate the free-free and ionizing luminosity for the sources lacking distance determinations. The total dust-corrected ionizing luminosity is Q = 2.9 \pm 0.5 x 10^53 photons s^-1, which implies a galactic star formation rate of 1.2 \pm 0.2 M_{\sun} yr^-1. We present the (ionizing) luminosity function of the SFCs, and show that 24 sources emit half the ionizing luminosity of the Galaxy. The SFCs appear as bubbles in GLIMPSE or MSX images; the radial velocities associated with the bubble walls allow us to infer the expansion velocity of the bubbles. We calculate the kinetic luminosity of the bubble expansion and compare it to the turbulent luminosity of the inner molecular disk. SFCs emitting 80% of the total galactic free-free luminosity produce a kinetic luminosity equal to 65% of the turbulent luminosity in the inner molecular disk. This suggests that the expansion of the bubbles is a major driver of the turbulent motion of the inner Milky Way molecular gas.

We demonstrate that the observed distribution of `Hot Neptune'/`Super-Earth' systems is well reproduced by a model in which planet assembly occurs in situ, with no significant migration post-assembly. This is achieved only if the amount of mass in rocky material is $\sim 50$--$100 M_{\oplus}$ interior to 1 AU. Such a reservoir of material implies that significant radial migration of solid material takes place, and that it occur before the stage of final planet assembly. The model not only reproduces the general distribution of mass versus period, but also the detailed statistics of multiple planet systems in the sample. We furthermore demonstrate that cores of this size are also likely to meet the criterion to gravitationally capture gas from the nebula, although accretion is rapidly limited by the opening of gaps in the gas disk. If the mass growth is limited by this tidal truncation, then the scenario sketched here naturally produces Neptune-mass objects with substantial components of both rock and gas, as is observed. The quantitative expectations of this scenario are that most planets in the `Hot Neptune/Super-Earth' class inhabit multiple-planet systems, with characteristic orbital spacings. The model also provides a natural division into gas-rich (Hot Neptune) and gas-poor (Super-Earth) classes at fixed period. The dividing mass ranges from $\sim 3 M_{\oplus}$ at 10 day orbital periods to $\sim 10 M_{\oplus}$ at 100 day orbital periods. For orbital periods $< 10$ days, the division is less clear because a gas atmosphere may be significantly eroded by stellar radiation.

Photoionization modeling of the low-ionization broad absorption lines of certain quasars, known as FeLoBALs, has recently revealed the number density of the wind absorbers and their distance from the central supermassive black hole. From these, the feedback efficiency of the quasars can in principle be derived. The implied properties of the FeLoBALs are, however, surprising, with the thickness of the absorbers relative to their distance from the black hole, Delta R/R, as small as ~10^-5. Such absorbers are unlikely to survive the journey from the supermassive black hole to their inferred location. We show that the observed FeLoBAL properties are readily explained if they are formed in situ in radiative shocks produced when a quasar blast wave impacts a moderately dense interstellar clump along the line of sight. This physical picture differs significantly from the thin shell approximation often assumed, and implies outflow rates, kinetic luminosities and momentum fluxes that differ correspondingly, in some cases at the order of magnitude level. Using the radiative shock model, we estimate the ratio of the outflow kinetic luminosity to bolometric luminosity for three bright FeLoBAL quasars in the literature. We find Edot/Lbol~2-5% (and corresponding momentum fluxes Pdot~2-15 Lbol/c), similar to what is adopted in models reproducing the M-sigma relation. These outflow properties are also comparable to those recently inferred for molecular outflows in local ultra-luminous infrared galaxies, suggesting that active galactic nuclei are capable of driving such outflows.

It has been suggested that the 170 day period in the light curve of the low mass X-ray binary 4U 1820-30 arises from the presence of a third body with a large inclination to the binary orbit. We show that this long period motion arises if the system is librating around the stable fixed point in a Kozai resonance. We demonstrate that mass transfer drives the system toward this fixed point, and calculate, both analytically and via numerical integrations, that the period of libration is of order 170 days when the mutual inclination is near the Kozai critical value. The non-zero eccentricity of the binary, combined with tidal dissipation, implies that the rate of change of the binary period would be slower than, or even of opposite sign to, that implied by standard mass transfer models. If the 170 day period results from libration, then, contrary to appearances, the orbital period of the inner binary is increasing with time; in that case, (e/0.009)^2Q/k_2 > 2.5 x 10^9, where k_2 = 0.01 is the tidal Love number and e = 0.009 is the fiducial eccentricity of the inner binary. It appears unlikely that the observed negative period derivative results from the smaller than expected (but positive) value of \dot P combined with the previously suggested acceleration of the system in the gravitational field of the host globular cluster NGC 6624. The discrepancy between the observed and expected period derivative requires further investigation.

Feedback from massive stars is believed to play a critical role in shaping the galaxy mass function, the structure of the interstellar medium (ISM), and the low efficiency of star formation, but the exact form of the feedback is uncertain. In this paper, the first in a series, we present and test a novel numerical implementation of stellar feedback resulting from momentum imparted to the ISM by radiation, supernovae, and stellar winds. We employ a realistic cooling function, and find that a large fraction of the gas cools to <100K, so that the ISM becomes highly inhomogeneous. Despite this, our simulated galaxies reach an approximate steady state, in which gas gravitationally collapses to form giant molecular clouds (GMCs), dense clumps, and stars; subsequently, stellar feedback disperses the GMCs, repopulating the diffuse ISM. This collapse and dispersal cycle is seen in models of SMC-like dwarfs, the Milky-Way, and z~2 clumpy disk analogues. The simulated global star formation efficiencies are consistent with the observed Kennicutt-Schmidt relation. Moreover, the star formation rates are nearly independent of the numerically imposed high-density star formation efficiency, density threshold, and density scaling. This is a consequence of the fact that, in our simulations, star formation is regulated by stellar feedback limiting the amount of very dense gas available for forming stars. In contrast, in simulations without stellar feedback, i.e. under the action of only gravity and gravitationally-induced turbulence, the ISM experiences runaway collapse to very high densities. In these simulations without feedback, the global star formation rates exceed observed galactic star formation rates by 1-2 orders of magnitude, demonstrating that stellar feedback is crucial to the regulation of star formation in galaxies.

Essentially all stars form in giant molecular clouds (GMCs). However, inside GMCs, most of the gas does not participate in star formation; rather, denser gas accumulates in clumps in the GMC, with the bulk of the stars in a given GMC forming in a few of the most massive clumps. In the Milky Way, these clumps have masses $M_{\rm cl}\lesssim 5\times 10^{-2}$ of the GMC, radii $r_{\rm cl} \sim 1$pc, and free-fall times $\tau_{\rm cl} \sim 2\times 10^5\yr$. We show that clumps inside giant molecular clouds should accrete at a modified Bondi accretion rate, which depends on clump mass as $\dot M_{\rm cl}\sim M_{\rm cl}^{5/4}$. This rate is initially rather slow, usually slower than the initial star formation rate inside the clump (we adopt the common assumption that inside the clump, $\dot M_*=\epsilon_{\rm ff} M_{\rm cl}/\tau_{\rm cl}$, with $\epsilon_{\rm ff} \approx 0.017$). However, after $\sim 2$ GMC free-fall times $\tau_{\rm GMC}$, the clump accretion rate accelerates rapidly; formally, the clump can accrete the entire GMC in $\sim 3\tau_{\rm GMC}$. At the same time, the star formation rate accelerates, tracking the Bondi accretion rate. If the GMC is disrupted by feedback from the largest clump, half the stars in that clump form in the final $\taug$ before the GMC is disrupted. The theory predicts that the distribution of effective star formation rates, measured per GMC free-fall time, is broad, ranging from $\sim 0.001$ up to 0.1 or larger and that the mass spectrum of star clusters is flatter than that of clumps, consistent with observations.

We used WIRCam on CFHT to observe four transits of the super-Earth GJ 1214b in the near-infrared. For each transit we observed in two bands nearly-simultaneously by rapidly switching the WIRCam filter wheel back and forth for the duration of the observations. By combining all our J-band (~1.25 microns) observations we find a transit depth in this band of 1.338\pm0.013% - a value consistent with the optical transit depth reported by Charbonneau and collaborators. However, our best-fit combined Ks-band (~2.15 microns) transit depth is deeper: 1.438\pm0.019%. Formally our Ks-band transits are deeper than the J-band transits observed simultaneously by a factor of 1.072\pm0.018 - a 4-sigma discrepancy. The most straightforward explanation for our deeper Ks-band depth is a spectral absorption feature from the limb of the atmosphere of the planet; for the spectral absorption feature to be this prominent the atmosphere of GJ 1214b must have a large scale height and a low mean molecular weight. That is, it would have to be hydrogen/helium dominated and this planet would be better described as a mini-Neptune. However, recently published observations from 0.78 - 1.0 microns, by Bean and collaborators, show a lack of spectral features and transit depths consistent with those obtained by Charbonneau and collaborators. The most likely atmospheric composition for GJ 1214b that arises from combining all these observations is less clear; if the atmosphere of GJ 1214b is hydrogen/helium dominated then it must have either a haze layer that is obscuring transit depth differences at shorter wavelengths, or significantly different spectral features than current models predict. Our observations disfavour a water-world composition, but such a composition will remain a possibility until observations reconfirm our deeper Ks-band transit depth or detect features at other wavelengths. [Abridged]

We present Ks, H & J-band photometry of the very highly irradiated hot Jupiter WASP-12b using the Wide-field Infrared Camera on the Canada-France-Hawaii telescope. Our photometry brackets the secondary eclipse of WASP-12b in the Ks and H-bands, and in J-band starts in mid-eclipse and continues until well after the end of the eclipse. We detect its thermal emission in all three near-infrared bands. Our secondary eclipse depths are 0.309 +/- 0.013% in Ks-band (24-sigma), 0.176 +/- 0.020% in H-band (9-sigma) and 0.131 +/- 0.028% in J-band (4-sigma). All three secondary eclipses are best-fit with a consistent phase that is compatible with a circular orbit. By combining our secondary eclipse times with others published in the literature, as well as the radial velocity and transit timing data for this system, we show that there is no evidence that WASP-12b is precessing at a detectable rate, and show that its orbital eccentricity is likely zero. Our thermal emission measurements also allow us to constrain the characteristics of the planet's atmosphere; our Ks-band eclipse depth argues in favour of inefficient day to nightside redistribution of heat and a low Bond albedo for this very highly irradiated hot Jupiter. The J and H-band brightness temperatures are slightly cooler than the Ks-band brightness temperature, and thus hint at the possibility of a modest temperature inversion deep in the atmosphere of WASP-12b; the high pressure, deep atmospheric layers probed by our J and H-band observations are likely more homogenized than the higher altitude layer. Lastly, our best-fit Ks-band eclipse has a marginally longer duration than would otherwise be expected; this may be tentative evidence for material being tidally stripped from the planet - as was predicted for this system by Li & collaborators, and for which observational confirmation was recently arguably provided by Fossati & collaborators.

We study the kinematically narrow, low-ionization line emission from a bright, starburst galaxy at z = 0.69 using slit spectroscopy obtained with Keck/LRIS. The spectrum reveals strong absorption in MgII and FeII resonance transitions with Doppler shifts of -200 to -300 km/s, indicating a cool gas outflow. Emission in MgII near and redward of systemic velocity, in concert with the observed absorption, yields a P Cygni-like line profile similar to those observed in the Ly alpha transition in Lyman Break Galaxies. Further, the MgII emission is spatially resolved, and extends significantly beyond the emission from stars and HII regions within the galaxy. Assuming the emission has a simple, symmetric surface brightness profile, we find that the gas extends to distances > ~7 kpc. We also detect several narrow FeII* fine-structure lines in emission near the systemic velocity, arising from energy levels which are radiatively excited directly from the ground state. We suggest that the MgII and FeII* emission is generated by photon scattering in the observed outflow, and emphasize that this emission is a generic prediction of outflows. These observations provide the first direct constraints on the minimum spatial extent and morphology of the wind from a distant galaxy. Estimates of these parameters are crucial for understanding the impact of outflows in driving galaxy evolution.

The increasing abundance of passive "red-sequence" galaxies since z=1-2 is mirrored by a coincident rise in the number of galaxies with spheroidal morphologies. In this paper, however, we show that in detail the correspondence between galaxy morphology and color is not perfect, providing insight into the physical origin of this evolution. Using the COSMOS survey, we study a significant population of red sequence galaxies with disk-like morphologies. These passive disks typically have Sa-Sb morphological types with large bulges, but they are not confined to dense environments. They represent nearly one-half of all red-sequence galaxies and dominate at lower masses (log Mstar < 10) where they are increasingly disk-dominated. As a function of time, the abundance of passive disks with log Mstar < 11 increases, but not as fast as red-sequence spheroidals in the same mass range. At higher mass, the passive disk population has declined since z~1, likely because they transform into spheroidals. We estimate that as much as 60% of galaxies transitioning onto the red sequence evolve through a passive disk phase. The origin of passive disks therefore has broad implications for understanding how star formation shuts down. Because passive disks tend to be more bulge-dominated than their star-forming counterparts, a simple fading of blue disks does not fully explain their origin. We explore several more sophisticated explanations, including environmental effects, internal stabilization, and disk regrowth during gas-rich mergers. While previous work has sought to explain color and morphological transformations with a single process, these observations open the way to new insight by highlighting the fact that galaxy evolution may actually proceed through several separate stages.

We use a sample of the 13 most luminous WMAP Galactic free-free sources, responsible for 33% of the free- free emission of the Milky Way, to investigate star formation. The sample contains 40 star forming complexes; we combine this sample with giant molecular cloud (GMC) catalogs in the literature, to identify the host GMCs of 32 of the complexes. We estimate the star formation efficiency epsilon_GMC and star formation rate per free-fall time epsilon_ff. We find that epsilon_GMC ranges from 0.002 to 0.2, with an ionizing luminosity-weighted average epsilon_GMC = 0.08, compared to the Galactic average = 0.005. Turning to the star formation rate per free-fall time, we find values that range up to epsilon_ff = 1. Weighting by ionizing luminosity, we find an average of epsilon_ff = 0.16 - 0.24 depending on the estimate of the age of the system. Once again, this is much larger than the Galaxy-wide average value epsilon_ff = 0.008. We show that the lifetimes of giant molecular clouds at the mean mass found in our sample is 17 plus or minus 4 Myr, about two free-fall times. The GMCs hosting the most luminous clusters are being disrupted by those clusters. Accordingly, we interpret the range in epsilon_ff as the result of a time-variable star formation rate; the rate of star formation increases with the age of the host molecular cloud, until the stars disrupt the cloud. These results are inconsistent with the notion that the star formation rate in Milky Way GMCs is determined by the properties of supersonic turbulence

Galactic outflows of low ionization, cool gas are ubiquitous in local starburst galaxies, and in the majority of galaxies at high redshift. How these cool outflows arise is still in question. Hot gas from supernovae has long been suspected as the primary driver, but this mechanism suffers from its tendency to destroy the cool gas as the latter is accelerated. We propose a modification of the supernova scenario that overcomes this difficulty. Star formation is observed to take place in clusters; in a given galaxy, the bulk of the star formation is found in the ~20 most massive clusters. We show that, for L* galaxies, the radiation pressure from clusters with M>10^6 M_sun is able to expel the surrounding gas at velocities in excess of the circular velocity of the disk galaxy. This cool gas can travel above the galactic disk in less than 2 Myr, well before any supernovae erupt in the driving cluster. Once above the disk, the cool outflowing gas is exposed to radiation, and supernovae induced hot gas outflows, from other clusters in the disk, which drive it to distances of several tens to hundreds of kpc. Because the radiatively driven clouds grow in size as they travel, and because the hot gas is more dilute at large distance, the clouds are less subject to destruction if they do eventually encounter hot gas. Therefore, unlike wind driven clouds, radiatively driven clouds can survive to distances ~50 kpc. We identify these cluster-driven winds with large-scale galactic outflows. Another implication of our model is that only starburst galaxies, where massive clusters reside, are able to drive winds cold outflows on galactic scales via this mechanism. We find that the critical star formation rates above which large scale cool outflows will be launched to be ~0.1 M_sun/yr/kpc^2, which is in good agreement with observations.

We examine the thirteen most luminous sources in the WMAP free-free map using the Spitzer GLIMPSE and MSX surveys to identify massive star formation regions, emitting one-third of the Galactic free-free luminosity. We identify star forming regions by a combination of bubble morphology in 8 $\micronm$ (PAH) emission and radio recombination line radial velocities. We find 40 star forming regions associated with our WMAP sources, and determine unique distances to 31. We interpret the bubbles as evidence for radial expansion. The radial velocity distribution for each source allows us to measure the intrinsic speed of a region's expansion. This speed is consistent with the size and age of the bubbles. The high free-free luminosities, combined with negligible synchrotron emission, demonstrate that the bubbles are not driven by supernovae. The kinetic energy of the largest bubbles is a substantial fraction of that measured in the older superbubbles found by Heiles. We find that the energy injected into the ISM by our bubbles is similar to that required to maintain the turbulent motion in the gas disk inside 8 kpc. We report a number of new star forming regions powered by massive ($\textrm{M}_{*} > 10^4 \textrm{M}_\sun$) star clusters. We measure the scale height of the Galactic O stars to be $h_{\textrm{*}} = 35 \pm 5 \pc$. We determine an empirical relationship between the PAH and free-free emission of the form $F_{\textrm{PAH}} \propto F^2_{\textrm{ff}}$. Finally, we find that the bubble geometry is more consistent with a spherical shell rather than a flattened disk.

We compile observations of the surface mass density profiles of dense stellar systems, including globular clusters in the Milky Way and nearby galaxies, massive star clusters in nearby starbursts, nuclear star clusters in dwarf spheroidals and late-type disks, ultra-compact dwarfs, and galaxy spheroids spanning the range from low-mass cusp bulges and ellipticals to massive core ellipticals. We show that in all cases the maximum stellar surface density attained in the central regions of these systems is similar, Sigma_max ~ 10^11 M_sun/kpc^2 (~20 g/cm^2), despite the fact that the systems span 7 orders of magnitude in total stellar mass M_star, 5 in effective radius R_e, and have a wide range in effective surface density M_star/R_e^2. The surface density limit is reached on a wide variety of physical scales in different systems and is thus not a limit on three-dimensional stellar density. Given the very different formation mechanisms involved in these different classes of objects, we argue that a single piece of physics likely determines Sigma_max. The radiation fields and winds produced by massive stars can have a significant influence on the formation of both star clusters and galaxies, while neither supernovae nor black hole accretion are important in star cluster formation. We thus conclude that feedback from massive stars likely accounts for the observed Sigma_max, plausibly because star formation reaches an Eddington-like flux that regulates the growth of these diverse systems. This suggests that current models of galaxy formation, which focus on feedback from supernovae and active galactic nuclei, are missing a crucial ingredient.

Supermassive black holes (BHs) obey tight scaling relations between their mass and their host galaxy properties such as total stellar mass, velocity dispersion, and potential well depth. This has led to the development of self-regulated models for BH growth, in which feedback from the central BH halts its own growth upon reaching a critical threshold. However, models have also been proposed in which feedback plays no role: so long as a fixed fraction of the host gas supply is accreted, relations like those observed can be reproduced. Here, we argue that the scatter in the observed BH-host correlations, and its run with scale, presents a demanding constraint on any model for these correlations, and that it favors self-regulated models of BH growth. We show that the scatter in the stellar mass fraction within a radius R in observed ellipticals and spheroids increases strongly at small R. At fixed total stellar mass (or host velocity dispersion), on very small scales near the BH radius of influence, there is an order-of-magnitude scatter in the amount of gas that must have entered and formed stars. In short, the BH appears to 'know more' about the global host galaxy potential on large scales than the stars and gas supply on small scales. This is predicted in self-regulated models; however, models where there is no feedback would generically predict order-of-magnitude scatter in the BH-host correlations. Likewise, models in which the BH feedback in the 'bright' mode does not regulate the growth of the BH itself, but sets the stellar mass of the galaxy by inducing star formation or blowing out a mass in gas much larger than the galaxy stellar mass, are difficult to reconcile with the scatter on small scales.

Observations suggest that effective radii of high-z massive spheroids are as much as a factor ~6 smaller than low-z galaxies of comparable mass. Given the apparent absence of low-z counterparts, this has often been interpreted as indicating that the high density, compact red galaxies must be 'puffed up' by some mechanism. We compare the ensemble of high-z observations with large samples of well-observed low-z ellipticals. At the same physical radii, the stellar surface mass densities of low and high-z systems are comparable. Moreover, the abundance of high surface density material at low redshift is comparable to or larger than that observed at z>1-2, consistent with the continuous buildup of spheroids over this time. The entire population of compact, high-z red galaxies may be the progenitors of the high-density cores of present-day ellipticals, with no need for a decrease in stellar density from z=2 to z=0. The primary difference between low and high-z systems is thus the observed low-density material at large radii in low-z spheroids (rather than the high-density material in high-z spheroids). Such low-density material may either (1) assemble at z<2 or (2) be present, but not yet detected, at z>2. Mock observations of low-z massive systems show that the high-z observations do not yet probe sufficiently low surface brightness material to detect the low surface density 'wings' (if present). Thus, if the high-z galaxies resemble the most massive systems today, their inferred effective radii could be under-estimated by factors ~2-4. This difference arises because massive systems at low redshift are not well-fit by single Sersic profiles. We discuss implications of our results for physical models of galaxy evolution.

Star formation is slow, in the sense that the gas consumption time is much longer than the dynamical time. It is also inefficient; essentially all star formation in local galaxies takes place in giant molecular clouds (GMCs), but the fraction of a GMC converted to stars is very small, ~5%. In the most luminous starbursts, the GMC lifetime is shorter than the main sequence lifetime of even the most massive stars, so that supernovae can play no role in GMC disruption. We investigate the disruption of GMCs across a wide range of galaxies, from normal spirals to the densest starbursts; we take into account the effects of HII gas pressure, shocked stellar winds, protostellar jets, and radiation pressure produced by the absorption and scattering of starlight on dust grains. In the Milky Way, we find that a combination of three mechanisms, jets, HII gas pressure, and radiation pressure, disrupts the clouds. In more rapidly star forming galaxies such as ``clump'' galaxies at high-redshift, ultra-luminous infrared galaxies (ULIRGs) and submillimeter galaxies, radiation pressure dominates natal cloud distribution. We predict the presence of 10-20 clusters with masses ~10^7 Msun in local ULIRGs such as Arp 220 and a similar number of clusters with M_* ~ 10^8 Msun in high redshift clump galaxies; submillimeter galaxies will have even more massive clusters. We find that the mass fraction of a GMC that ends up in stars is an increasing function of the gas surface density of a galaxy, reaching ~35% in the most luminous starbursts. Furthermore, the disruption of bubbles by radiation pressure stirs the interstellar medium to velocities of ~10 km/s in normal galaxies and to ~100 km/s in ULIRGs like Arp 220, consistent with observations. Thus, radiation pressure may play a dominant role in the ISM of star-forming galaxies.

We use the WMAP maximum entropy method foreground emission map combined with previously determined distances to giant HII regions to measure the free-free flux at Earth and the free-free luminosity of the galaxy. We find a total flux f_\nu=54211 Jy and a flux from 88 sources of f_\nu=36043 Jy. The bulk of the sources are at least marginally resolved, with mean radii ~60 pc, electron density n_e ~ 9 cm^{-3}, and filling factor \phi_{HII}=0.005 (over the Galactic gas disk). The total dust-corrected ionizing photon luminosity is Q=3.2x10^{53} photons/s, in good agreement with previous estimates. We use GLIMPSE and MSX 8 micron images to show that the bulk of the free-free luminosity is associated with bubbles having radii r~5-100 pc, with a mean ~20 pc. These bubbles are leaky, so that ionizing photons from inside the bubble excite free-free emission beyond the bubble walls, producing WMAP sources that are larger than the 8 micron bubbles. We suggest that the WMAP sources are the counterparts of the extended low density HII regions described by Mezger (1978). Half the ionizing luminosity from the sources is emitted by the nine most luminous objects, while the seventeen most luminous emit half the total Galactic ionizing flux. These 17 sources have 4x10^{51} < Q <1.6x10^{52}, corresponding to 6x10^4M_\odot < M_*< 2x10^5M_\odot; half to two thirds of this will be in the central massive star cluster. We convert the measurement of Q to a Galactic star formation rate dM/dt=1.3M_\odot/yr, but point out that this is highly dependent on the exponent \Gamma~1.35 of the high mass end of the stellar initial mass function.

Using observations from the literature, we show that the non-thermal radio luminosity (L) of supernova remnants (SNRs) is a strong function of the average gas surface density (Sigma) of the galaxy in which the remnants reside, from normal spirals to luminous starbursts. We combine a simple theory for electron cooling in SNRs with the observed radio luminosities to estimate the remnant magnetic field strength (B_SNR): the correlation between L and Sigma implies that B_SNR also increases with Sigma. We explore two interpretations of this correlation: (1) B_SNR is generated by post-shock magnetic field amplification, with B_SNR^2 proportional to Sigma and (2) B_SNR results from shock-compression of the ambient ISM magnetic field (B_ISM), with B_ISM being larger in denser galaxies. We find that shock compression is, on average, sufficient to produce the observed radio emission from SNRs in the densest starbursts; amplification of post-shock magnetic fields is not required. By contrast, in normal spirals post-shock field amplification (by a factor of a few - 10) is consistent with the data; we find tentative evidence that both the Alfven speed and the ratio of B_SNR^2 to the post-shock pressure ("epsilon_B") are constant in SNRs from galaxy to galaxy. We discuss observational tests that can be used to distinguish between these two interpretations of the radio luminosities of SNRs. Regardless of which is correct, the radio emission from SNRs provides an upper limit to B_ISM that is independent of the minimum energy assumption. For the densest starbursts, the ISM magnetic energy density is below that required for hydrostatic equilibrium; thus magnetic fields are not dynamically important on the largest scales in starbursts, in contrast with spiral galaxies like our own. This dichotomy may have implications for galactic dynamo theory.

We have tested the two main theoretical models of bubbles around massive star clusters, Castor et al. and Chevalier & Clegg, against observations of the well studied Carina Nebula. The Castor et al. theory over-predicts the X-ray luminosity in the Carina bubble by a factor of 60 and expands too rapidly, by a factor of 4; if the correct radius and age are used, the predicted X-ray luminosity is even larger. In contrast, the Chevalier & Clegg model under-predicts the X-ray luminosity by a factor of 10. We modify the Castor et al. theory to take into account lower stellar wind mass loss rates, radiation pressure, gravity, and escape of or energy loss from the hot shocked gas. We argue that energy is advected rather than radiated from the bubble. We undertake a parameter study for reduced stellar mass loss rates and for various leakage rates and are able to find viable models. The X-ray surface brightness in Carina is highest close to the bubble wall, which is consistent with conductive evaporation from cold clouds. The picture that emerges is one in which the hot gas pressure is far below that found by dividing the time-integrated wind luminosity by the bubble volume; rather, the pressure in the hot gas is set by pressure equilibrium with the photoionized gas at T=10^4 K. It follows that the shocked stellar winds are not dynamically important in forming the bubbles.

We describe a newly developed hydrodynamic code for studying accretion disk processes. The numerical method uses a finite volume, nonlinear, Total Variation Diminishing (TVD) scheme to capture shocks and control spurious oscillations. It is second-order accurate in time and space and makes use of a FARGO-type algorithm to alleviate Courant-Friedrichs-Lewy time step restrictions imposed by the rapidly rotating inner disk region. OpenMP directives are implemented enabling faster computations on shared-memory, multi-processor machines. The resulting code is simple, fast and memory efficient. We discuss the relevant details of the numerical method and provide results of the code's performance on standard test problems. We also include a detailed examination of the code's performance on planetary disk-planet interactions. We show that the results produced on the standard problem setup are consistent with a wide variety of other codes.

Photoionization heating from UV radiation incident on the atmospheres of hot Jupiters may drive planetary mass loss. We construct a model of escape that includes realistic heating and cooling, ionization balance, tidal gravity, and pressure confinement by the host star wind. We show that mass loss takes the form of a hydrodynamic ("Parker") wind, emitted from the planet's dayside during lulls in the stellar wind. When dayside winds are suppressed by the confining action of the stellar wind, nightside winds might pick up if there is sufficient horizontal transport of heat. A hot Jupiter loses mass at maximum rates of ~2 x 10^12 g/s during its host star's pre-main-sequence phase and ~2 x10^10 g/s during the star's main sequence lifetime, for total maximum losses of ~0.06% and ~0.6% of the planet's mass, respectively. For UV fluxes F_UV < 10^4 erg/cm^2/s, the mass loss rate is approximately energy-limited and is proportional to F_UV^0.9. For larger UV fluxes, such as those typical of T Tauri stars, radiative losses and plasma recombination force the mass loss rate to increase more slowly as F_UV^0.6. Dayside winds are quenched during the T Tauri phase because of confinement by overwhelming stellar wind pressure. We conclude that while UV radiation can indeed drive winds from hot Jupiters, such winds cannot significantly alter planetary masses during any evolutionary stage. They can, however, produce observable signatures. Candidates for explaining why the Lyman-alpha photons of HD 209458 are absorbed at Doppler-shifted velocities of +/- 100 km/s include charge-exchange in the shock between the planetary and stellar winds.

The masses of star clusters range over seven decades, from ten up to one hundred million solar masses. Remarkably, clusters with masses in the range 10^4 to 10^6 solar mases show no systematic variation of radius with mass. However, recent observations have shown that clusters with masses greater than 3x10^6 solar masses do show an increase in size with increasing mass. We point out that clusters with m>10^6 solar masses were optically thick to far infrared radiation when they formed, and explore the hypothesis that the size of clusters with m> 3x10^6 solar masses is set by a balance between accretion powered radiation pressure and gravity when the clusters formed, yielding a mass-radius relation r~0.3(m/10^6M_\odot)^{3/5} pc. We show that the Jeans mass in optically thick objects increases systematically with cluster mass. We argue, by assuming that the break in the stellar initial mass function is set by the Jeans mass, that optically thick clusters are born with top heavy initial mass functions; it follows that they are over-luminous compared to optically thin clusters when young, and have a higher mass to light ratio Upsilon_V=m/L_V when older than ~1 Gyr. Old, optically thick clusters have Upsilon_V~ mcl^{0.1-0.3}. It follows that L_V~\sigma^{\beta}, where \sigma is the cluster velocity dispersion, and \beta~4. It appears that Upsilon_V is an increasing function of cluster mass for compact clusters and ultra-compact dwarf galaxies. We show that this is unlikely to be due to the presence of non-baryonic dark matter, by comparing clusters to Milky Way satellite galaxies, which are dark matter dominated. The satellite galaxies appear to have a fixed mass inside a fiducial radius, M(r=r_0)=const.

Radial velocity surveys find Jupiter mass planets with semi-major axes a less than 0.1 AU around ~1% of solar-type stars; counting planets with $a$ as large as 5 AU, the fraction of stars having planets reaches ~ 10% {Marcy,Butler}. An examination of the distribution of semi-major axes shows that there is a clear excess of planets with orbital periods around 3 or 4 days, corresponding to a~0.03$ AU, with a sharp cutoff at shorter periods (see Figure 1). It is believed that Jupiter mass planets form at large distances from their parent stars; some fraction then migrate in to produce the short period objects. We argue that a significant fraction of the `hot Jupiters' (a<0.1 AU) may arise in binary star systems in which the orbit of the binary is highly inclined to the orbit of the planet. Mutual torques between the two orbits drive down the minimum separation or periapse r_p between the planet and its host star (the Kozai mechanism). This periapse collapse is halted when tidal friction on the planet circularizes the orbit faster than Kozai torque can excite it. The same friction then circularizes the planet orbit, producing hot Jupiters with the peak of the semimajor axis distribution lying around 3 days. For the observed distributions of binary separation, eccentricity and mass ratio, roughly 2.5% of planets with initial semimajor axis a_p ~ 5au will migrate to within 0.1au of their parent star. Kozai migration could account for 10% or more of the observed hot Jupiters.

We argue that the X-ray and UV flux illuminating the parsec-scale accretion disk around luminous active galactic nuclei (AGN) is super-Eddington with respect to the local far-infrared dust opacity. The far infrared opacity may be larger than in the interstellar medium of the Milky Way due to a combination of supersolar metallicity and the growth of dust grains in the dense accretion disk. Because of the irradiating flux, the outer accretion disk puffs up with a vertical thickness $h\sim R$. This provides a mechanism for generating a geometrically thick obscuring region from an intrinsically thin disk. We find obscuring columns $\sim 10^{22} - 10^{23} {\rm cm}^{-2}$, in reasonable agreement with observations.

In the core accretion model, gas giant formation is a race between growth and migration; for a core to become a jovian planet, it must accrete its envelope before it spirals into the host star. We use a multizone numerical model to extend our previous investigation of the "window of opportunity" for gas giant formation within a disk. When the collision cross-section enhancement due to core atmospheres is taken into account, we find that a broad range of protoplanetary disks posses such a window.

We build and test Parker-wind models to apply to observations of large-scale (of order 100 pc) outflows from Active Galactic Nuclei (AGNs). These models include detailed photoionization simulations, the observed radially varying mass profile, adiabatic cooling, and approximations for clouds dragged along in the wind and the interaction of the wind with the circumnuclear ISM of the galaxy. We test this model against recent HST/STIS observations of [O III] emission-line kinematics (in particular, we test against those observed in NGC 4151, but approximately the same kinematics is observed in NGC 1068 and Mrk 3) to constrain the viability of large-scale thermal winds in AGNs. We find that adiabatic cooling dominates in these outflows, decelerating Parker winds on large scales, making them highly unlikely as explanations of the observed kinematics.

Roughly 80% of Ultraluminous Infrared Galaxies (ULIRGs) show blue shifted absorption in the resonance lines of neutral sodium, indicating that cool winds are common in such objects, as shown by Rupke et al and by Martin. The neutral sodium (NaI) columns indicated by these absorption lines are ~ 10^{13}-3x10^{14}/cm^2, while the bolometric luminosity varies by a factor of only four. We show that the gas in ULIRG outflows is likely to be in photoionization equilibrium. The very small ULIRG sample of Goldader et al. demonstrates that the ratio of ultraviolet flux to far infrared flux varies by a factor $\sim100$ from object to object. While the Goldader sample does not overlap with those of Rupke et al. and Martin, we show that such a large variation in ultraviolet flux will produce a similar variation in the column of neutral sodium for a fixed mass flux and density. However, if the cold gas is in pressure equilibrium with a hot outflow with a mass loss rate similar to the star formation rate, the range of ionization state is significantly smaller. Measurements of the UV flux for objects in the Martin and Rupke et al. catalogs will definitively determine if photoionization effects are responsible for the wide variation seen in the sodium columns. If they are, a determination of the gas density and mass loss rate in the cool winds will follow, with attendant improvements in our understanding of wind driving mechanisms and of the effects of galaxies on their surroundings.

In the standard model of gas giant planet formation, a large solid core (~ 10 times the Earth's mass) forms first, then accretes its massive envelope (100 or more Earth masses) of gas. However, inward planet migration due to gravitational interaction with the proto-stellar gas disk poses a difficulty in this model. Core-sized bodies undergo rapid "Type I" migration; for typical parameters their migration timescale is much shorter than their accretion timescale. How, then, do growing cores avoid spiraling into the central star before they ever get the chance to become gas giants? Here, we present a simple model of core formation in a gas disk which is viscously evolving. As the disk dissipates, accretion and migration timescales eventually become comparable. If this happens while there is still enough gas left in the disk to supply a jovian atmosphere, then a window of opportunity for gas giant formation opens. We examine under what circumstances this happens, and thus, what predictions our model makes about the link between proto-stellar disk properties and the likelihood of forming giant planets.

As a result of deep hard X-ray observations by Chandra and XMM-Newton a significant fraction of the cosmic X-ray background (CXRB) has been resolved into individual sources. These objects are almost all active galactic nuclei (AGN) and optical followup observations find that they are mostly obscured Type 2 AGN, have Seyfert-like X-ray luminosities (i.e., L_X ~ 10^{43-44} ergs s^{-1}), and peak in redshift at z~0.7. Since this redshift is similar to the peak in the cosmic star-formation rate, this paper proposes that the obscuring material required for AGN unification is regulated by star-formation within the host galaxy. We test this idea by computing CXRB synthesis models with a ratio of Type 2/Type 1 AGN that is a function of both z and 2-10 keV X-ray luminosity, L_X. The evolutionary models are constrained by parameterizing the observed Type 1 AGN fractions from the recent work by Barger et al. The parameterization which simultaneously best accounts for Barger's data, the CXRB spectrum and the X-ray number counts has a local, low-L_X Type 2/Type 1 ratio of 4, and predicts a Type 2 AGN fraction which evolves as (1+z)^{0.3}. Models with no redshift evolution yielded much poorer fits to the Barger Type 1 AGN fractions. This particular evolution predicts a Type 2/Type 1 ratio of 1-2 for log L_X > 44, and thus the deep X-ray surveys are missing about half the obscured AGN with these luminosities. These objects are likely to be Compton thick. Overall, these calculations show that the current data strongly supports a change to the AGN unification scenario where the obscuration is connected with star formation in the host galaxy rather than a molecular torus alone. The evolution of the obscuration implies a close relationship between star formation and AGN fueling, most likely due to minor mergers or interactions.

We estimate minimum energy magnetic fields (B_min) for a sample of galaxies with measured gas surface densities, spanning from normal spirals to starbursts. We show that the ratio of the minimum energy magnetic pressure to the total pressure in the ISM decreases substantially with increasing surface density; for Arp 220 this ratio is ~10^-4. Therefore, if the minimum energy estimate is applicable, magnetic fields in starbursts are dynamically weak compared to gravity, in contrast to normal spiral galaxies. We argue, however, that rapid cooling of relativistic electrons in starbursts invalidates the minimum energy estimate. We critically assess a number of independent constraints on the magnetic field strength in starbursts. In particular, we argue that the existence of the FIR-radio correlation implies that the synchrotron cooling timescale for cosmic ray electrons is much shorter than their escape time from the galactic disk; this in turn implies that the true magnetic field in starbursts is significantly larger than B_min. The strongest argument against such large fields is that one might expect starbursts to have steep radio spectra indicative of strong synchrotron cooling, which is not observed. We show, however, that ionization and bremsstrahlung losses can flatten the nonthermal spectra of starburst galaxies even in the presence of rapid cooling, providing much better agreement with observed spectra. We further demonstrate that ionization and bremsstrahlung losses are likely to be important in shaping the radio spectra of most starbursts at GHz frequencies, thereby preserving the linearity of the FIR-radio correlation. We thus conclude that magnetic fields in starbursts are significantly larger than B_min. We highlight several observations that can test this conclusion.

Oxygen abundances have been derived from the near-IR, high-excitation Lambda 7774 O I triplet in high-resolution, high signal-to-noise spectra of 45 Hyades dwarfs using standard one dimensional, plane-parallel LTE models. Effective temperatures of the stellar sample range from 4319-6301 K, and the derived relative O abundances as a function of T_eff evince a trichotomous morphology. At T_eff > 6100 K, there is evidence of an increase in the O abundances with increasing T_eff, consistent with non-LTE (NLTE) predictions. At intermediate T_eff (5450 < T_eff < 6100 K), the O abundances are flat, and star-to-star values are in good agreement, having a mean value of [O/H] = +0.25 +/- 0.02; however, systematic errors at the ~0.10 dex level might exist. The O abundances for stars with T_eff < 5450 K show a striking increase with decreasing T_eff, in stark contrast to expectations and canonical NLTE calculations. The cool Hyades triplet results are compared to those recently reported for dwarfs in the Pleiades cluster and the UMa moving group; qualitative differences between the trends observed in these stellar aggregates point to a possible age-related diminution of triplet abundance trends in cool open cluster dwarfs. Correlations with age-related phenomena, i.e., chromospheric activity and photospheric spots, faculae, and/or plages, are investigated. No correlation with Ca II H+K chromospheric activity indicators is observed. Multi-component LTE ``toy'' models have been constructed in order to simulate photospheric temperature inhomogeneities that could arise from the presence of starspots, and we demonstrate that photospheric spots are a plausible source of the triplet trends among the cool dwarfs.

We consider the structure of marginally Toomre-stable starburst disks under the assumption that radiation pressure on dust grains provides the dominant vertical support against gravity. This is particularly appropriate when the disk is optically thick to its own IR radiation, as in the central regions of ULIRGs. Because the disk radiates at its Eddington limit, the Schmidt-law for star formation changes in the optically-thick limit, with the star formation rate per unit area scaling as Sigma_g/kappa, where Sigma_g is the gas surface density and kappa is the mean opacity. We show that optically thick starburst disks have a characteristic flux and dust effective temperature of F ~ 10^{13} L_sun/kpc^2 and T_eff ~ 90K, respectively. We compare our predictions with observations and find good agreement. We extend our model from many-hundred parsec scales to sub-parsec scales and address the problem of fueling AGN. We assume that angular momentum transport proceeds via global torques rather than a local viscosity. We account for the radial depletion of gas due to star formation and find a strong bifurcation between two classes of disk models: (1) solutions with a starburst on large scales that consumes all of the gas with little fueling of a central AGN and (2) models with an outer large-scale starburst accompanied by a more compact starburst on 1-10 pc scales and a bright central AGN. The luminosity of the latter models is in many cases dominated by the AGN. We show that the vertical thickness of the starburst disk on pc scales can approach h ~ r, perhaps accounting for the nuclear obscuration in some Type 2 AGN. We also argue that the disk of young stars in the Galactic Center may be the remnant of such a compact nuclear starburst.

The two planets about the star GJ 876 appear to have undergone extensive migration from their point of origin in the protoplanetary disk -- both because of their close proximity to the star (30 and 60 day orbital periods) and because of their occupying three stable orbital resonances at the 2:1 mean-motion commensurability. The resonances were most likely established by converging differential migration of the planets leading to capture into the resonances. A problem with this scenario is that continued migration of the system while it is trapped in the resonances leads to orbital eccentricities that rapidly exceed the observational upper limits of e_1 = 0.31 and e_2 = 0.05. As seen in forced 3-body simulations, lower eccentricities would persist during migration only for an applied eccentricity damping. Here we explore the evolution of the GJ 876 system using two-dimensional hydrodynamical simulations that include viscous heating and radiative effects. We find that a hydrodynamic evolution within the resonance, where only the outer planet interacts with the disk, always rapidly leads to large values of eccentricities that exceed those observed. Only if mass is removed from the disk on a time scale of the order of the migration time scale (before there has been extensive migration after capture), as might occur for photoevaporation in the late phases of planet formation, can we end up with eccentricities that are consistent with the observations.

Future surveys for transiting extrasolar planets, including the space-based mission Kepler (Borucki et al 2003), are expected to detect hundreds of Jovian mass planets and tens of terrestrial mass planets. For many of these newly discovered planets, the intervals between successive transits will be measured with an accuracy of 0.1--100 minutes. We show that these timing measurements will allow for the detection of additional planets in the system (not necessarily transiting), via their gravitational interaction with the transiting planet. The transit time variations depend on the mass of the additional planet, and in some cases Earth-mass planets will produce a measurable effect.

We investigate large-scale galactic winds driven by momentum deposition. Momentum injection is provided by (1) radiation pressure produced by the continuum absorption and scattering of UV photons on dust grains and (2) supernovae. UV radiation can be produced by a starburst or AGN activity. We argue that momentum-driven winds are an efficient mechanism for feedback during the formation of galaxies. We show that above a limiting luminosity, momentum deposition from star formation can expel a significant fraction of the gas in a galaxy. The limiting, Eddington-like luminosity is $L_{\rm M}\simeq(4f_g c/G) \sigma^4$, where $\sigma$ is the galaxy velocity dispersion and $f_g$ is the gas fraction. A starburst that attains $L_{\rm M}$ moderates its star formation rate and its luminosity does not increase significantly further. We argue that ellipticals attain this limit during their growth at $z \gtrsim 1$ and that this is the origin of the Faber-Jackson relation. We show that Lyman break galaxies and ultra-luminous infrared galaxies have luminosities near $L_{\rm M}$. Star formation is unlikely to efficiently remove gas from very small scales in galactic nuclei, i.e., scales much smaller than that of a nuclear starburst. This gas is available to fuel a central black hole (BH). We argue that a BH clears gas out of its galactic nucleus when the luminosity of the BH itself reaches $\approx L_{\rm M}$. This shuts off the fuel supply to the BH and may also terminate star formation in the surrounding galaxy. As a result, the BH mass is fixed to be $M_{\rm BH}\simeq (f_g \kappa_{\rm es}/\pi G^2)\sigma^4$, where $\kappa_{\rm es}$ is the electron scattering opacity. This limit is in accord with the observed $M_{\rm BH}-\sigma$ relation. (Abridged)

Using high-resolution, high quality spectra we investigate the presence of $^6$Li in two lithium-poor stars that host extrasolar planetary systems. We present improved atomic and molecular line lists for the region in the vicinity of the lithium line at 6707.8 \AA, and we produce an excellent fit to the solar spectrum. From line profile fitting, we find results consistent with no $^6$Li in either of the lithium-poor planet-bearing stars or in three comparison stars with and without planets, and 1-$\sigma$ upper limits of 0.04 for the isotopic ratios of the two lithium-poor stars give an upper limit of 0.3 Jupiter masses of material with primordial abundances that could have been recently deposited in their outer layers. Our results suggest that post-main sequence accretion of planets or planetary material that is undepleted in lithium is uncommon.

Micron size extrasolar dust particles have been convincingly detected by satellites. Larger extrasolar meteoroids (5-35 microns) have most likely been detected by ground based radar at Arecibo and New Zealand. We present estimates of the minimum detectable particle sizes and collecting areas for both radar systems. We show that particles larger than about 10 microns can propagate for tens of parsecs through the interstellar medium, opening up the possibility that ground based radar systems can detect AGB stars, young stellar objects such as T Tauri stars, and debris disks around Vega-like stars. We provide analytical and numerical estimates of the ejection velocity in the case of a debris disk interacting with a Jupiter mass planet. We provide rough estimates of the flux of large micrometeoroids from all three classes of sources. Current radar systems are unlikely to detect significant numbers of meteors from debris disks such as Beta Pictoris. However, we suggest improvements to radar systems that should allow for the detection of multiple examples of all three classes.

Hubble Space Telescope images of MCG-6-30-15 show a dust lane crossing the galaxy just below the nucleus. In this paper, we argue that this dust lane is responsible for the observed reddening of the nuclear emission and the Fe I edge hinted at in the Chandra spectrum of MCG-6-30-15. We further suggest that the gas within the dust lane can comprise much of the low ionization component (i.e., the one contributing the O VII edge) of the observed warm absorber. Moreover, placing the warm absorbing material at such distances (hundreds of pc) can account for the small outflow velocities of the low ionization absorption lines as well as the constancy of the O VIII edge. Photoionization models of a dusty interstellar gas cloud (with a column appropriate for the reddening toward MCG-6-30-15) using a toy Seyfert 1 spectral energy distribution show that it is possible to obtain a significant O VII edge (\tau~0.2) if the material is ~150 pc from the ionizing source. For MCG-6-30-15, such a distance is consistent with the observed dust lane. The current data on MCG-6-30-15 is unable to constrain the dust composition within the warm absorber. Astronomical silicate is a viable candidate, but there are indications of a very low O abundance in the dust, which is inconsistent with a silicate origin. If true, this may indicate that there were repeated cycles of grain destruction and growth from shocks in the interstellar medium of MCG-6-30-15. Pure iron grains are an unlikely dust constituent due to the limit on their abundance in the Galaxy, yet they cannot be ruled out. The high column densities inferred from the highly ionized zone of the warm absorber implies that this gas is dust-free.

The exo-solar planet HD 80606b has a highly eccentric (e=0.93) and tight (a=0.47 AU) orbit. We study how it might arrive at such an orbit and how it has avoided being tidally circularized until now. The presence of a stellar companion to the host star suggests the possibility that the Kozai mechanism and tidal dissipation combined to draw the planet inward well after it formed: Kozai oscillations produce periods of extreme eccentricity in the planet orbit, and the tidal dissipation that occurs during these periods of small pericentre distances leads to gradual orbital decay. We call this migration mechanism the 'Kozai migration'. It requires that the initial planet orbit is highly inclined relative to the binary orbit. For a companion at 1000 AU and an initial planet orbit at 5 AU, the minimum relative inclination required is ~ 85. We discuss the efficiency of tidal dissipation inferred from the observations of exo-planets. Moreover, we investigate possible explanations for the velocity residual (after the motion induced by the planet is removed) observed on the host star: a second planet in the system is excluded over a large extent of semi-major axis space if Kozai migration is to work, and the tide raised on the star by HD 80606b is likely too small in amplitude. Lastly, we discuss the relevance of Kozai migration for other planetary systems.

The recent discovery of Jupiter-mass planets orbiting at a few AU from their stars compliments earlier detections of massive planets on very small orbits. The short period orbits strongly suggest that planet migration has occurred, with the likely mechanism being tidal interactions between the planets and the gas disks out of which they formed. The newly discovered long period planets, together with the gas giant planets in our solar system, show that migration is either absent or rapidly halted in at least some systems. We propose a mechanism for halting type-II migration at several AU in a gas disk. Photoevaporation of the disk by irradiation from the central star can produce a gap in the disk at a few AU, preventing planets outside the gap from migrating down to the star. This would result in an excess of systems with planets at or just outside the photoevaporation radius.

Recently, Maiolino et al. (2001a, A&A, 365, 28) constructed a sample of active galactic nuclei for which both the reddening E(B-V) and the column density N_H to the nucleus could be determined. For most of the galaxies in their sample, they found that E(B-V)/N_H is substantially smaller than for the diffuse ISM of our Galaxy. They asserted that either the dust-to-gas ratio is lower than in the Galaxy or that the grains are so large that they do not extinct or redden efficiently in the optical. We show that there is no systematic increase in E(B-V) with N_H for the Maiolino et al. (2001a) galaxies, which suggests that the X-ray absorption and optical extinction occur in distinct media. Maiolino et al. (2001b, A&A, 365, 37) suggested that the observed lines of sight for the Maiolino et al. (2001a) galaxies pass through the ``torus'' that obscures the broad line region and nuclear continuum in Seyfert 2 galaxies and argued that the torus grains are larger than Galactic grains. There is no reason to believe that the lines of sight for these galaxies pass through the torus, since the observed column densities are lower than those typically observed in Seyfert 2 galaxies. We suggest instead that the X-ray absorption occurs in material located off the torus and/or accretion disk while the optical extinction occurs in material located beyond the torus. The X-ray absorbing material could either be dust-free or could contain large grains that do not extinct efficently in the optical. There is no conclusive evidence that the grains in active galactic nuclei are systematically larger than those in the diffuse ISM of our Galaxy.

The orbits of the outer two known planets orbiting Upsilon Andromedae are remarkably eccentric. Planet C possesses an orbital eccentricity of e1 = 0.253. For the more distant planet D, e2 = 0.308. Previous dynamical analyses strongly suggest that the two orbits are nearly co-planar and are trapped in an apsidal resonance in which the difference between their longitudes of periastron undergoes a bounded oscillation about 0 degrees. Here we elucidate the origin of these large eccentricities and of the apsidal alignment. Resonant interactions between a remnant circumstellar disk of gas lying exterior to the orbits of both planets can smoothly grow e2. Secular interactions between planets D and C can siphon off the eccentricity of the former to grow that of the latter. Externally amplifying e2 during the phase of the apsidal oscillation when e2/e1 is smallest drives the oscillation amplitude towards zero. Thus, the substantial eccentricity of planet C and the locking of orbital apsides are both consequences of externally pumping the eccentricity of planet D over timescales exceeding apsidal precession periods of order 1e4 yr. We explain why the recently detected stellar companion to Upsilon Andromedae is largely dynamically decoupled from the planetary system.

We study the structure and properties of hot MHD accretion onto a Kerr black hole. In such a system, the hole is magnetically coupled to the inflowing gas and exerts a torque onto the accretion flow. A hot settling flow can form around the hole and transport the angular momentum outward, to the outer edge of the flow. Unlike other hot flows, such as advection- and convection-dominated flows and inflow-outflow solutions (ADAFs, CDAFs, and ADIOS), the properties of the hot settling flow are determined by the spin of the central black hole, but are insensitive to the mass accretion rate. Therefore, it may be possible to identify rapidly spinning BHs simply from their broad-band spectra. Observationally, the hot settling flow around a Kerr hole is somewhat similar to other hot flows in that they all have hard, power-law spectra and relatively low luminosities. Thus, most black hole candidates in the low/hard and, perhaps, intermediate X-ray state may potentially accrete via the hot settling flow. However, a settling flow will be somewhat more luminous than ADAFs/CDAFs/ADIOS, will exhibit high variability in X-rays, and may have relativistic jets. This suggests that galactic microquasars and active galactic nuclei may be powered by hot settling flows. We identify several galactic X-ray sources as the best candidates.

Our understanding of the Solar System has been revolutionized over the past decade by the finding that the orbits of the planets are inherently chaotic. In extreme cases, chaotic motions can change the relative positions of the planets around stars, and even eject a planet from a system. Moreover, the spin axis of a planet-Earth's spin axis regulates our seasons-may evolve chaotically, with adverse effects on the climates of otherwise biologically interesting planets. Some of the recently discovered extrasolar planetary systems contain multiple planets, and it is likely that some of these are chaotic as well.

The physical basis of chaos in the solar system is now better understood: in all cases investigated so far, chaotic orbits result from overlapping resonances. Perhaps the clearest examples are found in the asteroid belt. Overlapping resonances account for its Kirkwood gaps and were used to predict and find evidence for very narrow gaps in the outer belt. Further afield, about one new ``short-period'' comet is discovered each year. They are believed to come from the ``Kuiper Belt'' (at 40 AU or more) via chaotic orbits produced by mean-motion and secular resonances with Neptune. Finally, the planetary system itself is not immune from chaos. In the inner solar system, overlapping secular resonances have been identified as the possible source of chaos. For example, Mercury, in 10^{12} years, may suffer a close encounter with Venus or plunge into the Sun. In the outer solar system, three-body resonances have been identified as a source of chaos, but on an even longer time scale of 10^9 times the age of the solar system. On the human time scale, the planets do follow their orbits in a stately procession, and we can predict their trajectories for hundreds of thousands of years. That is because the mavericks, with shorter instability times, have long since been ejected. The solar system is not stable; it is just old!

We describe deep (40 orbits) HST/STIS observations of the BALQSO PG 0946+301 and make them available to the community. These observations are the major part of a multi-wavelength campaign on this object aimed at determining the ionization equilibrium and abundances (IEA) in broad absorption line (BAL) QSOs. We present simple template fits to the entire data set, which yield firm identifications for more than two dozen BALs from 18 ions and give lower limits for the ionic column densities. We find that the outflow's metalicity is consistent with being solar, while the abundance ratio of phosphorus to other metals is at least ten times solar. These findings are based on diagnostics that are not sensitive to saturation and partial covering effects in the BALs, which considerably weakened previous claims for enhanced metalicity. Ample evidence for these effects is seen in the spectrum. We also discuss several options for extracting tighter IEA constraints in future analyses, and present the significant temporal changes which are detected between these spectra and those taken by the HST/FOS in 1992.

We compare the metallicities of stars with radial velocity planets to the metallicity of a sample of field dwarfs. We confirm recent work indicating that the stars-with-planet sample as a whole is iron rich. However, the lowest mass stars tend to be iron poor, with several having [Fe/H]<-0.2, demonstrating that high metallicity is not required for the formation of short period Jupiter-mass planets. We show that the average [Fe/H] increases with increasing stellar mass (for masses below 1.25 solar masses) in both samples, but that the increase is much more rapid in the stars-with-planet sample. The variation of metallicity with stellar age also differs between the two samples. We examine possible selection effects related to variations in the sensitivity of radial velocity surveys with stellar mass and metallicity, and identify a color cutoff (B-V>0.48) that contributes to but does not explain the mass-metallicity trend in the stars-with-planets sample. We use Monte Carlo models to show that adding an average of 6.5 Earth masses of iron to each star can explain both the mass-metallicity and the age-metallicity relations of the stars-with-planets sample. However, for at least one star, HD 38529, there is good evidence that the bulk metallicity is high. We conclude that the observed metallicities and metallicity trends are the result of the interaction of three effects; accretion of about 6 Earth masses of iron rich material, selection effects, and in some cases, high intrinsic metallicity.

We consider the amplification and transport of a magnetic field in the collapsed core of a massive star, including both the region between the neutrinosphere and the shock, and the central, opaque core. An analytical argument explains why rapid convective overturns persist within a newly formed neutron star for roughly 10 seconds ($> 10^3$ overturns), consistent with recent numerical models. A dynamical balance between turbulent and magnetic stresses within this convective layer corresponds to flux densities in excess of $10^{15}$G. Material accreting onto the core is heated by neutrinos and also becomes strongly convective. We compare the expected magnetic stresses in this convective `gain layer' with those deep inside the neutron core. Buoyant motions of magnetized fluid are greatly aided by the intense neutrino flux. We calculate the transport rate through a medium containing free neutrons protons, and electrons, in the limiting cases of degenerate or non-degenerate nucleons. Fields stronger than $\sim 10^{13}$ G are able to rise through the outer degenerate layers of the neutron core during the last stages of Kelvin-Helmholtz cooling (up to 10 seconds post-collapse), even though these layers have become stable to convection. We also find the equilibrium shape of a thin magnetic flux rope in the dense hydrostatic atmosphere of the neutron star, along with the critical separation of the footpoints above which the rope undergoes unlimited expansion against gravity. The implications of these results for pulsar magnetism are summarized, and applied to the case of late fallback over the first 1,000-10,000 s of the life of a neutron star

We examine the eccentricity evolution of a system of two planets locked in a mean motion resonance, in which the outer planet loses energy and angular momentum. The sink of energy and angular momentum could be either a gas or planetesimal disk. We show that the eccentricity of both planetary bodies can grow to large values, particularly if the inner body does not directly exchange energy or angular momentum with the disk. We analytically calculate the eccentricity damping rate in the case of a single planet migrating through a planetesimal disk. We present the results of numerical integrations of two resonant planets showing rapid growth of eccentricity. We also present integrations in which a Jupiter-mass planet is forced to migrate inward through a system of 5-10 roughly Earth mass planets. The migrating planet can eject or accrete the smaller bodies; roughly 5% of the mass (averaged over all the integrations) accretes onto the central star. The results are discussed in the context of the currently known extrasolar planetary systems.

We study spectroscopically determined iron abundances of 642 solar-type stars to search for the signature of accreted iron-rich material. We find that the metallicity [Fe/H] of a subset of 466 main sequence stars, when plotted as a function of stellar mass, mimics the pattern seen in lithium abundances in open clusters. Using Monte Carlo models we find that, on average, these stars have accreted about 0.4 Earth masses of iron while on the main sequence. A much smaller sample of 19 stars in the Hertzsprung gap, which are slightly evolved and whose convection zones are significantly more massive, have lower average [Fe/H], and their metallicity shows no clear variation with stellar mass. These findings suggest that terrestrial-type material is common around solar type stars.

We present simultaneous observations by EUVE, ASCA, and RXTE of the type~1 Seyfert galaxy NGC 5548. These data indicate that variations in the EUV emission (at $\sim 0.2 $keV) appear to lead similar modulations in higher energy ($\ga 1 $keV) X-rays by $\sim$10--30 ks. This is contrary to popular models which attribute the correlated variability of the EUV, UV and optical emission in type~1 Seyferts to reprocessing of higher energy radiation. This behavior instead suggests that the variability of the optical through EUV emission is an important driver for the variability of the harder X-rays which are likely produced by thermal Comptonization. We also investigate the spectral characteristics of the fluorescent iron K$\alpha$ line and Compton reflection emission. In contrast to prior measurements of these spectral features, we find that the iron K$\alpha$ line has a relatively small equivalent width ($W_{K\alpha} \sim 100 $eV) and that the reflection component is consistent with a covering factor which is significantly less than unity ($\Omega/2\pi \sim 0.4$--0.5). Notably, although the 2--10 keV X-ray flux varies by $\sim \pm 25$% and the derived reflection fraction appears to be constant throughout our observations, the flux in the Fe~K$\alpha$ line is also constant. This behavior is difficult to reconcile in the context of standard Compton reflection models.

Classical analytic theories of the solar system indicate that it is stable, but numerical integrations suggest that it is chaotic. This disagreement is resolved by a new analytic theory. The theory shows that the chaos among the Jovian planets results from the overlap of the components of a mean motion resonance among Jupiter, Saturn, and Uranus, and provides rough estimates of the Lyapunov time (10 million years) and the dynamical lifetime of Uranus (10^{18} years). The Jovian planets must have entered the resonance after all the gas and most of the planetesimals in the protoplanetary disk were removed.

A planet orbiting in a disk of planetesimals can experience an instability in which it migrates to smaller orbital radii. Resonant interactions between the planet and planetesimals remove angular momentum from the planetesimals, increasing their eccentricities. Subsequently, the planetesimals either collide with or are ejected by the planet, reducing the semimajor axis of the planet. If the surface density of planetesimals exceeds a critical value, corresponding to 0.03 solar masses of gas inside the orbit of Jupiter, the planet will migrate inward a large distance. This instability may explain the presence of Jupiter-mass objects in small orbits around nearby stars.

Using the disk wind model of Murray et al. (1995), we calculate line profiles and frequency-resolved response functions for broad line emission from the surface of an accretion disk in an AGN in the presence of a radiatively driven wind. We find that the combined effects of the shears in the wind and in the disk itself produce anisotropic line emission which solves several well-known problems connected with disk models of the broad line region. In particular, the broadening of resonance lines such as \Civ, \Lya, and \Nv\/ can be attributed to orbital motion of the disk gas at radii as close as $\sim 10^{16}$~cm in Seyferts without requiring unrealistically large emission regions in order to produce single-peaked profiles. Furthermore, the anisotropy of the line emission results in frequency-dependent response functions which are no longer red-blue symmetric so that the time delays inferred for the various red and blue components of the line agree qualitatively with recent reverberation mapping observations of NGC~5548.