The integrated X-ray luminosity ($L_{\mathrm{X}}$) of high-mass X-ray binaries (HMXBs) in a galaxy is correlated with its star formation rate (SFR), and the normalization of this correlation increases with redshift. Population synthesis models suggest that the redshift evolution of $L_{\mathrm{X}}$/SFR is driven by the metallicity ($Z$) dependence of HMXBs, and the first direct evidence of this connection was recently presented using galaxies at $z\sim2$. To confirm this result with more robust measurements and better constrain the $L_{\mathrm{X}}$-SFR-$Z$ relation, we have studied the $Z$ dependence of $L_{\mathrm{X}}$/SFR at lower redshifts. Using samples of star-forming galaxies at $z=0.1-0.9$ with optical spectra from the hCOSMOS and zCOSMOS surveys, we stacked \textit{Chandra} data from the COSMOS Legacy survey to measure the average $L_{\mathrm{X}}$/SFR as a function of $Z$ in three redshift ranges: $z=0.1-0.25$, $0.25-0.4$, and $0.5-0.9$. We find no significant variation of the $L_{\mathrm{X}}$-SFR-$Z$ relation with redshift. Our results provide further evidence that the $Z$ dependence of HMXBs is responsible for the redshift evolution of $L_{\mathrm{X}}$/SFR. Combining all available $z>0$ measurements together, we derive a best-fitting $L_{\mathrm{X}}$-SFR-$Z$ relation and assess how different population synthesis models describe the data. These results provide the strongest constraints to date on the $L_{\mathrm{X}}$-SFR-$Z$ relation in the range of $8.0<$12+log(O/H)$<9.0$.

The International Gamma-Ray Astrophysics Laboratory (INTEGRAL) satellite has detected in excess of 1000 sources in the ~20-100 keV band during its surveys of the sky over the past 17 years. We obtained 5 ks observations of 15 unclassified IGR sources with the Chandra X-ray Observatory in order to localize them, to identify optical/IR counterparts, to measure their soft X-ray spectra, and to classify them. For 10 of the IGR sources, we detect Chandra sources that are likely (or in some cases certain) to be the counterparts. IGR J18007-4146 and IGR J15038-6021 both have Gaia parallax distances, placing them at 2.5+0.5-0.4 and 1.1+1.5-0.4 kpc, respectively. We tentatively classify both of them as intermediate polar-type Cataclysmic Variables. Also, IGR J17508-3219 is likely to be a Galactic source, but it is unclear if it is a Dwarf Nova or another type of transient. For IGR J17118-3155, we provide a Chandra localization, but it is unclear if the source is Galactic or extragalactic. Based on either near-IR/IR colors or the presence of extended near-IR emission, we classify four sources as Active Galactic Nuclei (IGR J16181-5407, IGR J16246-4556, IGR J17096-2527, and IGR J19294+1327), and IGR J20310+3835 and IGR J15541-5613 are AGN candidates. In addition, we identified an AGN in the INTEGRAL error circle of IGR J16120-3543 that is a possible counterpart.

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.

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.

We present a source catalog from the first deep hard X-ray ($E>10$ keV) survey of the Small Magellanic Cloud (SMC), the NuSTAR Legacy Survey of the SMC. We observed three fields, for a total exposure time of 1 Ms, along the bar of this nearby star-forming galaxy. Fields were chosen for their young stellar and accreting binary populations. We detected 10 sources above a 3$\sigma$ significance level (4$-$25 keV) and obtained upper limits on an additional 40 sources. We reached a 3$\sigma$ limiting luminosity in the 4$-$25 keV band of $\sim$ $10^{35}$ erg s$^{-1}$, allowing us to probe fainter X-ray binary (XRB) populations than has been possible with other extragalactic NuSTAR surveys. We used hard X-ray colors and luminosities to constrain the compact-object type, exploiting the spectral differences between accreting black holes and neutron stars at $E>10$ keV. Several of our sources demonstrate variability consistent with previously observed behavior. We confirmed pulsations for seven pulsars in our 3$\sigma$ sample. We present the first detection of pulsations from a Be-XRB, SXP305 (CXO J005215.4$-$73191), with an X-ray pulse period of $305.69\pm0.16$ seconds and a likely orbital period of $\sim$1160-1180 days. Bright sources ($\gtrsim 5\times 10^{36}$ erg s$^{-1}$) in our sample have compact-object classifications consistent with their previously reported types in the literature. Lower luminosity sources ($\lesssim 5\times 10^{36}$ erg s$^{-1}$) have X-ray colors and luminosities consistent with multiple classifications. We raise questions about possible spectral differences at low luminosity between SMC pulsars and the Galactic pulsars used to create the diagnostic diagrams.

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.

The cosmic X-ray background (CXB), which peaks at an energy of ~30 keV, is produced primarily by emission from accreting supermassive black holes (SMBHs). The CXB therefore serves as a constraint on the integrated SMBH growth in the Universe and the accretion physics and obscuration in active galactic nuclei (AGNs). This paper gives an overview of recent progress in understanding the high-energy (>~10 keV) X-ray emission from AGNs and the synthesis of the CXB, with an emphasis on results from NASA's NuSTAR hard X-ray mission. We then discuss remaining challenges and open questions regarding the nature of AGN obscuration and AGN physics. Finally, we highlight the exciting opportunities for a next-generation, high-resolution hard X-ray mission to achieve the long-standing goal of resolving and characterizing the vast majority of the accreting SMBHs that produce the CXB.

The discoveries made over the past 20 years by Chandra and XMM-Newton surveys in conjunction with multiwavelength imaging and spectroscopic data available in the same fields have significantly changed the view of the supermassive black hole (SMBH) and galaxy connection. These discoveries have opened up several exciting questions that are beyond the capabilities of current X-ray telescopes and will need to be addressed by observatories in the next two decades. As new observatories peer into the early Universe, we will begin to understand the physics and demographics of SMBH infancy (at $z>6$) and investigate the influence of their accretion on the formation of the first galaxies ($\S$ 2.1). We will also be able to understand the accretion and evolution over the cosmic history (at $z\sim$1-6) of the full population of black holes in galaxies, including low accretion rate, heavily obscured AGNs at luminosities beyond the reach of current X-ray surveys ($\S$2.2 and $\S$2.3), enabling us to resolve the connection between SMBH growth and their environment.

What are the most important conditions and processes governing the growth of stellar-origin compact objects? The identification of compact object type as either black hole (BH) or neutron star (NS) is fundamental to understanding their formation and evolution. To date, time-domain determination of compact object type remains a relatively untapped tool. Measurement of orbital periods, pulsations, and bursts will lead to a revolution in the study of the demographics of NS and BH populations, linking source phenomena to accretion and galaxy parameters (e.g., star formation, metallicity). To perform these measurements over sufficient parameter space, a combination of a wide-field (>5000 deg^2) transient X-ray monitor over a dynamic energy range (~1-100 keV) and an X-ray telescope for deep surveys with <5 arcsec PSF half-energy width (HEW) angular resolution are required. Synergy with multiwavelength data for characterizing the underlying stellar population will transform our understanding of the time domain properties of transient sources, helping to explain details of supernova explosions and gravitational wave event rates.

We write in response to the call from the 2020 Decadal Survey to submit white papers illustrating the most pressing scientific questions in astrophysics for the coming decade. We propose exploration as the central question for the Decadal Committee's discussions.The history of astronomy shows that paradigm changing discoveries are not driven by well formulated scientific questions, based on the knowledge of the time. They were instead the result of the increase in discovery space fostered by new telescopes and instruments. An additional tool for increasing the discovery space is provided by the analysis and mining of the increasingly larger amount of archival data available to astronomers. Revolutionary observing facilities, and the state of the art astronomy archives needed to support these facilities, will open up the universe to new discovery. Here we focus on exploration for compact objects and multi messenger science. This white paper includes science examples of the power of the discovery approach, encompassing all the areas of astrophysics covered by the 2020 Decadal Survey.

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.

The NuSTAR serendipitous survey has already uncovered a large number of Active Galactic Nuclei (AGN), providing new information about the composition of the Cosmic X-ray Background. For the AGN off the Galactic plane, it has been possible to use the existing X-ray archival data to improve source localizations, identify optical counterparts, and classify the AGN with optical spectroscopy. However, near the Galactic Plane, better X-ray positions are necessary to achieve optical or near-IR identifications due to the higher levels of source crowding. Thus, we have used observations with the Chandra X-ray Observatory to obtain the best possible X-ray positions. With eight observations, we have obtained coverage for 19 NuSTAR serendips within 12 deg of the plane. One or two Chandra sources are detected within the error circle of 15 of the serendips, and we report on these sources and search for optical counterparts. For one source (NuSTAR J202421+3350.9), we obtained a new optical spectrum and detected the presence of hydrogen emission lines. The source is Galactic, and we argue that it is likely a Cataclysmic Variable. For the other sources, the Chandra positions will enable future classifications in order to place limits on faint Galactic populations, including high-mass X-ray binaries and magnetars.

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.

Nearby galaxy surveys have long classified X-ray binaries (XRBs) by the mass category of their donor stars (high-mass and low-mass). The NuSTAR observatory, which provides imaging data at E $>10$ keV, has enabled the classification of extragalactic XRBs by their compact object type: neutron star (NS) or black hole (BH). We analyzed NuSTAR/Chandra/XMM-Newton observations from a NuSTAR-selected sample of 12 galaxies within 5 Mpc having stellar masses ($M_{\star}$) $10^{7-11}$ $M_{\odot}$ and star formation rates (SFR) $\approx0.01-15$ $M_{\odot}$ yr$^{-1}$. We detect 128 NuSTAR sources to a sensitivity of $\approx10^{38}$ erg s$^{-1}$. Using NuSTAR color-intensity and color-color diagrams we classify 43 of these sources as candidate NS and 47 as candidate BH. We further subdivide BH by accretion states (soft, intermediate, and hard) and NS by weak (Z/Atoll) and strong (accreting pulsar) magnetic field. Using 8 normal (Milky Way-type) galaxies in the sample, we confirm the relation between SFR and galaxy X-ray point source luminosity in the 4-25 and 12-25 keV energy bands. We also constrain galaxy X-ray point source luminosity using the relation $L_{\rm{X}}=\alpha M_{\star}+\beta\text{SFR}$, finding agreement with previous work. The XLF of all sources in the 4-25 and 12-25 keV energy bands matches with the $\alpha=1.6$ slope for high-mass XRBs. We find that NS XLFs suggest a decline beginning at the Eddington limit for a 1.4 $M_{\odot}$ NS, whereas the BH fraction shows an approximate monotonic increase in the 4-25 and 12-25keV energy bands. We calculate the overall ratio of BH to NS to be $\approx1$ for 4-25 keV and $\approx2$ for 12-25 keV.

We present an X-ray stacking analysis of $\sim$75,000 star-forming galaxies between $0.1<z<5.0$ using the Chandra COSMOS Legacy survey to study the X-ray emission of low-luminosity active galactic nuclei (AGN) and its connection to host galaxy properties. The stacks at $z<0.9$ have luminosity limits as low as $10^{40}-10^{41}$ erg s$^{-1}$, a regime in which X-ray binaries (XRBs) can dominate the X-ray emission. Comparing the measured luminosities to established XRB scaling relations, we find that the redshift evolution of the luminosity per star formation rate (SFR) of XRBs depends sensitively on the assumed obscuration and may be weaker than previously found. The XRB scaling relation based on stacks from the Chandra Deep Field South overestimates the XRB contribution to the COSMOS high specific SFR (sSFR) stacks, possibly due to a bias affecting the CDF-S stacks because of their small galaxy samples. After subtracting the estimated XRB contribution from the stacks, we find that most stacks at $z>1.3$ exhibit a significant X-ray excess indicating nuclear emission. The AGN emission is strongly correlated with stellar mass but does not exhibit an additional correlation with SFR. The hardness ratios of the high-redshift stacks indicate that the AGN are substantially obscured ($N_{\mathrm{H}}\sim10^{23}$ cm$^{-2}$). These obscured AGN are not identified by IRAC color selection and have $L_X\sim10^{41}-10^{43}$ erg s$^{-1}$, consistent with accretion at an Eddington rate of $\sim10^{-3}$ onto $10^7-10^8 M_{\odot}$ black holes. Combining our results with other X-ray studies suggests that AGN obscuration depends on stellar mass and an additional variable, possibly the Eddington rate.

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 the results and the source catalog of the NuSTAR survey in the UKIDSS Ultra Deep Survey (UDS) field, bridging the gap in depth and area between NuSTAR's ECDFS and COSMOS surveys. The survey covers a $\sim 0.6$ deg$^2$ area of the field for a total observing time of $\sim$ 1.75 Ms, to a half-area depth of $\sim$ 155 ks corrected for vignetting at $3-24$ keV, and reaching sensitivity limits at half-area in the full ($3-24$ keV), soft ($3-8$ keV) and hard ($8-24$ keV) bands of $2.2 \times 10^{-14}$ erg cm$^{-2}$ s$^{-1}$, $1.0 \times 10^{-14}$ erg cm$^{-2}$ s$^{-1}$, and $2.7 \times 10^{-14}$ erg cm$^{-2}$ s$^{-1}$, respectively. A total of 67 sources are detected in at least one of the three bands, 56 of which have a robust optical redshift with a median of $\langle z\rangle \sim 1.1$. Through a broadband ($0.5-24$ keV) spectral analysis of the whole sample combined with the NuSTAR hardness ratios, we compute the observed Compton-thick (CT; $N_{\rm H} > 10^{24}$ cm$^{-2}$) fraction. Taking into account the uncertainties on each $N_{\rm H}$ measurement, the final number of CT sources is $6.8\pm1.2$. This corresponds to an observed CT fraction of $(11.5\pm2.0)\%$, providing a robust lower limit to the intrinsic fraction of CT AGN and placing constraints on cosmic X-ray background synthesis models.

IGR J18214-1318, a Galactic source discovered by the International Gamma-Ray Astrophysics Laboratory, is a high-mass X-ray binary (HMXB) with a supergiant O-type stellar donor. We report on the XMM-Newton and NuSTAR observations that were undertaken to determine the nature of the compact object in this system. This source exhibits high levels of aperiodic variability, but no periodic pulsations are detected with a 90% confidence upper limit of 2% fractional rms between 0.00003-88 Hz, a frequency range that includes the typical pulse periods of neutron stars (NSs) in HMXBs (0.1-10$^3$ s). Although the lack of pulsations prevents us from definitively identifying the compact object in IGR J18214-1318, the presence of an exponential cutoff with e-folding energy $\lesssim30$ keV in its 0.3-79 keV spectrum strongly suggests that the compact object is an NS. The X-ray spectrum also shows a Fe K$\alpha$ emission line and a soft excess, which can be accounted for by either a partial-covering absorber with $N_{\mathrm{H}}\approx10^{23}$ cm$^{-2}$ which could be due to the inhomogeneous supergiant wind, or a blackbody component with $kT=1.74^{+0.04}_{-0.05}$ keV and $R_{BB}\approx0.3$ km, which may originate from NS hot spots. Although neither explanation for the soft excess can be excluded, the former is more consistent with the properties observed in other supergiant HMXBs. We compare IGR J18214-1318 to other HMXBs that lack pulsations or have long pulsation periods beyond the range covered by our observations.

The Nuclear Spectroscopic Telescope Array (NuSTAR) provides an improvement in sensitivity at energies above 10 keV by two orders of magnitude over non-focusing satellites, making it possible to probe deeper into the Galaxy and Universe. Lansbury and collaborators recently completed a catalog of 497 sources serendipitously detected in the 3-24 keV band using 13 deg2 of NuSTAR coverage. Here, we report on an optical and X-ray study of 16 Galactic sources in the catalog. We identify eight of them as stars (but some or all could have binary companions), and use information from Gaia to report distances and X-ray luminosities for three of them. There are four CVs or CV candidates, and we argue that NuSTAR J233426-2343.9 is a relatively strong CV candidate based partly on an X-ray spectrum from XMM-Newton. NuSTAR J092418-3142.2, which is the brightest serendipitous source in the Lansbury catalog, and NuSTAR J073959-3147.8 are LMXB candidates, but it is also possible that these two sources are CVs. One of the sources is a known HMXB, and NuSTAR J105008-5958.8 is a new HMXB candidate, which has strong Balmer emission lines in its optical spectrum and a hard X-ray spectrum. We discuss the implications of finding these HMXBs for the surface density (logN-logS) and luminosity function of Galactic HMXBs. We conclude that, with the large fraction of unclassified sources in the Galactic plane detected by NuSTAR in the 8-24 keV band, there could be a significant population of low luminosity HMXBs.

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

Studying the population of faint hard X-ray sources along the plane of the Galaxy is challenging because of high-extinction and crowding, which make the identification of individual sources more difficult. IGR J18293-1213 is part of the population of persistent sources which have been discovered by the INTEGRAL satellite. We report on NuSTAR and Swift/XRT observations of this source, performed on 2015 September 11. We detected three eclipsing intervals in the NuSTAR light curve, allowing us to constrain the duration of these eclipses, $\Delta t = 30.8^{+6.3}_{-0.0}$ min, and the orbital period of the system, $T = 6.92\pm0.01$ hr. Even though we only report an upper limit on the amplitude of a putative spin modulation, the orbital period and the hard thermal Bremsstrahlung spectrum of IGR J18293-1213 provide strong evidence that this source is a magnetic Cataclysmic Variable (CV). Our NuSTAR and Swift/XRT joint spectral analysis places strong constraints on the white dwarf mass $M_{\rm wd} = 0.78^{+0.10}_{-0.09}$ M$_\odot$. Assuming that the mass to radius ratio of the companion star $M_* / R_* = 1$ (solar units) and using $T$, $\Delta t$ and $M_{\rm wd}$, we derived the mass of the companion star $M_* = 0.82\pm0.01$ M$_\odot$, the orbital separation of the binary system $a=2.14\pm0.04$ R$_\odot$, and its orbital inclination compared to the line of sight $i=(72.2^{+2.4}_{-0.0})\pm1.0^\circ$.

Recently, unresolved hard (20-40 keV) X-ray emission has been discovered within the central 10 pc of the Galaxy, possibly indicating a large population of intermediate polars (IPs). Chandra and XMM-Newton measurements in the surrounding ~50 pc imply a much lighter population of IPs with $\langle M_{\rm WD} \rangle \approx 0.5 M_\odot$. Here we use broad-band NuSTAR observations of two IPs: TV Columbae, which has a fairly typical but widely varying reported mass of $M_{\rm WD} \approx 0.5-1.0 M_\odot$, and IGR J17303-0601, with a heavy reported mass of $M_{\rm WD} \approx 1.0-1.2 M_\odot$. We investigate how varying spectral models and observed energy ranges influence estimated white dwarf mass. Observations of the inner 10 pc can be accounted for by IPs with $\langle M_{\rm WD} \rangle \approx 0.9 M_\odot$, consistent with that of the CV population in general, and the X-ray observed field IPs in particular. The lower mass derived by Chandra and XMM-Newton appears to be an artifact of narrow energy band fitting. To explain the (unresolved) CHXE by IPs requires an X-ray (2-8 keV) luminosity function (XLF) extending down to at least $5\times10^{31}$ erg/s. The CHXE XLF, if extended to the surrounding ~50 pc observed by Chandra and XMM-Newton, requires at least ~20-40% of the $\sim$9000 point sources are IPs. If the XLF extends just a factor of a few lower in luminosity, then the vast majority of these sources are IPs. This is in contrast to recent observations of the Galactic ridge, where the bulk of the 2-8 keV emission is ascribed to dwarf novae.

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

The high-mass X-ray binary and accreting X-ray pulsar IGR J16393-4643 was observed by NuSTAR in the 3-79 keV energy band for a net exposure time of 50 ks. We present the results of this observation which enabled the discovery of a cyclotron resonant scattering feature with a centroid energy of 29.3(+1.1/-1.3) keV. This allowed us to measure the magnetic field strength of the neutron star for the first time: B = (2.5+/-0.1)e12 G. The known pulsation period is now observed at 904.0+/-0.1 s. Since 2006, the neutron star has undergone a long-term spin-up trend at a rate of P' = -2e-8 s/s (-0.6 s per year, or a frequency derivative of nu' = 3e-14 Hz/s ). In the power density spectrum, a break appears at the pulse frequency which separates the zero slope at low frequency from the steeper slope at high frequency. This addition of angular momentum to the neutron star could be due to the accretion of a quasi-spherical wind, or it could be caused by the transient appearance of a prograde accretion disk that is nearly in corotation with the neutron star whose magnetospheric radius is around 2e8 cm.

We report on 0.3-10 keV observations with the Chandra X-ray Observatory of eight hard X-ray sources discovered within 8 degrees of the Galactic plane by the INTEGRAL satellite. The short (5 ks) Chandra observations of the IGR source fields have yielded very likely identifications of X-ray counterparts for three of the IGR sources: IGR J14091-6108, IGR J18088-2741, and IGR J18381-0924. The first two have very hard spectra in the Chandra band that can be described by a power-law with photon indices of Gamma = 0.6+/-0.4 and -0.7(+0.4)(-0.3), respectively (90% confidence errors are given), and both have a unique near-IR counterpart consistent with the Chandra position. IGR J14091-6108 also displays a strong iron line and a relatively low X-ray luminosity, and we argue that the most likely source type is a Cataclysmic Variable (CV), although we do not completely rule out the possibility of a High Mass X-ray Binary. IGR J18088-2741 has an optical counterpart with a previously measured 6.84 hr periodicity, which may be the binary orbital period. We also detect five cycles of a possible 800-950 s period in the Chandra light curve, which may be the compact object spin period. We suggest that IGR J18088-2741 is also most likely a CV. For IGR J18381-0924, the spectrum is intrinsically softer with Gamma = 1.5(+0.5)(-0.4), and it is moderately absorbed, nH = (4+/-1)e22 cm-2. There are two near-IR sources consistent with the Chandra position, and they are both classified as galaxies, making it likely that IGR J18381-0924 is an Active Galactic Nucleus (AGN). For the other five IGR sources, we provide lists of nearby Chandra sources, which may be used along with further observations to identify the correct counterparts, and we discuss the implications of the low inferred Chandra count rates for these five sources.

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

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

We report on CTIO/NEWFIRM and CTIO/OSIRIS photometric and spectroscopic observations of 20 new X-ray (0.5-10 keV) emitters discovered in the Norma Arm Region Chandra Survey (NARCS). NEWFIRM photometry was obtained to pinpoint the near-infrared counterparts of NARCS sources, while OSIRIS spectroscopy was used to help identify 20 sources with possible high mass X-ray binary properties. We find that (1) two sources are WN8 Wolf-Rayet stars, maybe in colliding wind binaries, part of the massive star cluster Mercer 81; (2) two are emission-line stars, possibly in X-ray binaries, that exhibit near- and mid-infrared excesses either due to free-free emission from the decretion discs of Be stars or warm dust in the stellar winds of peculiar massive stars such as B[e] supergiants or luminous blue variables; (3) one is a B8-A3 IV-V star that could be in a quiescent high mass X-ray binary system; (4) two are cataclysmic variables including one intermediate polar; (5) three may be neutron star symbiotic binaries; (6) five are most likely white dwarf symbiotic binaries; and (7) five exhibit properties more consistent with isolated giant/dwarf stars. The possible detection of one to three high mass X-ray binaries is in good agreement with our predictions. However, our study illustrates the difficulty of clearly differentiating quiescent or intermediate X-ray luminosity systems from isolated massive stars, which may lead to an underestimation of the number of known high mass X-ray binaries.

During hard X-ray observations of the Norma spiral arm region by the Nuclear Spectroscopic Telescope Array (NuSTAR) in 2013 February, a new transient source, NuSTAR J163433-4738.7, was detected at a significance level of 8-sigma in the 3-10 keV bandpass. The source is consistent with having a constant NuSTAR count rate over a period of 40 ks and is also detected simultaneously by Swift at lower significance. The source is not significantly detected by NuSTAR, Swift, or Chandra in the days before or weeks after the discovery of the transient, indicating that the strong X-ray activity lasted for between ~0.5 and 1.5 days. Near-IR imaging observations were carried out before and after the X-ray activity, but we are not able to identify the counterpart. The combined NuSTAR and Swift energy spectrum is consistent with a power-law with a photon index of Gamma = 4.1(+1.5)(-1.0) (90% confidence errors), a blackbody with kT = 1.2+/-0.3 keV, or a bremsstrahlung model with kT = 3.0(+2.1)(-1.2) keV. The reduced-chi2 values for the three models are not significantly different, ranging from 1.23 to 1.44 for 8 degrees of freedom. The spectrum is strongly absorbed with NH = 2.8(+2.3)(-1.4)e23 cm-2, 9(+15)(-7)e22 cm-2, and 1.7(+1.7)(-0.9)e23 cm-2, for the power-law, blackbody, and bremsstrahlung models, respectively. Although the high column density could be due to material local to the source, it is consistent with absorption from interstellar material along the line of sight at a distance of 11 kpc, which would indicate an X-ray luminosity >1e34 erg/s. Although we do not reach a definitive determination of the nature of NuSTAR J163433-4738.7, we suggest that it may be an unusually bright active binary or a magnetar.

We report on Chandra X-ray and Parkes radio observations of IGR J11014-6103, which is a possible pulsar wind nebula with a complex X-ray morphology and a likely radio counterpart. With the superb angular resolution of Chandra, we find evidence that a portion of the extended emission may be related to a bow shock due to the putative pulsar moving through the interstellar medium. The inferred direction of motion is consistent with IGR J11014-6103 having been born in the event that produced the supernova remnant (SNR) MSH 11-61A. If this association is correct, then previous constraints on the expansion of MSH 11-61A imply a transverse velocity for IGR J11014-6103 of 2,400-2,900 km/s, depending on the SNR model used. This would surpass the kick velocities of any known pulsars and rival or surpass the velocities of any compact objects that are associated with SNRs. While it is important to confirm the nature of the source, our radio pulsation search did not yield a detection.