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.

In this paper I describe numerical calculations of the motion of particles in a disk about a solar-mass object perturbed by a planet on a circular orbit with mass greater than 0.001 of the stellar mass. A simple algorithm for simulating bulk viscosity is added to the ensemble of particles, and the response of the disk is followed for several planet orbital periods. A two-arm spiral structure forms near the inner resonance (2-1) and extends to the planetary orbit radius (corotation). In the same way for gaseous disks on a galactic scale perturbed by a weak rotating bar-like distortion, this is shown to be related to the appearance of two perpendicular families of periodic orbits near the resonance combined with dissipation which inhibits the crossing of streamlines. Spiral density enhancements result from the crowding of streamlines due to the gradual shift between families. The results, such as the dependence of pitch-angle on radius and the asymmetry of the spiral features, resemble those of sophisticated calculations that include more physical effects. The morphology of structure generated in this way clearly resembles that observed in objects with well-defined two-arm spirals, such as SAO 206462. This illustrates that the process of spiral formation via interaction with planets in such disks can be due to orbital motion in a perturbed Keplerian field combined with kinematic viscosity.}

Physical properties of galaxies at z>7 are of interest for understanding both the early phases of star formation and the process of cosmic reionization. Chemical abundance measurements offer valuable information on the integrated star formation history, and hence ionizing photon production, as well as the rapid gas accretion expected at such high redshifts. We use reported measurements of [O III] 88$\mu$m emission and star formation rate to estimate gas-phase oxygen abundances in five galaxies at z=7.1-9.1 using the direct T_e method. We find typical abundances 12+log(O/H) = 7.9 ($\sim$0.2 times the solar value) and an evolution of 0.9$\pm$0.5 dex in oxygen abundance at fixed stellar mass from z$\simeq$8 to 0. These results are compatible with theoretical predictions, albeit with large (conservative) uncertainties in both mass and metallicity. We assess both statistical and systematic uncertainties to identify promising means of improvement with the Atacama Large Millimeter Array (ALMA) and the James Webb Space Telescope (JWST). In particular we highlight [O III] 52$\mu$m as a valuable feature for robust metallicity measurements. Precision of 0.1-0.2 dex in T_e-based O/H abundance can be reasonably achieved for galaxies at z$\approx$5-8 by combining [O III] 52$\mu$m with rest-frame optical strong lines. It will also be possible to probe gas mixing and mergers via resolved T_e-based abundances on kpc scales. With ALMA and JWST, direct metallicity measurements will thus be remarkably accessible in the reionization epoch.

Galaxy evolution is thought to be driven in large part by the flow of gas between galaxies and the circumgalactic medium (CGM), a halo of metal-enriched gas extending out to $\gtrsim100$ kpc from each galaxy. Studying the spatial structure of the CGM holds promise for understanding these gas flow mechanisms; however, the common method using background quasar sightlines provides minimal spatial information. Recent works have shown the utility of extended background sources such as giant gravitationally lensed arcs. Using background lensed arcs from the CSWA 38 lens system, we continuously probed, at a resolution element of about 15 kpc$^2$, the spatial and kinematic distribution of MgII absorption in a star-forming galaxy at $z=0.77$ (stellar mass $\approx 10^{9.7}$ M$_\odot$, star formation rate $\approx 10$ M$_\odot$ yr$^{-1}$) at impact parameters $D=5-40$ kpc. Our results present an anisotropic, optically thick medium whose absorption strength decreases with increasing impact parameter, in agreement with the statistics towards quasars and other gravitational arcs. Furthermore, we detect low line-of-sight velocities ($v\approx-25-80$ km s$^{-1}$) and relatively high velocity dispersion ($\sigma\approx50\pm30$ km s$^{-1}$) in the MgII gas. These measures provide evidence of a mainly pressure-supported, metal-enriched diffuse gas recycling through the CGM rather than an active galactic outflow.

We present a joint analysis of rest-UV and rest-optical spectra obtained using Keck/LRIS and Keck/MOSFIRE for a sample of 62 star-forming galaxies at $z\sim2.3$. We divide our sample into 2 bins based on their location in the [OIII]/Hb vs. [NII]/Ha BPT diagram, and perform the first differential study of the rest-UV properties of massive ionizing stars as a function of rest-optical emission-line ratios. Fitting BPASS stellar population synthesis models, including nebular continuum emission, to our rest-UV spectra, we find that high-redshift galaxies offset towards higher [OIII]/Hb and [NII]/Ha have younger ages (log(Age/yr)=$7.20^{+0.57}_{-0.20}$) and lower stellar metallicities ($Z_*=0.0010^{+0.0011}_{-0.0003}$) resulting in a harder ionizing spectrum, compared to the galaxies in our sample that lie on the local BPT star-forming sequence (log(Age/yr)=$8.57^{+0.88}_{-0.84}$, $Z_*=0.0019\pm0.0006$). Additionally, we find that the offset galaxies have an ionization parameter of log(U)=$-3.04^{+0.06}_{-0.11}$ and nebular metallicity of 12+log(O/H)=$8.40^{+0.06}_{-0.07}$, and the non-offset galaxies have log(U)=$-3.11\pm0.08$ and 12+log(O/H)=$8.30^{+0.05}_{-0.06}$. The stellar and nebular metallicities derived for our sample imply that the galaxies offset from the local BPT relation are more alpha-enhanced ($7.28^{+2.52}_{-2.82}$ O/Fe$_\odot$) compared to those consistent with the local sequence ($3.04^{+0.95}_{-0.54}$ O/Fe$_\odot$). However, even galaxies that are entirely consistent with the local nebular excitation sequence appear to be alpha-enhanced -- in contrast with typical local systems. Such differences must be considered when estimating gas-phase O/H at high redshift based on strong emission-line ratios. Specifically, a similarity in the location of high-redshift and local galaxies in the BPT diagram may not be indicative of a similarity in their physical properties.

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 dark energy-cold dark matter paradigm ($\Lambda$CDM) has gained widespread acceptance because it explains the pattern of anisotropies observed in the cosmic microwave background radiation, the observed distribution of large scale inhomogeneities in detectable matter, and the perceived overall expansion history of the Universe. It is further {\it assumed} that the cosmic dark matter component clusters on the scale of bound astronomical systems and thereby accounts for the observed difference between the directly detectable (baryonic) mass and the total Newtonian dynamical mass. In this respect the paradigm fails; it is falsified by the existence of a simple algorithm, modified Newtonian dynamics (MOND), which explains, not only general scaling relations for astronomical systems, but quite precisely predicts the effective gravitational acceleration in such objects from the observed distribution of detectable baryonic matter -- all of this with one additional universal parameter having units of acceleration. On this sub-Hubble scale, the dark matter hypothesis is essentially reactive, while MOND is successfully predictive.

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.

We present the discovery of PS18kh, a tidal disruption event (TDE) discovered at the center of SDSS J075654.53+341543.6 ($d\simeq322$ Mpc) by the Pan-STARRS Survey for Transients. Our dataset includes pre-discovery survey data from Pan-STARRS, the All-Sky Automated Survey for Supernovae (ASAS-SN), and the Asteroid Terrestrial-impact Last Alert System (ATLAS) as well as high-cadence, multi-wavelength follow-up data from ground-based telescopes and Swift, spanning from 56 days before peak light until 75 days after. The optical/UV emission from PS18kh is well-fit as a blackbody with temperatures ranging from $T\simeq12000$ K to $T\simeq25000$ K and it peaked at a luminosity of $L\simeq8.8\times10^{43}$ ergs s$^{-1}$. PS18kh radiated $E=(3.45\pm0.22)\times10^{50}$ ergs over the period of observation, with $(1.42\pm0.20)\times10^{50}$ ergs being released during the rise to peak. Spectra of PS18kh show a changing, boxy/double-peaked H$\alpha$ emission feature, which becomes more prominent over time. We use models of non-axisymmetric accretion disks to describe the profile of the H$\alpha$ line and its evolution. We find that at early times the high accretion rate leads the disk to emit a wind which modifies the shape of the line profile and makes it bell-shaped. At late times, the wind becomes optically thin, allowing the non-axisymmetric perturbations to show up in the line profile. The line-emitting portion of the disk extends from $r_{\rm in}\sim60r_{\rm g}$ to an outer radius of $r_{\rm out}\sim1400r_{\rm g}$ and the perturbations can be represented either as an eccentricity in the outer rings of the disk or as a spiral arm in the inner disk.

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.

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.

(Abridged) Long gamma-ray bursts (LGRB) have been suggested as promising tracers of star formation owing to their association with the core-collapse of massive stars. The goal of this work is to characterise the population of host galaxies of LGRBs at 1 < z < 2, investigate the conditions in which LGRBs form at these redshifts and assess their use as tracers of star formation. We perform a spectro-photometric analysis to determine the stellar mass, star formation rate, specific star formation rate and metallicity of the complete, unbiased host galaxy sample of the Swift/BAT6 LGRB sample at 1 < z < 2. We compare the distribution of these properties to the ones of typical star-forming galaxies from the MOSDEF and COSMOS2015 Ultra Deep surveys, within the same redshift range. We find that, similarly to z < 1, LGRBs do not directly trace star formation at 1 < z < 2, and they tend to avoid high-mass, high-metallicity host galaxies. We also find evidence for an enhanced fraction of starbursts among the LGRB host sample with respect to the star-forming population of galaxies. Nonetheless we demonstrate that the driving factor ruling the LGRB efficiency is metallicity. The LGRB host distributions can be reconciled with the ones expected from galaxy surveys by imposing a metallicity upper limit of 12+logOH ~ 8.55. Metallicity rules the LGRB production efficiency, which is stifled at Z > 0.7 Zsun. Under this hypothesis we can expect LGRBs to trace star formation at z > 3, once the bulk of the star forming galaxy population are characterised by metallicities below this limit. The moderately high metallicity threshold found is in agreement with the conditions necessary to rapidly produce a fast-rotating Wolf-Rayet star a in close binary system, and could be accommodated by single star models under chemically homogeneous mixing with very rapid rotation and weak magnetic coupling.

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.

I consider the observed rotation curves of 12 gas-dominated low-surface-brightness galaxies -- objects in which the mass of gas ranges between 2.2 and 27 times the mass of the stellar disk (mean=9.4). This means that, in the usual decomposition of rotation curves into those resulting from various mass components, the mass-to-light ratio of the luminous stellar disk effectively vanishes as an additional adjustable parameter. It is seen that the observed rotation curves reflect the observed structure in gas surface density distribution often in detail. This fact is difficult to comprehend in the context of the dark matter paradigm where the dark halo completely dominates the gravitational potential in the low surface density systems; however it is expected result in the context of modified Newtonian dynamics (MOND) in which the baryonic matter is the only component. With MOND the calculated rotation curves are effectively parameter-free predictions.

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.

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.

The gravitational-wave event GW170817 and the near-simultaneous corresponding gamma-ray burst (GRB 170817 A) falsify modified gravity theories in which the gravitational geometry differs non-conformally from physical geometry. Thus, the observations of this event definitively rule out theories, such as TeVeS, a suggested relativistic extension of Milgrom's modified Newtonian dynamics (MOND), that predict a significantly different Shapiro delay for electromagnetic and gravitational radiation. While not falsifying MOND per se, GW170817 severely constrains relativistic extensions of MOND to theories that do not rely on additional matter-coupling fields but rather upon modified field equations for one universal gravitational and physical metric. Here I mention a simple preferred-frame theory as an example.

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.

Galaxy metallicity scaling relations provide a powerful tool for understanding galaxy evolution, but obtaining unbiased global galaxy gas-phase oxygen abundances requires proper treatment of the various line-emitting sources within spectroscopic apertures. We present a model framework that treats galaxies as ensembles of HII and diffuse ionized gas (DIG) regions of varying metallicities. These models are based upon empirical relations between line ratios and electron temperature for HII regions, and DIG strong-line ratio relations from SDSS-IV MaNGA IFU data. Flux-weighting effects and DIG contamination can significantly affect properties inferred from global galaxy spectra, biasing metallicity estimates by more than 0.3 dex in some cases. We use observationally-motivated inputs to construct a model matched to typical local star-forming galaxies, and quantify the biases in strong-line ratios, electron temperatures, and direct-method metallicities as inferred from global galaxy spectra relative to the median values of the HII region distributions in each galaxy. We also provide a generalized set of models that can be applied to individual galaxies or galaxy samples in atypical regions of parameter space. We use these models to correct for the effects of flux-weighting and DIG contamination in the local direct-method mass-metallicity and fundamental metallicity relations, and in the mass-metallicity relation based on strong-line metallicities. Future photoionization models of galaxy line emission need to include DIG emission and represent galaxies as ensembles of emitting regions with varying metallicity, instead of as single HII regions with effective properties, in order to obtain unbiased estimates of key underlying physical properties.

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 results on the variation of 7.7 micron Polycyclic Aromatic Hydrocarbon (PAH) emission in galaxies spanning a wide range in metallicity at z ~ 2. For this analysis, we use rest-frame optical spectra of 476 galaxies at 1.37 < z < 2.61 from the MOSFIRE Deep Evolution Field (MOSDEF) survey to infer metallicities and ionization states. Spitzer/MIPS 24 micron and Herschel/PACS 100 and 160 micron observations are used to derive rest-frame 7.7 micron luminosities (L(7.7)) and total IR luminosities (L(IR)), respectively. We find significant trends between the ratio of L(7.7) to L(IR) (and to dust-corrected SFR) and both metallicity and [OIII]/[OII] (O32) emission-line ratio. The latter is an empirical proxy for the ionization parameter. These trends indicate a paucity of PAH emission in low metallicity environments with harder and more intense radiation fields. Additionally, L(7.7)/L(IR) is significantly lower in the youngest quartile of our sample (ages of 500 Myr) compared to older galaxies, which may be a result of the delayed production of PAHs by AGB stars. The relative strength of L(7.7) to L(IR) is also lower by a factor of ~ 2 for galaxies with masses $M_* < 10^{10}M_{\odot}$, compared to the more massive ones. We demonstrate that commonly-used conversions of L(7.7) (or 24 micron flux density; f(24)) to L(IR) underestimate the IR luminosity by more than a factor of 2 at $M_*$ ~ $10^{9.6-10.0} M_{\odot}$. We adopt a mass-dependent conversion of L(7.7) to L(IR) with L(7.7)/L(IR)= 0.09 and 0.22 for $M_* < 10^{10}$ and $> 10^{10} M_{\odot}$, respectively. Based on the new scaling, the SFR-$M_*$ relation has a shallower slope than previously derived. Our results also suggest a higher IR luminosity density at z ~ 2 than previously measured, corresponding to a ~ 30% increase in the SFR density.

We highlight phenomenological aspects of Verlinde's recent proposal to account for the mass anomalies in galactic systems without dark matter -- in particular in their relation to MOND. Welcome addition to the MOND lore as it is, this approach have reproduced, so far, only a small fraction of MOND phenomenology, and is still rather tentative, both in its theoretical foundations and in its phenomenology. What Verlinde has extracted from this approach, so far, is a formula -- of rather limited applicability, and with no road to generalization in sight -- for the effective gravitational field of a spherical, isolated, static baryonic system. This formula cannot be used to calculate the gravitational field of disk galaxies, with their rich MOND phenomenology. Notably, it cannot predict their rotation curves, except asymptotically. It does not apply to the few-, or many-body problem; so, it cannot give, e.g., the two-body force between two galaxies, or be used to conduct N-body calculations of galaxy formation, evolution, and interactions. The formula cannot be applied to the internal dynamics of a system embedded in an external field, where MOND predicts important consequences. etc. MOND is backed by full-fledged, Lagrangian theories that can be, and are, routinely applied to all the above phenomena, and more. Verlinde's formula, as it now stands, strongly conflicts with solar-system and possibly earth-surface constraints, and cannot fully account for the mass anomalies in the cores of galaxy clusters (a standing conundrum in MOND). The recent weak-lensing test of the formula is, in fact, testing a cornerstone prediction of MOND, one that the formula does reproduce, and which has been tested before in the very same way.

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 present measurements of the electron-temperature based oxygen abundance for a highly star-forming galaxy at z=3.08, COSMOS-1908. This is the highest redshift at which [OIII]$\lambda$4363 has been detected, and the first time that this line has been measured at z>2. We estimate an oxygen abundance of 12+log(O/H)$=8.00^{+0.13}_{-0.14}$. This galaxy is a low-mass ($10^{9.3}$ M$_{\odot}$), highly star-forming ($\sim50$ M$_{\odot}$ yr$^{-1}$) system that hosts a young stellar population ($\sim160$ Myr). We investigate the physical conditions of the ionized gas in COSMOS-1908 and find that this galaxy has a high ionization parameter, little nebular reddening ($E(B-V)_{\rm gas}<0.14$), and a high electron density ($n_e\sim500$ cm$^{-3}$). We compare the ratios of strong oxygen, neon, and hydrogen lines to the direct-method oxygen abundance for COSMOS-1908 and additional star-forming galaxies at z=0-1.8 with [OIII]$\lambda$4363 measurements, and show that galaxies at z$\sim$1-3 follow the same strong-line correlations as galaxies in the local universe. This agreement suggests that the relationship between ionization parameter and O/H is similar for z$\sim$0 and high-redshift galaxies. These results imply that metallicity calibrations based on lines of oxygen, neon, and hydrogen do not strongly evolve with redshift and can reliably estimate abundances out to z$\sim$3, paving the way for robust measurements of the evolution of the mass-metallicity relation to high redshift.

We present the first direct comparison between Balmer line and panchromatic SED-based SFRs for z~2 galaxies. For this comparison we used 17 star-forming galaxies selected from the MOSFIRE Deep Evolution Field (MOSDEF) survey, with $3\sigma$ detections for H$\alpha$ and at least two IR bands (Spitzer/MIPS 24$\mu$m and Herschel/PACS 100 and 160$\mu$m, and in some cases Herschel/SPIRE 250, 350, and 500$\mu$m). The galaxies have total IR (8-1000$\mu$m) luminosities of $\sim10^{11.4}-10^{12.4}\,\textrm{L}_\odot$ and star-formation rates (SFRs) of $\sim30-250\,\textrm{M}_\odot\,\mathrm{yr^{-1}}$. We fit the UV-to-far-IR SEDs with flexible stellar population synthesis (FSPS) models - which include both stellar and dust emission - and compare the inferred SFRs with the SFR(H$\alpha$,H$\beta$) values corrected for dust attenuation using Balmer decrements. The two SFRs agree with a scatter of 0.17 dex. Our results imply that the Balmer decrement accurately predicts the obscuration of the nebular lines and can be used to robustly calculate SFRs for star-forming galaxies at z~2 with SFRs up to $\sim200\,\textrm{M}_\odot\,\mathrm{yr^{-1}}$. We also use our data to assess SFR indicators based on modeling the UV-to-mid-IR SEDs or by adding SFR(UV) and SFR(IR), for which the latter is based on the mid-IR only or on the full IR SED. All these SFRs show a poorer agreement with SFR(H$\alpha$,H$\beta$) and in some cases large systematic biases are observed. Finally, we show that the SFR and dust attenuation derived from the UV-to-near-IR SED alone are unbiased when assuming a delayed exponentially declining star-formation history.

We present H$\alpha$ gas kinematics for 178 star-forming galaxies at z~2 from the MOSFIRE Deep Evolution Field survey. We have developed models to interpret the kinematic measurements from fixed-angle multi-object spectroscopy, using structural parameters derived from CANDELS HST/F160W imaging. For 35 galaxies we measure resolved rotation with a median $(V/\sigma_{V,0})_{R_E}=2.1$. We derive dynamical masses from the kinematics and sizes and compare them to baryonic masses, with gas masses estimated from dust-corrected H$\alpha$ star formation rates (SFRs) and the Kennicutt-Schmidt relation. When assuming that galaxies with and without observed rotation have the same median $(V/\sigma_{V,0})_{R_E}$, we find good agreement between the dynamical and baryonic masses, with a scatter of $\sigma_{RMS}=0.34$dex and a median offset of $\Delta\log_{10}M=0.04$dex. This comparison implies a low dark matter fraction (8% within an effective radius) for a Chabrier initial mass function (IMF), and disfavors a Salpeter IMF. Moreover, the requirement that $M_{dyn}/M_{baryon}$ should be independent of inclination yields a median value of $(V/\sigma_{V,0})_{R_E}=2.1$ for galaxies without observed rotation. If instead we treat the galaxies without detected rotation as early-type galaxies, the masses are also in reasonable agreement ($\Delta\log_{10}M=-0.07$dex, $\sigma_{RMS}=0.37$dex). The inclusion of gas masses is critical in this comparison; if gas masses are excluded there is an increasing trend of $M_{dyn}/M_{*}$ with higher specific SFR (SSFR). Furthermore, we find indications that $V/\sigma$ decreases with increasing H$\alpha$ SSFR for our full sample, which may reflect disk settling. We also study the Tully-Fisher relation and find that at fixed stellar mass $S_{0.5}=(0.5V_{2.2}^2+\sigma_{V,0}^2)^{1/2}$ was higher at earlier times. At fixed baryonic mass, we observe the opposite trend. [abridged]

We present results on the star-formation rate (SFR) versus stellar mass ($M_*$) relation (i.e., the "main sequence") among star-forming galaxies at 1.37$\leq$$z$$\leq$2.61 using the MOSFIRE Deep Evolution Field (MOSDEF) survey. Based on a sample of 261 galaxies with H$\alpha$ and H$\beta$ spectroscopy, we have estimated robust dust-corrected instantaneous SFRs over a large range in $M_*$ ($\sim10^{9.5}-10^{11.5}M_\odot$). We find a correlation between log(SFR(H$\alpha$)) and log($M_*$) with a slope of 0.65$\pm$0.08 (0.58$\pm$0.10) at 1.4<z<2.6 (2.1<z<2.6). We find that different assumptions for the dust correction, such as using the color-excess of the stellar continuum to correct the nebular lines, sample selection biases against red star-forming galaxies, and not accounting for Balmer absorption can yield steeper slopes of the log(SFR)-log($M_*$) relation. Our sample is immune from these biases as it is rest-frame optically selected, H$\alpha$ and H$\beta$ are corrected for Balmer absorption, and the H$\alpha$ luminosity is dust-corrected using the nebular color-excess computed from the Balmer decrement. The scatter of the log(SFR(H$\alpha$))-log($M_*$) relation, after accounting for the measurement uncertainties, is 0.31 dex at 2.1<z<2.6, which is 0.05 dex larger than the scatter in log(SFR(UV))-log($M_*$). Based on comparisons to a simulated SFR-$M_*$ relation with some intrinsic scatter, we argue that in the absence of direct measurements of galaxy-to-galaxy variations in the attenuation/extinction curves and the IMF, one cannot use the difference in the scatter of the SFR(H$\alpha$)- and SFR(UV)-$M_*$ relations to constrain the stochasticity of star formation in high-redshift galaxies.

Using observations from the MOSFIRE Deep Evolution Field (MOSDEF) survey, we investigate the physical conditions of star-forming regions in $z\sim2.3$ galaxies, specifically the electron density and ionization state. From measurements of the [O II]$\lambda\lambda$3726,3729 and [S II]$\lambda\lambda$6716,6731 doublets, we find a median electron density of $\sim250$ cm$^{-3}$ at $z\sim2.3$, an increase of an order of magnitude compared to measurements of galaxies at $z\sim0$. While $z\sim2.3$ galaxies are offset towards significantly higher O$_{32}$ values relative to local galaxies at fixed stellar mass, we find that the high-redshift sample follows a similar distribution to the low-metallicity tail of the local distribution in the O$_{32}$ vs. R$_{23}$ and O3N2 diagrams. Based on these results, we propose that $z\sim2.3$ star-forming galaxies have the same ionization parameter as local galaxies at fixed metallicity. In combination with simple photoionization models, the position of local and $z\sim2.3$ galaxies in excitation diagrams suggests that there is no significant change in the hardness of the ionizing spectrum at fixed metallicity from $z\sim0$ to $z\sim2.3$. We find that $z\sim2.3$ galaxies show no offset compared to low-metallicity local galaxies in emission line ratio diagrams involving only lines of hydrogen, oxygen, and sulfur, but show a systematic offset in diagrams involving [N II]$\lambda$6584. We conclude that the offset of $z\sim2.3$ galaxies from the local star-forming sequence in the [N II] BPT diagram is primarily driven by elevated N/O at fixed O/H compared to local galaxies. These results suggest that the local gas-phase and stellar metallicity sets the ionization state of star-forming regions at $z\sim0$ and $z\sim2$.

In this paper we present the MOSFIRE Deep Evolution Field (MOSDEF) survey. The MOSDEF survey aims to obtain moderate-resolution (R=3000-3650) rest-frame optical spectra (~3700-7000 Angstrom) for ~1500 galaxies at 1.37<z<3.80 in three well-studied CANDELS fields: AEGIS, COSMOS, and GOODS-N. Targets are selected in three redshift intervals: 1.37<z<1.70, 2.09<z<2.61, and 2.95<z<3.80, down to fixed H_AB (F160W) magnitudes of 24.0, 24.5 and 25.0, respectively, using the photometric and spectroscopic catalogs from the 3D-HST survey. We target both strong nebular emission lines (e.g., [OII], Hbeta, [OIII], 5008, Halpha, [NII], and [SII]) and stellar continuum and absorption features (e.g., Balmer lines, Ca-II H and K, Mgb, 4000 Angstrom break). Here we present an overview of our survey, the observational strategy, the data reduction and analysis, and the sample characteristics based on spectra obtained during the first 24 nights. To date, we have completed 21 masks, obtaining spectra for 591 galaxies. For ~80% of the targets we derive a robust redshift from either emission or absorption lines. In addition, we confirm 55 additional galaxies, which were serendipitously detected. The MOSDEF galaxy sample includes unobscured star-forming, dusty star-forming, and quiescent galaxies and spans a wide range in stellar mass (~10^9-10^11.5 Msol) and star formation rate (~10^0-10^3 Msol/yr). The spectroscopically confirmed sample is roughly representative of an H-band limited galaxy sample at these redshifts. With its large sample size, broad diversity in galaxy properties, and wealth of available ancillary data, MOSDEF will transform our understanding of the stellar, gaseous, metal, dust, and black hole content of galaxies during the time when the universe was most active.

We present results on the dust attenuation curve of z~2 galaxies using early observations from the MOSFIRE Deep Evolution Field (MOSDEF) survey. Our sample consists of 224 star-forming galaxies with nebular spectroscopic redshifts in the range z= 1.36-2.59 and high S/N measurements of, or upper limits on, the H-alpha and H-beta emission lines obtained with Keck/MOSFIRE. We construct composite SEDs of galaxies in bins of specific SFR and Balmer optical depth in order to directly constrain the dust attenuation curve from the UV through near-IR for typical star-forming galaxies at high redshift. Our results imply an attenuation curve that is very similar to the SMC extinction curve at wavelengths redward of 2500 Angstroms. At shorter wavelengths, the shape of the curve is identical to that of the Calzetti relation, but with a lower normalization (R_V). Hence, the new attenuation curve results in SFRs that are ~20% lower, and log stellar masses that are 0.16 dex lower, than those obtained with the Calzetti attenuation curve. Moreover, we find that the difference in the reddening---and the total attenuation---of the ionized gas and stellar continuum correlates strongly with SFR, such that for dust-corrected SFRs larger than 20 Msun/yr assuming a Chabrier IMF, the nebular emission lines suffer an increasing degree of obscuration relative to the continuum. A simple model that can account for these trends is one in which the UV through optical stellar continuum is dominated by a population of less reddened stars, while the nebular line and bolometric luminosities become increasingly dominated by dustier stellar populations for galaxies with large SFRs, as a result of the increased dust enrichment that accompanies such galaxies. Consequently, UV- and SED-based SFRs may underestimate the total SFR at even modest levels of ~20 Msun/yr. [Abridged]

We present results on the excitation properties of z~2.3 galaxies using early observations from the MOSFIRE Deep Evolution Field (MOSDEF) Survey. With its coverage of the full suite of strong rest-frame optical emission lines, MOSDEF provides an unprecedented view of the rest-frame optical spectra of a representative sample of distant star-forming galaxies. We investigate the locations of z~2.3 MOSDEF galaxies in multiple emission-line diagnostic diagrams. These include the [OIII]/Hb vs. [NII]/Ha and [OIII]/Hb vs. [SII]/Ha "BPT" diagrams, as well as the O_32 vs. R_23 excitation diagram. We recover the well-known offset in the star-forming sequence of high-redshift galaxies in the [OIII]/Hb vs. [NII]/Ha BPT diagram relative to SDSS star-forming galaxies. However, the shift for our rest-frame optically selected sample is less significant than for rest-frame-UV selected and emission-line selected galaxies at z~2. Furthermore, we find that the offset is mass-dependent, only appearing within the low-mass half of the z~2.3 MOSDEF sample, where galaxies are shifted towards higher [NII]/Ha at fixed [OIII]/Hb. Within the [OIII]/Hb vs. [SII]/Ha and O_32 vs. R_23 diagrams, we find that z~2.3 galaxies are distributed like local ones, and therefore attribute the shift in the [OIII]/Hb vs. [NII]/Ha BPT diagram to elevated N/O abundance ratios among lower-mass (M_*<10^10 M_sun) high-redshift galaxies. The variation in N/O ratios calls into question the use at high redshift of oxygen abundance indicators based on nitrogen lines, but the apparent invariance with redshift of the excitation sequence in the O_32 vs. R_23 diagram paves the way for using the combination of O_32 and R_23 as an unbiased metallicity indicator over a wide range in redshift. This indicator will allow for an accurate characterization of the shape and normalization of the mass-metallicity relationship over more than 10 Gyr.

We present results on the z~2.3 mass-metallicity relation (MZR) using early observations from the MOSFIRE Deep Evolution Field (MOSDEF) survey. We use an initial sample of 87 star-forming galaxies with spectroscopic coverage of H\beta, [OIII]\lambda 5007, H\alpha, and [NII]\lambda 6584 rest-frame optical emission lines, and estimate the gas-phase oxygen abundance based on the N2 and O3N2 strong-line indicators. We find a positive correlation between stellar mass and metallicity among individual z~2.3 galaxies using both the N2 and O3N2 indicators. We also measure the emission-line ratios and corresponding oxygen abundances for composite spectra in bins of stellar mass. Among composite spectra, we find a monotonic increase in metallicity with increasing stellar mass, offset ~0.15-0.3 dex below the local MZR. When the sample is divided at the median star-formation rate (SFR), we do not observe significant SFR dependence of the z~2.3 MZR among either individual galaxies or composite spectra. We furthermore find that z~2.3 galaxies have metallicities ~0.1 dex lower at a given stellar mass and SFR than is observed locally. This offset suggests that high-redshift galaxies do not fall on the local "fundamental metallicity relation" among stellar mass, metallicity, and SFR, and may provide evidence of a phase of galaxy growth in which the gas reservoir is built up due to inflow rates that are higher than star-formation and outflow rates. However, robust conclusions regarding the gas-phase oxygen abundances of high-redshift galaxies await a systematic reappraisal of the application of locally calibrated metallicity indicators at high redshift.

I review the history and development of Modified Newtonian Dynamics (MOND) beginning with the phenomenological basis as it existed in the early 1980s. I consider Milgrom's papers of 1983 introducing the idea and its consequences for galaxies and galaxy groups, as well as the initial reactions, both negative and positive. The early criticisms were primarily on matters of principle, such as the absence of conservation laws and perceived cosmological problems; an important step in addressing these issues was the development of the Lagrangian-based non-relativistic theory of Bekenstein and Milgrom. This theory led to the development of a tentative relativistic theory that formed the basis for later multi-field theories of gravity. On an empirical level the predictive success of the idea with respect to the phenomenology of galaxies presents considerable challenges for cold dark matter. For MOND the essential challenge remains the absence of a generally accepted theoretical underpinning of the idea and, thus, cosmological predictions. I briefly review recent progress in this direction. Finally I discuss the role and sociology of unconventional ideas in astronomy in the presence of a strongly entrenched standard paradigm.

I show that the lensing masses of the SLACS sample of strong gravitational lenses are consistent with the stellar masses determined from population synthesis models using the Salpeter IMF. This is true in the context of both General Relativity and modified Newtonian dynamics, and is in agreement with the expectation of MOND that there should be little classical discrepancy within the high surface brightness regions probed by strong gravitational lensing. There is also dynamical evidence from this sample supporting the claim that the mass-to-light ratio of the stellar component increases with the velocity dispersion.

The LUX experimental group has just announced the most stringent upper limits so far obtained on the cross section of WIMP-nucleon elastic scattering [1]. This result is a factor of two to five below the previous best upper limit [2] and effectively rules out earlier suggestions of low mass WIMP detection signals. The experimental expertise exhibited by this group is extremely impressive, but the fact of continued negative results raises the more basic question of whether or not this is the right approach to solving the dark matter problem. Here I comment upon this question, using as a basis the final chapter of my book on dark matter [3], somewhat revised and extended. I muse on dark matter and its alternative, modified Newtonian dynamics, or MOND.

We have analyzed OGLE-III photometry for stars in the SMC cluster NGC 299. Two eclipsing binaries and one intrinsic variable (most likely a Be star) are identified. Unlike other young SMC clusters, no low-amplitude pulsating variables are present.

I argue that, despite repeated claims of Ibata et al., the globular cluster NGC 2419 does not pose a problem for modified Newtonian dynamics (MOND). I present a new polytropic model with a running polytropic index. This model provides an improved representation of the radial distribution of surface brightness while maintaining a reasonable fit to the velocity dispersion profile. Although it may be argued that the differences with these observations remain large compared to the reported random errors, there are several undetectable systematic effects which render a formal likelihood analysis irrelevant. I comment generally upon these effects and upon the intrinsic limitations of pressure supported objects as tests of gravity.

Wojtak, Hansen and Hjorth have recently claimed to confirm general relativity and to rule out the tensor-vector-scalar (TeVeS) gravitational theory based on an analysis of the gravitational redshifts of galaxies in 7800 clusters. But their ubiquitous modeling of the sources of cluster gravitational fields with Navarro-Frenk-White mass profiles is neither empirically justified out to the necessary radii in clusters, nor germane in the case of TeVeS. Using MONDian isothermal sphere models consistently constructed within MOND (equivalent to TeVeS models), we can fit the determined redshifts no worse than does general relativity with dark halos. Wojtak, Hansen and Hjorth's work is further marred by confusion between the primitive mu-function of TeVeS and the MOND interpolation function.

Ho\v{r}ava gravity is a attempt to construct a renormalizable theory of gravity by breaking the Lorentz Invariance of the gravitational action at high energies. The underlying principle is that Lorentz Invariance is an approximate symmetry and its violation by gravitational phenomena is somehow hidden to present limits of observational precision. Here I point out that a simple modification of the low energy limit of Ho\v{r}ava gravity in its non-projectable form can effectively camouflage the presence of a preferred frame in regions where the Newtonian gravitational field gradient is higher than $cH_0$; this modification results in the phenomenology of MOND at lower accelerations.

I show that, in the context of MOND, non-isothermal models, approximated by high order polytropic spheres, are consistent with the observations of the radial distribution of the line-of-sight velocity dispersion in the distant globular cluster, NGC 2419. This calls into question the claim by Ibata et al. that the object constitutes a severe challenge for MOND. In general, the existence and properties of globular clusters are more problematic for LCDM than for MOND.

Dwarf and low surface brightness galaxies are ideal objects to test modified Newtonian dynamics (MOND), because in most of these galaxies the accelerations fall below the threshold below where MOND supposedly applies. We have selected from the literature a sample of 27 dwarf and low surface brightness galaxies. MOND is successful in explaining the general shape of the observed rotation curves for roughly three quarters of the galaxies in the sample presented here. However, for the remaining quarter, MOND does not adequately explain the observed rotation curves. Considering the uncertainties in distances and inclinations for the galaxies in our sample, a small fraction of poor MOND predictions is expected and is not necessarily a problem for MOND. We have also made fits taking the MOND acceleration constant, a_0, as a free parameter in order to identify any systematic trends. We find that there appears to be a correlation between central surface brightness and the best-fit value of a_0, in the sense that lower surface brightness galaxies tend to have lower a_0. However, this correlation depends strongly on a small number of galaxies whose rotation curves might be uncertain due to either bars or warps. Without these galaxies, there is less evidence of a trend, but the average value we find for a_0 ~ 0.7*10^-8 cm s^-2 is somewhat lower than derived from previous studies. Such lower fitted values of a_0 could occur if external gravitational fields are important.

In the context of modified Newtonian dynamics, the fundamental plane, as the observational signature of the Newtonian virial theorem, is defined by high surface brightness objects that deviate from being purely isothermal: the line-of-sight velocity dispersion should slowly decline with radius as observed in luminous elliptical galaxies. All high surface brightness objects (e.g. globular clusters, ultra-compact dwarfs) will lie, more or less, on the fundamental plane defined by elliptical galaxies, but low surface brightness objects (dwarf spheroidals) would be expected to deviate from this relation. This is borne out by observations. With MOND, the Faber-Jackson relation (the power-law relation between luminosity and velocity dispersion), ranging from globular clusters to clusters of galaxies and including both high and low surface brightness objects, is the more fundamental and universal scaling relation in spite of its larger scatter. Faber-Jackson reflects the presence of an additional dimensional constant (the MOND acceleration) in the structure equation.

MOND-- modified Newtonian dynamics-- may be viewed as an algorithm for calculating the distribution of force in an astronomical object from the observed distribution of baryonic matter. The fact that it works for galaxies is quite problematic for Cold Dark Matter. Moreover, MOND explains or subsumes systematic aspects of galaxy photometry and kinematics-- aspects that CDM does not address or gets wrong. I will present evidence here in support of these assertions and claim that this is effectively a falsification of dark matter that is dynamically important on the scale of galaxies.

Bolton et al. (2007) have derived a mass-based fundamental plane using photometric and spectroscopic observations of 36 strong gravitational lenses. The lensing allows a direct determination of the mass-surface density and so avoids the usual dependence on mass-to-light ratio. We consider this same sample in the context of modified Newtonian dynamics (MOND) and demonstrate that the observed mass-based fundamental plane coincides with the MOND fundamental plane determined previously for a set of high-order polytropic spheres chosen to match the observed range of effective radii and velocity dispersions in elliptical galaxies. Moreover, the observed projected mass within one-half an effective radius is consistent with the mass in visible stars plus a small additional component of ``phantom dark matter'' resulting from the MOND contribution to photon deflection.

MOND predicts that a mass, M, contained within its transition radius r_t=(MG/a0)^{1/2}, may exhibit a feature at about that radius in the form of a shell, or projected ring, in the deduced distribution of its phantom dark matter. This occurs despite the absence of any underlying feature in the true (``baryon'') source distribution itself. The phenomenon is similar to the appearance of an event horizon and other unusual physics ``in the middle of nothing'' near the transition radius of General Relativity $MG/c^2$. We consider the possibility that this pure MOND phenomenon is in the basis of the recent finding of such a ring in the galaxy cluster Cl 0024+17 by Jee et al. We find that the parameters of the observed ring can be naturally explained in this way; this feature may therefore turn out to be a direct evidence for MOND. We study this phenomenon in simple, axisymmetric configurations aligned with the line of sight: spherical masses, a dumbbell of spherical masses, and an elongated, thin structure. The properties of the apparent ring: its radius, surface density, and contrast, depend on the form of the MOND interpolating function and on the exact three dimensional distribution of the sources (the thin-lens approximation is quite invalid in MOND). We also comment on the possible appearance of orphan features, marking the Newtonian-to-MOND transition, in high surface brightness galaxies. In particular, we find that previously unexplained structure in the rotation curves of some galaxies may be evidence for such features.

Beginning with a simple model for the growth of structure, I consider the dissipationless evolution of a MOND-dominated region in an expanding Universe by means of a spherically symmetric N-body code. I demonstrate that the final virialized objects resemble elliptical galaxies with well-defined relationships between the mass, radius, and velocity dispersion. These calculations suggest that, in the context of MOND, massive elliptical galaxies may be formed early (z > 10) as a result of monolithic dissipationless collapse. Then I reconsider the classic argument that a galaxy of stars results from cooling and fragmentation of a gas cloud on a time scale shorter than that of dynamical collapse. Qualitatively, the results are similar to that of the traditional picture; moreover, the existence, in MOND, of a density-temperature relation for virialized, near isothermal objects as well as a mass-temperature relation implies that there is a definite limit to the mass of a gas cloud where this condition can be met-- an upper limit corresponding to that of presently observed massive galaxies.

The dynamical mass of clusters of galaxies, calculated in terms of modified Newtonian dynamics, is a factor of two or three times smaller than the Newtonian dynamical mass but remains significantly larger than the observed baryonic mass in the form of hot gas and stars in galaxies. Here I consider further the suggestion that the undetected matter might be in the form of cosmological neutrinos with mass on the order of 2 eV. If the neutrinos and baryons have comparable velocity dispersions and if the two components maintain their cosmological density ratio, then the electron density in the cores of clusters should be proportional to temperature to the 3/2 power, as appears to be true in non-cooling flow clusters. This is equivalent to the ``entropy floor'' proposed to explain the steepness of the observed luminosity-temperature relation, but here preheating of the medium is not required. Two fluid (neutrino-baryon) hydrostatic models of clusters, in the context of MOND, reproduce the observed luminosity-temperature relation of clusters. If the beta law is imposed the gas density distribution, then the self-consistent models predict the general form of the observed temperature profile in both cooling and non-cooling flow clusters.

We extend the MOND analysis to a sample of 17 high surface brightness, early-type disc galaxies with rotation curves derived from a combination of 21cm HI line observations and optical spectroscopic data. A number of these galaxies have asymptotic rotation velocities between 250 and 350 km/s making them among the most massive systems (in terms of baryonic mass) considered in the context of MOND. We find that the general MOND prediction for such galaxies -- a rotation curve which gradually declines to the asymptotic value -- is confirmed, and in most cases the MOND rotation curve, determined from the mean radial light and gas distribution, agrees in detail with the observed rotation curve. In the few cases where MOND appears not to work well, the discrepancies can generally be understood in terms of various observational errors -- such as incorrect orientation angles and/or distances -- or of unmodelled physical effects -- such as non-circular motions. The implied mass-to-light ratios for the stellar disc and bulge constrain the MOND interpolating function; the form recently suggested by Zhao & Famaey (2005) yields more sensible values than the one traditionally used in MOND determinations of galaxy rotation curves.

Below a threshold energy, gas in the constant density core of a triaxial galaxy can find no simple non-intersecting periodic orbit to act as an attractor for its trajectory (El-Zant et al. 2003). If a disc of gas arriving from further out in the galaxy dissipates sufficient energy to fall below this threshold, it will thereafter collapse into the very centre. Such a mechanism may be relevant to the early growth of super-massive black holes at the Eddington limit and the appearance of the quasar phenomenon at high redshift. This process is self-limiting in the sense that, when the black hole mass has grown to a significant fraction of the core mass, simple angular momentum conserving orbits are restored and accretion reverts to the slow viscous mode. The mechanism depends upon the pre-existence of constant density cores in triaxial spheroidal galaxies.

We present MOND analysis for several of the lowest mass disc galaxies currently amenable to such analysis--with (baryonic) masses below 4x10^8 solar masses. The agreement is good, extending the validity of MOND and its predicted mass velocity relation, to such low masses.

Any viable theory of modified Newtonian dynamcs (MOND) as modified gravity is likely to require fields in addition to the usual tensor field of General Relativity. For such theories the MOND phenomenology emerges from an effective fifth force probably associated with a scalar field. Here I consider the constraints imposed upon such theories by solar system phenomenology, primarily by the absence of significant deviations from inverse square attraction in the inner solar system as well as detectable local preferred frame effects. The current examples of multi-field theories can be constructed to satisfy these constraints and such theories lead inevitably to an anomalous non inverse-square force in the outer solar system.

On an empirical level, the most successful alternative to dark matter in bound gravitational systems is the modified Newtonian dynamics, or MOND, proposed by Milgrom. Here I discuss the attempts to formulate MOND as a modification of General Relativity. I begin with a summary of the phenomenological successes of MOND and then discuss the various covariant theories that have been proposed as a basis for the idea. I show why these proposals have led inevitably to a multi-field theory. I describe in some detail TeVeS, the tensor-vector-scalar theory proposed by Bekenstein, and discuss its successes and shortcomings. This lecture is primarily pedagogical and directed to those with some, but not a deep, background in General Relativity

I review various ideas on MOND cosmology and structure formation beginning with non-relativistic models in analogy with Newtonian cosmology. I discuss relativistic MOND cosmology in the context of Bekenstein's theory and propose an alternative biscalar effective theory of MOND in which the acceleration parameter is identified with the cosmic time derivative of a matter coupling scalar field. Cosmic CDM appears in this theory as scalar field oscillations of the auxiliary "coupling strength" field.

We first review the nonrelativistic lagrangian theory as a framework for the MOND equation. Obstructions to a relativistic version of it are discussed leading up to TeVeS, a relativistic tensor-vector-scalar field theory which displays both MOND and Newtonian limits. The whys for its particular structure are discussed and its achievements so far are summarized.

I describe a tensor-vector-scalar theory that reconciles the galaxy scale success of modified Newtonian dynamics (MOND) with the cosmological scale evidence for CDM. The theory provides a cosmological basis for MOND in the sense that the predicted phenomenology only arises in a cosmological background. The theory contains an evolving effective potential, and scalar field oscillations in this potential comprise the cold dark matter; the de Broglie wavelength of these soft bosons, however, is sufficiently large that they cannot accumulate in galaxies. The theory predicts, inevitably, a constant anomalous acceleration in the outer solar system which, depending upon the choice of parameters, can be consistent with that detected by the Pioneer spacecrafts.

We examine two corollaries of MOND pertaining to the properties of the equivalent dark-matter halo. MOND predicts for pure exponential discs a tight relation involving the halo and disc scale lengths and the mean acceleration in the disc, which we find to test favorably against the Verheijen sample of Ursa Major galaxies. A correlation between halo and disc length scales is also apparent when the "maximum disc" contribution is assumed, but we demonstrate that this follows from the more general MOND prediction. The second MOND prediction involves the existence of a maximum halo acceleration, which also tests favorably against the Ursa Major sample for different assumptions on the disc contribution.

I discuss the classical cosmological tests-- angular size-redshift, flux-redshift, and galaxy number counts-- in the light of the cosmology prescribed by the intepretation of CMB anisotropies. The discussion is somewhat of a primer for physicists, with emphasis upon the possible systematic uncertainties in the observations and their interpretation. Given the curious composition of the Universe inherent in the emerging cosmological model, I stress the value of searching for inconsistencies rather than concordance, and suggest that the prevailing mood of triumphalism in cosmology is premature.

The recent findings of Romanowsky et al., of an ``unexpectedly'' small mass discrepancy within 5 effective radii in several elliptical galaxies, are not surprising in the context of MOND. As we show here, they are, in fact, in full concordance with its predictions. One is dealing with high-surface-density galaxies with mean accelerations rather larger than the acceleration constant of MOND. These findings continue, and are now the extreme examples of, the trend predicted by MOND: the mass discrepancy sets in at larger and larger scaled radii in galaxies with larger and larger mean surface densities; or, equivalently, mean accelerations.

X-ray emitting clusters of galaxies are considered in the context of modified Newtonian dynamics (MOND). I show that self-gravitating isothermal gas spheres are not good representations of rich clusters with respect to the radial gas density distribution as indicated by the X-ray surface brightness. Pure gas spheres with a density distribution described by a ``beta'' model, as observed, also fail because, with MOND, these objects are far from isothermal. These problems may be resolved by adding an additional dark mass component in the central regions; a constant density sphere contained within two core radii and having a total mass of one to two times that in the gas. With this additional component, the observed luminosity-temperature relation for clusters of galaxies is reproduced. When the observed X-ray surface brightness distribution in actual clusters is modeled by such a two-component structure, the typical mass discrepancy is three to four times smaller than with Newtonian dynamics. Thus while MOND significantly reduces the mass of the dark component in clusters it does not remove it completely. I speculate on the nature of the dark component and argue that this is not a fundamental problem for MOND.

Rotation curves for four spiral galaxies with recently determined Cepheid-based distances are reconsidered in terms of modified Newtonian dynamics (MOND). For two of the objects, NGC 2403 and NGC 7331, the rotation curves predicted by MOND are compatible with the observed curves when these galaxies are taken to be at the Cepheid distance. For NGC 3198, the largest distance for which reasonable agreement is obtained is 10% smaller than the Cepheid-based distance; i.e., MOND clearly prefers a smaller distance. This conclusion is unaltered when new near-infrared photometry of NGC 3198 is taken as the tracer of the stellar mass distribution. For the large Sc spiral, NGC 2841, MOND requires a distance which is at least 20% larger than the Cepheid-based distance. However, the discrepancy of the Tully-Fisher and SNIa distances with the Cepheid determination casts some doubt upon the Cepheid method in this case.

Modified Newtonian dynamics (MOND) is an empirically motivated modification of Newtonian gravity or inertia suggested by Milgrom as an alternative to cosmic dark matter. The basic idea is that at accelerations below a0 ~ 10^{-8} cm/s^2 ~ cH0/6 the effective gravitational attraction approaches sqrt{gN*a0} where gN is the usual Newtonian acceleration. This simple algorithm yields flat rotation curves for spiral galaxies and a mass-rotation velocity relation of the form M ~ V^4 that forms the basis for the observed luminosity-rotation velocity relation-- the Tully-Fisher law. We review the phenomenological success of MOND on scales ranging from dwarf spheroidal galaxies to superclusters, and demonstrate that the evidence for dark matter can be equally well interpreted as evidence for MOND. We discuss the possible physical basis for an acceleration-based modification of Newtonian dynamics as well as the extension of MOND to cosmology and structure formation.

Milgrom has proposed that the appearance of discrepancies between the Newtonian dynamical mass and the directly observable mass in astronomical systems could be due to a breakdown of Newtonian dynamics in the limit of low accelerations rather than the presence of unseen matter. Milgrom's hypothesis, modified Newtonian dynamics or MOND, has been remarkably successful in explaining systematic properties of spiral and elliptical galaxies and predicting in detail the observed rotation curves of spiral galaxies with only one additional parameter-- a critical acceleration which is on the order of the cosmologically interesting value of $cH_o$. Here I review the empirical successes of this idea and discuss its possible extention to cosmology and structure formation.

I consider the growth of inhomogeneities in a low-density baryonic, vacuum energy-dominated universe in the context of modified Newtonian dynamics (MOND). I first write down a two-field Langrangian-based theory of MOND (non-relativistic), which embodies several assumptions such as constancy of the MOND acceleration parameter, association of a MOND force with peculiar accelerations only, and the deceleration of the Hubble flow as a background field which influences the dynamics of a finite size region. In the context of this theory, the equation for the evolution of spherically symmetric over-densities is non-linear and implies very rapid growth even in a low-density background, particularly at the epoch when the putative cosmological constant begins to dominate the Hubble expansion. Small comoving scales enter the MOND regime earlier than larger scales and therefore evolve to large over-densities sooner. Taking the initial COBE-normalized power spectrum provided by CMBFAST (Seljak & Zeldarriaga 1996), I find that the final power-spectrum resembles that of the standard LCDM universe and thus retains the empirical successes of that model.

In a significant recent paper A. Edery undertakes a new study of light deflection in generalizations of general relativity (GR). He claims to prove that any metric-based gravitational theory that proposes to explain the flat rotation curves of disk galaxies without postulating dark matter halos must conflict with observations of gravitational lensing by galaxies and clusters of galaxies because any such theory inevitably make a negative contribution to light deflection. Here we show that some of the basic steps of Edery's argument are invalid, and no such general result obtains.

The modified Newtonian dynamics (MOND), suggested by Milgrom as an alternative to dark matter, implies that isothermal spheres with a fixed anisotropy parameter should exhibit a near perfect relation between the mass and the fourth power of the velocity dispersion. This is consistent with the observed Faber-Jackson relation for elliptical galaxies-- a luminosity-velocity dispersion relation with large scatter. However, the observable global properties of elliptical galaxies comprise a three parameter family; they lie on a ``fundamental plane'' in a logarithmic space consisting of central velocity dispersion, effective radius, and luminosity. The scatter perpendicular to this plane is significantly less than that about the Faber-Jackson relation. I show here that, in order to match the observed global properties of elliptical galaxies with MOND, models must deviate from being strictly isothermal and isotropic; such objects can be approximated by high-order polytropic spheres with a radial orbit anisotropy in the outer regions. MOND imposes boundary conditions on the inner Newtonian regions which restrict these models to a dynamical fundamental plane which may differ from that implied by the traditional virial theorem. Scatter about this plane is relatively insensitive to the necessary deviations from homology.

Sticky particle calculations indicate that a coherent structure, a dispersion ring, forms when a cloud on a low angular momentum orbit passes close to the dynamical center of an isothermal sphere containing a central point mass. The cloud is tidally stretched and differentially wrapped, and dissipation in shocks organizes the gas into a precessing off-set elliptical ring which can persist for many rotation periods. The morphology and kinematics of the circumnuclear disk (CND) between 2 and 5 pc and the Northern arm in the inner 1 pc are well-represented by such structures. In the case of the Northern Arm, strong shocks which arise during the formation of the dispersion ring can lead to star formation even in the near tidal field of a massive black hole.

A sample of 197 X-ray emitting clusters of galaxies is considered in the context of Milgrom's modified Newtonian dynamics (MOND). It is shown that the gas mass, extrapolated via an assumed $\beta$ model to a fixed radius of 3 Mpc, is correlated with the gas temperature as predicted by MOND ($M_g \propto T^2$). The observed temperatures are generally consistent with the inferred mass of hot gas; no substantial quantity of additional unseen matter is required in the context of MOND. However, modified dynamics cannot resolve the strong lensing discrepancy in those clusters where this phenomenon occurs. The prediction is that additional baryonic matter may be detected in the central regions of rich clusters.

This is the third in a series of papers in which spiral galaxy rotation curves are considered in the context of Milgrom's modified dynamics (MOND). The present sample of 30 objects is drawn from a complete sample of galaxies in the Ursa Major cluster with photometric data by Tully et al. (1996) and 21 cm line data by Verheijen (1997). The galaxies are roughly all at the same distance (15 to 16 Mpc). The radio observations are made with the Westerbork Synthesis Array which means that the linear resolution of all rotation curves is comparable. The greatest advantage of this sample is the existance of K'-band surface photometry for all galaxies; the near-infrared emission, being relatively free of the effects of dust absorption and less sensitive to recent star formation, is a more precise tracer of the mean radial distribution of the dominant stellar population. The predicted rotation curves are calculated from the K'-band surface photometry and the observed distribution of neutral hydrogen using the simple MOND prescription where the one adjustable parameter is the mass of the stellar disk or the implied mass-to-light ratio. The predicted rotation curves generally agree with the observed curves and the mean M/L in the near-infrared is about 0.9 with a small dispersion. The fitted M/L in the B-band is correlated with B-V color in the sense predicted by population synthesis models. Including earlier work, about 80 galaxy rotation curves are now well-reproduced from the observed distribution of detectable matter using the MOND formula to calculate the gravitational acceleration; this lends considerable observational support to Milgrom's unconventional hypothesis.

It is well-known that the application of Newtonian dynamics to an expanding spherical region leads to the correct relativistic expression (the Friedmann equation) for the evolution of the cosmic scale factor. Here, the cosmological implications of Milgrom's modified Newtonian dynamics (MOND) are considered by means of a similar procedure. Earlier work by Felten demonstrated that in a region dominated by modified dynamics, the expansion cannot be uniform (separations cannot be expressed in terms of a scale factor) and that any such region will eventually re-collapse regardless of the initial expansion velocity and mean density. Here I show that, because of the acceleration threshold for the MOND phenomenology, a region dominated by MOND will have a finite size which, in the early Universe ($z > 3$), is smaller than the horizon scale. Therefore, uniform expansion and homogeneity on horizon scale are consistent with MOND-dominated non-uniform expansion and the development of inhomogeneities on smaller scale. In the radiation-dominated era, the amplitude of MOND-induced inhomogeneities is much smaller than that implied by observations of the cosmic background radiation, and the thermal and dynamical history of the Universe is identical to that of the standard Big Bang. When matter first dominates the energy density of the Universe, the cosmology diverges from that of the standard model. Objects of galaxy mass are the first virialized objects to form (by z=10) and larger structure develops rapidly. At the present epoch, the Universe would be inhomogeneous out to a substantial fraction of the Hubble radius.

On the basis of ``sticky particle'' calculations, it is argued that the gas features observed within 10 pc of the Galactic Centre-- the circumnuclear disk (CND) and the ionized gas filaments-- as well as the newly formed stars in the inner one parsec can be understood in terms of tidal capture and disruption of gas clouds on low angular momentum orbits in a potential containing a point mass. The calculations demonstrate that a dissipative component forms a ``dispersion ring'', an asymmetric elliptical torus precessing counter to the direction of rotation, and that this shape can be maintained for many orbital periods. For a range of plausible initial conditions, such a sturcture can explain the morphology and kinematics of the CND and of the most conspicuous ionized filament. While forming the dispersion ring, a small cloud with low specific angular momentum is drawn into a long filament which repeatedly collides with itself at high velocity. The compression in strong shocks is likely to lead to star formation even in the near tidal field of the point mass. This process may have general relevance to accretion onto massive black holes in normal and active galactic nuclei.

Although the modified dynamics (MOND) proposed by Milgrom successfully accounts for the systematics of galaxy rotation curves and cluster dynamics without invoking dark matter, the idea remains a largely ad hoc modification of Newtonian dynamics with no basis in deeper theory. Non-standard scalar-tensor theories have been suggested as a theoretical basis for MOND; however, any such theory with the usual conformal relation between the Einstein and physical metrics fails to predict the degree of light deflection observed in distant clusters of galaxies. In the present paper, I demonstrate that one can write down a framework for scalar-tensor theories which predict the MOND phenomenology for the low velocity (v<<c) dynamics of galaxies and clusters of galaxies and which are consistent with observations of extragalactic gravitational lenses provided that one drops the requirement of the Lorentz Invariance of gravitational dynamics. This leads to ``preferred frame'' theories characterized by a non-conformal relation between the two metrics. I describe a toy theory in which the local environment (the solar system, binary pulsars) is protected from detectable preferred frame effects by the very same non-standard (aquadratic) scalar Lagrangian which gives rise to the MOND phenomenology. Although this particular theory is also contrived, it represents a limiting case for two-field theories of MOND, and is consistent with a wide range of gravitational phenomena.

A sample of 22 spiral galaxy rotation curves, measured in the 21 cm line of neutral hydrogen, is considered in the context of Milgrom's modified dynamics (MOND). Combined with the previous highly selected sample of Begeman et al. (1990), this comprises the current total sample of galaxies with published (or available) extended rotation curves and photometric observations of the light distribution. This is the observational basis of present quantitative understanding of the discrepancy between the visible mass and classical dynamical mass in galaxies. It is found that the gravitational force calculated from the observed distribution of luminous material and gas using the simple MOND formula can account for the overall shape and amplitude of these 22 rotation curves, and in some cases, the predicted curve agrees with the observed rotation curve in detail. The fitted rotation curves have, in 13 cases, only one free parameter which is the mass-to-light ratio of the luminous disk; in nine cases, there is an additional free parameter which is M/L of a central bulge or light concentration. The values of the global M/L (bulge plus disk) are reasonable and, when the gas mass is also included, show a scatter which is consistent with that in the Tully-Fisher relation. The success of the MOND prescription in predicting the rotation curves in this larger, less stringently selected sample, lends further support to the idea that dynamics or gravity is non-Newtonian in the limit of low accelerations and that it is unnecessary to invoke the presence of large quantities of unseen matter.

We study gravitational lensing by clusters of galaxies in the context of the generic class of unconventional gravity theories of the scalar--tensor type. For positive energy scalar fields with any dynamics, the bending of light by a weakly gravitating system (galaxy or cluster) is always smaller than the bending predicted by general relativity for the mass of visible and hitherto undetected matter (but excluding the scalar field energy). The same conclusion obtains within general relativity if a nonnegligible part of the mass in clusters is in coherent scalar fields, {\it i.e.\/} Higgs fields. Thus the observational claim that clusters of galaxies deflect light much more strongly than would be expected from the observable matter contained by them, cannot be interpreted in terms of some scalar-tensor unconventional gravity theory with no dark matter. If the observations eventually show that the matter distribution inferred via general relativity from the lensing is very much like that determined from the dynamics of test objects, then scalar--tensor unconventional gravity will be irrelevant for understanding the mass discrepancy in clusters. However, even a single system in which the dynamically determined mass significantly exceeds the lensing mass suggested by general relativity, would be very problematic for the dark matter picture, but would be entirely consistent with unconventional scalar--tensor gravity.