-
TESS photometry of the nova eruption in V606 Vul: asymmetric photosphere and multiple ejections?
Authors:
Kirill V. Sokolovsky,
Elias Aydi,
Konstantin Malanchev,
Colin J. Burke,
Koji Mukai,
J. L. Sokoloski,
Brian D. Metzger,
Kirill E. Atapin,
Aleksandre A. Belinski,
Yu-Ching Chen,
Laura Chomiuk,
Pavol A. Dubovsky,
Claude-Andre Faucher-Giguere,
Rebekah A. Hounsell,
Natalia P. Ikonnikova,
Vsevolod Yu. Lander,
Junyao Li,
Justin D. Linford,
Amy J. Mioduszewski,
Isabella Molina,
Ulisse Munari,
Sergey A. Potanin,
Robert M. Quimby,
Michael P. Rupen,
Simone Scaringi
, et al. (46 additional authors not shown)
Abstract:
Lightcurves of many classical novae deviate from the canonical "fast rise - smooth decline" pattern and display complex variability behavior. We present the first TESS-space-photometry-based investigation of this phenomenon. We use TESS Sector 41 full-frame images to extract a lightcurve of the slow Galactic nova V606 Vul that erupted nine days prior to the start of the TESS observations. The ligh…
▽ More
Lightcurves of many classical novae deviate from the canonical "fast rise - smooth decline" pattern and display complex variability behavior. We present the first TESS-space-photometry-based investigation of this phenomenon. We use TESS Sector 41 full-frame images to extract a lightcurve of the slow Galactic nova V606 Vul that erupted nine days prior to the start of the TESS observations. The lightcurve covers the first of two major peaks of V606 Vul that was reached 19 days after the start of the eruption. The nova reached its brightest visual magnitude V=9.9 in its second peak 64 days after the eruption onset, following the completion of Sector 41 observations. To increase the confidence level of the extracted lightcurve, we performed the analysis using four different codes implementing the aperture photometry (Lightkurve, VaST) and image subtraction (TESSreduce, tequila_shots) and find good agreement between them. We performed ground-based photometric and spectroscopic monitoring to complement the TESS data. The TESS lightcurve reveals two features: periodic variations (0.12771d, 0.01mag average peak-to-peak amplitude) that disappeared when the source was within 1mag of peak optical brightness and a series of isolated mini-flares (with peak-to-peak amplitudes of up to 0.5mag) appearing at seemingly random times. We interpret the periodic variations as the result of azimuthal asymmetry of the photosphere engulfing the nova-hosting binary that was distorted by and rotating with the binary. Whereas we use spectra to associate the two major peaks in the nova lightcurve with distinct episodes of mass ejection, the origin of mini-flares remains elusive.
△ Less
Submitted 8 November, 2023;
originally announced November 2023.
-
Frequency-dependent Discrete Implicit Monte-Carlo Scheme for the Radiative Transfer Equation
Authors:
Elad Steinberg,
Shay I. Heizler
Abstract:
This work generalizes the discrete implicit Monte-Carlo (DIMC) method for modeling the radiative transfer equation from a gray treatment to an frequency-dependent one. The classic implicit Monte-Carlo (IMC) algorithm, that has been used for several decades, suffers from a well-known numerical problem, called teleportation, where the photons might propagate faster than the exact solution due to the…
▽ More
This work generalizes the discrete implicit Monte-Carlo (DIMC) method for modeling the radiative transfer equation from a gray treatment to an frequency-dependent one. The classic implicit Monte-Carlo (IMC) algorithm, that has been used for several decades, suffers from a well-known numerical problem, called teleportation, where the photons might propagate faster than the exact solution due to the finite size of the spatial and temporal resolution. The Semi-analog Monte-Carlo algorithm proposed the use of two kinds of particles, photons and material particles that are born when a photon is absorbed. The material particle can `propagate' only by transforming into a photon, due to black-body emission. While this algorithm produces a teleportation-free result, it is noisier results compared to IMC due to the discrete nature of the absorption-emission process. In a previous work [Steinberg and Heizler, ApJS, 258:14 (2022)], proposed a gray version of DIMC, that makes use of two kinds of particles, and therefore has teleportation-free results, but also uses the continuous absorption algorithm of IMC, yielding smoother results. This work is a direct frequency-dependent (energy-dependent) generalization of the DIMC algorithm. We find in several one and two dimensional benchmarks, that the new frequency-dependent DIMC algorithm yields teleportation-free results on one hand, and smooth results with IMC-like noise level.
△ Less
Submitted 12 March, 2023;
originally announced March 2023.
-
The multi-wavelength view of shocks in the fastest nova V1674 Her
Authors:
K. V. Sokolovsky,
T. J. Johnson,
S. Buson,
P. Jean,
C. C. Cheung,
K. Mukai,
L. Chomiuk,
E. Aydi,
B. Molina,
A. Kawash,
J. D. Linford,
A. J. Mioduszewski,
M. P. Rupen,
J. L. Sokoloski,
M. N. Williams,
E. Steinberg,
I. Vurm,
B. D. Metzger,
K. L. Page,
M. Orio,
R. M. Quimby,
A. W. Shafter,
H. Corbett,
S. Bolzoni,
J. DeYoung
, et al. (19 additional authors not shown)
Abstract:
Classical novae are shock-powered multi-wavelength transients triggered by a thermonuclear runaway on an accreting white dwarf. V1674 Her is the fastest nova ever recorded (time to declined by two magnitudes is t_2=1.1 d) that challenges our understanding of shock formation in novae. We investigate the physical mechanisms behind nova emission from GeV gamma-rays to cm-band radio using coordinated…
▽ More
Classical novae are shock-powered multi-wavelength transients triggered by a thermonuclear runaway on an accreting white dwarf. V1674 Her is the fastest nova ever recorded (time to declined by two magnitudes is t_2=1.1 d) that challenges our understanding of shock formation in novae. We investigate the physical mechanisms behind nova emission from GeV gamma-rays to cm-band radio using coordinated Fermi-LAT, NuSTAR, Swift and VLA observations supported by optical photometry. Fermi-LAT detected short-lived (18 h) 0.1-100 GeV emission from V1674 Her that appeared 6 h after the eruption began; this was at a level of (1.6 +/- 0.4)x10^-6 photons cm^-2 s^-1. Eleven days later, simultaneous NuSTAR and Swift X-ray observations revealed optically thin thermal plasma shock-heated to kT_shock = 4 keV. The lack of a detectable 6.7 keV Fe K_alpha emission suggests super-solar CNO abundances. The radio emission from V1674 Her was consistent with thermal emission at early times and synchrotron at late times. The radio spectrum steeply rising with frequency may be a result of either free-free absorption of synchrotron and thermal emission by unshocked outer regions of the nova shell or the Razin-Tsytovich effect attenuating synchrotron emission in dense plasma. The development of the shock inside the ejecta is unaffected by the extraordinarily rapid evolution and the intermediate polar host of this nova.
△ Less
Submitted 21 March, 2023; v1 submitted 6 February, 2023;
originally announced February 2023.
-
Stream-Disk Shocks as the Origins of Peak Light in Tidal Disruption Events
Authors:
Elad Steinberg,
Nicholas C. Stone
Abstract:
Tidal disruption events occur when stars are ripped apart by massive black holes, and result in highly luminous, multi-wavelength flares. Optical/UV observations of tidal disruption events (TDEs) contradict simple models of TDE emission, but the debate between alternative models (e.g. shock power or reprocessed accretion power remains unsettled, as the dynamic range of the problem has so far preve…
▽ More
Tidal disruption events occur when stars are ripped apart by massive black holes, and result in highly luminous, multi-wavelength flares. Optical/UV observations of tidal disruption events (TDEs) contradict simple models of TDE emission, but the debate between alternative models (e.g. shock power or reprocessed accretion power remains unsettled, as the dynamic range of the problem has so far prevented ab initio hydrodynamical simulations. Consequently, past simulations have resorted to unrealistic parameter choices, artificial mass injection schemes or very short run-times. Here we present a 3D radiation-hydrodynamic simulation of a TDE flare from disruption to peak emission, with typical astrophysical parameters. At early times, shocks near pericenter power the light curve and a novel source of X-ray emission, but circularization and outflows are inefficient. Near peak light, stream-disk shocks efficiently circularize returning debris, power stronger outflows, and reproduce observed peak optical/UV luminosities. Peak emission in this simulation is shock-powered, but upper limits on accretion power become competitive near peak light as circularization runs away. This simulation shows how deterministic predictions of TDE light curves and spectra can be calculated using moving-mesh hydrodynamics algorithms.
△ Less
Submitted 25 December, 2023; v1 submitted 21 June, 2022;
originally announced June 2022.
-
A New Discrete Implicit Monte Carlo Scheme for Simulating Radiative Transfer Problems
Authors:
Elad Steinberg,
Shay I. Heizler
Abstract:
We present a new algorithm for radiative transfer-based on a statistical Monte Carlo approach-that does not suffer from teleportation effects, on the one hand, and yields smooth results, on the other hand. Implicit Monte Carlo (IMC) techniques for modeling radiative transfer have existed from the 1970s. When they are used for optically thick problems, however, the basic algorithm suffers from "tel…
▽ More
We present a new algorithm for radiative transfer-based on a statistical Monte Carlo approach-that does not suffer from teleportation effects, on the one hand, and yields smooth results, on the other hand. Implicit Monte Carlo (IMC) techniques for modeling radiative transfer have existed from the 1970s. When they are used for optically thick problems, however, the basic algorithm suffers from "teleportation" errors, where the photons propagate faster than the exact physical behavior, due to the absorption-blackbody emission processes. One possible solution is to use semianalog Monte Carlo, in its new implicit form (ISMC), which uses two kinds of particles, photons and discrete material particles. This algorithm yields excellent teleportation-free results, but it also produces noisier solutions (relative to classic IMC), due to its discrete nature. Here, we derive a new Monte Carlo algorithm, Discrete Implicit Monte Carlo (DIMC), which also uses the idea of two kinds of discrete particles, and thus does not suffer from teleportation errors. DIMC implements the IMC discretization and creates new radiation photons for each time step, unlike ISMC. Using the continuous absorption technique, DIMC yields smooth results like classic IMC. One of the main elements of the algorithm is the avoidance of the explosion of the particle population, by using particle merging. We test the new algorithm on 1D and 2D cylindrical problems, and show that it yields smooth, teleportation-free results. We finish by demonstrating the power of the new algorithm on a classic radiative hydrodynamic problem-an opaque radiative shock wave. This demonstrates the power of the new algorithm for astrophysical scenarios.
△ Less
Submitted 11 January, 2022; v1 submitted 30 August, 2021;
originally announced August 2021.
-
The first nova eruption in a novalike variable: YZ Ret as seen in X-rays and gamma-rays
Authors:
Kirill V. Sokolovsky,
Kwan-Lok Li,
Raimundo Lopes de Oliveira,
Jan-Uwe Ness,
Koji Mukai,
Laura Chomiuk,
Elias Aydi,
Elad Steinberg,
Indrek Vurm,
Brian D. Metzger,
Aliya-Nur Babul,
Adam Kawash,
Justin D. Linford,
Thomas Nelson,
Kim L. Page,
Michael P. Rupen,
Jennifer L. Sokoloski,
Jay Strader,
David Kilkenny
Abstract:
Peaking at 3.7 mag on 2020 July 11, YZ Ret was the second-brightest nova of the decade. The nova's moderate proximity (2.7 kpc from Gaia) provided an opportunity to explore its multi-wavelength properties in great detail. Here we report on YZ Ret as part of a long-term project to identify the physical mechanisms responsible for high-energy emission in classical novae. We use simultaneous Fermi/LAT…
▽ More
Peaking at 3.7 mag on 2020 July 11, YZ Ret was the second-brightest nova of the decade. The nova's moderate proximity (2.7 kpc from Gaia) provided an opportunity to explore its multi-wavelength properties in great detail. Here we report on YZ Ret as part of a long-term project to identify the physical mechanisms responsible for high-energy emission in classical novae. We use simultaneous Fermi/LAT and NuSTAR observations complemented by XMM-Newton X-ray grating spectroscopy to probe the physical parameters of the shocked ejecta and the nova-hosting white dwarf. The XMM-Newton observations revealed a super-soft X-ray emission which is dominated by emission lines of CV, CVI, NVI, NVII, and OVIII rather than a blackbody-like continuum, suggesting CO-composition of the white dwarf in a high-inclination binary system. Fermi/LAT detected YZ Ret for 15 days with the gamma-ray spectrum best described by a power law with an exponential cut-off at 1.9 +/-0.6 GeV. In stark contrast with theoretical predictions and in keeping with previous NuSTAR observations of Fermi-detected classical novae (V5855 Sgr and V906 Car), the 3.5-78 keV X-ray emission is found to be two orders of magnitude fainter than the GeV emission. The X-ray emission observed by NuSTAR is consistent with a single-temperature thermal plasma. We detect no non-thermal tail of the GeV emission expected to extend down to the NuSTAR band. NuSTAR observations continue to challenge theories of high-energy emission from shocks in novae.
△ Less
Submitted 21 May, 2022; v1 submitted 6 August, 2021;
originally announced August 2021.
-
Multi-Frequency Implicit Semi-analog Monte-Carlo (ISMC) Radiative Transfer Solver in Two-Dimensions (without Teleportation)
Authors:
Elad Steinberg,
Shay I. Heizler
Abstract:
We study the multi-dimensional radiative transfer phenomena using the ISMC scheme, in both gray and multi-frequency problems. Implicit Monte-Carlo (IMC) schemes have been in use for five decades. The basic algorithm yields teleportation errors, where photons propagate faster than the correct heat front velocity. Recently [Poëtte and Valentin, J. Comp. Phys., 412, 109405 (2020)], a new implicit sch…
▽ More
We study the multi-dimensional radiative transfer phenomena using the ISMC scheme, in both gray and multi-frequency problems. Implicit Monte-Carlo (IMC) schemes have been in use for five decades. The basic algorithm yields teleportation errors, where photons propagate faster than the correct heat front velocity. Recently [Poëtte and Valentin, J. Comp. Phys., 412, 109405 (2020)], a new implicit scheme based on the semi-analog scheme was presented and tested in several one-dimensional gray problems. In this scheme, the material energy of the cell is carried by material-particles, and the photons are produced only from existing material particles. As a result, the teleportation errors vanish, due to the infinite discrete spatial accuracy of the scheme. We examine the validity of the new scheme in two-dimensional problems, both in Cartesian and Cylindrical geometries. Additionally, we introduce an expansion of the new scheme for multi-frequency problems. We show that the ISMC scheme presents excellent results without teleportation errors in a large number of benchmarks, especially against the slow classic IMC convergence.
△ Less
Submitted 5 January, 2022; v1 submitted 5 August, 2021;
originally announced August 2021.
-
Early spectral evolution of classical novae: consistent evidence for multiple distinct outflows
Authors:
E. Aydi,
L. Chomiuk,
L. Izzo,
E. J. Harvey,
J. Leahy-McGregor,
J. Strader,
D. A. H. Buckley,
K. V. Sokolovsky,
A. Kawash,
C. S. Kochanek,
J. D. Linford,
B. D. Metzger,
K. Mukai,
M. Orio,
B. J. Shappee,
L. Shishkovsky,
E. Steinberg,
S. J. Swihart,
J. L. Sokoloski,
F. M. Walter,
P. A. Woudt
Abstract:
The physical mechanism driving mass ejection during a nova eruption is still poorly understood. Possibilities include ejection in a single ballistic event, a common envelope interaction, a continuous wind, or some combination of these processes. Here we present a study of 12 Galactic novae, for which we have pre-maximum high-resolution spectroscopy. All 12 novae show the same spectral evolution. B…
▽ More
The physical mechanism driving mass ejection during a nova eruption is still poorly understood. Possibilities include ejection in a single ballistic event, a common envelope interaction, a continuous wind, or some combination of these processes. Here we present a study of 12 Galactic novae, for which we have pre-maximum high-resolution spectroscopy. All 12 novae show the same spectral evolution. Before optical peak, they show a slow P Cygni component. After peak a fast component quickly arises, while the slow absorption remains superimposed on top of it, implying the presence of at least two physically distinct flows. For novae with high-cadence monitoring, a third, intermediate-velocity component is also observed.
These observations are consistent with a scenario where the slow component is associated with the initial ejection of the accreted material and the fast component with a radiation-driven wind from the white dwarf. When these flows interact, the slow flow is swept up by the fast flow, producing the intermediate component. These colliding flows may produce the gamma-ray emission observed in some novae. Our spectra also show that the transient heavy element absorption lines seen in some novae have the same velocity structure and evolution as the other lines in the spectrum, implying an association with the nova ejecta rather than a pre-existing circumbinary reservoir of gas or material ablated from the secondary. While this basic scenario appears to qualitatively reproduce multi-wavelength observations of classical novae, substantial theoretical and observational work is still needed to untangle the rich diversity of nova properties.
△ Less
Submitted 14 October, 2020;
originally announced October 2020.
-
X-ray Spectroscopy of the gamma-ray Brightest Nova V906 Car (ASASSN-18fv)
Authors:
Kirill V. Sokolovsky,
Koji Mukai,
Laura Chomiuk,
Raimundo Lopes de Oliveira,
Elias Aydi,
Kwan-Lok Li,
Elad Steinberg,
Indrek Vurm,
Brian D. Metzger,
Adam Kawash,
Justin D. Linford,
Amy J. Mioduszewski,
Thomas Nelson,
Jan-Uwe Ness,
Kim L. Page,
Michael P. Rupen,
Jennifer L. Sokoloski,
Jay Strader
Abstract:
Shocks in gamma-ray emitting classical novae are expected to produce bright thermal and non-thermal X-rays. We test this prediction with simultaneous NuSTAR and Fermi/LAT observations of nova V906 Car, which exhibited the brightest GeV gamma-ray emission to date. The nova is detected in hard X-rays while it is still gamma-ray bright, but contrary to simple theoretical expectations, the detected 3.…
▽ More
Shocks in gamma-ray emitting classical novae are expected to produce bright thermal and non-thermal X-rays. We test this prediction with simultaneous NuSTAR and Fermi/LAT observations of nova V906 Car, which exhibited the brightest GeV gamma-ray emission to date. The nova is detected in hard X-rays while it is still gamma-ray bright, but contrary to simple theoretical expectations, the detected 3.5-78 keV emission of V906 Car is much weaker than the simultaneously observed >100 MeV emission. No non-thermal X-ray emission is detected, and our deep limits imply that the gamma-rays are likely hadronic. After correcting for substantial absorption (N_H ~ 2 x 10^23 cm^-2), the thermal X-ray luminosity (from a 9 keV optically-thin plasma) is just ~2% of the gamma-ray luminosity. We consider possible explanations for the low thermal X-ray luminosity, including the X-rays being suppressed by corrugated, radiative shock fronts or the X-rays from the gamma-ray producing shock are hidden behind an even larger absorbing column (N_H >10^25 cm^-2). Adding XMM-Newton and Swift/XRT observations to our analysis, we find that the evolution of the intrinsic X-ray absorption requires the nova shell to be expelled 24 days after the outburst onset. The X-ray spectra show that the ejecta are enhanced in nitrogen and oxygen, and the nova occurred on the surface of a CO-type white dwarf. We see no indication of a distinct super-soft phase in the X-ray lightcurve, which, after considering the absorption effects, may point to a low mass of the white dwarf hosting the nova.
△ Less
Submitted 15 July, 2020;
originally announced July 2020.
-
Direct evidence for shock-powered optical emission in a nova
Authors:
Elias Aydi,
Kirill V. Sokolovsky,
Laura Chomiuk,
Elad Steinberg,
Kwan Lok Li,
Indrek Vurm,
Brian D. Metzger,
Jay Strader,
Koji Mukai,
Ondřej Pejcha,
Ken J. Shen,
Gregg A. Wade,
Rainer Kuschnig,
Anthony F. J. Moffat,
Herbert Pablo,
Andrzej Pigulski,
Adam Popowicz,
Werner Weiss,
Konstanze Zwintz,
Luca Izzo,
Karen R. Pollard,
Gerald Handler,
Stuart D. Ryder,
Miroslav D. Filipović,
Rami Z. E. Alsaberi
, et al. (17 additional authors not shown)
Abstract:
Classical novae are thermonuclear explosions that occur on the surfaces of white dwarf stars in interacting binary systems (Bode & Evans 2008). It has long been thought that the luminosity of classical novae is powered by continued nuclear burning on the surface of the white dwarf after the initial runaway (Gallaher & Starrfield 1978). However, recent observations of GeV $γ$-rays from classical no…
▽ More
Classical novae are thermonuclear explosions that occur on the surfaces of white dwarf stars in interacting binary systems (Bode & Evans 2008). It has long been thought that the luminosity of classical novae is powered by continued nuclear burning on the surface of the white dwarf after the initial runaway (Gallaher & Starrfield 1978). However, recent observations of GeV $γ$-rays from classical novae have hinted that shocks internal to the nova ejecta may dominate the nova emission. Shocks have also been suggested to power the luminosity of events as diverse as stellar mergers (Metzger & Pejcha 2017), supernovae (Moriya et al. 2018), and tidal disruption events (Roth et al. 2016), but observational confirmation has been lacking. Here we report simultaneous space-based optical and $γ$-ray observations of the 2018 nova V906 Carinae (ASASSN-18fv), revealing a remarkable series of distinct correlated flares in both bands. The optical and $γ$-ray flares occur simultaneously, implying a common origin in shocks. During the flares, the nova luminosity doubles, implying that the bulk of the luminosity is shock-powered. Furthermore, we detect concurrent but weak X-ray emission from deeply embedded shocks, confirming that the shock power does not appear in the X-ray band and supporting its emergence at longer wavelengths. Our data, spanning the spectrum from radio to $γ$-ray, provide direct evidence that shocks can power substantial luminosity in classical novae and other optical transients.
△ Less
Submitted 12 April, 2020;
originally announced April 2020.
-
Internal Shocks from Variable Outflows in Classical Novae
Authors:
Elad Steinberg,
Brian D. Metzger
Abstract:
We present one-dimensional hydrodynamical simulations including radiative losses, of internal shocks in the outflows from classical novae, to explore the role of shocks in powering multi-wavelength emission from radio to gamma-ray wavelengths. Observations support a picture in which the initial phases of some novae generate a slow, equatorially-focused outflow (directly from the outer Lagrange poi…
▽ More
We present one-dimensional hydrodynamical simulations including radiative losses, of internal shocks in the outflows from classical novae, to explore the role of shocks in powering multi-wavelength emission from radio to gamma-ray wavelengths. Observations support a picture in which the initial phases of some novae generate a slow, equatorially-focused outflow (directly from the outer Lagrange point, or from a circumbinary disk), which then transitions to, or is overtaken by, a faster more isotropic outflow from the white dwarf which collides and shocks the slower flow, powering gamma-ray and optical emission through reprocessing by the ejecta. However, the common occurence of multiple peaks in nova light curves suggests that the outflow's acceleration need not be monotonic, but instead can involve successive transitions between "fast" and "slow" modes. Such a time-fluctuating outflow velocity naturally can reproduce several observed properties of nova, such as correlated gamma-ray and optical flares, expansion of the photosphere coincident with (though lagging slightly) the peak flare luminosity, and complex time-evolution of spectral lines (including accelerating, decelerating, and merging velocity components). While the shocks are still deeply embedded during the gamma-ray emission, the onset of $\sim$ keV X-ray and $\sim 10$ GHz radio synchrotron emission is typically delayed until the forward shock of the outermost monolithic shell (created by merger of multiple internal shock-generated shells) reaches a sufficiently low column through the dense external medium generated by the earliest phase of the outburst.
△ Less
Submitted 21 November, 2019; v1 submitted 5 August, 2019;
originally announced August 2019.
-
Astro2020 Science White Paper: A Shocking Shift in Paradigm for Classical Novae
Authors:
Laura Chomiuk,
Elias Aydi,
Aliya-Nur Babul,
Andrea Derdzinski,
Adam Kawash,
Kwan-Lok Li,
Justin Linford,
Brian D. Metzger,
Koji Mukai,
Michael P. Rupen,
Jennifer Sokoloski,
Kirill Sokolovsky,
Elad Steinberg
Abstract:
The discovery of GeV gamma-rays from classical novae has led to a reassessment of these garden-variety explosions, and highlighted their importance for understanding radiative shocks, particle acceleration, and dust formation in more exotic, distant transients. Recent collaboration between observers and theorists has revealed that shocks in novae are energetically important, and can even dominate…
▽ More
The discovery of GeV gamma-rays from classical novae has led to a reassessment of these garden-variety explosions, and highlighted their importance for understanding radiative shocks, particle acceleration, and dust formation in more exotic, distant transients. Recent collaboration between observers and theorists has revealed that shocks in novae are energetically important, and can even dominate their bolometric luminosity. Shocks may also explain long-standing mysteries in novae such as dust production, super-Eddington luminosities, and `flares' in optical light curves. Here, we highlight the multi-wavelength facilities of the next decade that will further test our nova shock model and fulfill the promise of novae as powerful astrophysical laboratories.
△ Less
Submitted 19 March, 2019;
originally announced March 2019.
-
Thawing the frozen-in approximation: implications for self-gravity in deeply plunging tidal disruption events
Authors:
Elad Steinberg,
Eric R. Coughlin,
Nicholas C. Stone,
Brian D. Metzger
Abstract:
The tidal destruction of a star by a massive black hole, known as a tidal disruption event (TDE), is commonly modeled using the "frozen-in" approximation. Under this approximation, the star maintains exact hydrostatic balance prior to entering the tidal sphere (radius $r_{\rm t}$), after which point its internal pressure and self-gravity become instantaneously negligible and the debris undergoes b…
▽ More
The tidal destruction of a star by a massive black hole, known as a tidal disruption event (TDE), is commonly modeled using the "frozen-in" approximation. Under this approximation, the star maintains exact hydrostatic balance prior to entering the tidal sphere (radius $r_{\rm t}$), after which point its internal pressure and self-gravity become instantaneously negligible and the debris undergoes ballistic free fall. We present a suite of hydrodynamical simulations of TDEs with high penetration factors $β\equiv r_{\rm t}/r_{\rm p} = 5-7$, where $r_{\rm p}$ is the pericenter of the stellar center of mass, calculated using a Voronoi-based moving-mesh technique. We show that basic assumptions of the frozen-in model, such as the neglect of self-gravity inside $r_{\rm t}$, are violated. Indeed, roughly equal fractions of the final energy spread accumulate exiting and entering the tidal sphere, though the frozen-in prediction is correct at the order-of-magnitude level. We also show that an $\mathcal{O}(1)$ fraction of the debris mass remains transversely confined by self-gravity even for large $β$ which has implications for the radio emission from the unbound debris and, potentially, for the circularization efficiency of the bound streams.
△ Less
Submitted 2 April, 2019; v1 submitted 9 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.
-
The Multi-Dimensional Structure of Radiative Shocks: Suppressed Thermal X-rays and Relativistic Ion Acceleration
Authors:
Elad Steinberg,
Brian D. Metzger
Abstract:
Radiative shocks, behind which gas cools faster than the dynamical time, play a key role in many astrophysical transients, including classical novae and young supernovae interacting with circumstellar material. The dense layer behind high Mach number $\mathcal{M} \gg 1$ radiative shocks is susceptible to thin-shell instabilities, creating a "corrugated" shock interface. We present two and three-di…
▽ More
Radiative shocks, behind which gas cools faster than the dynamical time, play a key role in many astrophysical transients, including classical novae and young supernovae interacting with circumstellar material. The dense layer behind high Mach number $\mathcal{M} \gg 1$ radiative shocks is susceptible to thin-shell instabilities, creating a "corrugated" shock interface. We present two and three-dimensional hydrodynamical simulations of optically-thin radiative shocks to study their thermal radiation and acceleration of non-thermal relativistic ions. We employ a moving-mesh code and a specialized numerical technique to eliminate artificial heat conduction across grid cells. The fraction of the shock's luminosity $L_{\rm tot}$ radiated at X-ray temperatures $kT_{\rm sh} \approx (3/16)μm_p v_{\rm sh}^{2}$ expected from a one-dimensional analysis is suppressed by a factor $L(>T_{\rm sh}/3)/L_{\rm tot} \approx 4.5/\mathcal{M}^{4/3}$ for $\mathcal{M} \approx 4-36$. This suppression results in part from weak shocks driven into under-pressured cold filaments by hot shocked gas, which sap thermal energy from the latter faster than it is radiated. Combining particle-in-cell simulation results for diffusive shock acceleration with the inclination angle distribution across the shock (relative to an upstream magnetic field in the shock plane$-$the expected geometry for transient outflows), we predict the efficiency and energy spectrum of ion acceleration. Though negligible acceleration is predicted for adiabatic shocks, the corrugated shock front enables local regions to satisfy the quasi-parallel magnetic field geometry required for efficient acceleration, resulting in an average acceleration efficiency of $ε_{\rm nth} \sim 0.005-0.02$ for $\mathcal{M} \approx 12-36$, in agreement with modeling of the gamma-ray nova ASASSN-16ma.
△ Less
Submitted 8 June, 2018; v1 submitted 8 May, 2018;
originally announced May 2018.
-
Instability of Supersonic Cold Streams Feeding Galaxies II. Nonlinear Evolution of Surface and Body Modes of Kelvin-Helmholtz Instability
Authors:
Dan Padnos,
Nir Mandelker,
Yuval Birnboim,
Avishai Dekel,
Mark R. Krumholz,
Elad Steinberg
Abstract:
As part of our long-term campaign to understand how cold streams feed massive galaxies at high redshift, we study the Kelvin-Helmholtz instability (KHI) of a supersonic, cold, dense gas stream as it penetrates through a hot, dilute circumgalactic medium (CGM). A linear analysis (Paper I) showed that, for realistic conditions, KHI may produce nonlinear perturbations to the stream during infall. The…
▽ More
As part of our long-term campaign to understand how cold streams feed massive galaxies at high redshift, we study the Kelvin-Helmholtz instability (KHI) of a supersonic, cold, dense gas stream as it penetrates through a hot, dilute circumgalactic medium (CGM). A linear analysis (Paper I) showed that, for realistic conditions, KHI may produce nonlinear perturbations to the stream during infall. Therefore, we proceed here to study the nonlinear stage of KHI, still limited to a two-dimensional slab with no radiative cooling or gravity. Using analytic models and numerical simulations, we examine stream breakup, deceleration and heating via surface modes and body modes. The relevant parameters are the density contrast between stream and CGM ($δ$), the Mach number of the stream velocity with respect to the CGM ($M_{\rm b}$) and the stream radius relative to the halo virial radius ($R_{\rm s}/R_{\rm v}$). We find that sufficiently thin streams disintegrate prior to reaching the central galaxy. The condition for breakup ranges from $R_{\rm s} < 0.03 R_{\rm v}$ for $(M_{\rm b} \sim 0.75, δ\sim 10)$ to $R_{\rm s} < 0.003 R_{\rm v}$ for $(M_{\rm b} \sim 2.25, δ\sim 100)$. However, due to the large stream inertia, KHI has only a small effect on the stream inflow rate and a small contribution to heating and subsequent Lyman-$α$ cooling emission.
△ Less
Submitted 24 March, 2018;
originally announced March 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.
-
Probing the Gas Density in our Galactic Center: Moving Mesh Simulations of G2
Authors:
Elad Steinberg,
Re'em Sari,
Orly Gnat,
Stefan Gillessen,
Philipp Plewa,
Reinhard Genzel,
Frank Eisenhauer,
Thomas Ott,
Oliver Pfuhl,
Maryam Habibi,
Idel Waisberg,
Sebastiano von Fellenberg,
Jason Dexter,
Michi Bauböck
Abstract:
The G2 object has recently passed its pericenter passage in our Galactic Center. While the $Br_γ$ emission shows clear signs of tidal interaction, the change in the observed luminosity is only of about a factor of 2, in contention with all previous predictions.
We present high resolution simulations performed with the moving mesh code, RICH, together with simple analytical arguments that reprodu…
▽ More
The G2 object has recently passed its pericenter passage in our Galactic Center. While the $Br_γ$ emission shows clear signs of tidal interaction, the change in the observed luminosity is only of about a factor of 2, in contention with all previous predictions.
We present high resolution simulations performed with the moving mesh code, RICH, together with simple analytical arguments that reproduce the observed $Br_γ$ emission. In our model, G2 is a gas cloud that undergoes tidal disruption in a dilute ambient medium. We find that during pericenter passage, the efficient cooling of the cloud results in a vertical collapse, compressing the cloud by a factor of $\sim5000$. By properly taking into account the ionization state of the gas, we find that the cloud is UV starved and are able to reproduce the observed $Br_γ$ luminosity.
For densities larger than $\approx500\;\mathrm{cm}^{-3}$ at pericenter, the cloud fragments, due to cooling instabilities and the emitted radiation is inconsistent with observations. For lower densities, the cloud survives the pericenter passage intact and its emitted radiation matches the observed lightcurve.
From the duration of $Br_γ$ emission which contains both redshifted and blueshifted components, we show that the cloud is not spherical but rather elongated with a size ratio of 4 at year 2001. The simulated cloud's elongation grows as it travels towards pericenter and is consistent with observations, due to viewing angles. The simulation is also consistent with having a spherical shape at apocenter.
△ Less
Submitted 28 September, 2017; v1 submitted 29 May, 2017;
originally announced May 2017.
-
Instability of Supersonic Cold Streams Feeding Galaxies I: Linear Kelvin-Helmholtz Instability with Body Modes
Authors:
Nir Mandelker,
Dan Padnos,
Avishai Dekel,
Yuval Birnboim,
Andreas Burkert,
Mark R. Krumholz,
Elad Steinberg
Abstract:
Massive galaxies at high redshift are predicted to be fed from the cosmic web by narrow, dense, cold streams. These streams penetrate supersonically through the hot medium encompassed by a stable shock near the virial radius of the dark-matter halo. Our long-term goal is to explore the heating and dissipation rate of the streams and their fragmentation and possible breakup, in order to understand…
▽ More
Massive galaxies at high redshift are predicted to be fed from the cosmic web by narrow, dense, cold streams. These streams penetrate supersonically through the hot medium encompassed by a stable shock near the virial radius of the dark-matter halo. Our long-term goal is to explore the heating and dissipation rate of the streams and their fragmentation and possible breakup, in order to understand how galaxies are fed, and how this affects their star-formation rate and morphology. We present here the first step, where we analyze the linear Kelvin-Helmholtz instability (KHI) of a cold, dense slab or cylinder flowing through a hot, dilute medium in the transonic regime. The current analysis is limited to the adiabatic case with no gravity and assuming equal pressure in the stream and the medium. By analytically solving the linear dispersion relation, we find a transition from a dominance of the familiar rapidly growing surface modes in the subsonic regime to more slowly growing body modes in the supersonic regime. The system is parameterized by three parameters: the density contrast between the stream and the medium, the Mach number of stream velocity with respect to the medium, and the stream width with respect to the halo virial radius. We find that a realistic choice for these parameters places the streams near the mode transition, with the KHI exponential-growth time in the range 0.01-10 virial crossing times for a perturbation wavelength comparable to the stream width. We confirm our analytic predictions with idealized hydrodynamical simulations. Our linear-KHI estimates thus indicate that KHI may in principle be effective in the evolution of streams by the time they reach the galaxy. More definite conclusions await the extension of the analysis to the nonlinear regime and the inclusion of cooling, thermal conduction, the halo potential well, self-gravity and magnetic fields.
△ Less
Submitted 18 September, 2016; v1 submitted 20 June, 2016;
originally announced June 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.
-
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.
-
The Galactic Center cloud G2 and its gas streamer
Authors:
Oliver Pfuhl,
Stefan Gillessen,
Frank Eisenhauer,
Reinhard Genzel,
Philipp M. Plewa,
Thomas Ott,
Alessandro Ballone,
Marc Schartmann,
Andreas Burkert,
Tobias K. Fritz,
Re'em Sari,
Elad Steinberg,
Ann-Marie Madigan
Abstract:
We present new, deep near-infrared SINFONI @ VLT integral field spectroscopy of the gas cloud G2 in the Galactic Center, from late summer 2013 and spring 2014. G2 is visible in recombination line emission. The spatially resolved kinematic data track the ongoing tidal disruption. As expected for an observation near pericenter passage, roughly half of the gas in 2014 is found at the redshifted, pre-…
▽ More
We present new, deep near-infrared SINFONI @ VLT integral field spectroscopy of the gas cloud G2 in the Galactic Center, from late summer 2013 and spring 2014. G2 is visible in recombination line emission. The spatially resolved kinematic data track the ongoing tidal disruption. As expected for an observation near pericenter passage, roughly half of the gas in 2014 is found at the redshifted, pre-pericenter side of the orbit, while the other half is at the post-pericenter, blueshifted side. We also present an orbital solution for the gas cloud G1, which was discovered a decade ago in L'-band images when it was spatially almost coincident with Sgr A*. The orientation of the G1 orbit in the three angles is almost identical to the one of G2, but it has a lower eccentricity and smaller semi-major axis. We show that the observed astrometric positions and radial velocities of G1 are compatible with the G2 orbit, assuming that (i) G1 was originally on the G2 orbit preceding G2 by 13 years and (ii) a simple drag force acted on it during pericenter passage. Taken together with the previously described tail of G2, which we detect in recombination line emission and thermal broadband emission, we propose that G2 may be a bright knot in a much more extensive gas streamer. This matches purely gaseous models for G2, such as a stellar wind clump or the tidal debris from a partial disruption of a star.
△ Less
Submitted 5 November, 2014; v1 submitted 16 July, 2014;
originally announced July 2014.
-
Spins of large asteroids: Hint on a primordial distribution of their spin rates
Authors:
Elad Steinberg,
Re'em Sari
Abstract:
The Asteroid Belt and the Kuiper Belt are relics from the formation of our solar system. Understanding the size and spin distribution of the two belts is crucial for a deeper understanding of the formation of our solar system and the dynamical process that govern it. In this paper, we investigate the effect of collisions on the evolution of the spin distribution of asteroids and KBO's.
We find t…
▽ More
The Asteroid Belt and the Kuiper Belt are relics from the formation of our solar system. Understanding the size and spin distribution of the two belts is crucial for a deeper understanding of the formation of our solar system and the dynamical process that govern it. In this paper, we investigate the effect of collisions on the evolution of the spin distribution of asteroids and KBO's.
We find that the power law nature of the impactors' size distribution leads to a Lévy distribution of the spin rates. This results in a power law tail of the spin distribution, in stark contrast to the usually quoted Maxwellian distribution. We show that for bodies larger than 10 km, collisions alone lead to spin rates peaking at 0.15-0.5 revolutions per day. Comparing that to the observed spin rates of large asteroids ($R>50$ km), we find that the spins of large asteroids, peaking at $\sim1-2$ revolutions per day, are dominated by a primordial component that reflects the formation mechanism of the asteroids. Similarly, the Kuiper Belt has undergone virtually no collisional spin evolution, assuming current day density. Collisions contribute a spin rate of $\sim0.01$ revolutions per day, thus the observed fast spin rates of KBOs are also primordial in nature.
△ Less
Submitted 19 February, 2015; v1 submitted 22 June, 2014;
originally announced June 2014.
-
Sterile and Fertile Planetary Systems - Statistical Analysis of Multi-Planet Systems in Kepler's data
Authors:
Amir Weissbein,
Elad Steinberg,
Re'em Sari
Abstract:
The Kepler mission has discovered a large number of planetary systems. We analyze the implications of the discovered single/multi-exoplanet systems from Kepler's data.
As done in previous works, we test a simple model in which the intrinsic occurrence of planet is an independent process, and with equal probability around all planet producing stars. This leads to a Poisson distribution for the in…
▽ More
The Kepler mission has discovered a large number of planetary systems. We analyze the implications of the discovered single/multi-exoplanet systems from Kepler's data.
As done in previous works, we test a simple model in which the intrinsic occurrence of planet is an independent process, and with equal probability around all planet producing stars. This leads to a Poisson distribution for the intrinsic number of planets around each host. However, the possibility of zero/low mutual inclination is taken into account, creating a correlation between detecting different planets in a given stellar system, leading to a non Poisson distribution for the number of transiting planets per system. Comparing the model's predictions with the observations made by Kepler, we find that the correlation produced by planarity is insufficient and a higher correlation is needed; either the formation of one planet in the system enhances the likelihood of other planets to form, and/or that some stars are considerably more fertile than others. Kepler's data presents evidences that both correlations might play a part, in particular a significant dependency in the radial distribution of planets in multi-planet systems is shown. Followup observations on Kepler planet's hosts can help pinpoint the physical nature of this correlation.
△ Less
Submitted 20 July, 2012; v1 submitted 27 March, 2012;
originally announced March 2012.
-
Binary YORP and Evolution of Binary Asteroids
Authors:
Elad Steinberg,
Re'em Sari
Abstract:
The rotation states of kilometer sized near earth asteroids are known to be affected by the YORP effect. In a related effect, Binary YORP (BYORP) the orbital properties of a binary asteroid evolves under a radiation effect mostly acting on a tidally locked secondary. The BYORP effect can alter the orbital elements in $\sim 10^{4-5}$ years for a $D_{p}=2\;km$ primary with a $D_{s}=0.4\; km$ seconda…
▽ More
The rotation states of kilometer sized near earth asteroids are known to be affected by the YORP effect. In a related effect, Binary YORP (BYORP) the orbital properties of a binary asteroid evolves under a radiation effect mostly acting on a tidally locked secondary. The BYORP effect can alter the orbital elements in $\sim 10^{4-5}$ years for a $D_{p}=2\;km$ primary with a $D_{s}=0.4\; km$ secondary at $1\; AU$. It can either separate the binary components or cause them to collide. In this paper we devise a simple approach to calculate the YORP effect on asteroids and BYORP effect on binaries including $J_2$ effects due to primary oblateness and the sun. We apply this to asteroids with known shapes as well as a set of randomly generated bodies with various degrees of smoothness. We find a strong correlation between the strengths of an asteroid's YORP and BYORP effects. Therefore, a statistical knowledge on one, could be used to estimate the effect of the other. We show that the action of BYORP preferentially shrinks rather than expands the binary orbit and that YORP preferentially slows down asteroids. This conclusion holds for the two extremes of thermal conductivities studied in this work and assuming the asteroid reaches a stable point, but may break down for moderate thermal conductivity. The YORP and BYORP effects are shown to be smaller than what could be naively expected due to near cancellation of the effects on small scales. Taking this near cancellation into account, a simple order of magnitude estimate of the YORP and BYORP effects as function of the sizes and smoothness of the bodies is calculated. Finally, we provide a simple proof showing that there is no secular effect due to absorption of radiation in BYORP.
△ Less
Submitted 21 December, 2010; v1 submitted 13 October, 2010;
originally announced October 2010.
