-
Observability of flashes from ejecta crashes in aspherical supernovae, with application to SN 2008D
Authors:
Benjamin Scully,
Christopher D. Matzner,
Almog Yalinewich
Abstract:
A new class of transient, which has been hypothesized to accompany the explosion of an aspherical compact supernova, would arise when streams of ejecta collide outside the star. However, conditions that favour the prompt release of radiation from the collision, such as a diffuse stellar envelope, disfavour the creation of non-radial ejecta in the first place. To determine whether the collision can…
▽ More
A new class of transient, which has been hypothesized to accompany the explosion of an aspherical compact supernova, would arise when streams of ejecta collide outside the star. However, conditions that favour the prompt release of radiation from the collision, such as a diffuse stellar envelope, disfavour the creation of non-radial ejecta in the first place. To determine whether the collision can both occur and be visible, we simulate aspherical explosions using the HUJI-RICH moving-mesh hydrodynamics code and analyze them in terms of diffusion measures defined for individual fluid elements. While our simulations are highly idealized, they connect to realistic explosions via a single dimensionless parameter. Defining two measures of the importance of diffusivity (two versions of the inverse Péclet number), we find that one varies in a way that indicates colliding ejecta can release a photon flash, while the other does not. Examining the x-ray transient XT 080109 associated with supernova SN 2008D, we find that its fluence and duration are consistent with the properties of an ejecta collision in the aspherical model that is most likely to emit a flash. Our results give tentative evidence for the possibility of collision-induced flashes for a narrow and radius-dependent range of asphericity, and motivate future radiation hydrodynamics simulations.
△ Less
Submitted 17 August, 2023; v1 submitted 28 July, 2023;
originally announced July 2023.
-
The Moving Mirror model for Fast Radio Bursts
Authors:
Almog Yalinewich,
Ue-Li Pen
Abstract:
Recent observations of coherent radiation from the Crab pulsar (Bij et al 2021) suggest the emission is driven by an ultra - relativistic ($γ\sim 10^4$), cold plasma flow. A relativistically expanding plasma shell can compress the ambient magnetic field, like a moving mirror, and thus produce coherent radiation whose wavelength is shorter than that of the ambient medium by $γ^2$. This mechanism ha…
▽ More
Recent observations of coherent radiation from the Crab pulsar (Bij et al 2021) suggest the emission is driven by an ultra - relativistic ($γ\sim 10^4$), cold plasma flow. A relativistically expanding plasma shell can compress the ambient magnetic field, like a moving mirror, and thus produce coherent radiation whose wavelength is shorter than that of the ambient medium by $γ^2$. This mechanism has been studied in the past by Colgate and Noerdelinger (1971), in the context of radio loud supernova explosions. In this work we propose that a similar mechanism drives the coherent emission in fast radio bursts. The high Lorenz factors dramatically lower the implied energy and magnetic field requirements, allowing the spin down energy of regular (or even recycled), fast spinning pulsars, rather than slow spinning magnetars, to explain FRBs. We show that this model can explain the frequency and the time evolution of observed FRBs, as well as their duration, energetics and absence of panchromatic counterparts. We also predict that the peak frequency of sub pulses decline with observation time as $ω_{\rm obs} \propto t_{\rm obs}^{-1/2}$. Unfortunately, with current capabilities it is not possible to constrain the shape of the curve $ω_{\rm obs} \left(t_{\rm obs} \right)$. Finally, we find that a variation of this model can explain weaker radio transients, such as the one observed from a galactic magnetar. In this variant, the shock wave produces low frequency photons which are then Compton scattered to the GHz range.
△ Less
Submitted 25 April, 2022;
originally announced April 2022.
-
Plug Disintegration in GRB Jet Eruption
Authors:
Almog Yalinewich,
Paz Beniamini
Abstract:
In this work we consider the eruption of a tenuous relativistic hydrodynamic jet from a dense baryonic envelope. As the jet moves out and away, it carries along and continues to accelerate a layer of baryonic material which we refer to as the plug. We solve the relativistic equations of motion for the trajectory of the plug, and verify it using a relativistic hydrodynamic simulation. We show that…
▽ More
In this work we consider the eruption of a tenuous relativistic hydrodynamic jet from a dense baryonic envelope. As the jet moves out and away, it carries along and continues to accelerate a layer of baryonic material which we refer to as the plug. We solve the relativistic equations of motion for the trajectory of the plug, and verify it using a relativistic hydrodynamic simulation. We show that under these conditions, the plug breaks up at a radius larger by a factor of a few from the radius of the envelope, due to the onset of the Rayleigh Taylor instability. After breakup the jet continues to accelerate to higher Lorentz factors while the plug fragments maintain a moderate Lorentz factor. The presence of slower moving ejecta can explain late time features of GRBs such as X ray flares without recourse to a long lived engine.
△ Less
Submitted 26 January, 2022;
originally announced January 2022.
-
The Propagation of Strong Shocks into Planetary and Stellar Atmospheres
Authors:
Almog Yalinewich,
Andrey Remorov
Abstract:
In this work we present a mathematical model for the propagation of the shock waves that occur in graded density profiles. These waves can occur in a wide range of astrophysical events, such as collisions in planetary and stellar atmospheres, common envelope explosions and peculiar type Ia supernovae. The behaviour of the shock wave and its evolution can be modelled using type II self similar solu…
▽ More
In this work we present a mathematical model for the propagation of the shock waves that occur in graded density profiles. These waves can occur in a wide range of astrophysical events, such as collisions in planetary and stellar atmospheres, common envelope explosions and peculiar type Ia supernovae. The behaviour of the shock wave and its evolution can be modelled using type II self similar solutions. In such solutions the evolution of the shock wave is determined by boundary conditions at the shock front and a singular point in the shocked region. We show how the evolution can be determined for different equations of state and density profiles, and compare these results to numerical simulations. These findings are also applied to a variety of astrophysical phenomena to further test their validity.
△ Less
Submitted 4 July, 2021;
originally announced July 2021.
-
Crater Morphology of Primordial Black Hole Impacts
Authors:
Almog Yalinewich,
Matthew E. Caplan
Abstract:
In this work we propose a novel campaign for constraining relativistically compact MACHO dark matter, such as primordial black holes (PBHs), using the moon as a detector. PBHs of about $10^{19} \textrm{ g}$ to $10^{22} \textrm{ g}$ may be sufficiently abundant to have collided with the moon in the history of the solar system. We show that the crater profiles of a PBH collision differ from traditio…
▽ More
In this work we propose a novel campaign for constraining relativistically compact MACHO dark matter, such as primordial black holes (PBHs), using the moon as a detector. PBHs of about $10^{19} \textrm{ g}$ to $10^{22} \textrm{ g}$ may be sufficiently abundant to have collided with the moon in the history of the solar system. We show that the crater profiles of a PBH collision differ from traditional impactors and may be detectable in high resolution lunar surface scans now available. Any candidates may serve as sites for in situ measurements to identify high pressure phases of matter which may have formed near the PBH during the encounter. While we primarily consider PBH dark matter, the discussion generalises to the entire family of MACHO candidates with relativistic compactness. Moreover, we focus on the Moon since it has been studied well, but the same principles can be applied to other rocky bodies in our solar system without an atmosphere.
△ Less
Submitted 31 March, 2021;
originally announced April 2021.
-
Self Similar Adiabatic Strong Explosion in a Medium Gravitationally Free Falling to a Point Mass
Authors:
Almog Yalinewich
Abstract:
We develop a generalisation to the classical Sedov Taylor explosion where the medium free falls to a point mass at the centre of the explosion. To verify our analytic results, we compare them to a suite of numerical simulations. We find that there exists a critical energy below which, instead of propagating outward the shock stalls and collapses under gravity. Furthermore, we find that the value o…
▽ More
We develop a generalisation to the classical Sedov Taylor explosion where the medium free falls to a point mass at the centre of the explosion. To verify our analytic results, we compare them to a suite of numerical simulations. We find that there exists a critical energy below which, instead of propagating outward the shock stalls and collapses under gravity. Furthermore, we find that the value of the critical energy threshold decreases when the adiabatic index increases and material is more evenly distributed within the shocked region. We apply this model to the problem of a shock bounce in core collapse supernova, in which the proto neutron star serves as the point mass. The relation between the threshold energy and the distribution of mass in the shock might help explain how turbulence prevents shock stalling and recession in a core collapse supernova explosion.
△ Less
Submitted 25 August, 2020;
originally announced August 2020.
-
Self similar Shocks in Atmospheric Mass Loss due to Planetary Collisions
Authors:
Almog Yalinewich,
Andrey Remorov
Abstract:
We present a mathematical model for the propagation of the shock waves that occur during planetary collisions. Such collisions are thought to occur during the formation of terrestrial planets, and they have the potential to erode the planet's atmosphere. We show that under certain assumptions, this evolution of the shock wave can be determined using the method of self similar solutions. This self…
▽ More
We present a mathematical model for the propagation of the shock waves that occur during planetary collisions. Such collisions are thought to occur during the formation of terrestrial planets, and they have the potential to erode the planet's atmosphere. We show that under certain assumptions, this evolution of the shock wave can be determined using the method of self similar solutions. This self similar solution is of type II, which means that it only applies to a finite region behind the shock front. This region is bounded by the shock front and the sonic point. Energy and matter continuously flow through the sonic point, so that energy in the self similar region is not conserved, as is the case for type I solutions. Instead, the evolution of the shock wave is determined by boundary conditions at the shock front and at the sonic point. We show how the evolution can be determined for different equations of state, allowing these results to be readily used to calculate the atmospheric mass loss from planetary cores made of different materials.
△ Less
Submitted 15 November, 2019;
originally announced November 2019.
-
Optical Transient from an Explosion Close to the Stellar Surface
Authors:
Almog Yalinewich,
Christopher D. Matzner
Abstract:
We study the hydrodynamic evolution of an explosion close to the stellar surface, and give predictions for the radiation from such an event. We show that such an event will give rise to a multi-wavelength transient. We apply this model to describe a precursor burst to the peculiar supernova iPTF14hls, which occurred in 1954, sixty year before the supernova. We propose that the new generation of op…
▽ More
We study the hydrodynamic evolution of an explosion close to the stellar surface, and give predictions for the radiation from such an event. We show that such an event will give rise to a multi-wavelength transient. We apply this model to describe a precursor burst to the peculiar supernova iPTF14hls, which occurred in 1954, sixty year before the supernova. We propose that the new generation of optical surveys might detect similar transients, and they can be used to identify supernova progenitors well before the explosion.
△ Less
Submitted 17 July, 2019;
originally announced July 2019.
-
Nekhoroshev Estimates for the Survival Time of Tightly Packed Planetary Systems
Authors:
Almog Yalinewich,
Cristobal Petrovich
Abstract:
$N…
▽ More
$N$-body simulations of non-resonant tightly-packed planetary systems have found that their survival time (i.e. time to first close encounter) grows exponentially with their interplanetary spacing and planetary masses. Although this result has important consequences for the assembly of planetary systems by giants collisions and their long-term evolution, this underlying exponential dependence is not understood from first principles, and previous attempts based on orbital diffusion have only yielded power-law scalings. We propose a different picture, where large deviations of the system from its initial conditions is due to few slowly developing high-order resonances. Thus, we show that the survival time of the system $T$ can be estimated using a heuristic motivated by Nekhoroshev's theorem, and obtain a formula for systems away from overlapping two-body mean-motion resonances as: $T/P=c_1 \frac{a}{Δa} \exp \left(c_2 \frac{Δa}{a} /μ^{1/4}\right)$, where $P$ is the average Keplerian period, $a$ is the average semi major axis, $Δa\ll a$ is the difference between the semi major axes of neighbouring planets, $μ$ is the planet to star mass ratio, and $c_1$ and $c_2$ are dimensionless constants. We show that this formula is in good agreement with numerical N-body experiments for $c_1=5 \cdot 10^{-4}$ and $c_2=8$.
△ Less
Submitted 16 January, 2020; v1 submitted 15 July, 2019;
originally announced July 2019.
-
Fast Radio Bursts from Terraformation
Authors:
Almog Yalinewich,
Mubdi Rahman,
Alysa Obertas,
Patrick C. Breysse
Abstract:
Fast radio bursts (FRBs) are, as the name implies, short and intense pulses of radiation at wavelengths of roughly one metre. FRBs have extremely high brightness temperatures, which points to a coherent source of radiation. The energy of a single burst ranges from $10^{36}$ to $10^{39}$ erg. At the high end of the energy range, FRBs have enough energy to unbind an earth-sized planet, and even at t…
▽ More
Fast radio bursts (FRBs) are, as the name implies, short and intense pulses of radiation at wavelengths of roughly one metre. FRBs have extremely high brightness temperatures, which points to a coherent source of radiation. The energy of a single burst ranges from $10^{36}$ to $10^{39}$ erg. At the high end of the energy range, FRBs have enough energy to unbind an earth-sized planet, and even at the low end, there is enough energy to vaporise and unbind the atmosphere and the oceans. We therefore propose that FRBs are signatures of an artificial terraformer, capable of eradicating life on another planet, or even destroy the planet entirely. The necessary energy can be harvested from Wolf-Rayet stars with a Dyson sphere ($\sim 10^{38}$ erg s$^{-1}$) , and the radiation can be readily produced by astrophysical masers. We refer to this mechanism as Volatile Amplification of a Destructive Emission of Radiation (VADER). We use the observational information to constrain the properties of the apparatus. We speculate that the non-repeating FRBs are low-energy pulses used to exterminate life on a single planet, but leaving it otherwise intact, and that the stronger repeating FRB is part of an effort to destroy multiple objects in the same solar system, perhaps as a preventative measure against panspermia. In this picture, the persistent synchrotron source associated with the first repeating FRB arises from the energy harvesting process. Finally we propose that Oumuamua might have resulted from a destruction of a planet in this manner.
△ Less
Submitted 28 March, 2019;
originally announced March 2019.
-
Radio Emission from the unbound Debris of Tidal Disruption Events
Authors:
Almog Yalinewich,
Elad Steinberg,
Tsvi Piran,
Julian H. Krolik
Abstract:
When a star gets too close to a supermassive black hole, it is torn apart by the tidal forces. Roughly half of the stellar mass becomes unbound and flies away at tremendous velocities - around $10^4$ km/s. In this work we explore the idea that the shock produced by the interaction of the unbound debris with the ambient medium gives rise to the synchrotron radio emission observed in several TDEs. W…
▽ More
When a star gets too close to a supermassive black hole, it is torn apart by the tidal forces. Roughly half of the stellar mass becomes unbound and flies away at tremendous velocities - around $10^4$ km/s. In this work we explore the idea that the shock produced by the interaction of the unbound debris with the ambient medium gives rise to the synchrotron radio emission observed in several TDEs. We use a moving mesh numerical simulation to study the evolution of the unbound debris and the bow shock around it. We find that as the periapse distance of the star decreases, the outflow becomes faster and wider. A tidal disruption event whose periapse distance is a factor of 7 smaller than the tidal radius can account for the radio emission observed in ASASSN-14li. This model also allows us to obtain a more accurate estimate for the gas density around the centre of the host galaxy of ASASSN-14li.
△ Less
Submitted 4 June, 2019; v1 submitted 6 March, 2019;
originally announced March 2019.
-
Statistical inference of the distance to ASKAP FRBs
Authors:
Dongzi Li,
Almog Yalinewich,
Patrick C. Breysse
Abstract:
The distances to fast radio bursts (FRBs) are crucial for understanding their underlying engine, and for their use as cosmological probes. In this paper, we provide three statistical estimates of the distance to ASKAP FRBs. First, we show that the number of events of similar luminosity in ASKAP does not scale as distance cubed, as one would expect, when directly using the observed dispersion measu…
▽ More
The distances to fast radio bursts (FRBs) are crucial for understanding their underlying engine, and for their use as cosmological probes. In this paper, we provide three statistical estimates of the distance to ASKAP FRBs. First, we show that the number of events of similar luminosity in ASKAP does not scale as distance cubed, as one would expect, when directly using the observed dispersion measure (DM) to infer distance. Second, by comparing the average DMs of FRBs observed with different instruments, we estimated the average redshift of ASKAP FRBs to be $z\sim 0.01$ using CHIME and ASKAP, and $z\lesssim0.07$ using Parkes and ASKAP. Both values are much smaller than the upper limit $z\sim0.3$ estimated directly from the DM. Third, we cross-correlate the locations of the ASKAP FRBs with existing large-area redshift surveys, and see a 3$σ$ correlation with the 2MASS Redshift Survey and a 5$σ$ correlation with the HI Parkes All Sky Survey at $z\sim0.007$. This corresponds well with the redshift of the most likely host galaxy of ASKAP FRB 171020, which is at $z=0.00867$. These arguments combined suggest an extremely nearby origin of ASKAP FRBs and a local environment with accumulated electrons that contribute a DM of several hundred pc/cm$^3$, which should be accounted for in theoretical models.
△ Less
Submitted 26 February, 2019;
originally announced February 2019.
-
Atmospheric Mass Loss from High Velocity Giant Impacts
Authors:
Almog Yalinewich,
Hilke E. Schlichting
Abstract:
Using moving mesh hydrodynamic simulations, we determine the shock propagation and resulting ground velocities for a planet hit by a high velocity impactor. We use our results to determine the atmospheric mass loss caused by the resulting ground motion due to the impact shock wave. We find that there are two distinct shock propagation regimes: In the limit in which the impactor is significantly sm…
▽ More
Using moving mesh hydrodynamic simulations, we determine the shock propagation and resulting ground velocities for a planet hit by a high velocity impactor. We use our results to determine the atmospheric mass loss caused by the resulting ground motion due to the impact shock wave. We find that there are two distinct shock propagation regimes: In the limit in which the impactor is significantly smaller than the target ($R_i<< R_t$), the solutions are self-similar and the shock velocity at a fixed point on the target scale as $m_i^{2/3}$, where $m_i$ is the mass of the impactor. In addition, the ground velocities follow a universal profile given by $v_g/v_i=(14.2x^2-25.3x+11.3)/(x^2-2.5x+1.9) +2\ln{R_i/R_t}$, where $x=\sin\left(θ/2\right)$, $θ$ is the latitude on the target measured from the impact site, and $v_g$ and $v_i$ are the ground velocity and impact velocity, respectively. In contrast, in the limit in which the impactor is comparable to the size of the target ($R_i \sim R_t$), we find that shock velocities decline with the mass of the impactor significantly more weakly than $m_i^{2/3}$. We use the resulting surface velocity profiles to calculate the atmospheric mass loss for a large range of impactor masses and impact velocities and apply them to the Kepler-36 system and the Moon forming impact. Finally, we present and generalise our results in terms of the $v_g/v_i$ and the impactor to target size ratio ($R_i/R_t$) such that they can easily be applied to other collision scenarios.
△ Less
Submitted 16 April, 2019; v1 submitted 28 November, 2018;
originally announced November 2018.
-
The Signature of a Windy Radio Supernova Progenitor in a Binary System
Authors:
Almog Yalinewich,
Simon Portegies Zwart
Abstract:
Type II supernova progenitors are expected to emit copious amounts of mass in a dense stellar wind prior to the explosion. When the progenitor is a member of a binary, the orbital motion modulates the density of this wind. When the progenitor explodes, the high-velocity ejecta collides with the modulated wind, which in turn produces a modulated radio signal. In this work we derive general analytic…
▽ More
Type II supernova progenitors are expected to emit copious amounts of mass in a dense stellar wind prior to the explosion. When the progenitor is a member of a binary, the orbital motion modulates the density of this wind. When the progenitor explodes, the high-velocity ejecta collides with the modulated wind, which in turn produces a modulated radio signal. In this work we derive general analytic relations between the parameters of the radio signal modulations and binary parameter in the limit of large member mass ratio. We use these relations to infer the semi major axis of SN1979c and a lower bound for the mass of the companion. We further constrain the analytic estimates by numerical simulations using the AMUSE framework. In these calculations we simulate the progenitor binary system including the wind and the gravitational effect of a companion star. The simulation output is compared to the observed radio signal in supernova SN1979C. We find that it must have been a binary with an orbital period of about 2000 year. If the exploding star evolved from a $\sim 18 M_{\odot}$ zero-age main-sequence at solar metalicity, we derive a companion mass of $5$ to $12 M_{\odot}$ in an orbit with an eccentricity lower than about 0.8.
△ Less
Submitted 29 January, 2019; v1 submitted 27 November, 2018;
originally announced November 2018.
-
Shock Breakouts from Tidal Disruption Events
Authors:
Almog Yalinewich,
James Guillochon,
Re'em Sari,
Abraham Loeb
Abstract:
Tidal disruption events of stars by supermassive black holes have so far been discovered months to years after the fact. In this paper we explore the short, faint and hard burst of radiation is emitted at maximum compression, as a result of shock breakout. The detection of this burst can be used to capture tidal disruption events in real time. We verify that shock breakout from main sequence stars…
▽ More
Tidal disruption events of stars by supermassive black holes have so far been discovered months to years after the fact. In this paper we explore the short, faint and hard burst of radiation is emitted at maximum compression, as a result of shock breakout. The detection of this burst can be used to capture tidal disruption events in real time. We verify that shock breakout from main sequence stars produces radiation in the X-ray range, but find that it is difficult to detect using all sky X-ray surveying telescopes. In the case of shock breakout from red giants, most of the radiation is emitted in the UV and visible range, which is significantly easier to detect. A similar burst of UV/optical radiation will also be emitted by stars puffed by tidal heating from a previous passage close to the central black hole. This radiation can be detected by surveys like ZTF and LSST. We calculate detection rates for different types of galactic nuclei. For the case of a very full loss cone we predict a detection rate of once per month with LSST, whereas for the case of a very empty loss cone we predict a rate of once per year with LSST. Evidence from a recent tidal disruption event, ASASSN-14li, seems to favour a very full loss cone, in which case LSST is expected to detect one such event every month.
△ Less
Submitted 6 September, 2018; v1 submitted 30 August, 2018;
originally announced August 2018.
-
Dark Passengers in Stellar Surveys
Authors:
Almog Yalinewich,
Paz Beniamini,
Kenta Hotokezaka,
Wei Zhu
Abstract:
We develop stellar population models to predict the number of binaries with a single luminous member is Gaia and Hipparcos. Our models yield dozens of detections of black hole - luminous companion binaries (BHLC) and hundreds to thousands of neutron star - luminous companion binaries (NSLC) with Gaia. Interestingly, our models also yield a single detection of BHLC binary with Hipparcos, and a few…
▽ More
We develop stellar population models to predict the number of binaries with a single luminous member is Gaia and Hipparcos. Our models yield dozens of detections of black hole - luminous companion binaries (BHLC) and hundreds to thousands of neutron star - luminous companion binaries (NSLC) with Gaia. Interestingly, our models also yield a single detection of BHLC binary with Hipparcos, and a few NSLC binaries. We also show how the statistical distribution of detected binaries with a single luminous companion can be used to constrain the formation process of neutron stars and black holes.
△ Less
Submitted 22 August, 2018; v1 submitted 2 July, 2018;
originally announced July 2018.
-
A Generalised Bondi Accretion Model for the Galactic Centre
Authors:
Almog Yalinewich,
Re'em Sari,
Aleksey Generozov,
Nicholas Chamberlain Stone,
Brian David Metzger
Abstract:
We develop an analytic, steady-state model for the gas environment in quiescent galactic nuclei. We assume that the mass is constantly supplied by a spherically symmetric distribution of wind emitting stars, and that gravity is solely due to a central supermassive black hole. We show that at some finite radius, where the Keplerian velocity is comparable to the wind velocity, the bulk velocity vani…
▽ More
We develop an analytic, steady-state model for the gas environment in quiescent galactic nuclei. We assume that the mass is constantly supplied by a spherically symmetric distribution of wind emitting stars, and that gravity is solely due to a central supermassive black hole. We show that at some finite radius, where the Keplerian velocity is comparable to the wind velocity, the bulk velocity vanishes. Matter generated below that radius will be accreted onto the black hole, while matter outside it will escape the system. Under certain conditions, the flow may become supersonic at both domains. We obtain radial profiles of the hydrodynamic variables and verify them using a time-dependent hydrodynamic simulation. We delineate the conditions under which radiative cooling can be neglected, and predict the luminosity and spectrum of the free-free X-ray emission from such a system. We discuss applications of our solution to our own Galactic Centre and other quiescent galactic nuclei.
△ Less
Submitted 12 May, 2018; v1 submitted 9 May, 2018;
originally announced May 2018.
-
Ignition of detonation in accreted helium
Authors:
Ami Glasner,
Eli Livne,
Elad Steinberg,
Almog Yalinewich,
James W. Truran
Abstract:
Sub-Chandrasekhar CO white dwarfs accreting helium have been considered as candidates for Type Ia supernova(SNIa) progenitors since the early 1980s (helium shell mass $> 0.1 M_\odot $). These models, once detonated did not fit the observed spectra and light curve of typical SNIa observations. New theoretical work examined detonations on much less massive ($< 0.05 M_\odot $) envelopes. They find st…
▽ More
Sub-Chandrasekhar CO white dwarfs accreting helium have been considered as candidates for Type Ia supernova(SNIa) progenitors since the early 1980s (helium shell mass $> 0.1 M_\odot $). These models, once detonated did not fit the observed spectra and light curve of typical SNIa observations. New theoretical work examined detonations on much less massive ($< 0.05 M_\odot $) envelopes. They find stable detonations that lead to light curves, spectra and abundances that compare relatively well with the observational data. The exact mechanism leading to the ignition of helium detonation is a key issue, since it is a mandatory first step for the whole scenario. As the flow of the accreted envelope is unstable to convection long before any hydrodynamic phenomena develops, a multidimensional approach is needed in order to study the ignition process. The complex convective reactive flow is challenging to any hydrodynamical solver. To the best of our knowledge, all previous 2D studies ignited the detonation artificially. We present here, for the first time, fully consistent results from two hydrodynamical 2D solvers that adopt two independent accurate schemes. For both solvers an effort was made to overcome the problematics raised by the finite resolution and numerical diffusion by the advective terms. Our best models lead to the ignition of a detonation in a convective cell. Our results are robust and the agreement between the two different numerical approaches is very good.
△ Less
Submitted 2 March, 2018;
originally announced March 2018.
-
Supernovae generated High Velocity Compact Clouds
Authors:
Almog Yalinewich,
Paz Beniamini
Abstract:
It has been proposed that an intermediate mass black hole can produce compact molecular clouds with large velocity dispersions. This model is called into question due to the discovery of several such short lived clouds, which suggests they are common, and hence it is unlikely that their formation involves an exotic object such as an intermediate mass black hole. In this paper we propose an alterna…
▽ More
It has been proposed that an intermediate mass black hole can produce compact molecular clouds with large velocity dispersions. This model is called into question due to the discovery of several such short lived clouds, which suggests they are common, and hence it is unlikely that their formation involves an exotic object such as an intermediate mass black hole. In this paper we propose an alternative model, where supernovae produce such clouds. We apply our model to a such a cloud called CO-0.40-0.22. A previous work pointed out the fact that the kinetic energy of this cloud is comparable to that from a supernova, but disqualified the supernova scenario due to lack of thermal emission. We explain the lack of thermal emission as a consequence of the large density of the medium. In such a medium, the supernova shock is radiation dominated, and hence the initial temperatures are lower by many orders of magnitude compared with a typical supernova. Finally, we discuss some of the consequences of having a significant population of intermediate mass black holes in the galactic centre.
△ Less
Submitted 17 September, 2017;
originally announced September 2017.
-
Evolution of Supernova Remnants near the Galactic Center
Authors:
Almog Yalinewich,
Tsvi Piran,
Re'em Sari
Abstract:
Supernovae near the galactic center evolve differently from regular galactic supernovae. This is mainly due to the environment into which the supernova remnants propagate. Instead of a static, uniform density medium, SNRs near the galactic center propagate into a wind swept environment with a velocity away from the galactic center, and a graded density profile. This causes these SNRs to be non - s…
▽ More
Supernovae near the galactic center evolve differently from regular galactic supernovae. This is mainly due to the environment into which the supernova remnants propagate. Instead of a static, uniform density medium, SNRs near the galactic center propagate into a wind swept environment with a velocity away from the galactic center, and a graded density profile. This causes these SNRs to be non - spherical, and to evolve faster than their galactic counterparts. We develop an analytic theory for the evolution of explosions within a stellar wind, and verify it using a hydrodynamic code. We show that such explosions can evolve in one of three possible morphologies. Using these results we discuss the association between the two SNRs (SGR East and SGR A's bipolar radio/X-ray Lobes) and the two neutron stars (the cannonball and SGR J1745-2900) near the galactic center. We show that, given the morphologies of the SNR and positions of the neutron stars, the only possible association is between SGR A's bipolar radio/X-ray Lobes and SGR J1745-2900. The compact object created in the explosion of SGR East remains undetected, and the SNR of the supernova that created the cannonball has already disappeared.
△ Less
Submitted 21 August, 2016;
originally announced August 2016.
-
Analytic Asymptotic Solution to Spherical Relativistic Shock Breakout
Authors:
Almog Yalinewich,
Re'em Sari
Abstract:
We investigate the relativistic breakout of a shock wave from the surface of a star. In this process, each fluid shell is endowed with some kinetic and thermal energy by the shock, and then continues to accelerate adiabatically by converting thermal energy into kinetic energy. This problem has been previously studied for a mildly relativistic breakout, where the acceleration ends close to the surf…
▽ More
We investigate the relativistic breakout of a shock wave from the surface of a star. In this process, each fluid shell is endowed with some kinetic and thermal energy by the shock, and then continues to accelerate adiabatically by converting thermal energy into kinetic energy. This problem has been previously studied for a mildly relativistic breakout, where the acceleration ends close to the surface of the star. The current work focuses on the case where the acceleration ends at distances much greater than the radius of the star. We derive an analytic description for the hydrodynamic evolution of the ejecta in this regime, and validate it using a numerical simulation. We also provide predictions for the expected light curves and spectra from such an explosion. The relevance to astrophysical explosions is discussed, and it is shown that such events require more energy than is currently believed to result from astrophysical explosions.
△ Less
Submitted 1 August, 2016;
originally announced August 2016.
-
Grid Noise in Moving Mesh Codes: Fixing The Volume Inconsistency Problem
Authors:
Elad Steinberg,
Almog Yalinewich,
Re'em Sari
Abstract:
Current Voronoi based moving mesh hydro codes suffer from "grid noise". We identify the cause of this noise as the volume inconsistency error, where the volume that is transferred between cells is inconsistent with the hydrodynamical calculations. As a result, the codes do not achieve second order convergence. In this paper we describe how a simple fix allows Voronoi based moving mesh codes to att…
▽ More
Current Voronoi based moving mesh hydro codes suffer from "grid noise". We identify the cause of this noise as the volume inconsistency error, where the volume that is transferred between cells is inconsistent with the hydrodynamical calculations. As a result, the codes do not achieve second order convergence. In this paper we describe how a simple fix allows Voronoi based moving mesh codes to attain second order convergence. The fix is based on the understanding that the volume exchanged between cells should be consistent with the hydrodynamical calculations. We benchmark our fix with three test problems and show that it can significantly improve the computational accuracy. We also examine the effect of initial mesh initialization and present an improved model for the Green-Gauss based gradient estimator.
△ Less
Submitted 18 March, 2016;
originally announced March 2016.
-
Asymptotic Steady State Solution to a Bow Shock with an Infinite Mach Number
Authors:
Almog Yalinewich,
Re'em Sari
Abstract:
The problem of a cold gas flowing past a stationary object is considered. It is shown that at large distances from the obstacle the shock front forms a parabolic solid of revolution. The interior of the shock front is obtained by solution of the hydrodynamic equations in parabolic coordinates. The results are verified with a hydrodynamic simulation. The drag force and expected spectra are calculat…
▽ More
The problem of a cold gas flowing past a stationary object is considered. It is shown that at large distances from the obstacle the shock front forms a parabolic solid of revolution. The interior of the shock front is obtained by solution of the hydrodynamic equations in parabolic coordinates. The results are verified with a hydrodynamic simulation. The drag force and expected spectra are calculated for such shock, both in case of an optically thin and thick media. Finally, relations to astrophysical bow shocks and other analytic works on oblique shocks are discussed.
△ Less
Submitted 26 October, 2015;
originally announced October 2015.
-
Rich: Open Source Hydrodynamic Simulation on a Moving Voronoi Mesh
Authors:
Almog Yalinewich,
Elad Steinberg,
Re'em Sari
Abstract:
We present here RICH, a state of the art 2D hydrodynamic code based on Godunov's method, on an unstructured moving mesh (the acronym stands for Racah Institute Computational Hydrodynamics). This code is largely based on the code AREPO. It differs from AREPO in the interpolation and time advancement scheme as well as a novel parallelization scheme based on Voronoi tessellation. Using our code we st…
▽ More
We present here RICH, a state of the art 2D hydrodynamic code based on Godunov's method, on an unstructured moving mesh (the acronym stands for Racah Institute Computational Hydrodynamics). This code is largely based on the code AREPO. It differs from AREPO in the interpolation and time advancement scheme as well as a novel parallelization scheme based on Voronoi tessellation. Using our code we study the pros and cons of a moving mesh (in comparison to a static mesh). We also compare its accuracy to other codes. Specifically, we show that our implementation of external sources and time advancement scheme is more accurate and robust than AREPO's, when the mesh is allowed to move. We performed a parameter study of the cell rounding mechanism (Llyod iterations) and it effects. We find that in most cases a moving mesh gives better results than a static mesh, but it is not universally true. In the case where matter moves in one way, and a sound wave is traveling in the other way (such that relative to the grid the wave is not moving) a static mesh gives better results than a moving mesh. Moreover, we show that Voronoi based moving mesh schemes suffer from an error, that is resolution independent, due to inconsistencies between the flux calculation and change in the area of a cell. Our code is publicly available as open source and designed in an object oriented, user friendly way that facilitates incorporation of new algorithms and physical processes.
△ Less
Submitted 13 October, 2014;
originally announced October 2014.
-
Balancing the load: A Voronoi based scheme for parallel computations
Authors:
Elad Steinberg,
Almog Yalinewich,
Re'em Sari,
Paul Duffell
Abstract:
The use of numerical simulations in science is ever increasing and with it the computational size. In many cases single processors are no longer adequate and simulations are run on multiple core machines or supercomputers. One of the key issues when running a simulation on multiple CPUs is maintaining a proper load balance throughout the run and minimizing communications between CPUs.
We propose…
▽ More
The use of numerical simulations in science is ever increasing and with it the computational size. In many cases single processors are no longer adequate and simulations are run on multiple core machines or supercomputers. One of the key issues when running a simulation on multiple CPUs is maintaining a proper load balance throughout the run and minimizing communications between CPUs.
We propose a novel method of utilizing a Voronoi diagram to achieve a nearly perfect load balance without the need of any global redistributions of data. As a show case, we implement our method in RICH, a 2D moving mesh hydrodynamical code, but it can be extended trivially to other codes in 2D or 3D. Our tests show that this method is indeed efficient and can be used in a large variety of existing hydrodynamical codes as well as other applications.
△ Less
Submitted 16 December, 2014; v1 submitted 14 August, 2014;
originally announced August 2014.
-
Discrete Self Similarity in Filled Type I Strong Explosions
Authors:
Almog Yalinewich,
Re'em Sari
Abstract:
We present new solutions to the strong explosion problem in a non power law density profi{}le. The unperturbed self similar solutions developed by Sedov, Taylor and Von Neumann describe strong Newtonian shocks propagating into a cold gas with a density profile falling off as $r^{-ω}$, where $ω\le\frac{7-γ}{γ+1}$ (filled type I solutions), and $γ$ is the adiabatic index of the gas. The perturbation…
▽ More
We present new solutions to the strong explosion problem in a non power law density profi{}le. The unperturbed self similar solutions developed by Sedov, Taylor and Von Neumann describe strong Newtonian shocks propagating into a cold gas with a density profile falling off as $r^{-ω}$, where $ω\le\frac{7-γ}{γ+1}$ (filled type I solutions), and $γ$ is the adiabatic index of the gas. The perturbations we consider are spherically symmetric and log periodic with respect to the radius. While the unperturbed solutions are continuously self similar, the log periodicity of the density perturbations leads to a discrete self similarity of the perturbations, i.e., the solution repeats itself up to a scaling at discrete time intervals. We discuss these solutions and verify them against numerical integrations of the time dependent hydrodynamic equations. This is an extension of a previous investigation on type II solutions and helps clarifying boundary conditions for perturbations to type I self similar solutions.
△ Less
Submitted 25 July, 2014;
originally announced July 2014.
-
Angular Density Perturbations to Filled Type I Strong Explosions
Authors:
Almog Yalinewich,
Re'em Sari
Abstract:
In this paper we extend the Sedov - Taylor - Von Neumann model for a strong explosion to account for small angular and radial variations in the density. We assume that the density profile is given by $ρ\left(r,θ,φ\right)=kr^{-ω}\left(1+\varepsilon\left(\frac{r}{r_{0}}\right)^{q}Y_{lm}\left(θ,φ\right)\right)$, where $\varepsilon\ll1$ and $ω\le\frac{7-γ}{γ+1}$. In order to verify our results we comp…
▽ More
In this paper we extend the Sedov - Taylor - Von Neumann model for a strong explosion to account for small angular and radial variations in the density. We assume that the density profile is given by $ρ\left(r,θ,φ\right)=kr^{-ω}\left(1+\varepsilon\left(\frac{r}{r_{0}}\right)^{q}Y_{lm}\left(θ,φ\right)\right)$, where $\varepsilon\ll1$ and $ω\le\frac{7-γ}{γ+1}$. In order to verify our results we compare them to analytical approximations and full hydrodynamic simulations. We demonstrate how this method can be used to describe arbitrary (not just self similar) angular perturbations. This work complements our previous analysis on radial, spherically symmetric perturbations, and allows one to calculate the response of an explosion to arbitrary perturbations in the upstream density. Together, they settle an age old controversy about the inner boundary conditions
△ Less
Submitted 25 July, 2014;
originally announced July 2014.
-
Atmospheric Mass Loss During Planet Formation: The Importance of Planetesimal Impacts
Authors:
Hilke Schlichting,
Re'em Sari,
Almog Yalinewich
Abstract:
We quantify the atmospheric mass loss during planet formation by examining the contributions to atmospheric loss from both giant impacts and planetesimal accretion. Giant impacts cause global motion of the ground. Using analytic self-similar solutions and full numerical integrations we find (for isothermal atmospheres with adiabatic index ($γ=5/3$) that the local atmospheric mass loss fraction for…
▽ More
We quantify the atmospheric mass loss during planet formation by examining the contributions to atmospheric loss from both giant impacts and planetesimal accretion. Giant impacts cause global motion of the ground. Using analytic self-similar solutions and full numerical integrations we find (for isothermal atmospheres with adiabatic index ($γ=5/3$) that the local atmospheric mass loss fraction for ground velocities $v_g < 0.25 v_{esc}$ is given by $χ_{loss}=(1.71 v_g/v_{esc})^{4.9}$, where $v_{esc}$ is the escape velocity from the target. Yet, the global atmospheric mass loss is a weaker function of the impactor velocity $v_{Imp}$ and mass $m_{Imp}$ and given by $X_{loss} ~ 0.4x+1.4x^2-0.8x^3$ (isothermal atmosphere) and $X_{loss} ~ 0.4x+1.8x^2-1.2x^3$ (adiabatic atmosphere), where $x=(v_{Imp}m/v_{esc}M)$. Atmospheric mass loss due to planetesimal impacts proceeds in two different regimes: 1) Large enough impactors $m > \sqrt{2} ρ_0 (πh R)^{3/2}$ (25~km for the current Earth), are able to eject all the atmosphere above the tangent plane of the impact site, which is $h/2R$ of the whole atmosphere, where $h$, $R$ and $ρ_0$ are the atmospheric scale height, radius of the target, and its atmospheric density at the ground. 2) Smaller impactors, but above $m>4 πρ_0 h^3$ (1~km for the current Earth) are only able to eject a fraction of the atmospheric mass above the tangent plane. We find that the most efficient impactors (per unit impactor mass) for atmospheric loss are planetesimals just above that lower limit and that the current atmosphere of the Earth could have resulted from an equilibrium between atmospheric erosion and volatile delivery to the atmosphere from planetesimals. We conclude that planetesimal impacts are likely to have played a major role in atmospheric mass loss over the formation history of the terrestrial planets. (Abridged)
△ Less
Submitted 26 August, 2014; v1 submitted 24 June, 2014;
originally announced June 2014.
-
Slightly Two or Three Dimensional Self-Similar Solutions
Authors:
Re'em Sari,
J. Nate Bode,
Almog Yalinewich,
Andrew MacFadyen
Abstract:
Self similarity allows for analytic or semi-analytic solutions to many hydrodynamics problems. Most of these solutions are one dimensional. Using linear perturbation theory, expanded around such a one-dimensional solution, we find self-similar hydrodynamic solutions that are two- or three-dimensional. Since the deviation from a one-dimensional solution is small, we call these slightly two-dimensio…
▽ More
Self similarity allows for analytic or semi-analytic solutions to many hydrodynamics problems. Most of these solutions are one dimensional. Using linear perturbation theory, expanded around such a one-dimensional solution, we find self-similar hydrodynamic solutions that are two- or three-dimensional. Since the deviation from a one-dimensional solution is small, we call these slightly two-dimensional and slightly three-dimensional self-similar solutions, respectively. As an example, we treat strong spherical explosions of the second type. A strong explosion propagates into an ideal gas with negligible temperature and density profile of the form rho(r,theta,phi)=r^{-omega}[1+sigma*F(theta,phi)], where omega>3 and sigma << 1. Analytical solutions are obtained by expanding the arbitrary function F(theta,phi) in spherical harmonics. We compare our results with two dimensional numerical simulations, and find good agreement.
△ Less
Submitted 15 September, 2011;
originally announced September 2011.