Researcher profile

Rodrigo Fernández

Rodrigo Fernández contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
8works
0followers
7topics
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

8 published item(s)

preprint2025arXiv

Collapsar Disk Outflows III: Detectable Neutrino and Gravitational Wave Signatures

We investigate the neutrino and gravitational wave (GW) signals from accretion disks formed during the failed collapse of a rotating massive star (a collapsar). Following black hole formation, a neutrino-cooled, shocked accretion disk forms, which displays non-spherical oscillations for a period of seconds before becoming advective and exploding the star. We compute the neutrino and GW signals (matter quadrupole, frequencies $\lesssim 100$ Hz) from collapsar disks using global axisymmetric, viscous hydrodynamic simulations. The neutrino signal with typical energies of O$(10)$ MeV is maximal during the neutrino-cooled (NDAF) phase that follows shock formation. This phase lasts for a few seconds and is easily detectable within O$(10-100)$ kpc by the IceCube Neutrino Telescope. Additional neutrino signatures from a precursor equatorial shock and by stochastic accretion plumes during the advective phase are detectable within the galaxy. The GW signal during the NDAF phase is detectable in the galaxy by current and next-generation ground-based observatories. The explosion (memory) GW signal is similar to that of standard core-collapse supernovae and can be probed with a deci-Hertz space-based detector. Shock oscillations during the NDAF phase impart time variations with frequency O$(10-100)$ Hz to the neutrino and GW signals, encoding information about the shock dynamics and inner disk. These time variations can be detectable in neutrinos by IceCube within O$(1-10)$ kpc depending on progenitor model, flavor transformation scenario, and detailed properties of the angular momentum transport mechanism.

preprint2022arXiv

Long-term 3D-MHD Simulations of Black Hole Accretion Disks formed in Neutron Star Mergers

We examine the long-term evolution of accretion tori around black hole (BH) remnants of compact object mergers involving at least one neutron star, to better understand their contribution to kilonovae and the synthesis of r-process elements. To this end, we modify the unsplit magnetohydrodynamic (MHD) solver in FLASH4.5 to work in non-uniform three-dimensional spherical coordinates, enabling more efficient coverage of a large dynamic range in length scales while exploiting symmetries in the system. This modified code is used to perform BH accretion disk simulations that vary the initial magnetic field geometry and disk compactness, utilizing a physical equation of state, a neutrino leakage scheme for emission and absorption, and modeling the BH's gravity with a pseudo-Newtonian potential. Simulations run for long enough to achieve a radiatively-inefficient state in the disk. We find robust mass ejection with both poloidal and toroidal initial field geometries, and suppressed outflow at high disk compactness. With the included physics, we obtain bimodal velocity distributions that trace back to mass ejection by magnetic stresses at early times, and to thermal processes in the radiatively-inefficient state at late times. The electron fraction distribution of the disk outflow is broad in all models, and the ejecta geometry follows a characteristic hourglass shape. We test the effect of removing neutrino absorption or nuclear recombination with axisymmetric models, finding $\sim 50\%$ less mass ejection and more neutron-rich composition without neutrino absorption, and a subdominant contribution from nuclear recombination. Tests of the MHD and neutrino leakage implementations are included.

preprint2021arXiv

Probing magnetar emission mechanisms with spectropolarimetry

Over the next year, a new era of observations of compact objects in X-ray polarization will commence. Among the key targets for the upcoming Imaging X-ray Polarimetry Explorer mission, will be the magnetars 4U 0142+61 and 1RXS J170849.0-400910. Here we present the first detailed predictions of the expected polarization from these sources that incorporate realistic models of emission physics at the surface (gaseous or condensed), the temperature distribution on the surface, general relativity, quantum electrodynamics and scattering in the magnetosphere, and also account for the broadband spectral energy distribution of these sources from below 1 keV to nearly 100 keV. We find that either atmospheres or condensed surfaces can account for the emission at a few keV; in both cases either a small hot polar cap or scattering is required to account for the emission at 5-10 keV, and above 10 keV scattering by a hard population of electrons can account for the rising power in the hard X-rays observed in many magnetars in quiescence. Although these different scenarios result in very similar spectral energy distributions, they generate dramatically different polarization signatures from 2-10 keV, which is the range of sensitivity of the Imaging X-ray Polarimetry Explorer. Observations of these sources in X-ray polarization will therefore probe the emission from magnetars in an essentially new way.

preprint2020arXiv

A Deep CFHT Optical Search for a Counterpart to the Possible Neutron Star -- Black Hole Merger GW190814

We present a wide-field optical imaging search for electromagnetic counterparts to the likely neutron star - black hole (NS-BH) merger GW190814/S190814bv. This compact binary merger was detected through gravitational waves by the LIGO/Virgo interferometers, with masses suggestive of a NS-BH merger. We imaged the LIGO/Virgo localization region using the MegaCam instrument on the Canada-France-Hawaii Telescope. We describe our hybrid observing strategy of both tiling and galaxy-targeted observations, as well as our image differencing and transient detection pipeline. Our observing campaign produced some of the deepest multi-band images of the region between 1.7 and 8.7 days post-merger, reaching a 5sigma depth of g > 22.8 (AB mag) at 1.7 days and i > 23.1 and i > 23.9 at 3.7 and 8.7 days, respectively. These observations cover a mean total integrated probability of 67.0% of the localization region. We find no compelling candidate transient counterparts to this merger in our images, which suggests that either the lighter object was tidally disrupted inside of the BH&#39;s innermost stable circular orbit, the transient lies outside of the observed sky footprint, or the lighter object is a low-mass BH. We use 5sigma source detection upper limits from our images in the NS-BH interpretation of this merger to constrain the mass of the kilonova ejecta to be Mej < 0.015Msun for a &#39;blue&#39; (kappa = 0.5 cm^2 g^-1) kilonova, and Mej < 0.04Msun for a &#39;red&#39; (kappa = 5-10 cm^2 g^-1) kilonova. Our observations emphasize the key role of large-aperture telescopes and wide-field imagers such as CFHT MegaCam in enabling deep searches for electromagnetic counterparts to gravitational wave events.

preprint2020arXiv

The landscape of disk outflows from black hole - neutron star mergers

We investigate mass ejection from accretion disks formed in mergers of black holes (BHs) and neutron stars (NSs). The third observing run of the LIGO/Virgo interferometers provided BH-NS candidate events that yielded no electromagnetic (EM) counterparts. The broad range of disk configurations expected from BH-NS mergers motivates a thorough exploration of parameter space to improve EM signal predictions. Here we conduct 27 high-resolution, axisymmetric, long-term hydrodynamic simulations of the viscous evolution of BH accretion disks that include neutrino emission/absorption effects and post-processing with a nuclear reaction network. In the absence of magnetic fields, these simulations provide a lower-limit to the fraction of the initial disk mass ejected. We find a nearly linear inverse dependence of this fraction on disk compactness (BH mass over initial disk radius). The dependence is related to the fraction of the disk mass accreted before the outflow is launched, which depends on the disk position relative to the innermost stable circular orbit. We also characterize a trend of decreasing ejected fraction and decreasing lanthanide/actinide content with increasing disk mass at fixed BH mass. This trend results from a longer time to reach weak freezout and an increasingly dominant role of neutrino absorption at higher disk masses. We estimate the radioactive luminosity from the disk outflow alone available to power kilonovae over the range of configurations studied, finding a spread of two orders of magnitude. For most of the BH-NS parameter space, the disk outflow contribution is well below the kilonova mass upper limits for GW190814.

preprint2019arXiv

LRP2020: The cosmic origin and evolution of the elements

The origin of many elements of the periodic table remains an unsolved problem. While many nucleosynthetic channels are broadly understood, significant uncertainties remain regarding certain groups of elements such as the intermediate and rapid neutron-capture processes, the p-process, or the origin of odd-Z elements in the most metal-poor stars. Canada has a long tradition of leadership in nuclear astrophysics, dating back to the work of Alastair Cameron in the 1950s. Recent faculty hires have further boosted activity in the field, including transient observation and theory, survey science on galactic nucleosynthesis, and nuclear experiments. This white paper contains a brief overview of recent activity in the community, highlighting strengths in each sub-field, and provides recommendations to improve interdisciplinary collaboration. Sustaining Canadian leadership in the next decade will require, on the observational side, access to transient and non-transient surveys like LSST, SKA, or MSE, support for target-of-opportunity observing in current and future Canadian telescopes, and participation in next-generation X-ray telescopes such as ATHENA. State-of-the-art theoretical predictions will require an ambitious succession plan for the Niagara supercomputer to support large parallel jobs. We propose a funding instrument for postdoctoral training that reflects the interdisciplinary nature of nuclear astrophysics research, and the creation of a national collaborative funding program that allows for joint projects and workshop organization.

preprint2009arXiv

Stability of a Spherical Accretion Shock with Nuclear Dissociation

We examine the stability of a standing shock wave within a spherical accretion flow onto a gravitating star, in the context of core-collapse supernova explosions. Our focus is on the effect of nuclear dissociation below the shock on the linear growth, and non-linear saturation, of non-radial oscillations of the shocked fluid. We combine two-dimensional, time-dependent hydrodynamic simulations using FLASH2.5 with a solution to the linear eigenvalue problem, and demonstrate the consistency of the two approaches. Previous studies of this `Standing Accretion Shock Instability&#39; (SASI) have focused either on zero-energy accretion flows without nuclear dissociation, or made use of a detailed finite-temperature nuclear equation of state and included strong neutrino heating. Our main goal in this and subsequent papers is to introduce equations of state of increasing complexity, in order to isolate the various competing effects. In this work we employ an ideal gas equation of state with a constant rate of nuclear dissociation below the shock, and do not include neutrino heating. We find that a negative Bernoulli parameter below the shock significantly lowers the real frequency, growth rate, and saturation amplitude of the SASI. A decrease in the adiabatic index has similar effects. The non-linear development of the instability is characterized by an expansion of the shock driven by turbulent kinetic energy at nearly constant internal energy. Our results also provide further insight into the instability mechanism: the rate of growth of a particular mode is fastest when the radial advection time from the shock to the accretor overlaps with the period of a standing lateral sound wave. The fastest-growing mode can therefore be modified by nuclear dissociation.

preprint2005arXiv

Rotochemical Heating in Millisecond Pulsars. Formalism and Non-superfluid case

Rotochemical heating originates in a departure from beta equilibrium due to spin-down compression in a rotating neutron star. The main consequence is that the star eventually arrives at a quasi-equilibrium state, in which the thermal photon luminosity depends only on the current value of the spin-down power, which is directly measurable. Only in millisecond pulsars the spin-down power remains high long enough for this state to be reached with a substantial luminosity. We report an extensive study of the effect of this heating mechanism on the thermal evolution of millisecond pulsars, developing a general formalism in the slow-rotation approximation of general relativity that takes the spatial structure of the star fully into account, and using a sample of realistic equations of state to solve the non-superfluid case numerically. We show that nearly all observed millisecond pulsars are very likely to be in the quasi-equilibrium state. Our predicted quasi-equilibrium temperatures for PSR J0437-4715 are only 20% lower than inferred from observations. Accounting for superfluidity should increase the predicted value.