Search

The radial gradient of gas-phase metallicity is a powerful probe of the chemical and structural evolution of star-forming galaxies, closely tied to disk formation and gas kinematics in the early universe. We present spatially resolved chemical and dynamical properties for a sample of 26 galaxies at $0.5 \lesssim z \lesssim 1.7$ from the MSA-3D survey. These innovative observations provide 3D spectroscopy of galaxies at a spatial resolution approaching JWST's diffraction limit and a high spectral resolution of $R\simeq2700$. The metallicity gradients measured in our galaxy sample range from $-$0.05 to 0.02 dex~kpc$^{-1}$. Most galaxies exhibit negative or flat radial gradients, indicating lower metallicity in the outskirts or uniform metallicity throughout the entire galaxy. We confirm a tight relationship between stellar mass and metallicity gradient at $z\sim1$ with small intrinsic scatter of 0.02 dex~kpc$^{-1}$. Our results indicate that metallicity gradients become increasingly negative as stellar mass increases, likely because the more massive galaxies tend to be more "disky". This relationship is consistent with the predictions from cosmological hydrodynamic zoom-in simulations with strong stellar feedback. This work presents the effort to harness the multiplexing capability of JWST NIRSpec/MSA in slit-stepping mode to map the chemical and kinematic profiles of high-redshift galaxies in large samples and at high spatial and spectral resolution., Comment: 23 pages, 6 figures, 1 table; submitted to ApJL

We show that a quantum spin system has an exact free-parafermion solution if its frustration graph is an oriented indifference graph. Further, we show that if the frustration graph of a model can be dipath oriented via switching operations, then the model is integrable in the sense that there is a family of commuting independent set charges. Additionally, we establish an efficient algorithm for deciding whether this is possible. Our characterisation extends that given for free-fermion solvability. Finally, we apply our results to solve three qudit spin models., Comment: 22 pages, 2 figures

We use JWST/NIRSpec observations from the Assembly of Ultradeep Rest-optical Observations Revealing Astrophysics (AURORA) survey to constrain the shape of the nebular attenuation curve of a star-forming galaxy at z=4.41, GOODSN-17940. We utilize 11 unblended HI recombination lines to derive the attenuation curve spanning optical to near-infrared wavelengths (3751-9550 \r{A}). We then leverage a high-S/N spectroscopic detection of the rest-frame ultraviolet continuum in combination with rest-UV photometric measurements to constrain the shape of the curve at ultraviolet wavelengths. While this UV constraint is predominantly based on stellar emission, the large measured equivalent widths of H$\alpha$ and H$\beta$ indicate that GOODSN-17940 is dominated by an extremely young stellar population <10 Myr in age such that the UV stellar continuum experiences the same attenuation as the nebular emission. The resulting combined nebular attenuation curve spans 1400-9550 \r{A} and has a shape that deviates significantly from commonly assumed dust curves in high-redshift studies. Relative to the Milky Way, SMC, and Calzetti curves, the new curve has a steeper slope at long wavelengths ($\lambda>5000$ \r{A}) while displaying a similar slope across blue-optical wavelengths ($\lambda=3750-5000$ \r{A}). In the ultraviolet, the new curve is shallower than the SMC and Calzetti curves and displays no significant 2175 \r{A} bump. This work demonstrates that the most commonly assumed dust curves are not appropriate for all high-redshift galaxies. These results highlight the ability to derive nebular attenuation curves for individual high-redshift sources with deep JWST/NIRSpec spectroscopy, thereby improving the accuracy of physical properties inferred from nebular emission lines., Comment: 21 pages, 6 figures, 1 table, submitted to ApJ

The metal content of galaxies is a direct probe of the baryon cycle. A hallmark example is the relationship between a galaxy's stellar mass, star formation rate (SFR), and gas-phase metallicity: the Fundamental Metallicity Relation (FMR). While low-redshift ($z\lesssim4$) observational studies suggest that the FMR is redshift-invariant, recent JWST data indicate deviations from this model. In this study, we utilize the FMR to predict the evolution of the normalisation of the mass-metallicity relation (MZR) using the cosmological simulations Illustris, IllustrisTNG, EAGLE, and SIMBA. Our findings demonstrate that a $z = 0$ calibrated FMR struggles to predict the evolution in the MZR of each simulation. To quantify the divergence of the predictions, we introduce the concepts of a ''static'' FMR, where the role of the SFR in setting the normalization of the MZR does not change with redshift, and a ''dynamic'' FMR, where the role of SFR evolves over time. We find static FMRs in Illustris and SIMBA and dynamic FMRs in IllustrisTNG and EAGLE. We suggest that the differences between these models likely points to the subtle differences in the implementation of the baryon cycle. Moreover, we echo recent JWST results at $z > 4$ by finding significant offsets from the FMR in IllustrisTNG and EAGLE, suggesting that the observed FMR may be dynamic as well. Overall, our findings imply that the current FMR framework neglects important variations in the baryon cycle through cosmic time., Comment: Submitted to MNRAS. 17 pages, 8 Figures. See also paper one of series (arXiv:2403.08856). Comments welcome

We present results on the emission-line properties of z=1.4-7.5 star-forming galaxies in the Assembly of Ultradeep Rest-optical Observations Revealing Astrophysics (AURORA) Cycle 1 JWST/NIRSpec program. Based on its depth, continuous wavelength coverage from 1--5 microns, and medium spectral resolution (R~1000), AURORA includes detections of a large suite of nebular emission lines spanning a broad range in rest wavelength. We investigate the locations of AURORA galaxies in multiple different emission-line diagrams, including traditional "BPT" diagrams of [OIII]/Hbeta vs. [NII]/Halpha, [SII]/Halpha, and [OI]/Halpha, and the "ionization-metallicity" diagram of [OIII]/[OII] (O32) vs. ([OIII]+[OII])/Hbeta (R23). We also consider a bluer rest-frame "ionization-metallicity" diagram introduced recently to characterize z>10 galaxies: [NeIII]/[OII] vs. ([NeIII]+[OII])/Hdelta; as well as longer-wavelength diagnostic diagrams extending into the rest-frame near-IR: [OIII]/Hbeta vs. [SIII]/[SII] (S32); and HeI/Pagamma and [SIII]/Pagamma vs. [FeII]/Pabeta. With a significant boost in signal-to-noise and large, representative samples of individual galaxy detections, the AURORA emission-line diagrams presented here definitively confirm a physical picture in which chemically-young, alpha-enhanced, massive stars photoionize the ISM in distant galaxies with a harder ionizing spectrum at fixed nebular metallicity than in their z~0 counterparts. We also uncover previously unseen evolution prior to z~2 in the [OIII]/Hbeta vs. [NII]/Halpha diagram, which motivates deep NIRSpec observations at even higher redshift. Finally, we present the first statistical sample of rest-frame near-IR emission-line diagnostics in star-forming galaxies at high redshift. In order to truly interpret rest-frame near-IR line ratios including [FeII], we must obtain better constraints on dust depletion in the high-redshift ISM., Comment: 17 pages, 8 figures, submitted to ApJ

We have used public JWST/NIRSpec and JWST/NIRCam observations from the CEERS and JADES surveys in order to analyze the star-forming main sequence (SFMS) over the redshift range $1.4 \leq z < 7$. We calculate the star-formation rates (SFRs) of the galaxy sample using three approaches: Balmer line luminosity, spectral energy distribution (SED) fitting, and UV luminosity. We find a larger degree of scatter about the SFMS using the Balmer-based SFRs compared to the UV-based SFRs. Because these SFR indicators are sensitive to star formation on different time scales, the difference in scatter may be evidence of bursty star-formation histories in the early universe. We additionally compare the H$\alpha$-to-UV luminosity ratio (L(H$\alpha$)/$\nu$L$_{\nu,1600}$) for individual galaxies in the sample and find that 29\%$-$52\% of the ratios across the sample are poorly described by predictions from a smooth star-formation history. Measuring the burstiness of star formation in the early universe has multiple significant implications, such as deriving accurate physical parameters from SED fitting, explaining the evolution of the UV luminosity function, and providing constraints for sub-grid models of feedback in simulations of galaxy formation and evolution., Comment: 21 pages, 7 figures, submitted to ApJ

Accurate chemical compositions of star-forming regions are a critical diagnostic tool to characterize the star formation history and gas flows which regulate galaxy formation. However, the abundance discrepancy factor (ADF) between measurements from the "direct" optical electron temperature ($T_e$) method and from the recombination lines (RL) represents $\sim0.2$ dex systematic uncertainty in oxygen abundance. The degree of uncertainty for other elements is unknown. We conduct a comprehensive analysis of O$^{++}$ and N$^+$ ion abundances using optical and far-infrared spectra of a star-forming region within the nearby dwarf galaxy Haro 3, which exhibits a typical ADF. Assuming homogeneous conditions, the far-IR emission indicates an O abundance which is higher than the $T_e$ method and consistent with the RL value, as would be expected from temperature fluctuations, whereas the N abundance is too large to be explained by temperature fluctuations. Instead a component of highly obscured gas is likely required to explain the high far-IR to optical flux ratios. Accounting for this obscured component reduces both the IR-based metallicities and the inferred magnitude of temperature fluctuations, such that they cannot fully explain the ADF in Haro 3. Additionally, we find potential issues when predicting the RL fluxes from current atomic data. Our findings underscore the critical importance of resolving the cause of abundance discrepancies and understanding the biases between different metallicity methods. This work represents a promising methodology, and we identify further approaches to address the current dominant uncertainties., Comment: 18 pages, 7 figures, and 2 tables. Submitted to ApJ

About 40 transneptunian binaries (TNBs) have fully determined orbits with about 10 others being solved except for breaking the mirror ambiguity. Despite decades of study almost all TNBs have only ever been analyzed with a model that assumes perfect Keplerian motion (e.g., two point masses). In reality, all TNB systems are non-Keplerian due to non-spherical shapes, possible presence of undetected system components, and/or solar perturbations. In this work, we focus on identifying candidates for detectable non-Keplerian motion based on sample of 45 well-characterized binaries. We use MultiMoon, a non-Keplerian Bayesian inference tool, to analyze published relative astrometry allowing for non-spherical shapes of each TNB system's primary. We first reproduce the results of previous Keplerian fitting efforts with MultiMoon, which serves as a comparison for the non-Keplerian fits and confirms that these fits are not biased by the assumption of a Keplerian orbit. We unambiguously detect non-Keplerian motion in 8 TNB systems across a range of primary radii, mutual orbit separations, and system masses. As a proof of concept for non-Keplerian fitting, we perform detailed fits for (66652) Borasisi-Pabu, possibly revealing a $J_2 \approx 0.44$, implying Borasisi (and/or Pabu) may be a contact binary or an unresolved compact binary. However, full confirmation of this result will require new observations. This work begins the next generation of TNB analyses that go beyond the point mass assumption to provide unique and valuable information on the physical properties of TNBs with implications for their formation and evolution., Comment: Accepted to AJ

We study the optimal scheduling of graph states in measurement-based quantum computation, establishing an equivalence between measurement schedules and path decompositions of graphs. We define the spatial cost of a measurement schedule based on the number of simultaneously active qubits and prove that an optimal measurement schedule corresponds to a path decomposition of minimal width. Our analysis shows that approximating the spatial cost of a graph is $\textsf{NP}$-hard, while for graphs with bounded spatial cost, we establish an efficient algorithm for computing an optimal measurement schedule., Comment: 5 pages, 1 figure

We study the spatially resolved outflow properties of CSWA13, an intermediate mass ($M_*=10^{9}~\mathrm{M}_{\odot}$), gravitationally lensed star-forming galaxy at $z=1.87$. We use Keck/KCWI to map outflows in multiple rest-frame ultraviolet ISM absorption lines, along with fluorescent Si II$^*$ emission, and nebular emission from C III] tracing the local systemic velocity. The spatial structure of outflow velocity mirrors that of the nebular kinematics, which we interpret to be a signature of a young galactic wind that is pressurizing the ISM of the galaxy but is yet to burst out. From the radial extent of Si II$^*$ emission, we estimate that the outflow is largely encapsulated within $3.5$ kpc. We explore the geometry (e.g., patchiness) of the outflow by measuring the covering fraction at different velocities, finding that the maximum covering fraction is at velocities $v\simeq-150$ km$\,$s$^{-1}$. Using the outflow velocity ($v_{out}$), radius ($R$), column density ($N$), and solid angle ($\Omega$) based on the covering fraction, we measure the mass loss rate $\log\dot{m}_{out}/(\mathrm{M}_{\odot}\text{yr}^{-1}) = 1.73\pm0.23$ and mass loading factor $\log\eta = 0.04\pm0.34$ for the low-ionization outflowing gas in this galaxy. These values are relatively large and the bulk of the outflowing gas is moving with speeds less than the escape velocity of the galaxy halo, suggesting that the majority of outflowing mass will remain in the circumgalactic medium and/or recycle back into the galaxy. The results support a picture of high outflow rates transporting mass and metals into the inner circumgalactic medium, providing the gas reservoir for future star formation., Comment: Submitted to ApJ. Comments welcome (17 pages, 10 figures, 1 table)

We generalise results by Lamberti and Lanza de Cristoforis (2005) concerning the continuity of projections onto eigenspaces of self-adjoint differential operators with compact inverses as the (spatial) domain of the functions is perturbed in $\mathbb{R}^2$. Our main case of interest is the Dirichlet Laplacian. We extend these results from bounds from $H_0^1$ to $H_0^1$ to bounds from $L^p$ to $L^p$, under the assumption that $(-\Delta^{-1}-z)^{-1}$ is $L^p$ bounded when $z$ lies outside of the spectrum of $-\Delta^{-1}$. We show that this assumption is met if the initial domain is a square or a rectangle., Comment: 27 pages, 2 figures

Decentralized autonomous organizations (DAOs) are a new, rapidly-growing class of organizations governed by smart contracts. Here we describe how researchers can contribute to the emerging science of DAOs and other digitally-constituted organizations. From granular privacy primitives to mechanism designs to model laws, we identify high-impact problems in the DAO ecosystem where existing gaps might be tackled through a new data set or by applying tools and ideas from existing research fields such as political science, computer science, economics, law, and organizational science. Our recommendations encompass exciting research questions as well as promising business opportunities. We call on the wider research community to join the global effort to invent the next generation of organizations., Comment: includes major coordination problems

The ratios of strong rest-frame optical emission lines are the dominant indicator of metallicities in high-redshift galaxies. Since typical strong-line based metallicity indicators are calibrated on auroral lines at $z=0$, their applicability for galaxies in the distant Universe is unclear. In this paper, we make use of mock emission line data from cosmological simulations to investigate the calibration of rest-frame optical emission lines as metallicity indicators at high redshift. Our model, which couples the SIMBA cosmological galaxy formation simulation with cloudy photoionization calculations, includes contributions from HII regions, post-AGB stars and Diffuse Ionized Gas (DIG). We find mild redshift evolution in the 12 indicators that we study, which implies that the dominant physical properties that evolve in our simulations do have a discernible impact on the metallicity calibrations at high redshifts. When comparing our calibrations with high redshift auroral line observations from James Webb Space Telescope we find a slight offset between our model results and the observations and find that a higher ionization parameter at high redshifts can be one of the possible explanations. We explore the physics that drives the shapes of strong-line metallicity relationships and propose calibrations for hitherto unexplored low-metallicity regimes. Finally, we study the contribution of DIG to total line fluxes. We find that the contribution of DIG increases with metallicity at z $\sim$ 0 for singly ionized oxygen and sulfur lines and can be as high as 70% making it crucial to include their contribution when modeling nebular emission., Comment: 28 pages, 6 figures, typo corrected: in polynomial fits, x represents metallicity and y represents the line ratio

Gait generation for soft robots is challenging due to the nonlinear dynamics and high dimensional input spaces of soft actuators. Limitations in soft robotic control and perception force researchers to hand-craft open loop controllers for gait sequences, which is a non-trivial process. Moreover, short soft actuator lifespans and natural variations in actuator behavior limit machine learning techniques to settings that can be learned on the same time scales as robot deployment. Lastly, simulation is not always possible, due to heterogeneity and nonlinearity in soft robotic materials and their dynamics change due to wear. We present a sample-efficient, simulation free, method for self-generating soft robot gaits, using very minimal computation. This technique is demonstrated on a motorized soft robotic quadruped that walks using four legs constructed from 16 "handed shearing auxetic" (HSA) actuators. To manage the dimension of the search space, gaits are composed of two sequential sets of leg motions selected from 7 possible primitives. Pairs of primitives are executed on one leg at a time; we then select the best-performing pair to execute while moving on to subsequent legs. This method -- which uses no simulation, sophisticated computation, or user input -- consistently generates good translation and rotation gaits in as low as 4 minutes of hardware experimentation, outperforming hand-crafted gaits. This is the first demonstration of completely autonomous gait generation in a soft robot., Comment: 7 Pages, 6 Figures, Published at IROS 2023

We report on a $\rm{[CII]}_{158\mu\rm{m}}$ search using the Atacama Large Millimeter/submillimeter Array (ALMA) on three lensed, confirmed {\lya} emitting galaxies at $z \sim 7$. Our targets are ultra-violet (UV) faint systems with stellar masses on the order of $M_{*} \sim 10^{9} M_{\odot}$. We detect a single [CII] line emission ($4\sigma$) from the brightest ($L \sim 2.4 \times 10^{10}L_{\odot}$) galaxy in our sample, MACS0454-1251. We determine a systemic redshift ($z_{\rm{[CII]}} = 6.3151 \pm 0.0005$) for MACS0454-1251 and measure a {\lya} velocity offset of $\Delta v \approx 300 \pm 70 \rm{km\,s}^{-1}$. The remaining two galaxies we detect no {\ct} but provide $3 \sigma$ upper limits on their {\ct} line luminosities which we use to investigate the $L_{\textrm{[CII]}} - \rm{SFR}$ relation. Overall our single {\ct} detection shows agreement with the relation for dwarf and local starburst galaxies. Our [CII] deficient galaxies could potentially be exhibiting low metallicities ($Z

Predicting sea surface temperature (SST) within the El Ni\~no-Southern Oscillation (ENSO) region has been extensively studied due to its significant influence on global temperature and precipitation patterns. Statistical models such as linear inverse model (LIM), analog forecasting (AF), and recurrent neural network (RNN) have been widely used for ENSO prediction, offering flexibility and relatively low computational expense compared to large dynamic models. However, these models have limitations in capturing spatial patterns in SST variability or relying on linear dynamics. Here we present a modified Convolutional Gated Recurrent Unit (ConvGRU) network for the ENSO region spatio-temporal sequence prediction problem, along with the Ni\~no 3.4 index prediction as a down stream task. The proposed ConvGRU network, with an encoder-decoder sequence-to-sequence structure, takes historical SST maps of the Pacific region as input and generates future SST maps for subsequent months within the ENSO region. To evaluate the performance of the ConvGRU network, we trained and tested it using data from multiple large climate models. The results demonstrate that the ConvGRU network significantly improves the predictability of the Ni\~no 3.4 index compared to LIM, AF, and RNN. This improvement is evidenced by extended useful prediction range, higher Pearson correlation, and lower root-mean-square error. The proposed model holds promise for improving our understanding and predicting capabilities of the ENSO phenomenon and can be broadly applicable to other weather and climate prediction scenarios with spatial patterns and teleconnections., Comment: 13 pages, 7 figures

We establish a general framework for developing approximation algorithms for a class of counting problems. Our framework is based on the cluster expansion of abstract polymer models formalism of Koteck\'y and Preiss. We apply our framework to obtain efficient algorithms for (1) approximating probability amplitudes of a class of quantum circuits close to the identity, (2) approximating expectation values of a class of quantum circuits with operators close to the identity, (3) approximating partition functions of a class of quantum spin systems at high temperature, and (4) approximating thermal expectation values of a class of quantum spin systems at high temperature with positive-semidefinite operators. Further, we obtain hardness of approximation results for approximating probability amplitudes of quantum circuits and partition functions of quantum spin systems. This establishes a computational complexity transition for these problems and shows that our algorithmic conditions are optimal under complexity-theoretic assumptions. Finally, we show that our algorithmic condition is almost optimal for expectation values and optimal for thermal expectation values in the sense of zero freeness., Comment: 22 pages, 0 figures, published version

We report on the discovery of cool gas inflows towards three star-forming galaxies at $\left\sim$ 2.30. Analysis of Keck Low-Resolution Imaging Spectrometer spectroscopy reveals redshifted low-ionisation interstellar (LIS) metal absorption lines with centroid velocities of 60 - 130 km $\rm{s}^{-1}$. These inflows represent some of the most robust detections of inflowing gas into isolated, star-forming galaxies at high redshift. Our analysis suggests that the inflows are due to recycling metal-enriched gas from previous ejections. Comparisons between the galaxies with inflows and a larger parent sample of 131 objects indicate that galaxies with detected inflows may have higher specific star-formation rates (sSFR) and star-formation-rate surface densities. However, when additional galaxies without robustly detected inflows based on centroid velocity but whose LIS absorption line profiles indicate large red-wing velocities are considered, galaxies with inflows do not show unique properties relative to those lacking inflows. Additionally, we calculate the covering fraction of cool inflowing gas as a function of red-wing inflow velocity, finding an enhancement in high sSFR binned galaxies, likely due to an increase in the amount of recycling gas. Together, these results suggest that the low detection rate of galaxies with cool inflows is primarily related to the viewing angle rather than the physical properties of the galaxies., Comment: 10 pages, 7 figures, accepted for publication in MNRAS

We show that a quantum spin system has an exact description by non-interacting fermions if its frustration graph is claw-free and contains a simplicial clique. The frustration graph of a spin model captures the pairwise anticommutation relations between Pauli terms of its Hamiltonian in a given basis. This result captures a vast family of known free-fermion solutions. In previous work, it was shown that a free-fermion solution exists if the frustration graph is either a line graph, or (even-hole, claw)-free. The former case generalizes the celebrated Jordan-Wigner transformation and includes the exact solution to the Kitaev honeycomb model. The latter case generalizes a non-local solution to the four-fermion model given by Fendley. Our characterization unifies these two approaches, extending generalized Jordan-Wigner solutions to the non-local setting and generalizing the four-fermion solution to models of arbitrary spatial dimension. Our key technical insight is the identification of a class of cycle symmetries for all models with claw-free frustration graphs. We prove that these symmetries commute, and this allows us to apply Fendley's solution method to each symmetric subspace independently. Finally, we give a physical description of the fermion modes in terms of operators generated by repeated commutation with the Hamiltonian. This connects our framework to the developing body of work on operator Krylov subspaces. Our results deepen the connection between many-body physics and the mathematical theory of claw-free graphs., Comment: 49 pages, 8 figures

Using deep near-infrared Keck/MOSFIRE observations, we analyze the rest-optical spectra of eight star-forming galaxies in the COSMOS and GOODS-N fields. We reach integration times of $\sim$10 hours in the deepest bands, pushing the limits on current ground-based observational capabilities. The targets fall into two redshift bins -- 5 galaxies at $z \sim 1.7$ and 3 at $z \sim 2.5$ -- and were selected as likely to yield significant auroral-line detections. Even with long integration times, detection of the auroral lines remains challenging. We stack the spectra together into subsets based on redshift, improving the signal-to-noise ratio on the [O III] $\lambda 4364$ auroral emission line and, in turn, enabling a direct measurement of the oxygen abundance for each stack. We compare these measurements to commonly-employed strong-line ratios alongside measurements from the literature. We find that the stacks fall within the distribution of $z>1$ literature measurements, but a larger sample size is needed to robustly constrain the relationships between strong-line ratios and oxygen abundance at high redshift. We additionally report detections of [O I] $\lambda6302$ for nine individual galaxies and composite spectra of 21 targets in the MOSFIRE pointings. We plot their line ratios on the [O III] $\lambda 5008$/H$\beta$ vs. [O I] $\lambda 6302$/H$\alpha$ diagnostic BPT diagram, comparing our targets to local galaxies and H II regions. We find that the [O I]/H$\alpha$ ratios in our sample of galaxies are consistent with being produced in gas ionized by $\alpha$-enhanced massive stars, as has been previously inferred for rapidly-forming galaxies at early cosmic times., Comment: 18 pages, 6 figures, submitted to ApJ

We investigate dust attenuation and its dependence on viewing angle for 308 star-forming galaxies at $1.3\leq z\leq2.6$ from the MOSFIRE Deep Evolution Field (MOSDEF) survey. We divide galaxies with a detected H$\alpha$ emission line and coverage of H$\beta$ into eight groups by stellar mass, star formation rate (SFR), and inclination (i.e., axis ratio), then stack their spectra. From each stack, we measure Balmer decrement and gas-phase metallicity, then we compute median \AV and UV continuum spectral slope ($\beta$). First, we find that none of the dust properties (Balmer decrement, \AV, $\beta$) vary with axis ratio. Second, both stellar and nebular attenuation increase with increasing galaxy mass, showing little residual dependence on SFR or metallicity. Third, nebular emission is more attenuated than stellar emission, and this difference grows even larger at higher galaxy masses and SFRs. Based on these results, we propose a three-component dust model where attenuation predominantly occurs in star-forming regions and large, dusty star-forming clumps, with minimal attenuation in the diffuse ISM. In this model, nebular attenuation primarily originates in clumps, while stellar attenuation is dominated by star-forming regions. Clumps become larger and more common with increasing galaxy mass, creating the above mass trends. Finally, we argue that a fixed metal yield naturally leads to mass regulating dust attenuation. Infall of low-metallicity gas increases SFR and lowers metallicity, but leaves dust column density mostly unchanged. We quantify this idea using the Kennicutt-Schmidt and fundamental metallicity relations, showing that galaxy mass is indeed the primary driver of dust attenuation., Comment: 19 pages, 13 figures, accepted for publication in ApJ