We derive the UV-optical stellar dust attenuation curve of galaxies at z=1.4-2.6 as a function of gas-phase metallicity. We use a sample of 218 star-forming galaxies, excluding those with very young or heavily obscured star formation, from the MOSFIRE Deep Evolution Field (MOSDEF) survey with H$\alpha$, H$\beta$, and [NII]$\lambda 6585$ spectroscopic measurements. We constrain the shape of the attenuation curve by comparing the average flux densities of galaxies sorted into bins of dust obscuration using Balmer decrements, i.e., H$\alpha$-to-H$\beta$ luminosities. The average attenuation curve for the high-metallicity sample (12+log(O/H)>8.5, corresponding to $M_*\gtrsim10^{10.4}\,M_{\odot}$) has a shallow slope, identical to that of the Calzetti local starburst curve, and a significant UV 2175A extinction bump that is $\sim 0.5\times$ the strength of the Milky Way bump. On the other hand, the average attenuation curve of the low-metallicity sample (12+log(O/H) $\sim 8.2-8.5$) has a steeper slope similar to that of the SMC curve, only consistent with the Calzetti slope at the $3\sigma$ level. The UV bump is not detected in the low-metallicity curve, indicating the relative lack of the small dust grains causing the bump at low metallicities. Furthermore, we find that on average the nebular reddening (E(B-V)) is a factor of 2 times larger than that of the stellar continuum for galaxies with low metallicities, while the nebular and stellar reddening are similar for galaxies with higher metallicities. The latter is likely due to a high surface density of dusty clouds embedding the star forming regions but also reddening the continuum in the high-metallicity galaxies.

We present ionized gas kinematics for 681 galaxies at $z\sim 1.4-3.8$ from the MOSFIRE Deep Evolution Field survey, measured using models which account for random galaxy-slit misalignments together with structural parameters derived from CANDELS Hubble Space Telescope (HST) imaging. Kinematics and sizes are used to derive dynamical masses. Baryonic masses are estimated from stellar masses and inferred gas masses from dust-corrected star formation rates (SFRs) and the Kennicutt-Schmidt relation. We measure resolved rotation for 105 galaxies. For the remaining 576 galaxies we use models based on HST imaging structural parameters together with integrated velocity dispersions and baryonic masses to statistically constrain the median ratio of intrinsic ordered to disordered motion, $V/\sigma_{V,0}$. We find that $V/\sigma_{V,0}$ increases with increasing stellar mass and decreasing specific SFR (sSFR). These trends may reflect marginal disk stability, where systems with higher gas fractions have thicker disks. For galaxies with detected rotation we assess trends between their kinematics and mass, sSFR, and baryon surface density ($\Sigma_{\mathrm{bar},e}$). Intrinsic dispersion correlates most with $\Sigma_{\mathrm{bar},e}$ and velocity correlates most with mass. By comparing dynamical and baryonic masses, we find that galaxies at $z\sim 1.4-3.8$ are baryon dominated within their effective radii ($R_E$), with Mdyn/Mbaryon increasing over time. The inferred baryon fractions within $R_E$, $f_{\mathrm{bar}}$, decrease over time, even at fixed mass, size, or surface density. At fixed redshift, $f_{\mathrm{bar}}$ does not appear to vary with stellar mass but increases with decreasing $R_E$ and increasing $\Sigma_{\mathrm{bar},e}$. For galaxies at $z\geq2$, the median inferred baryon fractions generally exceed 100%. We discuss possible explanations and future avenues to resolve this tension.

We present measurements of [SIII]$\lambda\lambda$9069,9531 for a sample of $z\sim1.5$ star-forming galaxies, the first sample with measurements of these lines at z>0.1. We employ the line ratio S$_{32}$$\equiv$[SIII]$\lambda\lambda$9069,9531/[SII]$\lambda\lambda$6716,6731 as a novel probe of evolving ISM conditions. Since this ratio includes the low-ionization line [SII], it is crucial that the effects of diffuse ionized gas (DIG) on emission-line ratios be accounted for in $z\sim0$ integrated galaxy spectra, or else that comparisons be made to samples of local HII regions in which DIG emission is not present. We find that S$_{32}$ decreases with increasing stellar mass at both $z\sim1.5$ and $z\sim0$, but that the dependence is weak suggesting S$_{32}$ has a very shallow anticorrelation with metallicity, in contrast with O$_{32}$ that displays a strong metallicity dependence. As a result, S$_{32}$ only mildly evolves with redshift at fixed stellar mass. The $z\sim1.5$ sample is systematicallty offset towards lower S$_{32}$ and higher [SII]/H$\alpha$ at fixed [OIII]/H$\beta$ relative to $z=0$ HII regions. By comparing to photoionization model grids, we find that such trends can be explained by a scenario in which the ionizing spectrum is harder at fixed O/H with increasing redshift, but are inconsistent with an increase in ionization parameter at fixed O/H. This analysis demonstrates the advantages of expanding beyond the strongest rest-optical lines for evolutionary studies, and the particular utility of [SIII] for characterizing evolving ISM conditions and stellar compositions. These measurements provide a basis for estimating [SIII] line strengths for high-redshift galaxies, a line that the James Webb Space Telescope will measure out to z~5.5.

The five globular clusters (GCs) of the Fornax dSph are puzzling for two reasons; the mass in GCs is high with respect to the galaxy's old stellar mass, and their survival and large distance (> 1 kpc) is at odds with naive expectations of dynamical friction. We present here a semi-analytic model, simultaneously addressing both problems in a comprehensive evolutionary framework for Fornax. Key to the model is inclusion of: 1) hydrodynamical constraints on the GC formation locations, 2) self-consistent velocity distribution functions in the dynamical friction calculations and 3) expansion of GC orbits due to a past dwarf-dwarf merger in the orbit integrations. The latter is crucial for reconciling the dynamical survival of the clusters, and their chemical properties with respect to the Fornax field stars. We find that in order for four of the GCs to survive at their observed projected location, a dark matter core of radius > 1.5 kpc and a dwarf merger with dynamical mass ratio of 1:5 to 1:2 with Fornax is required. We support the merger scenario by showing that aspects of the field star metallicity distribution function and anomalous chemical properties of GC5, are representative of a merging galaxy which is ~1/3 less massive than Fornax. Together the chemical and dynamical models suggest a scenario where three in-situ GCs in proto-Fornax were ejected to the outskirts during the merger, a GC4 formed during the merger at about 10 Gyrs ago, with GC5 being brought in by the merging galaxy to Fornax.

Ninety per cent of baryons are located outside galaxies, either in the circumgalactic or intergalactic medium. Theory points to galactic winds as the primary source of the enriched and massive circumgalactic medium. Winds from compact starbursts have been observed to flow to distances somewhat greater than ten kiloparsecs, but the circumgalactic medium typically extends beyond a hundred kiloparsecs. Here we report optical integral field observations of the massive but compact galaxy SDSS J211824.06+001729.4. The oxygen [O II] lines at wavelengths of 3726 and 3729 angstroms reveal an ionized outflow spanning 80 by 100 square kiloparsecs, depositing metal-enriched gas at 10,000 kelvin through an hourglass-shaped nebula that resembles an evacuated and limb-brightened bipolar bubble. We also observe neutral gas phases at temperatures of less than 10,000 kelvin reaching distances of 20 kiloparsecs and velocities of around 1,500 kilometres per second. This multi-phase outflow is probably driven by bursts of star formation, consistent with theory.

We present detections of [OIII]$\lambda$4363 and direct-method metallicities for star-forming galaxies at $z=1.7-3.6$. We combine new measurements from the MOSFIRE Deep Evolution Field (MOSDEF) survey with literature sources to construct a sample of 18 galaxies with direct-method metallicities at $z>1$, spanning $7.5<1$2+log(O/H$)<8.2$ and log(M$_*$/M$_{\odot})=7-10$. We find that strong-line calibrations based on local analogs of high-redshift galaxies reliably reproduce the metallicity of the $z>1$ sample on average. We construct the first mass-metallicity relation at $z>1$ based purely on direct-method O/H, finding a slope that is consistent with strong-line results. Direct-method O/H evolves by $\lesssim0.1$ dex at fixed M$_*$ and SFR from $z\sim0-2.2$. We employ photoionization models to constrain the ionization parameter and ionizing spectrum in the high-redshift sample. Stellar models with super-solar O/Fe and binary evolution of massive stars are required to reproduce the observed strong-line ratios. We find that the $z>1$ sample falls on the $z\sim0$ relation between ionization parameter and O/H, suggesting no evolution of this relation from $z\sim0$ to $z\sim2$. These results suggest that the offset of the strong-line ratios of this sample from local excitation sequences is driven primarily by a harder ionizing spectrum at fixed nebular metallicity compared to what is typical at $z\sim0$, naturally explained by super-solar O/Fe at high redshift caused by rapid formation timescales. Given the extreme nature of our $z>1$ sample, the implications for representative $z\sim2$ galaxy samples at $\sim10^{10}$ M$_{\odot}$ are unclear, but similarities to $z>6$ galaxies suggest that these conclusions can be extended to galaxies in the epoch of reionization.

Using data from the MOSFIRE Deep Evolution Field (MOSDEF) survey, we present a census of AGN-driven ionized outflows in a sample of 159 AGNs at $1.4 \le z \le 3.8$. The sample spans AGN bolometric luminosities of $10^{44-47} \mathrm{~erg~s}^{-1}$ and includes both quiescent and star-forming galaxies extending across three orders of magnitude in stellar mass. We identify and characterize outflows from the \hbeta, [OIII], \halpha ~and [NII] emission line spectra. We detect outflows in $17\%$ of the AGNs, seven times more often than in a mass-matched sample of inactive galaxies in MOSDEF. The outflows are fast and galaxy-wide, with velocities of $\sim 400-3500 ~\mathrm{km~s}^{-1}$ and spatial extents of $0.3-11.0$ kpc. The incidence of outflows among AGNs is independent of the stellar mass of the host galaxy, with outflows detected in both star-forming and quiescent galaxies. This suggests that outflows exist across different phases in galaxy evolution. We investigate relations between outflow kinematic, spatial, and energetic properties and both AGN and host galaxy properties. Our results show that AGN-driven outflows are widespread in galaxies along the star-forming main sequence. The mass-loading factors of the outflows are typically $0.1-1$ and increase with AGN luminosity, capable of exceeding unity at $L_\mathrm{AGN} \gtrsim 10^{46.3} \mathrm{~erg~s}^{-1}$. In these more luminous sources the ionized outflow alone is likely sufficient to regulate star formation, and when combined with outflowing neutral and molecular gas may be able to quench star formation in their host galaxies.

Population synthesis models predict that high-mass X-ray binary (HMXB) populations produced in low metallicity environments should be more X-ray luminous, a trend supported by studies of nearby galaxies. This trend may be responsible for the observed increase of the X-ray luminosity ($L_{\mathrm{X}}$) per star formation rate (SFR) with redshift due to the decrease of metallicity ($Z$) at fixed stellar mass as a function of redshift. To test this hypothesis, we use a sample of 79 $z\sim2$ star-forming galaxies with oxygen abundance measurements from the MOSDEF survey, which obtained rest-frame optical spectra for $\sim1500$ galaxies in the CANDELS fields at $1.37<z<3.80$. Using Chandra data from the AEGIS-X Deep, Deep Field North, and Deep Field South surveys, we stack the X-ray data at the galaxy locations in bins of redshift and $Z$ because the galaxies are too faint to be individually detected. In agreement with previous studies, the average $L_{\mathrm{X}}$/SFR of our $z\sim2$ galaxy sample is enhanced by $\approx0.4-0.8$ dex relative to local HMXB $L_{\mathrm{X}}$-SFR scaling relations. Splitting our sample by $Z$, we find that $L_{\mathrm{X}}$/SFR and $Z$ are anti-correlated with 97% confidence. This observed $Z$ dependence for HMXB-dominated galaxies is consistent both with the local $L_{\mathrm{X}}$-SFR-$Z$ relation and a subset of population synthesis models. Although the statistical significance of the observed trends is weak due to the low X-ray statistics, these results constitute the first direct evidence connecting the redshift evolution of $L_{\mathrm{X}}$/SFR and the $Z$ dependence of HMXBs.

We present results on the emission-line properties of 1.3<=z<=2.7 galaxies drawn from the complete MOSFIRE Deep Evolution Field (MOSDEF) survey. Specifically, we use observations of the emission-line diagnostic diagram of [OIII]5007/Hb vs. [SII]6717,6731/Ha, i.e., the "[SII] BPT diagram," to gain insight into the physical properties of high-redshift star-forming regions. High-redshift MOSDEF galaxies are offset towards lower [SII]6717,6731/Ha at fixed [OIII]5007/Hb, relative to local galaxies from the Sloan Digital Sky Survey (SDSS). Furthermore, at fixed [OIII]5007/Hb, local SDSS galaxies follow a trend of decreasing [SII]6717,6731/Ha as the surface density of star formation (Sigma_SFR) increases. We explain this trend in terms of the decreasing fractional contribution from diffuse ionized gas (f_DIG) as Sigma_SFR increases in galaxies, which causes galaxy-integrated line ratios to shift towards the locus of pure HII-region emission. The z~0 relationship between f_DIG and Sigma_SFR implies that high-redshift galaxies have lower f_DIG values than typical local systems, given their significantly higher typical Sigma_SFR. When an appropriate low-redshift benchmark with zero or minimal f_DIG is used, high-redshift MOSDEF galaxies appear offset towards higher [SII]6717,6731/Ha and/or [OIII]5007/Hb. The joint shifts of high-redshift galaxies in the [SII] and [NII] BPT diagrams are best explained in terms of the harder spectra ionizing their star-forming regions at fixed nebular oxygen abundance (expected for chemically-young galaxies), as opposed to large variations in N/O ratios or higher ionization parameters. The evolving mixture of HII regions and DIG is an essential ingredient to our description of the ISM over cosmic time.

Stellar feedback plays a significant role in modulating star formation, redistributing metals, and shaping the baryonic and dark structure of galaxies -- however, the efficiency of its energy deposition to the interstellar medium is challenging to constrain observationally. Here we leverage HST and ALMA imaging of a molecular gas and dust shell ($M_{H2} \sim 2\times 10^{5} ~{\rm M}_{\odot}$) in an outflow from the nuclear star forming ring of the galaxy NGC 3351, to serve as a boundary condition for a dynamical and energetic analysis of the outflowing ionised gas seen in our MUSE TIMER survey. We use \texttt{STARBURST99} models and prescriptions for feedback from simulations to demonstrate that the observed star formation energetics can reproduce the ionised and molecular gas dynamics -- provided a dominant component of the momentum injection comes from direct photon pressure from young stars, on top of supernovae, photoionisation heating and stellar winds. The mechanical energy budget from these sources is comparable to low luminosity AGN, suggesting that stellar feedback can be a relevant driver of bulk gas motions in galaxy centres - although here $\lesssim 10^{-3}$ of the ionized gas mass is escaping the galaxy. We test several scenarios for the survival/formation of the cold gas in the outflow, including in-situ condensation and cooling. Interestingly, the geometry of the molecular gas shell, observed magnetic field strengths and emission line diagnostics are consistent with a scenario where magnetic field lines aided survival of the dusty ISM as it was initially launched (with mass loading factor $\lesssim 1$) from the ring by stellar feedback. This system's unique feedback driven morphology can hopefully serve as a useful litmus test for feedback prescriptions in magnetohydrodynamical galaxy simulations.

We study the properties of 30 spectroscopically-identified pairs of galaxies observed during the peak epoch of star formation in the universe. These systems are drawn from the MOSFIRE Deep Evolution Field (MOSDEF) Survey at $1.4 \leq z \leq 3.8$, and are interpreted as early-stage galaxy mergers. Galaxy pairs in our sample are identified as two objects whose spectra were collected on the same Keck/MOSFIRE spectroscopic slit. Accordingly, all pairs in the sample have projected separations $R_{{\rm proj}}\leq 60$ kpc. The velocity separation for pairs was required to be $\Delta v \leq 500 \mbox{ km s}^{-1}$, which is a standard threshold for defining interacting galaxy pairs at low redshift. Stellar mass ratios in our sample range from 1.1 to 550, with 12 ratios closer than or equal to 3:1, the common definition of a "major merger." Studies of merging pairs in the local universe indicate anenhancement in star-formation activity and deficit in gas-phase oxygen abundance relative to isolated galaxies of the same mass. We compare the MOSDEF pairs sample to a control sample of isolated galaxies at the same redshift, finding no measurable SFR enhancement or metallicity deficit at fixed stellar mass for the pairs sample. The lack of significant difference between the average properties of pairs and control samples appears in contrast to results from low-redshift studies, although the small sample size and lower signal-to-noise of the high-redshift data limit definitive conclusions on redshift evolution. These results are consistent with some theoretical works suggesting a reduced differential effect of pre-coalescence mergers on galaxy properties at high redshift -- specifically that pre-coalescence mergers do not drive strong starbursts.

We use extensive spectroscopy from the MOSFIRE Deep Evolution Field (MOSDEF) survey to investigate the relationships between rest-frame optical emission line equivalent widths ($W$) and a number of galaxy and ISM characteristics for a sample of $1134$ star-forming galaxies at redshifts $1.4\lesssim z\lesssim 3.8$. We examine how the equivalent widths of [OII]$\lambda\lambda 3727, 3730$, H$\beta$, [OIII]$\lambda\lambda 4960, 5008$, [OIII]$+$H$\beta$, H$\alpha$, and H$\alpha$+[NII]$\lambda\lambda 6550, 6585$, depend on stellar mass, UV slope, age, star-formation rate (SFR) and specific SFR (sSFR), ionization parameter and excitation conditions (O32 and [OIII]/H$\beta$), gas-phase metallicity, and ionizing photon production efficiency ($\xi_{\rm ion}$). The trend of increasing $W$ with decreasing stellar mass is strongest for [OIII] (and [OIII]+H$\beta$). More generally, the equivalent widths of all the lines increase with redshift at a fixed stellar mass or fixed gas-phase metallicity, suggesting that high equivalent width galaxies are common at high redshift. This redshift evolution in equivalent widths can be explained by the increase in SFR and decrease in metallicity with redshift at a fixed stellar mass. Consequently, the dependence of $W$ on sSFR is largely invariant with redshift, particularly when examined for galaxies of a given metallicity. Our results show that high equivalent width galaxies, specifically those with high $W({\rm [OIII]})$, have low stellar masses, blue UV slopes, young ages, high sSFRs, ISM line ratios indicative of high ionization parameters, high $\xi_{\rm ion}$, and low metallicities. As these characteristics are often attributed to galaxies with high ionizing escape fractions, galaxies with high $W$ are likely candidates for the population that dominates cosmic reionization.

Using the MOSDEF rest-frame optical spectroscopic survey, we investigate the star-formation histories (SFHs) of different galaxy types, ranging from actively star forming to quiescent at $1.4\leq~z\leq2.6$. SFHs are constrained utilizing stellar continuum spectroscopy, specifically through a combination of Balmer absorption lines, the 4000$\AA$ break, and the equivalent width of the H$\alpha$ emission line. To attain a sufficiently high signal-to-noise ratio (S/N) to conduct these measurements we stack spectra of galaxies with similar spectral types, as determined from their rest-frame $U-V$ and $V-J$ colors. We bin the MOSDEF sample into five spectral types, subdividing the quiescent and star-forming bins to better explore galaxies transitioning between the two. We constrain the average SFHs for each type, finding that quiescent and transitional galaxies in the MOSDEF sample are dominated by an SFH with an average star-formation timescale of $\tau\sim0.1-0.2$~Gyr. These findings contrast with measurements from the low-redshift Universe where, on average, galaxies form their stars over a more extended time period ($\tau>1$~Gyr). Furthermore, our spectral index measurements correlate with mass surface density for all spectral types. Finally, we compare the average properties of the galaxies in our transitional bins to investigate possible paths to quiescence, and speculate on the viability of a dusty post-starburst phase.

The Time Inference with MUSE in Extragalactic Rings (TIMER) project is a survey with the VLT-MUSE integral-field spectrograph of 24 nearby barred galaxies with prominent central structures (e.g., nuclear rings or inner discs). The main goals of the project are: (i) estimating the cosmic epoch when discs of galaxies settle, leading to the formation of bars; (ii) testing the hypothesis whereby discs in more massive galaxies are assembled first; and (iii) characterising the history of external gas accretion in disc galaxies. We present details on the sample selection, observations, data reduction, and derivation of high-level data products, including stellar kinematics, ages and metallicities. We also derive star formation histories and physical properties and kinematics of ionised gas. We illustrate how this dataset can be used for a plethora of scientific applications, e.g., stellar feedback, outflows, nuclear and primary bars, stellar migration and chemical enrichment, and the gaseous and stellar dynamics of nuclear spiral arms, barlenses, box/peanuts and bulges. Amongst our first results - based on a few selected galaxies -, we show that the dynamics of nuclear rings and inner discs is consistent with the picture in which they are formed by bars, that the central few hundred parsecs in massive disc galaxies tend to show a pronounced peak in stellar metallicity, and that nuclear rings can efficiently prevent star formation in this region. Finally, we present evidence that star-bursting nuclear rings can be fed with low-metallicity gas from low-mass companions.

We investigate the nature of the relation among stellar mass, star-formation rate, and gas-phase metallicity (the M$_*$-SFR-Z relation) at high redshifts using a sample of 260 star-forming galaxies at $z\sim2.3$ from the MOSDEF survey. We present an analysis of the high-redshift M$_*$-SFR-Z relation based on several emission-line ratios for the first time. We show that a M$_*$-SFR-Z relation clearly exists at $z\sim2.3$. The strength of this relation is similar to predictions from cosmological hydrodynamical simulations. By performing a direct comparison of stacks of $z\sim0$ and $z\sim2.3$ galaxies, we find that $z\sim2.3$ galaxies have $\sim0.1$ dex lower metallicity at fixed M$_*$ and SFR. In the context of chemical evolution models, this evolution of the M$_*$-SFR-Z relation suggests an increase with redshift of the mass-loading factor at fixed M$_*$, as well as a decrease in the metallicity of infalling gas that is likely due to a lower importance of gas recycling relative to accretion from the intergalactic medium at high redshifts. Performing this analysis simultaneously with multiple metallicity-sensitive line ratios allows us to rule out the evolution in physical conditions (e.g., N/O ratio, ionization parameter, and hardness of the ionizing spectrum) at fixed metallicity as the source of the observed trends with redshift and with SFR at fixed M$_*$ at $z\sim2.3$. While this study highlights the promise of performing high-order tests of chemical evolution models at high redshifts, detailed quantitative comparisons ultimately await a full understanding of the evolution of metallicity calibrations with redshift.

We present an analysis using the MOSFIRE Deep Evolution Field (MOSDEF) survey on the nature of "MIR-excess" galaxies, which have star formation rates (SFR) inferred from mid-infrared (MIR) data that is substantially elevated relative to that estimated from dust-corrected UV data. We use a sample of $\sim$200 galaxies and AGN at $1.40<z<2.61$ with 24 $\mu$m detections (rest-frame 8$\mu$m) from MIPS/\textit{Spitzer}. We find that the identification of MIR-excess galaxies strongly depends on the methodologies used to estimate IR luminosity ($\rm L_{IR}$) and to correct the UV light for dust attenuation. We find that extrapolations of the SFR from the observed 24 $\mu$m flux, using luminosity-dependent templates based on local galaxies, substantially overestimate $\rm L_{IR}$ in $z\sim2$ galaxies. By including \textit{Herschel} observations and using a stellar mass-dependent, luminosity-independent $\rm L_{IR}$, we obtain more reliable estimates of the SFR and a lower fraction of MIR-excess galaxies. Once stellar mass selection biases are taken into account, we identify $\sim24\%$ of our galaxies as MIR-excess. However, $\rm SFR_{H\alpha}$ is not elevated in MIR-excess galaxies compared to MIR-normal galaxies, indicating that the intrinsic fraction of MIR-excess may be lower. Using X-ray, IR, and optically-selected AGN in MOSDEF, we do not find a higher prevalence for AGN in MIR-excess galaxies relative to MIR-normal galaxies. A stacking analysis of X-ray undetected galaxies does not reveal a harder spectrum in MIR-excess galaxies relative to MIR-normal galaxies. Our analysis indicates that AGN activity does not contribute substantially to the MIR excess and instead implies that it is likely due to the enhanced PAH emission.

We present a comparative study of molecular and ionized gas kinematics in nearby galaxies. These results are based on observations from the EDGE survey, which measured spatially resolved $^{12}$CO(J=1-0) in 126 nearby galaxies. Every galaxy in EDGE has corresponding resolved ionized gas measurements from CALIFA. Using a sub-sample of 17 rotation dominated, star-forming galaxies where precise molecular gas rotation curves could be extracted, we derive CO and H$\alpha$ rotation curves using the same geometric parameters out to $\gtrsim$1 $R_e$. We find that $\sim$75% of our sample galaxies have smaller ionized gas rotation velocities than the molecular gas in the outer part of the rotation curve. In no case is the molecular gas rotation velocity measurably lower than that of the ionized gas. We suggest that the lower ionized gas rotation velocity can be attributed to a significant contribution from extraplanar diffuse ionized gas in a thick, turbulence supported disk. Using observations of the H$\gamma$ transition also available from CALIFA, we measure ionized gas velocity dispersions and find that these galaxies have sufficiently large velocity dispersions to support a thick ionized gas disk. Kinematic simulations show that a thick disk with a vertical rotation velocity gradient can reproduce the observed differences between the CO and H$\alpha$ rotation velocities. Observed line ratios tracing diffuse ionized gas are elevated compared to typical values in the midplane of the Milky Way. In galaxies affected by this phenomenon, dynamical masses measured using ionized gas rotation curves will be systematically underestimated.

Deriving circular velocities of galaxies from stellar kinematics can provide an estimate of their total dynamical mass, provided a contribution from the velocity dispersion of the stars is taken into account. Molecular gas (e.g., CO) on the other hand, is a dynamically cold tracer and hence acts as an independent circular velocity estimate without needing such a correction. In this paper we test the underlying assumptions of three commonly used dynamical models, deriving circular velocities from stellar kinematics of 54 galaxies (S0-Sd) that have observations of both stellar kinematics from the CALIFA survey, and CO kinematics from the EDGE survey. We test the Asymmetric Drift Correction (ADC) method, as well as Jeans, and Schwarzschild models. The three methods each reproduce the CO circular velocity at 1Re to within 10%. All three methods show larger scatter (up to 20%) in the inner regions (R < 0.4Re) which may be due to an increasingly spherical mass distribution (which is not captured by the thin disk assumption in ADC), or non-constant stellar M/L ratios (for both the JAM and Schwarzschild models). This homogeneous analysis of stellar and gaseous kinematics validates that all three models can recover Mdyn at 1Re to better than 20%, but users should be mindful of scatter in the inner regions where some assumptions may break down.

We combine spectroscopic measurements of H$\alpha$ and H$\beta$ and UV continuum photometry for a sample of 673 galaxies from the MOSFIRE Deep Evolution Field survey to constrain hydrogen ionizing photon production efficiencies ($\xi_{\rm ion}$, xi_ion) at z=1.4-2.6. We find average log(xi_ion/[Hz erg$^{-1}$])=25.06 (25.34), assuming the Calzetti (SMC) curve for the UV dust correction and a scatter of 0.28 dex in xi_ion distribution. After accounting for observational uncertainties and variations in dust attenuation, we conclude that the remaining scatter in xi_ion is likely dominated by galaxy-to-galaxy variations in stellar populations, including the slope and upper-mass cutoff of the initial mass function, stellar metallicity, star-formation burstiness, and stellar evolution (e.g., single/binary star evolution). Moreover, xi_ion is elevated in galaxies with high ionization states (high [OIII]/[OII]) and low oxygen abundances (low [NII]/H$\alpha$ and high [OIII]/H$\beta$) in the ionized ISM. However, xi_ion does not correlate with the offset from the z~0 star-forming locus in the BPT diagram, suggesting no change in the hardness of ionizing radiation accompanying the offset from the z~0 sequence. We also find that galaxies with blue UV spectral slopes ($\langle\beta\rangle$=-2.1) have elevated xi_ion by a factor of ~2 relative to the average xi_ion of the sample ($\langle\beta\rangle$=-1.4). If these blue galaxies are similar to those at z > 6, our results suggest that a lower Lyman continuum escape fraction is required for galaxies to maintain reionization, compared to the canonical xi_ion predictions from stellar population models. Furthermore, we demonstrate that even with robustly dust-corrected H$\alpha$, the UV dust attenuation can cause on average a ~0.3dex systematic uncertainty in xi_ion calculations.

After changing optical AGN type from 1.9 to 1 in 1984, the AGN Mrk 1018 recently reverted back to its type 1.9 state. Our ongoing monitoring now reveals that the AGN has halted its dramatic dimming, reaching a minimum around October 2016. The minimum was followed by an outburst rising with $\sim$0.25 U-band mag/month. The rebrightening lasted at least until February 2017, as confirmed by joint Chandra and Hubble observations. Monitoring was resumed in July 2017 after the source emerged from sunblock, at which point the AGN was found only $\sim$0.4 mag brighter than its minimum. The intermittent outburst was accompanied by the appearance of a red wing asymmetry in broad-line shape, indicative of an inhomogeneous broad-line region. The current flickering brightness of Mrk 1018 following its rapid fading either suggests that the source has reignited, remains variable at a low level, or may continue dimming over the next few years. Distinguishing between these possibilities requires continuous multiwavelength monitoring.

We present results from the MOSFIRE Deep Evolution Field (MOSDEF) survey on broad flux from the nebular emission lines H$\alpha$, [NII], [OIII], H$\beta$, and [SII]. The sample consists of 127 star-forming galaxies at $1.37 < z < 2.61$ and 84 galaxies at $2.95 < z < 3.80$. We decompose the emission lines using narrow ($\text{FWHM} < 275 \ \text{km s}^{-1}$) and broad ($\text{FWHM} > 300 \ \text{km s}^{-1}$) Gaussian components for individual galaxies and stacks. Broad emission is detected at $>3\sigma$ in $<10$% of galaxies and the broad flux accounts for 10-70% of the total flux. We find a slight increase in broad to narrow flux ratio with mass but note that we cannot reliably detect broad emission with $\text{FWHM} < 275 \ \text{km s}^{-1}$, which may be significant at low masses. Notably, there is a correlation between higher signal-to-noise (S/N) spectra and a broad component detection indicating a S/N dependence in our ability to detect broad flux. When placed on the N2-BPT diagram ([OIII]/H$\beta$ vs. [NII]/H$\alpha$) the broad components of the stacks are shifted towards higher [OIII]/H$\beta$ and [NII]/$\alpha$ ratios compared to the narrow component. We compare the location of the broad components to shock models and find that the broad component could be due to shocks, but we do not rule out other possibilities such as the presence of an AGN. We estimate the mass loading factor (mass outflow rate/star formation rate) assuming the broad component is a photoionized outflow and find that the mass loading factor increases as a function of mass which agrees with previous studies. We show that adding emission from shocked gas to $z\sim0$ SDSS spectra shifts galaxies towards the location of $z\sim2$ galaxies on several emission line diagnostic diagrams.

Using observations from the first two years of the MOSFIRE Deep Evolution Field (MOSDEF) survey, we study 13 AGN-driven outflows detected from a sample of 67 X-ray, IR and/or optically-selected AGN at $z \sim 2$. The AGN have bolometric luminosities of $\sim10^{44}-10^{46} ~\mathrm{erg~s^{-1}}$, including both quasars and moderate-luminosity AGN. We detect blueshifted, ionized gas outflows in the H$\beta$ , [OIII], H$\alpha$ ~and/or [NII] emission lines of $19\%$ of the AGN, while only 1.8\% of the MOSDEF galaxies have similarly-detected outflows. The outflow velocities span $\sim$300 to 1000 km s$^{-1}$. Eight of the 13 outflows are spatially extended on similar scales as the host galaxies, with spatial extents of 2.5 to 11.0 kpc. Outflows are detected uniformly across the star-forming main sequence, showing little trend with the host galaxy SFR. Line ratio diagnostics indicate that the outflowing gas is photoionized by the AGN. We do not find evidence for positive AGN feedback, in either our small MOSDEF sample or a much larger SDSS sample, using the BPT diagram. Given that a galaxy with an AGN is ten times more likely to have a detected outflow, the outflowing gas is photoionzed by the AGN, and estimates of the mass and energy outflow rates indicate that stellar feedback is insufficient to drive at least some of these outflows, they are very likely to be AGN-driven. The outflows have mass-loading factors of the order of unity, suggesting that they help regulate star formation in their host galaxies, though they may be insufficient to fully quench it.

We present the first spectroscopic measurement of multiple rest-frame optical emission lines at $z>4$. During the MOSFIRE Deep Evolution Field (MOSDEF) survey, we observed the galaxy GOODSN-17940 with the Keck I/MOSFIRE spectrograph. The K-band spectrum of GOODSN-17940 includes significant detections of the [OII]$\lambda\lambda 3726,3729$, [NeIII]$\lambda3869$, and H$\gamma$ emission lines and a tentative detection of H$\delta$, indicating $z_{\rm{spec}}=4.4121$. GOODSN-17940 is an actively star-forming $z>4$ galaxy based on its K-band spectrum and broadband spectral energy distribution. A significant excess relative to the surrounding continuum is present in the Spitzer/IRAC channel 1 photometry of GOODSN-17940, due primarily to strong H$\alpha$ emission with a rest-frame equivalent width of $\mbox{EW(H}\alpha)=1200$ \AA. Based on the assumption of $0.5 Z_{\odot}$ models and the Calzetti attenuation curve, GOODSN-17940 is characterized by $M_*=5.0^{+4.3}_{-0.2}\times 10^9 M_{\odot}$. The Balmer decrement inferred from H$\alpha$/H$\gamma$ is used to dust correct the H$\alpha$ emission, yielding $\mbox{SFR(H}\alpha)=320^{+190}_{-140} M_{\odot}\mbox{ yr}^{-1}$. These $M_*$ and SFR values place GOODSN-17940 an order of magnitude in SFR above the $z\sim 4$ star-forming "main sequence." Finally, we use the observed ratio of [NeIII]/[OII] to estimate the nebular oxygen abundance in GOODSN-17940, finding $\mbox{O/H}\sim 0.2 \mbox{ (O/H)}_{\odot}$. Combining our new [NeIII]/[OII] measurement with those from stacked spectra at $z\sim 0, 2, \mbox{ and } 3$, we show that GOODSN-17940 represents an extension to $z>4$ of the evolution towards higher [NeIII]/[OII] (i.e., lower $\mbox{O/H}$) at fixed stellar mass. It will be possible to perform the measurements presented here out to $z\sim 10$ using the James Webb Space Telescope.

We present interferometric CO observations made with the Combined Array for Millimeter-wave Astronomy (CARMA) of galaxies from the Extragalactic Database for Galaxy Evolution survey (EDGE). These galaxies are selected from the Calar Alto Legacy Integral Field Area (CALIFA) sample, mapped with optical integral field spectroscopy. EDGE provides good quality CO data (3$\sigma$ sensitivity $\Sigma_{\rm mol}$ $\sim$ 11 M$_\odot$ pc$^{-2}$ before inclination correction, resolution $\sim1.4$ kpc) for 126 galaxies, constituting the largest interferometric CO survey of galaxies in the nearby universe. We describe the survey, the data characteristics, the data products, and present initial science results. We find that the exponential scale-lengths of the molecular, stellar, and star-forming disks are approximately equal, and galaxies that are more compact in molecular gas than in stars tend to show signs of interaction. We characterize the molecular to stellar ratio as a function of Hubble type and stellar mass, present preliminary results on the resolved relations between the molecular gas, stars, and star formation rate, and discuss the dependence of the resolved molecular depletion time on stellar surface density, nebular extinction, and gas metallicity. EDGE provides a key dataset to address outstanding topics regarding gas and its role in star formation and galaxy evolution, which will be publicly available on completion of the quality assessment.

We present results from the MOSFIRE Deep Evolution Field (MOSDEF) survey on the identification, selection biases, and host galaxy properties of 55 X-ray, IR and optically-selected active galactic nuclei (AGN) at $1.4 < z < 3.8$. We obtain rest-frame optical spectra of galaxies and AGN and use the BPT diagram to identify optical AGN. We examine the uniqueness and overlap of the AGN identified at different wavelengths. There is a strong bias against identifying AGN at any wavelength in low mass galaxies, and an additional bias against identifying IR AGN in the most massive galaxies. AGN hosts span a wide range of star formation rate (SFR), similar to inactive galaxies once stellar mass selection effects are accounted for. However, we find (at $\sim 2-3\sigma$ significance) that IR AGN are in less dusty galaxies with relatively higher SFR and optical AGN in dusty galaxies with relatively lower SFR. X-ray AGN selection does not display a bias with host galaxy SFR. These results are consistent with those from larger studies at lower redshifts. Within star-forming galaxies, once selection biases are accounted for, we find AGN in galaxies with similar physical properties as inactive galaxies, with no evidence for AGN activity in particular types of galaxies. This is consistent with AGN being fueled stochastically in any star-forming host galaxy. We do not detect a significant correlation between SFR and AGN luminosity for individual AGN hosts, which may indicate the timescale difference between the growth of galaxies and their supermassive black holes.

We explore the correlations between primordial non-Gaussianity and isocurvature perturbation. We sketch the generic relation between the bispectrum of the curvature perturbation and the cross-correlation power spectrum in the presence of explicit couplings between the inflaton and another light field which gives rise to isocurvature perturbation. Using a concrete model of a Peccei-Quinn type field with generic gravitational couplings, we illustrate explicitly how the primordial bispectrum correlates with the cross-correlation power spectrum. Assuming the resulting fnl ~ O(1), we find that the form of the correlation depends mostly upon the inflation model but only weakly on the axion parameters, even though fnl itself does depend heavily on the axion parameters.

We study the impacts of reheating temperature on the inflationary predictions of the spectral index and tensor-to-scalar ratio. Assuming sinusoidal oscillations and that reheating process is very fast, the reheating temperature can be constrained for sinusoidal oscillation within a factor of 10 - 100 or even better with the prospect of future observations. Beyond this, we find that the predictions can also be insensitive to the reheating temperature in certain models, including Higgs inflation.

Millisecond pulsars (MSPs) are thought to born in low-mass X-ray binaries when the neutron star has gained enough angular momentum from the accreting materials of its companion star. It is generally believed that a radio MSP is born when the neutron star stops accreting and enters a rotation-powered state. Exactly what happens during the transition time was poorly understood until a year ago. In the past year, observations have revealed a few objects that not only switched from one state to the other (as predicted in the above picture), but also have swung between the two states within weeks to years. In this work, we present observations of two of these transition objects (PSR J1023+0038 and XSS J12270-4859) and a theoretical framework that tries to explain their high-energy radiation.

We study a possible connection between processes of gamma-ray emission and hydrogen ionization in a few pc of central region around Sgr A*. Previous investigations showed there is a discrepancy between interpretation of gamma-ray and ionization data if gamma-rays are generated by proton-proton collisions. Here we provided analysis of processes of ionization and emission basing on analytical and numerical calculations of kinetic equations which describe processes of particle propagation and their energy losses. The origin of gamma rays could be either due to collisions of relativistic protons with the dense gas of the surrounding circumnuclear disk (CND) or bremsstrahlung and inverse Compton scattering of relativistic electrons. The hydrogen ionization in this case is produced by a low energy component of the CR spectrum. We found that if ionization is produced by protons the expected ionization rate of hydrogen in the CND is of the same order as derived from IR observations. So we do not see any discrepancy between the gamma-ray and ionization data for the hadronic model. In the case of ionization by electrons we obtained the ionization rate one order of magnitude higher than follows from the IR data. In principle, a selection between the leptonic and hadronic interpretations can be performed basing on measurements of radio and X-ray fluxes from this region because the leptonic and hadronic models give different values of the fluxes from there. We do not exclude that gamma-ray production and hydrogen ionization in the CND are due to a past activity of Sgr A* which occurred about 100 year ago. Then we hypothesize that there may be connection between a past proton eruption and a flux of hard X-rays emitted by Sgr A* hundred years ago as follows from the observed time variability of the iron line seen in the direction of GC molecular clouds.

The First \fermi-LAT Catalog of Sources Above 10 GeV reported evidence of pulsed emission above 25 GeV from 12 pulsars, including the Vela pulsar, which showed evidence of pulsation at $>37$ GeV energy bands. Using 62 months of \fermi-LAT data, we analyzed the gamma-ray emission from the Vela pulsar and searched for pulsed emission above 50 GeV. Having confirmed the significance of the pulsation in 30-50 GeV with the H-test (p-value $\sim10^{-77}$), we extracted its pulse profile using the Bayesian block algorithm and compared it with the distribution of the 5 observed photons above 50 GeV using the likelihood ratio test. Pulsation was significantly detected for photons above 50 GeV with p-value $=3\times10^{-5}$ ($4.2\sigma$). The detection of pulsation is significant above $4\sigma$ at $>79$ GeV and above $3\sigma$ at $>90$ GeV energy bands, making this the highest energy pulsation significantly detected by the LAT. We explore non-stationary outer gap scenario of the very high-energy emissions from the Vela pulsar.

There is a growing new class of young spin-down powered pulsars called GeV-quiet soft gamma-ray pulsar; (1) spectral turnover appears around~10MeV, (2) the X-ray spectra of below 20 keV can be described by power law with photon index around 1.2 and (3) the light curve in X-ray/soft gamma-ray bands shows single broad pulse. Their emission properties are distinct from the normal gamma-ray pulsars, for which the spectral peak in $\nu F_{\nu}$ appears in GeV energy bands and the X-ray/gamma-ray light curves show sharp and double (or more) peaks. In this paper, we discuss that X-ray/soft gamma-ray emissions of the GeV-quiet soft gamma-ray pulsars are caused bythe synchrotron radiation of the electron/positron pairs, which are created by the magnetic pair-creation process near the stellar surface. In our model, the viewing geometry is crucial factor to discriminate between the normal gamma-ray pulsars and soft gamma-ray pulsars. Our model suggests that the difference between the magnetic inclination angle ($\alpha$) and the Earth viewing angle ($\beta$) of the soft gamma-ray pulsars is small, so that the synchrotron emissions from the high magnetic field region around the polar cap region dominates in the observed emissions. Furthermore, the inclination angle of the soft gamma-ray pulsar is relatively small, $\alpha\leq 30$~degree, and our line of sight is out of the gamma-ray beam emitted via the curvature radiation process in the outer gap. We also analysis the six year $Fermi$ data for four soft gamma-ray pulsars to determine the upper limit of the GeV flux.

We study mechanisms of multi-wavelength emissions (X-ray, GeV and TeV gamma-rays) from the gamma-ray binary LS~5039. This paper is composed of two parts. In the first part, we report on results of observational analysis using four year data of \fermi\ Large Area Telescope. Due to the improvement of instrumental response function and increase of the statistics, the observational uncertainties of the spectrum in $\sim$100-300 MeV bands and $>10$GeV bands are significantly improved. The present data analysis suggests that the 0.1-100GeV emissions from LS~5039 contain three different components; (i) the first component contributes to $<$1GeV emissions around superior conjunction, (ii) the second component dominates in 1-10GeV energy bands and (iii) the third component is compatible to lower energy tail of the TeV emissions. In the second part, we develop an emission model to explain the properties of the phase-resolved emissions in multi-wavelength observations. Assuming that LS~5039 includes a pulsar, we argue that both emissions from magnetospheric outer gap and inverse-Compton scattering process of cold-relativistic pulsar wind contribute to the observed GeV emissions. We assume that the pulsar is wrapped by two kinds of termination shock; Shock-I due to the interaction between the pulsar wind and the stellar wind and Shock-II due to the effect of the orbital motion. We propose that the X-rays are produced by the synchrotron radiation at Shock-I region and the TeV gamma-rays are produced by the inverse-Compton scattering process at Shock-II region.

We report on X-ray and gamma-ray observations of the millisecond pulsar B1937+21 taken with the Chandra X-ray Observatory, XMM-Newton, and the Fermi Large Area Telescope. The pulsar X-ray emission shows a purely non-thermal spectrum with a hard photon index of 0.9+/-0.1, and is nearly 100% pulsed. We found no evidence of varying pulse profile with energy as previously claimed. We also analyzed 5.5 yr of Fermi survey data and obtained much improved constraints on the pulsar's timing and spectral properties in gamma-rays. The pulsed spectrum is adequately fitted by a simple power-law with a photon index of 2.38+/-0.07. Both the gamma-ray and X-ray pulse profiles show similar two-peak structure and generally align with the radio peaks. We found that the aligned profiles and the hard spectrum in X-rays seem to be common properties among millisecond pulsars with high magnetic fields at the light cylinder. We discuss a possible physical scenario that could give rise to these features.

Recent observations strongly suggest that the millisecond pulsar binary PSR J1023+0038 has developed an accretion disk since 2013 June. We present a multi-wavelength analysis of PSR J1023+0038, which reveals that 1) its gamma-rays suddenly brightened within a few days in June/July 2013 and has remained at a high gamma-ray state for several months; 2) both UV and X-ray fluxes have increased by roughly an order of magnitude, and 3) the spectral energy distribution has changed significantly after the gamma-ray sudden flux change. Time variabilities associated with UV and X-rays are on the order of 100-500 seconds and 50-100 seconds, respectively. Our model suggests that a newly formed accretion disk due to the sudden increase of the stellar wind could explain the changes of all these observed features. The increase of UV is emitted from the disk, and a new component in gamma-rays is produced by inverse Compton scattering between the new UV component and pulsar wind. The increase of X-rays results from the enhancement of injection pulsar wind energy into the intra-binary shock due to the increase of the stellar wind. We also predict that the radio pulses may be blocked by the evaporated winds from the disk and the pulsar is still powered by rotation.

We study the evolution of the non-linear curvature perturbation during perturbative reheating, and hence how observables evolve to their final values which we may compare against observations. Our study includes the evolution of the two trispectrum parameters, $\gnl$ and $\taunl$, as well as the scale dependence of both $\fnl$ and $\taunl$. In general the evolution is significant and must be taken into account, which means that models of multifield inflation cannot be compared to observations without specifying how the subsequent reheating takes place. If the trispectrum is large at the end of inflation, it normally remains large at the end of reheating. In the classes of models we study, it is very hard to generate $\taunl\gg\fnl^2$, regardless of the decay rates of the fields. Similarly, for the classes of models in which $\gnl\simeq\taunl$ during slow--roll inflation, we find the relation typically remains valid during reheating. Therefore it is possible to observationally test such classes of models without specifying the parameters of reheating, even though the individual observables are sensitive to the details of reheating. It is hard to generate an observably large $\gnl$ however. The runnings, $\nfnl$ and $\ntaunl$, tend to satisfy a consistency relation $\ntaunl=(3/2)\nfnl$, but are in general too small to be observed for the class of models considered regardless of reheating timescale.

Three transiting circumbinary planets (Kepler-16 b, Kepler-34 b, and Kepler-35 b) have recently been discovered from photometric data taken by the Kepler spacecraft. Their orbits are significantly non-Keplerian because of the large secondary-to-primary mass ratio and orbital eccentricity of the binaries, as well as the proximity of the planets to the binaries. We present an analytic theory, with the planet treated as a test particle, which shows that the planetary motion can be represented by the superposition of the circular motion of a guiding center, the forced oscillations due to the non-axisymmetric components of the binary's potential, the epicyclic motion, and the vertical motion. In this analytic theory, the periapse and ascending node of the planet precess at nearly equal rates in opposite directions. The largest forced oscillation term corresponds to a forced eccentricity (which is an explicit function of the parameters of the binary and of the guiding center radius of the planet), and the amplitude of the epicyclic motion (which is a free parameter of the theory) is the free eccentricity. Comparisons with direct numerical orbit integrations show that this analytic theory gives an accurate description of the planetary motion for all three Kepler systems. We find that all three Kepler circumbinary planets have nonzero free eccentricities.

We study the effects of perturbative reheating on the evolution of the curvature perturbation \zeta, in two-field inflation models. We use numerical methods to explore the sensitivity of f_NL, n_s and r to the reheating process, and present simple qualitative arguments to explain our results. In general, if a large non-Gaussian signal exists at the start of reheating, it will remain non zero at the end of reheating. Unless all isocurvature modes have completely decayed before the start of reheating, we find that the non-linearity parameter, f_NL, can be sensitive to the reheating timescale, and that this dependence is most appreciable for `runaway' inflationary potentials that only have a minimum in one direction. For potentials with a minimum in both directions, f_NL can also be sensitive to reheating if a mild hierarchy exists between the decay rates of each field. Within the class of models studied, we find that the spectral index n_s, is fairly insensitive to large changes in the field decay rates, indicating that n_s is a more robust inflationary observable, unlike the non-linearity parameter f_NL. Our results imply that the statistics of \zeta, especially f_NL, can only be reliably used to discriminate between models of two-field inflation if the physics of reheating are properly accounted for.