-
Nuclear uncertainties associated with the ejecta of a neutron-star black-hole accretion disk
Authors:
M. R. Mumpower,
T. M. Sprouse,
J. M. Miller,
K. A. Lund,
J. Cabrera Garcia,
N. Vassh,
G. C. McLaughlin,
R. Surman
Abstract:
The simulation of heavy element nucleosynthesis requires input from yet-to-be-measured nuclear properties. The uncertainty in the values of these off-stability nuclear properties propagates to uncertainties in the predictions of elemental and isotopic abundances. However, for any given astrophysical explosion, there are many different trajectories, i.e. temperature and density histories, experienc…
▽ More
The simulation of heavy element nucleosynthesis requires input from yet-to-be-measured nuclear properties. The uncertainty in the values of these off-stability nuclear properties propagates to uncertainties in the predictions of elemental and isotopic abundances. However, for any given astrophysical explosion, there are many different trajectories, i.e. temperature and density histories, experienced by outflowing material and thus different nuclear properties can come into play. We consider combined nucleosynthesis results from 460,000 trajectories from a neutron star-black hole accretion disk and the find spread in elemental predictions due solely to unknown nuclear properties to be a factor of a few. We analyze this relative spread in model predictions due to nuclear variations and conclude that the uncertainties can be attributed to a combination of properties in a given region of the abundance pattern. We calculate a cross-correlation between mass changes and abundance changes to show how variations among the properties of participating nuclei may be explored. Our results provide further impetus for measurements of multiple quantities on individual short-lived neutron-rich isotopes at modern experimental facilities.
△ Less
Submitted 3 April, 2024;
originally announced April 2024.
-
Magnetic Field Strength Effects on Nucleosynthesis from Neutron Star Merger Outflows
Authors:
Kelsey Alexandra Lund,
Gail McLaughlin,
Jonah Miller,
Matthew Mumpower
Abstract:
Magnetohydrodynamic turbulence drives the central engine of post-merger remnants, potentially powering both a nucleosynthetically active disk wind and the relativistic jet behind a short gamma ray burst. We explore the impact of the magnetic field on this engine by simulating three post-merger black hole accretion disks using general relativistic magnetohydrodynamics with Monte Carlo neutrino tran…
▽ More
Magnetohydrodynamic turbulence drives the central engine of post-merger remnants, potentially powering both a nucleosynthetically active disk wind and the relativistic jet behind a short gamma ray burst. We explore the impact of the magnetic field on this engine by simulating three post-merger black hole accretion disks using general relativistic magnetohydrodynamics with Monte Carlo neutrino transport, in each case varying the initial magnetic field strength. We find increasing ejecta masses associated with increasing magnetic field strength. We find that a fairly robust main r -process pattern is produced in all three cases, scaled by the ejected mass. Changing the initial magnetic field strength has a considerable effect on the geometry of the outflow and hints at complex central engine dynamics influencing lanthanide outflows. We find that actinide production is especially sensitive to magnetic field strength, with overall actinide mass fraction calculated at 1 Gyr post-merger increasing by more than a factor of six with a tenfold increase in magnetic field strength. This hints at a possible connection to the variability in actinide enhancements exhibited by metal poor, r -process-enhanced stars.
△ Less
Submitted 9 November, 2023;
originally announced November 2023.
-
Emergent nucleosynthesis from a 1.2 second long simulation of a black-hole accretion disk
Authors:
T. M. Sprouse,
K. A. Lund,
J. M. Miller,
G. C. McLaughlin,
M. R. Mumpower
Abstract:
We simulate a black-hole accretion disk system with full-transport general relativistic neutrino radiation magnetohydrodynamics (GR$ν$RMHD) for 1.2 seconds. This system is likely to form after the merger of two compact objects and is thought to be a robust site of $r$-process nucleosynthesis. We consider the case of a black-hole accretion disk arising from the merger of two neutron stars. Our simu…
▽ More
We simulate a black-hole accretion disk system with full-transport general relativistic neutrino radiation magnetohydrodynamics (GR$ν$RMHD) for 1.2 seconds. This system is likely to form after the merger of two compact objects and is thought to be a robust site of $r$-process nucleosynthesis. We consider the case of a black-hole accretion disk arising from the merger of two neutron stars. Our simulation time coincides with the nucleosynthesis timescale of the $r$ process ($\sim$ 1 second). Because these simulations are time consuming, it is common practice to run for `short' duration of approximately 0.1 to 0.3 seconds. We analyze the nucleosynthetic outflow from this system and compare the results between stopping at 0.12 and 1.2 seconds respectively. We find that the addition of mass ejected in the longer simulation as well as more favorable thermodynamic conditions from emergent viscous ejecta greatly impacts the nucleosynthetic outcome. We quantify the error in nucleosynthetic outcomes between short and long cuts.
△ Less
Submitted 14 September, 2023;
originally announced September 2023.
-
Superheavy Elements in Kilonovae
Authors:
Erika M. Holmbeck,
Jennifer Barnes,
Kelsey A. Lund,
Trevor M. Sprouse,
G. C. McLaughlin,
Matthew R. Mumpower
Abstract:
As LIGO-Virgo-KAGRA enters its fourth observing run, a new opportunity to search for electromagnetic counterparts of compact object mergers will also begin. The light curves and spectra from the first "kilonova" associated with a binary neutron star binary (NSM) suggests that these sites are hosts of the rapid neutron capture ("$r$") process. However, it is unknown just how robust elemental produc…
▽ More
As LIGO-Virgo-KAGRA enters its fourth observing run, a new opportunity to search for electromagnetic counterparts of compact object mergers will also begin. The light curves and spectra from the first "kilonova" associated with a binary neutron star binary (NSM) suggests that these sites are hosts of the rapid neutron capture ("$r$") process. However, it is unknown just how robust elemental production can be in mergers. Identifying signposts of the production of particular nuclei is critical for fully understanding merger-driven heavy-element synthesis. In this study, we investigate the properties of very neutron rich nuclei for which superheavy elements ($Z\geq 104$) can be produced in NSMs and whether they can similarly imprint a unique signature on kilonova light-curve evolution. A superheavy-element signature in kilonovae represents a route to establishing a lower limit on heavy-element production in NSMs as well as possibly being the first evidence of superheavy element synthesis in nature. Favorable NSMs conditions yield a mass fraction of superheavy elements is $X_{Z\geq 104}\approx 3\times 10^{-2}$ at 7.5 hours post-merger. With this mass fraction of superheavy elements, we find that the component of kilonova light curves possibly containing superheavy elements may appear similar to those arising from lanthanide-poor ejecta. Therefore, photometric characterizations of superheavy-element rich kilonova may possibly misidentify them as lanthanide-poor events.
△ Less
Submitted 1 June, 2023; v1 submitted 4 April, 2023;
originally announced April 2023.
-
The Influence of Beta Decay Rates on r-Process Observables
Authors:
Kelsey A. Lund,
J. Engel,
G. C. McLaughlin,
M. R. Mumpower,
E. M. Ney,
R. Surman
Abstract:
The rapid neutron capture process (r-process) is one of the main mechanisms whereby elements heavier than iron are synthesized, and is entirely responsible for the natural production of the actinides. Kilonova emissions are modeled as being largely powered by the radioactive decay of species synthesized via the r -process. Given that the r -process occurs far from nuclear stability, unmeasured bet…
▽ More
The rapid neutron capture process (r-process) is one of the main mechanisms whereby elements heavier than iron are synthesized, and is entirely responsible for the natural production of the actinides. Kilonova emissions are modeled as being largely powered by the radioactive decay of species synthesized via the r -process. Given that the r -process occurs far from nuclear stability, unmeasured beta decay rates play an essential role in setting the time scale for the r -process. In an effort to better understand the sensitivity of kilonova modeling to different theoretical global beta-decay descriptions, we incorporate these into nucleosynthesis calculations. We compare the results of these calculations and highlight differences in kilonova nuclear energy generation and light curve predictions, as well as final abundances and their implications for nuclear cosmochronometry. We investigate scenarios where differences in beta decay rates are responsible for increased nuclear heating on time scales of days that propagates into a significantly increased average bolometric luminosity between 1-10 days post-merger. We identify key nuclei, both measured and unmeasured, whose decay rates are directly impact nuclear heating generation on timescales responsible for light curve evolution. We also find that uncertainties in beta decay rates significantly impact ages estimates from cosmochronometry.
△ Less
Submitted 12 August, 2022;
originally announced August 2022.
-
Field linkage and magnetic helicity density
Authors:
K. Lund,
M. Jardine,
A. J. B. Russell,
J. -F. Donati,
R. Fares,
C. P. Folsom
Abstract:
The helicity of a magnetic field is a fundamental property that is conserved in ideal MHD. It can be explored in the stellar context by mapping large-scale magnetic fields across stellar surfaces using Zeeman-Doppler imaging. A recent study of 51 stars in the mass range 0.1-1.34 M$_\odot$ showed that the photospheric magnetic helicity density follows a single power law when plotted against the tor…
▽ More
The helicity of a magnetic field is a fundamental property that is conserved in ideal MHD. It can be explored in the stellar context by mapping large-scale magnetic fields across stellar surfaces using Zeeman-Doppler imaging. A recent study of 51 stars in the mass range 0.1-1.34 M$_\odot$ showed that the photospheric magnetic helicity density follows a single power law when plotted against the toroidal field energy, but splits into two branches when plotted against the poloidal field energy. These two branches divide stars above and below $\sim$ 0.5 M$_\odot$. We present here a novel method of visualising the helicity density in terms of the linkage of the toroidal and poloidal fields that are mapped across the stellar surface. This approach allows us to classify the field linkages that provide the helicity density for stars of different masses and rotation rates. We find that stars on the lower-mass branch tend to have toroidal fields that are non-axisymmetric and so link through regions of positive and negative poloidal field. A lower-mass star may have the same helicity density as a higher-mass star, despite having a stronger poloidal field. Lower-mass stars are therefore less efficient at generating large-scale helicity.
△ Less
Submitted 11 March, 2021; v1 submitted 22 February, 2021;
originally announced February 2021.
-
Kilonovae across the nuclear physics landscape: The impact of nuclear physics uncertainties on r-process-powered emission
Authors:
Jennifer Barnes,
Y. L. Zhu,
K. A. Lund,
T. M. Sprouse,
N. Vassh,
G. C. McLaughlin,
M. R. Mumpower,
R. Surman
Abstract:
Merging neutron stars produce "kilonovae"---electromagnetic transients powered by the decay of unstable nuclei synthesized via rapid neutron capture (the r-process) in material that is gravitationally unbound during inspiral and coalescence. Kilonova emission, if accurately interpreted, can be used to characterize the masses and compositions of merger-driven outflows, helping to resolve a long-sta…
▽ More
Merging neutron stars produce "kilonovae"---electromagnetic transients powered by the decay of unstable nuclei synthesized via rapid neutron capture (the r-process) in material that is gravitationally unbound during inspiral and coalescence. Kilonova emission, if accurately interpreted, can be used to characterize the masses and compositions of merger-driven outflows, helping to resolve a long-standing debate about the origins of r-process material in the Universe. We explore how the uncertain properties of nuclei involved in the r-process complicate the inference of outflow properties from kilonova observations. Using r-process simulations, we show how nuclear physics uncertainties impact predictions of radioactive heating and element synthesis. For a set of models that span a large range in both predicted heating and final abundances, we carry out detailed numerical calculations of decay product thermalization and radiation transport in a kilonova ejecta with a fixed mass and density profile. The light curves associated with our models exhibit great diversity in their luminosities, with peak brightness varying by more than an order of magnitude. We also find variability in the shape of the kilonova light curves and their color, which in some cases runs counter to the expectation that increasing levels of lanthanide and/or actinide enrichment will be correlated with longer, dimmer, redder emission.
△ Less
Submitted 21 October, 2020;
originally announced October 2020.
-
Modeling Kilonova Light Curves: Dependence on Nuclear Inputs
Authors:
Y. L. Zhu,
K. Lund,
J. Barnes,
T. M. Sprouse,
N. Vassh,
G. C. McLaughlin,
M. R. Mumpower,
R. Surman
Abstract:
The mergers of binary neutron stars, as well as black hole-neutron star systems, are expected to produce an electromagnetic counterpart that can be analyzed to infer the element synthesis that occurred in these events. We investigate one source of uncertainties pertinent to lanthanide-rich outflows: the nuclear inputs to rapid neutron capture nucleosynthesis calculations. We begin by examining thi…
▽ More
The mergers of binary neutron stars, as well as black hole-neutron star systems, are expected to produce an electromagnetic counterpart that can be analyzed to infer the element synthesis that occurred in these events. We investigate one source of uncertainties pertinent to lanthanide-rich outflows: the nuclear inputs to rapid neutron capture nucleosynthesis calculations. We begin by examining thirty-two different combinations of nuclear inputs: eight mass models, two types of spontaneous fission rates, and two types of fission daughter product distributions. We find that such nuclear physics uncertainties typically generate at least one order of magnitude uncertainty in key quantities such as the nuclear heating (one and a half orders of magnitude at one day post-merger), the bolometric luminosity (one order of magnitude at five days post-merger), and the inferred mass of material from the bolometric luminosity (factor of eight when considering the eight to ten days region). Since particular nuclear processes are critical for determining the electromagnetic signal, we provide tables of key nuclei undergoing $β$-decay, $α$-decay, and spontaneous fission important for heating at different times, identifying decays that are common among the many nuclear input combinations.
△ Less
Submitted 7 October, 2020;
originally announced October 2020.
-
Measuring stellar magnetic helicity density
Authors:
K. Lund,
M. Jardine,
L. T. Lehmann,
D. H. Mackay,
V. See,
A. A. Vidotto,
J. -F. Donati,
R. Fares,
C. P. Folsom,
S. V. Jeffers,
S. C. Marsden,
J. Morin,
P. Petit
Abstract:
Helicity is a fundamental property of a magnetic field but to date it has only been possible to observe its evolution in one star - the Sun. In this paper we provide a simple technique for mapping the large-scale helicity density across the surface of any star using only observable quantities: the poloidal and toroidal magnetic field components (which can be determined from Zeeman-Doppler imaging)…
▽ More
Helicity is a fundamental property of a magnetic field but to date it has only been possible to observe its evolution in one star - the Sun. In this paper we provide a simple technique for mapping the large-scale helicity density across the surface of any star using only observable quantities: the poloidal and toroidal magnetic field components (which can be determined from Zeeman-Doppler imaging) and the stellar radius. We use a sample of 51 stars across a mass range of 0.1-1.34 M$_\odot$ to show how the helicity density relates to stellar mass, Rossby number, magnetic energy and age. We find that the large-scale helicity density increases with decreasing Rossby number $R_o$, peaking at $R_o \simeq 0.1$, with a saturation or decrease below that. For both fully- and partially-convective stars we find that the mean absolute helicity density scales with the mean squared toroidal magnetic flux density according to the power law: $|\langle{h\,}\rangle|$ $\propto$ $\langle{\rm{B_{tor}}^2_{}\,\rangle}^{0.86\,\pm\,0.04}$. The scatter in this relation is consistent with the variation across a solar cycle, which we compute using simulations and observations across solar cycles 23 and 24 respectively. We find a significant decrease in helicity density with age.
△ Less
Submitted 31 January, 2020;
originally announced January 2020.
-
Radiation hydrodynamic simulations of massive star formation via gravitationally trapped HII regions - Spherically symmetric ionised accretion flows
Authors:
Kristin Lund,
Kenneth Wood,
Diego Falceta-Gonçalves,
Bert Vandenbroucke,
Nina Sartorio,
Ian Bonnell,
Katharine Johnston,
Eric Keto
Abstract:
This paper investigates the gravitational trapping of HII regions predicted by steady-state analysis using radiation hydrodynamical simulations. We present idealised spherically symmetric radiation hydrodynamical simulations of the early evolution of HII regions including the gravity of the central source. As with analytic steady state solutions of spherically symmetric ionised Bondi accretion flo…
▽ More
This paper investigates the gravitational trapping of HII regions predicted by steady-state analysis using radiation hydrodynamical simulations. We present idealised spherically symmetric radiation hydrodynamical simulations of the early evolution of HII regions including the gravity of the central source. As with analytic steady state solutions of spherically symmetric ionised Bondi accretion flows, we find gravitationally trapped HII regions with accretion through the ionisation front onto the source. We found that, for a constant ionising luminosity, fluctuations in the ionisation front are unstable. This instability only occurs in this spherically symmetric accretion geometry. In the context of massive star formation, the ionising luminosity increases with time as the source accretes mass. The maximum radius of the recurring HII region increases on the accretion timescale until it reaches the sonic radius, where the infall velocity equals the sound speed of the ionised gas, after which it enters a pressure-driven expansion phase. This expansion prevents accretion of gas through the ionisation front, the accretion rate onto the star decreases to zero, and it stops growing from accretion. Because of the time required for any significant change in stellar mass and luminosity through accretion our simulations keep both mass and luminosity constant and follow the evolution from trapped to expanding in a piecewise manner. Implications of this evolution of HII regions include a continuation of accretion of material onto forming stars for a period after the star starts to emit ionising radiation, and an extension of the lifetime of ultracompact HII regions.
△ Less
Submitted 1 March, 2019;
originally announced March 2019.
-
Testing the stability of supersonic ionized Bondi accretion flows with radiation hydrodynamics
Authors:
Bert Vandenbroucke,
Nina S. Sartorio,
Kenneth Wood,
Kristin Lund,
Diego Falceta-Gonçalves,
Thomas J. Haworth,
Ian Bonnell,
Eric Keto,
Daniel Tootill
Abstract:
We investigate the general stability of 1D spherically symmetric ionized Bondi accretion onto a massive object in the specific context of accretion onto a young stellar object. We first derive a new analytic expression for a steady state two temperature solution that predicts the existence of compact and hypercompact HII regions. We then show that this solution is only marginally stable if ionizat…
▽ More
We investigate the general stability of 1D spherically symmetric ionized Bondi accretion onto a massive object in the specific context of accretion onto a young stellar object. We first derive a new analytic expression for a steady state two temperature solution that predicts the existence of compact and hypercompact HII regions. We then show that this solution is only marginally stable if ionization is treated self-consistently. This leads to a recurring collapse of the HII region over time. We derive a semi-analytic model to explain this instability, and test it using spatially converged 1D radiation hydrodynamical simulations. We discuss the implications of the 1D instability on 3D radiation hydrodynamics simulations of supersonic accreting flows.
△ Less
Submitted 1 March, 2019;
originally announced March 2019.
-
The formation of high-mass binary star systems
Authors:
Kristin Lund,
Ian Bonnell
Abstract:
We develop a semi-analytic model to investigate how accretion onto wide low-mass binary stars can result in a close high-mass binary system. The key ingredient is to allow mass accretion while limiting the gain in angular momentum. We envision this process as being regulated by an external magnetic field during infall. Molecular clouds are made to collapse spherically with material either accretin…
▽ More
We develop a semi-analytic model to investigate how accretion onto wide low-mass binary stars can result in a close high-mass binary system. The key ingredient is to allow mass accretion while limiting the gain in angular momentum. We envision this process as being regulated by an external magnetic field during infall. Molecular clouds are made to collapse spherically with material either accreting onto the stars or settling in a disk. Our aim is to determine what initial conditions are needed for the resulting binary to be both massive and close. Whether material accretes, and what happens to the binary separation as a result, depends on the relative size of its specific angular momentum, compared to the specific angular momentum of the binary. When we add a magnetic field we are introducing a torque to the system which is capable of stripping the molecular cloud of some of its angular momentum, and consequently easing the formation of high-mass binaries. Our results suggest that clouds in excess of 1000 M$_\odot$ and radii of 0.5 pc or larger, can easily form binary systems with masses in excess of 25 M$_\odot$ and separations of order 10 R$_\odot$ with magnetic fields of order 100 μG (mass-to-flux ratios of order 5).
△ Less
Submitted 19 June, 2018;
originally announced June 2018.