Researcher profile

F. K. Roepke

F. K. Roepke contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
14works
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

14 published item(s)

preprint2025arXiv

NLTE spectral modelling for a carbon-oxygen and helium white-dwarf merger as a Ca-rich transient candidate

We carry out NLTE (non local thermodynamic equilibrium) radiative transfer simulations to determine whether explosion during the merger of a carbon-oxygen (CO) white dwarf (WD) with a helium (He) WD can reproduce the characteristic Ca II/[Ca II] and He I lines observed in Ca-rich transients. Our study is based on a 1D representation of a hydrodynamic simulation of a 0.6 $M_{\odot}$ CO + 0.4 $M_{\odot}$ He WD merger. We calculate both photospheric and nebular-phase spectra including treatment for non-thermal electrons, as is required for accurate modelling of He I and [Ca II]. Consistent with Ca-rich transients, our simulation predicts a nebular spectrum dominated by emission from [Ca II] 7291, 7324 angstrom and the Ca II near-infrared (NIR) triplet. The photospheric-phase synthetic spectrum also exhibits a strong Ca II NIR triplet, prominent optical absorption due to He I 5876 angstrom and He I 10830 angstrom in the NIR, as is commonly observed for Ca-rich transients. Overall, our results therefore suggest that CO+He WD mergers are a promising channel for Ca-rich transients. However, the current simulation overpredicts some He I features, in particular both He I 6678 and 7065 angstrom and shows a significant contribution from Ti II, which results in a spectral energy distribution that is substantially redder than most Ca-rich transients at peak. Additionally the Ca II nebular emission features are too broad. Future work should investigate if these discrepancies can be resolved by considering full 3D models and exploring a range of CO+He WD binary configurations.

preprint2022arXiv

Bipolar planetary nebulae from common envelope evolution of binary stars

Asymmetric shapes and evidence for binary central stars suggest a common-envelope origin for many bipolar planetary nebulae. The bipolar components of the nebulae are observed to expand faster than the rest and the more slowly expanding material has been associated with the bulk of the envelope ejected during the common-envelope phase of a stellar binary system. Common-envelope evolution in general remains one of the biggest uncertainties in binary star evolution and the origin of the fast outflow has not been explained satisfactorily. We perform three-dimensional magnetohydrodynamic simulations of common-envelope interaction with the moving-mesh code AREPO. Starting from the plunge-in of the companion into the envelope of an asymptotic giant branch star and covering hundreds of orbits of the binary star system, we are able to follow the evolution to complete envelope ejection. We find that magnetic fields are strongly amplified in two consecutive episodes. First, when the companion spirals in the envelope and, second, when it forms a contact binary with the core of the former giant star. In the second episode, a magnetically-driven, high-velocity outflow of gas is launched self-consistently in our simulations. The outflow is bipolar and the gas is additionally collimated by the ejected common envelope. The resulting structure reproduces typical morphologies and velocities observed in young planetary nebulae. We propose that the magnetic driving mechanism is a universal consequence of common envelope interaction responsible for a substantial fraction of observed planetary nebulae. Such a mechanism likely also exists in the common-envelope phase of other binary stars that lead to the formation of Type Ia supernovae, X-ray binaries and gravitational-wave merger events.

preprint2022arXiv

Dynamics in a stellar convective layer and at its boundary: Comparison of five 3D hydrodynamics codes

Our ability to predict the structure and evolution of stars is in part limited by complex, 3D hydrodynamic processes such as convective boundary mixing. Hydrodynamic simulations help us understand the dynamics of stellar convection and convective boundaries. However, the codes used to compute such simulations are usually tested on extremely simple problems and the reliability and reproducibility of their predictions for turbulent flows is unclear. We define a test problem involving turbulent convection in a plane-parallel box, which leads to mass entrainment from, and internal-wave generation in, a stably stratified layer. We compare the outputs from the codes FLASH, MUSIC, PPMSTAR, PROMPI, and SLH, which have been widely employed to study hydrodynamic problems in stellar interiors. The convection is dominated by the largest scales that fit into the simulation box. All time-averaged profiles of velocity components, fluctuation amplitudes, and fluxes of enthalpy and kinetic energy are within $\lesssim 3σ$ of the mean of all simulations on a given grid ($128^3$ and $256^3$ grid cells), where $σ$ describes the statistical variation due to the flow's time dependence. They also agree well with a $512^3$ reference run. The $128^3$ and $256^3$ simulations agree within $9\%$ and $4\%$, respectively, on the total mass entrained into the convective layer. The entrainment rate appears to be set by the amount of energy that can be converted to work in our setup and details of the small-scale flows in the boundary layer seem to be largely irrelevant. Our results lend credence to hydrodynamic simulations of flows in stellar interiors. We provide in electronic form all outputs of our simulations as well as all information needed to reproduce or extend our study.

preprint2022arXiv

Models of pulsationally assisted gravitationally confined detonations with different ignition conditions

Over the past decades, many explosion scenarios for Type Ia supernovae have been proposed and investigated including various combinations of deflagrations and detonations in white dwarfs of different masses up to the Chandrasekhar mass. One of these is the gravitationally confined detonation model. In this case a weak deflagration burns to the surface, wraps around the bound core, and collides at the antipode. A subsequent detonation is then initiated in the collision area. Since the parameter space for this scenario, that is, varying central densities and ignition geometries, has not been studied in detail, we used pure deflagration models of a previous parameter study dedicated to Type Iax supernovae as initial models to investigate the gravitationally confined detonation scenario. We aim to judge whether this channel can account for one of the many subgroups of Type Ia supernovae, or even normal events. To this end, we employed a comprehensive pipeline for three-dimensional Type Ia supernova modeling that consists of hydrodynamic explosion simulations, nuclear network calculations, and radiative transfer. The observables extracted from the radiative transfer are then compared to observed light curves and spectra. The study produces a wide range in masses of synthesized 56 Ni ranging from 0.257 to 1.057 $M_\odot$ , and, thus, can potentially account for subluminous as well as overluminous Type Ia supernovae in terms of brightness. However, a rough agreement with observed light curves and spectra can only be found for 91T-like objects. Although several discrepancies remain, we conclude that the gravitationally confined detonation model cannot be ruled out as a mechanism to produce 91T-like objects. However, the models do not provide a good explanation for either normal Type Ia supernovae or Type Iax supernovae.

preprint2021arXiv

Type Iax supernovae from deflagrations in Chandrasekhar mass white dwarfs

Due to the increasing number of observations Type Ia supernovae are nowadays regarded as a heterogeneous class of objects consisting of several subclasses. One of the largest of these is the class of Type Iax supernovae (SNe Iax) which have been suggested to originate from pure deflagrations in CO Chandrasekhar-mass white dwarfs (WDs). Although a few deflagration studies have been carried out, the full diversity of the class is not captured yet. We therefore present a parameter study of single-spot ignited deflagrations with varying ignition locations, central densities, metallicities and compositions. We also explore a rigidly rotating progenitor and carry out 3D hydrodynamic simulations, nuclear network calculations and radiative transfer. The new models extend the range in brightness covered by previous studies to the lower end. Our explosions produce $^{56}$Ni masses from $5.8 \times 10^{-3}$ to $9.2 \times 10^{-2}\,M_\odot$. In spite of the wide exploration of the parameter space the main characteristics of the models are primarily driven by the mass of $^{56}$Ni. Secondary parameters have too little impact to explain the observed trend among faint SNe~Iax. We report kick velocities of the bound explosion remnants from $6.9$ to $369.8\,$km$\,s^{-1}$. The wide exploration of the parameter space and viewing-angle effects in the radiative transfer lead to a significant spread in the synthetic observables. The trends towards the faint end of the class are, however, not reproduced. This motivates a quantification of the systematic uncertainties in the modeling procedure and the influence of the $^{56}$Ni-rich bound remnant. While the pure deflagration scenario remains a favorable explanation for bright and intermediate luminosity SNe~Iax, the possibility that SNe~Iax do not consist of a single explosion scenario needs to be considered.

preprint2020arXiv

Common-envelope evolution with an asymptotic giant branch star

Common-envelope phases are decisive for the evolution of many binary systems. Of particular interest are cases with asymptotic giant branch (AGB) primary stars, because they are thought to be progenitors of various astrophysical transients. In three-dimensional hydrodynamic simulations with the moving-mesh code AREPO, we study the common-envelope evolution of a $1.0\,M_{\odot}$ early-AGB star with companions of different masses. Although the stellar envelope of the AGB star is less tightly bound than that of a red giant, we find that the release of orbital energy of the core binary is insufficient to eject more than about twenty percent of the envelope mass. Ionization energy released in the expanding envelope, however, can lead to complete envelope ejection. Because recombination proceeds largely at high optical depths in our simulations, it is likely that this effect indeed plays a significant role in the considered systems. The efficiency of mass loss and the final orbital separation of the core binary system depend on the mass ratio between the companion and the primary star. Our results suggest a linear relation between the ratio of final to initial orbital separation and this parameter.

preprint2020arXiv

Fully compressible simulations of waves and core convection in main-sequence stars

Context. Recent, nonlinear simulations of wave generation and propagation in full-star models have been carried out in the anelastic approximation using spectral methods. Although it makes long time steps possible, this approach excludes the physics of sound waves completely and rather high artificial viscosity and thermal diffusivity are needed for numerical stability. Direct comparison with observations is thus limited. Aims. We explore the capabilities of our compressible multidimensional hydrodynamics code SLH to simulate stellar oscillations. Methods. We compare some fundamental properties of internal gravity and pressure waves in 2D SLH simulations to linear wave theory using two test cases: (1) an interval gravity wave packet in the Boussinesq limit and (2) a realistic $3\mathrm{M}_\odot$ stellar model with a convective core and a radiative envelope. Oscillation properties of the stellar model are also discussed in the context of observations. Results. Our tests show that specialized low-Mach techniques are necessary when simulating oscillations in stellar interiors. Basic properties of internal gravity and pressure waves in our simulations are in good agreement with linear wave theory. As compared to anelastic simulations of the same stellar model, we can follow internal gravity waves of much lower frequencies. The temporal frequency spectra of velocity and temperature are flat and compatible with observed spectra of massive stars. Conclusion. The low-Mach compressible approach to hydrodynamical simulations of stellar oscillations is promising. Our simulations are less dissipative and require less luminosity boosting than comparable spectral simulations. The fully-compressible approach allows the coupling of gravity and pressure waves to be studied too.

preprint2020arXiv

Long-term evolution of a magnetic massive merger product

About 10% of stars more massive than ${\approx}\,1.5\,\mathrm{M}_\odot$ have strong, large-scale surface magnetic fields and are being discussed as progenitors of highly-magnetic white dwarfs and magnetars. The origin of these fields remains uncertain. Recent 3D magnetohydrodynamical simulations have shown that strong magnetic fields can be generated in the merger of two massive stars. Here, we follow the long-term evolution of such a 3D merger product in a 1D stellar evolution code. During a thermal relaxation phase after the coalescence, the merger product reaches critical surface rotation, sheds mass and then spins down primarily because of internal mass readjustments. The spin of the merger product after thermal relaxation is mainly set by the co-evolution of the star-torus structure left after coalescence. This evolution is still uncertain, so we also consider magnetic braking and other angular-momentum-gain and -loss mechanisms that may influence the final spin of the merged star. Because of core compression and mixing of carbon and nitrogen in the merger, enhanced nuclear burning drives a transient convective core that greatly contributes to the rejuvenation of the star. Once the merger product relaxed back to the main sequence, it continues its evolution similar to that of a genuine single star of comparable mass. It is a slow rotator that matches the magnetic blue straggler $τ$ Sco. Our results show that merging is a promising mechanism to explain some magnetic massive stars and it may also be key to understand the origin of the strong magnetic fields of highly-magnetic white dwarfs and magnetars.

preprint2020arXiv

Nucleosynthesis imprints from different Type Ia Supernova explosion scenarios and implications for galactic chemical evolution

We analyze the nucleosynthesis yields of various Type Ia supernova explosion simulations including pure detonations in sub- Chandrasekhar mass white dwarfs, double detonations and pure helium detonations of sub-Chandrasekhar mass white dwarfs with an accreted helium envelope, a violent merger model of two white dwarfs and deflagrations as well as delayed detonations in Chandrasekhar mass white dwarfs. We focus on the iron peak elements Mn, Zn and Cu. To this end, we also briefly review the different burning regimes and production sites of these elements as well as the results of abundance measurements and several galactic chemical evolution studies. We find that super-solar values of [Mn/Fe] are not restricted to Chandrasekhar mass explosion models. Scenarios including a helium detonation can significantly contribute to the production of Mn, in particular the models proposed for calcium-rich transients. Although Type Ia supernovae are often not accounted for as production sites of Zn and Cu, our models involving helium shell detonations can produce these elements in super-solar ratios relative to Fe. Our results suggest a re-consideration of Type Ia supernova yields in galactic chemical evolution models. A detailed comparison with observations can provide new insight into the progenitor and explosion channels of these events.

preprint2020arXiv

Photometric detection of internal gravity waves in upper main-sequence stars. II. Combined TESS photometry and high-resolution spectroscopy

Context. Massive stars are predicted to excite internal gravity waves (IGWs) by turbulent core convection and from turbulent pressure fluctuations in their near-surface layers. These IGWs are extremely efficient at transporting angular momentum and chemical species within stellar interiors, but they remain largely unconstrained observationally. Aims. We aim to characterise the photometric detection of IGWs across a large number of O and early-B stars in the Hertzsprung-Russell diagram, and explain the ubiquitous detection of stochastic variability in the photospheres of massive stars. Methods. We combined high-precision time-series photometry from the NASA Transiting Exoplanet Survey Satellite with high-resolution ground-based spectroscopy of 70 stars with spectral types O and B to probe the relationship between the photometric signatures of IGWs and parameters such as spectroscopic mass, luminosity, and macroturbulence. Results. A relationship is found between the location of a star in the spectroscopic Hertzsprung-Russell diagram and the amplitudes and frequencies of stochastic photometric variability in the light curves of massive stars. Furthermore, the properties of the stochastic variability are statistically correlated with macroturbulent velocity broadening in the spectral lines of massive stars. Conclusions. The common ensemble morphology for the stochastic low-frequency variability detected in space photometry and its relationship to macroturbulence is strong evidence for IGWs in massive stars, since these types of waves are unique in providing the dominant tangential velocity field required to explain the observed spectroscopy.

preprint2020arXiv

White dwarf deflagrations for Type Iax supernovae: Polarisation signatures from the explosion and companion interaction

Growing evidence suggests that Type Iax supernovae might be the result of thermonuclear deflagrations of Chandrasekhar-mass white dwarfs in binary systems. We carry out Monte Carlo radiative transfer simulations and predict spectropolarimetric features originating from the supernova explosion and subsequent ejecta interaction with the companion star. Specifically, we calculate viewing-angle dependent flux and polarisation spectra for a 3D model simulating the deflagration of a Chandrasekhar-mass white dwarf and, for a second model, simulating the ejecta interaction with a main-sequence star. We find that the intrinsic signal is weakly polarised and only mildly viewing-angle dependent, owing to the overall spherical symmetry of the explosion and the depolarising contribution of iron-group elements dominating the ejecta composition. The interaction with the companion star carves out a cavity in the ejecta and produces a detectable, but modest signal that is significant only at relatively blue wavelengths ($\lesssim$ 5000 $\unicode{x212B}$). In particular, increasingly fainter and redder spectra are predicted for observer orientations further from the cavity, while a modest polarisation signal $P\sim0.2$ per cent is found at blue wavelengths for orientations 30$^\circ$ and 45$^\circ$ away from the cavity. We find a reasonable agreement between the interaction model viewed from these orientations and spectropolarimetric data of SN 2005hk and interpret the maximum-light polarisation signal seen at blue wavelengths for this event as a possible signature of the ejecta-companion interaction. We encourage further polarimetric observations of SNe Iax to test whether our results can be extended and generalised to the whole SN Iax class.

preprint2019arXiv

Discovery of an Exceptionally Strong $β$-Decay Transition of $^{20}$F and Implications for the Fate of Intermediate-Mass Stars

A significant fraction of stars between 7-11 solar masses are thought to become supernovae, but the explosion mechanism is unclear. The answer depends critically on the rate of electron capture on $^{20}$Ne in the degenerate oxygen-neon stellar core. However, due to the unknown strength of the transition between the ground states of $^{20}$Ne and $^{20}$F, it has not previously been possible to fully constrain the rate. By measuring the transition, we have established that its strength is exceptionally large and enhances the capture rate by several orders of magnitude. This has a decisive impact on the evolution of the core, increasing the likelihood that the star is (partially) disrupted by a thermonuclear explosion rather than collapsing to form a neutron star. Importantly, our measurement resolves the last remaining nuclear physics uncertainty in the final evolution of degenerate oxygen-neon stellar cores, allowing future studies to address the critical role of convection, which at present is poorly understood.

preprint2019arXiv

Monte Carlo radiative transfer for the nebular phase of Type Ia supernovae

We extend the range of validity of the ARTIS 3D radiative transfer code up to hundreds of days after explosion, when Type Ia supernovae are in their nebular phase. To achieve this, we add a non-local thermodynamic equilibrium (non-LTE) population and ionisation solver, a new multi-frequency radiation field model, and a new atomic dataset with forbidden transitions. We treat collisions with non-thermal leptons resulting from nuclear decays to account for their contribution to excitation, ionisation, and heating. We validate our method with a variety of tests including comparing our synthetic nebular spectra for the well-known one-dimensional W7 model with the results of other studies. As an illustrative application of the code, we present synthetic nebular spectra for the detonation of a sub-Chandrasekhar white dwarf in which the possible effects of gravitational settling of Ne22 prior to explosion have been explored. Specifically, we compare synthetic nebular spectra for a 1.06 M$_\odot$ white dwarf model obtained when 5.5 Gyr of very-efficient settling is assumed to a similar model without settling. We find that this degree of Ne22 settling has only a modest effect on the resulting nebular spectra due to increased Ni58 abundance. Due to the high ionisation in sub-Chandrasekhar models, the nebular [Ni II] emission remains negligible, while the [Ni III] line strengths are increased and the overall ionisation balance is slightly lowered in the model with Ne22 settling. In common with previous studies of sub-Chandrasekhar models at nebular epochs, these models overproduce [Fe III] emission relative to [Fe II] in comparison to observations of normal Type Ia supernovae.

preprint2018arXiv

The interaction of core-collapse supernova ejecta with a stellar companion

The progenitors of many core-collapse supernovae (CCSNe) are expected to be in binary systems. By performing a series of three-dimensional hydrodynamical simulations, we investigate how CCSN explosions affect their binary companion. We find that the amount of removed stellar mass, the resulting impact velocity, and the chemical contamination of the companion that results from the impact of the SN ejecta, strongly increases with decreasing binary separation and increasing explosion energy. Also, it is foud that the impact effects of CCSN ejecta on the structure of main-sequence (MS) companions, and thus their long term post-explosion evolution, is in general not be dramatic.