Showing 1–50 of 50 results for author: Jonsson, P

Extended atomic data for oxygen abundance analyses

Authors: W. Li, P. Jönsson, A. M. Amarsi, M. C. Li, J. Grumer

Abstract: As the most abundant element in the universe after hydrogen and helium, oxygen plays a key role in planetary, stellar, and galactic astrophysics. Its abundance is especially influential on stellar structure and evolution, and as the dominant opacity contributor at the base of the Sun's convection zone it is central to the discussion around the solar modelling problem. However, abundance analyses r… ▽ More As the most abundant element in the universe after hydrogen and helium, oxygen plays a key role in planetary, stellar, and galactic astrophysics. Its abundance is especially influential on stellar structure and evolution, and as the dominant opacity contributor at the base of the Sun's convection zone it is central to the discussion around the solar modelling problem. However, abundance analyses require complete and reliable sets of atomic data. We present extensive atomic data for O I, by using the multiconfiguration Dirac-Hartree-Fock and relativistic configuration interaction methods. Lifetimes and transition probabilities for radiative electric dipole transitions are given and compared with results from previous calculations and available measurements. The accuracy of the computed transition rates is evaluated by the differences between the transition rates in Babushkin and Coulomb gauges, as well as by a cancellation factor analysis. Out of the 989 computed transitions in this work, 205 are assigned to the accuracy classes AA-B, that is, with uncertainties less than 10%, following the criteria defined by the National Institute of Standards and Technology Atomic Spectra Database. We discuss the influence of the new log(gf) values on the solar oxygen abundance and ultimately advocate $\logε_{\mathrm{O}}=8.70\pm0.04$. △ Less

Submitted 5 April, 2023; originally announced April 2023.

Comments: 13 pages, 5 figures; Accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 674, A54 (2023)

Extended MCDHF calculations of energy levels and transition data for N I

Authors: Meichun Li, Wenxian Li, P. Jönsson, Anish M. Amarsi, Jon Grumer

Abstract: Accurate and extensive atomic data are essential for spectroscopic analyses of stellar atmospheres and other astronomical objects. We present energy levels, lifetimes, and transition probabilities for neutral nitrogen, the sixth most abundant element in the cosmos. The calculations employ the fully relativistic multiconfiguration Dirac-Hartree-Fock and relativistic configuration interaction method… ▽ More Accurate and extensive atomic data are essential for spectroscopic analyses of stellar atmospheres and other astronomical objects. We present energy levels, lifetimes, and transition probabilities for neutral nitrogen, the sixth most abundant element in the cosmos. The calculations employ the fully relativistic multiconfiguration Dirac-Hartree-Fock and relativistic configuration interaction methods, and span the 103 lowest states up to and including 2s$^2$2p$^2$5s. Our theoretical energies are in excellent agreement with the experimental data, with an average relative difference of 0.07%. In addition, our transition probabilities are in good agreement with available experimental and theoretical data. We further verify the agreement of our data with experimental results via a re-analysis of the solar nitrogen abundance, with the results from the Babushkin and Coulomb gauges consistent to 2% or 0.01 dex. We estimated the uncertainties of the computed transition data based on a statistical analysis of the differences between the transition rates in Babushkin and Coulomb gauges. Out of the 1701 computed electric dipole transitions in this work, 83 (536) are associated with uncertainties less than 5% (10%). △ Less

Submitted 31 January, 2023; originally announced February 2023.

Comments: 17 pages, 7 figures; Accepted for publication in The Astrophysical Journal Supplement Series

Extended theoretical transition data in C I - IV

Authors: W. Li, A. M. Amarsi, A. Papoulia, J. Ekman, P. Jönsson

Abstract: Accurate atomic data are essential for opacity calculations and for abundance analyses of the Sun and other stars. The aim of this work is to provide accurate and extensive results of energy levels and transition data for C I - IV. The Multiconfiguration Dirac-Hartree-Fock and relativistic configuration interaction methods were used in the present work. To improve the quality of the wave functions… ▽ More Accurate atomic data are essential for opacity calculations and for abundance analyses of the Sun and other stars. The aim of this work is to provide accurate and extensive results of energy levels and transition data for C I - IV. The Multiconfiguration Dirac-Hartree-Fock and relativistic configuration interaction methods were used in the present work. To improve the quality of the wave functions and reduce the relative differences between length and velocity forms for transition data involving high Rydberg states, alternative computational strategies were employed by imposing restrictions on the electron substitutions when constructing the orbital basis for each atom and ion. Transition data, e.g., weighted oscillator strengths and transition probabilities, are given for radiative electric dipole (E1) transitions involving levels up to 1s$^2$2s$^2$2p6s for C I, up to 1s$^2$2s$^2$7f for C II, up to 1s$^2$2s7f for C III, and up to 1s$^2$8g for CIV. Using the difference between the transition rates in length and velocity gauges as an internal validation, the average uncertainties of all presented E1 transitions are estimated to be 8.05%, 7.20%, 1.77%, and 0.28%, respectively, for C I - IV. Extensive comparisons with available experimental and theoretical results are performed and good agreement is observed for most of the transitions. In addition, the C I data were employed in a reanalysis of the solar carbon abundance. The new transition data give a line-by-line dispersion similar to the one obtained when using transition data that are typically used in stellar spectroscopic applications today. △ Less

Submitted 22 January, 2021; originally announced January 2021.

Comments: 21 pages, 4 figures

Supernova Model Discrimination with Hyper-Kamiokande

Authors: Hyper-Kamiokande Collaboration, :, K. Abe, P. Adrich, H. Aihara, R. Akutsu, I. Alekseev, A. Ali, F. Ameli, I. Anghel, L. H. V. Anthony, M. Antonova, A. Araya, Y. Asaoka, Y. Ashida, V. Aushev, F. Ballester, I. Bandac, M. Barbi, G. J. Barker, G. Barr, M. Batkiewicz-Kwasniak, M. Bellato, V. Berardi, M. Bergevin , et al. (478 additional authors not shown)

Abstract: Core-collapse supernovae are among the most magnificent events in the observable universe. They produce many of the chemical elements necessary for life to exist and their remnants -- neutron stars and black holes -- are interesting astrophysical objects in their own right. However, despite millennia of observations and almost a century of astrophysical study, the explosion mechanism of core-colla… ▽ More Core-collapse supernovae are among the most magnificent events in the observable universe. They produce many of the chemical elements necessary for life to exist and their remnants -- neutron stars and black holes -- are interesting astrophysical objects in their own right. However, despite millennia of observations and almost a century of astrophysical study, the explosion mechanism of core-collapse supernovae is not yet well understood. Hyper-Kamiokande is a next-generation neutrino detector that will be able to observe the neutrino flux from the next galactic core-collapse supernova in unprecedented detail. We focus on the first 500 ms of the neutrino burst, corresponding to the accretion phase, and use a newly-developed, high-precision supernova event generator to simulate Hyper-Kamiokande's response to five different supernova models. We show that Hyper-Kamiokande will be able to distinguish between these models with high accuracy for a supernova at a distance of up to 100 kpc. Once the next galactic supernova happens, this ability will be a powerful tool for guiding simulations towards a precise reproduction of the explosion mechanism observed in nature. △ Less

Submitted 20 July, 2021; v1 submitted 13 January, 2021; originally announced January 2021.

Comments: 21 pages, 7 figures. Article based on thesis published as arXiv:2002.01649. v2: added references and some explanations in response to reviewer comments

Journal ref: Astrophys.J. 916 (2021) 15

Large-scale multiconfiguration Dirac-Hartree-Fock calculations for astrophysics: n=4 levels in P-like ions from Mn~XI to Ni~XIV

Authors: C. X. Song, K. Wang, G. Del Zanna, P. Jönsson, R. Si, M. Godefroid, G. Gaigalas, L. Radži{ū}t{ė}, P. Rynkun, X. H. Zhao, J. Yan, C. Y. Chen

Abstract: Using the multiconfiguration Dirac-Hartree-Fock and the relativistic configuration interaction methods, a consistent set of transition energies and radiative transition data for the lowest 546 (623, 701, 745) states of the $3p^4 3d$, $3s 3p^2 3d^2$, $3s 3p^3 4p$, $3s 3p^4$, $3s^2 3d^3$, $3s^2 3p^2 3d$, $3s^2 3p^2 4d$, $3s^2 3p^2 4s$, $3p^3 3d^2$, $3p^5$, $3s 3p 3d^3$, $3s 3p^3 3d$, $3s 3p^3 4s$,… ▽ More Using the multiconfiguration Dirac-Hartree-Fock and the relativistic configuration interaction methods, a consistent set of transition energies and radiative transition data for the lowest 546 (623, 701, 745) states of the $3p^4 3d$, $3s 3p^2 3d^2$, $3s 3p^3 4p$, $3s 3p^4$, $3s^2 3d^3$, $3s^2 3p^2 3d$, $3s^2 3p^2 4d$, $3s^2 3p^2 4s$, $3p^3 3d^2$, $3p^5$, $3s 3p 3d^3$, $3s 3p^3 3d$, $3s 3p^3 4s$, $3s^2 3p 3d^2$, %$3s^2 3p^2 4f$, $3s^2 3p^2 4p$, $3s^2 3p^3$ configurations in Mn~XI (Fe~XII, Co~XIII, Ni~XIV) is provided. The comparison between calculated excitation energies for the $n=4$ states and available experimental values for Fe XII indicate that the calculations are highly accurate, with uncertainties of only a few hundred cm$^{-1}$. Lines from these states are prominent in the soft X-rays. With the present calculations, several recent new identifications are confirmed. Other identifications involving $3p^2 4d$ levels in Fe~XII that were found questionable are discussed and a few new assignments are recommended. As some $n=4$ states of the other ions also show large discrepancies between experimental and calculated energies, we reassess their identification. The present study provides highly accurate atomic data for the $n=4$ states of P-like ions of astrophysical interest, for which experimental data are scarce. △ Less

Submitted 28 February, 2020; originally announced March 2020.

Energy level structure and transition data of Er$^{2+}$

Authors: Gediminas Gaigalas, Pavel Rynkun, Laima Radžiūtė, Daiji Kato, Masaomi Tanaka, Per Jönsson

Abstract: The main aim of this paper is to present accurate energy levels of the ground [Xe]$4f^{12}$ and first excited [Xe]$4f^{11}5d$ configurations of Er$^{2+}$. The energy level structure of the Er$^{2+}$ ion was computed using the multiconfiguration Dirac-Hartree-Fock and relativistic configuration interaction (RCI) methods, as implemented in the GRASP2018 program package. The Breit interaction, self-e… ▽ More The main aim of this paper is to present accurate energy levels of the ground [Xe]$4f^{12}$ and first excited [Xe]$4f^{11}5d$ configurations of Er$^{2+}$. The energy level structure of the Er$^{2+}$ ion was computed using the multiconfiguration Dirac-Hartree-Fock and relativistic configuration interaction (RCI) methods, as implemented in the GRASP2018 program package. The Breit interaction, self-energy and vacuum polarization corrections were included in the RCI computations. The zero-first-order approach was used in the computations. Energy levels with the identification in $LS$ coupling for all (399) states belonging to the [Xe]$4f^{12}$ and [Xe]$4f^{11}5d$ configurations are presented. Electric dipole (E1) transition data between the levels of these two configurations are computed. The accuracy of the these data are evaluated by studying the behaviour of the transition rates as functions of the gauge parameter as well as by evaluating the cancellation factors. The core electron correlations were studied using different strategies. Root-mean-square deviations obtained in this study for states of the ground and excited configurations from the available experimental or semi-empirical data are 649 cm$^{-1}$, and 747 cm$^{-1}$, respectively. △ Less

Submitted 18 February, 2020; originally announced February 2020.

Comments: 20 pages, 9 tables, 17 figures

Large-scale multiconfiguration Dirac-Hartree-Fock calculations for astrophysics: Cl-like ions from Cr~VIII to Zn~XIV

Authors: K. Wang, P. Jönsson, G. Del Zanna, M. Godefroid, Z. B. Chen, C. Y. Chen, J. Yan

Abstract: We use the multiconfiguration Dirac-Hartree-Fock (MCDHF) method combined with the relativistic configuration interaction (RCI) approach (GRASP2K) to provide a consistent set of transition energies and radiative transition data for the lower $n =3$ states in all Cl-like ions of astrophysical importance, from \ion{Cr}{8} to \ion{Zn}{14}. We also provide excitation energies calculated for \mbox{Fe X}… ▽ More We use the multiconfiguration Dirac-Hartree-Fock (MCDHF) method combined with the relativistic configuration interaction (RCI) approach (GRASP2K) to provide a consistent set of transition energies and radiative transition data for the lower $n =3$ states in all Cl-like ions of astrophysical importance, from \ion{Cr}{8} to \ion{Zn}{14}. We also provide excitation energies calculated for \mbox{Fe X} using the many-body perturbation theory (MBPT, implemented within FAC). The comparison of the present MCDHF results with MBPT and with the available experimental energies indicates that the theoretical excitation energies are highly accurate, with uncertainties of only a few hundred cm$^{-1}$. Detailed comparisons for Fe~X and Ni~XII highlight discrepancies in the experimental energies found in the literature. Several new identifications are proposed. △ Less

Submitted 9 November, 2019; originally announced November 2019.

Hyper-Kamiokande Design Report

Authors: Hyper-Kamiokande Proto-Collaboration, :, K. Abe, Ke. Abe, H. Aihara, A. Aimi, R. Akutsu, C. Andreopoulos, I. Anghel, L. H. V. Anthony, M. Antonova, Y. Ashida, V. Aushev, M. Barbi, G. J. Barker, G. Barr, P. Beltrame, V. Berardi, M. Bergevin, S. Berkman, L. Berns, T. Berry, S. Bhadra, D. Bravo-Berguño, F. d. M. Blaszczyk , et al. (291 additional authors not shown)

Abstract: On the strength of a double Nobel prize winning experiment (Super)Kamiokande and an extremely successful long baseline neutrino programme, the third generation Water Cherenkov detector, Hyper-Kamiokande, is being developed by an international collaboration as a leading worldwide experiment based in Japan. The Hyper-Kamiokande detector will be hosted in the Tochibora mine, about 295 km away from th… ▽ More On the strength of a double Nobel prize winning experiment (Super)Kamiokande and an extremely successful long baseline neutrino programme, the third generation Water Cherenkov detector, Hyper-Kamiokande, is being developed by an international collaboration as a leading worldwide experiment based in Japan. The Hyper-Kamiokande detector will be hosted in the Tochibora mine, about 295 km away from the J-PARC proton accelerator research complex in Tokai, Japan. The currently existing accelerator will be steadily upgraded to reach a MW beam by the start of the experiment. A suite of near detectors will be vital to constrain the beam for neutrino oscillation measurements. A new cavern will be excavated at the Tochibora mine to host the detector. The experiment will be the largest underground water Cherenkov detector in the world and will be instrumented with new technology photosensors, faster and with higher quantum efficiency than the ones in Super-Kamiokande. The science that will be developed will be able to shape the future theoretical framework and generations of experiments. Hyper-Kamiokande will be able to measure with the highest precision the leptonic CP violation that could explain the baryon asymmetry in the Universe. The experiment also has a demonstrated excellent capability to search for proton decay, providing a significant improvement in discovery sensitivity over current searches for the proton lifetime. The atmospheric neutrinos will allow to determine the neutrino mass ordering and, together with the beam, able to precisely test the three-flavour neutrino oscillation paradigm and search for new phenomena. A strong astrophysical programme will be carried out at the experiment that will detect supernova neutrinos and will measure precisely solar neutrino oscillation. △ Less

Submitted 28 November, 2018; v1 submitted 9 May, 2018; originally announced May 2018.

Comments: 325 pages

Experimental and theoretical oscillator strengths of Mg I for accurate abundance analysis

Authors: A. Pehlivan Rhodin, H. Hartman, H. Nilsson, P. Jonsson

Abstract: Context. With the aid of stellar abundance analysis, it is possible to study the galactic formation and evolution. Magnesium is an important element to trace the alpha-element evolution in our Galaxy. For chemical abundance analysis, such as magnesium abundance, accurate and complete atomic data are essential. Inaccurate atomic data lead to uncertain abundances and prevent discrimination between d… ▽ More Context. With the aid of stellar abundance analysis, it is possible to study the galactic formation and evolution. Magnesium is an important element to trace the alpha-element evolution in our Galaxy. For chemical abundance analysis, such as magnesium abundance, accurate and complete atomic data are essential. Inaccurate atomic data lead to uncertain abundances and prevent discrimination between different evolution models. Aims. We study the spectrum of neutral magnesium from laboratory measurements and theoretical calculations. Our aim is to improve the oscillator strengths ( f -values) of Mg I lines and to create a complete set of accurate atomic data, particularly for the near-IR region. Methods. We derived oscillator strengths by combining the experimental branching fractions with radiative lifetimes reported in the literature and computed in this work. A hollow cathode discharge lamp was used to produce free atoms in the plasma and a Fourier transform spectrometer recorded the intensity-calibrated high-resolution spectra. In addition, we performed theoretical calculations using the multiconfiguration Hartree-Fock program ATSP2K. Results. This project provides a set of experimental and theoretical oscillator strengths. We derived 34 experimental oscillator strengths. Except from the Mg I optical triplet lines (3p 3^P^o_0,1,2 - 4s 3^S_1), these oscillator strengths are measured for the first time. The theoretical oscillator strengths are in very good agreement with the experimental data and complement the missing transitions of the experimental data up to n = 7 from even and odd parity terms. We present an evaluated set of oscillator strengths, gf, with uncertainties as small as 5%. The new values of the Mg I optical triplet line (3p 3^P^o_0,1,2 - 4s 3^S_1) oscillator strength values are ~0.08 dex larger than the previous measurements. △ Less

Submitted 21 November, 2016; originally announced November 2016.

Comments: Accepted by A&A

Journal ref: A&A 598, A102 (2017)

Abundances of disk and bulge giants from hi-res optical spectra: I. O, Mg, Ca, and Ti in the Solar neighborhood and Kepler field samples

Authors: H. Jönsson, N. Ryde, T. Nordlander, A. Pehlivan, H. Hartman, P. Jönsson, K. Eriksson

Abstract: The galactic bulge is a significant part of our galaxy, but it is hard to observe, being both distant and covered by dust in the disk. Therefore there do not exist many hi-res optical spectra of bulge stars with large wavelength coverage, whose determined abundances can be compared with nearby, similarly analyzed stellar samples. We aim to determine the, for chemical evolution models, so important… ▽ More The galactic bulge is a significant part of our galaxy, but it is hard to observe, being both distant and covered by dust in the disk. Therefore there do not exist many hi-res optical spectra of bulge stars with large wavelength coverage, whose determined abundances can be compared with nearby, similarly analyzed stellar samples. We aim to determine the, for chemical evolution models, so important alpha elements of a sample of bulge giants using hi-res optical spectra with large wavelength coverage. The abundances found will be compared to similarly derived abundances from similar spectra of similar stars in the local thin and thick disks. In this first paper we focus on the Solar neighborhood reference sample. We use spectral synthesis to derive the stellar parameters as well as the elemental abundances of both the local as well as the bulge samples of giants. Special care is taken to benchmark our method of determining stellar parameters against independent measurements of effective temperatures from angular diameter measurements and surface gravities from asteroseismology. In this first paper we present the method used to determine the stellar parameters as well as the elemental abundances, evaluate them, and present the results for our local disk sample of 291 giants. When comparing our determined spectroscopic temperatures to those derived from angular diameter measurements, we reproduce these with a systematic difference of +10 K and a standard deviation of 53 K. The spectroscopic gravities are reproducing the ones determined from asteroseismology with a systematic offset of +0.10 dex and a standard deviation of 0.12 dex. When it comes to the abundance trends, our sample of local disk giants is closely following that of other works analyzing solar neighborhood dwarfs, showing that the much brighter giant stars are as good abundance probes as the often used dwarfs. △ Less

Submitted 16 November, 2016; originally announced November 2016.

Comments: Accepted for publication in A&A

Journal ref: A&A 598, A100 (2017)

Plasma diagnostic potential of 2p4f in N$^+$ -- accurate wavelengths and oscillator strengths

Authors: Xiaozhi Shen, Jiguang Li, Per Jönsson, Jianguo Wang

Abstract: Radiative emission lines from nitrogen and its ions are often observed in nebulae spectra, where the N$^{2+}$ abundance can be inferred from lines of the 2p4f configuration. In addition, intensity ratios between lines of the 2p3p -- 2p3s and 2p4f -- 2p3d transition arrays can serve as temperature diagnostics. To aid abundance determinations and plasma diagnostics, wavelengths and oscillator streng… ▽ More Radiative emission lines from nitrogen and its ions are often observed in nebulae spectra, where the N$^{2+}$ abundance can be inferred from lines of the 2p4f configuration. In addition, intensity ratios between lines of the 2p3p -- 2p3s and 2p4f -- 2p3d transition arrays can serve as temperature diagnostics. To aid abundance determinations and plasma diagnostics, wavelengths and oscillator strengths were calculated with high-precision for electric-dipole (E1) transitions from levels in the 2p4f configuration of N$^{+}$. Electron correlation and relativistic effects, including the Breit interaction, were systematically taken into account within the framework of the multiconfiguration Dirac-Hartree-Fock (MCDHF) method. Except for the 2p4f - 2p4d transitions with quite large wavelengths and the two-electron-one-photon 2p4f -2s2p$^3$ transitions, the uncertainties of the present calculations were controlled to within 3% and 5% for wavelengths and oscillator strengths, respectively. We also compared our results with other theoretical and experimental values when available. Discrepancies were found between our calculations and previous calculations due to the neglect of relativistic effects in the latter. △ Less

Submitted 15 March, 2015; originally announced March 2015.

Hyperfine dependent $gf$-values of Mn I lines in the 1.49 $-$ 1.80 $μ$m H-band

Authors: Martin Andersson, Jon Grumer, Nils Ryde, Richard Blackwell-Whitehead, Roger Hutton, Yaming Zou, Per Jönsson, Tomas Brage

Abstract: The three Mn I lines at 17325, 17339 and 17349 Å~are among the 25 strongest lines ($\log(gf)>0.5$) in the H-band. They are all heavily broadened due to hyperfine structure and the profiles of these lines have so far not been understood. Earlier studies of these lines even suggested that they were blended. In this work, the profiles of these three infra-red (IR) lines have been studied theoreticall… ▽ More The three Mn I lines at 17325, 17339 and 17349 Å~are among the 25 strongest lines ($\log(gf)>0.5$) in the H-band. They are all heavily broadened due to hyperfine structure and the profiles of these lines have so far not been understood. Earlier studies of these lines even suggested that they were blended. In this work, the profiles of these three infra-red (IR) lines have been studied theoretically and compared to experimental spectra to assist in the complete understanding of the solar spectrum in the IR. It is shown that the structure of these lines can not be described in the conventional way by the diagonal $A$ and $B$ hyperfine interaction constants. The off-diagonal hyperfine interaction not only has large impact on the energies of the hyperfine levels, but also introduces a large intensity redistribution among the hyperfine lines, changing the line profiles dramatically. By performing large-scale calculations of the diagonal and off-diagonal hyperfine interaction and $gf$-values between the upper and lower hyperfine levels and using a semi-empirical fitting procedure, agreement between our synthetic and experimental spectra was achieved. Furthermore, we compare our results with observations of stellar spectra. The spectra of the Sun and the K1.5 III red giant star Arcturus were modelled in the relevant region, $1.73-1.74$ $μ$m using our theoretically predicted $gf$-values and energies for each individual hyperfine line. Satisfactory fits were obtained and clear improvements were found using our new data compared with the old available Mn I data. A complete list of energies and $gf$-values for all the $3d^54s({^7S})4d$ e$^{6}$D $-$ $3d^54s({^7S})4f$ w$^{6}$F hyperfine lines are available as supplementary material online, whereas only the stronger lines are presented and discussed in detail in this paper. △ Less

Submitted 11 June, 2014; originally announced June 2014.

Comments: 24 pages, 6 figures, 3 tables (more data tables will be available online when published in Astrophysical Journal Supplement Series)

The Spatial Extent and Distribution of Star Formation in 3D-HST Mergers at z~1.5

Authors: Kasper B. Schmidt, Hans-Walter Rix, Elisabete da Cunha, Gabriel B. Brammer, Thomas J. Cox, Pieter van Dokkum, Natascha M. Förster Schreiber, Marijn Franx, Mattia Fumagalli, Patrik Jonsson, Britt Lundgren, Michael V. Maseda, Ivelina Momcheva, Erica J. Nelson, Rosalind E. Skelton, Arjen van der Wel, Katherine E. Whitaker

Abstract: We present an analysis of the spatial distribution of star formation in a sample of 60 visually identified galaxy merger candidates at z>1. Our sample, drawn from the 3D-HST survey, is flux-limited and was selected to have high star formation rates based on fits of their broad-band, low spatial resolution spectral energy distributions. It includes plausible pre-merger (close pairs) and post-merger… ▽ More We present an analysis of the spatial distribution of star formation in a sample of 60 visually identified galaxy merger candidates at z>1. Our sample, drawn from the 3D-HST survey, is flux-limited and was selected to have high star formation rates based on fits of their broad-band, low spatial resolution spectral energy distributions. It includes plausible pre-merger (close pairs) and post-merger (single objects with tidal features) systems, with total stellar masses and star formation rates derived from multi-wavelength photometry. Here we use near-infrared slitless spectra from 3D-HST which produce Halpha or [OIII] emission line maps as proxies for star-formation maps. This provides a first comprehensive high-resolution, empirical picture of where star formation occurred in galaxy mergers at the epoch of peak cosmic star formation rate. We find that detectable star formation can occur in one or both galaxy centres, or in tidal tails. The most common case (58%) is that star formation is largely concentrated in a single, compact region, coincident with the centre of (one of) the merger components. No correlations between star formation morphology and redshift, total stellar mass, or star formation rate are found. A restricted set of hydrodynamical merger simulations between similarly massive and gas-rich objects implies that star formation should be detectable in both merger components, when the gas fractions of the individual components are the same. This suggests that z~1.5 mergers typically occur between galaxies whose gas fractions, masses, and/or star formation rates are distinctly different from one another. △ Less

Submitted 14 March, 2013; originally announced March 2013.

Comments: Accepted for publication in MNRAS, 16 pages, 10 figures

Global Star Formation Rates and Dust Emission Over the Galaxy Interaction Sequence

Authors: Lauranne Lanz, Andreas Zezas, Nicola Brassington, Howard A. Smith, Matthew L. N. Ashby, Elisabete da Cunha, Giovanni G. Fazio, Christopher C. Hayward, Lars Hernquist, Patrik Jonsson

Abstract: We measured and modeled the spectral energy distributions (SEDs) in 28 bands from the ultraviolet to the far-infrared (FIR) for 31 interacting galaxies in 14 systems. The sample is drawn from the Spitzer Interacting Galaxy Survey, which probes a range of galaxy interaction parameters at multiple wavelengths with an emphasis on the infrared bands. The subset presented in this paper consists of all… ▽ More We measured and modeled the spectral energy distributions (SEDs) in 28 bands from the ultraviolet to the far-infrared (FIR) for 31 interacting galaxies in 14 systems. The sample is drawn from the Spitzer Interacting Galaxy Survey, which probes a range of galaxy interaction parameters at multiple wavelengths with an emphasis on the infrared bands. The subset presented in this paper consists of all galaxies for which FIR Herschel SPIRE observations are publicly available. Our SEDs combine the Herschel photometry with multi-wavelength data from Spitzer, GALEX, Swift UVOT, and 2MASS. While the shapes of the SEDs are broadly similar across our sample, strongly interacting galaxies typically have more mid-infrared emission relative to their near-infrared and FIR emission than weakly or moderately interacting galaxies. We modeled the full SEDs to derive host galaxy star formation rates (SFR), specific star formation rates (sSFR), stellar masses, dust temperatures, dust luminosities, and dust masses. We find increases in the dust luminosity and mass, SFR, and cold (15-25 K) dust temperature as the interaction progresses from moderately to strongly interacting and between non- interacting and strongly interacting galaxies. We also find increases in the SFR between weakly and strongly interacting galaxies. In contrast, the sSFR remains unchanged across all the interaction stages. The SFR derived from the SED modeling agrees well with rates estimated by proportionality relations that depend on infrared emission. (abridged) △ Less

Submitted 20 February, 2013; originally announced February 2013.

Comments: 37 emulateapj pages, 11 tables, 25 figures, re-submitted to ApJ after revision per referee report

Modeling Mid-Infrared Diagnostics of Obscured Quasars and Starbursts

Authors: Gregory F. Snyder, Christopher C. Hayward, Anna Sajina, Patrik Jonsson, Thomas J. Cox, Lars Hernquist, Philip F. Hopkins, Lin Yan

Abstract: We analyze the link between active galactic nuclei (AGN) and mid-infrared flux using dust radiative transfer calculations of starbursts realized in hydrodynamical simulations. Focusing on the effects of galaxy dust, we evaluate diagnostics commonly used to disentangle AGN and star formation in ultraluminous infrared galaxies (ULIRGs). We examine these quantities as a function of time, viewing angl… ▽ More We analyze the link between active galactic nuclei (AGN) and mid-infrared flux using dust radiative transfer calculations of starbursts realized in hydrodynamical simulations. Focusing on the effects of galaxy dust, we evaluate diagnostics commonly used to disentangle AGN and star formation in ultraluminous infrared galaxies (ULIRGs). We examine these quantities as a function of time, viewing angle, dust model, AGN spectrum, and AGN strength in merger simulations representing two possible extremes of the ULIRG population: one is a typical gas-rich merger at z ~ 0, and the other is characteristic of extremely obscured starbursts at z ~ 2 to 4. This highly obscured burst begins star-formation-dominated with significant PAH emission, and ends with a ~10^9 yr period of red near-IR colors. At coalescence, when the AGN is most luminous, dust obscures the near-infrared AGN signature, reduces the relative emission from polycyclic aromatic hydrocarbons (PAHs), and enhances the 9.7 micron absorption by silicate grains. Although generally consistent with previous interpretations, our results imply none of these indicators can unambiguously estimate the AGN luminosity fraction in all cases. Motivated by the simulations, we show that a combination of the extinction feature at 9.7 micron, the PAH strength, and a near-infrared slope can simultaneously constrain the AGN fraction and dust grain distribution for a wide range of obscuration. We find that this indicator, accessible to the James Webb Space Telescope, may estimate the AGN power as tightly as the hard X-ray flux alone, thereby providing a valuable future cross-check and constraint for large samples of distant ULIRGs. △ Less

Submitted 11 April, 2013; v1 submitted 23 October, 2012; originally announced October 2012.

Comments: 21 pages, 13 figures, ApJ accepted version

Journal ref: Astrophys.J. 768 (2013) 168

Submillimetre galaxies in a hierarchical universe: number counts, redshift distribution, and implications for the IMF

Authors: Christopher C. Hayward, Desika Narayanan, Dušan Kereš, Patrik Jonsson, Philip F. Hopkins, T. J. Cox, Lars Hernquist

Abstract: High-redshift submillimetre galaxies (SMGs) are some of the most rapidly star-forming galaxies in the Universe. Historically, galaxy formation models have had difficulty explaining the observed number counts of SMGs. We combine a semi-empirical model with 3-D hydrodynamical simulations and 3-D dust radiative transfer to predict the number counts of unlensed SMGs. Because the stellar mass functions… ▽ More High-redshift submillimetre galaxies (SMGs) are some of the most rapidly star-forming galaxies in the Universe. Historically, galaxy formation models have had difficulty explaining the observed number counts of SMGs. We combine a semi-empirical model with 3-D hydrodynamical simulations and 3-D dust radiative transfer to predict the number counts of unlensed SMGs. Because the stellar mass functions, gas and dust masses, and sizes of our galaxies are constrained to match observations, we can isolate uncertainties related to the dynamical evolution of galaxy mergers and the dust radiative transfer. The number counts and redshift distributions predicted by our model agree well with observations. Isolated disc galaxies dominate the faint (S_1.1 < ~1 mJy, or S_850 < ~2 mJy) population. The brighter sources are a mix of merger-induced starbursts and galaxy-pair SMGs; the latter subpopulation accounts for ~30 - 50 per cent of all SMGs at all S_1.1 > ~0.5 mJy (S_850 > ~1 mJy). The mean redshifts are ~3.0 - 3.5, depending on the flux cut, and the brightest sources tend to be at higher redshifts. Because the galaxy-pair SMGs will be resolved into multiple fainter sources by ALMA, the bright ALMA counts should be as much as 2 times less than those observed using single-dish telescopes. The agreement between our model, which uses a Kroupa IMF, and observations suggests that the IMF in high-redshifts starbursts need not be top-heavy; if the IMF were top-heavy, our model would over-predict the number counts. We conclude that the difficulty some models have reproducing the observed SMG counts is likely indicative of more general problems -- such as an under-prediction of the abundance of massive galaxies or a star formation rate--stellar mass relation normalisation lower than that observed -- rather than a problem specific to the SMG population. △ Less

Submitted 15 October, 2012; v1 submitted 11 September, 2012; originally announced September 2012.

Comments: Accepted for publication in MNRAS. Insignificant changes from the previous version

Journal ref: MNRAS 428 (2013) 2529-2547

On Galaxies and Homology

Authors: Gregory S. Novak, Patrik Jonsson, Joel R. Primack, Thomas J. Cox, Avishai Dekel

Abstract: The definition of homology for single-component galaxies is clear, but for multi-component (luminous and dark matter) galaxies there is some ambiguity. We attempt to clarify the situation by carefully separating the different concepts of homology that have been used to date. We argue that the most useful definition is that a set of galaxies is homologous if they are the same in all respects up to… ▽ More The definition of homology for single-component galaxies is clear, but for multi-component (luminous and dark matter) galaxies there is some ambiguity. We attempt to clarify the situation by carefully separating the different concepts of homology that have been used to date. We argue that the most useful definition is that a set of galaxies is homologous if they are the same in all respects up to a set of three dimensional scaling constants which may differ from one galaxy to the next. Noting that we are free to choose the dimensional constants, we find that a set of hydrodynamic simulated galaxy merger remnants is significantly closer to homologous when the dimensional length constant is taken to be the radius containing equal amounts of dark and baryonic matter rather than the usual observationally motivated choice of the baryonic half-mass radius. Once the correct dimensional scaling constants are used, the stellar velocity dispersion anisotropy is essentially the sole source of the variation in the kinematic structure of these simulated merger remnants. In order to facilitate the use of these scaling constants to analyse observed galaxies, we calculated the relationship between our preferred dimensional scaling constants and the typical observationally accessible quantities. △ Less

Submitted 11 May, 2012; originally announced May 2012.

Comments: 15 pages, 9 figures, 3 tables. Accepted to MNRAS

How to distinguish starbursts and quiescently star-forming galaxies: The `bimodal' submillimetre galaxy population as a case study

Authors: Christopher C. Hayward, Patrik Jonsson, Dušan Kereš, Benjamin Magnelli, Lars Hernquist, T. J. Cox

Abstract: In recent work (arXiv:1101.0002) we have suggested that the high-redshift (z ~ 2-4) bright submillimetre galaxy (SMG) population is heterogeneous, with major mergers contributing both at early stages, where quiescently star-forming discs are blended into one submm source (`galaxy-pair SMGs'), and late stages, where mutual tidal torques drive gas inflows and cause strong starbursts. Here we combine… ▽ More In recent work (arXiv:1101.0002) we have suggested that the high-redshift (z ~ 2-4) bright submillimetre galaxy (SMG) population is heterogeneous, with major mergers contributing both at early stages, where quiescently star-forming discs are blended into one submm source (`galaxy-pair SMGs'), and late stages, where mutual tidal torques drive gas inflows and cause strong starbursts. Here we combine hydrodynamic simulations of major mergers with 3-D dust radiative transfer calculations to determine observational diagnostics that can distinguish between quiescently star-forming SMGs and starburst SMGs via integrated data alone. We fit the far-IR SEDs of the simulated galaxies with the optically thin single-temperature modified blackbody, the full form of the single-temperature modified blackbody, and a power-law temperature-distribution model. The effective dust temperature, T_dust, and power-law index of the dust emissivity in the far-IR, β, derived can significantly depend on the fitting form used, and the intrinsic β of the dust is not recovered. However, for all forms used here, there is a T_dust above which almost all simulated galaxies are starbursts, so a T_dust cut is very effective at selecting starbursts. Simulated merger-induced starbursts also have higher L_IR/M_gas and L_IR/L_FUV than quiescently star-forming galaxies and lie above the star formation rate-stellar mass relation. These diagnostics can be used to test our claim that the SMG population is heterogeneous and to observationally determine what star formation mode dominates a given galaxy population. We comment on applicability of these diagnostics to ULIRGs that would not be selected as SMGs. These `hot-dust ULIRGs' are typically starburst galaxies lower in mass than SMGs, but they can also simply be SMGs observed from a different viewing angle. △ Less

Submitted 9 May, 2012; v1 submitted 6 March, 2012; originally announced March 2012.

Comments: 21 pages, 11 figures. Accepted for publication in MNRAS. Minor changes to text but otherwise identical to v1

Journal ref: MNRAS 424 (2012) 951-970

The SEDs of Interacting Galaxies

Authors: Lauranne Lanz, Nicola Brassington, Andreas Zezas, Howard A. Smith, Matthew L. N. Ashby, Elisabete da Cunha, Christopher Klein, Patrik Jonsson, Christopher C. Hayward, Lars Hernquist, Giovanni Fazio

Abstract: The evolution of galaxies is greatly influenced by their interactions. As part of a program to study interacting galaxies, we have measured and modeled the spectral energy distri- butions (SEDs) from the ultraviolet (UV) to the far-infrared (FIR). We describe the constraints imposed on star formation histories by these SEDs, and the variations therein seen across the interaction sequence, and we c… ▽ More The evolution of galaxies is greatly influenced by their interactions. As part of a program to study interacting galaxies, we have measured and modeled the spectral energy distri- butions (SEDs) from the ultraviolet (UV) to the far-infrared (FIR). We describe the constraints imposed on star formation histories by these SEDs, and the variations therein seen across the interaction sequence, and we compare the results of different star formation rate prescriptions applied to the data. The sample itself is based on the Spitzer Interacting Galaxy Survey (SIGS) of 111 galaxies in 50 systems, a project designed to probe a range of galaxy interaction parameters in the infrared. Our SEDs combine the Spitzer results with multiwavelength data from other missions, in particular GALEX and Herschel. The subset presented here is the sample for which FIR Herschel observations are currently publicly available. △ Less

Submitted 18 November, 2011; originally announced November 2011.

Comments: 5 pages, 4 figures, to appear in the proceedings of "The Spectral Energy Distribution of Galaxies", Preston England, September 2011

The Major and Minor Galaxy Merger Rates at z < 1.5

Authors: Jennifer M. Lotz, Patrik Jonsson, T. J. Cox, Darren Croton, Joel R. Primack, Rachel S. Somerville, Kyle Stewart

Abstract: Calculating the galaxy merger rate requires both a census of galaxies identified as merger candidates, and a cosmologically-averaged `observability' timescale T_obs(z) for identifying galaxy mergers. While many have counted galaxy mergers using a variety of techniques, T_obs(z) for these techniques have been poorly constrained. We address this problem by calibrating three merger rate estimators wi… ▽ More Calculating the galaxy merger rate requires both a census of galaxies identified as merger candidates, and a cosmologically-averaged `observability' timescale T_obs(z) for identifying galaxy mergers. While many have counted galaxy mergers using a variety of techniques, T_obs(z) for these techniques have been poorly constrained. We address this problem by calibrating three merger rate estimators with a suite of hydrodynamic merger simulations and three galaxy formation models. We estimate T_obs(z) for (1) close galaxy pairs with a range of projected separations, (2) the morphology indicator G-M20, and (3) the morphology indicator asymmetry A. Then we apply these timescales to the observed merger fractions at z < 1.5 from the recent literature. When our physically-motivated timescales are adopted, the observed galaxy merger rates become largely consistent. The remaining differences between the galaxy merger rates are explained by the differences in the range of mass-ratio measured by different techniques and differing parent galaxy selection. The major merger rate per unit co-moving volume for samples selected with constant number density evolves much more strongly with redshift (~ (1+z)^(+3.0 \pm 1.1)) than samples selected with constant stellar mass or passively evolving luminosity (~ (1+z)^(+0.1 \pm 0.4)). We calculate the minor merger rate (1:4 < M_{sat}/M_{primary} <~ 1:10) by subtracting the major merger rate from close pairs from the `total' merger rate determined by G-M20. The implied minor merger rate is ~3 times the major merger rate at z ~ 0.7, and shows little evolution with redshift. △ Less

Submitted 11 August, 2011; originally announced August 2011.

Comments: accepted by ApJ; 25 pages; 13 figures; 4 tables

K+A Galaxies as the Aftermath of Gas-Rich Mergers: Simulating the Evolution of Galaxies as Seen by Spectroscopic Surveys

Authors: Gregory F. Snyder, Thomas J. Cox, Christopher C. Hayward, Lars Hernquist, Patrik Jonsson

Abstract: Models of poststarburst (or "K+A") galaxies are constructed by combining fully three-dimensional hydrodynamic simulations of galaxy mergers with radiative transfer calculations of dust attenuation. Spectral line catalogs are generated automatically from moderate-resolution optical spectra calculated as a function of merger progress in each of a large suite of simulations. The mass, gas fraction, o… ▽ More Models of poststarburst (or "K+A") galaxies are constructed by combining fully three-dimensional hydrodynamic simulations of galaxy mergers with radiative transfer calculations of dust attenuation. Spectral line catalogs are generated automatically from moderate-resolution optical spectra calculated as a function of merger progress in each of a large suite of simulations. The mass, gas fraction, orbital parameters, and mass ratio of the merging galaxies are varied systematically, showing that the lifetime and properties of the K+A phase are strong functions of merger scenario. K+A durations are generally less than ~0.1-0.3 Gyr, significantly shorter than the commonly assumed 1 Gyr, which is obtained only in rare cases, owing to a wide variation in star formation histories resulting from different orbital and progenitor configurations. Combined with empirical merger rates, the model lifetimes predict rapidly-rising K+A fractions as a function of redshift that are consistent with results of large spectroscopic surveys, resolving tension between the observed K+A abundance and that predicted when one assumes the K+A duration is the lifetime of A stars (~1 Gyr). The effects of dust attenuation, viewing angle, and aperture bias on our models are analyzed. In some cases, the K+A features are longer-lived and more pronounced when AGN feedback removes dust from the center, uncovering the young stars formed during the burst. In this picture, the K+A phase begins during or shortly after the bright starburst/AGN phase in violent mergers, and thus offers a unique opportunity to study the effects of quasar and star formation feedback on the gas reservoir and evolution of the remnant. Analytic fitting formulae are provided for the estimates of K+A incidence as a function of merger scenario. △ Less

Submitted 20 October, 2011; v1 submitted 17 February, 2011; originally announced February 2011.

Comments: 26 pages, 13 figures; ApJ; minor changes to reflect accepted version

Journal ref: Astrophys.J. 741 (2011) 77

What Does a Submillimeter Galaxy Selection Actually Select? The Dependence of Submillimeter Flux Density on Star Formation Rate and Dust Mass

Authors: Christopher C. Hayward, Dušan Kereš, Patrik Jonsson, Desika Narayanan, T. J. Cox, Lars Hernquist

Abstract: We perform 3-D dust radiative transfer (RT) calculations on hydrodynamic simulations of isolated and merging disk galaxies in order to quantitatively study the dependence of observed-frame submillimeter (submm) flux density on galaxy properties. We find that submm flux density and star formation rate (SFR) are related in dramatically different ways for quiescently star-forming galaxies and starbur… ▽ More We perform 3-D dust radiative transfer (RT) calculations on hydrodynamic simulations of isolated and merging disk galaxies in order to quantitatively study the dependence of observed-frame submillimeter (submm) flux density on galaxy properties. We find that submm flux density and star formation rate (SFR) are related in dramatically different ways for quiescently star-forming galaxies and starbursts. Because the stars formed in the merger-induced starburst do not dominate the bolometric luminosity and the rapid drop in dust mass and more compact geometry cause a sharp increase in dust temperature during the burst, starbursts are very inefficient at boosting submm flux density (e.g., a $\ga16$x boost in SFR yields a $\la 2$x boost in submm flux density). Moreover, the ratio of submm flux density to SFR differs significantly between the two modes; thus one cannot assume that the galaxies with highest submm flux density are necessarily those with the highest bolometric luminosity or SFR. These results have important consequences for the bright submillimeter-selected galaxy (SMG) population. Among them are: 1. The SMG population is heterogeneous. In addition to merger-driven starbursts, there is a subpopulation of galaxy pairs, where two disks undergoing a major merger but not yet strongly interacting are blended into one submm source because of the large ($\ga 15$", or $\sim 130$ kpc at $z = 2$) beam of single-dish submm telescopes. 2. SMGs must be very massive ($M_{\star} \ga 6 \times 10^{10} \msun$). 3. The infall phase makes the SMG duty cycle a factor of a few greater than what is expected for a merger-driven starburst. (Abridged.) △ Less

Submitted 12 September, 2011; v1 submitted 29 December, 2010; originally announced January 2011.

Comments: 15 pages, 3 figures. Accepted for publication in ApJ. Minor changes from v1

Journal ref: ApJ 743 (2011) 159

Interpreting the Evolution of the Size - Luminosity Relation for Disk Galaxies from Redshift 1 to the Present

Authors: Alyson Brooks, Adam Solomon, Fabio Governato, Jacqueline McCleary, Lauren MacArthur, Chris Brook, Patrik Jonsson, Tom Quinn, James Wadsley

Abstract: A sample of very high resolution cosmological disk galaxy simulations is used to investigate the evolution of galaxy disk sizes back to redshift 1 within the Lambda CDM cosmology. Artificial images in the rest frame B band are generated, allowing for a measurement of disk scale lengths using surface brightness profiles as observations would, and avoiding any assumption that light must follow mass… ▽ More A sample of very high resolution cosmological disk galaxy simulations is used to investigate the evolution of galaxy disk sizes back to redshift 1 within the Lambda CDM cosmology. Artificial images in the rest frame B band are generated, allowing for a measurement of disk scale lengths using surface brightness profiles as observations would, and avoiding any assumption that light must follow mass as previous models have assumed. We demonstrate that these simulated disks are an excellent match to the observed magnitude - size relation for both local disks, and for disks at z=1 in the magnitude/mass range of overlap. We disentangle the evolution seen in the population as a whole from the evolution of individual disk galaxies. In agreement with observations, our simulated disks undergo roughly 1.5 magnitudes/arcsec^2 of surface brightness dimming since z=1. We find evidence that evolution in the magnitude - size plane varies by mass, such that galaxies with M* > 10^9 M_sun undergo more evolution in size than luminosity, while dwarf galaxies tend to evolve potentially more in luminosity. The disks grow in such a way as to stay on roughly the same stellar mass - size relation with time. Finally, due to an evolving stellar mass - SFR relation, a galaxy at a given stellar mass (or size) at z=1 will reside in a more massive halo and have a higher SFR, and thus a higher luminosity, than a counterpart of the same stellar mass at z=0. △ Less

Submitted 1 November, 2010; originally announced November 2010.

Comments: Version resubmitted to ApJ, after referee's comments

Do sub-millimeter galaxy number counts provide evidence for a top-heavy IMF?

Authors: Christopher C. Hayward, Desika Narayanan, Patrik Jonsson, T. J. Cox, Dušan Kereš, Philip F. Hopkins, Lars Hernquist

Abstract: Theoretical models have had difficulty matching the observed number density of sub-millimeter galaxies (SMGs), causing some authors (e.g., Baugh et al. 2005) to suggest that SMGs provide evidence for a top-heavy initial mass function (IMF). To test this claim, we have, for the first time, combined high-resolution 3-D hydrodynamic simulations of isolated and merging massive, gas-rich galaxies, radi… ▽ More Theoretical models have had difficulty matching the observed number density of sub-millimeter galaxies (SMGs), causing some authors (e.g., Baugh et al. 2005) to suggest that SMGs provide evidence for a top-heavy initial mass function (IMF). To test this claim, we have, for the first time, combined high-resolution 3-D hydrodynamic simulations of isolated and merging massive, gas-rich galaxies, radiative transfer, and a semi-empirical merger rate model to predict the number density of SMGs. Our model can reproduce the observed SMG number density even when using a standard (Kroupa) IMF. The agreement is due to a combination of relatively long sub-mm duty cycles for mergers (a few times 10^8 years for our most massive models), which owe to our combination of high-resolution 3-D hydrodynamic simulations and dust radiative transfer; sufficient number densities of massive, gas-rich mergers; and the decrease in sub-mm counts observed by recent deep/wide surveys (e.g., Austermann et al. 2010) relative to previous surveys. Our results suggest that the observed SMG number counts do not provide evidence for a top-heavy IMF at high redshift. △ Less

Submitted 30 August, 2010; v1 submitted 26 August, 2010; originally announced August 2010.

Comments: Minor typo in the abstract corrected

On Sizes, Kinematics, M/L Gradients, and Light Profiles of Massive Compact Galaxies at z~2

Authors: Stijn Wuyts, Thomas J. Cox, Christopher C. Hayward, Marijn Franx, Lars Hernquist, Philip F. Hopkins, Patrik Jonsson, Pieter G. van Dokkum

Abstract: We present a detailed analysis of the structure and resolved stellar populations of simulated merger remnants, and compare them to observations of compact quiescent galaxies at z ~ 2. We find that major merging is a viable mechanism to produce systems of ~ 10^11 Msun and ~ 1 kpc size, provided the gas fraction at the time of final coalescence is high (~ 40%), and provided that the progenitors are… ▽ More We present a detailed analysis of the structure and resolved stellar populations of simulated merger remnants, and compare them to observations of compact quiescent galaxies at z ~ 2. We find that major merging is a viable mechanism to produce systems of ~ 10^11 Msun and ~ 1 kpc size, provided the gas fraction at the time of final coalescence is high (~ 40%), and provided that the progenitors are compact star-forming galaxies, as expected at high redshift. Their integrated spectral energy distributions and velocity dispersions are in good agreement with the observations, and their position in the (v_{maj}/sigma, ellipticity) diagram traces the upper envelope of the distribution of lower redshift early-type galaxies. The simulated merger remnants show time- and sightline-dependent M/L ratio gradients that result from a superposition of radially dependent stellar age, stellar metallicity, and extinction. The median ratio of effective radius in rest-frame V-band light to that in mass surface density is ~ 2 during the quiescent remnant phase. This is typically expressed by a negative color gradient (i.e., red core), which we expect to correlate with the integrated color of the system. Finally, the simulations differ from the observations in their surface brightness profile shape. The simulated remnants are typically best fit by high (n >> 4) Sersic indices, whereas observed quiescent galaxies at z ~ 2 tend to be less cuspy (median n ~ 2.3). Limiting early star formation in the progenitors may be required to prevent the simulated merger remnants from having extended wings. △ Less

Submitted 24 August, 2010; originally announced August 2010.

Comments: Accepted for publication in The Astrophysical Journal, 21 pages, 17 figures

An observer's view of simulated galaxies: disc-to-total ratios, bars, and (pseudo-)bulges

Authors: Cecilia Scannapieco, Dimitri A. Gadotti, Patrik Jonsson, Simon D. M. White

Abstract: We use cosmological hydrodynamical simulations of the formation of Milky Way-mass galaxies to study the relative importance of the main stellar components, i.e., discs, bulges, and bars, at redshift zero. The main aim of this work is to understand if estimates of the structural parameters of these components determined from kinematics (as is usually done in simulations) agree well with those obtai… ▽ More We use cosmological hydrodynamical simulations of the formation of Milky Way-mass galaxies to study the relative importance of the main stellar components, i.e., discs, bulges, and bars, at redshift zero. The main aim of this work is to understand if estimates of the structural parameters of these components determined from kinematics (as is usually done in simulations) agree well with those obtained using a photometric bulge/disc/bar decomposition (as done in observations). To perform such a comparison, we have produced synthetic observations of the simulation outputs with the Monte-Carlo radiative transfer code SUNRISE and used the BUDDA code to make 2D photometric decompositions of the resulting images (in the i and g bands). We find that the kinematic disc-to-total ratio (D/T) estimates are systematically and significantly lower than the photometric ones. While the maximum D/T ratios obtained with the former method are of the order of 0.2, they are typically >0.4, and can be as high as 0.7, according to the latter. The photometric decomposition shows that many of the simulated galaxies have bars, with Bar/T ratios in the range 0.2-0.4, and that bulges have in all cases low Sersic indices, resembling observed pseudo-bulges instead of classical ones. Simulated discs, bulges and bars generally have similar (g-i) colours, which are in the blue tail of the distribution of observed colours. This is not due to the presence of young stars, but rather to low metallicities and poor gas content in the simulated galaxies, which makes dust extinction low. Photometric decompositions thus match the component ratios usually quoted for spiral galaxies better than kinematic decompositions, but the shift is insufficient to make the simulations consistent with observed late-type systems. △ Less

Submitted 11 June, 2010; v1 submitted 27 January, 2010; originally announced January 2010.

Comments: MNRAS Letters accepted, revised version, conclusions unchanged

The Effect of Gas Fraction on the Morphology and Time-scales of Disc Galaxy Mergers

Authors: Jennifer M. Lotz, Patrik Jonsson, T. J. Cox, Joel R. Primack

Abstract: Gas-rich galaxy mergers are more easily identified by their disturbed morphologies than mergers with less gas. Because the typical gas fraction of galaxy mergers is expected to increase with redshift, the under-counting of low gas-fraction mergers may bias morphological estimates of the evolution of galaxy merger rate. To understand the magnitude of this bias, we explore the effect of gas fracti… ▽ More Gas-rich galaxy mergers are more easily identified by their disturbed morphologies than mergers with less gas. Because the typical gas fraction of galaxy mergers is expected to increase with redshift, the under-counting of low gas-fraction mergers may bias morphological estimates of the evolution of galaxy merger rate. To understand the magnitude of this bias, we explore the effect of gas fraction on the morphologies of a series of simulated disc galaxy mergers. With the resulting g-band images, we determine how the time-scale for identifying major and minor galaxy mergers via close projected pairs and quantitative morphology (the Gini coefficient G, the second-order moment of the brightest 20% of the light M20, and asymmetry A) depends on baryonic gas fraction f(gas). Strong asymmetries last significantly longer in high gas-fraction mergers of all mass ratios, with time-scales ranging from >= 300 Myr for f(gas) ~ 20% to >= 1 Gyr for f(gas) ~ 50%. Therefore the strong evolution with redshift observed in the fraction of asymmetric galaxies may reflect evolution in the gas properties of galaxies rather than the global galaxy merger rate. On the other hand, the time-scale for identifying a galaxy merger via G-M20 is weakly dependent on gas-fraction (~ 200-400 Myr), consistent with the weak evolution observed for G-M20 mergers. △ Less

Submitted 8 December, 2009; originally announced December 2009.

Comments: 15 pages; resubmitted to MNRAS

The Effect of Mass Ratio on the Morphology and Time-scales of Disc Galaxy Mergers

Authors: Jennifer M. Lotz, Patrik Jonsson, T. J. Cox, Joel R. Primack

Abstract: The majority of galaxy mergers are expected to be minor mergers. The observational signatures of minor mergers are not well understood, thus there exist few constraints on the minor merger rate. This paper seeks to address this gap in our understanding by determining if and when minor mergers exhibit disturbed morphologies and how they differ from the morphology of major mergers. We simulate a s… ▽ More The majority of galaxy mergers are expected to be minor mergers. The observational signatures of minor mergers are not well understood, thus there exist few constraints on the minor merger rate. This paper seeks to address this gap in our understanding by determining if and when minor mergers exhibit disturbed morphologies and how they differ from the morphology of major mergers. We simulate a series of unequal-mass moderate gas-fraction disc galaxy mergers. With the resulting g-band images, we determine how the time-scale for identifying galaxy mergers via projected separation and quantitative morphology (the Gini coefficient G, asymmetry A, and the second-order moment of the brightest 20% of the light M20) depends on the merger mass ratio, relative orientations and orbital parameters. We find that G-M20 is as sensitive to 9:1 baryonic mass ratio mergers as 1:1 mergers, with observability time-scales ~ 0.2-0.4 Gyr. In contrast, asymmetry finds mergers with baryonic mass ratios between 4:1 and 1:1 (assuming local disc galaxy gas-fractions). Asymmetry time-scales for moderate gas-fraction major disc mergers are ~ 0.2-0.4 Gyr, and less than 0.06 Gyr for moderate gas-fraction minor mergers. The relative orientations and orbits have little effect on the time-scales for morphological disturbances. Observational studies of close pairs often select major mergers by choosing paired galaxies with similar luminosities and/or stellar masses. Therefore, the various ways of finding galaxy mergers (G-M20, A, close pairs) are sensitive to galaxy mergers of different mass ratios. By comparing the frequency of mergers selected by different techniques, one may place empirical constraints on the major and minor galaxy merger rates. △ Less

Submitted 8 December, 2009; originally announced December 2009.

Comments: 16 pages; resubmitted to MNRAS

At the heart of the matter: the origin of bulgeless dwarf galaxies and Dark Matter cores

Authors: Fabio Governato, Chris Brook, Lucio Mayer, Alyson Brooks, George Rhee, James Wadsley, Patrik Jonsson, Beth Willman, Greg Stinson, Thomas Quinn, Piero Madau

Abstract: For almost two decades the properties of "dwarf" galaxies have challenged the Cold Dark Matter (CDM) paradigm of galaxy formation. Most observed dwarf galaxies consists of a rotating stellar disc embedded in a massive DM halo with a near constant-density core. Yet, models based on the CDM scenario invariably form galaxies with dense spheroidal stellar "bulges" and steep central DM profiles, as l… ▽ More For almost two decades the properties of "dwarf" galaxies have challenged the Cold Dark Matter (CDM) paradigm of galaxy formation. Most observed dwarf galaxies consists of a rotating stellar disc embedded in a massive DM halo with a near constant-density core. Yet, models based on the CDM scenario invariably form galaxies with dense spheroidal stellar "bulges" and steep central DM profiles, as low angular momentum baryons and DM sink to the center of galaxies through accretion and repeated mergers. Processes that decrease the central density of CDM halos have been identified, but have not yet reconciled theory with observations of present day dwarfs. This failure is potentially catastrophic for the CDM model, possibly requiring a different DM particle candidate. This Letter presents new hydrodynamical simulations in a Lambda$CDM framework where analogues of dwarf galaxies, bulgeless and with a shallow central DM profile, are formed. This is achieved by resolving the inhomogeneous interstellar medium, resulting in strong outflows from supernovae explosions which remove low angular momentum gas. This inhibits the formation of bulges and decreases the dark-matter density to less than half within the central kiloparsec. Realistic dwarf galaxies are thus shown to be a natural outcome of galaxy formation in the CDM scenario. △ Less

Submitted 11 November, 2009; originally announced November 2009.

Comments: To appear on Nature. Includes Supplemental Information. A movie can be downloaded at http://www.astro.washington.edu/users/fabio/movies/dwarfCDM2.mov

Journal ref: Nature, Volume 463, pp. 203-206 (2010)

A Physical Model for z~2 Dust Obscured Galaxies

Authors: Desika Narayanan, Arjun Dey, Christopher Hayward, Thomas J. Cox, R. Shane Bussmann, Mark Brodwin, Patrik Jonsson, Philip Hopkins, Brent Groves, Joshua D. Younger, Lars Hernquist

Abstract: We present a physical model for the origin of z~2 Dust-Obscured Galaxies (DOGs), a class of high-redshift ULIRGs selected at 24 micron which are particularly optically faint (24/R>1000). By combining N-body/SPH simulations of high redshift galaxy evolution with 3D polychromatic dust radiative transfer models, we find that luminous DOGs (with F24 > 0.3 mJy at z~2 are well-modeled as extreme gas-ric… ▽ More We present a physical model for the origin of z~2 Dust-Obscured Galaxies (DOGs), a class of high-redshift ULIRGs selected at 24 micron which are particularly optically faint (24/R>1000). By combining N-body/SPH simulations of high redshift galaxy evolution with 3D polychromatic dust radiative transfer models, we find that luminous DOGs (with F24 > 0.3 mJy at z~2 are well-modeled as extreme gas-rich mergers in massive (~5x10^12-10^13 Msun) halos, with elevated star formation rates (~500-1000 Msun/yr) and/or significant AGN growth (Mdot > 0.5 Msun/yr), whereas less luminous DOGs are more diverse in nature. At final coalescence, merger-driven DOGs transition from being starburst dominated to AGN dominated, evolving from a "bump" to a power-law shaped mid-IR (IRAC) spectral energy distribution (SED). After the DOG phase, the galaxy settles back to exhibiting a "bump" SED with bluer colors and lower star formation rates. While canonically power-law galaxies are associated with being AGN-dominated, we find that the power-law mid-IR SED can owe both to direct AGN contribution, as well as to a heavily dust obscured stellar bump at times that the galaxy is starburst dominated. Thus power-law galaxies can be either starburst or AGN dominated. Less luminous DOGs can be well-represented either by mergers, or by massive ($M_{\rm baryon} ~5x10^11 Msun) secularly evolving gas-rich disc galaxies (with SFR > 50 Msun/yr). By utilising similar models as those employed in the SMG formation study of Narayanan et al. (2010), we investigate the connection between DOGs and SMGs. We find that the most heavily star-forming merger driven DOGs can be selected as Submillimetre Galaxies (SMGs), while both merger-driven and secularly evolving DOGs typically satisfy the BzK selection criteria. △ Less

Submitted 11 May, 2010; v1 submitted 13 October, 2009; originally announced October 2009.

Comments: Accepted by MNRAS; major changes include better description of dependency on ISM specification and updated models allowing dust to evolve with metallicity.

Accelerating Dust Temperature Calculations with Graphics Processing Units

Authors: Patrik Jonsson, Joel Primack

Abstract: When calculating the infrared spectral energy distributions (SEDs) of galaxies in radiation-transfer models, the calculation of dust grain temperatures is generally the most time-consuming part of the calculation. Because of its highly parallel nature, this calculation is perfectly suited for massively parallel general-purpose Graphics Processing Units (GPUs). This paper presents an implementati… ▽ More When calculating the infrared spectral energy distributions (SEDs) of galaxies in radiation-transfer models, the calculation of dust grain temperatures is generally the most time-consuming part of the calculation. Because of its highly parallel nature, this calculation is perfectly suited for massively parallel general-purpose Graphics Processing Units (GPUs). This paper presents an implementation of the calculation of dust grain equilibrium temperatures on GPUs in the Monte-Carlo radiation transfer code Sunrise, using the CUDA API. The GPU can perform this calculation 69 times faster than the 8 CPU cores, showing great potential for accelerating calculations of galaxy SEDs. △ Less

Submitted 23 December, 2009; v1 submitted 21 July, 2009; originally announced July 2009.

Comments: 7 pages, 2 figures, accepted to New Astronomy. Minor updates to text and performance based on feedback from referee

High-Resolution Panchromatic Spectral Models of Galaxies including Photoionisation and Dust

Authors: Patrik Jonsson, Brent Groves, T. J. Cox

Abstract: An updated version of the dust radiation transfer code Sunrise, including models for star-forming regions and a self-consistent calculation of the spatially dependent dust and PAH emission, is presented. Given a hydrodynamic simulation of a galaxy, this model can calculate a realistic 2-dimensional ultraviolet--submillimeter spectral energy distribution of the galaxy, including emission lines fr… ▽ More An updated version of the dust radiation transfer code Sunrise, including models for star-forming regions and a self-consistent calculation of the spatially dependent dust and PAH emission, is presented. Given a hydrodynamic simulation of a galaxy, this model can calculate a realistic 2-dimensional ultraviolet--submillimeter spectral energy distribution of the galaxy, including emission lines from HII regions, from any viewpoint. To model the emission from star-forming regions, the MAPPINGSIII photoionization code is used. The high wavelength resolution (~ 1000 wavelengths) is made possible by the polychromatic Monte-Carlo algorithm employed by Sunrise. From the 2-D spectral energy distributions, images in any filter bands or integrated galaxy SEDs can be created. Using a suite of hydrodynamic simulations of disc galaxies, the output broad-band images and spectral energy distributions are compared with observed galaxies from the multiwavelength SINGS and SLUGS galaxy surveys. Overall, the output spectral energy distributions show a good match with observed galaxies in colours ranging from GALEX far-UV to SCUBA submillimeter wavelengths. The only possible exception is the 160 micron/850 micron colour, which the simulations underestimate by a factor "of order 5" compared to the SINGS sample. However, the simulations here agree with the SLUGS galaxies, which consistently have significantly larger amounts of cold dust than the SINGS galaxies. The Sunrise model can be used to generate simulated observations of arbitrary hydrodynamic galaxy simulations. In this way, predictions of galaxy formation theories can be directly tested against observations of galaxies. △ Less

Submitted 11 June, 2009; originally announced June 2009.

Comments: 32 pages, 22 figures, submitted to MNRAS

The merger-driven evolution of warm infrared luminous galaxies

Authors: J. D. Younger, C. C. Hayward, D. Narayanan, T. J. Cox, L. Hernquist, P. Jonsson

Abstract: We present a merger-driven evolutionary model for the production of luminous (LIRGs) and ultraluminous infrared galaxies (ULIRGs) with warm IR colours. Our results show that simulations of gas--rich major mergers including star formation, black hole growth, and feedback can produce warm (U)LIRGs. We also find that while the warm evolutionary phase is associated with increased AGN activity, star… ▽ More We present a merger-driven evolutionary model for the production of luminous (LIRGs) and ultraluminous infrared galaxies (ULIRGs) with warm IR colours. Our results show that simulations of gas--rich major mergers including star formation, black hole growth, and feedback can produce warm (U)LIRGs. We also find that while the warm evolutionary phase is associated with increased AGN activity, star formation alone may be sufficient to produce warm IR colours. However, the transition can be suppressed entirely - even when there is a significant AGN contribution - when we assume a single-phase ISM, which maximizes the attenuation. Finally, our evolutionary models are consistent with the 25-to-60 micron flux density ratio vs. L_HX/L_IR relation for local LIRGs and ULIRGs, and predict the observed scatter in IR color at fixed L_HX/L_IR. Therefore, our models suggest a cautionary note in the interpretation of warm IR colours: while associated with periods of active black hole growth, they are probably produced by a complex mix of star formation and AGN activity intermediate between the cold star formation dominated phase and the birth of a bright, unobscured quasar. △ Less

Submitted 1 April, 2009; originally announced April 2009.

Comments: 5 pages, accepted for publication in MNRAS Letters

The Formation of High Redshift Submillimeter Galaxies

Authors: Desika Narayanan, Christopher C. Hayward, Thomas J. Cox, Lars Hernquist, Patrik Jonsson, Joshua D. Younger, Brent Groves

Abstract: We describe a model for the formation of \zsim 2 Submillimeter Galaxies (SMGs) which simultaneously accounts for both average and bright SMGs while providing a reasonable match to their mean observed spectral energy distributions (SEDs). By coupling hydrodynamic simulations of galaxy mergers with the high resolution 3D polychromatic radiative transfer code Sunrise, we find that a mass sequence o… ▽ More We describe a model for the formation of \zsim 2 Submillimeter Galaxies (SMGs) which simultaneously accounts for both average and bright SMGs while providing a reasonable match to their mean observed spectral energy distributions (SEDs). By coupling hydrodynamic simulations of galaxy mergers with the high resolution 3D polychromatic radiative transfer code Sunrise, we find that a mass sequence of merger models which use observational constraints as physical input naturally yield objects which exhibit black hole, bulge, and H2 gas masses similar to those observed in SMGs. The dominant drivers behind the 850 micron flux are the masses of the merging galaxies and the stellar birthcloud covering fraction. The most luminous (S850 ~ 15 mJy) sources are recovered by ~10^13 Msun 1:1 major mergers with a birthcloud covering fraction close to unity, whereas more average SMGs ~5-7 mJy) may be formed in lower mass halos ~5x10^12 Msun. These models demonstrate the need for high spatial resolution hydrodynamic and radiative transfer simulations in matching both the most luminous sources as well as the full SEDs of SMGs. While these models suggest a natural formation mechanism for SMGs, they do not attempt to match cosmological statistics of galaxy populations; future efforts along this line will help ascertain the robustness of these models. △ Less

Submitted 8 October, 2009; v1 submitted 31 March, 2009; originally announced April 2009.

Comments: MNRAS Accepted; Revised version includes expanded discussion of simulated radio properties of SMGs

Astronomical Software Wants To Be Free: A Manifesto

Authors: Benjamin J. Weiner, Michael R. Blanton, Alison L. Coil, Michael C. Cooper, Romeel Davé, David W. Hogg, Bradford P. Holden, Patrik Jonsson, Susan A. Kassin, Jennifer M. Lotz, John Moustakas, Jeffrey A. Newman, J. X. Prochaska, Peter J. Teuben, Christy A. Tremonti, Christopher N. A. Willmer

Abstract: Astronomical software is now a fact of daily life for all hands-on members of our community. Purpose-built software for data reduction and modeling tasks becomes ever more critical as we handle larger amounts of data and simulations. However, the writing of astronomical software is unglamorous, the rewards are not always clear, and there are structural disincentives to releasing software publicl… ▽ More Astronomical software is now a fact of daily life for all hands-on members of our community. Purpose-built software for data reduction and modeling tasks becomes ever more critical as we handle larger amounts of data and simulations. However, the writing of astronomical software is unglamorous, the rewards are not always clear, and there are structural disincentives to releasing software publicly and to embedding it in the scientific literature, which can lead to significant duplication of effort and an incomplete scientific record. We identify some of these structural disincentives and suggest a variety of approaches to address them, with the goals of raising the quality of astronomical software, improving the lot of scientist-authors, and providing benefits to the entire community, analogous to the benefits provided by open access to large survey and simulation datasets. Our aim is to open a conversation on how to move forward. We advocate that: (1) the astronomical community consider software as an integral and fundable part of facility construction and science programs; (2) that software release be considered as integral to the open and reproducible scientific process as are publication and data release; (3) that we adopt technologies and repositories for releasing and collaboration on software that have worked for open-source software; (4) that we seek structural incentives to make the release of software and related publications easier for scientist-authors; (5) that we consider new ways of funding the development of grass-roots software; (6) and that we rethink our values to acknowledge that astronomical software development is not just a technical endeavor, but a fundamental part of our scientific practice. △ Less

Submitted 23 March, 2009; originally announced March 2009.

Comments: State of the Profession white paper submitted to the Astro2010 Decadal Survey. The authors welcome comments and discussion

Low-Energy Astrophysics: Stimulating the Reduction of Energy Consumption in the Next Decade

Authors: P. J. Marshall, N. Bennert, E. S. Rykoff, K. J. Shen, J. D. R. Steinfadt, J. Fregeau, R-R. Chary, K. Sheth, B. Weiner, K. B. Henisey, E. L. Quetin, R. Antonucci, D. Kaplan, P. Jonsson, M. W. Auger, C. Cardamone, T. Tao, D. E. Holz, M. Bradac, T. S. Metcalfe, S. McHugh, M. Elvis, B. J. Brewer, T. Urrutia, F. Guo , et al. (5 additional authors not shown)

Abstract: In this paper we address the consumption of energy by astronomers while performing their professional duties. Although we find that astronomy uses a negligible fraction of the US energy budget, the rate at which energy is consumed by an average astronomer is similar to that of a typical high-flying businessperson. We review some of the ways in which astronomers are already acting to reduce their… ▽ More In this paper we address the consumption of energy by astronomers while performing their professional duties. Although we find that astronomy uses a negligible fraction of the US energy budget, the rate at which energy is consumed by an average astronomer is similar to that of a typical high-flying businessperson. We review some of the ways in which astronomers are already acting to reduce their energy consumption. In the coming decades, all citizens will have to reduce their energy consumption to conserve fossil fuel reserves and to help avert a potentially catastrophic change in the Earth's climate. The challenges are the same for astronomers as they are for everyone: decreasing the distances we travel and investing in energy-efficient infrastructure. The high profile of astronomy in the media, and the great public interest in our field, can play a role in promoting energy-awareness to the wider population. Our specific recommendations are therefore to 1) reduce travel when possible, through efficient meeting organization, and by investing in high-bandwidth video conference facilities and virtual-world software, 2) create energy-efficient observatories, computing centers and workplaces, powered by sustainable energy resources, and 3) actively publicize these pursuits. △ Less

Submitted 19 March, 2009; originally announced March 2009.

Comments: "State of the Profession" white paper submitted to the 2010 Astronomy & Astrophysics Decadal Survey. Add your thoughts to the wiki at http://low-energy-astro.physics.ucsb.edu

Evolution of the Stellar Mass Tully-Fisher Relation in Disk Galaxy Merger Simulations

Authors: M. D. Covington, S. A. Kassin, A. A. Dutton, B. J. Weiner, T. J. Cox, P. Jonsson, J. R. Primack, S. M. Faber, D. C. Koo

Abstract: There is a large observational scatter toward low velocities in the stellar mass Tully-Fisher relation if disturbed and compact objects are included. However, this scatter can be eliminated if one replaces rotation velocity with $\rm S_{\rm 0.5}$, a quantity that includes a velocity dispersion term added in quadrature with the rotation velocity. In this work we use a large suite of hydrodynamic… ▽ More There is a large observational scatter toward low velocities in the stellar mass Tully-Fisher relation if disturbed and compact objects are included. However, this scatter can be eliminated if one replaces rotation velocity with $\rm S_{\rm 0.5}$, a quantity that includes a velocity dispersion term added in quadrature with the rotation velocity. In this work we use a large suite of hydrodynamic N-body galaxy merger simulations to explore a possible mechanism for creating the observed relations. Using mock observations of the simulations, we test for the presence of observational effects and explore the relationship between $\rm S_{\rm 0.5}$ and intrinsic properties of the galaxies. We find that galaxy mergers can explain the scatter in the TF as well as the tight $\rm S_{\rm 0.5}$-stellar mass relation. Furthermore, $\rm S_{\rm 0.5}$ is correlated with the total central mass of a galaxy, including contributions due to dark matter. △ Less

Submitted 3 December, 2009; v1 submitted 3 February, 2009; originally announced February 2009.

Comments: ApJ accepted

Journal ref: Astrophys.J.710:279-288,2010

Forming a Large Disc Galaxy from a z<1 Major Merger

Authors: F. Governato, C. B. Brook, A. M. Brooks, L. Mayer, B. Willman, P. Jonsson, A. M. Stilp, L. Pope, C. Christensen, J. Wadsley, T. Quinn

Abstract: Using high resolution SPH simulations in a fully cosmological Lambda CDM context we study the formation of a bright disk dominated galaxy that originates from a "wet" major merger at z=0.8. The progenitors of the disk galaxy are themselves disk galaxies that formed from early major mergers between galaxies with blue colors. A substantial thin stellar disk grows rapidly following the last major m… ▽ More Using high resolution SPH simulations in a fully cosmological Lambda CDM context we study the formation of a bright disk dominated galaxy that originates from a "wet" major merger at z=0.8. The progenitors of the disk galaxy are themselves disk galaxies that formed from early major mergers between galaxies with blue colors. A substantial thin stellar disk grows rapidly following the last major merger and the present day properties of the final remnant are typical of early type spiral galaxies, with an i band B/D ~0.65, a disk scale length of 7.2 kpc, g-r = 0.5 mag, an HI line width (W_{20}/2) of 238 km/sec and total magnitude i = -22.4. The key ingredients for the formation of a dominant stellar disk component after a major merger are: i) substantial and rapid accretion of gas through cold flows followed at late times by cooling of gas from the hot phase, ii) supernova feedback that is able to partially suppress star formation during mergers and iii) relative fading of the spheroidal component. The gas fraction of the progenitors' disks does not exceed 25% at z<3, emphasizing that the continuous supply of gas from the local environment plays a major role in the regrowth of disks and in keeping the galaxies blue. The results of this simulation alleviate the problem posed for the existence of disk galaxies by the high likelihood of interactions and mergers for galaxy sized halos at relatively low z. △ Less

Submitted 20 July, 2009; v1 submitted 1 December, 2008; originally announced December 2008.

Comments: MNRAS in press. Minor additions to text. Movie at: http://www.astro.washington.edu/fabio/movies/Merger.mpg

Journal ref: MNRAS, 398, pp. 312-320 (2009)

Galaxy Merger Morphologies and Time-Scales from Simulations of Equal-Mass Gas-Rich Disc Mergers

Authors: Jennifer M. Lotz, Patrik Jonsson, T. J. Cox, Joel R. Primack

Abstract: A key obstacle to understanding the galaxy merger rate and its role in galaxy evolution is the difficulty in constraining the merger properties and time-scales from instantaneous snapshots of the real universe.The most common way to identify galaxy mergers is by morphology, yet current theoretical calculations of the time-scales for galaxy disturbances are quite crude. We present a morphological… ▽ More A key obstacle to understanding the galaxy merger rate and its role in galaxy evolution is the difficulty in constraining the merger properties and time-scales from instantaneous snapshots of the real universe.The most common way to identify galaxy mergers is by morphology, yet current theoretical calculations of the time-scales for galaxy disturbances are quite crude. We present a morphological analysis of a large suite of GADGET N-Body/hydro-dynamical equal-mass gas-rich disc galaxy mergers which have been processed through the Monte-Carlo radiative transfer code SUNRISE. With the resulting images, we examine the dependence of quantitative morphology (G, M20, C, A) in the SDSS g-band on merger stage, dust, viewing angle, orbital parameters, gas properties, supernova feedback, and total mass. We find that mergers appear most disturbed in G-M20 and asymmetry at the first pass and at the final coalescence of their nuclei, but can have normal quantitative morphologies at other merger stages. The merger observability time-scales depend on the method used to identify the merger as well as the gas fraction, pericentric distance, and relative orientation of the merging galaxies. Enhanced star formation peaks after and lasts significantly longer than strong morphological disturbances. Despite their massive bulges, the majority of merger remnants appear disc-like and dusty in g-band light because of the presence of a low-mass star-forming disc. △ Less

Submitted 25 September, 2008; v1 submitted 8 May, 2008; originally announced May 2008.

Comments: accepted to MNRAS; 29 pages, 19 color figures; go to http://www.noao.edu/staff/lotz/pdfs/lotz_arXiv.0805.1246.pdf for full resolution

The Formation of Polar Disk Galaxies

Authors: Chris B. Brook, Fabio Governato, Thomas Quinn, James Wadsley, Alyson M. Brooks, Beth Willman, Adrienne Stilp, Patrik Jonsson

Abstract: Polar Ring Galaxies, such as NGC4650A, are a class of galaxy which have two kinematically distinct components that are inclined by almost 90 degrees to each other. These striking galaxies challenge our understanding of how galaxies form; the origin of their distinct components has remained uncertain, and the subject of much debate. We use high-resolution cosmological simulations of galaxy format… ▽ More Polar Ring Galaxies, such as NGC4650A, are a class of galaxy which have two kinematically distinct components that are inclined by almost 90 degrees to each other. These striking galaxies challenge our understanding of how galaxies form; the origin of their distinct components has remained uncertain, and the subject of much debate. We use high-resolution cosmological simulations of galaxy formation to show that Polar Ring Galaxies are simply an extreme example of the angular moment misalignment that occurs during the hierarchical structure formation characteristic of Cold Dark Matter cosmology. In our model, Polar Ring Galaxies form through the continuous accretion of gas whose angular momentum is misaligned with the central galaxy. △ Less

Submitted 8 February, 2008; v1 submitted 7 February, 2008; originally announced February 2008.

Comments: Submitted to ApJ. 8 pages in emulate ApJ style. 2 associated animations are found at: http://www.youtube.com/watch?v=c-H3WzaewdY http://www.youtube.com/watch?v=9Xf3fJkgWEg

Predicting the Properties of the Remnants of Dissipative Galaxy Mergers

Authors: M. Covington, A. Dekel, T. J. Cox, P. Jonsson, J. R. Primack

Abstract: We construct a physically motivated model for predicting the properties of the remnants of gaseous galaxy mergers, given the properties of the progenitors and the orbit. The model is calibrated using a large suite of SPH merger simulations. It implements generalized energy conservation while accounting for dissipative energy losses and star formation. The dissipative effects are evaluated from t… ▽ More We construct a physically motivated model for predicting the properties of the remnants of gaseous galaxy mergers, given the properties of the progenitors and the orbit. The model is calibrated using a large suite of SPH merger simulations. It implements generalized energy conservation while accounting for dissipative energy losses and star formation. The dissipative effects are evaluated from the initial gas fractions and from the orbital parameters via an "impulse" parameter, which characterizes the strength of the encounter. Given the progenitor properties, the model predicts the remnant stellar mass, half-mass radius, and velocity dispersion to an accuracy of 25%. The model is valid for both major and minor mergers. We provide an explicit recipe for semi-analytic models of galaxy formation. △ Less

Submitted 24 October, 2007; originally announced October 2007.

Comments: 14 pages, 9 figures. MNRAS accepted

The effect of galaxy mass ratio on merger--driven starbursts

Authors: T. J. Cox, Patrik Jonsson, Rachel S. Somerville, Joel R. Primack, Avishai Dekel

Abstract: We employ numerical simulations of galaxy mergers to explore the effect of galaxy mass ratio on merger--driven starbursts. Our numerical simulations include radiative cooling of gas, star formation, and stellar feedback to follow the interaction and merger of four disk galaxies. The galaxy models span a factor of 23 in total mass and are designed to be representative of typical galaxies in the l… ▽ More We employ numerical simulations of galaxy mergers to explore the effect of galaxy mass ratio on merger--driven starbursts. Our numerical simulations include radiative cooling of gas, star formation, and stellar feedback to follow the interaction and merger of four disk galaxies. The galaxy models span a factor of 23 in total mass and are designed to be representative of typical galaxies in the local Universe. We find that the merger--driven star formation is a strong function of merger mass ratio, with very little, if any, induced star formation for large mass ratio mergers. We define a burst efficiency that is useful to characterize the merger--driven star formation and test that it is insensitive to uncertainties in the feedback parameterization. In accord with previous work we find that the burst efficiency depends on the structure of the primary galaxy. In particular, the presence of a massive stellar bulge stabilizes the disk and suppresses merger--driven star formation for large mass ratio mergers. Direct, co--planar merging orbits produce the largest tidal disturbance and yield that most intense burst of star formation. Contrary to naive expectations, a more compact distribution of gas or an increased gas fraction both decrease the burst efficiency. Owing to the efficient feedback model and the newer version of SPH employed here, the burst efficiencies of the mergers presented here are smaller than in previous studies. △ Less

Submitted 21 September, 2007; originally announced September 2007.

Comments: 26 pages, 21 figures, submitted to MNRAS

Dust Attenuation in Hydrodynamic Simulations of Spiral Galaxies

Authors: Miguel Rocha, Patrik Jonsson, Joel R. Primack, T. J. Cox

Abstract: We study the effects of dust in hydrodynamic simulations of spiral galaxies when different radial metallicity gradients are assumed. SUNRISE, a Monte-Carlo radiative-transfer code, is used to make detailed calculations of the internal extinction of disk galaxies caused by their dust content. SUNRISE is used on eight different Smooth Particle Hydrodynamics (SPH) simulations of isolated spiral g… ▽ More We study the effects of dust in hydrodynamic simulations of spiral galaxies when different radial metallicity gradients are assumed. SUNRISE, a Monte-Carlo radiative-transfer code, is used to make detailed calculations of the internal extinction of disk galaxies caused by their dust content. SUNRISE is used on eight different Smooth Particle Hydrodynamics (SPH) simulations of isolated spiral galaxies. These galaxies vary mainly in mass and hence luminosity, spanning a range in luminosities from -16 to -22 magnitudes in the B band. We focus on the attenuation in different wavelength bands as a function of the disk inclination and the luminosity of the models, and compare this to observations. Observations suggest different metallicity gradients for galaxies of different luminosities. These metallicity gradients were explored in our different models, finding that the resulting dust attenuation matches observations for edge-on galaxies, but do not show a linear behaviour in log axis ratio as some observations have suggested. A quadratic law describing the dependence of attenuation on inclination, as proposed by more recent observations, reconciles the attenuation of the simulations at intermediate inclinations with observations. We also compare the total infrared-to-ultraviolet flux ratios for the simulated galaxies with those of the SINGS sample and find general agreement. Finally we compare our results with those from simpler models that do not take into account structure such as spiral arms, finding that the inclusion of sub-structure on the size scale of spiral arms does not change conclusions about the attenuation dependence on inclination or wavelength. △ Less

Submitted 1 February, 2008; v1 submitted 20 February, 2007; originally announced February 2007.

Comments: 12 pages, 8 figures, all figures require colors, accepted for publication in MNRAS. Major changes from v1 include: Comparison of the calculated infrared-to-ultraviolet ratios in our models to those observed by SIGNS, clearer fonts in our plots and clearer justification of the rationale behind the adopted gradients. Changes from v2: Added journal-ref

Journal ref: MNRAS. 383 (2008), Issue 4, pp. 1281-1291

Shapes of Stellar Systems and Dark Halos from Simulations of Galaxy Major Mergers

Authors: Gregory S. Novak, Thomas J. Cox, Joel R. Primack, Patrik Jonsson, Avishai Dekel

Abstract: Using a sample of 89 snapshots from 58 hydrodynamic binary galaxy major merger simulations, we find that stellar remnants are mostly oblate while dark matter halos are mostly prolate or triaxial. The stellar minor axis and the halo major axis are almost always nearly perpendicular. This can be understood by considering the influence of angular momentum and dissipation during the merger. If binar… ▽ More Using a sample of 89 snapshots from 58 hydrodynamic binary galaxy major merger simulations, we find that stellar remnants are mostly oblate while dark matter halos are mostly prolate or triaxial. The stellar minor axis and the halo major axis are almost always nearly perpendicular. This can be understood by considering the influence of angular momentum and dissipation during the merger. If binary mergers of spiral galaxies are responsible for the formation of elliptical galaxies or some subpopulation thereof, these galaxies can be expected to be oblate and inhabit their halos with the predicted shapes and orientations. These predictions are potentially relevant to observational studies of weak gravitational lensing, where one must stack many optically aligned galaxies in order to determine the shape of the resulting stacked mass distribution. The simple relationship between the dark and luminous matter presented here can be used to guide the stacking of galaxies to minimize the information lost. △ Less

Submitted 21 July, 2006; v1 submitted 6 April, 2006; originally announced April 2006.

Comments: 4 pages, 5 figures. Minor changes to match published version

Journal ref: Astrophys.J.646:L9-L12,2006

Sunrise: Polychromatic Dust Radiative Transfer in Arbitrary Geometries

Authors: Patrik Jonsson

Abstract: This paper describes Sunrise, a parallel, free Monte-Carlo code for the calculation of radiation transfer through astronomical dust. Sunrise uses an adaptive-mesh refinement grid to describe arbitrary geometries of emitting and absorbing/scattering media, with spatial dynamical range exceeding 10^4, and it can efficiently generate images of the emerging radiation at arbitrary points in space. In… ▽ More This paper describes Sunrise, a parallel, free Monte-Carlo code for the calculation of radiation transfer through astronomical dust. Sunrise uses an adaptive-mesh refinement grid to describe arbitrary geometries of emitting and absorbing/scattering media, with spatial dynamical range exceeding 10^4, and it can efficiently generate images of the emerging radiation at arbitrary points in space. In addition to the monochromatic radiative transfer typically used by Monte-Carlo codes, Sunrise is capable of propagating a range of wavelengths simultaneously. This ``polychromatic'' algorithm gives significant improvements in efficiency and accuracy when spectral features are calculated. Sunrise is used to study the effects of dust in hydrodynamic simulations of interacting galaxies, and the procedure for this is described. The code is tested against previously published results. △ Less

Submitted 21 July, 2006; v1 submitted 5 April, 2006; originally announced April 2006.

Comments: 22 pages, 11 figures, submitted to MNRAS. Changed based on referee response. Section on benchmarks expanded

Journal ref: Mon.Not.Roy.Astron.Soc.372:2-20,2006

Feedback in simulations of disc-galaxy major mergers

Authors: T. J. Cox, Patrik Jonsson, Joel R. Primack, Rachel S. Somerville

Abstract: Using hydrodynamic simulations of disc-galaxy major mergers, we investigate the star formation history and remnant properties when various parametrizations of a simple stellar feedback model are implemented. The simulations include radiative cooling, a density-dependent star formation recipe and a model for feedback from massive stars. The feedback model stores supernova feedback energy within i… ▽ More Using hydrodynamic simulations of disc-galaxy major mergers, we investigate the star formation history and remnant properties when various parametrizations of a simple stellar feedback model are implemented. The simulations include radiative cooling, a density-dependent star formation recipe and a model for feedback from massive stars. The feedback model stores supernova feedback energy within individual gas particles and dissipates this energy on a time-scale specified by two free parameters; tau_fb, which sets the dissipative time-scale, and n, which sets the effective equation of state in star-forming regions. Using a self-consistent disc galaxy, modelled after a local Sbc spiral, in both isolated and major-merger simulations, we investigate parametrizations of the feedback model that are selected with respect to the quiescent disc stability. These models produce a range of star formation histories that are consistent with the star formation relation found by Kennicutt. All major mergers produce a population of new stars that is highly centrally concentrated, demonstrating a distinct break in the r1/4 surface density profile, consistent with previous findings. The half-mass radius and one-dimensional velocity dispersion are affected by the feedback model used. Finally, we compare our results to those of previous simulations of star formation in disc-galaxy major mergers, addressing the effects of star formation normalization, the version of smoothed particle hydrodynamics (SPH) employed and assumptions about the interstellar medium. △ Less

Submitted 25 September, 2007; v1 submitted 9 March, 2005; originally announced March 2005.

Comments: update to match published version

Journal ref: 2006 MNRAS, 373, 1013

Simulations of Dust in Interacting Galaxies I: Dust Attenuation

Authors: Patrik Jonsson, T. J. Cox, Joel R. Primack, Rachel S. Somerville

Abstract: A new Monte-Carlo radiative-transfer code, Sunrise, is used in conjunction with hydrodynamic simulations of major galaxy mergers to calculate the effects of dust in such systems. The simulations are in good agreement with observations of dust absorption in starburst galaxies, and the dust has a profound effect on their appearance. The dust attenuation increases with luminosity such that at peak… ▽ More A new Monte-Carlo radiative-transfer code, Sunrise, is used in conjunction with hydrodynamic simulations of major galaxy mergers to calculate the effects of dust in such systems. The simulations are in good agreement with observations of dust absorption in starburst galaxies, and the dust has a profound effect on their appearance. The dust attenuation increases with luminosity such that at peak luminosities ~90% of the bolometric luminosity is absorbed by dust. In general, the detailed appearance of the merging event depends on the stage of the merger and the geometry of the encounter. The fraction of bolometric energy absorbed by the dust, however, is a robust quantity that can be predicted from the intrinsic properties bolometric luminosity, baryonic mass, star-formation rate, and metallicity of the system. This paper presents fitting formulae, valid over a wide range of masses and metallicities, from which the absorbed fraction of luminosity (and consequently also the infrared dust luminosity) can be predicted. The attenuation of the luminosity at specific wavelengths can also be predicted, albeit with a larger scatter due to the variation with viewing angle. These formulae for dust attenuation appear to be valid for both isolated and interacting galaxies, are consistent with earlier studies, and would be suitable for inclusion in theoretical models, e.g. semi-analytic models of galaxy formation. △ Less

Submitted 6 March, 2005; originally announced March 2005.

Comments: 12 pages, 10 figures, submitted to ApJ

Journal ref: Astrophys.J.637:255-268,2006

Simulations of Dust in Interacting Galaxies

Authors: Patrik Jonsson, T. J. Cox, Joel R. Primack

Abstract: A new Monte-Carlo radiative-transfer code, Sunrise, is used to study the effects of dust in N-body/hydrodynamic simulations of interacting galaxies. Dust has a profound effect on the appearance of the simulated galaxies. At peak luminosities, about 90% of the bolometric luminosity is absorbed, and the dust obscuration scales with luminosity in such a way that the brightness at UV/visual waveleng… ▽ More A new Monte-Carlo radiative-transfer code, Sunrise, is used to study the effects of dust in N-body/hydrodynamic simulations of interacting galaxies. Dust has a profound effect on the appearance of the simulated galaxies. At peak luminosities, about 90% of the bolometric luminosity is absorbed, and the dust obscuration scales with luminosity in such a way that the brightness at UV/visual wavelengths remains roughly constant. A general relationship between the fraction of energy absorbed and the ratio of bolometric luminosity to baryonic mass is found. Comparing to observations, the simulations are found to follow a relation similar to the observed IRX-Beta relation found by Meurer et al (1999) when similar luminosity objects are considered. The highest-luminosity simulated galaxies depart from this relation and occupy the region where local (U)LIRGs are found. This agreement is contingent on the presence of Milky-Way-like dust, while SMC-like dust results in far too red a UV continuum slope to match observations. The simulations are used to study the performance of star-formation indicators in the presence of dust. The far-infrared luminosity is found to be reliable. In contrast, the H-alpha and far-UV luminosity suffer severely from dust attenuation, and dust corrections can only partially remedy the situation. △ Less

Submitted 5 December, 2004; originally announced December 2004.

Comments: 4 pages, 5 figures, to appear in the proceedings of the conference "The Spectral Energy Distribution of Gas-Rich Galaxies", eds. C.C. Popescu & R.J. Tuffs (Heidelberg, October 2004)

Journal ref: AIP Conf.Proc.761:320-323,2005

Generating Hot Gas in Simulations of Disk-Galaxy Major Mergers

Authors: T. J. Cox, Joel Primack, Patrik Jonsson, Rachel Somerville

Abstract: We report on the merger-induced generation of a shock-heated gas wind and formation of a remnant gas halo in simulations of colliding disk galaxies. The simulations use cosmologically motivated initial conditions and include the effects of radiative cooling, star formation, stellar feedback and the non-adiabatic heating of gas. The non-adiabatic heating, i.e. shocks, generated in the final merge… ▽ More We report on the merger-induced generation of a shock-heated gas wind and formation of a remnant gas halo in simulations of colliding disk galaxies. The simulations use cosmologically motivated initial conditions and include the effects of radiative cooling, star formation, stellar feedback and the non-adiabatic heating of gas. The non-adiabatic heating, i.e. shocks, generated in the final merger forces gas out of the central region of the merger remnant and into the dark-matter halo. We demonstrate that the amount of heating depends on the size of the progenitor disk galaxy as well as the initial orbit the galaxies are placed on. Based upon these dependencies, we motivate a possible recipe for including this effect in semi-analytic models of galaxy formation. △ Less

Submitted 27 February, 2004; originally announced February 2004.

Comments: 4 pages, 2 figures, uses emulateapj.cls, submitted to ApJ Letters

Journal ref: Astrophys.J.607:L87-L90,2004

Young galaxies: what turns them on?

Authors: Tsafrir S. Kolatt, James S. Bullock, Rachel S. Somerville, Yair Sigad, Patrik Jonsson, Andrey V. Kravtsov, Anatoly A. Klypin, Joel R. Primack, Sandra M. Faber, Avishai Dekel, .

Abstract: Lyman break galaxies (LBGs) at $z \sim 3$ exhibit number densities and clustering similar to local $L_*$ galaxies; however, their star formation rates (SFRs) are much higher. We explore the scenario in which LBGs are starburst galaxies triggered by collisions, and thus provide an explanation for these key properties. The relative paucity of starburst galaxies at low redshift can be attributed to… ▽ More Lyman break galaxies (LBGs) at $z \sim 3$ exhibit number densities and clustering similar to local $L_*$ galaxies; however, their star formation rates (SFRs) are much higher. We explore the scenario in which LBGs are starburst galaxies triggered by collisions, and thus provide an explanation for these key properties. The relative paucity of starburst galaxies at low redshift can be attributed to a much higher collision rate in the past. We use high-resolution cosmological $N$-body simulations and a hierarchical halo finder to estimate the galaxy collision rate as a function of time in the popular $Λ$CDM cosmological model. We find that bright collisional starbursts are frequent enough to account for most of the high-$z$ ($2.5-4.5$) LBGs. Moreover, many of the objects are of relatively small mass, but they cluster about large-mass halos. They therefore exhibit strong clustering, similar to that observed and stronger than that of the relevant massive halos. △ Less

Submitted 7 June, 1999; originally announced June 1999.

Comments: 7 pages, three figures (one jpg), submitted to the Astrophysical Journal Letters