Showing 1–50 of 103 results for author: Snyder, G

The role of lattice thermal conductivity suppression by dopants from a holistic perspective

Authors: Shengnan Dai, Shijie Zhang, Ye Sheng, Erting Dong, Sheng Sun, Lili Xi, G. Jeffrey Snyder, Jinyang Xi, Jiong Yang

Abstract: Dopants play an important role in improving electrical and thermal transport. In the traditional perspective, a dopant suppresses lattice thermal conductivity kL by adding point defect (PD) scattering term to the phonon relaxation time, which has been adopted for decades. In this study, we propose an innovative perspective to solve the kL of defective systems-the holistic approach, i.e., treating… ▽ More Dopants play an important role in improving electrical and thermal transport. In the traditional perspective, a dopant suppresses lattice thermal conductivity kL by adding point defect (PD) scattering term to the phonon relaxation time, which has been adopted for decades. In this study, we propose an innovative perspective to solve the kL of defective systems-the holistic approach, i.e., treating dopant and matrix as a holism. This approach allows us to handle the influences from defects explicitly by the calculations of defective systems, about their changed phonon dispersion, phonon-phonon and electron-phonon interaction, etc, due to the existence of dopants. The kL reduction between defective MxNb1-xFeSb (M=V, Ti) and NbFeSb is used as an example for the holistic approach, and comparable results with experiments are obtained. It is notable that light elemental dopants also induced the avoided-crossing behavior. It can be further rationalized by a one-dimensional atomic chain model. The mass and force constant imbalance generally generates the avoided-crossing phonons, mathematically in a similar way as the coefficients in traditional PD scattering, but along a different direction in kL reduction. Our work provides another perspective for understanding the mechanism of dopants influence in material's thermal transport. △ Less

Submitted 29 June, 2024; originally announced July 2024.

UVCANDELS: The role of dust on the stellar mass-size relation of disk galaxies at 0.5 $\leq z \leq$ 3.0

Authors: Kalina V. Nedkova, Marc Rafelski, Harry I. Teplitz, Vihang Mehta, Laura DeGroot, Swara Ravindranath, Anahita Alavi, Alexander Beckett, Norman A. Grogin, Boris Häußler, Anton M. Koekemoer, Grecco A. Oyarzún, Laura Prichard, Mitchell Revalski, Gregory F. Snyder, Ben Sunnquist, Xin Wang, Rogier A. Windhorst, Nima Chartab, Christopher J. Conselice, Yicheng Guo, Nimish Hathi, Matthew J. Hayes, Zhiyuan Ji, Keunho J. Kim , et al. (8 additional authors not shown)

Abstract: We use the Ultraviolet Imaging of the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey fields (UVCANDELS) to measure half-light radii in the rest-frame far-UV for $\sim$16,000 disk-like galaxies over $0.5\leq z \leq 3$. We compare these results to rest-frame optical sizes that we measure in a self-consistent way and find that the stellar mass-size relation of disk galaxies is steeper… ▽ More We use the Ultraviolet Imaging of the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey fields (UVCANDELS) to measure half-light radii in the rest-frame far-UV for $\sim$16,000 disk-like galaxies over $0.5\leq z \leq 3$. We compare these results to rest-frame optical sizes that we measure in a self-consistent way and find that the stellar mass-size relation of disk galaxies is steeper in the rest-frame UV than in the optical across our entire redshift range. We show that this is mainly driven by massive galaxies ($\gtrsim10^{10}$M$_\odot$), which we find to also be among the most dusty. Our results are consistent with the literature and have commonly been interpreted as evidence of inside-out growth wherein galaxies form their central structures first. However, they could also suggest that the centers of massive galaxies are more heavily attenuated than their outskirts. We distinguish between these scenarios by modeling and selecting galaxies at $z=2$ from the VELA simulation suite in a way that is consistent with UVCANDELS. We show that the effects of dust alone can account for the size differences we measure at $z=2$. This indicates that, at different wavelengths, size differences and the different slopes of the stellar mass-size relation do not constitute evidence for inside-out growth. △ Less

Submitted 28 June, 2024; v1 submitted 17 May, 2024; originally announced May 2024.

Comments: Accepted for publication in ApJ. 22 pages, 12 figures, and 4 tables

Why does silicon have an indirect band gap?

Authors: Emily Oliphant, Veda Mantena, Madison Brod, G. Jeffrey Snyder, Wenhao Sun

Abstract: Density functional theory (DFT) is a powerful tool for predicting the electronic band structures of functional materials. However, it is often difficult to intuit how major band structure features$-$such as band gap magnitude, location of band extrema, effective masses, etc.$-$arise from the underlying crystal chemistry of a material. Band structure is represented in reciprocal space, but arises f… ▽ More Density functional theory (DFT) is a powerful tool for predicting the electronic band structures of functional materials. However, it is often difficult to intuit how major band structure features$-$such as band gap magnitude, location of band extrema, effective masses, etc.$-$arise from the underlying crystal chemistry of a material. Band structure is represented in reciprocal space, but arises from the orbital and bonding interactions between atoms in real space. Here, we present a conceptual and computable framework to extract chemical bonding origins of DFT-calculated band structure features. As a key example, we explain here the indirect band gap of silicon, which profoundly impacts its properties for use in photovoltaics and electronics. Even in this basic semiconductor material, our approach leads to new insights to understand and engineer its conduction band minimum position. These calculation techniques can be broadly applied to reveal the crystal chemistry origins of electronic structure features in other optical, electronic, and magnetic materials. △ Less

Submitted 9 April, 2024; v1 submitted 8 April, 2024; originally announced April 2024.

The Next Generation Deep Extragalactic Exploratory Public Near-Infrared Slitless Survey Epoch 1 (NGDEEP-NISS1): Extra-Galactic Star-formation and Active Galactic Nuclei at 0.5 < z < 3.6

Authors: Nor Pirzkal, Barry Rothberg, Casey Papovich, Lu Shen, Gene C. K. Leung, Micaela B. Bagley, Steven L. Finkelstein, Brittany N. Vanderhoof, Jennifer M. Lotz, Anton M. Koekemoer, Nimish P. Hathi, Yingjie Cheng, Nikko J. Cleri, Norman A. Grogin, L. Y. Aaron Yung, Mark Dickinson, Henry C. Ferguson, Jonathan P. Gardner, Intae Jung, Jeyhan S. Kartaltepe, Russell Ryan, Raymond C. Simons, Swara Ravindranath, Danielle A. Berg, Bren E. Backhaus , et al. (26 additional authors not shown)

Abstract: The Next Generation Deep Extragalactic Exploratory Public (NGDEEP) survey program was designed specifically to include Near Infrared Slitless Spectroscopic observations (NGDEEP-NISS) to detect multiple emission lines in as many galaxies as possible and across a wide redshift range using the Near Infrared Imager and Slitless Spectrograph (NIRISS). We present early results obtained from the the firs… ▽ More The Next Generation Deep Extragalactic Exploratory Public (NGDEEP) survey program was designed specifically to include Near Infrared Slitless Spectroscopic observations (NGDEEP-NISS) to detect multiple emission lines in as many galaxies as possible and across a wide redshift range using the Near Infrared Imager and Slitless Spectrograph (NIRISS). We present early results obtained from the the first set of observations (Epoch 1, 50$\%$ of the allocated orbits) of this program (NGDEEP-NISS1). Using a set of independently developed calibration files designed to deal with a complex combination of overlapping spectra, multiple position angles, and multiple cross filters and grisms, in conjunction with a robust and proven algorithm for quantifying contamination from overlapping dispersed spectra, NGDEEP-NISS1 has achieved a 3$σ$ sensitivity limit of 2 $\times$ 10$^{-18}$ erg/s/cm$^2$. We demonstrate the power of deep wide field slitless spectroscopy (WFSS) to characterize the star-formation rates, and metallicity ([OIII]/H$β$), and dust content, of galaxies at $1<z<3.5$. The latter showing intriguing initial results on the applicability and assumptions made regarding the use of Case B recombination. Further, we identify the presence of active galactic nuclei (AGN) and infer the mass of their supermassive black holes (SMBHs) using broadened restframe MgII and H$β$ emission lines. The spectroscopic results are then compared with the physical properties of galaxies extrapolated from fitting spectral energy distribution (SED) models to photometry alone. The results clearly demonstrate the unique power and efficiency of WFSS at near-infrared wavelengths over other methods to determine the properties of galaxies across a broad range of redshifts. △ Less

Submitted 20 April, 2024; v1 submitted 15 December, 2023; originally announced December 2023.

Comments: 33 pages, 28 Figures, Accepted (ApJ)

Coronal Heating as Determined by the Solar Flare Frequency Distribution Obtained by Aggregating Case Studies

Authors: James Paul Mason, Alexandra Werth, Colin G. West, Allison A. Youngblood, Donald L. Woodraska, Courtney Peck, Kevin Lacjak, Florian G. Frick, Moutamen Gabir, Reema A. Alsinan, Thomas Jacobsen, Mohammad Alrubaie, Kayla M. Chizmar, Benjamin P. Lau, Lizbeth Montoya Dominguez, David Price, Dylan R. Butler, Connor J. Biron, Nikita Feoktistov, Kai Dewey, N. E. Loomis, Michal Bodzianowski, Connor Kuybus, Henry Dietrick, Aubrey M. Wolfe , et al. (977 additional authors not shown)

Abstract: Flare frequency distributions represent a key approach to addressing one of the largest problems in solar and stellar physics: determining the mechanism that counter-intuitively heats coronae to temperatures that are orders of magnitude hotter than the corresponding photospheres. It is widely accepted that the magnetic field is responsible for the heating, but there are two competing mechanisms th… ▽ More Flare frequency distributions represent a key approach to addressing one of the largest problems in solar and stellar physics: determining the mechanism that counter-intuitively heats coronae to temperatures that are orders of magnitude hotter than the corresponding photospheres. It is widely accepted that the magnetic field is responsible for the heating, but there are two competing mechanisms that could explain it: nanoflares or Alfvén waves. To date, neither can be directly observed. Nanoflares are, by definition, extremely small, but their aggregate energy release could represent a substantial heating mechanism, presuming they are sufficiently abundant. One way to test this presumption is via the flare frequency distribution, which describes how often flares of various energies occur. If the slope of the power law fitting the flare frequency distribution is above a critical threshold, $α=2$ as established in prior literature, then there should be a sufficient abundance of nanoflares to explain coronal heating. We performed $>$600 case studies of solar flares, made possible by an unprecedented number of data analysts via three semesters of an undergraduate physics laboratory course. This allowed us to include two crucial, but nontrivial, analysis methods: pre-flare baseline subtraction and computation of the flare energy, which requires determining flare start and stop times. We aggregated the results of these analyses into a statistical study to determine that $α= 1.63 \pm 0.03$. This is below the critical threshold, suggesting that Alfvén waves are an important driver of coronal heating. △ Less

Submitted 9 May, 2023; originally announced May 2023.

Comments: 1,002 authors, 14 pages, 4 figures, 3 tables, published by The Astrophysical Journal on 2023-05-09, volume 948, page 71

The Physical Thickness of Stellar Disks to z ~ 2

Authors: Kathleen A. Hamilton-Campos, Raymond C. Simons, Molly S. Peeples, Gregory F. Snyder, Timothy M. Heckman

Abstract: In local disk galaxies such as our Milky Way, older stars generally inhabit a thicker disk than their younger counterparts. Two competing models have attempted to explain this result: one in which stars first form in thin disks that gradually thicken with time through dynamical heating, and one in which stars form in thick disks at early times and in progressively thinner disks at later times. We… ▽ More In local disk galaxies such as our Milky Way, older stars generally inhabit a thicker disk than their younger counterparts. Two competing models have attempted to explain this result: one in which stars first form in thin disks that gradually thicken with time through dynamical heating, and one in which stars form in thick disks at early times and in progressively thinner disks at later times. We use a direct measure of the thicknesses of stellar disks at high redshift to discriminate between these scenarios. Using legacy HST imaging from the CANDELS and GOODS surveys, we measure the rest-optical scale heights of 491 edge-on disk galaxies spanning 0.4 < z < 2.5. We measure a median intrinsic scale height for the full sample of 0.74 +/- 0.03 kpc, with little redshift evolution of both the population median and scatter. The median is consistent with the thick disk of the Milky Way today (0.6 - 1.1 kpc), but is smaller than the median scale height of local disks (~1.5 kpc) which are matched to our high-redshift sample by descendant mass. These findings indicate that (1) while disks as thick as the Milky Way's thick disk were in place at early times, (2) to explain the full disk galaxy population today, the stellar disks in galaxies need to on average physically thicken after formation. △ Less

Submitted 7 March, 2023; originally announced March 2023.

Comments: 14 pages, 8 figures, submitted to ApJ

The Next Generation Deep Extragalactic Exploratory Public (NGDEEP) Survey

Authors: Micaela B. Bagley, Nor Pirzkal, Steven L. Finkelstein, Casey Papovich, Danielle A. Berg, Jennifer M. Lotz, Gene C. K. Leung, Henry C. Ferguson, Anton M. Koekemoer, Mark Dickinson, Jeyhan S. Kartaltepe, Dale D. Kocevski, Rachel S. Somerville, L. Y. Aaron Yung, Bren E. Backhaus, Caitlin M. Casey, Marco Castellano, Óscar A. Chávez Ortiz, Katherine Chworowsky, Isabella G. Cox, Romeel Davé, Kelcey Davis, Vicente Estrada-Carpenter, Adriano Fontana, Seiji Fujimoto , et al. (23 additional authors not shown)

Abstract: We present the Next Generation Deep Extragalactic Exploratory Public (NGDEEP) Survey, a deep slitless spectroscopic and imaging Cycle 1 JWST treasury survey designed to constrain feedback mechanisms in low-mass galaxies across cosmic time. NGDEEP targets the Hubble Ultra Deep Field (HUDF) with NIRISS slitless spectroscopy (f~1.2e-18 erg/s/cm^2, 5sigma) to measure metallicities and star-formation r… ▽ More We present the Next Generation Deep Extragalactic Exploratory Public (NGDEEP) Survey, a deep slitless spectroscopic and imaging Cycle 1 JWST treasury survey designed to constrain feedback mechanisms in low-mass galaxies across cosmic time. NGDEEP targets the Hubble Ultra Deep Field (HUDF) with NIRISS slitless spectroscopy (f~1.2e-18 erg/s/cm^2, 5sigma) to measure metallicities and star-formation rates (SFRs) for low-mass galaxies through the peak of the cosmic SFR density (0.5<z<4). In parallel, NGDEEP targets the HUDF-Par2 parallel field with NIRCam (m=30.6-30.9, 5sigma) to discover galaxies to z>12, constraining the slope of the faint-end of the rest-ultraviolet luminosity function. NGDEEP overlaps with the deepest HST ACS optical imaging in the sky: F435W in the HUDF (m=29.6), and F814W in HUDF-Par2 (m=30), making this a premier HST+JWST Deep Field. As a treasury survey, NGDEEP data is public immediately, and we will rapidly release data products and catalogs in the spirit of previous deep field initiatives. In this paper we present the NGDEEP survey design, summarize the science goals, and detail plans for the public release of NGDEEP reduced data products. △ Less

Submitted 10 February, 2023; originally announced February 2023.

Comments: 21 pages, 9 figures, submitted to AAS Journals

Mock Galaxy Surveys for HST and JWST from the IllustrisTNG Simulations

Authors: Gregory F. Snyder, Theodore Pena, L. Y. Aaron Yung, Caitlin Rose, Jeyhan Kartaltepe, Harry Ferguson

Abstract: We present and analyze a series of synthetic galaxy survey fields based on the IllustrisTNG Simulation suite. With the Illustris public data release and JupyterLab service, we generated a set of twelve lightcone catalogs covering areas from 5 to 365 square arcminutes, similar to several JWST Cycle 1 programs, including JADES, CEERS, PRIMER, and NGDEEP. From these catalogs, we queried the public AP… ▽ More We present and analyze a series of synthetic galaxy survey fields based on the IllustrisTNG Simulation suite. With the Illustris public data release and JupyterLab service, we generated a set of twelve lightcone catalogs covering areas from 5 to 365 square arcminutes, similar to several JWST Cycle 1 programs, including JADES, CEERS, PRIMER, and NGDEEP. From these catalogs, we queried the public API to generate simple mock images in a series of broadband filters used by JWST-NIRCam and the Hubble Space Telescope cameras. This procedure generates wide-area simulated mosaic images that can support investigating the predicted evolution of galaxies alongside real data. Using these mocks, we demonstrate a few simple science cases, including morphological evolution and close pair selection. We publicly release the catalogs and mock images through MAST, along with the code used to generate these projects, so that the astrophysics community can make use of these products in their scientific analyses of JWST deep field observations. △ Less

Submitted 17 November, 2022; originally announced November 2022.

Comments: Accepted to MNRAS

Morphological signatures of mergers in the TNG50 simulation and the Kilo-Degree Survey: the merger fraction from dwarfs to Milky Way-like galaxies

Authors: Alejandro Guzmán-Ortega, Vicente Rodriguez-Gomez, Gregory F. Snyder, Katie Chamberlain, Lars Hernquist

Abstract: Using the TNG50 cosmological simulation and observations from the Kilo-Degree Survey (KiDS), we investigate the connection between galaxy mergers and optical morphology in the local Universe over a wide range of galaxy stellar masses ($8.5\leqslant\log(M_\ast/\text{M}_\odot)\leqslant11$). To this end, we have generated over 16,000 synthetic images of TNG50 galaxies designed to match KiDS observati… ▽ More Using the TNG50 cosmological simulation and observations from the Kilo-Degree Survey (KiDS), we investigate the connection between galaxy mergers and optical morphology in the local Universe over a wide range of galaxy stellar masses ($8.5\leqslant\log(M_\ast/\text{M}_\odot)\leqslant11$). To this end, we have generated over 16,000 synthetic images of TNG50 galaxies designed to match KiDS observations, including the effects of dust attenuation and scattering, and used the $\mathrm{\mathtt{statmorph}}$ code to measure various image-based morphological diagnostics in the $r$-band for both data sets. Such measurements include the Gini-$M_{20}$ and concentration-asymmetry-smoothness statistics. Overall, we find good agreement between the optical morphologies of TNG50 and KiDS galaxies, although the former are slightly more concentrated and asymmetric than their observational counterparts. Afterwards, we trained a random forest classifier to identify merging galaxies in the simulation (including major and minor mergers) using the morphological diagnostics as the model features, along with merger statistics from the merger trees as the ground truth. We find that the asymmetry statistic exhibits the highest feature importance of all the morphological parameters considered. Thus, the performance of our algorithm is comparable to that of the more traditional method of selecting highly asymmetric galaxies. Finally, using our trained model, we estimate the galaxy merger fraction in both our synthetic and observational galaxy samples, finding in both cases that the galaxy merger fraction increases steadily as a function of stellar mass. △ Less

Submitted 10 November, 2022; originally announced November 2022.

Comments: 18 pages, 12 figures. Accepted for publication in MNRAS

Effects of feedback on galaxies in the VELA simulations: elongation, clumps and compaction

Authors: Daniel Ceverino, Nir Mandelker, Gregory F. Snyder, Sharon Lapiner, Avishai Dekel, Joel Primack, Omri Ginzburg, Sean Larkin

Abstract: The evolution of star-forming galaxies at high redshifts is very sensitive to the strength and nature of stellar feedback. Using two sets of cosmological, zoom-in simulations from the VELA suite, we compare the effects of two different models of feedback: with and without kinetic feedback from the expansion of supernovae shells and stellar winds. At a fixed halo mass and redshift, the stellar mass… ▽ More The evolution of star-forming galaxies at high redshifts is very sensitive to the strength and nature of stellar feedback. Using two sets of cosmological, zoom-in simulations from the VELA suite, we compare the effects of two different models of feedback: with and without kinetic feedback from the expansion of supernovae shells and stellar winds. At a fixed halo mass and redshift, the stellar mass is reduced by a factor of 1-3 in the models with stronger feedback, so the stellar-mass-halo-mass relation is in better agreement with abundance matching results. On the other hand, the three-dimensional shape of low-mass galaxies is elongated along a major axis in both models. At a fixed stellar mass, Ms<10^10 Msun, galaxies are more elongated in the strong-feedback case. More massive, star-forming discs with high surface densities form giant clumps. However, the population of round, compact, old (age_c > 300 Myr), quenched, stellar (or gas-poor) clumps is absent in the model with strong feedback. On the other hand, giant star-forming clumps with intermediate ages (age_c= 100 - 300 Myr) can survive for several disc dynamical times, independently of feedback strength. The evolution through compaction followed by quenching in the plane of central surface density and specific star-formation rate is similar under the two feedback models. △ Less

Submitted 25 April, 2023; v1 submitted 27 October, 2022; originally announced October 2022.

Comments: 14 pages, 12 figures, accepted version at MNRAS

CEERS Key Paper III: The Diversity of Galaxy Structure and Morphology at z=3-9 with JWST

Authors: Jeyhan S. Kartaltepe, Caitlin Rose, Brittany N. Vanderhoof, Elizabeth J. McGrath, Luca Costantin, Isabella G. Cox, L. Y. Aaron Yung, Dale D. Kocevski, Stijn Wuyts, Henry C. Ferguson Brett H. Andrews, Micaela B. Bagley, Steven L. Finkelstein, Ricardo O. Amorin, Pablo Arrabal Haro, Bren E. Backhaus, Peter Behroozi, Laura Bisigello, Antonello Calabro, Caitlin M. Casey, Rosemary T. Coogan, Darren Croton, Alexander de la Vega, Mark Dickinson, M. C. Cooper, Adriano Fontana , et al. (36 additional authors not shown)

Abstract: We present a comprehensive analysis of the evolution of the morphological and structural properties of a large sample of galaxies at z=3-9 using early JWST CEERS NIRCam observations. Our sample consists of 850 galaxies at z>3 detected in both CANDELS HST imaging and JWST CEERS NIRCam images to enable a comparison of HST and JWST morphologies. Our team conducted a set of visual classifications, wit… ▽ More We present a comprehensive analysis of the evolution of the morphological and structural properties of a large sample of galaxies at z=3-9 using early JWST CEERS NIRCam observations. Our sample consists of 850 galaxies at z>3 detected in both CANDELS HST imaging and JWST CEERS NIRCam images to enable a comparison of HST and JWST morphologies. Our team conducted a set of visual classifications, with each galaxy in the sample classified by three different individuals. We also measure quantitative morphologies using the publicly available codes across all seven NIRCam filters. Using these measurements, we present the fraction of galaxies of each morphological type as a function of redshift. Overall, we find that galaxies at z>3 have a wide diversity of morphologies. Galaxies with disks make up a total of 60\% of galaxies at z=3 and this fraction drops to ~30% at z=6-9, while galaxies with spheroids make up ~30-40% across the whole redshift range and pure spheroids with no evidence for disks or irregular features make up ~20%. The fraction of galaxies with irregular features is roughly constant at all redshifts (~40-50%), while those that are purely irregular increases from ~12% to ~20% at z>4.5. We note that these are apparent fractions as many selection effects impact the visibility of morphological features at high redshift. The distributions of Sérsic index, size, and axis ratios show significant differences between the morphological groups. Spheroid Only galaxies have a higher Sérsic index, smaller size, and higher axis ratio than Disk/Irregular galaxies. Across all redshifts, smaller spheroid and disk galaxies tend to be rounder. Overall, these trends suggest that galaxies with established disks and spheroids exist across the full redshift range of this study and further work with large samples at higher redshift is needed to quantify when these features first formed. △ Less

Submitted 13 January, 2023; v1 submitted 26 October, 2022; originally announced October 2022.

Comments: Accepted for publication to ApJL, 24 pages, 14 figures

High Resolution HST Imaging Survey of Local Star-Forming Galaxies I: Spatially-Resolved Obscured Star Formation with H$α$ and Paschen-$β$ Recombination Lines

Authors: Clara Giménez-Arteaga, Gabriel B. Brammer, Danilo Marchesini, Luis Colina, Varun Bajaj, Malte Brinch, Daniela Calzetti, Daniel Lange-Vagle, Eric J. Murphy, Michele Perna, Javier Piqueras-López, Gregory F. Snyder

Abstract: We present a sample of 24 local star-forming galaxies observed with broad- and narrow-band photometry from the Hubble Space Telescope, that are part of the GOALS survey of local luminous and ultra-luminous infrared galaxies. With narrow-band filters around the emission lines H$α$ (and [NII]) and Pa$β$, we obtain robust estimates of the dust attenuation affecting the gas in each galaxy, probing hig… ▽ More We present a sample of 24 local star-forming galaxies observed with broad- and narrow-band photometry from the Hubble Space Telescope, that are part of the GOALS survey of local luminous and ultra-luminous infrared galaxies. With narrow-band filters around the emission lines H$α$ (and [NII]) and Pa$β$, we obtain robust estimates of the dust attenuation affecting the gas in each galaxy, probing higher attenuation than can be traced by the optical Balmer decrement H$α$/H$β$ alone by a factor of $>1$ mag. We also infer the dust attenuation towards the stars via a spatially-resolved SED-fitting procedure that uses all available HST imaging filters. We use various indicators to obtain the star formation rate (SFR) per spatial bin, and find that Pa$β$ traces star-forming regions where the H$α$ and the optical stellar continuum are heavily obscured. The dust-corrected Pa$β$ SFR recovers the 24$μ$m-inferred SFR with a ratio $-0.14\pm0.32$ dex and the SFR inferred from the $8\mathrm{-}1000\,μ\mathrm{m}$ infrared luminosity at $-0.04\pm0.23$ dex. Both in a spatially-resolved and integrated sense, rest-frame near infrared recombination lines can paint a more comprehensive picture of star formation across cosmic time, particularly with upcoming JWST observations of Paschen-series line emission in galaxies as early as the epoch of reionization. △ Less

Submitted 30 September, 2022; originally announced October 2022.

Comments: Accepted for publication in ApJS

Identifying Galaxy Mergers in Simulated CEERS NIRCam Images using Random Forests

Authors: Caitlin Rose, Jeyhan S. Kartaltepe, Gregory F. Snyder, Vicente Rodriguez-Gomez, L. Y. Aaron Yung, Pablo Arrabal Haro, Micaela B. Bagley, Antonello Calabrò, Nikko J. Cleri, M. C. Cooper, Luca Costantin, Darren Croton, Mark Dickinson, Steven L. Finkelstein, Boris Häußler, Benne W. Holwerda, Anton M. Koekemoer, Peter Kurczynski, Ray A. Lucas, Kameswara Bharadwaj Mantha, Casey Papovich, Pablo G. Pérez-González, Nor Pirzkal, Rachel S. Somerville, Amber N. Straughn , et al. (1 additional authors not shown)

Abstract: Identifying merging galaxies is an important - but difficult - step in galaxy evolution studies. We present random forest classifications of galaxy mergers from simulated JWST images based on various standard morphological parameters. We describe (a) constructing the simulated images from IllustrisTNG and the Santa Cruz SAM, and modifying them to mimic future CEERS observations as well as nearly n… ▽ More Identifying merging galaxies is an important - but difficult - step in galaxy evolution studies. We present random forest classifications of galaxy mergers from simulated JWST images based on various standard morphological parameters. We describe (a) constructing the simulated images from IllustrisTNG and the Santa Cruz SAM, and modifying them to mimic future CEERS observations as well as nearly noiseless observations, (b) measuring morphological parameters from these images, and (c) constructing and training the random forests using the merger history information for the simulated galaxies available from IllustrisTNG. The random forests correctly classify $\sim60\%$ of non-merging and merging galaxies across $0.5 < z < 4.0$. Rest-frame asymmetry parameters appear more important for lower redshift merger classifications, while rest-frame bulge and clump parameters appear more important for higher redshift classifications. Adjusting the classification probability threshold does not improve the performance of the forests. Finally, the shape and slope of the resulting merger fraction and merger rate derived from the random forest classifications match with theoretical Illustris predictions, but are underestimated by a factor of $\sim 0.5$. △ Less

Submitted 23 August, 2022; originally announced August 2022.

Comments: 21 pages, 14 figures, submitted to ApJ

Dusty Starbursts Masquerading as Ultra-high Redshift Galaxies in JWST CEERS Observations

Authors: Jorge A. Zavala, Veronique Buat, Caitlin M. Casey, Denis Burgarella, Steven L. Finkelstein, Micaela B. Bagley, Laure Ciesla, Emanuele Daddi, Mark Dickinson, Henry C. Ferguson, Maximilien Franco, E. F. Jim'enez-Andrade, Jeyhan S. Kartaltepe, Anton M. Koekemoer, Aurélien Le Bail, E. J. Murphy, Casey Papovich, Sandro Tacchella, Stephen M. Wilkins, Itziar Aretxaga, Peter Behroozi, Jaclyn B. Champagne, Adriano Fontana, Mauro Giavalisco, Andrea Grazian , et al. (99 additional authors not shown)

Abstract: Lyman Break Galaxy (LBG) candidates at z>10 are rapidly being identified in JWST/NIRCam observations. Due to the (redshifted) break produced by neutral hydrogen absorption of rest-frame UV photons, these sources are expected to drop out in the bluer filters while being well detected in redder filters. However, here we show that dust-enshrouded star-forming galaxies at lower redshifts (z<7) may als… ▽ More Lyman Break Galaxy (LBG) candidates at z>10 are rapidly being identified in JWST/NIRCam observations. Due to the (redshifted) break produced by neutral hydrogen absorption of rest-frame UV photons, these sources are expected to drop out in the bluer filters while being well detected in redder filters. However, here we show that dust-enshrouded star-forming galaxies at lower redshifts (z<7) may also mimic the near-infrared (near-IR) colors of z>10 LBGs, representing potential contaminants in LBG candidate samples. First, we analyze CEERS-DSFG-1, a NIRCam dropout undetected in the F115W and F150W filters but detected at longer wavelengths. Combining the JWST data with (sub)millimeter constraints, including deep NOEMA interferometric observations, we show that this source is a dusty star-forming galaxy (DSFG) at z~5.1. We also present a tentative 2.6sigma SCUBA-2 detection at 850um around a recently identified z~16 LBG candidate in the same field and show that, if the emission is real and associated with this candidate, the available photometry is consistent with a z~5 dusty galaxy with strong nebular emission lines despite its blue near-IR colors. Further observations on this candidate are imperative to mitigate the low confidence of this tentative submillimeter emission and its positional uncertainty. Our analysis shows that robust (sub)millimeter detections of NIRCam dropout galaxies likely imply z=4-6 redshift solutions, where the observed near-IR break would be the result of a strong rest-frame optical Balmer break combined with high dust attenuation and strong nebular line emission, rather than the rest-frame UV Lyman break. This provides evidence that DSFGs may contaminate searches for ultra high-redshift LBG candidates from JWST observations. △ Less

Submitted 30 January, 2023; v1 submitted 2 August, 2022; originally announced August 2022.

Comments: Published in The Astrophysical Journal Letters (updated to match the published version)

Journal ref: ApJL 943 L9 (2023)

Probing the earliest phases in the formation of massive galaxies with simulated HST+JWST imaging data from Illustris

Authors: Ángela García-Argumánez, Pablo G. Pérez-González, Armando Gil de Paz, Gregory F. Snyder, Pablo Arrabal Haro, Micaela B. Bagley, Steven L. Finkelstein, Jeyhan S. Kartaltepe, Anton Koekemoer, Casey Papovich, Nor Pirzkal, Harry C. Ferguson, L. Y. Aaron Yung, Marianna Annunziatella, Nikko J. Cleri, M. C. Cooper, Luca Costantin, Benne W. Holwerda, Rosa María Mérida González, Caitlin Rose, Mauro Giavalisco, Norman A. Grogin, Dale D. Kocevski

Abstract: We use the Illustris-1 simulation to explore the capabilities of the $\textit{Hubble}$ and $\textit{James Webb Space Telescope}$ data to analyze the stellar populations in high-redshift galaxies, taking advantage of the combined depth, spatial resolution, and wavelength coverage. For that purpose, we use simulated broad-band ACS, WFC3 and NIRCam data and 2-dimensional stellar population synthesis… ▽ More We use the Illustris-1 simulation to explore the capabilities of the $\textit{Hubble}$ and $\textit{James Webb Space Telescope}$ data to analyze the stellar populations in high-redshift galaxies, taking advantage of the combined depth, spatial resolution, and wavelength coverage. For that purpose, we use simulated broad-band ACS, WFC3 and NIRCam data and 2-dimensional stellar population synthesis (2D-SPS) to derive the integrated star formation history (SFH) of massive (M$_{\ast}>10^{10}\,$M$_{\odot}$) simulated galaxies at $1<z<4$ that evolve into a local M$_{\ast}>10^{11}\,$M$_{\odot}$ galaxy. In particular, we explore the potential of HST and JWST datasets reaching a depth similar to those of the CANDELS and ongoing CEERS observations, respectively, and concentrate on determining the capabilities of this dataset for characterizing the first episodes in the SFH of local M$_{\ast}>10^{11}\,$M$_{\odot}$ galaxies by studying their progenitors at $z>1$. The 2D-SPS method presented in this paper has been calibrated to robustly recover the cosmic times when the first star formation episodes occurred in massive galaxies, i.e., the first stages in their integrated SFHs. In particular, we discuss the times when the first 1% to 50% of their total stellar mass formed in the simulation. We demonstrate that we can recover these ages with typical median systematic offset of less than 5% and scatter around 20%-30%. According to our measurements on Illustris data, we are able to recover that local M$_{\ast}>10^{11}\,$M$_{\odot}$ galaxies would have started their formation by $z=16$, forming the first 5% of their stellar mass present at $z \sim 1$ by $z=4.5$, 10% by $z=3.7$, and 25% by $z=2.7$. △ Less

Submitted 19 December, 2022; v1 submitted 28 July, 2022; originally announced July 2022.

Comments: 28 pages, 13 figures, 4 tables. ApJ in press. Summary of changes from original submission: the major change is that we now include in Sec. 6 the comparison of the results obtained for our sample of massive 1 < z < 4 progenitors with those obtained by considering all massive galaxies at 1 < z < 4 in the simulated images. Several figures and sections have been updated

A Long Time Ago in a Galaxy Far, Far Away: A Candidate z ~ 12 Galaxy in Early JWST CEERS Imaging

Authors: Steven L. Finkelstein, Micaela B. Bagley, Pablo Arrabal Haro, Mark Dickinson, Henry C. Ferguson, Jeyhan S. Kartaltepe, Casey Papovich, Denis Burgarella, Dale D. Kocevski, Marc Huertas-Company, Kartheik G. Iyer, Rebecca L. Larson, Pablo G. Pérez-González, Caitlin Rose, Sandro Tacchella, Stephen M. Wilkins, Katherine Chworowsky, Aubrey Medrano, Alexa M. Morales, Rachel S. Somerville, L. Y. Aaron Yung, Adriano Fontana, Mauro Giavalisco, Andrea Grazian, Norman A. Grogin , et al. (95 additional authors not shown)

Abstract: We report the discovery of a candidate galaxy with a photo-z of z~12 in the first epoch of the JWST Cosmic Evolution Early Release Science (CEERS) Survey. Following conservative selection criteria we identify a source with a robust z_phot = 11.8^+0.3_-0.2 (1-sigma uncertainty) with m_F200W=27.3, and >7-sigma detections in five filters. The source is not detected at lambda < 1.4um in deep imaging f… ▽ More We report the discovery of a candidate galaxy with a photo-z of z~12 in the first epoch of the JWST Cosmic Evolution Early Release Science (CEERS) Survey. Following conservative selection criteria we identify a source with a robust z_phot = 11.8^+0.3_-0.2 (1-sigma uncertainty) with m_F200W=27.3, and >7-sigma detections in five filters. The source is not detected at lambda < 1.4um in deep imaging from both HST and JWST, and has faint ~3-sigma detections in JWST F150W and HST F160W, which signal a Ly-alpha break near the red edge of both filters, implying z~12. This object (Maisie's Galaxy) exhibits F115W-F200W > 1.9 mag (2-sigma lower limit) with a blue continuum slope, resulting in 99.6% of the photo-z PDF favoring z > 11. All data quality images show no artifacts at the candidate's position, and independent analyses consistently find a strong preference for z > 11. Its colors are inconsistent with Galactic stars, and it is resolved (r_h = 340 +/- 14 pc). Maisie's Galaxy has log M*/Msol ~ 8.5 and is highly star-forming (log sSFR ~ -8.2 yr^-1), with a blue rest-UV color (beta ~ -2.5) indicating little dust though not extremely low metallicity. While the presence of this source is in tension with most predictions, it agrees with empirical extrapolations assuming UV luminosity functions which smoothly decline with increasing redshift. Should followup spectroscopy validate this redshift, our Universe was already aglow with galaxies less than 400 Myr after the Big Bang. △ Less

Submitted 7 September, 2022; v1 submitted 25 July, 2022; originally announced July 2022.

Comments: 17 pages, 6 figures, 3 tables, ApJL in press. Summary of changes from original submission: Improvements in astrometry generated a weak detection in F150W that reduces the photo-z to 11.8 but does not increase the likelihood of lower-z solutions. A full discussion of changes from the original version is available at: https://web.corral.tacc.utexas.edu/ceersdata/papers/Maisie_update.pdf

The Science Performance of JWST as Characterized in Commissioning

Authors: Jane Rigby, Marshall Perrin, Michael McElwain, Randy Kimble, Scott Friedman, Matt Lallo, René Doyon, Lee Feinberg, Pierre Ferruit, Alistair Glasse, Marcia Rieke, George Rieke, Gillian Wright, Chris Willott, Knicole Colon, Stefanie Milam, Susan Neff, Christopher Stark, Jeff Valenti, Jim Abell, Faith Abney, Yasin Abul-Huda, D. Scott Acton, Evan Adams, David Adler , et al. (601 additional authors not shown)

Abstract: This paper characterizes the actual science performance of the James Webb Space Telescope (JWST), as determined from the six month commissioning period. We summarize the performance of the spacecraft, telescope, science instruments, and ground system, with an emphasis on differences from pre-launch expectations. Commissioning has made clear that JWST is fully capable of achieving the discoveries f… ▽ More This paper characterizes the actual science performance of the James Webb Space Telescope (JWST), as determined from the six month commissioning period. We summarize the performance of the spacecraft, telescope, science instruments, and ground system, with an emphasis on differences from pre-launch expectations. Commissioning has made clear that JWST is fully capable of achieving the discoveries for which it was built. Moreover, almost across the board, the science performance of JWST is better than expected; in most cases, JWST will go deeper faster than expected. The telescope and instrument suite have demonstrated the sensitivity, stability, image quality, and spectral range that are necessary to transform our understanding of the cosmos through observations spanning from near-earth asteroids to the most distant galaxies. △ Less

Submitted 10 April, 2023; v1 submitted 12 July, 2022; originally announced July 2022.

Comments: 5th version as accepted to PASP; 31 pages, 18 figures; https://iopscience.iop.org/article/10.1088/1538-3873/acb293

Journal ref: PASP 135 048001 (2023)

3D-DASH: The Widest Near-Infrared Hubble Space Telescope Survey

Authors: Lamiya A. Mowla, Sam E. Cutler, Gabriel B. Brammer, Ivelina G. Momcheva, Katherine E. Whitaker, Pieter G. van Dokkum, Rachel S. Bezanson, Natascha M. Forster Schreiber, Marijn Franx, Kartheik G. Iyer, Danilo Marchesini, Adam Muzzin, Erica J. Nelson, Rosalind E. Skelton, Gregory F. Snyder, David A. Wake, Stijn Wuyts, Arjen van der Wel

Abstract: The 3D-Drift And SHift (3D-DASH) program is a \textit{Hubble Space Telescope} WFC3 F160W imaging and G141 grism survey of the equatorial COSMOS field. 3D-DASH extends the legacy of HST near-infrared imaging and spectroscopy to degree-scale swaths of the sky, enabling the identification and study of distant galaxies ($z>2$) that are rare or in short-lived phases of galaxy evolution at rest-frame op… ▽ More The 3D-Drift And SHift (3D-DASH) program is a \textit{Hubble Space Telescope} WFC3 F160W imaging and G141 grism survey of the equatorial COSMOS field. 3D-DASH extends the legacy of HST near-infrared imaging and spectroscopy to degree-scale swaths of the sky, enabling the identification and study of distant galaxies ($z>2$) that are rare or in short-lived phases of galaxy evolution at rest-frame optical wavelengths. Furthermore, when combined with existing ACS/F814W imaging, the program facilitates spatially-resolved studies of the stellar populations and dust content of intermediate-redshift ($0.5<z<2$) galaxies. Here we present the reduced F160W imaging mosaic available to the community. Observed with the efficient DASH technique, the mosaic comprises 1256 individual WFC3 pointings, corresponding to an area of 1.35 deg$^2$ (1.43 deg$^2$ in 1912 when including archival data). The median $5σ$ point-source limit in $H_{160}$ is 24.74 mag. We also provide tools to determine the local point spread function (PSF), create cutouts, and explore the image at any location within the 3D-DASH footprint. 3D-DASH is the widest \textit{HST}/WFC3 imaging survey in the F160W filter to date, increasing the existing extragalactic survey area in the near-infrared at HST resolution by an order of magnitude. △ Less

Submitted 2 June, 2022; originally announced June 2022.

Comments: The 3D-DASH mosaic, PSF and cutout generator tool, and the image explorer are available at archive.stsci.edu/hlsp/3d-dash and at www.lamiyamowla.com/3d-dash. Accepted for publication in Astrophysical Journal. arXiv admin note: text overlap with arXiv:1808.04379

Mapping Thermoelectric Transport in a Multicomponent Alloy Space

Authors: Ramya Gurunathan, Suchismita Sarker, Christopher K. H. Borg, James Saal, Logan Ward, Apurva Mehta, G. Jeffrey Snyder

Abstract: Interest in high entropy alloy thermoelectric materials is predicated on achieving ultralow lattice thermal conductivity $κ\sub{L}$ through large compositional disorder. However, here we show that for a given mechanism, such as mass contrast phonon scattering, $κ\sub{L}$ will be minimized along the binary alloy with the highest mass contrast, such that adding an intermediate-mass atom to increase… ▽ More Interest in high entropy alloy thermoelectric materials is predicated on achieving ultralow lattice thermal conductivity $κ\sub{L}$ through large compositional disorder. However, here we show that for a given mechanism, such as mass contrast phonon scattering, $κ\sub{L}$ will be minimized along the binary alloy with the highest mass contrast, such that adding an intermediate-mass atom to increase atomic disorder can increase thermal conductivity. Only when each component adds an independent scattering mechanism (such as adding strain fluctuation to an existing mass fluctuation) is there a benefit. In addition, both charge carriers and heat-carrying phonons are known to experience scattering due to alloying effects, leading to a trade-off in thermoelectric performance. We apply analytic transport models, based on perturbation and effective medium theories, to predict how alloy scattering will affect the thermal and electronic transport across the full compositional range of several pseudo-ternary and pseudo-quaternary alloy systems. To do so, we demonstrate a multicomponent extension to both thermal and electronic binary alloy scattering models based on the virtual crystal approximation. Finally, we show that common functional forms used in computational thermodynamics can be applied to this problem to further generalize the scattering behavior that is modeled. △ Less

Submitted 3 May, 2022; originally announced May 2022.

Mode- and Space- Resolved Thermal Transport of Alloy Nanostructures

Authors: S. Aria Hosseini, Sarah Khanniche, G. Jeffrey Snyder, Samuel Huberman, P. Alex Greaney, Giuseppe Romano

Abstract: Nanostructured semiconducting alloys obtain ultra-low thermal conductivity as a result of the scattering of phonons with a wide range of mean-free-paths (MFPs). In these materials, long-MFP phonons are scattered at the nanoscale boundaries whereas short-MFP high-frequency phonons are impeded by disordered point defects introduced by alloying. While this trend has been validated by simplified analy… ▽ More Nanostructured semiconducting alloys obtain ultra-low thermal conductivity as a result of the scattering of phonons with a wide range of mean-free-paths (MFPs). In these materials, long-MFP phonons are scattered at the nanoscale boundaries whereas short-MFP high-frequency phonons are impeded by disordered point defects introduced by alloying. While this trend has been validated by simplified analytical and numerical methods, an ab-initio space-resolved approach remains elusive. To fill this gap, we calculate the thermal conductivity reduction in porous alloys by solving the mode-resolved Boltzmann transport equation for phonons using the finite-volume approach. We analyze different alloys, length-scales, concentrations, and temperatures, obtaining a very large reduction in the thermal conductivity over the entire configuration space. For example, a ~97% reduction is found for Al$_{0.8}$In$_{0.2}$As with 25% porosity. Furthermore, we employ these simulations to validate our recently introduced "Ballistic Correction Model" (BCM), an approach that estimates the effective thermal conductivity using the characteristic MFP of the bulk alloy and the length-scale of the material. The BCM is then used to provide guiding principles in designing alloy-based nanostructures. Notably, it elucidates how porous alloys such as Si$_{x}$Ge$_{1-x}$ obtain larger thermal conductivity reduction compared to porous Si or Ge, while also explaining why we should not expect similar behavior in alloys such as Al$_{x}$In$_{1-x}$As. By taking into account the synergy from scattering at different scales, we provide a route for the design of materials with ultra-low thermal conductivity. △ Less

Submitted 24 March, 2022; originally announced March 2022.

DeepAdversaries: Examining the Robustness of Deep Learning Models for Galaxy Morphology Classification

Authors: Aleksandra Ćiprijanović, Diana Kafkes, Gregory Snyder, F. Javier Sánchez, Gabriel Nathan Perdue, Kevin Pedro, Brian Nord, Sandeep Madireddy, Stefan M. Wild

Abstract: With increased adoption of supervised deep learning methods for processing and analysis of cosmological survey data, the assessment of data perturbation effects (that can naturally occur in the data processing and analysis pipelines) and the development of methods that increase model robustness are increasingly important. In the context of morphological classification of galaxies, we study the eff… ▽ More With increased adoption of supervised deep learning methods for processing and analysis of cosmological survey data, the assessment of data perturbation effects (that can naturally occur in the data processing and analysis pipelines) and the development of methods that increase model robustness are increasingly important. In the context of morphological classification of galaxies, we study the effects of perturbations in imaging data. In particular, we examine the consequences of using neural networks when training on baseline data and testing on perturbed data. We consider perturbations associated with two primary sources: 1) increased observational noise as represented by higher levels of Poisson noise and 2) data processing noise incurred by steps such as image compression or telescope errors as represented by one-pixel adversarial attacks. We also test the efficacy of domain adaptation techniques in mitigating the perturbation-driven errors. We use classification accuracy, latent space visualizations, and latent space distance to assess model robustness. Without domain adaptation, we find that processing pixel-level errors easily flip the classification into an incorrect class and that higher observational noise makes the model trained on low-noise data unable to classify galaxy morphologies. On the other hand, we show that training with domain adaptation improves model robustness and mitigates the effects of these perturbations, improving the classification accuracy by 23% on data with higher observational noise. Domain adaptation also increases by a factor of ~2.3 the latent space distance between the baseline and the incorrectly classified one-pixel perturbed image, making the model more robust to inadvertent perturbations. △ Less

Submitted 6 July, 2022; v1 submitted 28 December, 2021; originally announced December 2021.

Comments: 20 pages, 6 figures, 5 tables; accepted in MLST

Report number: FERMILAB-PUB-21-767-SCD

Robustness of deep learning algorithms in astronomy -- galaxy morphology studies

Authors: A. Ćiprijanović, D. Kafkes, G. N. Perdue, K. Pedro, G. Snyder, F. J. Sánchez, S. Madireddy, S. M. Wild, B. Nord

Abstract: Deep learning models are being increasingly adopted in wide array of scientific domains, especially to handle high-dimensionality and volume of the scientific data. However, these models tend to be brittle due to their complexity and overparametrization, especially to the inadvertent adversarial perturbations that can appear due to common image processing such as compression or blurring that are o… ▽ More Deep learning models are being increasingly adopted in wide array of scientific domains, especially to handle high-dimensionality and volume of the scientific data. However, these models tend to be brittle due to their complexity and overparametrization, especially to the inadvertent adversarial perturbations that can appear due to common image processing such as compression or blurring that are often seen with real scientific data. It is crucial to understand this brittleness and develop models robust to these adversarial perturbations. To this end, we study the effect of observational noise from the exposure time, as well as the worst case scenario of a one-pixel attack as a proxy for compression or telescope errors on performance of ResNet18 trained to distinguish between galaxies of different morphologies in LSST mock data. We also explore how domain adaptation techniques can help improve model robustness in case of this type of naturally occurring attacks and help scientists build more trustworthy and stable models. △ Less

Submitted 2 November, 2021; v1 submitted 1 November, 2021; originally announced November 2021.

Comments: Accepted in: Fourth Workshop on Machine Learning and the Physical Sciences (35th Conference on Neural Information Processing Systems; NeurIPS2021); final version

Report number: FERMILAB-CONF-21-561-SCD

Path Planning for Optimal Coverage of Areas with Nonuniform Importance

Authors: Gregory Snyder, Sachin Shriwastav, Dylan Morrison-Fogel, Zhuoyuan Song

Abstract: Coverage of an inaccessible or challenging region with potential health and safety hazards, such as in a volcanic region, is difficult yet crucial from scientific and meteorological perspectives. Areas contained within the region often provide valuable information of varying importance. We present an algorithm to optimally cover a volcanic region in Hawai`i with an unmanned aerial vehicle (UAV). T… ▽ More Coverage of an inaccessible or challenging region with potential health and safety hazards, such as in a volcanic region, is difficult yet crucial from scientific and meteorological perspectives. Areas contained within the region often provide valuable information of varying importance. We present an algorithm to optimally cover a volcanic region in Hawai`i with an unmanned aerial vehicle (UAV). The target region is assigned with a nonuniform coverage importance score distribution. For a specified battery capacity of the UAV, the optimization problem seeks the path that maximizes the total coverage area and the accumulated importance score while penalizing the revisiting of the same area. Trajectories are generated offline for the UAV based on the available power and coverage information map. The optimal trajectory minimizes the unspent battery power while enforcing that the UAV returns to its starting location. This multi-objective optimization problem is solved by using sequential quadratic programming. The details of the competitive optimization problem are discussed along with the analysis and simulation results to demonstrate the applicability of the proposed algorithm. △ Less

Submitted 10 December, 2021; v1 submitted 19 October, 2021; originally announced October 2021.

Comments: 10 pages, 7 figures

Dynamic Compressed Sensing of Unsteady Flows with a Mobile Robot

Authors: Sachin Shriwastav, Gregory Snyder, Zhuoyuan Song

Abstract: Large-scale environmental sensing with a finite number of mobile sensors is a challenging task that requires a lot of resources and time. This is especially true when features in the environment are spatiotemporally changing with unknown or partially known dynamics. Fortunately, these dynamic features often evolve in a low-dimensional space, making it possible to capture their dynamics sufficientl… ▽ More Large-scale environmental sensing with a finite number of mobile sensors is a challenging task that requires a lot of resources and time. This is especially true when features in the environment are spatiotemporally changing with unknown or partially known dynamics. Fortunately, these dynamic features often evolve in a low-dimensional space, making it possible to capture their dynamics sufficiently well with only one or several properly planned mobile sensors. This paper investigates the problem of dynamic compressed sensing of an unsteady flow field, which takes advantage of the inherently low dimensionality of the underlying flow dynamics to reduce number of waypoints for a mobile sensing robot. The optimal sensing waypoints are identified by an iterative compressed sensing algorithm that optimizes the flow reconstruction based on the proper orthogonal decomposition modes. An optimal sampling trajectory is then found to traverse these waypoints while minimizing the energy consumption, time, and flow reconstruction error. Simulation results in an unsteady double gyre flow field is presented to demonstrate the efficacy of the proposed algorithms. Experimental results with an indoor quadcopter are presented to show the feasibility of the resulting trajectory. △ Less

Submitted 2 March, 2022; v1 submitted 16 October, 2021; originally announced October 2021.

Comments: 8 pages, 7 figures

MSC Class: 49M37 (Primary) 68T40 (Secondary)

Koopman Operator Theory for Nonlinear Dynamic Modeling using Dynamic Mode Decomposition

Authors: Gregory Snyder, Zhuoyuan Song

Abstract: The Koopman operator is a linear operator that describes the evolution of scalar observables (i.e., measurement functions of the states) in an infinitedimensional Hilbert space. This operator theoretic point of view lifts the dynamics of a finite-dimensional nonlinear system to an infinite-dimensional function space where the evolution of the original system becomes linear. In this paper, we provi… ▽ More The Koopman operator is a linear operator that describes the evolution of scalar observables (i.e., measurement functions of the states) in an infinitedimensional Hilbert space. This operator theoretic point of view lifts the dynamics of a finite-dimensional nonlinear system to an infinite-dimensional function space where the evolution of the original system becomes linear. In this paper, we provide a brief summary of the Koopman operator theorem for nonlinear dynamics modeling and focus on analyzing several data-driven implementations using dynamical mode decomposition (DMD) for autonomous and controlled canonical problems. We apply the extended dynamic mode decomposition (EDMD) to identify the leading Koopman eigenfunctions and approximate a finite-dimensional representation of the discovered linear dynamics. This allows us to apply linear control approaches towards nonlinear systems without linearization approximations around fixed points. We can then examine the fidelity of using a linear controller based on a Koopman operator approximated system on under-actuated systems with basic maneuvers. We demonstrate the effectiveness of this theory through numerical simulation on two classic dynamical systems are used to show DMD methods of evaluating and approximating the Koopman operator and its effectiveness at linearizing these systems. △ Less

Submitted 15 October, 2021; originally announced October 2021.

Comments: 8 pages, 16 figures

The mass-metallicity relation at z~1-2 and its dependence on star formation rate

Authors: Alaina Henry, Marc Rafelski, Ben Sunnquist, Norbert Pirzkal, Camilla Pacifici, Hakim Atek, Micaela Bagley, Ivano Baronchelli, Guillermo Barro, Andrew J. Bunker, James Colbert, Y. Sophia Dai, Bruce G. Elmegreen, Debra Meloy Elmegreen, Steven Finkelstein, Dale Kocevski, Anton Koekemoer, Matthew Malkan, Crystal L. Martin, Vihang Mehta, Anthony Pahl, Casey Papovich, Michael Rutkowski, Jorge Sanchez Almeida, Claudia Scarlata , et al. (2 additional authors not shown)

Abstract: We present a new measurement of the gas-phase mass-metallicity relation (MZR), and its dependence on star formation rates (SFRs) at 1.3 < z < 2.3. Our sample comprises 1056 galaxies with a mean redshift of z = 1.9, identified from the Hubble Space Telescope Wide Field Camera 3 (WFC3) grism spectroscopy in the Cosmic Assembly Near-Infrared Deep Extragalactic Survey (CANDELS) and the WFC3 Infrared S… ▽ More We present a new measurement of the gas-phase mass-metallicity relation (MZR), and its dependence on star formation rates (SFRs) at 1.3 < z < 2.3. Our sample comprises 1056 galaxies with a mean redshift of z = 1.9, identified from the Hubble Space Telescope Wide Field Camera 3 (WFC3) grism spectroscopy in the Cosmic Assembly Near-Infrared Deep Extragalactic Survey (CANDELS) and the WFC3 Infrared Spectroscopic Parallel Survey (WISP). This sample is four times larger than previous metallicity surveys at z ~ 2, and reaches an order of magnitude lower in stellar mass (10^8 M_sun). Using stacked spectra, we find that the MZR evolves by 0.3 dex relative to z ~ 0.1. Additionally, we identify a subset of 49 galaxies with high signal-to-noise (SNR) spectra and redshifts between 1.3 < z < 1.5, where H-alpha emission is observed along with [OIII] and [OII]. With accurate measurements of SFR in these objects, we confirm the existence of a mass-metallicity-SFR (M-Z-SFR) relation at high redshifts. These galaxies show systematic differences from the local M-Z-SFR relation, which vary depending on the adopted measurement of the local relation. However, it remains difficult to ascertain whether these differences could be due to redshift evolution, as the local M-Z-SFR relation is poorly constrained at the masses and SFRs of our sample. Lastly, we reproduced our sample selection in the IllustrisTNG hydrodynamical simulation, demonstrating that our line flux limit lowers the normalization of the simulated MZR by 0.2 dex. We show that the M-Z-SFR relation in IllustrisTNG has an SFR dependence that is too steep by a factor of around three. △ Less

Submitted 1 July, 2021; originally announced July 2021.

Comments: Accepted for publication in ApJ; 41 pages, 20 figures

Dopant-segregation to grain boundaries controls electrical conductivity of n-type NbCo(Pt)Sn half-Heusler alloy mediating thermoelectric performance

Authors: Ting Luo, Federico Serrano-Sánchez, Hanna Bishara, Siyuan Zhang, Ruben Bueno Villoro, Jimmy Jiahong Kuo, Claudia Felser, Christina Scheu, G. Jeffrey Snyder, James P. Best, Gerhard Dehm, Yuan Yu, Dierk Raabe, Chenguang Fu, Baptiste Gault

Abstract: Science-driven design of future thermoelectric materials requires a deep understanding of the fundamental relationships between microstructure and transport properties. Grain boundaries in polycrystalline materials influence the thermoelectric performance through the scattering of phonons or the trapping of electrons due to space-charge effects. Yet, the current lack of careful investigations on g… ▽ More Science-driven design of future thermoelectric materials requires a deep understanding of the fundamental relationships between microstructure and transport properties. Grain boundaries in polycrystalline materials influence the thermoelectric performance through the scattering of phonons or the trapping of electrons due to space-charge effects. Yet, the current lack of careful investigations on grain boundary-associated features hinders further optimization of properties. Here, we study n-type NbCo1-xPtxSn half-Heusler alloys, which were synthesized by ball milling and spark plasma sintering (SPS). Post-SPS annealing was performed on one sample, leading to improved low-temperature electrical conductivity. The microstructure of both samples was examined by electron microscopy and atom probe tomography. The grain size increases from ~230 nm to ~2.38 μm upon annealing. Pt is found within grains and at grain boundaries, where it locally reduces the resistivity, as assessed by in situ four-point-probe electrical conductivity measurement. Our work showcases the correlation between microstructure and electrical conductivity, providing opportunities for future microstructural optimization by tuning the chemical composition at grain boundaries. △ Less

Submitted 1 July, 2021; originally announced July 2021.

A Census of the Bright z=8.5-11 Universe with the Hubble and Spitzer Space Telescopes in the CANDELS Fields

Authors: Steven L. Finkelstein, Micaela Bagley, Mimi Song, Rebecca Larson, Casey Papovich, Mark Dickinson, Keely Finkelstein, Anton M. Koekemoer, Norbert Pirzkal, Rachel S. Somerville, L. Y. Aaron Yung, Peter Behroozi, Harry Ferguson, Mauro Giavalisco, Norman Grogin, Nimish Hathi, Taylor Hutchison, Intae Jung, Dale Kocevski, Lalitwadee Kawinwanichakij, Sofia Rojas-Ruiz, Russell Ryan Jr., Gregory F. Snyder, Sandro Tacchella

Abstract: We present the results from a new search for candidate galaxies at z ~ 8.5-11 discovered over the 850 arcmin^2 area probed by the Cosmic Assembly Near-Infrared Deep Extragalactic Legacy Survey (CANDELS). We use a photometric redshift selection including both Hubble and Spitzer Space Telescope photometry to robustly identify galaxies in this epoch at F160W < 26.6. We use a detailed vetting procedur… ▽ More We present the results from a new search for candidate galaxies at z ~ 8.5-11 discovered over the 850 arcmin^2 area probed by the Cosmic Assembly Near-Infrared Deep Extragalactic Legacy Survey (CANDELS). We use a photometric redshift selection including both Hubble and Spitzer Space Telescope photometry to robustly identify galaxies in this epoch at F160W < 26.6. We use a detailed vetting procedure, including screening for persistence, stellar contamination, inclusion of ground-based imaging, and followup space-based imaging to build a robust sample of 11 candidate galaxies, three presented here for the first time. The inclusion of Spitzer/IRAC photometry in the selection process reduces contamination, and yields more robust redshift estimates than Hubble alone. We constrain the evolution of the rest-frame ultraviolet luminosity function via a new method of calculating the observed number densities without choosing a prior magnitude bin size. We find that the abundance at our brightest probed luminosities (M_UV=-22.3) is consistent with predictions from simulations which assume that galaxies in this epoch have gas depletion times at least as short as those in nearby starburst galaxies. Due to large Poisson and cosmic variance uncertainties we cannot conclusively rule out either a smooth evolution of the luminosity function continued from z=4-8, or an accelerate decline at z > 8. We calculate that the presence of seven galaxies in a single field (EGS) is an outlier at the 2-sigma significance level, implying the discovery of a significant overdensity. These scenarios will be imminently testable to high confidence within the first year of observations of the James Webb Space Telescope. △ Less

Submitted 19 November, 2021; v1 submitted 25 June, 2021; originally announced June 2021.

Comments: 47 figures, 25 pages, 10 tables, Accepted to the Astrophysical Journal. Photometric catalogs available via email request to lead author

Are all post-starbursts mergers? HST reveals hidden disturbances in the majority of PSBs

Authors: Elizaveta Sazonova, Katherine Alatalo, Kate Rowlands, Susana E. Deustua, Decker French, Timothy M. Heckman, Lauranne Lanz, Ute Lisenfeld, Yuanze Luo, Anne M. Medling, Kristina Nyland, Justin A. Otter, Andreea Petric, Gregory F. Snyder, Claudia M. Urry

Abstract: How do galaxies transform from blue, star-forming spirals to red, quiescent early-type galaxies? To answer this question, we analyzed a set of 26 gas-rich, shocked post-starburst galaxies with Hubble Space Telescope (HST) imaging in B, I, and H bands, and Sloan Digital Sky Survey (SDSS) i-band imaging of similar depth but lower resolution. We found that post-starbursts in our sample have intermedi… ▽ More How do galaxies transform from blue, star-forming spirals to red, quiescent early-type galaxies? To answer this question, we analyzed a set of 26 gas-rich, shocked post-starburst galaxies with Hubble Space Telescope (HST) imaging in B, I, and H bands, and Sloan Digital Sky Survey (SDSS) i-band imaging of similar depth but lower resolution. We found that post-starbursts in our sample have intermediate morphologies between disk- and bulge-dominated (Sérsic n$=1.7^{+0.3}_{-0.0}$) and have red bulges, likely due to dust obscuration in the cores. Majority of galaxies in our sample are more morphologically disturbed than regular galaxies (88%, corresponding to >3$σ$ significance) when observed with HST, with asymmetry and Sérsic residual flux fraction being the most successful measures of disturbance. Most disturbances are undetected at the lower resolution of SDSS imaging. Although ~27% galaxies are clear merger remnants, we found that disturbances in another ~30% of the sample are internal, caused by small-scale perturbations or dust substructures rather than tidal features, and require high-resolution imaging to detect. We found a 2.8$σ$ evidence that asymmetry features fade on timescales ~200 Myr, and may vanish entirely after ~750 Myr, so we do not rule out a possible merger origin of all post-starbursts given that asymmetric features may have already faded. This work highlights the importance of small-scale disturbances, detected only in high-resolution imaging, in understanding structural evolution of transitioning galaxies. △ Less

Submitted 20 May, 2021; originally announced May 2021.

Comments: 31 pages, 19 figures, 2 figure sets, 1 machine-readable table; accepted to ApJ

Comparison of the Tetrahedron Method to Smearing Methods for the Electronic Density of States

Authors: M. Y. Toriyama, A. M. Ganose, M. Dylla, S. Anand, J. Park, M. K. Brod, J. Munro, K. A. Persson, A. Jain, G. J. Snyder

Abstract: The electronic density of states (DOS) highlights fundamental properties of materials that oftentimes dictate their properties, such as the band gap and Van Hove singularities. In this short note, we discuss how sharp features of the density of states can be obscured by smearing methods (such as the Gaussian and Fermi smearing methods) when calculating the DOS. While the common approach to reach a… ▽ More The electronic density of states (DOS) highlights fundamental properties of materials that oftentimes dictate their properties, such as the band gap and Van Hove singularities. In this short note, we discuss how sharp features of the density of states can be obscured by smearing methods (such as the Gaussian and Fermi smearing methods) when calculating the DOS. While the common approach to reach a "converged" density of states of a material is to increase the discrete k-point mesh density, we show that the DOS calculated by smearing methods can appear to converge but not to the correct DOS. Employing the tetrahedron method for Brillouin zone integration resolves key features of the density of states far better than smearing methods. △ Less

Submitted 4 March, 2021; originally announced March 2021.

DeepMerge II: Building Robust Deep Learning Algorithms for Merging Galaxy Identification Across Domains

Authors: A. Ćiprijanović, D. Kafkes, K. Downey, S. Jenkins, G. N. Perdue, S. Madireddy, T. Johnston, G. F. Snyder, B. Nord

Abstract: In astronomy, neural networks are often trained on simulation data with the prospect of being used on telescope observations. Unfortunately, training a model on simulation data and then applying it to instrument data leads to a substantial and potentially even detrimental decrease in model accuracy on the new target dataset. Simulated and instrument data represent different data domains, and for a… ▽ More In astronomy, neural networks are often trained on simulation data with the prospect of being used on telescope observations. Unfortunately, training a model on simulation data and then applying it to instrument data leads to a substantial and potentially even detrimental decrease in model accuracy on the new target dataset. Simulated and instrument data represent different data domains, and for an algorithm to work in both, domain-invariant learning is necessary. Here we employ domain adaptation techniques$-$ Maximum Mean Discrepancy (MMD) as an additional transfer loss and Domain Adversarial Neural Networks (DANNs)$-$ and demonstrate their viability to extract domain-invariant features within the astronomical context of classifying merging and non-merging galaxies. Additionally, we explore the use of Fisher loss and entropy minimization to enforce better in-domain class discriminability. We show that the addition of each domain adaptation technique improves the performance of a classifier when compared to conventional deep learning algorithms. We demonstrate this on two examples: between two Illustris-1 simulated datasets of distant merging galaxies, and between Illustris-1 simulated data of nearby merging galaxies and observed data from the Sloan Digital Sky Survey. The use of domain adaptation techniques in our experiments leads to an increase of target domain classification accuracy of up to ${\sim}20\%$. With further development, these techniques will allow astronomers to successfully implement neural network models trained on simulation data to efficiently detect and study astrophysical objects in current and future large-scale astronomical surveys. △ Less

Submitted 1 March, 2021; originally announced March 2021.

Comments: Submitted to MNRAS; 21 pages, 9 figures, 9 tables

Report number: FERMILAB-PUB-21-072-SCD

Journal ref: MNRAS, Volume 506, Issue 1, September 2021, Page 677

Phonon Scattering in the Complex Strain Field of a Dislocation

Authors: Yandong Sun, Yanguang Zhou, Ramya Gurunathan, Jin-Yu Zhang, Ming Hu, Wei Liu, Ben Xu, G. Jeffrey Snyder

Abstract: Strain engineering is critical to the performance enhancement of electronic and thermoelectric devices because of its influence on the material thermal conductivity. However, current experiments cannot probe the detailed physics of the phonon-strain interaction due to the complex, inhomogeneous, and long-distance features of the strain field in real materials. Dislocations provide us with an excel… ▽ More Strain engineering is critical to the performance enhancement of electronic and thermoelectric devices because of its influence on the material thermal conductivity. However, current experiments cannot probe the detailed physics of the phonon-strain interaction due to the complex, inhomogeneous, and long-distance features of the strain field in real materials. Dislocations provide us with an excellent model to investigate these inhomogeneous strain fields. In this study, non-equilibrium molecular dynamics simulations were used to study the lattice thermal conductivity of PbTe under different strain status tuned by dislocation densities. The extended 1D McKelvey-Shockley flux method was used to analyze the frequency dependence of phonon scattering in the inhomogeneously strained regions of dislocations. A spatially resolved phonon dislocation scattering process was shown, where the unequal strain in different regions affected the magnitude and frequency-dependence of the scattering rate. Our study not only advances the knowledge of strain scattering of phonon propagation but offers fundamental guidance on optimizing thermal management by structure design. △ Less

Submitted 26 January, 2021; originally announced January 2021.

Comments: 31 pages,13 figures

The Connection between Mergers and AGN Activity in Simulated and Observed Massive Galaxies

Authors: Ray S. Sharma, Ena Choi, Rachel S. Somerville, Gregory F. Snyder, Dale D. Kocevski, Michaela Hirschmann, Benjamin P. Moster, Thorsten Naab, Desika Narayanan, Jeremiah P. Ostriker, David J. Rosario

Abstract: We analyze a suite of $30$ high resolution zoom-in cosmological hydrodynamic simulations of massive galaxies with stellar masses $M_{\ast} > 10^{10.9} M_\odot$, with the goal of better understanding merger activity in AGN, AGN activity in merging systems, SMBH growth during mergers, and the role of gas content. Using the radiative transfer code \textsc{Powderday}, we generate HST-WFC3 F160W synthe… ▽ More We analyze a suite of $30$ high resolution zoom-in cosmological hydrodynamic simulations of massive galaxies with stellar masses $M_{\ast} > 10^{10.9} M_\odot$, with the goal of better understanding merger activity in AGN, AGN activity in merging systems, SMBH growth during mergers, and the role of gas content. Using the radiative transfer code \textsc{Powderday}, we generate HST-WFC3 F160W synthetic observations of redshift $0.5 < z < 3$ central galaxies, add noise properties similar to the CANDELS survey, and measure morphological properties from the synthetic images using commonly adopted non-parametric statistics. We compare the distributions of morphological properties measured from the synthetic images with a sample of inactive galaxies and X-ray selected AGN hosts from CANDELS. We study the connection between mergers and AGN activity in the simulations, the synthetic images, and the observed CANDELS sample. We find that, in both the simulations and CANDELS, even the most luminous $(L_{\rm bol} > 10^{45}$ erg s$^{-1})$ AGN in our sample are no more likely than inactive galaxies $(L_{\rm bol} < 10^{43}$ erg s$^{-1})$ to be found in merging systems. We also find that AGN activity is not overall enhanced by mergers, nor enhanced at any specific time in the $1$ Gyr preceding and following a merger. Even gas rich major mergers (stellar mass ratio $>$1:4) do not necessarily enhance AGN activity or significantly grow the central SMBH. We conclude that in the simulated massive galaxies studied here, mergers are not the primary drivers of AGN. △ Less

Submitted 18 December, 2023; v1 submitted 5 January, 2021; originally announced January 2021.

Thermal Resistance at a Twist Boundary and Semicoherent Heterointerface

Authors: Ramya Gurunathan, Riley Hanus, Samuel Graham, Anupam Garg, G. Jeffrey Snyder

Abstract: Traditional models of interfacial phonon scattering, including the acoustic mismatch model (AMM) and diffuse mismatch model (DMM), take into account the bulk properties of the material surrounding the interface, but not the atomic structure and properties of the interface itself. Here, we derive a theoretical formalism for the phonon scattering at a dislocation grid, or two interpenetrating orthog… ▽ More Traditional models of interfacial phonon scattering, including the acoustic mismatch model (AMM) and diffuse mismatch model (DMM), take into account the bulk properties of the material surrounding the interface, but not the atomic structure and properties of the interface itself. Here, we derive a theoretical formalism for the phonon scattering at a dislocation grid, or two interpenetrating orthogonal arrays of dislocations, as this is the most stable structure of both the symmetric twist boundary and semicoherent heterointerface. With this approach, we are able to separately examine the contribution to thermal resistance due to the step function change in acoustic properties and due to interfacial dislocation strain fields, which induces diffractive scattering. Both low-angle Si-Si twist boundaries and the Si-Ge heterointerfaces are considered here and compared to previous experimental and simulation results. This work indicates that scattering from misfit dislocation strain fields doubles the thermal boundary resistance of Si-Ge heterointerfaces compared to scattering due to acoustic mismatch alone. Scattering from grain boundary dislocation strain fields is predicted to dominate the thermal boundary resistance of Si-Si twist boundaries. This physical treatment can guide the thermal design of devices by quantifying the relative importance of interfacial strain fields, which can be engineered via fabrication and processing methods, versus acoustic mismatch, which is fixed for a given interface. Additionally, this approach captures experimental and simulation trends such as the dependence of thermal boundary resistance on the grain boundary angle and interfacial strain energy. △ Less

Submitted 13 April, 2021; v1 submitted 4 January, 2021; originally announced January 2021.

Comments: 11 figures

Journal ref: Phys. Rev. B 103, 144302 (2021)

When Band Convergence is Not Beneficial for Thermoelectrics

Authors: Junsoo Park, Max Dylla, Yi Xia, Max Wood, G. Jeffrey Snyder, Anubhav Jain

Abstract: Band convergence is considered a clear benefit to thermoelectric performance because it increases the charge carrier concentration for a given Fermi level, which typically enhances charge conductivity while preserving the Seebeck coefficient. However, this advantage hinges on the assumption that interband scattering of carriers is weak or insignificant. With first-principles treatment of electron-… ▽ More Band convergence is considered a clear benefit to thermoelectric performance because it increases the charge carrier concentration for a given Fermi level, which typically enhances charge conductivity while preserving the Seebeck coefficient. However, this advantage hinges on the assumption that interband scattering of carriers is weak or insignificant. With first-principles treatment of electron-phonon scattering in CaMg$_{2}$Sb$_{2}$-CaZn$_{2}$Sb$_{2}$ Zintl system and full Heusler Sr$_{2}$SbAu, we demonstrate that the benefit of band convergence can be intrinsically negated by interband scattering depending on the manner in which bands converge. In the Zintl alloy, band convergence does not improve weighted mobility or the density-of-states effective mass. We trace the underlying reason to the fact that the bands converge at one k-point, which induces strong interband scattering of both the deformation-potential and the polar-optical kinds. The case contrasts with band convergence at distant k-points (as in the full Heusler), which better preserves the single-band scattering behavior thereby successfully leading to improved performance. Therefore, we suggest that band convergence as thermoelectric design principle is best suited to cases in which it occurs at distant k-points. △ Less

Submitted 26 May, 2021; v1 submitted 3 December, 2020; originally announced December 2020.

Morphological evolution of supermassive black hole merger hosts and multimessenger signatures

Authors: Colin DeGraf, Debora Sijacki, Tiziana Di Matteo, Kelly Holley-Bockelmann, Greg Snyder, Volker Springel

Abstract: With projects such as Laser Interferometer Space Antenna (LISA) and Pulsar Timing Arrays expected to detect gravitational waves from supermassive black hole mergers in the near future, it is key that we understand what we expect those detections to be, and maximize what we can learn from them. To address this, we study the mergers of supermassive black holes in the Illustris simulation, the overal… ▽ More With projects such as Laser Interferometer Space Antenna (LISA) and Pulsar Timing Arrays expected to detect gravitational waves from supermassive black hole mergers in the near future, it is key that we understand what we expect those detections to be, and maximize what we can learn from them. To address this, we study the mergers of supermassive black holes in the Illustris simulation, the overall rate of mergers, and the correlation between merging black holes and their host galaxies. We find that these mergers occur in typical galaxies along the $M_{\rm{BH}}-M_*$ relation, and that between LISA and PTAs we expect to probe the full range of galaxy masses. As galaxy mergers can trigger increased star formation, we find that galaxies hosting low-mass black hole mergers tend to show a slight increase in star formation rates compared to a mass-matched sample. However, high-mass merger hosts have typical star formation rates, due to a combination of low gas fractions and powerful AGN feedback. Although minor black hole mergers do not correlate with disturbed morphologies, major mergers (especially at high-masses) tend to show morphological evidence of recent galaxy mergers which survives for ~500 Myr. This is on the same scale as the infall/hardening time of the merging black holes, suggesting that electromagnetic followups to gravitational wave signals may not be able to observe this correlation. We further find that incorporating a realistic timescale delay for the black hole mergers could shift the distribution of merger masses toward higher-masses, decreasing the rate of LISA detections while increasing the rate of PTA detections. △ Less

Submitted 28 July, 2021; v1 submitted 1 December, 2020; originally announced December 2020.

Comments: 14 pages, 17 figures. Published in MNRAS

Journal ref: MNRAS (2021), 503, 3629D

The nature of giant clumps in high-z discs: a deep-learning comparison of simulations and observations

Authors: Omri Ginzburg, Marc Huertas-Company, Avishai Dekel, Nir Mandelker, Gregory Snyder, Daniel Ceverino, Joel Primack

Abstract: We use deep learning to explore the nature of observed giant clumps in high-redshift disc galaxies, based on their identification and classification in cosmological simulations. Simulated clumps are detected using the 3D gas and stellar densities in the VELA zoom-in cosmological simulation suite, with $\sim \!\! 25\ \!\rm{pc}$ maximum resolution, targeting main sequence galaxies at… ▽ More We use deep learning to explore the nature of observed giant clumps in high-redshift disc galaxies, based on their identification and classification in cosmological simulations. Simulated clumps are detected using the 3D gas and stellar densities in the VELA zoom-in cosmological simulation suite, with $\sim \!\! 25\ \!\rm{pc}$ maximum resolution, targeting main sequence galaxies at $1\!<\!z\!<\!3$. The clumps are classified as long-lived clumps (LLCs) or short-lived clumps (SLCs) based on their longevity in the simulations. We then train neural networks to detect and classify the simulated clumps in mock, multi-color, dusty and noisy HST-like images. The clumps are detected using an encoder-decoder convolutional neural network (CNN), and are classified according to their longevity using a vanilla CNN. Tests using the simulations show our detector and classifier to be $\sim80\%$ complete and $\sim80\%$ pure for clumps more massive than $\sim10^{7.5}\rm{M_\odot}$. When applied to observed galaxies in the CANDELS/GOODS S+N fields, we find both types of clumps to appear in similar abundances in the simulations and the observations. LLCs are, on average, more massive than SLCs by $\sim 0.5\ \rm{dex}$, and they dominate the clump population above $M_{\rm c}\gtrsim 10^{7.6}\ \rm{M_\odot}$. LLCs tend to be found closer to the galactic centre, indicating clump migration to the centre or preferential formation at smaller radii. The LLCs are found to reside in high mass galaxies, indicating better clump survivability under supernova feedback there, due to clumps being more massive in these galaxies. We find the clump properties in the simulations and the observations to agree within a factor of $\sim\! 2$. △ Less

Submitted 1 November, 2021; v1 submitted 12 November, 2020; originally announced November 2020.

Comments: 17 pages, 13 figures. Accepted for publication in MNRAS

Probing the Phonon Mean Free Paths in Dislocation Core by Molecular Dynamics Simulation

Authors: Yandong Sun, Yanguang Zhou, Ming Hu, G. Jeffrey Snyder, Ben Xu, Wei Liu

Abstract: Thermal management is extremely important for designing high-performance devices. The lattice thermal conductivity of materials is strongly dependent on the structural defects at different length scales, particularly point defects like vacancies, line defects like dislocations, and planar defects such as grain boundaries. Traditionally, the McKelvey-Shockley phonon Boltzmann's transport equation (… ▽ More Thermal management is extremely important for designing high-performance devices. The lattice thermal conductivity of materials is strongly dependent on the structural defects at different length scales, particularly point defects like vacancies, line defects like dislocations, and planar defects such as grain boundaries. Traditionally, the McKelvey-Shockley phonon Boltzmann's transport equation (BTE) method combined with molecular dynamics simulations has been widely used to evaluate the phonon mean free paths (MFPs) in defective systems. However, this method can only provide the aggregate MFPs of the whole sample. It is, therefore, challenging to extract the MFPs in the different regions with different thermal properties. In this study, the 1D McKelvey-Shockley phonon BTE method was extended to model inhomogeneous materials, where the effect of defects on the phonon MFPs is explicitly obtained. Then, the method was used to study the phonon interactions with the core structure of an edge dislocation. The phonon MFPs in the dislocation core were obtained and consistent with the analytical model such that high frequency phonons are likely to be scattered in this area. This method not only advances the knowledge of phonon-dislocation scattering but also shows the potential to investigate phonon transport behaviors in more complicated materials. △ Less

Submitted 26 January, 2021; v1 submitted 17 October, 2020; originally announced October 2020.

Comments: 21pages, 8 figures

MSC Class: 80A05 ACM Class: I.6.0

Investigating the Effect of Galaxy Interactions on AGN Enhancement at $0.5<z<3.0$

Authors: Ekta A. Shah, Jeyhan S. Kartaltepe, Christina T. Magagnoli, Isabella G. Cox, Caleb T. Wetherell, Brittany N. Vanderhoof, Antonello Calabro, Nima Chartab, Christopher J. Conselice, Darren J. Croton, Jennifer Donley, Laura de Groot, Alexander de la Vega, Nimish P. Hathi, Olivier Ilbert, Hanae Inami, Dale D. Kocevski, Anton M. Koekemoer, Brian C. Lemaux, Kameswara Bharadwaj Mantha, Stefano Marchesi, Marie Martig, Daniel C. Masters, Elizabeth J. McGrath, Daniel H. McIntosh , et al. (8 additional authors not shown)

Abstract: Galaxy interactions and mergers are thought to play an important role in the evolution of galaxies. Studies in the nearby universe show a higher AGN fraction in interacting and merging galaxies than their isolated counterparts, indicating that such interactions are important contributors to black hole growth. To investigate the evolution of this role at higher redshifts, we have compiled the large… ▽ More Galaxy interactions and mergers are thought to play an important role in the evolution of galaxies. Studies in the nearby universe show a higher AGN fraction in interacting and merging galaxies than their isolated counterparts, indicating that such interactions are important contributors to black hole growth. To investigate the evolution of this role at higher redshifts, we have compiled the largest known sample of major spectroscopic galaxy pairs (2381 with $ΔV <5000$ km s$^{-1}$) at $0.5<z<3.0$ from observations in the COSMOS and CANDELS surveys. We identify X-ray and IR AGN among this kinematic pair sample, a visually identified sample of mergers and interactions, and a mass-, redshift-, and environment-matched control sample for each in order to calculate AGN fractions and the level of AGN enhancement as a function of relative velocity, redshift, and X-ray luminosity. While we see a slight increase in AGN fraction with decreasing projected separation, overall, we find no significant enhancement relative to the control sample at any separation. In the closest projected separation bin ($<25$ kpc, $ΔV <1000$ km s$^{-1}$), we find enhancements of a factor of 0.94$^{+0.21}_{-0.16}$ and 1.00$^{+0.58}_{-0.31}$ for X-ray and IR-selected AGN, respectively. While we conclude that galaxy interactions do not significantly enhance AGN activity on average over $0.5<z<3.0$ at these separations, given the errors and the small sample size at the closest projected separations, our results would be consistent with the presence of low-level AGN enhancement. △ Less

Submitted 6 October, 2020; originally announced October 2020.

Comments: Accepted for publication in ApJ, 28 pages, 15 figures

The Morphology-Density relationship in 1<z<2 clusters

Authors: Elizaveta Sazonova, Katherine Alatalo, Jennifer Lotz, Kate Rowlands, Gregory F. Snyder, Kyle Boone, Mark Brodwin, Brian Hayden, Lauranne Lanz, Saul Perlmutter, Vicente Rodriguez-Gomez

Abstract: The morphology-density relationship states that dense cosmic environments such as galaxy clusters have an overabundance of quiescent elliptical galaxies, but it is unclear at which redshift this relationship is first established. We study the morphology of 4 clusters with $1.2<z<1.8$ using HST imaging and the morphology computation code statmorph. By comparing median morphology of cluster galaxies… ▽ More The morphology-density relationship states that dense cosmic environments such as galaxy clusters have an overabundance of quiescent elliptical galaxies, but it is unclear at which redshift this relationship is first established. We study the morphology of 4 clusters with $1.2<z<1.8$ using HST imaging and the morphology computation code statmorph. By comparing median morphology of cluster galaxies to CANDELS field galaxies using Monte Carlo analysis, we find that 2 out of 4 clusters (at z=1.19 and z=1.75) have an established morphology-density relationship with more than $3σ$ significance. $\sim$50% of galaxies in these clusters are bulge-dominated compared to $\sim$30% in the field, and they are significantly more compact. This result is more significant for low-mass galaxies with $\log M/M_\odot \lessapprox 10.5$, showing that low-mass galaxies are affected the most in clusters. We also find an intriguing system of two z $\approx$ 1.45 clusters at a unusually small separation 2D separation of $3'$ and 3D separation of $\approx73$ Mpc that exhibit no morphology-density relationship but have enhanced merger signatures. We conclude that the environmental mechanism responsible for the morphology-density relationship is 1) already active as early as z=1.75, 2) forms compact, bulge-dominated galaxies and 3) affects primarily low-mass galaxies. However, there is a significant degree of intracluster variance that may depend on the larger cosmological environment in which the cluster is embedded. △ Less

Submitted 9 July, 2020; v1 submitted 7 July, 2020; originally announced July 2020.

Comments: 27 pages, 14 figures, Accepted for publication in ApJ

Stellar Masses of Giant Clumps in CANDELS and Simulated Galaxies Using Machine Learning

Authors: M. Huertas-Company, Y. Guo, O. Ginzburg, C. T. Lee, N. Mandelker, M. Metter, J. R. Primack, A. Dekel, D. Ceverino, S. M. Faber, D. C. Koo, A. Koekemoer, G. Snyder, M. Giavalisco, H. Zhang

Abstract: A significant fraction of high redshift star-forming disc galaxies are known to host giant clumps, whose nature and role in galaxy evolution are yet to be understood. In this work we first present a new method based on neural networks to detect clumps in galaxy images. We use this method to detect clumps in the rest-frame optical and UV images of a complete sample of $\sim1500$ star forming galaxi… ▽ More A significant fraction of high redshift star-forming disc galaxies are known to host giant clumps, whose nature and role in galaxy evolution are yet to be understood. In this work we first present a new method based on neural networks to detect clumps in galaxy images. We use this method to detect clumps in the rest-frame optical and UV images of a complete sample of $\sim1500$ star forming galaxies at $1<z<3$ in the CANDELS survey as well as in images from the VELA zoom-in cosmological simulations. We show that observational effects have a dramatic impact on the derived clump properties leading to an overestimation of the clump mass up to a factor of 10, which highlights the importance of fair comparisons between observations and simulations and the limitations of current HST data to study the resolved structure of distant galaxies. After correcting for these effects with a mixture density network, we estimate that the clump stellar mass function follows a power-law down to the completeness limit ($10^{7}$ solar masses) with the majority of the clumps being less massive than $10^9$ solar masses. This is in better agreement with recent gravitational lensing based measurements. The simulations explored in this work overall reproduce the shape of the observed clump stellar mass function and clumpy fractions when confronted under the same conditions, although they tend to lie in the lower limit of the confidence intervals of the observations. This agreement suggests that most of the observed clumps are formed in-situ. △ Less

Submitted 8 September, 2020; v1 submitted 25 June, 2020; originally announced June 2020.

Comments: Accepted for publication in MNRAS - This is the final version

Powderday: Dust Radiative Transfer for Galaxy Simulations

Authors: Desika Narayanan, Matthew J. Turk, Thomas Robitaille, Ashley J. Kelly, B. Connor McClellan, Ray S. Sharma, Prerak Garg, Matthew Abruzzo, Ena Choi, Charlie Conroy, Benjamin D. Johnson, Benjamin Kimock, Qi Li, Christopher C. Lovell, Sidney Lower, George C. Privon, Jonathan Roberts, Snigdaa Sethuram, Gregory F. Snyder, Robert Thompson, John H. Wise

Abstract: We present Powderday, a flexible, fast, open-source dust radiative transfer package designed to interface with galaxy formation simulations. Powderday builds on FSPS population synthesis models, Hyperion dust radiative transfer, and employs yt to interface between different software packages. We include our stellar population synthesis modeling on the fly, which allows for significant run-time fle… ▽ More We present Powderday, a flexible, fast, open-source dust radiative transfer package designed to interface with galaxy formation simulations. Powderday builds on FSPS population synthesis models, Hyperion dust radiative transfer, and employs yt to interface between different software packages. We include our stellar population synthesis modeling on the fly, which allows for significant run-time flexibility in the assumed stellar physics. We include a model for nebular line emission that can employ either precomputed Cloudy lookup tables (for efficiency), or direct photoionization calculations for all young stars (for flexibility). The dust content follows either observationally-motivated prescriptions, direct modeling from galaxy formation simulations, or a novel approach that includes the dust content via learning-based algorithms from the SIMBA cosmological galaxy formation simulation. AGN can additionally be included via a range of prescriptions. The output of these models are broadband SEDs, as well as filter-convolved images. Powderday is designed to eliminate last-mile efforts by researchers that employ different hydrodynamic galaxy formation models, and seamlessly interfaces with GIZMO, AREPO, GASOLINE, CHANGA, and ENZO. We demonstrate the capabilities of the code via three applications: a model for the star formation rate (SFR) - infrared luminosity relation in galaxies (including the impact of AGN); the impact of circumstellar dust around AGB stars on the mid-infrared emission from galaxy SEDs; and the impact of galaxy inclination angle on dust attenuation laws. △ Less

Submitted 18 June, 2020; originally announced June 2020.

Comments: Submitted to ApJ; comments welcome

On the Dopability of Semiconductors and Governing Material Properties

Authors: Anuj Goyal, Prashun Gorai, Shashwat Anand, Eric S. Toberer, G. Jeffrey Snyder, Vladan Stevanovic

Abstract: To be practical, semiconductors need to be doped. Sometimes, to nearly degenerate levels, e.g. in applications such as thermoelectric, transparent electronics or power electronics. However, many materials with finite band gaps are not dopable at all, while many others exhibit strong preference toward allowing either p- or n-type doping, but not both. In this work, we develop a model description of… ▽ More To be practical, semiconductors need to be doped. Sometimes, to nearly degenerate levels, e.g. in applications such as thermoelectric, transparent electronics or power electronics. However, many materials with finite band gaps are not dopable at all, while many others exhibit strong preference toward allowing either p- or n-type doping, but not both. In this work, we develop a model description of semiconductor dopability and formulate design principles in terms of governing materials properties. Our approach, which builds upon the semiconductor defect theory applied to a suitably devised (tight-binding) model system, reveals analytic relationships between intrinsic materials properties and the semiconductor dopability, and elucidates the role and the insufficiency of previously suggested descriptors such as the absolute band edge positions. We validate our model against a number of classic binary semiconductors and discuss its extension to more complex chemistries and the utility in large-scale material searches. △ Less

Submitted 18 May, 2020; originally announced May 2020.

Journal ref: Chem. Mater. 2020

Figuring Out Gas & Galaxies In Enzo (FOGGIE). IV. The Stochasticity of Ram Pressure Stripping in Galactic Halos

Authors: Raymond C. Simons, Molly S. Peeples, Jason Tumlinson, Brian W. O'Shea, Britton D. Smith, Lauren Corlies, Cassandra Lochhaas, Yong Zheng, Ramona Augustin, Deovrat Prasad, Gregory F. Snyder, Erik Tollerud

Abstract: We study ram pressure stripping in simulated Milky Way-like halos at z>=2 from the Figuring Out Gas & Galaxies In Enzo (FOGGIE) project. These simulations reach exquisite resolution in their circumgalactic medium (CGM) gas owing to FOGGIE's novel refinement scheme. The CGM of each halo spans a wide dynamic range in density and velocity over its volume---roughly 6 dex and 1000 km/s, respectively---… ▽ More We study ram pressure stripping in simulated Milky Way-like halos at z>=2 from the Figuring Out Gas & Galaxies In Enzo (FOGGIE) project. These simulations reach exquisite resolution in their circumgalactic medium (CGM) gas owing to FOGGIE's novel refinement scheme. The CGM of each halo spans a wide dynamic range in density and velocity over its volume---roughly 6 dex and 1000 km/s, respectively---translating into a 5 dex range in ram pressure imparted to interacting satellites. The ram pressure profiles of the simulated CGM are highly stochastic, owing to kpc-scale variations of the density and velocity fields of the CGM gas. As a result, the efficacy of ram pressure stripping depends strongly on the specific path a satellite takes through the CGM. The ram-pressure history of a single satellite is generally unpredictable and not well correlated with its approach vector with respect to the host galaxy. The cumulative impact of ram pressure on the simulated satellites is dominated by only a few short strong impulses---on average, 90% of the total surface momentum gained through ram pressure is imparted in 20% or less of the total orbital time. These results reveal an erratic mode of ram pressure stripping in Milky-Way like halos at high redshift---one that is not captured by a smooth spherically-averaged model of the circumgalactic medium. △ Less

Submitted 29 April, 2020; originally announced April 2020.

Comments: 18 pages, 10 figures. Submitted to ApJ

DeepMerge: Classifying High-redshift Merging Galaxies with Deep Neural Networks

Authors: A. Ćiprijanović, G. F. Snyder, B. Nord, J. E. G. Peek

Abstract: We investigate and demonstrate the use of convolutional neural networks (CNNs) for the task of distinguishing between merging and non-merging galaxies in simulated images, and for the first time at high redshifts (i.e. $z=2$). We extract images of merging and non-merging galaxies from the Illustris-1 cosmological simulation and apply observational and experimental noise that mimics that from the H… ▽ More We investigate and demonstrate the use of convolutional neural networks (CNNs) for the task of distinguishing between merging and non-merging galaxies in simulated images, and for the first time at high redshifts (i.e. $z=2$). We extract images of merging and non-merging galaxies from the Illustris-1 cosmological simulation and apply observational and experimental noise that mimics that from the Hubble Space Telescope; the data without noise form a "pristine" data set and that with noise form a "noisy" data set. The test set classification accuracy of the CNN is $79\%$ for pristine and $76\%$ for noisy. The CNN outperforms a Random Forest classifier, which was shown to be superior to conventional one- or two-dimensional statistical methods (Concentration, Asymmetry, the Gini, $M_{20}$ statistics etc.), which are commonly used when classifying merging galaxies. We also investigate the selection effects of the classifier with respect to merger state and star formation rate, finding no bias. Finally, we extract Grad-CAMs (Gradient-weighted Class Activation Mapping) from the results to further assess and interrogate the fidelity of the classification model. △ Less

Submitted 24 April, 2020; originally announced April 2020.

Comments: 17 pages, 8 figures, submitted to Astronomy & Computing

Analytical Models of Phonon-Point Defect Scattering

Authors: Ramya Gurunathan, Riley Hanus, Maxwell Dylla, Ankita Katre, G. Jeffrey Snyder

Abstract: Point defects exist widely in engineering materials and are known to scatter vibrational modes to reduce thermal conductivity. The Klemens description of point defect scattering is the most prolific analytical model for this effect. This work reviews the essential physics of the model and compares its predictions to first principles results for isotope and alloy scattering, demonstrating the model… ▽ More Point defects exist widely in engineering materials and are known to scatter vibrational modes to reduce thermal conductivity. The Klemens description of point defect scattering is the most prolific analytical model for this effect. This work reviews the essential physics of the model and compares its predictions to first principles results for isotope and alloy scattering, demonstrating the model as a useful materials design metric. A treatment of the scattering parameter ($Γ$) for a multiatomic lattice is recommended and compared to other treatments presented in literature, which have been at times misused to yield incomplete conclusions about the system's scattering mechanisms. Additionally, we demonstrate a reduced sensitivity of the model to the full phonon dispersion and discuss its origin. Finally, a simplified treatment of scattering in alloy systems with vacancies and interstitial defects is demonstrated to suitably describe the potent scattering strength of these off-stoichiometric defects. △ Less

Submitted 8 October, 2019; originally announced October 2019.

Comments: 10 pages, 7 pages of supplemental, 5 figures total

Journal ref: Phys. Rev. Applied 13, 034011 (2020)

On the need for synthetic data and robust data simulators in the 2020s

Authors: Molly S. Peeples, Bjorn Emonts, Mark Kyprianou, Matthew T. Penny, Gregory F. Snyder, Christopher C. Stark, Michael Troxel, Neil T. Zimmerman, John ZuHone

Abstract: As observational datasets become larger and more complex, so too are the questions being asked of these data. Data simulations, i.e., synthetic data with properties (pixelization, noise, PSF, artifacts, etc.) akin to real data, are therefore increasingly required for several purposes, including: (1) testing complicated measurement methods, (2) comparing models and astrophysical simulations to obse… ▽ More As observational datasets become larger and more complex, so too are the questions being asked of these data. Data simulations, i.e., synthetic data with properties (pixelization, noise, PSF, artifacts, etc.) akin to real data, are therefore increasingly required for several purposes, including: (1) testing complicated measurement methods, (2) comparing models and astrophysical simulations to observations in a manner that requires as few assumptions about the data as possible, (3) predicting observational results based on models and astrophysical simulations for, e.g., proposal planning, and (4) mitigating risk for future observatories and missions by effectively priming and testing pipelines. We advocate for an increase in using synthetic data to plan for and interpret real observations as a matter of routine. This will require funding for (1) facilities to provide robust data simulators for their instruments, telescopes, and surveys, and (2) making synthetic data publicly available in archives (much like real data) so as to lower the barrier of entry to all. △ Less

Submitted 16 July, 2019; originally announced July 2019.

Comments: An Astro2020 APC white paper

Journal ref: Astro2020: Decadal Survey on Astronomy and Astrophysics, APC white papers, no. 210; Bulletin of the American Astronomical Society, Vol. 51, Issue 7, id. 210 (2019)

Critical Behavior and Anisotropy in Single Crystal SrRuO$_3$

Authors: G. Jeffrey Snyder

Abstract: The magnetization of single crystal SrRuO3 is studied as a function of temperature along different crystallographic directions. The magnetocrystalline anisotropy and behavior near the critical transition temperature are analyzed in detail. The magnetization vs temperature is found to vary more like $T^2$ rather than $T^{3/2}$ expected for spin waves. The magnetization of single crystal SrRuO3 is studied as a function of temperature along different crystallographic directions. The magnetocrystalline anisotropy and behavior near the critical transition temperature are analyzed in detail. The magnetization vs temperature is found to vary more like $T^2$ rather than $T^{3/2}$ expected for spin waves. △ Less

Submitted 27 April, 2019; originally announced April 2019.

Comments: This work was first published in a thesis in 1997, reference [1]

Observational constraints on the merger history of galaxies since $z\approx6$: Probabilistic galaxy pair counts in the CANDELS fields

Authors: Kenneth Duncan, Christopher J. Conselice, Carl Mundy, Eric Bell, Jennifer Donley, Audrey Galametz, Yicheng Guo, Norman A. Grogin, Nimish P. Hathi, Jeyhan Kartaltepe, Dale Kocevski, Anton M. Koekemoer, Pablo G. Pérez-González, Kameswara B. Mantha, Gregory F. Snyder, Mauro Stefanon

Abstract: Galaxy mergers are expected to have a significant role in the mass assembly of galaxies in the early Universe, but there are very few observational constraints on the merger history of galaxies at $z>2$. We present the first study of galaxy major mergers (mass ratios $>$ 1:4) in mass-selected samples out to $z\approx6$. Using all five fields of the HST/CANDELS survey and a probabilistic pair count… ▽ More Galaxy mergers are expected to have a significant role in the mass assembly of galaxies in the early Universe, but there are very few observational constraints on the merger history of galaxies at $z>2$. We present the first study of galaxy major mergers (mass ratios $>$ 1:4) in mass-selected samples out to $z\approx6$. Using all five fields of the HST/CANDELS survey and a probabilistic pair count methodology that incorporates the full photometric redshift posteriors and corrections for stellar mass completeness, we measure galaxy pair-counts for projected separations between 5 and 30 kpc in stellar mass selected samples at $9.7 < \log_{10}(\rm{M}_{*}/\rm{M}_{\odot}) < 10.3$ and $\log_{10}(\rm{M}_{*}/\rm{M}_{\odot}) > 10.3$. We find that the major merger pair fraction rises with redshift to $z\approx6$ proportional to $(1+z)^{m}$, with $m = 0.8\pm0.2$ ($m = 1.8\pm0.2$) for $\log_{10}(\rm{M}_{*} / \rm{M}_{\odot}) > 10.3$ ($9.7 < \log_{10}(\rm{M}_{*}/\rm{M}_{\odot}) < 10.3$). Investigating the pair fraction as a function of mass ratio between 1:20 and 1:1, we find no evidence for a strong evolution in the relative numbers of minor to major mergers out to $z<3$. Using evolving merger timescales we find that the merger rate per galaxy ($\mathcal{R}$) rises rapidly from $0.07\pm 0.01$ Gyr$^{-1}$ at $z < 1$ to $7.6\pm 2.7$ Gyr$^{-1}$ at $z = 6$ for galaxies at $\log_{10}(\rm{M}_{*}/\rm{M}_{\odot}) > 10.3$. The corresponding co-moving major merger rate density remains roughly constant during this time, with rates of $Γ\approx 10^{-4}$ Gyr$^{-1}$ Mpc$^{-3}$. Based on the observed merger rates per galaxy, we infer specific mass accretion rates from major mergers that are comparable to the specific star-formation rates for the same mass galaxies at $z>3$ - observational evidence that mergers are as important a mechanism for building up mass at high redshift as in-situ star-formation. △ Less

Submitted 28 March, 2019; originally announced March 2019.

Comments: 36 pages, 17 figures. Accepted for publication in ApJ

Studying the Physical Properties of Tidal Features I. Extracting Morphological Substructure in CANDELS Observations and VELA Simulations

Authors: Kameswara Bharadwaj Mantha, Daniel H. McIntosh, Cody P. Ciaschi, Rubyet Evan, Henry C. Ferguson, Logan B. Fries, Yicheng Guo, Anton M. Koekemoer, Luther D. Landry, Elizabeth J. McGrath, Raymond C. Simons, Gregory F. Snyder, Scott E. Thompson, Eric F. Bell, Daniel Ceverino, Nimish P. Hathi, Camilla Pacifici, Joel R. Primack, Marc Rafelski, Vicente Rodriguez-Gomez

Abstract: The role of major mergers in galaxy evolution remains a key open question. Existing empirical merger identification methods use non-parametric and subjective visual classifications which can pose systematic challenges to constraining merger histories. As a first step towards overcoming these challenges, we develop and share publicly a new Python-based software tool that identifies and extracts the… ▽ More The role of major mergers in galaxy evolution remains a key open question. Existing empirical merger identification methods use non-parametric and subjective visual classifications which can pose systematic challenges to constraining merger histories. As a first step towards overcoming these challenges, we develop and share publicly a new Python-based software tool that identifies and extracts the flux-wise and area-wise significant contiguous regions from the model-subtracted "residual" images produced by popular parametric light-profile fitting tools (e.g., GALFIT). Using Hubble Space Telescope ($HST$) $H$-band single-Sérsic residual images of $17$ CANDELS galaxies, we demonstrate the tool's ability to measure the surface brightness and improve the qualitative identification of a variety of common residual features (disk structures, spiral substructures, plausible tidal features, and strong gravitational arcs). We test our method on synthetic $HST$ observations of a $z\sim 1.5$ major merger from the VELA hydrodynamic simulations. We extract $H$-band residual features corresponding to the birth, growth, and fading of tidal features during different stages and viewing orientations at CANDELS depths and resolution. We find that the extracted features at shallow depths have noisy visual appearance and are susceptible to viewing angle effects. For a VELA $z\sim 3$ major merger, we find that James Webb Space Telescope NIRCam observations can probe high-redshift tidal features with considerable advantage over existing $HST$ capabilities. Further quantitative analysis of plausible tidal features extracted with our new software hold promise for the robust identification of hallmark merger signatures and corresponding improvements to merger rate constraints. △ Less

Submitted 26 March, 2019; originally announced March 2019.

Comments: 18 pages, 11 Figures, and accepted for publication in MNRAS. The data products discussed in this paper along with the software are publicly available at https://github.com/AgentM-GEG/residual_feature_extraction