BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Oleg Korobkin

Oleg Korobkin appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
astro-ph.HEastro-ph.SRgr-qcnucl-thastro-ph.IMDistributed, Parallel, and Cluster ComputingSoftware Engineering

Catalog footprint

What is connected

16works
7topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

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

Published work

16 published item(s)

preprint2022arXiv

Realistic Kilonova Up Close

Neutron star mergers are cosmic catastrophes that produce some of the most energetic observed phenomena: short gamma-ray bursts, gravitational wave signals, and kilonovae. The latter are optical transients, powered by radioactive nuclides which are synthesized when the neutron-rich ejecta of a disrupted neutron star undergoes decompression. We model this decompression phase using data from simulations of post-merger accretion disk winds. We use smoothed particle hydrodynamics with realistic nuclear heating to model the expansion over multiple scales, from initially several thousand km to billions of km. We then render a realistic image of a kilonova ejecta as it would appear for a nearby observer. This is the first time such a visualization is performed using input from state-of-the-art accretion disk simulations, nuclear physics and atomic physics. The volume rendering in our model computes an opacity transfer function on the basis of the physical opacity, varying significantly with the inhomogeneity of the neutron richness in the ejecta. Other physical quantities such as temperature or electron fraction can be visualized using an independent color transfer function. We discuss several difficulties with the ParaView application that we encountered during the visualization process, and give descriptions of our solutions and workarounds which could be used for future improvements.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Kilonova Detectability with Wide-Field Instruments

Kilonovae are ultraviolet, optical, and infrared transients powered by the radioactive decay of heavy elements following a neutron star merger. Joint observations of kilonovae and gravitational waves can offer key constraints on the source of Galactic $r$-process enrichment, among other astrophysical topics. However, robust constraints on heavy element production requires rapid kilonova detection (within $\sim 1$ day of merger) as well as multi-wavelength observations across multiple epochs. In this study, we quantify the ability of 13 wide field-of-view instruments to detect kilonovae, leveraging a large grid of over 900 radiative transfer simulations with 54 viewing angles per simulation. We consider both current and upcoming instruments, collectively spanning the full kilonova spectrum. The Roman Space Telescope has the highest redshift reach of any instrument in the study, observing kilonovae out to $z \sim 1$ within the first day post-merger. We demonstrate that BlackGEM, DECam, GOTO, the Vera C. Rubin Observatory's LSST, ULTRASAT, and VISTA can observe some kilonovae out to $z \sim 0.1$ ($\sim$475 Mpc), while DDOTI, MeerLICHT, PRIME, $Swift$/UVOT, and ZTF are confined to more nearby observations. Furthermore, we provide a framework to infer kilonova ejecta properties following non-detections and explore variation in detectability with these ejecta parameters.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Gamma-rays from kilonova: a potential probe of r-process nucleosynthesis

The mergers of compact binaries with at least one neutron star component are the potential leading sites of the production and ejection of $r$-process elements. Discoveries of galactic binary pulsars, short gamma-ray bursts, and gravitational-wave detections have all been constraining the rate of these events, while the gravitational wave plus broadband electromagnetic coverage of binary neutron star merger (GW170817) has also placed constraints on the properties (mass and composition) of the merger ejecta. But uncertainties and ambiguities in modeling the optical and infrared emission make it difficult to definitively measure the distribution of heavy isotopes in these mergers. In contrast, gamma rays emitted in the decay of these neutron-rich ejecta may provide a more direct measurement of the yields. We calculate the gamma production in remnants of neutron star mergers, considering two epochs: a kilonova epoch, lasting about two weeks, and a much later epoch of tens and hundreds of thousands of years after the merger. For the kilonova epoch, when the expanding ejecta is still only partially transparent to gamma radiation, we use 3D radiative transport simulations to produce the spectra. We show that the gamma-ray spectra associated with beta- and alpha-decay provide a fingerprint of the ejecta properties and, for a sufficiently nearby remnant, may be detectable, even for old remnants. We compare our gamma spectra with the potential detection limits of next generation detectors, including the Lunar Occultation Explorer (LOX), the All-sky Medium Energy Gamma-ray Observatory (AMEGO), and the Compton Spectrometer and Imager (COSI). We show that fission models can be discriminated via the presence of short-lived fission fragments in the remnant spectra.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Gamma-rays from kilonovae and the cosmic gamma-ray background

The recent detection of the gravitational wave event GW170817, produced by the coalescence of two neutron stars, and of its optical-infrared counterpart, powered by the radioactive decay of r-process elements, has opened a new window to gamma-ray astronomy: the direct detection of photons coming from such decays. Here we calculate the contribution of kilonovae to the diffuse gamma-ray background in the MeV range, using recent results on the spectra of the gamma-rays emitted in individual events, and we compare it with that from other sources. We find that the contribution from kilonovae is not dominant in such energy range, but within current uncertainties, its addition to other sources might help to fit the observational data.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Composition Effects on Kilonova Spectra and Light Curves: I

The merger of neutron star binaries is believed to eject a wide range of heavy elements into the universe. By observing the emission from this ejecta, scientists can probe the ejecta properties (mass, velocity and composition distributions). The emission (a.k.a. kilonova) is powered by the radioactive decay of the heavy isotopes produced in the merger and this emission is reprocessed by atomic opacities to optical and infra-red wavelengths. Understanding the ejecta properties requires calculating the dependence of this emission on these opacities. The strong lines in the optical and infra-red in lanthanide opacities have been shown to significantly alter the light-curves and spectra in these wavelength bands, arguing that the emission in these wavelengths can probe the composition of this ejecta. Here we study variations in the kilonova emission by varying individual lanthanide (and the actinide uranium) concentrations in the ejecta. The broad forest of lanthanide lines makes it difficult to determine the exact fraction of individual lanthanides. Nd is an exception. Its opacities above 1 micron are higher than other lanthanides and observations of kilonovae can potentially probe increased abundances of Nd. Similarly, at early times when the ejecta is still hot (first day), the U opacity is strong in the 0.2-1 micron wavelength range and kilonova observations may also be able to constrain these abundances.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Neutrino-driven winds in the aftermath of a neutron star merger: nucleosynthesis and electromagnetic transients

We present a comprehensive nucleosynthesis study of the neutrino-driven wind in the aftermath of a binary neutron star merger. Our focus is the initial remnant phase when a massive central neutron star is present. Using tracers from a recent hydrodynamical simulation, we determine total masses and integrated abundances to characterize the composition of unbound matter. We find that the nucleosynthetic yields depend sensitively on both the life time of the massive neutron star and the polar angle. Matter in excess of up to $9 \cdot 10^{-3} M_\odot$ becomes unbound until $\sim 200~{\rm ms}$. Due to electron fractions of $Y_{\rm e} \approx 0.2 - 0.4$ mainly nuclei with mass numbers $A < 130$ are synthesized, complementing the yields from the earlier dynamic ejecta. Mixing scenarios with these two types of ejecta can explain the abundance pattern in r-process enriched metal-poor stars. Additionally, we calculate heating rates for the decay of the freshly produced radioactive isotopes. The resulting light curve peaks in the blue band after about $4~{\rm h}$. Furthermore, high opacities due to heavy r-process nuclei in the dynamic ejecta lead to a second peak in the infrared after $3-4~{\rm d}$.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Nucleosynthesis in the Ejecta of Neutron Star Mergers

Heavy elements like gold, platinum or uranium are produced in the r-process, which needs neutron-rich and explosive environments. Neutron star mergers are a promising candidate for an r-process site. They exhibit three different channels for matter ejection fulfilling these conditions: dynamic ejecta due to tidal torques, neutrino-driven winds and evaporating matter from the accretion disk. We present a first study of the integrated nucleosynthesis for a neutrino-driven wind from a neutron star merger with a hyper-massive neutron star. Trajectories from a recent hydrodynamical simulation are divided into four different angle regions and post-processed with a reaction network. We find that the electron fraction varies around $Y_e \approx 0.1 - 0.4$, but its distribution differs for every angle of ejection. Hence, the wind ejecta do not undergo a robust r-process, but rather possess distinct nucleosynthesis yields depending on the angle range. Compared to the dynamic ejecta, a smaller amount of neutron-rich matter gets unbound, but the production of lighter heavy elements with $A \lesssim 130$ in the neutrino-driven wind can complement the strong r-process of the dynamic ejecta.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Relativistic effects on tidal disruption kicks of solitary stars

Solitary stars that wander too close to their galactic centres can become tidally disrupted, if the tidal forces due to the supermassive black hole (SMBH) residing there overcome the self-gravity of the star. If the star is only partially disrupted, so that a fraction survives as a self-bound object, this remaining core will experience a net gain in specific orbital energy, which translates into a velocity "kick" of up to $\sim 10^3$ km/s. In this paper, we present the result of smoothed particle hydrodynamics (SPH) simulations of such partial disruptions, and analyse the velocity kick imparted on the surviving core. We compare $γ$ = 5/3 and $γ$ = 4/3 polytropes disrupted in both a Newtonian potential, and a generalized potential that reproduces most relativistic effects around a Schwarzschild black hole either exactly or to excellent precision. For the Newtonian case, we confirm the results of previous studies that the kick velocity of the surviving core is virtually independent of the ratio of the black hole to stellar mass, and is a function of the impact parameter $β$ alone, reaching at most the escape velocity of the original star. For a given $β$, relativistic effects become increasingly important for larger black hole masses. In particular, we find that the kick velocity increases with the black hole mass, making larger kicks more common than in the Newtonian case, as low-$β$ encounters are statistically more likely than high-$β$ encounters. The analysis of the tidal tensor for the generalized potential shows that our results are robust lower limits on the true relativistic kick velocities, and are generally in very good agreement with the exact results.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The Role of Fission in Neutron Star Mergers and its Impact on the r-Process Peaks

Comparing observational abundance features with nucleosynthesis predictions of stellar evolution or explosion simulations can scrutinize two aspects: (a) the conditions in the astrophysical production site and (b) the quality of the nuclear physics input utilized. We test the abundance features of r-process nucleosynthesis calculations for the dynamical ejecta of neutron star merger simulations based on three different nuclear mass models: The Finite Range Droplet Model (FRDM), the (quenched version of the) Extended Thomas Fermi Model with Strutinsky Integral (ETFSI-Q), and the Hartree-Fock-Bogoliubov (HFB) mass model. We make use of corresponding fission barrier heights and compare the impact of four different fission fragment distribution models on the final r-process abundance distribution. In particular, we explore the abundance distribution in the second r-process peak and the rare-earth sub-peak as a function of mass models and fission fragment distributions, as well as the origin of a shift in the third r-process peak position. The latter has been noticed in a number of merger nucleosynthesis predictions. We show that the shift occurs during the r-process freeze-out when neutron captures and β-decays compete and an (n,γ)-(γ,n) equilibrium is not maintained anymore. During this phase neutrons originate mainly from fission of material above A = 240. We also investigate the role of β-decay half-lives from recent theoretical advances, which lead either to a smaller amount of fissioning nuclei during freeze-out or a faster (and thus earlier) release of fission neutrons, which can (partially) prevent this shift and has an impact on the second and rare-earth peak as well.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Implications of GRB 130603B and its macronova for r-process nucleosynthesis

The tentative identification of a Li-Paczynski macronova following the short GRB 130603B indicated that a few hundredths of a solar mass of neutron star matter were ejected and that this ejected mass has radioactively decayed into heavy r-process elements. If correct, this confirms long standing predictions (Eichler, Livio, Piran and Schramm 1989) that on the one hand, sGRBs are produced in compact binary mergers (CBMs) and on the other hand that these events are significant and possibly dominant sources of the heavy (A>130) r-process nuclei. Assuming that this interpretation is correct we obtain a lower limit of 0.02 m_sun on the ejected mass. Using the current estimates of the rate of sGRBs and with a beaming factor of 50, mergers associated with sGRBs can produce all the observed heavy r-process material in the Universe. We confront this conclusion with cosmochemistry and show that even though such events are rare, mixing is sufficient to account for the current homogeneous distribution of r-process material in the Galaxy. However, the appearance of significant amounts of Eu in some very low metallicity stars requires that some mergers took place very early on, namely with a very short time delay after the earliest star formation episodes. Alternatively, an additional early r-processsource may have contributed at that early stage. Finally, we note that evidence for short lived {244}^Pu in the very early solar system suggests that a merger of this kind took place within the vicinity of the solar system shortly (a few hundred million years) before its formation.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Neutrino-driven winds from neutron star merger remnants

We present a detailed, 3D hydrodynamics study of the neutrino-driven winds that emerge from the remnant of a NS merger. Our simulations are performed with the Newtonian, Eulerian code FISH, augmented by a detailed, spectral neutrino leakage scheme that accounts for heating due to neutrino absorption in optically thin conditions. Consistent with the 2D study of Dessart et al. (2009), we find that a strong baryonic wind is blown out along the original binary rotation axis within $100$ ms after the merger. We compute a lower limit on the expelled mass of $3.5 \times 10^{-3} M_{\odot}$, large enough to be relevant for heavy element nucleosynthesis. The physical properties vary significantly between different wind regions. For example, due to stronger neutrino irradiation, the polar regions show substantially larger $Y_e$ than those at lower latitudes. This has its bearings on the nucleosynthesis: the polar ejecta produce interesting r-process contributions from $A\sim 80$ to about 130, while the more neutron-rich, lower-latitude parts produce also elements up to the third r-process peak near $A\sim 195$. We also calculate the properties of electromagnetic transients that are powered by the radioactivity in the wind, in addition to the macronova transient that stems from the dynamic ejecta. The high-latitude (polar) regions produce UV/optical transients reaching luminosities up to $10^{41} {\rm erg \, s^{-1}}$, which peak around 1 day in optical and 0.3 days in bolometric luminosity. The lower-latitude regions, due to their contamination with high-opacity heavy elements, produce dimmer and more red signals, peaking after $\sim 2$ days in optical and infrared. Our numerical experiments indicate that it will be difficult to infer the collapse time-scale of the HMNS to a BH based on the wind electromagnetic transient, at least for collapse time-scales larger than the wind production time-scale.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The long-term evolution of neutron star merger remnants - II. Radioactively powered transients

We use 3D hydrodynamic simulations of the long-term evolution of neutron star merger ejecta to predict the light curves of electromagnetic transients that are powered by the decay of freshly produced r-process nuclei. For the dynamic ejecta that are launched by tidal and hydrodynamic interaction, we adopt grey opacities of 10 cm$^2$/g, as suggested by recent studies. For our reference case of a 1.3-1.4 $M_\odot$ merger, we find a broad IR peak 2-4 d after the merger. The peak luminosity is $\approx 2\times 10^{40}$ erg/s for an average orientation, but increased by up to a factor of 4 for more favourable binary parameters and viewing angles. These signals are rather weak and hardly detectable within the large error box (~100 deg$^2$) of a gravitational wave trigger. A second electromagnetic transient results from neutrino-driven winds. These winds produce `weak' r-process material with $50 < A < 130$ and abundance patterns that vary substantially between different merger cases. For an adopted opacity of 1 cm$^2$/g, the resulting transients peak in the UV/optical about 6 h after the merger with a luminosity of $\approx 10^{41}$ erg/s (for a wind of 0.01 $M_\odot$) These signals are marginally detectable in deep follow-up searches (e.g. using Hypersuprime camera on Subaru). A subsequent detection of the weaker but longer lasting IR signal would allow an identification of the merger event. We briefly discuss the implications of our results to the recent detection of an nIR transient accompanying GRB 130603B.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Ensuring Correctness at the Application Level: a Software Framework Approach

As scientific applications extend to the simulation of more and more complex systems, they involve an increasing number of abstraction levels, at each of which errors can emerge and across which they can propagate; tools for correctness evaluation and enforcement at every level (from the code level to the application level) are therefore necessary. Whilst code-level debugging tools are already a well established standard, application-level tools are lagging behind, possibly due to their stronger dependence on the application's details. In this paper, we describe the programming model introduced by the Cactus framework, review the High Performance Computing (HPC) challenges that Cactus is designed to address, and illustrate the correctness strategies that are currently available in Cactus at the code, component, and application level.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

On the frequency band of the f-mode CFS instability

Rapidly rotating neutron stars can be unstable to the gravitational-wave-driven CFS mechanism if they have a neutral point in the spectrum of nonaxisymmetric f-modes. We investigate the frequencies of these modes in two sequences of uniformly rotating polytropes using nonlinear simulations in full general relativity, determine the approximate locations of the neutral points, and derive limits on the observable frequency band available to the instability in these sequences. We find that general relativity enhances the detectability of a CFS-unstable neutron star substantially, both by widening the instability window and enlarging the band into the optimal range for interferometric detectors like LIGO, VIRGO, and GEO-600.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Stability of general-relativistic accretion disks

Self-gravitating relativistic disks around black holes can form as transient structures in a number of astrophysical scenarios such as binary neutron star and black hole-neutron star coalescences, as well as the core-collapse of massive stars. We explore the stability of such disks against runaway and non-axisymmetric instabilities using three-dimensional hydrodynamics simulations in full general relativity using the THOR code. We model the disk matter using the ideal fluid approximation with a $Γ$-law equation of state with $Γ=4/3$. We explore three disk models around non-rotating black holes with disk-to-black hole mass ratios of 0.24, 0.17 and 0.11. Due to metric blending in our initial data, all of our initial models contain an initial axisymmetric perturbation which induces radial disk oscillations. Despite these oscillations, our models do not develop the runaway instability during the first several orbital periods. Instead, all of the models develop unstable non-axisymmetric modes on a dynamical timescale. We observe two distinct types of instabilities: the Papaloizou-Pringle and the so-called intermediate type instabilities. The development of the non-axisymmetric mode with azimuthal number m = 1 is accompanied by an outspiraling motion of the black hole, which significantly amplifies the growth rate of the m = 1 mode in some cases. Overall, our simulations show that the properties of the unstable non-axisymmetric modes in our disk models are qualitatively similar to those in Newtonian theory.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Solving the Einstein constraint equations on multi-block triangulations using finite element methods

In order to generate initial data for nonlinear relativistic simulations, one needs to solve the Einstein constraints, which can be cast into a coupled set of nonlinear elliptic equations. Here we present an approach for solving these equations on three-dimensional multi-block domains using finite element methods. We illustrate our approach on a simple example of Brill wave initial data, with the constraints reducing to a single linear elliptic equation for the conformal factor $ψ$. We use quadratic Lagrange elements on semi-structured simplicial meshes, obtained by triangulation of multi-block grids. In the case of uniform refinement the scheme is superconvergent at most mesh vertices, due to local symmetry of the finite element basis with respect to local spatial inversions. We show that in the superconvergent case subsequent unstructured mesh refinements do not improve the quality of our initial data. As proof of concept that this approach is feasible for generating multi-block initial data in three dimensions, after constructing the initial data we evolve them in time using a high order finite-differencing multi-block approach and extract the gravitational waves from the numerical solution.

0 citations0 reviewsNo code linkedOpen access