-
Massive Star Cluster Formation II. Runaway Stars as Fossils of Sub-Cluster Mergers
Authors:
Brooke Polak,
Mordecai-Mark Mac Low,
Ralf S. Klessen,
Simon Portegies Zwart,
Eric P. Andersson,
Sabrina M. Appel,
Claude Cournoyer-Cloutier,
Simon C. O. Glover,
Stephen L. W. McMillan
Abstract:
Two main mechanisms have classically been proposed for the formation of runaway stars. The binary supernova scenario (BSS) suggests that a massive star in a binary explodes as a supernova, ejecting its companion. The dynamical ejection scenario suggests that a star is ejected through a strong dynamical encounter between multiple stars. We propose a third mechanism for the formation of runaway star…
▽ More
Two main mechanisms have classically been proposed for the formation of runaway stars. The binary supernova scenario (BSS) suggests that a massive star in a binary explodes as a supernova, ejecting its companion. The dynamical ejection scenario suggests that a star is ejected through a strong dynamical encounter between multiple stars. We propose a third mechanism for the formation of runaway stars: the sub-cluster ejection scenario (SCES), where an infalling sub-cluster of stars is ejected out of the cluster by a slingshot interaction with the contracting gravitational potential of the assembling cluster. We demonstrate the SCES in a star-by-star simulation of a young massive cluster forming from a $10^6\rm~M_\odot$ gas cloud using the Torch framework. This star cluster forms hierarchically through a sequence of sub-cluster mergers, determined by the initial turbulent, spherical initial conditions of the gas. We find that these mergers drive the formation of runaway stars in our model. Late-forming sub-clusters fall into the central potential, where they are ejected on a slingshot trajectory, forming groups of runaway stars that are distributed highly anisotropically. Runaways formed by the same SCES share similar ages, velocities, and ejection directions. Surveying observations, we identify several SCES candidate groups with anisotropic ejection directions. The SCES is capable of producing runaway binaries: two wide dynamical binaries in infalling sub-clusters were tightened through ejection. This allows for another velocity kick via subsequent BSS, which is a promising avenue for producing hypervelocity stars unbound to the Galaxy. The SCES occurs when sub-cluster formation is resolved. We expect non-spherical initial gas distributions to increase runaway star numbers. The observation of groups of runaway stars formed via SCES thus reveals the assembly history of their natal clusters.
△ Less
Submitted 20 May, 2024;
originally announced May 2024.
-
Massive Star Cluster Formation I. High Star Formation Efficiency While Resolving Feedback of Individual Stars
Authors:
Brooke Polak,
Mordecai-Mark Mac Low,
Ralf S. Klessen,
Jia Wei Teh,
Claude Cournoyer-Cloutier,
Eric P. Andersson,
Sabrina M. Appel,
Aaron Tran,
Sean C. Lewis,
Maite J. C. Wilhelm,
Simon Portegies Zwart,
Simon C. O. Glover,
Long Wang,
Stephen L. W. McMillan
Abstract:
The mode of star formation that results in the formation of globular clusters and young massive clusters is difficult to constrain through observations. We present models of massive star cluster formation using the Torch framework, which uses AMUSE to couple distinct multi-physics codes that handle star formation, stellar evolution and dynamics, radiative transfer, and magnetohydrodynamics. We upg…
▽ More
The mode of star formation that results in the formation of globular clusters and young massive clusters is difficult to constrain through observations. We present models of massive star cluster formation using the Torch framework, which uses AMUSE to couple distinct multi-physics codes that handle star formation, stellar evolution and dynamics, radiative transfer, and magnetohydrodynamics. We upgrade Torch by implementing the N-body code PeTar, thereby enabling Torch to handle massive clusters forming from $10^6\rm\, M_\odot$ clouds with $\ge10^5$ individual stars. We present results from Torch simulations of star clusters forming from $10^4, 10^5$, and $10^6\rm M_\odot$ turbulent, spherical gas clouds (named M4, M5, M6) of radius $R=11.7$ pc. We find that star formation is highly efficient and becomes more so at higher cloud mass and surface density. For M4, M5, and M6 with initial surface densities $2.325\times 10^{1,2,3}\rm\, M_\odot\, pc^{-2}$, after a free-fall time of $t_{ff}=6.7,2.1,0.67$ Myr, we find that $\sim30\%$, 40%, and 60% of the cloud mass has formed into stars, respectively. The final integrated star formation efficiency is $32\%,\, 65\%$, and 85\% for M4, M5, and M6. Observations of nearby clusters similar to M4 have similar integrated star formation efficiencies of $\leq30\%$. The M5 and M6 models represent a different regime of cluster formation that is more appropriate for the conditions in starburst galaxies and gas-rich galaxies at high redshift, and that leads to a significantly higher efficiency of star formation. We argue that young massive clusters build up through short efficient bursts of star formation in regions that are sufficiently dense ($\ge 10^2 \rm\,M_\odot\,pc^{-2}$) and massive ($\ge10^5\rm\, M_\odot$). In such environments, the dynamical time of the cloud becomes short enough that stellar feedback cannot act quickly enough to slow star formation.
△ Less
Submitted 11 December, 2023;
originally announced December 2023.
-
Pre-supernova feedback sets the star cluster mass function to a power law and reduces the cluster formation efficiency
Authors:
Eric P. Andersson,
Mordecai-Mark Mac Low,
Oscar Agertz,
Florent Renaud,
Hui Li
Abstract:
The star cluster initial mass function is observed to have an inverse power law exponent around 2, yet there is no consensus on what determines this distribution, and why some variation is observed in different galaxies. Furthermore, the cluster formation efficiency covers a range of values, particularly when considering different environments. These clusters are often used to empirically constrai…
▽ More
The star cluster initial mass function is observed to have an inverse power law exponent around 2, yet there is no consensus on what determines this distribution, and why some variation is observed in different galaxies. Furthermore, the cluster formation efficiency covers a range of values, particularly when considering different environments. These clusters are often used to empirically constrain star formation and as fundamental units for stellar feedback models. Detailed galaxy models must therefore accurately capture the basic properties of observed clusters to be considered predictive. We use hydrodynamical simulations of a dwarf galaxy as a laboratory to study star cluster formation. We test different combinations of stellar feedback mechanisms, including stellar winds, ionizing radiation, and supernovae. Each feedback mechanism affects the cluster formation efficiency and cluster mass function. Increasing the feedback budget by combining the different types of feedback decreases the cluster formation efficiency by reducing the number of massive clusters. Ionizing radiation is found to be especially influential. This effect depends on the timing of feedback initiation, as shown by comparing early and late feedback. Early feedback occurs from ionizing radiation and stellar winds with onset immediately after a massive star is formed. Late feedback occurs when energy injection only starts after the main-sequence lifetime of the most massive SN progenitor, a timing that is further influenced by the choice of the most massive SN progenitor. Late feedback alone results in a broad, flat mass function, approaching a log-normal shape in the complete absence of feedback. Early feedback, on the other hand, produces a power-law cluster mass function with lower formation efficiency, albeit with a steeper slope than that usually observed.
△ Less
Submitted 23 August, 2023;
originally announced August 2023.
-
EDGE: The direct link between mass growth history and the extended stellar haloes of the faintest dwarf galaxies
Authors:
Alex Goater,
Justin I. Read,
Noelia E. D. Noël,
Matthew D. A. Orkney,
Stacy Y. Kim,
Martin P. Rey,
Eric P. Andersson,
Oscar Agertz,
Andrew Pontzen,
Roberta Vieliute,
Dhairya Kataria,
Kiah Jeneway
Abstract:
Ultra-faint dwarf galaxies (UFDs) are commonly found in close proximity to the Milky Way and other massive spiral galaxies. As such, their projected stellar ellipticity and extended light distributions are often thought to owe to tidal forces. In this paper, we study the projected stellar ellipticities and faint stellar outskirts of tidally isolated ultra-faints drawn from the 'Engineering Dwarfs…
▽ More
Ultra-faint dwarf galaxies (UFDs) are commonly found in close proximity to the Milky Way and other massive spiral galaxies. As such, their projected stellar ellipticity and extended light distributions are often thought to owe to tidal forces. In this paper, we study the projected stellar ellipticities and faint stellar outskirts of tidally isolated ultra-faints drawn from the 'Engineering Dwarfs at Galaxy Formation's Edge' (EDGE) cosmological simulation suite. Despite their tidal isolation, our simulated dwarfs exhibit a wide range of projected ellipticities ($0.03 < \varepsilon < 0.85$), with many possessing anisotropic extended stellar haloes that mimic tidal tails, but owe instead to late-time accretion of lower mass companions. Furthermore, we find a strong causal relationship between ellipticity and formation time of an UFD, which is robust to a wide variation in the feedback model. We show that the distribution of projected ellipticities in our suite of simulated EDGE dwarfs matches well with that of 21 Local Group dwarf galaxies. Given the ellipticity in EDGE arises from an ex-situ accretion origin, the agreement in shape indicates the ellipticities of some observed dwarfs may also originate from a similar non-tidal scenario. The orbital parameters of these observed dwarfs further support that they are not currently tidally disrupting. If the baryonic content in these galaxies is still tidally intact, then the same may be true for their dark matter content, making these galaxies in our Local Group pristine laboratories for testing dark matter and galaxy formation models.
△ Less
Submitted 11 July, 2023;
originally announced July 2023.
-
Mortality Forecasting using Variational Inference
Authors:
Patrik Andersson,
Mathias Lindholm
Abstract:
This paper considers the problem of forecasting mortality rates. A large number of models have already been proposed for this task, but they generally have the disadvantage of either estimating the model in a two-step process, possibly losing efficiency, or relying on methods that are cumbersome for the practitioner to use. We instead propose using variational inference and the probabilistic progr…
▽ More
This paper considers the problem of forecasting mortality rates. A large number of models have already been proposed for this task, but they generally have the disadvantage of either estimating the model in a two-step process, possibly losing efficiency, or relying on methods that are cumbersome for the practitioner to use. We instead propose using variational inference and the probabilistic programming library Pyro for estimating the model. This allows for flexibility in modelling assumptions while still being able to estimate the full model in one step.
The models are fitted on Swedish mortality data and we find that the in-sample fit is good and that the forecasting performance is better than other popular models.
Code is available at https://github.com/LPAndersson/VImortality.
△ Less
Submitted 25 May, 2023;
originally announced May 2023.
-
Variational inference for max-stable processes
Authors:
Patrik Andersson,
Alexander Engberg
Abstract:
Max-stable processes provide natural models for the modelling of spatial extreme values observed at a set of spatial sites. Full likelihood inference for max-stable data is, however, complicated by the form of the likelihood function as it contains a sum over all partitions of sites. As such, the number of terms to sum over grows rapidly with the number of sites and quickly becomes prohibitively b…
▽ More
Max-stable processes provide natural models for the modelling of spatial extreme values observed at a set of spatial sites. Full likelihood inference for max-stable data is, however, complicated by the form of the likelihood function as it contains a sum over all partitions of sites. As such, the number of terms to sum over grows rapidly with the number of sites and quickly becomes prohibitively burdensome to compute.
We propose a variational inference approach to full likelihood inference that circumvents the problematic sum. To achieve this, we first posit a parametric family of partition distributions from which partitions can be sampled. Second, we optimise the parameters of the family in conjunction with the max-stable model to find the partition distribution best supported by the data, and to estimate the max-stable model parameters.
In a simulation study we show that our method enables full likelihood inference in higher dimensions than previous methods, and is readily applicable to data sets with a large number of observations. Furthermore, our method can easily be extended to a Bayesian setting. Code is available at https://github.com/LPAndersson/MaxStableVI.jl.
△ Less
Submitted 14 December, 2022;
originally announced December 2022.
-
PRISM: A Non-Equilibrium, Multiphase Interstellar Medium Model for Radiation Hydrodynamics Simulations of Galaxies
Authors:
Harley Katz,
Shenghua Liu,
Taysun Kimm,
Martin P. Rey,
Eric P. Andersson,
Alex J. Cameron,
Francisco Rodriguez-Montero,
Oscar Agertz,
Julien Devriendt,
Adrianne Slyz
Abstract:
We introduce the PRISM interstellar medium (ISM) model for thermochemistry and its implementation in the RAMSES-RTZ code. The model includes a non-equilibrium primordial, metal, and molecular chemistry network for 115 species coupled to on-the-fly multifrequency radiation transport. PRISM accurately accounts for the dominant ISM cooling and heating processes in the low-density regime (i.e.…
▽ More
We introduce the PRISM interstellar medium (ISM) model for thermochemistry and its implementation in the RAMSES-RTZ code. The model includes a non-equilibrium primordial, metal, and molecular chemistry network for 115 species coupled to on-the-fly multifrequency radiation transport. PRISM accurately accounts for the dominant ISM cooling and heating processes in the low-density regime (i.e. $ρ<10^5\ {\rm cm^{-3}}$), including photoheating, photoelectric heating, H$_2$ heating/cooling, cosmic-ray heating, H/He cooling, metal-line cooling, CO cooling, and dust cooling (recombination and gas-grain collisions). We validate the model by comparing 1D equilibrium simulations across six dex in metallicity to existing 1D ISM models in the literature. We apply PRISM to high-resolution (4.5 pc) isolated dwarf galaxy simulations that include state-of-the-art models for star formation and stellar feedback to take an inventory of which cooling and heating processes dominate each different gas phase of a galaxy and to understand the importance of non-equilibrium effects. We show that most of the ISM gas is either close to thermal equilibrium or exhibits a slight cooling instability, while from a chemical perspective, the non-equilibrium electron fraction is often more than three times higher or lower than the equilibrium value, which impacts cooling, heating, and observable emission lines. Electron enhancements are attributed to recombination lags while deficits are shown to be due to rapid cosmic-ray heating. The PRISM model and its coupling to RAMSES-RTZ is applicable to a wide variety of astrophysical scenarios, from cosmological simulations to isolated giant molecular clouds, and is particularly useful for understanding how changes to ISM physics impact observable quantities such as metallic emission lines.
△ Less
Submitted 8 November, 2022;
originally announced November 2022.
-
INFERNO: Galactic winds in dwarf galaxies with star-by-star simulations including runaway stars
Authors:
Eric P. Andersson,
Oscar Agertz,
Florent Renaud,
Romain Teyssier
Abstract:
The formation and evolution of galaxies have proved sensitive to the inclusion of stellar feedback, which is therefore crucial to any successful galaxy model. We present INFERNO, a new model for hydrodynamic simulations of galaxies, which incorporates resolved stellar objects with star-by-star calculations of when and where the injection of enriched material, momentum, and energy takes place. INFE…
▽ More
The formation and evolution of galaxies have proved sensitive to the inclusion of stellar feedback, which is therefore crucial to any successful galaxy model. We present INFERNO, a new model for hydrodynamic simulations of galaxies, which incorporates resolved stellar objects with star-by-star calculations of when and where the injection of enriched material, momentum, and energy takes place. INFERNO treats early stellar kinematics to include phenomena such as walkaway and runaway stars. We employ this innovative model on simulations of a dwarf galaxy and demonstrate that our physically motivated stellar feedback model can drive vigorous galactic winds. This is quantified by mass and metal loading factors in the range of $10-100$, and an energy loading factor close to unity. Outflows are established close to the disc, are highly multi-phase, spanning almost $8$ orders of magnitude in temperature, and with a clear dichotomy between mass ejected in cold, slow-moving ($T\lesssim5\times10^4\,{\rm K}$, $v<100\,{\rm km\,s}^{-1}$) gas and energy ejected in hot, fast-moving ($T>10^6\,{\rm K}$, $v>100\,{\rm km\,s}^{-1}$) gas. In contrast to massive disc galaxies, we find a surprisingly weak impact of the early stellar kinematics, with runaway stars having little to no effect on our results, despite exploding in diffuse gas outside the dense star-forming gas, as well as outside the galactic disc entirely. We demonstrate that this weak impact in dwarf galaxies stems from a combination of strong feedback and a porous interstellar medium, which obscure any unique signatures that runaway stars provide.
△ Less
Submitted 3 March, 2023; v1 submitted 13 September, 2022;
originally announced September 2022.
-
Quantum Computing: Fundamentals, Trends and Perspectives for Chemical and Biochemical Engineers
Authors:
Amirhossein Nourbakhsh,
Mark Nicholas Jones,
Kaur Kristjuhan,
Deborah Carberry,
Jay Karon,
Christian Beenfeldt,
Kyarash Shahriari,
Martin P. Andersson,
Mojgan A. Jadidi,
Seyed Soheil Mansouri
Abstract:
We use the benefits and components of classical computers every day. However, there are many types of problems which, as they grow in size, their computational complexity grows larger than classical computers will ever be able to solve. Quantum computing (QC) is a computation model that uses quantum physical properties to solve such problems. QC is at the early stage of large-scale adoption in var…
▽ More
We use the benefits and components of classical computers every day. However, there are many types of problems which, as they grow in size, their computational complexity grows larger than classical computers will ever be able to solve. Quantum computing (QC) is a computation model that uses quantum physical properties to solve such problems. QC is at the early stage of large-scale adoption in various industry domains to take advantage of the algorithmic speed-ups it has to offer. It can be applied in a variety of areas, such as computer science, mathematics, chemical and biochemical engineering, and the financial industry. The main goal of this paper is to give an overview to chemical and biochemical researchers and engineers who may not be familiar with quantum computation. Thus, the paper begins by explaining the fundamental concepts of QC. The second contribution this publication tries to tackle is the fact that the chemical engineering literature still lacks a comprehensive review of the recent advances of QC. Therefore, this article reviews and summarizes the state of the art to gain insight into how quantum computation can benefit and optimize chemical engineering issues. A bibliography analysis covers the comprehensive literature in QC and analyzes quantum computing research in chemical engineering on various publication topics, using Clarivate analytics covering the years 1990 to 2020. After the bibliographic analysis, relevant applications of QC in chemical and biochemical engineering are highlighted and a conclusion offers an outlook of future directions within the field.
△ Less
Submitted 8 January, 2022;
originally announced January 2022.
-
EDGE: The sensitivity of ultra-faint dwarfs to a metallicity-dependent initial mass function
Authors:
Mateo Prgomet,
Martin P. Rey,
Eric P. Andersson,
Alvaro Segovia Otero,
Oscar Agertz,
Florent Renaud,
Andrew Pontzen,
Justin I. Read
Abstract:
Motivated by the observed bottom-light initial mass function (IMF) in faint dwarfs, we study how a metallicity-dependent IMF affects the feedback budget and observables of an ultra-faint dwarf galaxy. We model the evolution of a low-mass ($\approx 8 \, \times \, 10^{8} \, \rm M_{\odot}$) dark matter halo with cosmological, zoomed hydrodynamical simulations capable of resolving individual supernova…
▽ More
Motivated by the observed bottom-light initial mass function (IMF) in faint dwarfs, we study how a metallicity-dependent IMF affects the feedback budget and observables of an ultra-faint dwarf galaxy. We model the evolution of a low-mass ($\approx 8 \, \times \, 10^{8} \, \rm M_{\odot}$) dark matter halo with cosmological, zoomed hydrodynamical simulations capable of resolving individual supernovae explosions, which we complement with an empirically motivated subgrid prescription for systematic IMF variations. In this framework, at the low gas metallicities typical of faint dwarfs, the IMF of newborn stellar populations becomes top-heavy, increasing the efficiency of supernova and photoionization feedback in regulating star formation. This results in a 100-fold reduction of the final stellar mass of the dwarf compared to a canonical IMF, at fixed dynamical mass. The increase in the feedback budget is nonetheless met by increased metal production from more numerous massive stars, leading to nearly constant iron content at $z=0$. A metallicity-dependent IMF therefore provides a mechanism to produce low-mass ($\rm M_{\star}\sim 10^3 \rm M_{\odot}$), yet enriched ($\rm [Fe/H]\approx -2$) field dwarf galaxies, thus opening a self-consistent avenue to populate the plateau in $\rm [Fe/H]$ at the faintest end of the mass--metallicity relation.
△ Less
Submitted 2 May, 2022; v1 submitted 1 July, 2021;
originally announced July 2021.
-
Runaway stars masquerading as star formation in galactic outskirts
Authors:
Eric P. Andersson,
Florent Renaud,
Oscar Agertz
Abstract:
In the outskirts of nearby spiral galaxies, star formation is observed in extremely low gas surface densities. Star formation in these regions, where the interstellar medium is dominated by diffuse atomic hydrogen, is difficult to explain with classic star formation theories. In this letter, we introduce runaway stars as an explanation for this observation. Runaway stars, produced by collisional d…
▽ More
In the outskirts of nearby spiral galaxies, star formation is observed in extremely low gas surface densities. Star formation in these regions, where the interstellar medium is dominated by diffuse atomic hydrogen, is difficult to explain with classic star formation theories. In this letter, we introduce runaway stars as an explanation for this observation. Runaway stars, produced by collisional dynamics in young stellar clusters, can travel kiloparsecs during their main-sequence lifetime. Using galactic-scale hydrodynamic simulations including a treatment of individual stars, we demonstrate that this mechanism enables the ejection of young massive stars into environments where the gas is not dense enough to trigger star formation. This results in the appearance of star formation in regions where it ought to be impossible. We conclude that runaway stars are a contributing, if not dominant, factor to the observations of star formation in the outskirts of spiral galaxies.
△ Less
Submitted 8 September, 2021; v1 submitted 22 October, 2020;
originally announced October 2020.
-
VINTERGATAN III: how to reset the metallicity of the Milky Way
Authors:
Florent Renaud,
Oscar Agertz,
Eric P. Andersson,
Justin I. Read,
Nils Ryde,
Thomas Bensby,
Martin P. Rey,
Diane K. Feuillet
Abstract:
Using the cosmological zoom simulation VINTERGATAN, we present a new scenario for the onset of star formation at the metal-poor end of the low-[$α$/Fe] sequence in a Milky Way-like galaxy. In this scenario, the galaxy is fueled by two distinct gas flows. One is enriched by outflows from massive galaxies, but not the other. While the former feeds the inner galactic region, the latter fuels an outer…
▽ More
Using the cosmological zoom simulation VINTERGATAN, we present a new scenario for the onset of star formation at the metal-poor end of the low-[$α$/Fe] sequence in a Milky Way-like galaxy. In this scenario, the galaxy is fueled by two distinct gas flows. One is enriched by outflows from massive galaxies, but not the other. While the former feeds the inner galactic region, the latter fuels an outer gas disk, inclined with respect to the main galactic plane, and with a significantly poorer chemical content. The first passage of the last major merger galaxy triggers tidal compression in the outer disk, which increases the gas density and eventually leads to star formation, at a metallicity 0.75 dex lower than the inner galaxy. This forms the first stars of the low-[$α$/Fe] sequence. These in situ stars have halo-like kinematics, similarly to what is observed in the Milky Way, due to the inclination of the outer disk which eventually aligns with the inner one via gravitational torques. We show that this tilting disk scenario is likely to be common in Milky-Way like galaxies. This process implies that the low-[$α$/Fe] sequence is populated in situ, simultaneously from two formation channels, in the inner and the outer galaxy, with distinct metallicities. This contrasts with purely sequential scenarios for the assembly of the Milky Way disk and could be tested observationally.
△ Less
Submitted 29 March, 2021; v1 submitted 10 June, 2020;
originally announced June 2020.
-
VINTERGATAN II: the history of the Milky Way told by its mergers
Authors:
Florent Renaud,
Oscar Agertz,
Justin I. Read,
Nils Ryde,
Eric P. Andersson,
Thomas Bensby,
Martin P. Rey,
Diane K. Feuillet
Abstract:
Using the VINTERGATAN cosmological zoom simulation, we explore the contributions of the in situ and accreted material, and the effect of galaxy interactions and mergers in the assembly of a Milky Way-like galaxy. We find that the initial growth phase of galaxy evolution, dominated by repeated major mergers, provides the necessary physical conditions for the assembly of a thick, kinematically hot d…
▽ More
Using the VINTERGATAN cosmological zoom simulation, we explore the contributions of the in situ and accreted material, and the effect of galaxy interactions and mergers in the assembly of a Milky Way-like galaxy. We find that the initial growth phase of galaxy evolution, dominated by repeated major mergers, provides the necessary physical conditions for the assembly of a thick, kinematically hot disk populated by high-[$α$/Fe] stars, formed both in situ and in accreted satellite galaxies. We find that the diversity of evolutionary tracks followed by the simulated galaxy and its progenitors leads to very little overlap of the in situ and accreted populations for any given chemical composition. At a given age, the spread in [$α$/Fe] abundance ratio results from the diversity of physical conditions in VINTERGATAN and its satellites, with an enhancement in [$α$/Fe] found in stars formed during starburst episodes. Later, the cessation of the merger activity promotes the in situ formation of stars in the low-[$α$/Fe] regime, in a radially extended, thin and overall kinematically colder disk, thus establishing chemically bimodal thin and thick disks, in line with observations. We draw links between notable features in the [Fe/H] - [$α$/Fe] plane with their physical causes, and propose a comprehensive formation scenario explaining self-consistently, in the cosmological context, the main observed properties of the Milky Way.
△ Less
Submitted 29 March, 2021; v1 submitted 10 June, 2020;
originally announced June 2020.
-
VINTERGATAN I: The origins of chemically, kinematically and structurally distinct discs in a simulated Milky Way-mass galaxy
Authors:
Oscar Agertz,
Florent Renaud,
Sofia Feltzing,
Justin I. Read,
Nils Ryde,
Eric P. Andersson,
Martin P. Rey,
Thomas Bensby,
Diane K. Feuillet
Abstract:
Spectroscopic surveys of the Milky Way's stars have revealed spatial, chemical and kinematical structures that encode its history. In this work, we study their origins using a cosmological zoom simulation, VINTERGATAN, of a Milky Way-mass disc galaxy. We find that in connection to the last major merger at $z\sim 1.5$, cosmological accretion leads to the rapid formation of an outer, metal-poor, low…
▽ More
Spectroscopic surveys of the Milky Way's stars have revealed spatial, chemical and kinematical structures that encode its history. In this work, we study their origins using a cosmological zoom simulation, VINTERGATAN, of a Milky Way-mass disc galaxy. We find that in connection to the last major merger at $z\sim 1.5$, cosmological accretion leads to the rapid formation of an outer, metal-poor, low-[$α$/Fe] gas disc around the inner, metal-rich galaxy containing the old high-[$α$/Fe] stars. This event leads to a bimodality in [$α$/Fe] over a range of [Fe/H]. A detailed analysis of how the galaxy evolves since $z\sim 1$ is presented. We demonstrate the way in which inside-out growth shapes the radial surface density and metallicity profile and how radial migration preferentially relocates stars from the inner to the outer disc. Secular disc heating is found to give rise to increasing velocity dispersions and scaleheights with stellar age, which together with disc flaring explains several trends observed in the Milky Way, including shallower radial [Fe/H]-profiles above the midplane. We show how the galaxy formation scenario imprints non-trivial mappings between structural associations (i.e. thick and thin discs), velocity dispersions, $α$-enhancements, and ages of stars, e.g. the most metal-poor stars in the low-[$α$/Fe] sequence are found to have a scaleheight comparable to old high-[$α$/Fe] stars. Finally, we illustrate how at low spatial resolution, comparable to the thickness of the galaxy, the proposed pathway to distinct sequences in [$α$/Fe]-[Fe/H] cannot be captured.
△ Less
Submitted 12 June, 2020; v1 submitted 10 June, 2020;
originally announced June 2020.
-
How runaway stars boost galactic outflows
Authors:
Eric P. Andersson,
Oscar Agertz,
Florent Renaud
Abstract:
Roughly ten per cent of OB stars are kicked out of their natal clusters before ending their life as supernovae. These so-called runaway stars can travel hundreds of parsecs into the low-density interstellar medium, where momentum and energy from stellar feedback is efficiently deposited. In this work we explore how this mechanism affects large scale properties of the galaxy, such as outflows. To d…
▽ More
Roughly ten per cent of OB stars are kicked out of their natal clusters before ending their life as supernovae. These so-called runaway stars can travel hundreds of parsecs into the low-density interstellar medium, where momentum and energy from stellar feedback is efficiently deposited. In this work we explore how this mechanism affects large scale properties of the galaxy, such as outflows. To do so we use a new model which treats OB stars and their associated feedback processes on a star-by-star basis. With this model we compare two hydrodynamical simulations of Milky Way-like galaxies, one where we include runaways, and one where we ignore them. Including runaway stars leads to twice as many supernovae explosions in regions with gas densities ranging from 1e-5 cm^-3 to 1e-3 cm^-3. This results in more efficient heating of the inter-arm regions, and drives strong galactic winds with mass loading factors boosted by up to one order of magnitude. These outflows produce a more massive and extended multi-phase circumgalactic medium, as well as a population of dense clouds in the halo. Conversely, since less energy and momentum is released in the dense star forming regions, the cold phase of the interstellar medium is less disturbed by feedback effects.
△ Less
Submitted 1 May, 2020; v1 submitted 27 March, 2020;
originally announced March 2020.
-
Tidal stripping as a mechanism for placing globular clusters on wide orbits: the case of MGC1 in M31
Authors:
Eric P. Andersson,
Melvyn B. Davies
Abstract:
The globular clusters of large spiral galaxies can be divided into two populations: one which formed in-situ and one which comprises clusters tidally stripped away from other galaxies. In this paper we investigate the contribution to the outer globular cluster population in the M31 galaxy through donation of clusters from dwarf galaxies. We test this numerically by comparing the contribution of gl…
▽ More
The globular clusters of large spiral galaxies can be divided into two populations: one which formed in-situ and one which comprises clusters tidally stripped away from other galaxies. In this paper we investigate the contribution to the outer globular cluster population in the M31 galaxy through donation of clusters from dwarf galaxies. We test this numerically by comparing the contribution of globular clusters from simulated encounters to the observed M31 globular cluster population. To constrain our simulations, we specifically investigate the outermost globular cluster in the M31 system, MGC1. The remote location of MGC1 favours the idea of it being captured, however, the cluster is devoid of features associated with tidal interactions. Hence we separate simulations where tidal features are present and where they are hidden. We find that our simulated encounters can place clusters on MGC1-like orbits. In addition, we find that tidal stripping of clusters from dwarf galaxies leaves them on orbits having a range of separations, broadly matching those observed in M31. We find that the specific energies of globular clusters captured by M31 closely matches those of the incoming host dwarf galaxies. Furthermore, in our simulations we find an equal number of accreted clusters on co-rotating and counter-rotating orbits within M31 and use this to infer the fraction of clusters that has been accreted. We find that even close in roughly 50% of the clusters are accreted, whilst this figure increases to over 80% further out.
△ Less
Submitted 11 March, 2019; v1 submitted 24 August, 2018;
originally announced August 2018.
-
ToxTrac: a fast and robust software for tracking organisms
Authors:
Alvaro Rodriquez,
Hanqing Zhang,
Jonatan Klaminder,
Tomas Brodin,
Patrik L. Andersson,
Magnus Andersson
Abstract:
1. Behavioral analysis based on video recording is becoming increasingly popular within research fields such as; ecology, medicine, ecotoxicology, and toxicology. However, the programs available to analyze the data, which are; free of cost, user-friendly, versatile, robust, fast and provide reliable statistics for different organisms (invertebrates, vertebrates and mammals) are significantly limit…
▽ More
1. Behavioral analysis based on video recording is becoming increasingly popular within research fields such as; ecology, medicine, ecotoxicology, and toxicology. However, the programs available to analyze the data, which are; free of cost, user-friendly, versatile, robust, fast and provide reliable statistics for different organisms (invertebrates, vertebrates and mammals) are significantly limited.
2. We present an automated open-source executable software (ToxTrac) for image-based tracking that can simultaneously handle several organisms monitored in a laboratory environment. We compare the performance of ToxTrac with current accessible programs on the web.
3. The main advantages of ToxTrac are: i) no specific knowledge of the geometry of the tracked bodies is needed; ii) processing speed, ToxTrac can operate at a rate >25 frames per second in HD videos using modern desktop computers; iii) simultaneous tracking of multiple organisms in multiple arenas; iv) integrated distortion correction and camera calibration; v) robust against false positives; vi) preservation of individual identification if crossing occurs; vii) useful statistics and heat maps in real scale are exported in: image, text and excel formats.
4. ToxTrac can be used for high speed tracking of insects, fish, rodents or other species, and provides useful locomotor information. We suggest using ToxTrac for future studies of animal behavior independent of research area. Download ToxTrac here: https://toxtrac.sourceforge.io
△ Less
Submitted 8 June, 2017;
originally announced June 2017.
-
A robust particle detection algorithm based on symmetry
Authors:
Alvaro Rodriguez,
Hanqing Zhang,
Krister Wiklund,
Tomas Brodin,
Jonatan Klaminder,
Patrik Andersson,
Magnus Andersson
Abstract:
Particle tracking is common in many biophysical, ecological, and micro-fluidic applications. Reliable tracking information is heavily dependent on of the system under study and algorithms that correctly determines particle position between images. However, in a real environmental context with the presence of noise including particular or dissolved matter in water, and low and fluctuating light con…
▽ More
Particle tracking is common in many biophysical, ecological, and micro-fluidic applications. Reliable tracking information is heavily dependent on of the system under study and algorithms that correctly determines particle position between images. However, in a real environmental context with the presence of noise including particular or dissolved matter in water, and low and fluctuating light conditions, many algorithms fail to obtain reliable information. We propose a new algorithm, the Circular Symmetry algorithm (C-Sym), for detecting the position of a circular particle with high accuracy and precision in noisy conditions. The algorithm takes advantage of the spatial symmetry of the particle allowing for subpixel accuracy. We compare the proposed algorithm with four different methods using both synthetic and experimental datasets. The results show that C-Sym is the most accurate and precise algorithm when tracking micro-particles in all tested conditions and it has the potential for use in applications including tracking biota in their environment.
△ Less
Submitted 11 May, 2016;
originally announced May 2016.
-
Light-ion production in the interaction of 175 MeV quasi-mono-energetic neutrons with iron and with bismuth
Authors:
R. Bevilacqua,
K. Jansson,
S. Pomp,
P. Andersson,
J. Blomgren,
C. Gustavsson,
A. Hjalmarsson,
V. D. Simutkin,
M. Österlund,
A. J. Koning,
A. V. Prokofiev,
M. Hayashi,
S. Hirayama,
Y. Naitou,
Y. Watanabe,
U. Tippawan,
S. G. Mashnik,
L. M. Kerby,
F. -R. Lecolley,
N. Marie,
J. -C. David,
S. Leray
Abstract:
Nuclear data for neutron-induced reactions in the intermediate energy range of 20 to 200 MeV are of great importance for the development of nuclear reaction codes since little data exist in that range. Also several different applications benefit from such data, notably accelerator-driven incineration of nuclear waste. The Medley setup was used for a series of measurements of p, d, t, $^3$He and…
▽ More
Nuclear data for neutron-induced reactions in the intermediate energy range of 20 to 200 MeV are of great importance for the development of nuclear reaction codes since little data exist in that range. Also several different applications benefit from such data, notably accelerator-driven incineration of nuclear waste. The Medley setup was used for a series of measurements of p, d, t, $^3$He and $α$-particle production by 175 MeV quasi-mono-energetic neutrons on various target nuclei. The measurements were performed at the The Svedberg Laboratory in Uppsala, Sweden. Eight detector telescopes placed at angles between 20$^\circ$ and 160$^\circ$ were used. Medley uses the $ΔE$-$ΔE$-$E$ technique to discriminate among the particle types and is able to measure double-differential cross sections over a wide range of particle energies. This paper briefly describes the experimental setup, summarizes the data analysis and reports on recent changes in the previously reported preliminary data set on bismuth. Experimental data are compared with INCL4.5-Abla07, MCNP6 using CEM03.03, TALYS and PHITS model calculations as well as with nuclear data evaluations. The models agree fairly well overall but in some cases systematic differences are found.
△ Less
Submitted 12 November, 2014;
originally announced November 2014.
-
Ice Formation via Deposition Mode Nucleation on Bare and Alcohol-covered Graphite Surfaces
Authors:
Xiangrui Kong,
Patrik U. Andersson,
Erik S. Thomson,
Jan B. C. Pettersson
Abstract:
Deposition of water on aerosol particles contributes to ice cloud formation in the atmosphere with implications for the water cycle and climate on Earth. The heterogeneous ice nucleation process is influenced by physico-chemical properties of the substrate, but the mechanisms remain incompletely understood. Here, we report on ice formation on bare and alcohol-covered graphite at temperatures from…
▽ More
Deposition of water on aerosol particles contributes to ice cloud formation in the atmosphere with implications for the water cycle and climate on Earth. The heterogeneous ice nucleation process is influenced by physico-chemical properties of the substrate, but the mechanisms remain incompletely understood. Here, we report on ice formation on bare and alcohol-covered graphite at temperatures from 175 to 213 K, probed by elastic helium and light scattering. Water has a low wettability on bare and butanol-covered graphite resulting in the growth of rough ice surfaces. In contrast, pre-adsorbed methanol provides hydrophilic surface sites and results in the formation of smooth crystalline ice; an effect that is pronounced also for sub-monolayer methanol coverages. The alcohols primarily reside at the ice surface and at the ice-graphite interface with a minor fraction being incorporated into the growing ice structures. Methanol has no observable effect on gas/solid water vapor exchange whereas butanol acts as a transport barrier for water resulting in a reduction in ice evaporation rate at 185 K. Implications for the description of deposition mode freezing are discussed.
△ Less
Submitted 12 September, 2013;
originally announced September 2013.
-
Characterization of a Be(p,xn) neutron source for fission yields measurements
Authors:
A. Mattera,
P. Andersson,
A. Hjalmarsson,
M. Lantz,
S. Pomp,
V. Rakopoulos,
A. Solders,
J. Valldor-Blücher,
D. Gorelov,
H. Penttilä,
S. Rinta-Antila,
A. V. Prokofiev,
E. Passoth,
R. Bedogni,
A. Gentile,
D. Bortot,
A. Esposito,
M. V. Introini,
A. Pola
Abstract:
We report on measurements performed at The Svedberg Laboratory (TSL) to characterize a proton-neutron converter for independent fission yield studies at the IGISOL-JYFLTRAP facility (Jyväskylä, Finland). A 30 MeV proton beam impinged on a 5 mm water-cooled Beryllium target. Two independent experimental techniques have been used to measure the neutron spectrum: a Time of Flight (TOF) system used to…
▽ More
We report on measurements performed at The Svedberg Laboratory (TSL) to characterize a proton-neutron converter for independent fission yield studies at the IGISOL-JYFLTRAP facility (Jyväskylä, Finland). A 30 MeV proton beam impinged on a 5 mm water-cooled Beryllium target. Two independent experimental techniques have been used to measure the neutron spectrum: a Time of Flight (TOF) system used to estimate the high-energy contribution, and a Bonner Sphere Spectrometer able to provide precise results from thermal energies up to 20 MeV. An overlap between the energy regions covered by the two systems will permit a cross-check of the results from the different techniques. In this paper, the measurement and analysis techniques will be presented together with some preliminary results.
△ Less
Submitted 2 April, 2013;
originally announced April 2013.
-
A microfluidic device for the study of the orientational dynamics of microrods
Authors:
Y. N. Mishra,
J. Einarsson,
O. A. John,
P. Andersson,
B. Mehlig,
D. Hanstorp
Abstract:
We describe a microfluidic device for studying the orientational dynamics of microrods. The device enables us to experimentally investigate the tumbling of microrods immersed in the shear flow in a microfluidic channel with a depth of 400 mu and a width of 2.5 mm. The orientational dynamics was recorded using a 20 X microscopic objective and a CCD camera. The microrods were produced by shearing mi…
▽ More
We describe a microfluidic device for studying the orientational dynamics of microrods. The device enables us to experimentally investigate the tumbling of microrods immersed in the shear flow in a microfluidic channel with a depth of 400 mu and a width of 2.5 mm. The orientational dynamics was recorded using a 20 X microscopic objective and a CCD camera. The microrods were produced by shearing microdroplets of photocurable epoxy resin. We show different examples of empirically observed tumbling. On the one hand we find that short stretches of the experimentally determined time series are well described by fits to solutions of Jeffery's approximate equation of motion [Jeffery, Proc. R. Soc. London. 102 (1922), 161-179]. On the other hand we find that the empirically observed trajectories drift between different solutions of Jeffery's equation. We discuss possible causes of this orbit drift.
△ Less
Submitted 11 February, 2012;
originally announced February 2012.
-
Collision Dynamics and Solvation of Water Molecules in a Liquid Methanol Film
Authors:
Erik S. Thomson,
Xiangrui Kong,
Patrik U. Andersson,
Nikola Marković,
Jan B. C. Pettersson
Abstract:
Environmental molecular beam experiments are used to examine water interactions with liquid methanol films at temperatures from 170 K to 190 K. We find that water molecules with 0.32 eV incident kinetic energy are efficiently trapped by the liquid methanol. The scattering process is characterized by an efficient loss of energy to surface modes with a minor component of the incident beam that is in…
▽ More
Environmental molecular beam experiments are used to examine water interactions with liquid methanol films at temperatures from 170 K to 190 K. We find that water molecules with 0.32 eV incident kinetic energy are efficiently trapped by the liquid methanol. The scattering process is characterized by an efficient loss of energy to surface modes with a minor component of the incident beam that is inelastically scattered. Thermal desorption of water molecules has a well characterized Arrhenius form with an activation energy of 0.47{\pm}0.11 eV and pre-exponential factor of 4.6 {\times} 10^(15{\pm}3) s^(-1). We also observe a temperature dependent incorporation of incident water into the methanol layer. The implication for fundamental studies and environmental applications is that even an alcohol as simple as methanol can exhibit complex and temperature dependent surfactant behavior.
△ Less
Submitted 8 September, 2011;
originally announced September 2011.
-
Rydberg Matter clusters of alkali metal atoms: the link between meteoritic matter, polar mesosphere summer echoes (PMSE), sporadic sodium layers, polar mesospheric clouds (PMCs, NLCs), and ion chemistry
Authors:
Frans Olofson,
Patrik U. Andersson,
Leif Holmlid
Abstract:
A material exists which links together the influx of meteoritic matter from interplanetary space, the polar mesosphere summer echoes (PMSE), the sporadic sodium layers, the polar mesospheric clouds (PMCs, NLCs), and the observed ion chemistry in the mesosphere. The evidence in these research fields is here analyzed and found to agree well with the properties of Rydberg Matter (RM). This material…
▽ More
A material exists which links together the influx of meteoritic matter from interplanetary space, the polar mesosphere summer echoes (PMSE), the sporadic sodium layers, the polar mesospheric clouds (PMCs, NLCs), and the observed ion chemistry in the mesosphere. The evidence in these research fields is here analyzed and found to agree well with the properties of Rydberg Matter (RM). This material has been studied with numerous methods in the laboratory. Alkali atoms, mainly Na, reach the mesosphere in the form of interplanetary (meteoritic, cometary) dust. The planar RM clusters NaN usually contain N = 19, 37 or 61 atoms, and have the density of air at 90 km altitude where they float. The diameters of the clusters are 10-100 nm from laboratory high precision radio frequency spectroscopic studies. Such experiments show that RM clusters interact strongly with radar frequencies: this explains the radio frequency heating and reflection studies of PMSE layers. The clusters give the low temperature in the mesosphere by efficient selective radiation at long wavelengths, which is observed in RF emission experiments. The lowest possible stable temperature of the mesopause is calculated for the first time to be 121 K in agreement with measurements, based on the strong optical activity at long wavelengths in RM. Sporadic sodium layers are explained in a unique way as due to shockwaves in the RM layers. Due to the high electronic excitation energy in RM clusters, they induce efficient reactions forming ions of all atoms and molecules in the atmosphere thus providing condensation nuclei for water vapour. This finally gives the visible part of the PMC structure. The present contribution fills the gap between and partially replaces the separate theories used to describe the various aspects of these intriguing phenomena.
△ Less
Submitted 8 February, 2010;
originally announced February 2010.
-
Dynamical stabilization of the bcc phase in lanthanum and thorium by phonon-phonon interaction
Authors:
P. Souvatzis,
T. Bjorkman,
O. Eriksson,
P Andersson,
M. I. Katsnelson S- P Rudin
Abstract:
The recently developed self consistent {\it ab initio} lattice dynamical method (SCAILD) has been applied to the high temperature bcc phase of La and Th which are dynamically unstable at low temperatures. The bcc phase of these metals is found to be stabilized by phonon-phonon interactions. The calculated high temperature phonon frequencies for La are found to be in good agreement with the corre…
▽ More
The recently developed self consistent {\it ab initio} lattice dynamical method (SCAILD) has been applied to the high temperature bcc phase of La and Th which are dynamically unstable at low temperatures. The bcc phase of these metals is found to be stabilized by phonon-phonon interactions. The calculated high temperature phonon frequencies for La are found to be in good agreement with the corresponding experimental data.
△ Less
Submitted 25 March, 2009;
originally announced March 2009.
-
On The Nature of the 5f States in Delta Plutonium
Authors:
J. M. Wills,
O. Eriksson,
A. Delin,
P. H. Andersson,
J. J. Joyce,
T. Durakiewicz,
M. T. Butterfield,
L. A. Morales,
D. P. Moore,
A. J. Arko
Abstract:
We present a theoretical model of the electronic structure of delta-Pu that is consistent with many of the electronic structure related properties of this complex metal. In particular we show that the theory is capable of reproducing the valence band photoelectron spectrum of delta-Pu. We report new experimental photoelectron spectra at several photon energies and present evidence that the elect…
▽ More
We present a theoretical model of the electronic structure of delta-Pu that is consistent with many of the electronic structure related properties of this complex metal. In particular we show that the theory is capable of reproducing the valence band photoelectron spectrum of delta-Pu. We report new experimental photoelectron spectra at several photon energies and present evidence that the electronic structure of delta-Pu is unique among the elements, involving a 5f shell with four 5f electrons in a localized multiplet, hybridizing with valence states, and approximately one 5f electron forming a completely delocalized band state.
△ Less
Submitted 31 July, 2003;
originally announced July 2003.
-
Spin-orbit coupling in the actinide elements: a critical evaluation of theoretical equilibrium volumes
Authors:
J. M. Wills,
P. H. Andersson,
L. Nordstrom,
P. Soderlind,
O. Eriksson
Abstract:
In a recent paper by Jones et al., it is argued, based on FP-LAPW band-structure calculations, that previous calculations of ground-state properties for actinides, using the FP-LMTO method implemented by J. M. Wills, are in error. We demonstrate in this paper that the conclusions of Jones et. al. are unfounded. Calculations using the FP-LMTO method are compared with calculations performed with t…
▽ More
In a recent paper by Jones et al., it is argued, based on FP-LAPW band-structure calculations, that previous calculations of ground-state properties for actinides, using the FP-LMTO method implemented by J. M. Wills, are in error. We demonstrate in this paper that the conclusions of Jones et. al. are unfounded. Calculations using the FP-LMTO method are compared with calculations performed with the FP-LAPW method and, in contradiction to statements by Jones et. al., they agree very well. The limitation of the present implementation of the spin-orbit coupling is discussed, where especially the 6p states pose a difficulty. The observed discrepancy (0-10%) between the FP-LMTO calculations and the recent FP-LAPW calculation of Jones et. al. is shown to be due to the choice of muffin-tin radius in the calculations. We argue that the choice of a constant muffin-tin radius, common for all calculated volumes, is less satisfactory compared to the choice of a muffin-tin radius that scales with the volume, especially when a large volume interval and open packed structures are considered. The conclusion of Jones et al. that alpha-Pu has partially delocalized 5f electrons is argued to be erroneous.
△ Less
Submitted 24 August, 1999;
originally announced August 1999.
-
Search for charginos, neutralinos, and gravitinos at LEP
Authors:
P. Andersson,
A. De Min,
J. J. Hernandez,
K. Hultqvist,
A. Lipniacka,
C. Matteuzzi,
S. Navas,
S. Ortuno,
P. Rebecchi,
A. Trombini
Abstract:
The hep-ex data base was decided not to be an appropriate place to make DELPHI notes public. Sorry for the inconvenience.
The hep-ex data base was decided not to be an appropriate place to make DELPHI notes public. Sorry for the inconvenience.
△ Less
Submitted 9 May, 1997; v1 submitted 30 April, 1997;
originally announced April 1997.