Source author record

Ondrej Pejcha

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

13works
4topics
4close collaborators

Actions

Connect this record

Log in to claim

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

13 published item(s)

preprint2022arXiv

The luminous red nova AT 2018bwo in NGC 45 and its binary yellow supergiant progenitor

Luminous Red Novae (LRNe) are astrophysical transients associated with the partial ejection of a binary system's common envelope (CE) shortly before its merger. Here we present the results of our photometric and spectroscopic follow-up campaign of AT2018bwo (DLT18x), a LRN discovered in NGC45, and investigate its progenitor system using binary stellar-evolution models. The transient reached a peak magnitude of $M_r=-10.97\pm0.11$ and maintained this brightness during its optical plateau of $t_p = 41\pm5$days. During this phase, it showed a rather stable photospheric temperature of ~3300K and a luminosity of ~$10^{40}$erg/s. The photosphere of AT2018bwo at early times appeared larger and cooler than other similar LRNe, likely due to an extended mass-loss episode before the merger. Towards the end of the plateau, optical spectra showed a reddened continuum with strong molecular absorption bands. The reprocessed emission by the cooling dust was also detected in the mid-infrared bands ~1.5 years after the outburst. Archival Spitzer and Hubble Space Telescope data taken 10-14 years before the transient event suggest a progenitor star with $T_{prog}\sim 6500$K, $R_{prog}\sim 100R_{\odot}$ and $L_{prog}\sim 2\times10^4L_{\odot}$, and an upper limit for optically thin warm (1000 K) dust mass of $M_d<10^{-6}M_{\odot}$. Using stellar binary-evolution models, we determined the properties of binary systems consistent with the progenitor parameter space. For AT2018bwo, we infer a primary mass of 12-16 $M_{\odot}$, which is 9-45% larger than the ~11$M_{\odot}$ obtained using single-star evolution models. The system, consistent with a yellow-supergiant primary, was likely in a stable mass-transfer regime with -2.4<log ($\dot{M}/M_{\odot}$/yr)<-1.2 a decade before the main instability occurred. During the dynamical merger, the system would have ejected 0.15-0.5$M_{\odot}$ with a velocity of ~500 km/s.

preprint2016arXiv

Binary Stellar Mergers with Marginally-Bound Ejecta: Excretion Disks, Inflated Envelopes, Outflows, and their Luminous Transients

We study mass loss from the outer Lagrange point (L2) in binary stellar mergers and their luminous transients by means of radiative hydrodynamical simulations. Previously, we showed that for binary mass ratios 0.06 < q < 0.8, synchronous L2 mass loss results in a radiatively inefficient, dust-forming unbound equatorial outflow. A similar outflow exists irrespective of q if the ratio of the sound speed to the orbital speed at the injection point is sufficiently large, ε= c_T/v_orb > 0.15. By contrast, for cold L2 mass-loss (ε< 0.15) from binaries with q < 0.06 or q > 0.8, the equatorial outflow instead remains marginally-bound and falls back to the binary over tens to hundreds of binary orbits, where it experiences additional tidal torqueing and shocking. As the bound gas becomes virialized with the binary, the luminosity of the system increases slowly at approximately constant photosphere radius, causing the temperature to rise. Subsequent evolution depends on the efficiency of radiative cooling. If the bound atmosphere is able to cool efficiently, as quantified by radiative diffusion time being shorter than the advection time (t_diff/t_adv < 1), then the virialized gas collapses to an excretion disk, while for t_diff/t_adv > 1 an isotropic wind is formed. Between these two extremes, an inflated envelope transports the heat generated near the binary to the surface by meridional flows. In all cases, the radiated luminosity reaches a fraction ~0.01 to 0.1 of Mdot v_orb^2/2, where Mdot is the mass outflow rate. We discuss the implications of our results for transients in the luminosity gap between classical novae and supernovae, such as V1309 Sco and V838 Mon.

preprint2015arXiv

Cool and Luminous Transients from Mass-Losing Binary Stars

We study transients produced by equatorial disk-like outflows from catastrophically mass-losing binary stars with an asymptotic velocity and energy deposition rate near the inner edge which are proportional to the binary escape velocity v_esc. As a test case, we present the first smoothed-particle radiation-hydrodynamics calculations of the mass loss from the outer Lagrange point with realistic equation of state and opacities. The resulting spiral stream becomes unbound for binary mass ratios 0.06 < q < 0.8. For synchronous binaries with non-degenerate components, the spiral-stream arms merge at a radius of ~10a, where a is the binary semi-major axis, and the accompanying shock thermalizes about 10% of the kinetic power of the outflow. The mass-losing binary outflows produce luminosities reaching up to ~10^6 L_Sun and effective temperatures spanning 500 < T_eff < 6000 K, which is compatible with many of the class of recently-discovered red transients such as V838 Mon and V1309 Sco. Dust readily forms in the outflow, potentially in a catastrophic global cooling transition. The appearance of the transient is viewing angle-dependent due to vastly different optical depths parallel and perpendicular to the binary plane. We predict a correlation between the peak luminosity and the outflow velocity, which is roughly obeyed by the known red transients. Outflows from mass-losing binaries can produce luminous (10^5 L_Sun) and cool (T_eff < 1500 K) transients lasting a year or longer, as has potentially been detected by Spitzer surveys of nearby galaxies.

preprint2015arXiv

On The Intrinsic Diversity of Type II-Plateau Supernovae

Hydrogen-rich Type II-Plateau supernovae exhibit correlations between the plateau luminosity $L_{\rm pl}$, the nickel mass $M_{\rm Ni}$, the explosion energy $E_{\rm exp}$, and the ejecta mass $M_{\rm ej}$. Using our global, self-consistent, multi-band model of nearby well-observed supernovae, we find that the covariances of these quantities are strong and that the confidence ellipsoids are oriented in the direction of the correlations, which reduces their significance. By proper treatment of the covariance matrix of the model, we discover a significant intrinsic width to the correlations between $L_{\rm pl}$, $E_{\rm exp}$, and $M_{\rm Ni}$, where the uncertainties due to the distance and the extinction dominate. For fixed $E_{\rm exp}$, the spread in $M_{\rm Ni}$ is about 0.25 dex, which we attribute to the differences in the progenitor internal structure. We argue that the effects of incomplete $γ$-ray trapping are not important in our sample. Similarly, the physics of the Type II-Plateau supernova light curves leads to inherently degenerate estimates of $E_{\rm exp}$ and $M_{\rm ej}$, which makes their observed correlation weak. Ignoring the covariances of supernova parameters or the intrinsic width of the correlations causes significant biases in the slopes of the fitted relations. Our results imply that Type II-Plateau supernova explosions are not described by a single physical parameter or a simple one-dimensional trajectory through the parameter space, but instead reflect the diversity of the core and surface properties of their progenitors. We discuss the implications for the physics of the explosion mechanism and possible future observational constraints.

preprint2015arXiv

The Landscape of the Neutrino Mechanism of Core-Collapse Supernovae: Neutron Star and Black Hole Mass Functions, Explosion Energies and Nickel Yields

If the neutrino luminosity from the proto-neutron star formed during a massive star core collapse exceeds a critical threshold, a supernova (SN) results. Using spherical quasi-static evolutionary sequences for hundreds of progenitors over a range of metallicities, we study how the explosion threshold maps onto observables, including the fraction of successful explosions, the neutron star (NS) and black hole (BH) mass functions, the explosion energies (E_SN) and nickel yields (M_Ni), and their mutual correlations. Successful explosions are intertwined with failures in a complex pattern that is not simply related to initial progenitor mass or compactness. We predict that progenitors with initial masses of 15 +/- 1, 19 +/- 1, and 21-26 M_Sun are most likely to form BHs, that the BH formation probability is non-zero at solar-metallicity and increases significantly at low metallicity, and that low luminosity, low Ni-yield SNe come from progenitors close to success/failure interfaces. We qualitatively reproduce the observed E_SN-M_Ni correlation, we predict a correlation between the mean and width of the NS mass and E_SN distributions, and that the means of the NS and BH mass distributions are correlated. We show that the observed mean NS mass of ~1.33 M_Sun implies that the successful explosion fraction is higher than 0.35. Overall, we show that the neutrino mechanism can in principle explain the observed properties of SNe and their compact objects. We argue that the rugged landscape of progenitors and outcomes mandates that SN theory should focus on reproducing the wide ranging distributions of observed SN properties.

preprint2014arXiv

A Global Model of The Light Curves and Expansion Velocities of Type II-Plateau Supernovae

We present a new self-consistent and versatile method that derives photospheric radius and temperature variations of Type II-Plateau supernovae based on their expansion velocities and photometric measurements. We apply the method to a sample of 26 well-observed, nearby supernovae with published light curves and velocities. We simultaneously fit ~230 velocity and ~6800 magnitude measurements distributed over 21 photometric passbands spanning wavelengths from 0.19 to 2.2 microns. The light curve differences among the Type II-Plateau supernovae are well-modeled by assuming different rates of photospheric radius expansion, which we explain as different density profiles of the ejecta and we argue that steeper density profiles result in flatter plateaus, if everything else remains unchanged. The steep luminosity decline of Type II-Linear supernovae is due to fast evolution of the photospheric temperature, which we verify with a successful fit of SN1980K. Eliminating the need for theoretical supernova atmosphere models, we obtain self-consistent relative distances, reddenings, and nickel masses fully accounting for all internal model uncertainties and covariances. We use our global fit to estimate the time evolution of any missing band tailored specifically for each supernova and we construct spectral energy distributions and bolometric light curves. We produce bolometric corrections for all filter combinations in our sample. We compare our model to the theoretical dilution factors and find good agreement for the B and V filters. Our results differ from the theory when the I, J, H, or K bands are included. We investigate the reddening law towards our supernovae and find reasonable agreement with standard R_V ~ 3.1 reddening law in UBVRI bands. Results for other bands are inconclusive. We make our fitting code publicly available.

preprint2014arXiv

Burying a Binary: Dynamical Mass Loss and a Continuous Optically-Thick Outflow Explain the Candidate Stellar Merger V1309 Scorpii

V1309 Sco was proposed to be a stellar merger and a common envelope transient based on the pre-outburst light curve of a contact eclipsing binary with a rapidly decaying orbital period. Using published data, I show that the period decay timescale P/Pdot of V1309 Sco decreased from ~1000 to ~170 years in less than about 6 years, which implies a very high value of second period derivative. I argue that V1309 Sco experienced an onset of dynamical mass loss through the outer Lagrange point, which eventually obscured the binary. The photosphere of the resulting continuous optically-thick outflow expands as the mass-loss rate increases, explaining the ~200 day rise to optical maximum. The model yields the mass-loss rate of the binary star as a function of time and fits the observed light curve remarkably well. It is also possible to observationally constrain the properties of the surface layers undergoing the dynamical mass loss. V1309 Sco is thus a prototype of a new class of stellar transients distinguished by a slow rise to optical maximum that are driven by dynamical mass loss from a binary. I discuss implications of these findings for stellar transients and other suggested common envelope events.

preprint2013arXiv

Greatly Enhanced Eccentricity Oscillations in Quadruple Systems Composed of Two Binaries: Implications for Stars, Planets, and Transients

We study the orbital evolution of hierarchical quadruple systems composed of two binaries on a long mutual orbit, where each binary acts as a Kozai-Lidov (KL) perturber on the other. We find that the coupling between the two binaries qualitatively changes the behavior of their KL cycles. The binaries can experience coherent eccentricity oscillations as well as excursions to very high eccentricity that occur over a much larger fraction of the parameter space than in triple systems. For a ratio of outer to inner semi-major axes of 10 to 20, about 30 to 50% of equal-mass quadruples reach eccentricity 1-e < 10^{-3} in one of the binaries. This is about 4 to 12 times more than for triples with equivalent parameters. Orbital "flips" and collisions without previous tidal interaction are similarly enhanced in quadruples relative to triples. We argue that the frequency of evolutionary paths influenced by KL cycles is comparable in the triple and quadruple populations even though field quadruples are a factor of ~5 less frequent than triples. Additionally, quadruples might be a non-negligible source of triples and provide fundamentally new evolutionary outcomes involving close binaries, mergers, collisions, and associated transients, which occur without any fine tuning of parameters. Finally, we study the perturbations to a planetary orbit due to a distant binary and we find that the fraction of orbital flips is a factor of 3 to 4 higher than for the corresponding triple system given our fiducial parameters with implications for hot Jupiters and star-planet collisions.

preprint2012arXiv

A Global Physical Model for Cepheids

We perform a global fit to ~5,000 radial velocity and ~177,000 magnitude measurements in 29 photometric bands covering 0.3 to 8.0 microns distributed among 287 Galactic, LMC, and SMC Cepheids with P > 10 days. We assume that the Cepheid light curves and radial velocities are fully characterized by distance, reddening, and time-dependent radius and temperature variations. We construct phase curves of radius and temperature for periods between 10 and 100 days, which yield light curve templates for all our photometric bands and can be easily generalized to any additional band. With only 4 to 6 parameters per Cepheid, depending on the existence of velocity data and the amount of freedom in the distance, the models have typical rms light and velocity curve residuals of 0.05 mag and 3.5 km/s. The model derives the mean Cepheid spectral energy distribution and its derivative with respect to temperature, which deviate from a black body in agreement with metal-line and molecular opacity effects. We determine a mean reddening law towards the Cepheids in our sample, which is not consistent with standard assumptions in either the optical or near-IR. Based on stellar atmosphere models we predict the biases in distance, reddening, and temperature determinations due to the metallicity and we quantify the metallicity signature expected for our fit residuals. The observed residuals as a function of wavelength show clear differences between the individual galaxies, which are compatible with these predictions. In particular, we find that metal-poor Cepheids appear hotter. Finally, we provide a framework for optimally selecting filters that yield the smallest overall errors in Cepheid parameter determination, or filter combinations for suppressing or enhancing the metallicity effects on distance determinations. We make our templates publicly available.

preprint2012arXiv

Effect of Collective Neutrino Oscillations on the Neutrino Mechanism of Core-Collapse Supernovae

In the seconds after collapse of a massive star, the newborn proto-neutron star (PNS) radiates neutrinos of all flavors. The absorption of electron-type neutrinos below the radius of the stalled shockwave may drive explosions (the "neutrino mechanism"). Because the heating rate is proportional to the square of neutrino energy, flavor conversion of mu and tau neutrinos to electron-type neutrinos via collective neutrino oscillations (CnuO) may in principle increase the heating rate and drive explosions. In order to assess the potential importance of CnuO for the shock revival, we solve the steady-state boundary value problem of spherically-symmetric accretion between the PNS surface (r_nu) and the shock (r_S), including a scheme for flavor conversion via CnuO. For a given r_nu, PNS mass (M), accretion rate (Mdot), and assumed values of the neutrino energies from the PNS, we calculate the critical neutrino luminosity above which accretion is impossible and explosion results. We show that CnuO can decrease the critical luminosity by a factor of at most ~1.5, but only if the flavor conversion is fully completed inside r_S and if there is no matter suppression. The magnitude of the effect depends on the model parameters (M, Mdot, and r_nu) through the shock radius and the physical scale for flavor conversion. We quantify these dependencies and find that CnuO could lower the critical luminosity only for small M and Mdot, and large r_nu. However, for these parameter values CnuO are suppressed due to matter effects. By quantifying the importance of CnuO and matter suppression at the critical neutrino luminosity for explosion, we show in agreement with previous studies that CnuO are unlikely to affect the neutrino mechanism of core-collapse supernovae significantly.

preprint2012arXiv

The observed neutron star mass distribution as a probe of the supernova explosion mechanism

The observed distribution of neutron star (NS) masses reflects the physics of core-collapse supernova explosions and the structure of the massive stars that produce them at the end of their evolution. We present a Bayesian analysis that directly compares the NS mass distribution observed in double NS systems to theoretical models of NS formation. We find that models with standard binary mass ratio distributions are strongly preferred over independently picking the masses from the initial mass function, although the strength of the inference depends on whether current assumptions for identifying the remnants of the primary and secondary stars are correct. Second, NS formation models with no mass fallback are favored because they reduce the dispersion in NS masses. The double NS system masses thus directly point to the mass coordinate where the supernova explosion was initiated, making them an excellent probe of the supernova explosion mechanism. If we assume no fallback and simply vary the mass coordinate separating the remnant and the supernova ejecta, we find that for solar metallicity stars the explosion most likely develops at the edge of the iron core at a specific entropy of about 2.8 k_B. The primary limitations of our study are the poor knowledge of the supernova explosion mechanism and the lack of broad range of SN model explosions of LMC to solar metallicity.

preprint2011arXiv

The Physics of the Neutrino Mechanism of Core-Collapse Supernovae

(Abridged) Neutrino heating may drive core-collapse supernova explosions. Although it is known that the stalled accretion shock turns into explosion when the neutrino luminosity from the collapsed core exceeds a critical value (L_crit) (the "neutrino mechanism"), the physics of L_crit, as well as its dependence on the properties of the proto-neutron star (PNS) and changes to the microphysics has never been systematically explored. We solve the one-dimensional steady-state accretion problem between the PNS surface and the accretion shock. We quantify the deep connection between the solution space of steady-state accretion flows with bounding shocks and the neutrino mechanism. We show that there is a maximum, critical sound speed above which it is impossible to maintain accretion with a standoff shock, because the shock jump conditions cannot be satisfied. The physics of this critical sound speed is general and does not depend on a specific heating mechanism. For the simple model of pressure-less free-fall onto a shock bounding an isothermal accretion flow with sound speed c_T, we show that if c_T^2/v_escape^2 > 3/16 explosion results. We generalize this result to the more complete supernova problem, showing explicitly that the same physics determines L_crit. We find that the critical condition for explosion can be written as c_S^2/v_escape^2 = 0.19, where c_S is the adiabatic sound speed. This "antesonic" condition describes L_crit over a broad range in accretion rate and microphysics. We show that the addition of the accretion luminosity (L_acc) reduces L_crit non-trivially. As in previous work, we find that L_crit is always significantly higher than the maximum possible value of L_acc. Finally, we provide evidence that the reduction in L_crit seen in recent multi-dimensional simulations results from a reduction in the efficiency of cooling, rather than an increase in the heating rate.

preprint2010arXiv

An R- and I-Band Photometric Variability Survey of the Cygnus OB2 Association

We present a catalog of photometrically variable stars discovered within two 21'.3 X 21'.3 fields centered on the Cygnus OB2 association. There have hitherto been no deep optical variability studies of Cyg OB2 despite it being replete with early-type massive stars, perhaps due to the high and variable extinction (up to A_V ~ 20) that permeates much of the region. Here we provide results of the first variability study with this combination of spatial coverage (~ 0.5 deg) and photometric depth (R ~ 21 mag). We find 121 stars to be variable in both R- and I-band, 116 of them newly discovered. Of the 121 variables, we identify 27 eclipsing binaries (EBs) and eclipsing binary candidates, 20 potential Herbig Ae/Be stars, and 52 pulsating variables. Confirming both the status and the cluster membership of the Herbig Ae/Be stars would address the uncertainty regarding the age and star formation history of Cyg OB2. We match our catalog to known variables and binaries in the region, 2MASS near-IR (NIR) data, and Chandra X-ray observations to find counterparts to new variables in other wavelengths.