-
A 9-Month Hubble Space Telescope Near-UV Survey of M87. I. Light and Color Curves of 94 Novae, and a Re-determination of the Nova Rate
Authors:
Michael M. Shara,
Alec M. Lessing,
Rebekah Hounsell,
Shifra Mandel,
David Zurek,
Matthew J. Darnley,
Or Graur,
Yael Hillman,
Eileen T. Meyer,
Joanna Mikolajewska,
James D. Neill,
Dina Prialnik,
William Sparks
Abstract:
M87 has been monitored with a cadence of 5 days over a 9 month-long span through the near-ultraviolet (NUV:F275W) and optical (F606W) filters of the Wide Field Camera 3 (WFC3) of the $\textit{Hubble Space Telescope}$. This unprecedented dataset yields the NUV and optical light and color curves of 94 M87 novae, characterizing the outburst and decline properties of the largest extragalactic nova dat…
▽ More
M87 has been monitored with a cadence of 5 days over a 9 month-long span through the near-ultraviolet (NUV:F275W) and optical (F606W) filters of the Wide Field Camera 3 (WFC3) of the $\textit{Hubble Space Telescope}$. This unprecedented dataset yields the NUV and optical light and color curves of 94 M87 novae, characterizing the outburst and decline properties of the largest extragalactic nova dataset in the literature (after M31 and M81). We test and confirm nova modelers' prediction that recurrent novae cannot erupt more frequently that once every 45 days; show that there are zero rapidly recurring novae in the central $\sim$ 1/3 of M87 with recurrence times $ < $ 130 days; demonstrate that novae closely follow the K-band light of M87 to within a few arcsec of the galaxy nucleus; show that nova NUV light curves are as heterogeneous as their optical counterparts, and usually peak 5 to 30 days after visible light maximum; determine our observations' annual detection completeness to be 71 - 77\%; and measure the rate Rnova of nova eruptions in M87 as $352_{-37}^{+37}$/yr. The corresponding luminosity-specific classical nova rate for this galaxy is $7.91_{-1.20}^{+1.20}/yr/10^{10}L_\odot,_{K}$. These rates confirm that ground-based observations of extragalactic novae miss most faint, fast novae and those near the centers of galaxies. An annual M87 nova rate of 300 or more seems inescapable. A luminosity-specific nova rate of $\sim$ $7 - 10/yr/10^{10}L_\odot,_{K}$ in ${\it all}$ types of galaxies is indicated by the data available in 2023.
△ Less
Submitted 9 October, 2023; v1 submitted 29 August, 2023;
originally announced August 2023.
-
Sublimation of ices during the early evolution of Kuiper belt objects
Authors:
Adam Parhi,
Dina Prialnik
Abstract:
Kuiper belt objects, such as Arrokoth, the probable progenitors of short-period comets, formed and evolved at large heliocentric distances, where the ambient temperatures appear to be sufficiently low for preserving volatile ices. By detailed numerical simulations, we follow the long-term evolution of small bodies, composed of amorphous water ice, dust, and ices of other volatile species that are…
▽ More
Kuiper belt objects, such as Arrokoth, the probable progenitors of short-period comets, formed and evolved at large heliocentric distances, where the ambient temperatures appear to be sufficiently low for preserving volatile ices. By detailed numerical simulations, we follow the long-term evolution of small bodies, composed of amorphous water ice, dust, and ices of other volatile species that are commonly observed in comets. The heat sources are solar radiation and the decay of short-lived radionuclides. The bodies are highly porous and gases released in the interior flow through the porous medium. The most volatile ices, CO and CH$_4$ , are found to be depleted down to the center over a time scale on the order of 100 Myr. Sublimation fronts advance from the surface inward, and when the temperature in the inner part rises sufficiently, bulk sublimation throughout the interior reduces gradually the volatile ices content until they are completely lost. All the other ices survive, which is compatible with data collected by New Horizons on Arrokoth, showing the presence of methanol, and possibly, H$_2$O, CO$_2$, NH$_3$ and C$_2$H$_6$, but no hypervolatiles. The effect of short-lived radionuclides is to increase the sublimation equilibrium temperatures and reduce volatile depletion times. We consider the effect of the bulk density, abundance ratios and heliocentric distance. At 100~au, CO is depleted, but CH$_4$ survives to present time, except for a thin outer layer. Since CO is abundantly detected in comets, we conclude that the source of highly volatile species in active comets must be gas trapped in amorphous ice.
△ Less
Submitted 10 April, 2023;
originally announced April 2023.
-
About the loss of a primordial atmosphere of super-Earths by planetesimal impacts
Authors:
Michael Lozovsky,
Dina Prialnik,
Morris Podolak
Abstract:
We consider planets composed of water ice and rock, located far from a central star. In an earlier study, computing the growth of planets by continuous accretion, we found that a large fraction of the ice evaporates upon accretion, creating a water vapor atmosphere. Here we consider accretion as a discrete series of planetesimal impacts (of order $10^8$), at the same time-dependent accretion rate,…
▽ More
We consider planets composed of water ice and rock, located far from a central star. In an earlier study, computing the growth of planets by continuous accretion, we found that a large fraction of the ice evaporates upon accretion, creating a water vapor atmosphere. Here we consider accretion as a discrete series of planetesimal impacts (of order $10^8$), at the same time-dependent accretion rate, and investigate the fate of the vapor, as a result of its interaction with the accreting planetesimals. We find that a large fraction of the vapor escapes. The remaining fraction may form an outer layer of ice after the termination of accretion and cooling of the surface. The escaped water mass may significantly alter the ice-to-rock ratio of the planet. We investigate the effect of different choices of parameters such as the ice-to-rock ratio, the planetesimal size distribution, and the impact velocities. We find that the planetesimal size distribution has a negligible effect and explain why. By contrast, the ice-to-rock ratio and impact velocities affect the fraction of retained water masses considerably.
△ Less
Submitted 12 February, 2023;
originally announced February 2023.
-
29P/Schwassmann-Wachmann: A Rosetta Stone for Amorphous Water Ice and CO <-> CO2 Conversion in Centaurs and Comets?
Authors:
C. M. Lisse,
J. K. Steckloff,
D. Prialnik,
M. Womack,
O. Harrington-Pinto,
G. Sarid,
Y. R. Fernandez,
C. A. Schambeau,
T. Kareta,
N. H. Samarasinha,
W. Harris,
K. Volk,
L. M. Woodney,
D. P. Cruikshank,
S. A. Sandford
Abstract:
Centaur 29P/Schwassmann-Wachmann 1 (SW1) is a highly active object orbiting in the transitional Gateway region (Sarid et al. 2019) between the Centaur and Jupiter Family Comet regions. SW1 is unique among the Centaurs in that it experiences quasi-regular major outbursts and produces CO emission continuously; however, the source of the CO is unclear. We argue that due to its very large size (approx…
▽ More
Centaur 29P/Schwassmann-Wachmann 1 (SW1) is a highly active object orbiting in the transitional Gateway region (Sarid et al. 2019) between the Centaur and Jupiter Family Comet regions. SW1 is unique among the Centaurs in that it experiences quasi-regular major outbursts and produces CO emission continuously; however, the source of the CO is unclear. We argue that due to its very large size (approx. 32 km radius), SW1 is likely still responding, via amorphous water ice (AWI) conversion to crystalline water ice (CWI), to the rapid change in its external thermal environment produced by its dynamical migration from the Kuiper belt to the Gateway Region at the inner edge of the Centaur region at 6 au. It is this conversion process that is the source of the abundant CO and dust released from the object during its quiescent and outburst phases. If correct, these arguments have a number of important predictions testable via remote sensing and in situ spacecraft characterization, including: the quick release on Myr timescales of CO from AWI conversion for any few km-scale scattered disk KBO transiting into the inner system; that to date SW1 has only converted between 50 to 65% of its nuclear AWI to CWI; that volume changes upon AWI conversion could have caused subsidence and cave-ins, but not significant mass wasting or crater loss on SW1; that SW1s coma should contain abundant amounts of CWI CO2-rich icy dust particles; and that when SW1 transits into the inner system within the next 10,000 years, it will be a very different kind of JFC comet.
△ Less
Submitted 27 October, 2022; v1 submitted 19 September, 2022;
originally announced September 2022.
-
Amorphous ice in comets: evidence and consequences
Authors:
Dina Prialnik,
David Jewitt
Abstract:
Ice naturally forms in the disordered or ``amorphous'' state when accreted from vapor at temperatures and pressures found in the interstellar medium and in the frigid, low density outer regions of the Sun's protoplanetary disk. It is therefore the expected form of ice in comets and other primitive bodies that have escaped substantial heating since formation. Despite expectations, however, the obse…
▽ More
Ice naturally forms in the disordered or ``amorphous'' state when accreted from vapor at temperatures and pressures found in the interstellar medium and in the frigid, low density outer regions of the Sun's protoplanetary disk. It is therefore the expected form of ice in comets and other primitive bodies that have escaped substantial heating since formation. Despite expectations, however, the observational evidence for amorphous ice in comets remains largely indirect. This is both because the spectral features of amorphous ice are subtle and because the solar system objects for which we possess high quality data are mostly too close to the Sun and too hot for amorphous ice to survive near the surface, where it can be detected. This chapter reviews the properties of amorphous ice, the evidence for its existence and its consequences for the behavior of comets.
△ Less
Submitted 13 September, 2022;
originally announced September 2022.
-
The Effect of Accretion Rate and Composition on the Structure of Ice-rich Super-Earths
Authors:
Michael Lozovsky,
Dina Prialnik,
Morris Podolak
Abstract:
It is reasonable to assume that the structure of a planet and the interior distribution of its components are determined by its formation history. We thus follow the growth of a planet from a small embryo through its subsequent evolution. We estimate the accretion rate range based on a protoplanetary disk model at a large enough distance from the central star, for water ice to be a major component…
▽ More
It is reasonable to assume that the structure of a planet and the interior distribution of its components are determined by its formation history. We thus follow the growth of a planet from a small embryo through its subsequent evolution. We estimate the accretion rate range based on a protoplanetary disk model at a large enough distance from the central star, for water ice to be a major component. We assume the accreted material to be a mixture of silicate rock and ice, with no H-He envelope, as the accretion timescale is much longer than the time required for the nebular gas to dissipate. We adopt a thermal evolution model that includes accretional heating, radioactive energy release, and separation of ice and rock. Taking the Safronov parameter and the ice-to-rock ratio as free parameters, we compute growth and evolutionary sequences for different parameter combinations, for 4.6 Gyr.
We find the final structure to depend significantly on both parameters. Low initial ice to rock ratios and high accretion rates, each resulting in increased heating rate, lead to the formation of extended rocky cores, while the opposite conditions leave the composition almost unchanged and result in relatively low internal temperatures. When rocky cores form, the ice-rich outer mantles still contain rock mixed with the ice. We find that a considerable fraction of the ice evaporates upon accretion, depending on parameters, and assume it is lost, thus the final surface composition and bulk density of the planet do not necessarily reflect the protoplanetary disk composition.
△ Less
Submitted 26 June, 2022;
originally announced June 2022.
-
On the structure and long-term evolution of ice-rich bodies
Authors:
Stephan Loveless,
Dina Prialnik,
Morris Podolak
Abstract:
The interest in the structure of ice-rich planetary bodies, in particular the differentiation between ice and rock, has grown due to the discovery of Kuiper belt objects and exoplanets. We thus carry out a parameter study for a range of planetary masses $M$, yielding radii $50 \aplt R \aplt 3000$~km, and for rock/ice mass ratios between 0.25 and 4, evolving them for 4.5~Gyr in a cold environment,…
▽ More
The interest in the structure of ice-rich planetary bodies, in particular the differentiation between ice and rock, has grown due to the discovery of Kuiper belt objects and exoplanets. We thus carry out a parameter study for a range of planetary masses $M$, yielding radii $50 \aplt R \aplt 3000$~km, and for rock/ice mass ratios between 0.25 and 4, evolving them for 4.5~Gyr in a cold environment, to obtain the present structure. We use a thermal evolution model that allows for liquid and vapor flow in a porous medium, solving mass and energy conservation equations under hydrostatic equilibrium for a spherical body in orbit around a central star. The model includes the effect of pressure on porosity and on the melting temperature, heating by long-lived radioactive isotopes, and temperature-dependent serpentinization and dehydration. We obtain the boundary in parameter space [size, rock-content] between bodies that differentiate, forming a rocky core, and those which remain undifferentiated: small bodies, bodies with a low rock content, and the largest bodies considered, which develop high internal pressures and barely attain the melting temperature. The final differentiated structure comprises a rocky core, an ice-rich mantle, and a thin dense crust below the surface. We obtain and discuss the bulk density-radius relationship. The effect of a very cold environment is investigated and we find that at an ambient temperature of $\sim$20~K, small bodies preserve the ice in amorphous form to the present.
△ Less
Submitted 9 January, 2022;
originally announced January 2022.
-
A unified theory of cataclysmic variables from self-consistent numerical simulations
Authors:
Yael Hillman,
Michael M. Shara,
Dina Prialnik,
Attay Kovetz
Abstract:
The hydrogen-rich envelopes accreted by white dwarf stars from their red dwarf companions lead to thermonuclear runaways observed as classical nova eruptions peaking at up to 1 Million solar luminosities. Virtually all nova progenitors are novalike binaries exhibiting high rates of mass transfer to their white dwarfs before and after an eruption. It is a puzzle that binaries indistinguishable from…
▽ More
The hydrogen-rich envelopes accreted by white dwarf stars from their red dwarf companions lead to thermonuclear runaways observed as classical nova eruptions peaking at up to 1 Million solar luminosities. Virtually all nova progenitors are novalike binaries exhibiting high rates of mass transfer to their white dwarfs before and after an eruption. It is a puzzle that binaries indistinguishable from novalikes, but with much lower mass transfer rates, and resulting dwarf nova outbursts, co-exist at the same orbital periods. Nova shells surrounding several dwarf novae demonstrate that at least some novae become dwarf novae between successive nova eruptions, though the mechanisms and timescales governing mass transfer rate variations are poorly understood. Here we report simulations of the multiGyr evolution of novae which self-consistently model every eruption's thermonuclear runaway, mass and angular momentum losses, feedback due to irradiation and variable mass transfer, and orbital size and period changes. The simulations reproduce the observed wide range of mass transfer rates at a given orbital period, with large and cyclic changes in white dwarf-red dwarf binaries emerging on kyr to Myr timescales. They also demonstrate that deep hibernation, (complete stoppage of mass transfer for long periods), occurs only in short-period binaries; that initially very different binaries converge to become nearly identical systems; that while almost all prenovae should be novalike binaries, dwarf novae should also occasionally be observed to give rise to novae; and that the masses of white dwarfs decrease only slightly while their red dwarf companions are consumed.
△ Less
Submitted 14 January, 2020;
originally announced January 2020.
-
The supersoft X-ray transient ASASSN-16oh as a thermonuclear runaway without mass ejection
Authors:
Yael Hillman,
Marina Orio,
Dina Prialnik,
Michael Shara,
Pavol Bezak,
Andrej Dobrotka
Abstract:
The supersoft X-ray and optical transient ASASSN-16oh has been interpreted by Maccarone et al. (2019) as having being induced by an accretion event on a massive white dwarf, resembling a dwarf nova super-outburst. These authors argued that the supersoft X-ray spectrum had a different origin than in an atmosphere heated by shell nuclear burning, because no mass was ejected. We find instead that the…
▽ More
The supersoft X-ray and optical transient ASASSN-16oh has been interpreted by Maccarone et al. (2019) as having being induced by an accretion event on a massive white dwarf, resembling a dwarf nova super-outburst. These authors argued that the supersoft X-ray spectrum had a different origin than in an atmosphere heated by shell nuclear burning, because no mass was ejected. We find instead that the event's timescale and other characteristics are typical of non-mass ejecting thermonuclear runaways, as already predicted by Shara et al. (1977) and the extensive grid of nova models by Yaron et al. (2005). We suggest that the low X-ray and bolometric luminosity in comparison to the predictions of the models of nuclear burning are due to an optically thick accretion disk, hiding most of the white dwarf surface. If this is the case, we calculated that the optical transient can be explained as a non-ejective thermonuclear event on a WD of $\simeq$1.1M$_\odot$ accreting at the rate of $\simeq3.5{-}5{\times}10^{-7}$M$_\odot$yr$^{-1}$. We make predictions that should prove whether the nature of the transient event was due to thermonuclear burning or to accretion; observational proof should be obtained in the next few years, because a new outburst should occur within $\simeq$10-15 years of the event.
△ Less
Submitted 30 June, 2019; v1 submitted 27 June, 2019;
originally announced June 2019.
-
The Masses and Accretion Rates of White Dwarfs in Classical and Recurrent Novae
Authors:
Michael M. Shara,
Dina Prialnik,
Yael Hillman,
Attay Kovetz
Abstract:
Models have long predicted that the frequency-averaged masses of white dwarfs in Galactic classical novae are twice as large as those of field white dwarfs. Only a handful of dynamically well-determined nova white dwarf masses have been published, leaving the theoretical predictions poorly tested. The recurrence time distributions and mass accretion rate distributions of novae are even more poorly…
▽ More
Models have long predicted that the frequency-averaged masses of white dwarfs in Galactic classical novae are twice as large as those of field white dwarfs. Only a handful of dynamically well-determined nova white dwarf masses have been published, leaving the theoretical predictions poorly tested. The recurrence time distributions and mass accretion rate distributions of novae are even more poorly known. To address these deficiencies, we have combined our extensive simulations of nova eruptions with the Strope et al (2010) and Schaefer et al (2010) databases of outburst characteristics of Galactic classical and recurrent novae to determine the masses of 92 white dwarfs in novae. We find that the mean mass (frequency averaged mean mass) of 82 Galactic classical novae is 1.06 (1.13) Msun, while the mean mass of 10 recurrent novae is 1.31 Msun. These masses, and the observed nova outburst amplitude and decline time distributions allow us to determine the long-term mass accretion rate distribution of classical novae. Remarkably, that value is just 1.3 x 10^{-10} Msun/yr, which is an order of magnitude smaller than that of cataclysmic binaries in the decades before and after classical nova eruptions. This predicts that old novae become low mass transfer rate systems, and hence dwarf novae, for most of the time between nova eruptions. We determine the mass accretion rates of each of the 10 known Galactic RN, finding them to be in the range 10^{-7} - 10^{-8} $ Msun/yr. We are able to predict the recurrence time distribution of novae and compare it with the predictions of population synthesis models.
△ Less
Submitted 23 April, 2018; v1 submitted 18 April, 2018;
originally announced April 2018.
-
The Castalia Mission to Main Belt Comet 133P/Elst-Pizarro
Authors:
C. Snodgrass,
G. H. Jones,
H. Boehnhardt,
A. Gibbings,
M. Homeister,
N. Andre,
P. Beck,
M. S. Bentley,
I. Bertini,
N. Bowles,
M. T. Capria,
C. Carr,
M. Ceriotti,
A. J. Coates,
V. Della Corte,
K. L. Donaldson Hanna,
A. Fitzsimmons,
P. J. Gutierrez,
O. R. Hainaut,
A. Herique,
M. Hilchenbach,
H. H. Hsieh,
E. Jehin,
O. Karatekin,
W. Kofman
, et al. (19 additional authors not shown)
Abstract:
We describe Castalia, a proposed mission to rendezvous with a Main Belt Comet (MBC), 133P/Elst-Pizarro. MBCs are a recently discovered population of apparently icy bodies within the main asteroid belt between Mars and Jupiter, which may represent the remnants of the population which supplied the early Earth with water. Castalia will perform the first exploration of this population by characterisin…
▽ More
We describe Castalia, a proposed mission to rendezvous with a Main Belt Comet (MBC), 133P/Elst-Pizarro. MBCs are a recently discovered population of apparently icy bodies within the main asteroid belt between Mars and Jupiter, which may represent the remnants of the population which supplied the early Earth with water. Castalia will perform the first exploration of this population by characterising 133P in detail, solving the puzzle of the MBC's activity, and making the first in situ measurements of water in the asteroid belt. In many ways a successor to ESA's highly successful Rosetta mission, Castalia will allow direct comparison between very different classes of comet, including measuring critical isotope ratios, plasma and dust properties. It will also feature the first radar system to visit a minor body, mapping the ice in the interior. Castalia was proposed, in slightly different versions, to the ESA M4 and M5 calls within the Cosmic Vision programme. We describe the science motivation for the mission, the measurements required to achieve the scientific goals, and the proposed instrument payload and spacecraft to achieve these.
△ Less
Submitted 11 September, 2017;
originally announced September 2017.
-
An Extensive Grid of Models Producing Extreme Horizontal Branch Stars
Authors:
Ofer Yaron,
Dina Prialnik,
Attay Kovetz,
Michael M. Shara
Abstract:
Horizontal branch (HB) morphology is a complex multiple-parameter problem. Besides the metallicity, two other leading parameters are the mass loss rate (MLR) and the initial He abundance of the HB progenitors. Using the STAREV stellar evolution code, we produce a wide array of Extreme Horizontal Branch (EHB) stars and also examine their post-HB evolution. EHB stars are produced in our calculations…
▽ More
Horizontal branch (HB) morphology is a complex multiple-parameter problem. Besides the metallicity, two other leading parameters are the mass loss rate (MLR) and the initial He abundance of the HB progenitors. Using the STAREV stellar evolution code, we produce a wide array of Extreme Horizontal Branch (EHB) stars and also examine their post-HB evolution. EHB stars are produced in our calculations by the so called `delayed (late) hot core flash' scenario. The MLR is increased on the red giant branch (RGB) to the extent that, prior to reaching core flash conditions, only a very thin H-rich envelope remains and helium ignition takes place at hotter positions on the HRD. We perform an extensive, self-consistent parameter study, covering populations I and II (Z=0.0001-0.03), for both normal initial helium abundances and He-enriched models (up to Y=0.40). For completeness of the study and in order to point to complete trends, we chose NOT to cut out several combinations (or results of) that may extend beyond realistic limits. We present results and complete evolutionary tracks for the covered parameter space, showing in particular that: a) Increased He abundance ON ITS OWN -- without having a significant-enough MLR on the RGB -- DOES NOT lead to the production of EHB stars; however, b) The bluest (hottest) HB positions do result from the COMBINED effect of He-enhancement and increased MLR; c) The general trend is that the effective temperature on the HB increases with decreasing metallicity, but there is an indication for a halt, or even a reversal of this trend, as Z further decreases below $10^{-3}$; d) EHB stars can serve as major contributors to the UV flux emanating from their host system. Thus, the present comprehensive study both complements and lends support to previous, more restricted studies of the HB phase, and adds results for unexplored regions of the parameter space.
△ Less
Submitted 7 September, 2017;
originally announced September 2017.
-
H-Atmospheres of Icy Super-Earths Formed in situ in the Outer Solar System: An Application to a Possible Planet Nine
Authors:
Amit Levi,
Scott J. Kenyon,
Morris Podolak,
Dina Prialnik
Abstract:
We examine the possibility that icy super-Earth mass planets, formed over long time scales (0.1--1~Gyr) at large distances ($\sim$ 200--1000~AU) from their host stars, will develop massive H-rich atmospheres. Within the interior of these planets, high pressure converts CH$_4$ into ethane, butane, or diamond and releases H$_2$. Using simplified models which capture the basic physics of the internal…
▽ More
We examine the possibility that icy super-Earth mass planets, formed over long time scales (0.1--1~Gyr) at large distances ($\sim$ 200--1000~AU) from their host stars, will develop massive H-rich atmospheres. Within the interior of these planets, high pressure converts CH$_4$ into ethane, butane, or diamond and releases H$_2$. Using simplified models which capture the basic physics of the internal structure, we show that the physical properties of the atmosphere depend on the outflux of H$_2$ from the mantle. When this outflux is $\lesssim 10^{10}$ [molec cm$^{-2}$ s$^{-1}$], the outgassed atmosphere has base pressure $\lesssim$ 1 bar. Larger outflows result in a substantial atmosphere where the base pressure may approach $10^3 - 10^4$ bar. For any pressure, the mean density of these planets, 2.4--3 [g cm$^{-3}$], is much larger than the mean density of Uranus and Neptune, 1.3--1.6 [g cm$^{-3}$]. Thus, observations can distinguish between a Planet Nine with a primordial H/He-rich atmosphere accreted from the protosolar nebula and one with an atmosphere outgassed from the core.
△ Less
Submitted 2 April, 2017;
originally announced April 2017.
-
A Hubble Space Telescope survey for novae in M87. III. Are novae good standard candles 15 days after maximum brightness?
Authors:
Michael M. Shara,
Trisha F. Doyle,
Ashley Pagnotta,
James T. Garland,
Tod R. Lauer,
David Zurek,
Edward A. Baltz,
Ariel Goerl,
Attay Kovetz,
Tamara Machac,
Juan Madrid,
Joanna Mikolajewska,
J. D. Neill,
Dina Prialnik,
Doug L. Welch,
Ofer Yaron
Abstract:
Ten weeks of daily imaging of the giant elliptical galaxy M87 with the Hubble Space Telescope (HST) has yielded 41 nova light curves of unprecedented quality for extragalactic cataclysmic variables. We have recently used these light curves to demonstrate that the observational scatter in the so-called Maximum-Magnitude Rate of Decline (MMRD) relation for classical novae is so large as to render th…
▽ More
Ten weeks of daily imaging of the giant elliptical galaxy M87 with the Hubble Space Telescope (HST) has yielded 41 nova light curves of unprecedented quality for extragalactic cataclysmic variables. We have recently used these light curves to demonstrate that the observational scatter in the so-called Maximum-Magnitude Rate of Decline (MMRD) relation for classical novae is so large as to render the nova-MMRD useless as a standard candle. Here we demonstrate that a modified Buscombe - de Vaucouleurs hypothesis, namely that novae with decline times t2 > 10 days converge to nearly the same absolute magnitude about two weeks after maximum light in a giant elliptical galaxy, is supported by our M87 nova data. For 13 novae with daily-sampled light curves, well determined times of maximum light in both the F606W and F814W filters, and decline times $t2 > 10 days we find that M87 novae display M(606W,15) = -6.37 +/- 0.46 and M(814W,15) = -6.11 +/- 0.43. If very fast novae with decline times t2 < 10 days are excluded, the distances to novae in elliptical galaxies with stellar binary populations similar to those of M87 should be determinable with 1 sigma accuracies of +/-20% with the above calibrations.
△ Less
Submitted 15 November, 2017; v1 submitted 22 February, 2017;
originally announced February 2017.
-
A Hubble Space Telescope Survey for Novae in M87. II. Snuffing out the Maximum Magnitude - Rate of Decline Relation for Novae as a Non-Standard Candle, and a Prediction of the Existence of Ultrafast Novae
Authors:
Michael M. Shara,
Trisha Doyle,
Tod R. Lauer,
David Zurek,
Edward A. Baltz,
Attay Kovetz,
Juan P. Madrid,
Joanna Mikolajewska,
J. D. Neill,
Dina Prialnik,
Doug L. Welch,
Ofer Yaron
Abstract:
The extensive grid of numerical simulations of nova eruptions of Yaron et al.(2005) first predicted that some classical novae might deviate significantly from the Maximum Magnitude - Rate of Decline (MMRD) relation, which purports to characterise novae as standard candles. Kasliwal et al. (2011) have announced the observational detection of an apparently new class of faint, fast classical novae in…
▽ More
The extensive grid of numerical simulations of nova eruptions of Yaron et al.(2005) first predicted that some classical novae might deviate significantly from the Maximum Magnitude - Rate of Decline (MMRD) relation, which purports to characterise novae as standard candles. Kasliwal et al. (2011) have announced the observational detection of an apparently new class of faint, fast classical novae in the Andromeda galaxy. These objects deviate strongly from the MMRD relationship, exactly as predicted by Yaron et al. (2005). Shara et al. (2016) recently reported the first detections of faint, fast novae in M87. These previously overlooked objects are as common in the giant elliptical galaxy M87 as they are in the giant spiral M31; they comprise about 40% of all classical nova eruptions and greatly increase the observational scatter in the MMRD relation. We use the extensive grid of nova simulations of Yaron et al. (2005) to identify the underlying causes of the existence of faint, fast novae. These are systems which have accreted, and can thus eject, only very low mass envelopes, of order 10^-7 - 10^-8 Msun, on massive white dwarfs. Such binaries include, but are not limited to, the recurrent novae. These same models predict the existence of ultrafast novae which display decline times t2 as short as five hours. We outline a strategy for their future detection.
△ Less
Submitted 19 February, 2017;
originally announced February 2017.
-
A 1-D evolutionary model for icy satellites, applied to Enceladus
Authors:
Uri Malamud,
Dina Prialnik
Abstract:
We develop a long-term 1-D evolution model for icy satellites that couples multiple processes: water migration and differentiation, geochemical reactions and silicate phase transitions, compaction by self-gravity, and ablation. The model further considers the following energy sources and sinks: tidal heating, radiogenic heating, geochemical energy released by serpentinization or absorbed by minera…
▽ More
We develop a long-term 1-D evolution model for icy satellites that couples multiple processes: water migration and differentiation, geochemical reactions and silicate phase transitions, compaction by self-gravity, and ablation. The model further considers the following energy sources and sinks: tidal heating, radiogenic heating, geochemical energy released by serpentinization or absorbed by mineral dehydration, gravitational energy and insolation, and heat transport by conduction, convection, and advection. We apply the model to Enceladus, by guessing the initial conditions that would render a structure compatible with present-day observations, assuming the initial structure to have been homogeneous. Assuming the satellite has been losing water continually along its evolution, we postulate that it was formed as a more massive, more icy and more porous satellite, and gradually transformed into its present day state due to sustained long-term tidal heating. We consider several initial compositions and evolution scenarios and follow the evolution for the age of the Solar System, testing the present day model results against the available observational constraints. Our model shows the present configuration to be differentiated into a pure icy mantle, several tens of km thick, overlying a rocky core, composed of dehydrated rock at the center and hydrated rock in the outer part. For Enceladus, it predicts a higher rock/ice mass ratio than previously assumed and a thinner ice mantle, compatible with recent estimates based on gravity field measurements. Although, obviously, the model cannot be used to explain local phenomena, it sheds light on the internal structure invoked in explanations of localized features and activities.
△ Less
Submitted 25 October, 2015;
originally announced October 2015.
-
Growing White Dwarfs to the Chandrasekhar Limit: The Parameter Space of the Single Degenerate SNIa Channel
Authors:
Yael Hillman,
Dina Prialnik,
Attay Kovetz,
Michael M. Shara
Abstract:
Can a white dwarf, accreting hydrogen-rich matter from a non-degenerate companion star, ever exceed the Chandrasekhar mass and explode as a type Ia supernova? We explore the range of accretion rates that allow a white dwarf (WD) to secularly grow in mass, and derive limits on the accretion rate and on the initial mass that will allow it to reach $1.4M_\odot$ --- the Chandrasekhar mass. We follow t…
▽ More
Can a white dwarf, accreting hydrogen-rich matter from a non-degenerate companion star, ever exceed the Chandrasekhar mass and explode as a type Ia supernova? We explore the range of accretion rates that allow a white dwarf (WD) to secularly grow in mass, and derive limits on the accretion rate and on the initial mass that will allow it to reach $1.4M_\odot$ --- the Chandrasekhar mass. We follow the evolution through a long series of hydrogen flashes, during which a thick helium shell accumulates. This determines the effective helium mass accretion rate for long-term, self-consistent evolutionary runs with helium flashes. We find that net mass accumulation always occurs despite helium flashes. Although the amount of mass lost during the first few helium shell flashes is a significant fraction of that accumulated prior to the flash, that fraction decreases with repeated helium shell flashes. Eventually no mass is ejected at all during subsequent flashes. This unexpected result occurs because of continual heating of the WD interior by the helium shell flashes near its surface. The effect of heating is to lower the electron degeneracy throughout the WD, and especially in the outer layers. This key result yields helium burning that is quasi-steady state, instead of explosive. We thus find a remarkably large parameter space within which long-term, self-consistent simulations show that a WD can grow in mass and reach the Chandrasekhar limit, despite its helium flashes.
△ Less
Submitted 13 August, 2015;
originally announced August 2015.
-
HST Images Flash Ionization of Old Ejecta by the 2011 Eruption of Recurrent Nova T Pyxidis
Authors:
Michael M. Shara,
David Zurek,
Bradley E. Schaefer,
Howard E. Bond,
Patrick Godon,
Mordecai-Mark Mac Low,
Ashley Pagnotta,
Dina Prialnik,
Edward M. Sion,
Jayashree Toraskar,
Robert E. Williams
Abstract:
T Pyxidis is the only recurrent nova surrounded by knots of material ejected in previous outbursts. Following the eruption that began on 2011 April 14.29, we obtained seven epochs (from 4 to 383 days after eruption) of Hubble Space Telescope narrowband Ha images of T Pyx . The flash of radiation from the nova event had no effect on the ejecta until at least 55 days after the eruption began. Photoi…
▽ More
T Pyxidis is the only recurrent nova surrounded by knots of material ejected in previous outbursts. Following the eruption that began on 2011 April 14.29, we obtained seven epochs (from 4 to 383 days after eruption) of Hubble Space Telescope narrowband Ha images of T Pyx . The flash of radiation from the nova event had no effect on the ejecta until at least 55 days after the eruption began. Photoionization of hydrogen located north and south of the central star was seen 132 days after the beginning of the eruption. That hydrogen recombined in the following 51 days, allowing us to determine a hydrogen atom density of at least 7e5 cm^-3 - at least an order of magnitude denser than the previously detected, unresolved [NII] knots surrounding T Pyx. Material to the northwest and southeast was photoionized between 132 and 183 days after the eruption began. 99 days later that hydrogen had recombined. Both then (282 days after outburst) and 101 days later, we detected almost no trace of hydrogen emission around T Pyx. There is a large reservoir of previously unseen, cold diffuse hydrogen overlapping the previously detected, [NII] - emitting knots of T Pyx ejecta. The mass of this newly detected hydrogen is probably an order of magnitude larger than that of the [NII] knots. We also determine that there is no significant reservoir of undetected ejecta from the outer boundaries of the previously detected ejecta out to about twice that distance, near the plane of the sky. The lack of distant ejecta is consistent with the Schaefer et al (2010) scenario for T Pyx, in which the star underwent its first eruption within five years of 1866 after many millennia of quiescence, followed by the six observed recurrent nova eruptions since 1890. This lack of distant ejecta is not consistent with scenarios in which T Pyx has been erupting continuously as a recurrent nova for many centuries or millennia.
△ Less
Submitted 30 March, 2015;
originally announced March 2015.
-
Observational Signatures of SNIa Progenitors, as Predicted by Models
Authors:
Yael Hillman,
Dina Prialnik,
Attay Kovetz,
Michael M. Shara
Abstract:
A definitive determination of the progenitors of type Ia supernovae (SNIa) has been a conundrum for decades. The single degenerate scenario $-$ a white dwarf (WD) in a semi-detached binary system accreting mass from its secondary $-$ is a plausible path; however, no simulation to date has shown that such an outcome is possible. In this study, we allowed a WD with a near Chandrasekhar mass of…
▽ More
A definitive determination of the progenitors of type Ia supernovae (SNIa) has been a conundrum for decades. The single degenerate scenario $-$ a white dwarf (WD) in a semi-detached binary system accreting mass from its secondary $-$ is a plausible path; however, no simulation to date has shown that such an outcome is possible. In this study, we allowed a WD with a near Chandrasekhar mass of $1.4M_\odot$ to evolve over tens of thousands of nova cycles, accumulating mass secularly while undergoing periodic nova eruptions. We present the mass accretion limits within which a SNIa can possibly occur. The results showed, for each parameter combination within the permitted limits, tens of thousands of virtually identical nova cycles where the accreted mass exceeded the ejected mass, i.e. the WD grew slowly but steadily in mass. Finally, the WD became unstable, the maximal temperature rose by nearly two orders of magnitude, heavy element production was enhanced by orders of magnitude and the nuclear and neutrino luminosities became enormous. We also found that this mechanism leading to WD collapse is robust, with WDs in the range $1.0 - 1.38M_\odot$, and an accretion rate of $5*10^{-7}M_\odot/yr$, all growing steadily in mass. These simulations of the onset of a SNIa event make observationally testable predictions about the light curves of pre-SN stars, and about the chemistry of SNIa ejecta.
△ Less
Submitted 3 November, 2014;
originally announced November 2014.
-
Asteroid rotation periods from the Palomar Transient Factory survey
Authors:
D. Polishook,
E. O. Ofek,
A. Waszczak,
S. R. Kulkarni,
A. Gal-Yam,
O. Aharonson,
R. Laher,
J. Surace,
C. Klein,
J. Bloom,
N. Brosch,
D. Prialnik,
C. Grillmair,
S. B. Cenko,
M. Kasliwal,
N. Law,
D. Levitan,
P. Nugent,
D. Poznanski,
R. Quimby
Abstract:
The Palomar Transient Factory (PTF) is a synoptic survey designed to explore the transient and variable sky in a wide variety of cadences. We use PTF observations of fields that were observed multiple times (>=10) per night, for several nights, to find asteroids, construct their lightcurves and measure their rotation periods. Here we describe the pipeline we use to achieve these goals and present…
▽ More
The Palomar Transient Factory (PTF) is a synoptic survey designed to explore the transient and variable sky in a wide variety of cadences. We use PTF observations of fields that were observed multiple times (>=10) per night, for several nights, to find asteroids, construct their lightcurves and measure their rotation periods. Here we describe the pipeline we use to achieve these goals and present the results from the first four (overlapping) PTF fields analyzed as part of this program. These fields, which cover an area of 21 deg^2, were observed on four nights with a cadence of ~20 min. Our pipeline was able to detect 624 asteroids, of which 145 (~20%) were previously unknown. We present high quality rotation periods for 88 main-belt asteroids and possible period or lower limit on the period for an additional 85 asteroids. For the remaining 451 asteroids, we present lower limits on their photometric amplitudes. Three of the asteroids have lightcurves that are characteristic of binary asteroids. We estimate that implementing our search for all existing high-cadence PTF data will provide rotation periods for about 10,000 asteroids mainly in the magnitude range ~14 to ~20.
△ Less
Submitted 9 January, 2012;
originally announced January 2012.
-
Rotation Periods of Binary Asteroids with Large Separations - Confronting the Escaping Ejecta Binaries Model with Observations
Authors:
D. Polishook,
N. Brosch,
D. Prialnik
Abstract:
Durda et al. (2004), using numerical models, suggested that binary asteroids with large separation, called Escaping Ejecta Binaries (EEBs), can be created by fragments ejected from a disruptive impact event. It is thought that six binary asteroids recently discovered might be EEBs because of the high separation between their components (~100 > a/Rp > ~20). However, the rotation periods of four out…
▽ More
Durda et al. (2004), using numerical models, suggested that binary asteroids with large separation, called Escaping Ejecta Binaries (EEBs), can be created by fragments ejected from a disruptive impact event. It is thought that six binary asteroids recently discovered might be EEBs because of the high separation between their components (~100 > a/Rp > ~20). However, the rotation periods of four out of the six objects measured by our group and others and presented here show that these suspected EEBs have fast rotation rates of 2.5 to 4 hours. Because of the small size of the components of these binary asteroids, linked with this fast spinning, we conclude that the rotational-fission mechanism, which is a result of the thermal YORP effect, is the most likely formation scenario. Moreover, scaling the YORP effect for these objects shows that its timescale is shorter than the estimated ages of the three relevant Hirayama families hosting these binary asteroids. Therefore, only the largest (D~19 km) suspected asteroid, (317) Roxane, could be, in fact, the only known EEB. In addition, our results confirm the triple nature of (3749) Balam by measuring mutual events on its lightcurve that match the orbital period of a nearby satellite in addition to its distant companion. Measurements of (1509) Esclangona at different apparitions show a unique shape of the lightcurve that might be explained by color variations.
△ Less
Submitted 21 December, 2010;
originally announced December 2010.
-
An Extended Grid of Nova Models. III. Very Luminous, Red Novae
Authors:
Michael M. Shara,
Ofer Yaron,
Dina Prialnik,
Attay Kovetz,
David Zurek
Abstract:
Extremely luminous, red eruptive variables like RV in M31 are being suggested as exemplars of a new class of astrophysical object. Our greatly extended series of nova simulations shows that classical nova models can produce very red, luminous eruptions. In a poorly studied corner of 3-D nova parameter space (very cold, low-mass white dwarfs, accreting at very low rates) we find bona fide classical…
▽ More
Extremely luminous, red eruptive variables like RV in M31 are being suggested as exemplars of a new class of astrophysical object. Our greatly extended series of nova simulations shows that classical nova models can produce very red, luminous eruptions. In a poorly studied corner of 3-D nova parameter space (very cold, low-mass white dwarfs, accreting at very low rates) we find bona fide classical novae that are very luminous and red because they eject very slowly moving, massive envelopes. A crucial prediction of these nova models - in contrast to the predictions of merging star ("mergeburst") models - is that a hot remnant, the underlying white dwarf, will emerge after the massive ejected envelope has expanded enough to become optically thin. This blue remnant must fade on a timescale of decades - much faster than a "mergeburst", which must fade on timescales of millennia or longer. Furthermore, the cooling nova white dwarf and its expanding ejecta must become redder in the years after eruption, while a contracting mergeburst must become hotter and bluer. We predict that red novae will always brighten to L~1000 Lsun for about 1 year before rising to maximum luminosity at L~10^6 - 10^7 Lsun. The maximum luminosity attainable by a nova is likely to be L~10^7 Lsun, corresponding to M-12. In an accompanying paper we describe a fading, luminous blue candidate for the remnant of M31-RV; it is observed with HST to be compatible only with the nova model.
△ Less
Submitted 20 September, 2010;
originally announced September 2010.
-
The Red Nova-like Variable in M31 - A Blue Candidate in Quiescence
Authors:
Michael M. Shara,
David Zurek,
Dina Prialnik,
Ofer Yaron,
Attay Kovetz
Abstract:
M31-RV was an extraordinarily luminous (~10^6 Lsun) eruptive variable, displaying very cool temperatures (roughly 1000 Kelvins) as it faded. The photometric behavior of M31-RV (and several other very red novae, i.e. luminous eruptive red variables) has led to several models of this apparently new class of astrophysical object. One of the most detailed models is that of "mergebursts": hypothetical…
▽ More
M31-RV was an extraordinarily luminous (~10^6 Lsun) eruptive variable, displaying very cool temperatures (roughly 1000 Kelvins) as it faded. The photometric behavior of M31-RV (and several other very red novae, i.e. luminous eruptive red variables) has led to several models of this apparently new class of astrophysical object. One of the most detailed models is that of "mergebursts": hypothetical mergers of close binary stars. These are predicted to rival or exceed the brightest classical novae in luminosity, but to be much cooler and redder than classical novae, and to become slowly hotter and bluer as they age. This prediction suggests two stringent and definitive tests of the mergeburst hypothesis. First, there should always be a cool red remnant, and NOT a hot blue remnant at the site of such an outburst. Second, the inflated envelope of a mergeburst event should be slowly contracting, hence it must display a slowly rising effective temperature. We have located a luminous, UV-bright object within 0.4 arcsec (1.5 sigma of the astrometric position) of M31-RV in archival WFPC2 images taken 10 years after the outburst: it resembles an old nova. Twenty years after the outburst, the object remains much too hot to be a mergeburst. Its behavior remains consistent with that of theoretical nova models which erupt on a low mass white dwarf. Future Hubble UV and visible images could determine if the M31-RV analogs (in M85 and in M99) are also behaving like old novae.
△ Less
Submitted 15 September, 2010;
originally announced September 2010.
-
Non-Equipartition of Energy, Masses of Nova Ejecta, and Type Ia Supernovae
Authors:
Michael M. Shara,
Ofer Yaron,
Dina Prialnik,
Attay Kovetz
Abstract:
The total masses ejected during classical nova eruptions are needed to answer two questions with broad astrophysical implications: Can accreting white dwarfs be pushed towards the Chandrasekhar mass limit to yield type Ia supernovae? Are Ultra-luminous red variables a new kind of astrophysical phenomenon, or merely extreme classical novae? We review the methods used to determine nova ejecta mass…
▽ More
The total masses ejected during classical nova eruptions are needed to answer two questions with broad astrophysical implications: Can accreting white dwarfs be pushed towards the Chandrasekhar mass limit to yield type Ia supernovae? Are Ultra-luminous red variables a new kind of astrophysical phenomenon, or merely extreme classical novae? We review the methods used to determine nova ejecta masses. Except for the unique case of BT Mon (nova 1939), all nova ejecta mass determinations depend on untested assumptions and multi-parameter modeling. The remarkably simple assumption of equipartition between kinetic and radiated energy (E_kin and E_rad, respectively) in nova ejecta has been invoked as a way around this conundrum for the ultra-luminous red variable in M31. The deduced mass is far larger than that produced by any classical nova model. Our nova eruption simulations show that radiation and kinetic energy in nova ejecta are very far from being in energy equipartition, with variations of four orders of magnitude in the ratio E_kin/E_rad being commonplace. The assumption of equipartition must not be used to deduce nova ejecta masses; any such "determinations" can be overestimates by a factor of up to 10,000. We data-mined our extensive series of nova simulations to search for correlations that could yield nova ejecta masses. Remarkably, the mass ejected during a nova eruption is dependent only on (and is directly proportional to) E_rad. If we measure the distance to an erupting nova and its bolometric light curve then E_rad and hence the mass ejected can be directly measured.
△ Less
Submitted 11 February, 2010;
originally announced February 2010.
-
From KBOs to Centaurs: The Thermal Connection
Authors:
G. Sarid,
D. Prialnik
Abstract:
We present results of thermal evolution calculations for objects originating in the Kuiper belt and transferring inwards, to the region of the outer planets. Kuiper belt objects (KBOs) are considered to be part of a reservoir that supplies the flux of small icy bodies, mainly Centaurs and Jupiter-family comets, to regions interior to the orbit of Neptune. We study the internal thermal evolution,…
▽ More
We present results of thermal evolution calculations for objects originating in the Kuiper belt and transferring inwards, to the region of the outer planets. Kuiper belt objects (KBOs) are considered to be part of a reservoir that supplies the flux of small icy bodies, mainly Centaurs and Jupiter-family comets, to regions interior to the orbit of Neptune. We study the internal thermal evolution, for yr, of three typical KBOs and use the end state of the simulation as initial conditions for evolutionary calculations of two typical Centaurs. Some evolutionary trends can be identified for the KBOs, depending on key physical parameters, such as size and composition. The subsequent evolution in the Centaur region results in both specific features for each modeled object (mainly surface and sub-surface composition) and common characteristics of thermally evolved Centaurs.
△ Less
Submitted 14 September, 2009;
originally announced September 2009.
-
Simultaneous Spectroscopic and Photometric Observations of Binary Asteroids
Authors:
D. Polishook,
N. Brosch,
D. Prialnik,
S. Kaspi
Abstract:
We present results of visible wavelengths spectroscopic measurements (0.45 to 0.72 microns) of two binary asteroids, obtained with the 1-m telescope at the Wise Observatory on January 2008. The asteroids (90) Antiope and (1509) Esclangona were observed to search for spectroscopic variations correlated with their rotation while presenting different regions of their surface to the viewer. Simultan…
▽ More
We present results of visible wavelengths spectroscopic measurements (0.45 to 0.72 microns) of two binary asteroids, obtained with the 1-m telescope at the Wise Observatory on January 2008. The asteroids (90) Antiope and (1509) Esclangona were observed to search for spectroscopic variations correlated with their rotation while presenting different regions of their surface to the viewer. Simultaneous photometric observations were performed with the Wise Observatory's 0.46-m telescope, to investigate the rotational phase behavior and possible eclipse events. (90) Antiope displayed an eclipse event during our observations. We could not measure any slope change of the spectroscopic albedo within the error range of 3%, except for a steady decrease in the total light flux while the eclipse took place. We conclude that the surface compositions of the two components do not differ dramatically, implying a common origin and history. (1509) Esclangona did not show an eclipse, but rather a unique lightcurve with three peaks and a wide and flat minimum, repeating with a period of 3.2524 hours. Careful measurements of the spectral albedo slopes reveal a color variation of 7 to 10 percent on the surface of (1509) Esclangona, which correlates with a specific region in the photometric lightcurve. This result suggests that the different features on the lightcurve are at least partially produced by color variations and could perhaps be explained by the existence of an exposed fresh surface on (1509) Esclangona.
△ Less
Submitted 2 September, 2009;
originally announced September 2009.
-
A New, Efficient Stellar Evolution Code for Calculating Complete Evolutionary Tracks
Authors:
Attay Kovetz,
Ofer Yaron,
Dina Prialnik
Abstract:
We present a new stellar evolution code and a set of results, demonstrating its capability at calculating full evolutionary tracks for a wide range of masses and metallicities. The code is fast and efficient, and is capable of following through all evolutionary phases, without interruption or human intervention. It is meant to be used also in the context of modeling the evolution of dense stella…
▽ More
We present a new stellar evolution code and a set of results, demonstrating its capability at calculating full evolutionary tracks for a wide range of masses and metallicities. The code is fast and efficient, and is capable of following through all evolutionary phases, without interruption or human intervention. It is meant to be used also in the context of modeling the evolution of dense stellar systems, for performing live calculations for both normal star models and merger-products.
The code is based on a fully implicit, adaptive-grid numerical scheme that solves simultaneously for structure, mesh and chemical composition. Full details are given for the treatment of convection, equation of state, opacity, nuclear reactions and mass loss.
Results of evolutionary calculations are shown for a solar model that matches the characteristics of the present sun to an accuracy of better than 1%; a 1 Msun model for a wide range of metallicities; a series of models of stellar populations I and II, for the mass range 0.25 to 64 Msun, followed from pre-main-sequence to a cool white dwarf or core collapse. An initial final-mass relationship is derived and compared with previous studies. Finally, we briefly address the evolution of non-canonical configurations, merger-products of low-mass main-sequence parents.
△ Less
Submitted 5 March, 2009; v1 submitted 24 September, 2008;
originally announced September 2008.
-
A fully 3-dimensional thermal model of a comet nucleus
Authors:
E. D. Rosenberg,
D. Prialnik
Abstract:
A 3-D numerical model of comet nuclei is presented. An implicit numerical scheme was developed for the thermal evolution of a spherical nucleus composed of a mixture of ice and dust. The model was tested against analytical solutions, simplified numerical solutions, and 1-D thermal evolution codes. The 3-D code was applied to comet 67P/Churyumov-Gerasimenko; surface temperature maps and the inter…
▽ More
A 3-D numerical model of comet nuclei is presented. An implicit numerical scheme was developed for the thermal evolution of a spherical nucleus composed of a mixture of ice and dust. The model was tested against analytical solutions, simplified numerical solutions, and 1-D thermal evolution codes. The 3-D code was applied to comet 67P/Churyumov-Gerasimenko; surface temperature maps and the internal thermal structure was obtained as function of depth, longitude and hour angle. The effect of the spin axis tilt on the surface temperature distribution was studied in detail. It was found that for small tilt angles, relatively low temperatures may prevail on near-pole areas, despite lateral heat conduction. A high-resolution run for a comet model of 67P/Churyumov-Gerasimenko with low tilt angle, allowing for crystallization of amorphous ice, showed that the amorphous/crystalline ice boundary varies significantly with depth as a function of cometary latitude.
△ Less
Submitted 11 March, 2007;
originally announced March 2007.
-
A Thousand and One Nova Outbursts
Authors:
N. Epelstain,
O. Yaron,
A. Kovetz,
D. Prialnik
Abstract:
Multicycle nova evolution models have been calculated over the past twenty years, the number being limited by numerical constraints. Here we present a long-term evolution code that enables a continuous calculation through an unlimited number of nova cycles for an unlimited evolution time, even up to (or exceeding) a Hubble time. Starting with two sets of the three independent nova parameters --…
▽ More
Multicycle nova evolution models have been calculated over the past twenty years, the number being limited by numerical constraints. Here we present a long-term evolution code that enables a continuous calculation through an unlimited number of nova cycles for an unlimited evolution time, even up to (or exceeding) a Hubble time. Starting with two sets of the three independent nova parameters -- the white dwarf mass, the temperature of its isothermal core, and the rate of mass transfer on to it -- we have followed the evolution of two models, with initial masses of 1 and 0.65 solar masses, accretion rates (constant throughout each calculation) of 1e-11 and 1e-9 solar-masses/yr, and relatively high initial temperatures (as they are likely to be at the onset of the outburst phase), through over 1000 and over 3000 cycles, respectively. The results show that although on the short-term consecutive outbursts are almost identical, on the long-term scale the characteristics change. This is mainly due to the changing core temperature, which decreases very similarly to that of a cooling white dwarf for a time, but at a slower rate thereafter. As the white dwarf's mass continually decreases, since both models lose more mass than they accrete, the central pressure decreases accordingly. The outbursts on the massive white dwarf change gradually from fast to moderately fast, and the other characteristics (velocity, abundance ratios, isotopic ratios) change, too. Very slowly, a steady state is reached, where all characteristics, both in quiescence and in outburst, remain almost constant. For the less massive white dwarf accreting at a high rate, outbursts are similar throughout the evolution.
△ Less
Submitted 5 November, 2006;
originally announced November 2006.
-
Thermal evolution and activity of Comet 9P/Tempel 1 and simulation of a deep impact
Authors:
G. Sarid,
D. Prialnik,
K. J. Meech,
J. Pittichova,
T. L. Farnham
Abstract:
We use a quasi 3-D thermal evolution model for a spherical comet nucleus, which takes into account the diurnal and latitudinal variation of the solar flux, but neglects lateral heat conduction. We model the thermal evolution and activity of Comet 9P/Tempel 1, in anticipation of the Deep Impact mission encounter with the comet. We also investigate the possible outcome of a projectile impact, assu…
▽ More
We use a quasi 3-D thermal evolution model for a spherical comet nucleus, which takes into account the diurnal and latitudinal variation of the solar flux, but neglects lateral heat conduction. We model the thermal evolution and activity of Comet 9P/Tempel 1, in anticipation of the Deep Impact mission encounter with the comet. We also investigate the possible outcome of a projectile impact, assuming that all the energy is absorbed as thermal energy. An interesting result of this investigation, is that the estimated amount of dust ejected due to the impact is equivalent to 2--2.6 days of activity, during "quiet" conditions, at perihelion.
We show that production rates of volatiles that are released in the interior of the nucleus depend strongly on the porous structure, in particular on the surface to volume ratio of the pores. We develop a more accurate model for calculating this parameter, based on a distribution of pore sizes, rather than a single, average pore size.
△ Less
Submitted 1 July, 2005;
originally announced July 2005.
-
An Extended Grid of Nova Models: II. The Parameter Space of Nova Outbursts
Authors:
O. Yaron,
D. Prialnik,
M. M. Shara,
A. Kovetz
Abstract:
This paper is a sequel to an earlier paper devoted to multiple, multicycle nova evolution models (Prialnik & Kovetz 1995, first paper of the series), which showed that the different characteristics of nova outbursts can be reproduced by varying the values of three basic and independent parameters: the white dwarf mass-M_{WD}, the temperature of its isothermal core-T_{WD} and the mass transfer ra…
▽ More
This paper is a sequel to an earlier paper devoted to multiple, multicycle nova evolution models (Prialnik & Kovetz 1995, first paper of the series), which showed that the different characteristics of nova outbursts can be reproduced by varying the values of three basic and independent parameters: the white dwarf mass-M_{WD}, the temperature of its isothermal core-T_{WD} and the mass transfer rate-Mdot. Apart from being the largest computational classical novae parameter-space survey, we show here that the parameter space is constrained by several analytical considerations and find its limiting surfaces. Consequently, we extend the grid of multicycle nova evolution models presented in Paper I almost to its limits, adding multicycle nova outburst calculations for a considerable number of new parameter combinations. In particular, the extended parameter pace that produces nova eruptions includes low mass transfer rates down to 5e-13 Msun/yr, and more models for low T_{WD}. Resulting characteristics of these runs are added to the former parameter combination results, to provide a full grid spanning the entire parameter space for Carbon-Oxygen white dwarfs. The full grid covers the entire range of observed nova characteristics, even those of peculiar objects, which have not been numerically reproduced until now. Most remarkably, runs for very low Mdot lead to very high values for some characteristics, such as outburst amplitude A>~20, high super-Eddington luminosities at maximum, heavy element abundance of the ejecta Z_{ej}~0.63 and high ejected masses m_{ej}~7e-4 Msun.
△ Less
Submitted 7 March, 2005;
originally announced March 2005.
-
FAUST Observations in the Fourth Galactic Quadrant
Authors:
N. Brosch,
A. Brook,
L. Wisotzki,
E. Almoznino,
Dina Prialnik,
S. Bowyer,
M. Lampton
Abstract:
We analyze UV observations with FAUST of four sky fields in the general direction of the fourth Galactic quadrant, in which we detect 777 UV sources. This is ~50% more than detected originally by Bowyer et al (1995). We discuss the source detection process and the identification of UV sources with optical counterparts. For the first time in this project we use ground-based objective-prism inform…
▽ More
We analyze UV observations with FAUST of four sky fields in the general direction of the fourth Galactic quadrant, in which we detect 777 UV sources. This is ~50% more than detected originally by Bowyer et al (1995). We discuss the source detection process and the identification of UV sources with optical counterparts. For the first time in this project we use ground-based objective-prism information for two of the fields, to select the best-matching optical objects with which to identify the UV sources. Using this, and correlations with existing catalogs, we present reliable identifications for \~75% of the sources. Most of the remaining sources have assigned optical counterparts, but lacking additional information we offer only plausible identification. We discuss the types of objects found, and compare the observed population with predictions of our UV Galaxy model.
△ Less
Submitted 22 December, 1999;
originally announced December 1999.