Researcher profile

Enrico Ramirez-Ruiz

Enrico Ramirez-Ruiz contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
26works
0followers
5topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

26 published item(s)

preprint2026arXiv

Cascading from the winds to the disk: the universality of supernovae-driven turbulence in different galactic interstellar media

Star-forming galaxies are in a state of turbulence, with one of the principle components of the turbulence sourced by the constant injection of momentum from supernovae (SNe) explosions. Utilizing high-resolution stratified, gravito-hydrodynamical models of SNe-driven turbulence with interstellar medium (ISM) cooling and heating, we explore how SNe-driven turbulence changes across different galactic conditions, parameterized by the galactic mass and potential, SNe-driving rate, and seeding functions. We show that even though the underlying ISM changes between starburst and Milky Way analogue models, the velocity fluctuations in the turbulence of both models, but not the kinetic energy fluctuations, can be normalized into a universal, single cascade, $du^2(k)/dk \propto k^{-3/2}$, where $u$ is the velocity and $k$ is the wavemode, indicating that the structure of the turbulence is robust to significant changes in the ISM and SNe seeding. Moreover, the cascades connect smoothly from the winds into the galactic disk, pushing the outer-scale of the turbulence, $\ell_{\rm cor}$, to over $\ell_{\rm cor} \approx 6 \ell_0$, where $\ell_0$ is the gaseous scale-height. By providing an analytical model for the sound speed spectrum, $dc_s^2(k)/dk$, in the weak-cooling, adiabatic limit, we show that it is the compressible turbulent modes, $u_c$, that control the volume-filling phase structure of the galactic disks in our models, with $dc_s^2(k)/dk \propto k^{-2} \propto du_c^2(k)/dk$. This may indicate that galactic turbulence does not only have highly-universal features across different galaxies, but also directly sets the volume-filling hot and warm phase structure of the underlying galactic ISM through turbulent compressible modes.

preprint2022arXiv

An asymmetric electron-scattering photosphere around optical tidal disruption events

A star crossing the tidal radius of a supermassive black hole will be spectacularly ripped apart with an accompanying burst of radiation. A few tens of such tidal disruption events (TDEs) have now been identified in the optical wavelengths, but the exact origin of the strong optical emission remains inconclusive. Here we report polarimetric observations of three TDEs. The continuum polarization is independent of wavelength, while emission lines are partially depolarized. These signatures are consistent with optical photons being scattered and polarized in an envelope of free electrons. An almost axisymmetric photosphere viewed from different angles is in broad agreement with the data, but there is also evidence for deviations from axial symmetry before the peak of the flare and significant time evolution at early times, compatible with the rapid formation of an accretion disk. By combining a super-Eddington accretion model with a radiative transfer code we generate predictions for the degree of polarization as a function of disk mass and viewing angle, and we show that the predicted levels are compatible with the observations, for extended reprocessing envelopes of $\sim$1000 gravitational radii. Spectropolarimetry therefore constitutes a new observational test for TDE models, and opens an important new line of exploration in the study of TDEs.

preprint2022arXiv

Evidence for the preferential disruption of moderately massive stars by supermassive black holes

Tidal disruption events (TDEs) provide a unique opportunity to probe the stellar populations around supermassive black holes (SMBHs). By combining light curve modeling with spectral line information and knowledge about the stellar populations in the host galaxies, we are able to constrain the properties of the disrupted star for three TDEs. The TDEs in our sample have UV spectra, and measurements of the UV N III to C III line ratios enabled estimates of the nitrogen-to-carbon abundance ratios for these events. We show that the measured nitrogen line widths are consistent with originating from the disrupted stellar material dispersed by the central SMBH. We find that these nitrogen-to-carbon abundance ratios necessitate the disruption of moderately massive stars ($\gtrsim 1 - 2 M_\odot$). We determine that these moderately massive disruptions are over-represented by a factor of $\gtrsim 10^2$ when compared to the overall stellar population of the post-starburst galaxy hosts. This implies that SMBHs are preferentially disrupting higher mass stars, possibly due to ongoing top-heavy star formation in nuclear star clusters or to dynamical mechanisms that preferentially transport higher mass stars to their tidal radii.

preprint2022arXiv

HARM3D+NUC: A new method for simulating the post-merger phase of binary neutron star mergers with GRMHD, tabulated EOS and neutrino leakage

The first binary neutron star merger has already been detected in gravitational waves. The signal was accompanied by an electromagnetic counterpart including a kilonova component powered by the decay of radioactive nuclei, as well as a short $γ$-ray burst. In order to understand the radioactively-powered signal, it is necessary to simulate the outflows and their nucleosynthesis from the post-merger disk. Simulating the disk and predicting the composition of the outflows requires general relativistic magnetohydrodynamical (GRMHD) simulations that include a realistic, finite-temperature equation of state (EOS) and self-consistently calculating the impact of neutrinos. In this work, we detail the implementation of a finite-temperature EOS and the treatment of neutrinos in the GRMHD code HARM3D+NUC, based on HARM3D. We include formal tests of both the finite-temperature EOS and the neutrino leakage scheme. We further test the code by showing that, given conditions similar to those of published remnant disks following neutron star mergers, it reproduces both recombination of free nucleons to a neutron-rich composition and excitation of a thermal wind.

preprint2022arXiv

Mergers prompted by dynamics in compact, multiple-star systems: a stellar-reduction case for the massive triple TIC 470710327

TIC 470710327, a massive compact hierarchical triple-star system, was recently identified by NASA's Transiting Exoplanet Survey Satellite (TESS). TIC 470710327 is comprised of a compact (1.10 d) circular eclipsing binary, with total mass $\approx 10.9-13.2\ \rm{M_{\odot}}$, and a more massive ($\approx 14-17\ \rm{M_{\odot}}$) eccentric non-eclipsing tertiary in a $52.04$ d orbit. Here we present a progenitor scenario for TIC 470710327 in which '2+2' quadruple dynamics result in Zeipel-Lidov-Kozai (ZLK) resonances that lead to a contact phase of the more massive binary. In this scenario, the two binary systems should form in a very similar manner, and dynamics trigger the merger of the more massive binary either during late phases of star formation or several Myr after the zero-age main sequence (ZAMS), when the stars begin to expand. Any evidence that the tertiary is a highly-magnetised ($\sim 1-10$ kG), slowly-rotating blue main-sequence star would hint towards a quadruple origin. Finally, our scenario suggests that the population of inclined, compact multiple-stellar systems is reduced into co-planar systems, via mergers, late during star formation or early in the main sequence. The elucidation of the origin of TIC 470710327 is crucial in our understanding of multiple massive-star formation and evolution.

preprint2022arXiv

Successful Common Envelope Ejection and Binary Neutron Star Formation in 3D Hydrodynamics

A binary neutron star merger has been observed in a multi-messenger detection of gravitational wave (GW) and electromagnetic (EM) radiation. Binary neutron stars that merge within a Hubble time, as well as many other compact binaries, are expected to form via common envelope evolution. Yet five decades of research on common envelope evolution have not yet resulted in a satisfactory understanding of the multi-spatial multi-timescale evolution for the systems that lead to compact binaries. In this paper, we report on the first successful simulations of common envelope ejection leading to binary neutron star formation in 3D hydrodynamics. We simulate the dynamical inspiral phase of the interaction between a 12$M_\odot$ red supergiant and a 1.4$M_\odot$ neutron star for different initial separations and initial conditions. For all of our simulations, we find complete envelope ejection and final orbital separations of $a_{\rm f} \approx 1.3$-$5.1 R_\odot$ depending on the simulation and criterion, leading to binary neutron stars that can merge within a Hubble time. We find $α_{\rm CE}$-equivalent efficiencies of $\approx 0.1$-$2.7$ depending on the simulation and criterion, but this may be specific for these extended progenitors. We fully resolve the core of the star to $\lesssim 0.005 R_\odot$ and our 3D hydrodynamics simulations are informed by an adjusted 1D analytic energy formalism and a 2D kinematics study in order to overcome the prohibitive computational cost of simulating these systems. The framework we develop in this paper can be used to simulate a wide variety of interactions between stars, from stellar mergers to common envelope episodes leading to GW sources.

preprint2022arXiv

Supernova-driven turbulent metal mixing in high redshift galactic disks: metallicity fluctuations in the interstellar medium and its imprints on metal poor stars in the Milky Way

The extent to which turbulence mixes gas in the face of recurrent infusions of fresh metals by supernovae (SN) could help provide important constraints on the local star formation conditions. This includes predictions of the metallicity dispersion amongst metal poor stars, which suggests that the interstellar medium was not very well mixed at these early times. The purpose of this {\it Letter} is to help isolate, via a series of numerical experiments, some of the key processes that regulate turbulent mixing of SN elements in galactic disks. We study the gas interactions in small simulated patches of a galaxy disk with the goal of resolving the small-scale mixing effects of metals at pc scales, which enables us to measure the turbulent diffusion coefficient in various galaxy environments. By investigating the statistics of variations of $α$ elements in these simulations, we are able to derive constraints not only on the allowed range of intrinsic yield variations in SN explosions but also on the star formation history of the Milky Way. We argue that the observed dispersion of [Mg/Fe] in metal poor halo stars is compatible with the star-forming conditions expected in dwarf satellites or in an early low star-forming Milky Way progenitor. In particular, metal variations in stars that have not been phase-mixed can be used to infer the star-forming conditions of disrupted dwarf satellites.

preprint2022arXiv

Testing the Momentum-driven Supernova Feedback Paradigm in M31

Momentum feedback from isolated supernova remnants (SNRs) have been increasingly recognized by modern cosmological simulations as a resolution-independent means to implement the effects of feedback in galaxies, such as turbulence and winds. However, the integrated momentum yield from SNRs is uncertain due to the effects of SN clustering and interstellar medium (ISM) inhomogeneities. In this paper, we use spatially-resolved observations of the prominent 10-kpc star-forming ring of M31 to test models of mass-weighted ISM turbulence driven by momentum feedback from isolated, non-overlapping SNRs. We use a detailed stellar-age distribution (SAD) map from the Panchromatic Hubble Andromeda Treasury (PHAT) survey, observationally-constrained SN delay-time distributions, and maps of the atomic and molecular hydrogen to estimate the mass-weighted velocity dispersion using the Martizzi et al. ISM turbulence model. Our estimates are within a factor of 2 of the observed mass-weighted velocity dispersion in most of the ring, but exceed observations at densities $\lesssim 0.2$ cm$^{-3}$ and SN rates $>2.1\times 10^{-4}$ SN yr$^{-1}$ kpc$^{-2}$, even after accounting for plausible variations in stellar-age distribution models and ISM scale height assumptions. We conclude that at high SN rates the momentum deposited is most likely suppressed by the non-linear effects of SN clustering, while at low densities, SNRs reach pressure equilibrium before the cooling phase. These corrections should be introduced in models of momentum-driven feedback and ISM turbulence.

preprint2022arXiv

The Combined Effects of Two-Body Relaxation Processes and the Eccentric Kozai-Lidov Mechanism on the EMRI Rate

Gravitational wave (GW) emissions from extreme-mass-ratio inspirals (EMRIs) are promising sources for low-frequency GW-detectors. They result from a compact object, such as a stellar-mass black-hole (BH), captured by a supermassive black hole (SMBH). Several physical processes have been proposed to form EMRIs. In particular, weak two-body interactions over a long time scale (i.e., relaxation processes) have been proposed as a likely mechanism to drive the BH orbit to high eccentricity. Consequently, it is captured by the SMBH and becomes an EMRI. Here we demonstrate that EMRIs are naturally formed in SMBH binaries. Gravitational perturbations from an SMBH companion, known as the eccentric Kozai-Lidov (EKL) mechanism, combined with relaxation processes, yield a significantly more enhanced rate than any of these processes operating alone. Since EKL is sensitive to the orbital configuration, two-body relaxation can alter the orbital parameters, rendering the system in a more EKL-favorable regime. As SMBH binaries are expected to be prevalent in the Universe, this process predicts a substantially high EMRI rate.

preprint2021arXiv

A Carbon/Oxygen-dominated Atmosphere Days After Explosion for the "Super-Chandrasekhar" Type Ia SN 2020esm

Seeing pristine material from the donor star in a Type Ia supernova (SN Ia) explosion can reveal the nature of the binary system. In this paper, we present photometric and spectroscopic observations of SN 2020esm, one of the best-studied SNe of the class of "super-Chandrasekhar" SNe Ia (SC SNe Ia), with data obtained $-12$ to +360 days relative to peak brightness, obtained from a variety of ground- and space-based telescopes. Initially misclassified as a Type II supernova, SN 2020esm peaked at $M_{B} = -19.9$ mag, declined slowly ($Δm_{15}(B) = 0.92$ mag), and had particularly blue UV and optical colors at early times. Photometrically and spectroscopically, SN 2020esm evolved similarly to other SC SNe Ia, showing the usual low ejecta velocities, weak intermediate mass elements (IMEs), and the enhanced fading at late times, but its early spectra are unique. Our first few spectra (corresponding to a phase of $\gtrsim$10~days before peak) reveal a nearly-pure carbon/oxygen atmosphere during the first days after explosion. This composition can only be produced by pristine material, relatively unaffected by nuclear burning. The lack of H and He may further indicate that SN 2020esm is the outcome of the merger of two carbon/oxygen white dwarfs (WDs). Modeling its bolometric light curve, we find a $^{56}$Ni mass of $1.23^{+0.14}_{-0.14}$ M$_{\odot}$ and an ejecta mass of $1.75^{+0.32}_{-0.20}$ M$_{\odot}$, in excess of the Chandrasekhar mass. Finally, we discuss possible progenitor systems and explosion mechanisms of SN 2020esm and, in general, the SC SNe Ia class.

preprint2021arXiv

Follow-up Observations of the Prolonged, super-Eddington, Tidal Disruption Event Candidate 3XMM~J150052.0+015452: the Slow Decline Continues

The X-ray source 3XMM~J150052.0+015452 was discovered as a spectacular tidal disruption event candidate during a prolonged ($>11$ yrs) outburst (Lin et al. 2017). It exhibited unique quasi-soft X-ray spectra of characteristic temperature $kT\sim0.3$ keV for several years at the peak, but in a recent Chandra observation (10 yrs into the outburst) a super-soft X-ray spectrum of $kT\sim0.15$ keV was detected. Such dramatic spectral softening could signal the transition from the super-Eddington to thermal state or the temporary presence of a warm absorber. Here we report on our study of four new XMM-Newton follow-up observations of the source. We found that they all showed super-soft spectra, suggesting that the source had remained super-soft for $>5$ yrs. Then its spectral change is best explained as due to the super-Eddington to thermal spectral state transition. The fits to the thermal state spectra suggested a smaller absorption toward the source than that obtained in Lin et al. (2017). This led us to update the modeling of the event as due to the disruption of a 0.75 msun star by a massive black hole of a few$\times10^5$ msun. We also obtained two HST images in the F606W and F814W filters and found that the dwarf star-forming host galaxy can be resolved into a dominant disk and a smaller bulge. No central point source was clearly seen in either filter, ruling out strong optical emission associated with the X-ray activity.

preprint2021arXiv

Probing the progenitors of spinning binary black-hole mergers with long gamma-ray bursts

Long-duration gamma-ray bursts are thought to be associated with the core-collapse of massive, rapidly spinning stars and the formation of black holes. However, efficient angular momentum transport in stellar interiors, currently supported by asteroseismic and gravitational-wave constraints, leads to predominantly slowly-spinning stellar cores. Here, we report on binary stellar evolution and population synthesis calculations, showing that tidal interactions in close binaries not only can explain the observed sub-population of spinning, merging binary black holes but also lead to long gamma-ray bursts at the time of black-hole formation. Given our model calibration against the distribution of isotropic-equivalent energies of luminous long gamma-ray bursts, we find that ~10% of the GWTC-2 reported binary black holes had a luminous long gamma-ray burst associated with their formation, with GW190517 and GW190719 having a probability of ~85% and ~60%, respectively, being among them. Moreover, given an assumption about their average beaming fraction, our model predicts the rate density of long gamma-ray bursts, as a function of redshift, originating from this channel. For a constant beaming fraction $f_\mathrm{B}\sim 0.05$ our model predicts a rate density comparable to the observed one, throughout the redshift range, while, at redshift $z \in [0,2.5]$, a tentative comparison with the metallicity distribution of observed LGRB host galaxies implies that between 20% to 85% of the observed long gamma-ray bursts may originate from progenitors of merging binary black holes. The proposed link between a potentially significant fraction of observed, luminous long gamma-ray bursts and the progenitors of spinning binary black-hole mergers allows us to probe the latter well outside the horizon of current-generation gravitational wave observatories, and out to cosmological distances.

preprint2021arXiv

Radio and X-ray observations of the luminous Fast Blue Optical Transient AT2020xnd

We present deep X-ray and radio observations of the Fast Blue Optical Transient (FBOT) AT2020xnd/ZTF20acigmel at $z=0.2433$ from $13$d to $269$d after explosion. AT2020xnd belongs to the category of optically luminous FBOTs with similarities to the archetypal event AT2018cow. AT2020xnd shows luminous radio emission reaching $L_ν\approx8\times10^{29}$ergs$^{-1}$Hz$^{-1}$ at 20GHz and $75$d post explosion, accompanied by luminous and rapidly fading soft X-ray emission peaking at $L_{X}\approx6\times10^{42}$ergs$^{-1}$. Interpreting the radio emission in the context of synchrotron radiation from the explosion's shock interaction with the environment we find that AT2020xnd launched a high-velocity outflow ($v\sim$0.1-0.2$c$) propagating into a dense circumstellar medium (effective $\dot M\approx10^{-3}M_{\rm{sol}}$yr$^{-1}$ for an assumed wind velocity of $v_w=1000$kms$^{-1}$). Similar to AT2018cow, the detected X-ray emission is in excess compared to the extrapolated synchrotron spectrum and constitutes a different emission component, possibly powered by accretion onto a newly formed black hole or neutron star. These properties make AT2020xnd a high-redshift analog to AT2018cow, and establish AT2020xnd as the fourth member of the class of optically-luminous FBOTs with luminous multi-wavelength counterparts.

preprint2021arXiv

The Effects of Metallicity and Abundance Pattern of the ISM on Supernova Feedback

Supernova (SN) feedback plays a vital role in the evolution of galaxies. While modern cosmological simulations capture the leading structures within galaxies, they struggle to provide sufficient resolution to study small-scale stellar feedback, such as the detailed evolution of SN remnants. It is thus common practice to assume subgrid models that are rarely extended to low metallicities, and which routinely use the standard solar abundance pattern. With the aid of 1-d hydrodynamical simulations, we extend these models to consider low metallicities and non-solar abundance patterns as derived from spectra of Milky Way stars. For that purpose, a simple, yet effective framework has been developed to generate non-solar abundance pattern cooling functions. We find that previous treatments markedly over-predict SN feedback at low metallicities and show that non-negligible changes in the evolution of SN remnants of up to $\approx 50\%$ in $cooling\; mass$ and $\approx 27\%$ in $momentum\; injection\; from\; SN\; remnants$ arise from non-solar abundance patterns. We use our simulations to quantify these results as a function of metallicity and abundance pattern variations and present analytic formulae to accurately describe the trends. These formulae have been designed to serve as subgrid models for SN feedback in cosmological hydrodynamical simulations.

preprint2021arXiv

The Gravity Collective: A Search for the Electromagnetic Counterpart to the Neutron Star-Black Hole Merger GW190814

We present optical follow-up imaging obtained with the Katzman Automatic Imaging Telescope, Las Cumbres Observatory Global Telescope Network, Nickel Telescope, Swope Telescope, and Thacher Telescope of the LIGO/Virgo gravitational wave (GW) signal from the neutron star-black hole (NSBH) merger GW190814. We searched the GW190814 localization region (19 deg$^{2}$ for the 90th percentile best localization), covering a total of 51 deg$^{2}$ and 94.6% of the two-dimensional localization region. Analyzing the properties of 189 transients that we consider as candidate counterparts to the NSBH merger, including their localizations, discovery times from merger, optical spectra, likely host-galaxy redshifts, and photometric evolution, we conclude that none of these objects are likely to be associated with GW190814. Based on this finding, we consider the likely optical properties of an electromagnetic counterpart to GW190814, including possible kilonovae and short gamma-ray burst afterglows. Using the joint limits from our follow-up imaging, we conclude that a counterpart with an $r$-band decline rate of 0.68 mag day$^{-1}$, similar to the kilonova AT 2017gfo, could peak at an absolute magnitude of at most $-17.8$ mag (50% confidence). Our data are not constraining for &#39;&#39;red&#39;&#39; kilonovae and rule out &#39;&#39;blue&#39;&#39; kilonovae with $M>0.5 M_{\odot}$ (30% confidence). We strongly rule out all known types of short gamma-ray burst afterglows with viewing angles $<$17$^{\circ}$ assuming an initial jet opening angle of $\sim$$5.2^{\circ}$ and explosion energies and circumburst densities similar to afterglows explored in the literature. Finally, we explore the possibility that GW190814 merged in the disk of an active galactic nucleus, of which we find four in the localization region, but we do not find any candidate counterparts among these sources.

preprint2020arXiv

A trend in the effective spin distribution of LIGO binary black holes with mass

Binary black holes (BBHs) detected by gravitational wave (GW) observations could be broadly divided into two formation channels: those formed through field binary evolution and those assembled dynamically in dense stellar systems. Each of these formation channels, and their sub-channels, populate a distinct region in the effective spin-mass ($χ_{\rm eff}-M$) plane. Depending on the branching ratio of different channels, an ensemble of BBHs could show a trend in this plane. Here we fit a mass-dependent distribution for $χ_{\rm eff}$ to the GWTC-1 BBHs from the first and second observing runs of Advanced LIGO and Advanced Virgo. We find a negative correlation between mass and the mean effective spin ($\barχ_{\mathrm{eff}}$), and positive correlation with its dispersion ($σ_{χ_\mathrm{eff}}$) at 75\% and 80\% confidence. This trend is robust against the choice of mass variable, but most pronounced when the mass variable is taken to be the chirp mass of the binary. The result is consistent with significant contributions from both dynamically assembled and field binaries in the GWTC-1 catalog. The upcoming LIGO O3a data release will critically test this interpretation.

preprint2020arXiv

An Energy Inventory of Tidal Disruption Events

Tidal disruption events (TDEs) offer a unique opportunity to study a single super-massive black hole (SMBH) under feeding conditions that change over timescales of days or months. However, the primary mechanism for generating luminosity during the flares remains debated. Despite the increasing number of observed TDEs, it is unclear whether most of the energy in the initial flare comes from accretion near the gravitational radius or from circularizing debris at larger distances from the SMBH. The energy dissipation efficiency increases with decreasing radii, therefore by measuring the total energy emitted and estimating the efficiency we can derive clues about the nature of the emission mechanism. Here we calculate the integrated energy, emission timescales, and average efficiencies for the TDEs using the Modular Open Source Fitter for Transients ({\tt MOSFiT}). Our calculations of the total energy generally yield higher values than previous estimates. This is predominantly because, if the luminosity follows the mass fallback rate, TDEs release a significant fraction of their energy long after their light curve peaks. We use {\tt MOSFiT} to calculate the conversion efficiency from mass to radiated energy, and find that for many of the events it is similar to efficiencies inferred for active galactic nuclei. There are, however, large systematic uncertainties in the measured efficiency due to model degeneracies between the efficiency and the mass of the disrupted star, and these must be reduced before we can definitively resolve the emission mechanism of individual TDEs.

preprint2020arXiv

Carbon star formation as seen through the non-monotonic initial-final mass relation

The initial-final mass relation (IFMR) links the birth mass of a star to the mass of the compact remnant left at its death. While the relevance of the IFMR across astrophysics is universally acknowledged, not all of its fine details have yet been resolved. A new analysis of a few carbon-oxygen white dwarfs in old open clusters of the Milky Way led us to identify a kink in the IFMR, located over a range of initial masses, $1.65 \lesssim M_{\rm i}/M_{\odot} \lesssim 2.10$. The kink&#39;s peak in WD mass of $\approx 0.70-0.75 \, M_{\odot}$ is produced by stars with $M_{\rm i} \simeq 1.8 - 1.9 \, M_{\odot}$, corresponding to ages of about $1.8 - 1.7 $ Gyr. Interestingly, this peak coincides with the initial mass limit between low-mass stars that develop a degenerate helium core after central hydrogen exhaustion, and intermediate-mass stars that avoid electron degeneracy. We interpret the IFMR kink as the signature of carbon star formation in the Milky Way. This finding is critical to constraining the evolution and chemical enrichment of low-mass stars, and their impact on the spectrophotometric properties of galaxies.

preprint2020arXiv

Common Envelope Wind Tunnel: Range of Applicability and Self-Similarity in Realistic Stellar Envelopes

Common envelope evolution, the key orbital tightening phase of the traditional formation channel for close binaries, is a multistage process that presents many challenges to the establishment of a fully descriptive, predictive theoretical framework. In an approach complementary to global 3D hydrodynamical modeling, we explore the range of applicability for a simplified drag formalism that incorporates the results of local hydrodynamic &#34;wind tunnel&#34; simulations into a semi-analytical framework in the treatment of the common envelope dynamical inspiral phase using a library of realistic giant branch stellar models across the low, intermediate, and high mass regimes. In terms of a small number of key dimensionless parameters, we characterize a wide range of common envelope events, revealing the broad range of applicability of the drag formalism as well its self-similar nature across mass regimes and ages. Limitations arising from global binary properties and local structural quantities are discussed together with the opportunity for a general prescriptive application for this formalism.

preprint2020arXiv

Does GW190425 require an alternative formation pathway than a fast-merging channel?

The LIGO/Virgo Scientific Collaboration (LVC) recently announced the detection of a compact object binary merger, GW190425, with a total mass of $3.4^{+0.3}_{-0.1}$ M$_{\odot}$, and individual component masses in the range of about 1.1 to 2.5 $M_{\odot}$. If the constituent compact objects are neutron stars, then the total mass is five standard deviations higher than the mean of $2.66\pm 0.12$ M$_{\odot}$ for Galactic binary neutron stars (BNSs). LVC suggests such massive BNS systems are born from a fast-merging channel, and therefore, their non-detection in the Galaxy to be due to a selection effect. However, we are unable to reconcile the inferred formation efficiency from the reported merger rate, $\mathcal{R}_{\rm GW190425}=460^{+1050}_{-390}$ yr$^{-1}$ Gpc$^{-3}$, with predictions from our own study for fast-merging BNS systems. Moreover, the comparable merger rates of GW190425, and GW170817 are possibly in tension with our results for two reasons: (i) more massive systems are expected to have a lower formation rate, and (ii) fast merging channels should constitute $\lesssim 10\%$ of the total BNS systems if case BB unstable mass transfer is permitted to take place as a formation pathway. We argue that to account for the high merger rate of GW190425 as a BNS system requires: (i) revisiting our understanding of NS formation in supernova explosions, or (ii): that more massive NSs need to be preferentially born with either very weak or very high magnetic fields so that they would be undetectable in the radio surveys. Perhaps massive NSs detected in NS-white dwarf binaries are our clues to the formation path of GW190425 systems.

preprint2020arXiv

Massive Stellar Triples Leading to Sequential Binary Black-Hole Mergers in the Field

Stellar triples with massive stellar components are common, and can lead to sequential binary black-hole mergers. Here, we outline the evolution towards these sequential mergers, and explore these events in the context of gravitational-wave astronomy and the pair-instability mass gap. We find that binary black-hole mergers in the pair-instability mass gap can be of triple origin and therefore are not exclusively formed in dense dynamical environments. We discuss the sequential merger scenario in the context of the most massive gravitational-wave sources detected to date: GW170729 and GW190521. We propose that the progenitor of GW170729 is a low-metallicity field triple. We support the premise that GW190521 could not have been formed in the field. We conclude that triple stellar evolution is fundamental in the understanding of gravitational-wave sources, and likely, other energetic transientsas well.

preprint2020arXiv

Stellar Tidal Disruption Events with Abundances and Realistic Structures (STARS): Library of Fallback Rates

We present the STARS library, a grid of tidal disruption event (TDE) simulations interpolated to provide the mass fallback rate ($dM/dt$) to the black hole for a main-sequence star of any stellar mass, stellar age, and impact parameter. We use a one-dimensional stellar evolution code to construct stars with accurate stellar structures and chemical abundances, then perform tidal disruption simulations in a three-dimensional adaptive-mesh hydrodynamics code with a Helmholtz equation of state, in unprecedented resolution: from 131 to 524 cells across the diameter of the star. The interpolated library of fallback rates is available on GitHub (https://github.com/jamielaw-smith/STARS_library) and version 1.0.0 is archived on Zenodo; one can query the library for any stellar mass, stellar age, and impact parameter. We provide new fitting formulae for important disruption quantities ($β_{\rm crit}, ΔM, \dot M_{\rm peak}, t_{\rm peak}, n_\infty$) as a function of stellar mass, stellar age, and impact parameter. Each of these quantities vary significantly with stellar mass and stellar age, but we are able to reduce all of our simulations to a single relationship that depends only on stellar structure, characterized by a single parameter $ρ_c/\barρ$, and impact parameter $β$. We also find that, in general, more centrally concentrated stars have steeper $dM/dt$ rise slopes and shallower decay slopes. For the same $ΔM$, the $dM/dt$ shape varies significantly with stellar mass, promising the potential determination of stellar properties from the TDE light curve alone. The $dM/dt$ shape depends strongly on stellar structure and to a certain extent stellar mass, meaning that fitting TDEs using this library offers a better opportunity to determine the nature of the disrupted star and the black hole.

preprint2020arXiv

The Fate of the Merger Remnant in GW170817 and its Imprint on the Jet Structure

The first neutron star binary merger detected in gravitational waves, GW170817 and the subsequent detection of its emission across the electromagnetic spectrum showed that these systems are viable progenitors of short $γ$-ray bursts (sGRB). The afterglow signal of GW170817 has been found to be consistent with a structured GRB jet seen off-axis, requiring significant amounts of relativistic material at large angles. This trait can be attributed to the interaction of the relativistic jet with the external wind medium. Here we perform numerical simulations of relativistic jets interacting with realistic wind environments in order to explore how the properties of the wind and central engine affect the structure of successful jets. We find that the angular energy distribution of the jet depends primarily on the ratio between the lifetime of the jet and the time it takes the merger remnant to collapse. We make use of these simulations to constrain the time it took for the merger remnant in GW170817 to collapse into a black hole based on the angular structure of the jet as inferred from afterglow observations. We conclude that the lifetime of the merger remnant in GW170817 was $\approx 1-1.7$s, which, after collapse, triggered the formation of the jet.

preprint2020arXiv

The Stars in M15 Were Born with the r-process

High-resolution spectroscopy of stars on the red giant branch (RGB) of the globular cluster M15 has revealed a large (~1 dex) dispersion in the abundances of r-process elements, like Ba and Eu. Neutron star mergers (NSMs) have been proposed as a major source of the r-process. However, most NSM models predict a delay time longer than the timescale for cluster formation. One possibility is that a NSM polluted the surfaces of stars in M15 long after the cluster finished forming. In this case, the abundances of the polluting elements would decrease in the first dredge-up as stars turn on to the RGB. We present Keck/DEIMOS abundances of Ba in 66 stars along the entire RGB and the top of the main sequence. The Ba abundances have no trend with stellar luminosity (evolutionary phase). Therefore, the stars were born with the Ba they have today, and Ba did not originate in a source with a delay time longer than the timescale for cluster formation. In particular, if the source of Ba was a neutron star merger, it would have had a very short delay time. Alternatively, if Ba enrichment took place before the formation of the cluster, an inhomogeneity of a factor of 30 in Ba abundance needs to be able to persist over the length scale of the gas cloud that formed M15, which is unlikely.

preprint2020arXiv

Updated Parameter Estimates for GW190425 Using Astrophysical Arguments and Implications for the Electromagnetic Counterpart

The progenitor system of the compact binary merger GW190425 had a total mass of $3.4^{+0.3}_{-0.1}$ M$_\odot$ (90th-percentile confidence region) as measured from its gravitational wave signal. This mass is significantly different from the Milky Way (MW) population of binary neutron stars (BNSs) that are expected to merge in a Hubble time and from that of the first BNS merger, GW170817. Here we explore the expected electromagnetic signatures of such a system. We make several astrophysically motivated assumptions to further constrain the parameters of GW190425. By simply assuming that both components were NSs, we reduce the possible component masses significantly, finding $m_1 = 1.85^{+0.27}_{-0.19}$ M$_\odot$ and $m_2 = 1.47^{+0.16}_{-0.18}$ M$_\odot$. However if the GW190425 progenitor system was a NS-black hole merger, we find best-fitting parameters $m_1 = 2.19^{+0.21}_{-0.17}$ M$_\odot$ and $m_2 = 1.26^{+0.10}_{-0.08}$ M$_\odot$. For a well-motivated BNS system where the lighter NS has a mass similar to the mass of non-recycled NSs in MW BNS systems, we find $m_1 = 2.03^{+0.15}_{-0.14}$ M$_\odot$ and $m_2 = 1.35 \pm 0.09$ M$_\odot$, corresponding to only 7% mass uncertainties. For all scenarios, we expect a prompt collapse of the resulting remnant to a black hole. Examining detailed models with component masses similar to our best-fitting results, we find the electromagnetic counterpart to GW190425 is expected to be significantly redder and fainter than that of GW170817. We find that almost all reported observations used to search for an electromagnetic counterpart for GW190425 were too shallow to detect the expected counterpart. If the LIGO-Virgo Collaboration promptly provides the chirp mass, the astronomical community can adapt their observations to improve the likelihood of detecting a counterpart for similarly &#34;high-mass&#34; BNS systems. (abridged)

preprint2017arXiv

Low-mass White Dwarfs with Hydrogen Envelopes as a Missing Link in the Tidal Disruption Menu

We construct a menu of objects that can give rise to bright flares when disrupted by massive black holes (BHs), ranging from planets to evolved stars. Through their tidal disruption, main sequence and evolved stars can effectively probe the existence of otherwise quiescent supermassive BHs and white dwarfs can probe intermediate mass BHs. Many low-mass white dwarfs possess extended hydrogen envelopes, which allow for the production of prompt flares in disruptive encounters with moderately massive BHs of $10^5$ to $10^7~M_\odot$--masses that may constitute the majority of massive BHs by number. These objects are a missing link in two ways: (1) for probing moderately massive BHs and (2) for understanding the hydrodynamics of the disruption of objects with tenuous envelopes. A flare arising from the tidal disruption of a $0.17~M_\odot$ white dwarf by a $10^5~M_\odot$ BH reaches a maximum between 0.6 and 11 days, with a peak fallback rate that is usually super-Eddington and results in a flare that is likely brighter than a typical tidal disruption event. Encounters stripping only the envelope can provide hydrogen-only fallback, while encounters disrupting the core evolve from H- to He-rich fallback. While most tidal disruption candidates observed thus far are consistent with the disruptions of main sequence stars, the rapid timescales of nuclear transients such as Dougie and PTF10iya are naturally explained by the disruption of low-mass white dwarfs. As the number of observed flares continues to increase, the menu presented here will be essential for characterizing nuclear BHs and their environments through tidal disruptions.