Study Of Galaxy Formation Research Articles

The baryon cycle in modern cosmological hydrodynamical simulations

ABSTRACT In recent years, cosmological hydrodynamical simulations have proven their utility as key interpretative tools in the study of galaxy formation and evolution. In this work, we present a comparative analysis of the baryon cycle in three publicly available, leading cosmological simulation suites: EAGLE, IllustrisTNG, and SIMBA. While these simulations broadly agree in terms of their predictions for the stellar mass content and star formation rates of galaxies at $z\approx 0$, they achieve this result for markedly different reasons. In EAGLE and SIMBA, we demonstrate that at low halo masses ($M_{\rm 200c}\lesssim 10^{11.5}\, \mathrm{M}_{\odot }$), stellar feedback (SF)-driven outflows can reach far beyond the scale of the halo, extending up to $2\!-\!3\times R_{\rm 200c}$. In contrast, in TNG, SF-driven outflows, while stronger at the scale of the interstellar medium, recycle within the circumgalactic medium (within $R_{\rm 200c}$). We find that active galactic nucleus (AGN)-driven outflows in SIMBA are notably potent, reaching several times $R_{\rm 200c}$ even at halo masses up to $M_{\rm 200c}\approx 10^{13.5}\, \mathrm{M}_{\odot }$. In both TNG and EAGLE, AGN feedback can eject gas beyond $R_{\rm 200c}$ at this mass scale, but seldom beyond $2\!-\!3\times R_{\rm 200c}$. We find that the scale of feedback-driven outflows can be directly linked with the prevention of cosmological inflow, as well as the total baryon fraction of haloes within $R_{\rm 200c}$. This work lays the foundation to develop targeted observational tests that can discriminate between feedback scenarios, and inform subgrid feedback models in the next generation of simulations.

CAVITY, Calar Alto Void Integral-field Treasury surveY and project extension

We have learnt in the last decades that the majority of galaxies belong to high density regions interconnected in a sponge-like fashion. This large-scale structure is characterised by clusters, filaments, and walls, where most galaxies concentrate, but also under-dense regions called voids. The void regions and the galaxies within represent an ideal place for the study of galaxy formation and evolution, as they are largely unaffected by the complex physical processes that transform galaxies in high-density environments. The void galaxies may hold the key to answer current challenges to the Lambda CDM paradigm as well. The CAVITY survey is a Legacy project approved by the Calar Alto Observatory to obtain spatially resolved spectroscopic information of $ void galaxies in the Local Universe (0.005 < z < 0.050), covering -17.0 to -21.5 in $ r$ band absolute magnitude. It officially started in January 2021 and has been awarded 110 useful dark observing nights at the 3.5 m telescope using the PMAS spectrograph. Complementary follow-up projects, including deep optical imaging, integrated as well as resolved CO data, and integrated HI spectra, have joined the PMAS observations and naturally complete the scientific aim of characterising galaxies in cosmic voids. The extension data has been named CAVITY+. The data will be available to the whole community in different data releases, the first of which is planned for July 2024, and it will provide the community with PMAS datacubes for around 100 void galaxies through a user friendly and well documented database platform. Here, we present the survey, sample selection, data reduction, quality control schemes, science goals, and some examples of the scientific power of the CAVITY and CAVITY+ data.

UPCluster-SZ: The Updated Catalog of Galaxy Clusters from the List of Planck Sunyaev–Zel’dovich Sources

We present the updated galaxy cluster catalog of the second Planck catalog of Sunyaev–Zel’dovich sources (PSZ2) through the compilation of the data for clusters and galaxies with spectroscopically measured redshifts in the literature. The original version of PSZ2 comprises 1653 Sunyaev–Zel’dovich (SZ) sources, of which 1203 have been validated as genuine galaxy clusters, while the remaining 450 sources are yet to be validated. To increase the number of genuine clusters in PSZ2, we first update the validations of the cluster candidates and their redshift information using the data compiled for the confirmed clusters and the member galaxies in the literature. We then use the galaxy redshift data in the fields of the remaining cluster candidates by searching for possible member galaxies with measured spectroscopic redshifts around the SZ centroids. In this search process, we classify clusters as strong candidates if they contain more than nine galaxies within a 4500 km s−1 velocity range and within 15′ around the SZ centroids. This process results in the validation of 139 new genuine clusters, the update of redshift information on 399 clusters, and the identification of 10 strong candidates, which increases the number of validated clusters up to 1334 among the 1653 SZ sources. Our updated galaxy cluster catalog will be very useful for studies of galaxy formation and cosmology through a combination with other all-sky surveys including the Wide-field Infrared Survey Explorer and SPHEREx.

Emulating the interstellar medium chemistry with neural operators

Context. The study of galaxy formation and evolution critically depends on our understanding of the complex photo-chemical processes that govern the evolution and thermodynamics of the interstellar medium (ISM). In a computational sense, resolving the chemistry is among the weightiest tasks in cosmological and astrophysical simulations. Aims. Astrophysical simulations can include photo-chemical models that allow for a wide range of densities (n), abundances of different species (ni/n) and temperature (T), and plausible evolution scenarios of the ISM under the action of a radiation field (F) with different spectral shapes and intensities. The evolution of such a non-equilibrium photo-chemical network relies on implicit, precise, computationally costly, ordinary differential equations (ODE) solvers. Here, we aim to substitute such procedural solvers with fast, pre-trained emulators based on neural operators. Methods. We emulated a non-equilibrium chemical network up to H2 formation (9 species, 52 reactions) by adopting the DeepONet formalism, namely: by splitting the ODE solver operator that maps the initial conditions and time evolution into a tensor product of two neural networks (named branch and trunk). We used KROME to generate a training set, spanning −2 < log(n/cm−3) ≤ 3.5, log(20) ≤ log(T/K) ≤ 5.5, −6 ≤ log(ni/n) < 0, and adopting an incident radiation field, F, sampled in 10 energy bins with a continuity prior. We separately trained the solver for T and each ni for ≃4.34 GPUhrs. Results. Compared with the reference solutions obtained by KROME for single-zone models, the typical precision obtained is of the order of 10−2, that is, it is 10 times better when using a training that is 40 times less costly, with respect to previous emulators that only considered a fixed F. DeepONet also performs well for T and ni outside the range of the training sample. Furthermore, the emulator aptly reproduces the ion and temperature profiles of photo dissociation regions as well; namely, by giving errors that are comparable to the typical difference between various photo-ionization codes. The present model achieves a speed-up of a factor of 128× with respect to stiff ODE solvers. Conclusions. Our neural emulator represents a significant leap forward in the modelling of ISM chemistry, offering a good balance of precision, versatility, and computational efficiency. Nevertheless, further work is required to address the challenges represented by the extrapolation beyond the training time domain and the removal of potential outliers.

Leaving No Branches Behind: Predicting Baryonic Properties of Galaxies from Merger Trees

Galaxies play a key role in our endeavor to understand how structure formation proceeds in the Universe. For any precision study of cosmology or galaxy formation, there is a strong demand for huge sets of realistic mock galaxy catalogs, spanning cosmologically significant volumes. For such a daunting task, methods that can produce a direct mapping between dark matter halos from dark matter-only simulations and galaxies are strongly preferred, as producing mocks from full-fledged hydrodynamical simulations or semi-analytical models is too expensive. Here, we present a graph-neural-network-based model that is able to accurately predict key properties of galaxies such as stellar mass, g − r color, star formation rate, gas mass, stellar metallicity, and gas metallicity, purely from dark matter properties extracted from halos along the full assembly history of the galaxies. Tests based on the TNG300 simulation of the IllustrisTNG project show that our model can recover the baryonic properties of galaxies to high accuracy, over a wide redshift range (z = 0–5), for all galaxies with stellar masses more massive than 109 M ⊙ and their progenitors, with strong improvements over the state-of-the-art methods. We further show that our method makes substantial strides toward providing an understanding of the implications of the IllustrisTNG galaxy formation model.

A Large Sample of Extremely Metal-poor Galaxies at z < 1 Identified from the DESI Early Data

Extremely metal-poor galaxies (XMPGs) at relatively low redshift are excellent laboratories for studying galaxy formation and evolution in the early universe. Much effort has been spent on identifying them from large-scale spectroscopic surveys or spectroscopic follow-up observations. Previous work has identified a few hundred XMPGs. In this work, we obtain a large sample of 223 XMPGs at z < 1 from the early data of the Dark Energy Spectroscopic Instrument (DESI). The oxygen abundance is determined using the direct T e method based on the detection of the [O iii]λ4363 line. The sample includes 95 confirmed XMPGs based on the oxygen abundance uncertainty; the remaining 128 galaxies are regarded as XMPG candidates. These XMPGs are only 0.01% of the total DESI observed galaxies. Their coordinates and other properties are provided in the paper. The most XMPGs have an oxygen abundance of ∼1/34 Z ⊙, a stellar mass of about 1.5 × 107 M ⊙, and a star formation rate of 0.22 M ⊙ yr−1. The two most XMPGs present distinct morphologies suggesting different formation mechanisms. The local environmental investigation shows that XMPGs preferentially reside in relatively low-density regions. Many of them fall below the stellar mass–metallicity relations (MZRs) of normal star-forming galaxies. From a comparison of the MZR with theoretical simulations, it appears that XMPGs are good analogs to high-redshift star-forming galaxies. The nature of these XMPG populations will be further investigated in detail with larger and more complete samples from the ongoing DESI survey.

Galaxy Interactions in Filaments and Sheets: Insights from EAGLE Simulations

We study the color and star formation rates of paired galaxies in filaments and sheets using the EAGLE simulations. We find that the major pairs with pair separation <50 kpc are bluer and more star-forming in filamentary environments compared to those hosted in sheet-like environments. This trend reverses beyond a pair separation of ∼50 kpc. The interacting pairs with larger separations (>50 kpc) in filaments are on average redder and low-star-forming compared to those embedded in sheets. The galaxies in filaments and sheets may have different stellar mass and cold gas mass distributions. Using a KS test, we find that for paired galaxies with pair separation <50 kpc, there are no significant differences in these properties in sheets and filaments. The filaments transport gas toward the cluster of galaxies. Some earlier studies find preferential alignment of galaxy pairs with the filament axis. Such alignment of galaxy pairs may lead to different gas accretion efficiency in galaxies residing in filaments and sheets. We propose that the enhancement of star formation rate at smaller pair separation in filaments is caused by the alignment of galaxy pairs. A recent study with SDSS data reports the same findings. The confirmation of these results by the EAGLE simulations suggests that the hydrodynamical simulations are powerful theoretical tools for studying galaxy formation and evolution in the cosmic web.

Thesan-hr: how does reionization impact early galaxy evolution?

ABSTRACT The feedback loop between the galaxies producing the background radiation field for reionization and their growth is crucial, particularly for low-mass haloes. Despite this, the vast majority of galaxy formation studies employ a spatially uniform, time-varying reionizing background, with the majority of reionization studies employing galaxy formation models only required to work at high redshift. This paper uses the well-studied TNG galaxy formation model, calibrated at low redshift, coupled to the arepo-rt code, to self-consistently solve the coupled problems of galaxy evolution and reionization, evaluating the impact of patchy (and slow) reionization on early galaxies. thesan-hr is an extension of the thesan project to higher resolution (a factor of 50 increase, with a baryonic mass of mb ≈ 104 M⊙), to additionally enable the study of ‘mini-haloes’ with virial temperatures Tvir &amp;lt; 104 K. Comparing the self-consistent model to a uniform UV background, we show that galaxies in thesan-hr are predicted to be larger in physical extent (by a factor ∼2), less metal enriched (by ∼0.2 dex), and less abundant (by a factor ∼10 at M1500 = − 10) by z = 5. We show that differences in star formation and enrichment patterns lead to significantly different predictions for star formation in low mass haloes, low-metallicity star formation, and even the occupation fraction of haloes. We posit that cosmological galaxy formation simulations aiming to study early galaxy formation (z ≳ 3) must employ a spatially inhomogeneous UV background to accurately reproduce galaxy properties.

The Maunakea Spectroscopic Explorer: Thousands of fibers, infinite possibilities

The Maunakea Spectroscopic Explorer (MSE) is a massively multiplexed spectroscopic survey facility that will replace the Canada–France–Hawai'i Telescope over the next two decades. This 12.5‐m telescope, with its 1.5 square degree field of view, will observe 18,000–20,000 astronomical targets in every pointing from 0.36 to 1.80 m at low/moderate resolution (R3,000, 6,000) and from 0.36 to 0.90 m at high resolution (R30,000). Parallel positioning of all fibers in the field will occur, providing simultaneous full‐field coverage for both resolution modes. Unveiling the composition and dynamics of the faint Universe, MSE will impact nearly every field of astrophysics across all spatial scales, from individual stars to the largest scale structures in the Universe, including (i) the ultimate Gaia follow‐up facility for understanding the chemistry and dynamics of the distant Milky Way, including the distant halo at high spectral resolution, (ii) the unparalleled study of galaxy formation and evolution at cosmic noon, (iii) the determination of the neutrino mass, and (iv) the generation of insights into inflationary physics through a cosmological redshift survey that probes a large volume of the Universe with a high galaxy density. Initially, CFHT will build a Pathfinder instrument to fast‐track the development of MSE technology while providing multi‐object and IFU spectroscopic capability.

Negligible Effects of Baryons on the Angular Momentum Scaling Relations of Galactic Dark Matter Halos

In cosmological simulations without baryons, the relation between the specific angular momentum j h and mass M h of galactic dark matter halos has the well-established form . This is invariably adopted as the starting point in efforts to understand the analogous relation between the specific angular momentum j * and mass M * of the stellar parts of galaxies, which are often re-expressed relative to the corresponding halo properties through the retention fractions f j = j */j h and f M = M */M h. An important caveat here is that the adopted relation could, in principle, be modified by the gravitational back-reaction of baryons on dark matter (DM). We have tested for this possibility by comparing the j h–M h relations in the IllustrisTNG100 and TNG50 simulations that include baryons (full-physics runs) with their counterparts that do not (DM-only runs). In all cases, we find scaling relations of the form , with α ≈ 2/3 over the ranges of mass and redshift studied here: M h ≥ 1010 M ⊙ and 0 ≤ z ≤ 2. The values of α are virtually identical in the full-physics and DM-only runs at the same redshift. The only detectable effect of baryons on the j h–M h relation is a slightly higher normalization, by 12%–15% at z = 0 and by 5% at z = 2. This implies that existing estimates of f j based on DM-only simulations should be adjusted downward by similar amounts. Finally, we discuss briefly some implications of this work for studies of galaxy formation.

Stellar Initial Mass Function (IMF) Probed with Supernova Rates and Neutrino Background: Cosmic-average IMF Slope Is ≃2–3 Similar to the Salpeter IMF

The stellar initial mass function (IMF) is expressed by ϕ(m) ∝ m − α with the slope α, and known as a poorly constrained but very important function in studies of star and galaxy formation. There are no sensible observational constraints on the IMF slopes beyond the Milky Way and nearby galaxies. Here we combine two sets of observational results, (1) cosmic densities of core-collapse supernova (CCSN) explosion rates and (2) cosmic far-UV radiation (and infrared reradiation) densities, which are sensitive to massive (≃8–50 M ⊙) and moderately massive (≃2.5–7 M ⊙) stars, respectively, and constrain the IMF slope at m > 1 M ⊙ with a freedom of redshift evolution. Although no redshift evolution is identified beyond the uncertainties, we find that the cosmic-average IMF slope at z = 0 is α = 1.8–3.2 at the 95% confidence level that is comparable with the Salpeter IMF, α = 2.35, which marks the first constraint on the cosmic-average IMF. We show a forecast for the Nancy Grace Roman Space Telescope supernova survey that will provide significantly strong constraints on the IMF slope with δ α ≃ 0.5 over z = 0–2. Moreover, as for an independent IMF probe instead of (1), we suggest to use diffuse supernovae neutrino background (DSNB), relic neutrinos from CCSNe. We expect that the Hyper-Kamiokande neutrino observations over 20 yr will improve the constraints on the IMF slope and the redshift evolution significantly better than those obtained today, if the systematic uncertainties of DSNB production physics are reduced in the future numerical simulations.

Predictions on the stellar-to-halo mass relation in the dwarf regime using the empirical model for galaxy formation Emerge

ABSTRACT One of the primary goals when studying galaxy formation is to understand how the luminous component of the Universe, galaxies, relate to the growth of structure which is dominated by the gravitational collapse of dark matter haloes. The stellar-to-halo mass relation probes how galaxies occupy dark matter haloes and what that entails for their star formation history. We deliver the first self-consistent empirical model that can place constraints on the stellar-to-halo mass relation down to log stellar mass log10(m*/M⊙) ≤ 5.0 by fitting our model directly to Local Group dwarf data. This is accomplished by penalizing galaxy growth in late-forming, low-mass haloes by mimicking the effects of reionization. This process serves to regulate the number density of galaxies by altering the scatter in halo peak mass $M^{\mathrm{peak}}_{h}$ at fixed stellar mass, creating a tighter scatter than would otherwise exist without a high-z quenching mechanism. Our results indicate that the previously established double-power law stellar-to-halo mass relation can be extended to include galaxies with $\log _{10}(M^{\mathrm{peak}}_{\mathrm{h}}/{\rm M}_{\odot })\gtrsim 10.0$. Furthermore, we show that haloes with $\log _{10}(M^{\mathrm{peak}}_{\mathrm{h}}/{\rm M}_{\odot })\lesssim 9.3$ by z = 4 are unlikely to host a galaxy with log10(m*/M⊙) &amp;gt; 5.0.

Structural Parameter of a Merging Compact Dwarf Galaxy CG 0315

Structural parameters like size, Sersic index, total luminosity and ellipticity are important tools to study galaxy formation and evolution. This work studies a compact dwarf galaxy CG 0325, originated from a merger, and performs a detailed analysis of morphological parameters. The galaxy size is derived by using two different methods, i.e., parametric and non-parametric. The parametric uses the Sersic modeling, but it is found that the observed light profile of CG 0325 cannot be modeled with a simple Sersic function. This is, however, common in merging feature galaxies. Therefore, for further analysis, the half-light radius derived from the non-parametric method can be used as the size of the galaxy. CG 0315 is far more compact than a typical galaxy. We measured its effective radius of 4.05 arcseconds and half-light radius of 2.60 arcseconds that corresponds to a physical value of 221 parsec. Its ratio between the minor and major axis is 0.82, and the position angle is 179 degrees to the North. The star formation rate of the galaxy is 0.051Mʘ year-1. We conclude from these observational data that CG 0315 will continue to be a compact dwarf galaxy when star formation stops and transforms into a well-known compact elliptical galaxy. Therefore, CG 0315 might be a progenitor of an elliptical galaxy.

The “External” Shears In Strong Lens Models

The “External” Shears In Strong Lens Models

PROBES. I. A Compendium of Deep Rotation Curves and Matched Multiband Photometry

We present the Photometry and Rotation Curve Observations from Extragalactic Surveys (PROBES) compendium of extended rotation curves for 3163 late-type spirals, with matching homogeneous multiband photometry for 1677 of them. PROBES rotation curves originally extracted from Hα long-slit spectra and aperture synthesis H i (21cm) velocity maps typically extend out to a median 2 R e (or 1 R 23.5,r ). Our uniform photometry takes advantage of GALEX, DESI-LIS, and WISE images and the software AutoProf to yield multiband azimuthally averaged surface brightness profiles that achieve depths greater than 25 mag arcsec−2 (FUV, NUV), 27 mag arcsec−2 (g, r), and 26 mag arcsec−2 (z, W1, and W2). With its library of spatially resolved profiles and an extensive table of structural parameters, the versatile PROBES data set will benefit studies of galaxy structure and formation.

Modelling the galaxy–halo connection with machine learning

ABSTRACT To extract information from the clustering of galaxies on non-linear scales, we need to model the connection between galaxies and haloes accurately and in a flexible manner. Standard halo occupation distribution (HOD) models make the assumption that the galaxy occupation in a halo is a function of only its mass, however, in reality; the occupation can depend on various other parameters including halo concentration, assembly history, environment, and spin. Using the IllustrisTNG hydrodynamical simulation as our target, we show that machine learning tools can be used to capture this high-dimensional dependence and provide more accurate galaxy occupation models. Specifically, we use a random forest regressor to identify which secondary halo parameters best model the galaxy–halo connection and symbolic regression to augment the standard HOD model with simple equations capturing the dependence on those parameters, namely the local environmental overdensity and shear, at the location of a halo. This not only provides insights into the galaxy formation relationship but also, more importantly, improves the clustering statistics of the modelled galaxies significantly. Our approach demonstrates that machine learning tools can help us better understand and model the galaxy–halo connection, and are therefore useful for galaxy formation and cosmology studies from upcoming galaxy surveys.

Galmask: A Python Package for Unsupervised Galaxy Masking

Galaxy morphological classification is a fundamental aspect of galaxy formation and evolution studies. Various machine learning tools have been developed for automated pipeline analysis of large-scale surveys, enabling a fast search for objects of interest. However, crowded regions in the image may pose a challenge as they can lead to bias in the learning algorithm. In this Research Note, we present galmask, an open-source package for unsupervised galaxy masking to isolate the central object of interest in the image. galmask is written in Python and can be installed from PyPI via the pip command.

The thesan project: predictions for multitracer line intensity mapping in the epoch of reionization

ABSTRACT Line intensity mapping (LIM) is rapidly emerging as a powerful technique to study galaxy formation and cosmology in the high-redshift Universe. We present LIM estimates of select spectral lines originating from the interstellar medium (ISM) of galaxies and 21 cm emission from neutral hydrogen gas in the Universe using the large volume, high resolution thesan reionization simulations. A combination of subresolution photoionization modelling for H ii regions and Monte Carlo radiative transfer calculations is employed to estimate the dust-attenuated spectral energy distributions (SEDs) of high-redshift galaxies (z ≳ 5.5). We show that the derived photometric properties such as the ultraviolet (UV) luminosity function and the UV continuum slopes match observationally inferred values, demonstrating the accuracy of the SED modelling. We provide fits to the luminosity–star formation rate relation (L–SFR) for the brightest emission lines and find that important differences exist between the derived scaling relations and the widely used low-z ones because the ISM of reionization era galaxies is generally less metal enriched than in their low-redshift counterparts. We use these relations to construct line intensity maps of nebular emission lines and cross-correlate with the 21 cm emission. Interestingly, the wavenumber at which the correlation switches sign (ktransition) depends heavily on the reionization model and to a lesser extent on the targeted emission line, which is consistent with the picture that ktransition probes the typical sizes of ionized regions. The derived scaling relations and intensity maps represent a timely state-of-the-art framework for forecasting and interpreting results from current and upcoming LIM experiments.

Galaxy Formation and Reionization: Key Unknowns and Expected Breakthroughs by the James Webb Space Telescope

The launch of the James Webb Space Telescope (JWST) in late 2021 marks a new start for studies of galaxy formation at high redshift ( z ≳ 6) during the era of cosmic reionization. JWST can capture sensitive, high-resolution images and multiobject spectroscopy in the IR that will transform our view of galaxy formation during the first billion years of cosmic history. This review summarizes our current knowledge of the role of galaxies in reionizing intergalactic hydrogen ahead of JWST, achieved through observations with the Hubble Space Telescope and ground-based facilities including Keck, the Very Large Telescope, Subaru, and the Atacama Large Millimeter/Submillimeter Array. We identify outstanding questions in the field that JWST can address during its mission lifetime, including with the planned JWST Cycle 1 programs. These findings include the following: ▪ Surveys with JWST have sufficient sensitivity and area to complete the census of galaxy formation at the current redshift frontier ( z ∼ 8–10). ▪ Rest-frame optical spectroscopy with JWST of galaxies will newly enable measures of star-formation rate, metallicity, and ionization at z ∼ 8–9, allowing for the astrophysics of early galaxies to be constrained. ▪ The presence of evolved stellar populations at z ∼ 8–10 can be definitively tested by JWST, which would provide evidence of star formation out to z ∼ 15.

ALMA Detections of [O iii] and [C ii] Emission Lines From A1689-zD1 at z = 7.13

A1689-zD1 is one of the most distant galaxies, discovered with the aid of gravitational lensing, providing us with an important opportunity to study galaxy formation in the very early universe. In this study, we report the detection of [C ii]158 μm and [O iii]88 μm emission lines of A1689-zD1 in the Atacama Large Millimeter/submillimeter Array (ALMA) Bands 6 and 8. We measure the redshift of this galaxy as z sys = 7.133 ± 0.005 based on the [C ii] and [O iii] emission lines, consistent with that adopted by Bakx et al. The observed L [O III]/L [C II] ratio is 2.09 ± 0.09, higher than that of most of the local galaxies, but consistent with other z ∼ 7 galaxies. The moderate spatial resolution of ALMA data provided us with a precious opportunity to investigate spatial variation of L [O III]/L [C II]. In contrast to the average value of 2.09, we find a much higher L [O III]/L [C II] of ∼7 at the center of the galaxy. This spatial variation of L [O III]/L [C II] was seldom reported for other high-z galaxies. It is also interesting that the peak of the ratio does not overlap with optical peaks. Possible physical reasons include a central active galactic nucleus, shock heating from merging, and a starburst. Our moderate spatial resolution data also reveal that in addition to the observed two clumps shown in previous Hubble Space Telescope images, there is a redshifted segment to the west of the northern optical clump. This structure is consistent with previous claims that A1689-zD1 is a merging galaxy, but with the northern redshifted part being some ejected material, or that the northern redshifted material stems from a third more highly obscured region of the galaxy.