-
Young asteroid families as the primary source of meteorites
Authors: M. Brož, P. Vernazza, M. Marsset, F. E. DeMeo, R. P. Binzel, D. Vokrouhlický, D. Nesvorný
Abstract: Understanding the origin of bright shooting stars and their meteorite samples is among the most ancient astronomy-related questions that at larger scales has human consequences [1-3]. As of today, only ${\sim}\,6\%$ of meteorite falls have been firmly linked to their sources (Moon, Mars, and asteroid (4) Vesta [4-6]). Here, we show that ${\sim}\,70\%$ of meteorites originate from three recent brea… ▽ More Understanding the origin of bright shooting stars and their meteorite samples is among the most ancient astronomy-related questions that at larger scales has human consequences [1-3]. As of today, only ${\sim}\,6\%$ of meteorite falls have been firmly linked to their sources (Moon, Mars, and asteroid (4) Vesta [4-6]). Here, we show that ${\sim}\,70\%$ of meteorites originate from three recent breakups of $D > 30\,{\rm km}$ asteroids that occurred 5.8, 7.5 and less than ${\sim}\,40$ million years ago. These breakups, including the well-known Karin family [7], took place in the prominent yet old Koronis and Massalia families and are at the origin of the dominance of H and L ordinary chondrites among meteorite falls. These young families distinguish themselves amidst all main belt asteroids by having a uniquely high abundance of small fragments. Their size-frequency distribution remains steep for a few tens of millions of years, exceeding temporarily the production of metre-sized fragments by the largest old asteroid families (e.g., Flora, Vesta). Supporting evidence includes the existence of associated dust bands [8-10], the cosmic-ray exposure ages of H-chondrite meteorites [11,12], or the distribution of pre-atmospheric orbits of meteorites [13-15]. △ Less
Submitted 13 March, 2024; originally announced March 2024.
-
The Massalia asteroid family as the origin of ordinary L chondrites
Authors: Michaël Marsset, Pierre Vernazza, Miroslav Brož, Cristina A. Thomas, Francesca E. DeMeo, Brian Burt, Richard P. Binzel, Vishnu Reddy, Allison McGraw, Chrysa Avdellidou, Benoit Carry, Stephen M. Slivan, David Polishook
Abstract: Studies of micrometeorites in mid-Ordovician limestones and Earth's impact craters indicate that our planet witnessed a massive infall of ordinary L chondrite material 466 million years (My) ago (Heck et al. 2017, Schmieder & Kring 2020, Kenkmann 2021) that may have been at the origin of the first major mass extinction event (Schmitz et al. 2019). The breakup of a large asteroid in the main belt i… ▽ More Studies of micrometeorites in mid-Ordovician limestones and Earth's impact craters indicate that our planet witnessed a massive infall of ordinary L chondrite material 466 million years (My) ago (Heck et al. 2017, Schmieder & Kring 2020, Kenkmann 2021) that may have been at the origin of the first major mass extinction event (Schmitz et al. 2019). The breakup of a large asteroid in the main belt is the likely cause of this massive infall. In modern times, material originating from this breakup still dominates meteorite falls (>20% of all falls) (Swindle et al. 2014). Here, we provide spectroscopic observations and dynamical evidence that the Massalia collisional family is the only plausible source of this catastrophic event and of the most abundant class of meteorites falling on Earth today. It is suitably located in the inner belt, at low-inclination orbits, which corresponds to the observed distribution of L-chondrite-like near-Earth objects (NEOs) and of interplanetary dust concentrated at 1.4 degrees (Sykes 1990, Reach et al. 1997). △ Less
Submitted 13 March, 2024; originally announced March 2024.
-
Orbital and absolute magnitude distribution of Jupiter Trojans
Authors: David Vokrouhlický, David Nesvorný, Miroslav Brož, William F. Bottke, Rogerio Deienno, Carson D. Fuls, Frank C. Shelly
Abstract: Jupiter Trojans (JTs) librate about the Lagrangian stationary centers L4 and L5 associated with this planet on a typically small-eccentricity and moderate-inclination heliocentric orbits. The physical and orbital properties of JTs provide important clues about the dynamical evolution of the giant planets in the early Solar System, as well as populations of planetesimals in their source regions. He… ▽ More Jupiter Trojans (JTs) librate about the Lagrangian stationary centers L4 and L5 associated with this planet on a typically small-eccentricity and moderate-inclination heliocentric orbits. The physical and orbital properties of JTs provide important clues about the dynamical evolution of the giant planets in the early Solar System, as well as populations of planetesimals in their source regions. Here we use decade long observations from the Catalina Sky Survey (station G96) to determine the bias-corrected orbital and magnitude distributions of JTs. We distinguish the background JT population, filling smoothly the long-term stable orbital zone about L4 and L5 points, and collisional families. We find that the cumulative magnitude distribution of JTs (the background population in our case) has a steep slope for $H\leq 9$, followed with a moderately shallow slope till $H\simeq 14.5$, beyond which the distribution becomes even shallower. At $H=15$ we find a local power-law exponent $0.38\pm 0.01$. We confirm the asymmetry between the magnitude limited background populations in L4 and L5 clouds characterized by a ratio $1.45\pm 0.05$ for $H<15$. Our analysis suggests an asymmetry in the inclination distribution of JTs, with the L4 population being tighter and the L5 population being broader. We also provide a new catalog of the synthetic proper elements for JTs with an updated identification of statistically robust families (9 at L4, and 4 at L5). The previously known Ennomos family is found to consist of two, overlapping Deiphobus and Ennomos families. △ Less
Submitted 27 January, 2024; originally announced January 2024.
-
Debiased population of very young asteroid families
Authors: David Vokrouhlický, David Nesvorný, Miroslav Brož, William F. Bottke
Abstract: We use observations from the Catalina Sky Survey (CSS) to determine the bias-corrected population of small members in four very young families down to sizes equivalent to several hundred meters. Using the most recent catalog of known asteroids, we identified members from four young families for which the population has grown appreciably over recent times. A large fraction of these bodies have also… ▽ More We use observations from the Catalina Sky Survey (CSS) to determine the bias-corrected population of small members in four very young families down to sizes equivalent to several hundred meters. Using the most recent catalog of known asteroids, we identified members from four young families for which the population has grown appreciably over recent times. A large fraction of these bodies have also been detected by CSS. We used synthetic populations of asteroids, with their magnitude distribution controlled by a small number of parameters, as a template for the bias-corrected model of these families. Applying the known detection probability of the CSS observations, we could adjust these model parameters to match the observed (biased) populations in the young families. In the case of three families, Datura, Adelaide, and Rampo, we find evidence that the magnitude distribution transitions from steep to shallow slopes near $300$ to $400$ meters. Conversely, the Hobson family population may be represented by a single power-law model. The Lucascavin family has a limited population; no new members have been discovered over the past two decades. We consider a model of parent body rotational fission with the escaping secondary tidally split into two components (thereby providing three members within this family). In support of this idea, we find that no other asteroid with absolute magnitude $H\leq 18.3$ accompanies the known three members in the Lucascavin family. A similar result is found for the archetypal asteroid pair Rheinland--Kurpfalz. △ Less
Submitted 27 October, 2023; originally announced October 2023.
-
An advanced multipole model of the (130) Elektra quadruple system
Authors: M. Fuksa, M. Brož, J. Hanuš, M. Ferrais, P. Fatka, P. Vernazza
Abstract: The Ch-type asteroid (130) Elektra is orbited by three moons, making it the first quadruple system in the main asteroid belt. We aim to characterise the irregular shape of Elektra and construct a complete orbital model of its unique moon system. We applied the All-Data Asteroid Modelling (ADAM) algorithm to 60 light curves of Elektra, including our new measurements, 46 adaptive-optics (AO) ima… ▽ More The Ch-type asteroid (130) Elektra is orbited by three moons, making it the first quadruple system in the main asteroid belt. We aim to characterise the irregular shape of Elektra and construct a complete orbital model of its unique moon system. We applied the All-Data Asteroid Modelling (ADAM) algorithm to 60 light curves of Elektra, including our new measurements, 46 adaptive-optics (AO) images obtained by the VLT/SPHERE and Keck/Nirc2 instruments, and two stellar occultation profiles. For the orbital model, we used an advanced $N$-body integrator, which includes a multipole expansion of the central body (with terms up to the order $\ell = 6$), mutual perturbations, internal tides, as well as the external tide of the Sun acting on the orbits. We fitted the astrometry measured with respect to the central body and also relatively, with respect to the moons themselves. We obtained a revised shape model of Elektra with the volume-equivalent diameter $(201\pm 2)\,{\rm km}$. Out of two pole solutions, $(λ, β) = (189; -88)\,{\rm deg}$ is preferred, because the other one leads to an incorrect orbital evolution of the moons. We also identified the true orbital period of the third moon S/2014 (130) 2 as $P_2 = (1.642112 \pm 0.000400)\,{\rm d}$, which is in between the other periods, $P_1 \simeq 1.212\,{\rm d}$, $P_3 \simeq 5.300\,{\rm d}$, of S/2014 (130) 1 and S/2003 (130) 1, respectively. The resulting mass of Elektra, $(6.606 \substack{+0.007 \\ -0.013}) \times 10^{18}\,{\rm kg}$, is precisely constrained by all three orbits. Its bulk density is then $(1.536 \pm 0.038)\,{\rm g\,cm}^{-3}$. The expansion with the assumption of homogeneous interior leads to the oblateness $J_2 = -C_{20} \simeq 0.16$. However, the best-fit precession rates indicate a slightly higher value, ${\simeq}\,0.18$. △ Less
Submitted 6 September, 2023; originally announced September 2023.
-
Science with a small two-band UV-photometry mission II: Observations of stars and stellar systems
Authors: J. Krtička, J. Benáček, J. Budaj, D. Korčáková, A. Pál, M. Piecka, M. Zejda, V. Bakış, M. Brož, Hsiang-Kuang Chang, N. Faltová, R. Gális, D. Jadlovský, J. Janík, J. Kára, J. Kolář, I. Krtičková, J. Kubát, B. Kubátová, P. Kurfürst, M. Labaj, J. Merc, Z. Mikulášek, F. Münz, E. Paunzen , et al. (10 additional authors not shown)
Abstract: We outline the impact of a small two-band UV-photometry satellite mission on the field of stellar physics, magnetospheres of stars, binaries, stellar clusters, interstellar matter, and exoplanets. On specific examples of different types of stars and stellar systems, we discuss particular requirements for such satellite missions in terms of specific mission parameters such as bandpass, precision, c… ▽ More We outline the impact of a small two-band UV-photometry satellite mission on the field of stellar physics, magnetospheres of stars, binaries, stellar clusters, interstellar matter, and exoplanets. On specific examples of different types of stars and stellar systems, we discuss particular requirements for such satellite missions in terms of specific mission parameters such as bandpass, precision, cadence, and mission duration. We show that such a mission may provide crucial data not only for hot stars that emit most of their light in UV, but also for cool stars, where UV traces their activity. This is important, for instance, for exoplanetary studies, because the level of stellar activity influences habitability. While the main asset of the two-band UV mission rests in time-domain astronomy, an example of open clusters proves that such a mission would be important also for the study of stellar populations. Properties of the interstellar dust are best explored when combining optical and IR information with observations in UV. It is well known that dust absorbs UV radiation efficiently. Consequently, we outline how such a UV mission can be used to detect eclipses of sufficiently hot stars by various dusty objects and study disks, rings, clouds, disintegrating exoplanets or exoasteroids. Furthermore, UV radiation can be used to study the cooling of neutron stars providing information about the extreme states of matter in the interiors of neutron stars and used for mapping heated spots on their surfaces. △ Less
Submitted 21 February, 2024; v1 submitted 26 June, 2023; originally announced June 2023.
-
NEOMOD: A New Orbital Distribution Model for Near Earth Objects
Authors: David Nesvorny, Rogerio Deienno, William F. Bottke, Robert Jedicke, Shantanu Naidu, Steven R. Chesley, Paul W. Chodas, Mikael Granvik, David Vokrouhlicky, Miroslav Broz, Alessandro Morbidelli, Eric Christensen, Bryce T. Bolin
Abstract: Near Earth Objects (NEOs) are a transient population of small bodies with orbits near or in the terrestrial planet region. They represent a mid-stage in the dynamical cycle of asteroids and comets, which starts with their removal from the respective source regions -- the main belt and trans-Neptunian scattered disk -- and ends as bodies impact planets, disintegrate near the Sun, or are ejected fro… ▽ More Near Earth Objects (NEOs) are a transient population of small bodies with orbits near or in the terrestrial planet region. They represent a mid-stage in the dynamical cycle of asteroids and comets, which starts with their removal from the respective source regions -- the main belt and trans-Neptunian scattered disk -- and ends as bodies impact planets, disintegrate near the Sun, or are ejected from the Solar System. Here we develop a new orbital model of NEOs by numerically integrating asteroid orbits from main belt sources and calibrating the results on observations of the Catalina Sky Survey. The results imply a size-dependent sampling of the main belt with the $ν_6$ and 3:1 resonances producing $\simeq 30$\% of NEOs with absolute magnitudes $H = 15$ and $\simeq 80$\% of NEOs with $H = 25$. Hence, the large and small NEOs have different orbital distributions. The inferred flux of $H<18$ bodies into the 3:1 resonance can be sustained only if the main-belt asteroids near the resonance drift toward the resonance at the maximal Yarkovsky rate ($\simeq 2 \times 10^{-4}$ au Myr$^{-1}$ for diameter $D=1$ km and semimajor axis $a=2.5$~au). This implies obliquities $θ\simeq 0^\circ$ for $a<2.5$~au and $θ\simeq 180^\circ$ for $a>2.5$~au, both in the immediate neighborhood of the resonance (the same applies to other resonances as well). We confirm the size-dependent disruption of asteroids near the Sun found in previous studies. An interested researcher can use the publicly available NEOMOD Simulator to generate user-defined samples of NEOs from our model. △ Less
Submitted 15 June, 2023; originally announced June 2023.
-
2021 occultations and transits of Linus orbiting (22) Kalliope: I. Polygonal and `cliptracing' algorithm
Authors: M. Brož, J. Ďurech, M. Ferrais, H. -J. Lee, M. -J. Kim, D. -G. Roh, H. -S. Yim, E. Jehin, A. Burdanov, J. de Wit, P. Fatka, J. Hanuš, B. Carry
Abstract: The satellite Linus orbiting the main-belt asteroid (22) Kalliope exhibited occultation and transit events in late 2021. A photometric campaign was organized and observations were taken by the TRAPPIST-South, SPECULOOS-Artemis, OWL-Net, and BOAO telescopes, with the goal to constrain models of this system. Our dynamical model is complex, with multipoles (up to the order $\ell = 2$), internal tides… ▽ More The satellite Linus orbiting the main-belt asteroid (22) Kalliope exhibited occultation and transit events in late 2021. A photometric campaign was organized and observations were taken by the TRAPPIST-South, SPECULOOS-Artemis, OWL-Net, and BOAO telescopes, with the goal to constrain models of this system. Our dynamical model is complex, with multipoles (up to the order $\ell = 2$), internal tides, and external tides. The model was constrained by astrometry (spanning 2001--2021), occultations, adaptive-optics imaging, calibrated photometry, as well as relative photometry. Our photometric model was substantially improved. A new precise (${<}\,0.1\,{\rm mmag}$) light curve algorithm was implemented, based on polygon intersections, which are computed exactly -- by including partial eclipses and partial visibility of polygons. Moreover, we implemented a `cliptracing' algorithm, based again on polygon intersections, in which partial contributions to individual pixels are computed exactly. Both synthetic light curves and synthetic images are then very smooth. Based on our combined solution, we confirmed the size of Linus, $(28\pm 1)\,{\rm km}$. However, this solution exhibits some tension between the light curves and the PISCO speckle-interferometry dataset. In most solutions, Linus is darker than Kalliope, with the albedos $A_{\rm w} = 0.40$ vs. $0.44$. This is confirmed on deconvolved images. A~detailed revision of astrometric data allowed us to revise also the $J_2 \equiv -C_{20}$ value of Kalliope. Most importantly, a~homogeneous body is excluded. For a differentiated body, two solutions exist: low-oblateness ($C_{20} \simeq -0.12$), with a~spherical iron core, and alternatively, high-oblateness ($C_{20} \simeq -0.22$) with an elongated iron core. These correspond to the low- and high-energy collisions, respectively, studied by means of SPH simulations in our previous work. △ Less
Submitted 7 June, 2023; originally announced June 2023.
-
V994 Her: A Unique Triply Eclipsing Sextuple Star System
Authors: P. Zasche, T. Borkovits, R. Jayaraman, S. A. Rappaport, M. Brož, D. Vokrouhlický, I. B. Bíró, T. Hegedüs, Z. T. Kiss, R. Uhlař, H. M. Schwengeler, A. Pál, M. Mašek, S. B. Howell, S. Dallaporta, U. Munari, R. Gagliano, T. Jacobs, M. H. Kristiansen, D. LaCourse, M. Omohundro, I. Terentev, A. Vanderburg, Z. Henzl, B. P. Powell , et al. (1 additional authors not shown)
Abstract: We report the discovery with $TESS$ of a third set of eclipses from V994 Herculis (TIC 424508303), previously only known as a doubly-eclipsing system. The key implication of this discovery and our analyses is that V994 Her is the second fully-characterized (2+2) + 2 sextuple system, in which all three binaries eclipse. In this work, we use a combination of ground-based observations and $TESS$ data… ▽ More We report the discovery with $TESS$ of a third set of eclipses from V994 Herculis (TIC 424508303), previously only known as a doubly-eclipsing system. The key implication of this discovery and our analyses is that V994 Her is the second fully-characterized (2+2) + 2 sextuple system, in which all three binaries eclipse. In this work, we use a combination of ground-based observations and $TESS$ data to analyze the eclipses of binaries A and B in order to update the parameters of the inner quadruple's orbit (with a derived period of 1062 $\pm$ 2d). The eclipses of binary C that were detected in the $TESS$ data were also found in older ground-based observations, as well as in more recently obtained observations. The eclipse timing variations of all three pairs were studied in order to detect the mutual perturbations of their constituent stars, as well as those of the inner pairs in the (2+2) core. At the longest periods they arise from apsidal motion, which may help constraining parameters of the component stars' internal structure. We also discuss the relative proximity of the periods of binaries A and B to a 3:2 mean motion resonance. This work represents a step forward in the development of techniques to better understand and characterize multiple star systems, especially those with multiple eclipsing components. △ Less
Submitted 31 January, 2023; originally announced January 2023.
-
Spectrum of the secondary component and new orbital elements of the massive triple star Delta Ori A
Authors: A. Oplištilová, P. Mayer, P. Harmanec, M. Brož, A. Pigulski, H. Božić, P. Zasche, M. Šlechta, H. Pablo, P. A. Kołaczek-Szymański, A. F. J. Moffat, C. C. Lovekin, G. A. Wade, K. Zwintz, A. Popowicz, W. W. Weiss
Abstract: $δ… ▽ More $δ$ Orionis is the closest massive multiple stellar system and one of the brightest members of the Orion OB association. The primary (Aa1) is a unique evolved O star. In this work, we applied a two-step disentangling method to a series of spectra in the blue region (430 to 450 nm), and we detected spectral lines of the secondary (Aa2). For the first time, we were able to constrain the orbit of the tertiary (Ab) - to 55 450 d or 152 yr - using variable $γ$ velocities and new speckle interferometric measurements, which have been published in the Washington Double Star Catalogue. In addition, the Gaia DR3 parallax of the faint component (Ca+Cb) constrains the distance of the system to (381 $\pm$ 8) pc, which is just in the centre of the Orion OB1b association, at (382 $\pm$ 1) pc. Consequently, we found that the component masses according to the three-body model are 17.8, 8.5, and 8.7 M$_{\odot}$, for Aa1, Aa2, and Ab, respectively, with the uncertainties of the order of $1$ M$_{\odot}$. We used new photometry from the BRITE satellites together with astrometry, radial velocities, eclipse timings, eclipse duration, spectral line profiles, and spectral energy distribution to refine radiative properties. The components, classified as O9.5 II + B2 V + B0 IV, have radii of 13.1, 4.1, and 12.0 R$_{\odot}$, which means that $δ$ Ori A is a pre-mass-transfer object. The frequency of 0.478 cycles per day, known from the Fourier analysis of the residual light curve and X-ray observations, was identified as the rotation frequency of the tertiary. $δ$ Ori could be related to other bright stars in Orion, in particular, $ζ$ Ori, which has a similar architecture, or $\varepsilon$ Ori, which is a single supergiant, and possibly a post-mass-transfer object. △ Less
Submitted 24 January, 2023; originally announced January 2023.
-
The equilibrium shape of (65) Cybele: primordial or relic of a large impact?
Authors: M. Marsset, M. Brož, J. Vermersch, N. Rambaux, M. Ferrais, M. Viikinkoski, J. Hanuš, E. Jehin, E. Podlewska-Gaca, P. Bartczak, G. Dudzinski, B. Carry, P. Vernazza, R. Szakáts, R. Duffard, A. Jones, D. Molina, T. Santana-Ros, Z. Benkhaldoun, M. Birlan, C. Dumas, R. Fétick, T. Fusco, L. Jorda, F. Marchis , et al. (2 additional authors not shown)
Abstract: Cybele asteroids constitute an appealing reservoir of primitive material genetically linked to the outer Solar System, and the physical properties of the largest members can be readily accessed by large telescopes. We took advantage of the bright apparition of (65) Cybele in July and August 2021 to acquire high-angular-resolution images and optical light curves of the asteroid with which we aim to… ▽ More Cybele asteroids constitute an appealing reservoir of primitive material genetically linked to the outer Solar System, and the physical properties of the largest members can be readily accessed by large telescopes. We took advantage of the bright apparition of (65) Cybele in July and August 2021 to acquire high-angular-resolution images and optical light curves of the asteroid with which we aim to analyse its shape and bulk properties. 7 series of images acquired with VLT/SPHERE were combined with optical light curves to reconstruct the shape of the asteroid using the ADAM, MPCD, and SAGE algorithms. The origin of the shape was investigated by means of N-body simulations. Cybele has a volume-equivalent diameter of 263+/-3km and a bulk density of 1.55+/-0.19g.cm-3. Notably, its shape and rotation state are closely compatible with those of a Maclaurin equilibrium figure. The lack of a collisional family associated with Cybele and the higher bulk density of that body with respect to other large P-type asteroids suggest that it never experienced any large disruptive impact followed by rapid re-accumulation. This would imply that its present-day shape represents the original one. However, numerical integration of the long-term dynamical evolution of a hypothetical family shows that it is dispersed by gravitational perturbations and chaotic diffusion over Gyrs of evolution. The very close match between Cybele and an equilibrium figure opens up the possibility that D>260km small bodies from the outer Solar System all formed at equilibrium. However, we cannot rule out an old impact as the origin of the equilibrium shape. Cybele itself is found to be dynamically unstable, implying that it was recently (<1Ga) placed on its current orbit either through slow diffusion from a relatively stable orbit in the Cybele region or, less likely, from an unstable, JFC orbit in the planet-crossing region. △ Less
Submitted 1 December, 2022; originally announced December 2022.
-
Towards a consistent model of the hot quadruple system HD 93206 = QZ Carinæ: II. N-body model
Authors: M. Brož, P. Harmanec, P. Zasche, R. Catalan-Hurtado, B. N. Barlow, W. Frondorf, M. Wolf, H. Drechsel, R. Chini, A. Nasseri, J. Labadie-Bartz, G. W. Christie, W. S. G. Walker, M. Blackford, D. Blane, A. A. Henden, T. Bohlsen, H. Božić, J. Jonák
Abstract: HD 93206 is early-type massive stellar system, composed of components resolved by direct imaging (Ab, Ad, B, C, D) as well as a compact sub-system (Aa1, Aa2, Ac1, Ac2). Its geometry was already determined on the basis of extensive photometric, spectroscopic and interferometric observations. However, the fundamental absolute parameters are still not known precisely enough. We use an advanced N-body… ▽ More HD 93206 is early-type massive stellar system, composed of components resolved by direct imaging (Ab, Ad, B, C, D) as well as a compact sub-system (Aa1, Aa2, Ac1, Ac2). Its geometry was already determined on the basis of extensive photometric, spectroscopic and interferometric observations. However, the fundamental absolute parameters are still not known precisely enough. We use an advanced N-body model to account for all mutual gravitational perturbations among the four close components, and all observational data types, including: astrometry, radial velocities, eclipse timing variations, squared visibilities, closure phases, triple products, normalized spectra, and spectral-energy distribution (SED). The respective model has 38 free parameters, namely three sets of orbital elements, component masses, and their basic radiative properties ($T$, $\log g$, $v_{\rm rot}$). We revised the fundamental parameters of QZ Car as follows. For a model with the nominal extinction coefficient $R_V \equiv A_V/E(B-V) = 3.1$, the best-fit masses are $m_1 = 26.1\,M_{\rm S}$, $m_2 = 32.3\,M_{\rm S}$, $m_3 = 70.3\,M_{\rm S}$, $m_4 = 8.8\,M_{\rm S}$, with uncertainties of the order of $2\,M_{\rm S}$, and the system distance $d = (2800\pm 100)\,{\rm pc}$. In an alternative model, where we increased the weights of RV and TTV observations and relaxed the SED constraints, because extinction can be anomalous with $R_V \sim 3.4$, the distance is smaller, $d = (2450\pm 100)\,{\rm pc}$. This would correspond to that of Collinder 228 cluster. Independently, this is confirmed by dereddening of the SED, which is only then consistent with the early-type classification (O9.7Ib for Aa1, O8III for Ac1). Future modelling should also account for an accretion disk around Ac2 component. △ Less
Submitted 24 June, 2022; originally announced June 2022.
-
Discovery of an asteroid family linked to (22) Kalliope and its moon Linus
Authors: M. Brož, M. Ferrais, P. Vernazza, P. Ševeček, M. Jutzi
Abstract: According to adaptive-optics observations by Ferrais et al., (22) Kalliope is a 150-km, dense and differentiated body. Here, we interpret (22) Kalliope in the context of bodies in its surroundings. While there is a known moon Linus, with a 5:1 size ratio, no family has been reported in the literature, which is in contradiction with the existence of the moon. Using the hierarchical clustering metho… ▽ More According to adaptive-optics observations by Ferrais et al., (22) Kalliope is a 150-km, dense and differentiated body. Here, we interpret (22) Kalliope in the context of bodies in its surroundings. While there is a known moon Linus, with a 5:1 size ratio, no family has been reported in the literature, which is in contradiction with the existence of the moon. Using the hierarchical clustering method (HCM) along with physical data, we identified the Kalliope family. Previously, it was associated to (7481) San Marcello. We then used various models (N-body, Monte-Carlo, SPH) of its orbital and collisional evolution, including the break-up of the parent body, to estimate the dynamical age of the family and address its link to Linus. The best-fit age is (900+-100) My according to our collisional model, in agreement with the position of (22) Kalliope, which was modified by chaotic diffusion due to 4-1-1 three-body resonance with Jupiter and Saturn. It seems possible to create Linus and the Kalliope family at the same time, although our SPH simulations show a variety of outcomes, for both satellite size and the family size-frequency distribution. The shape of (22) Kalliope itself was most likely affected by gravitational reaccumulation of `streams', which creates characteristic hills observed on the surface. If the body was differentiated, its internal structure is surely asymmetric. △ Less
Submitted 31 May, 2022; originally announced May 2022.
-
SPH simulations of high-speed collisions between asteroids and comets
Authors: J. Rozehnal, M. Brož, D. Nesvorný, K. J. Walsh, D. D. Durda, D. C. Richardson, E. Asphaug
Abstract: We studied impact processes by means of smoothed-particle hydrodynamics (SPH) simulations. The method was applied to modeling formation of main-belt families during the cometary bombardment (either early or late, ${\sim}\,3.85\,{\rm Gy}$ ago). If asteroids were bombarded by comets, as predicted by the Nice model, hundreds of asteroid families (catastrophic disruptions of diameter… ▽ More We studied impact processes by means of smoothed-particle hydrodynamics (SPH) simulations. The method was applied to modeling formation of main-belt families during the cometary bombardment (either early or late, ${\sim}\,3.85\,{\rm Gy}$ ago). If asteroids were bombarded by comets, as predicted by the Nice model, hundreds of asteroid families (catastrophic disruptions of diameter $D \ge 100\,{\rm km}$ bodies) should have been created, but the observed number is only 20. Therefore we computed a standard set of 125 simulations of collisions between representative $D = 100\,{\rm km}$ asteroids and high-speed icy projectiles (comets). According to our results, the largest remnant mass $M_{\rm lr}$ is similar as in low-speed collisions, due to appropriate scaling with the effective strength $Q_{\rm eff}$, but the largest fragment mass $M_{\rm lf}$ exhibits systematic differences - it is typically smaller for craterings and bigger for super-catastrophic events. This trend does not, however, explain the non-existence of old families. The respective parametric relations can be used in other statistical (Monte-Carlo) models to better understand collisions between asteroidal and cometary populations. △ Less
Submitted 28 May, 2022; originally announced May 2022.
-
Towards a consistent model of the hot quadruple system HD 93206 = QZ Carinae - I. Observations and their initial analyses
Authors: P. Harmanec, P. Zasche, M. Brož, R. Catalan-Hurtado, B. N. Barlow, W. Frondorf, M. Wolf, H. Drechsel, R. Chini, A. Nasseri, A. Pigulski, J. Labadie-Bartz, G. W. Christie, W. S. G. Walker, M. Blackford, D. Blane, A. A. Henden, T. Bohlsen, H. Božić, J. Jonák
Abstract: The hot nine-component system HD 93206, which contains a gravitationally bounded eclipsing Ac1+Ac2 binary ($P=5.9987$~d) and a spectroscopic Aa1+Aa2 ($P=20.734$~d) binary can provide~important insights into the origin and evolution of massive stars. Using archival and new spectra, and a~rich collection of ground-based and space photometric observations, we carried out a detailed study of this obje… ▽ More The hot nine-component system HD 93206, which contains a gravitationally bounded eclipsing Ac1+Ac2 binary ($P=5.9987$~d) and a spectroscopic Aa1+Aa2 ($P=20.734$~d) binary can provide~important insights into the origin and evolution of massive stars. Using archival and new spectra, and a~rich collection of ground-based and space photometric observations, we carried out a detailed study of this object. We provide a much improved description of both short orbits and a good estimate of the mutual period of both binaries of about 14500~d (i.e. 40 years). For the first time, we detected weak lines of the fainter component of the 6.0~d eclipsing binary in the optical region of the spectrum, measured their radial velocities, and derived a mass ratio of $M_{\rm Ac2}/M_{\rm Ac1}=1.29$, which is the opposite of what was estimated from the International Ultraviolet explorer (IUE) spectra. We confirm that the eclipsing subsystem Ac is semi-detached and is therefore in a phase of large-scale mass transfer between its components. The Roche-lobe filling and spectroscopically brighter component Ac1 is the less massive of the two and is eclipsed in the secondary minimum. We show that the bulk of the \ha emission, so far believed to be associated with the eclipsing system, moves with the primary O9.7I component Aa1 of the 20.73~d spectroscopic binary. However, the weak emission in the higher Balmer lines seems to be associated with the accretion disc around component Ac2. We demonstrate that accurate masses and other basic physical properties including the distance of this unique system can be obtained but require a more sophisticated modelling. A~first step in this direction is presented in the accompanying Paper~II (Brož et al.). △ Less
Submitted 14 April, 2022; originally announced April 2022.
-
Observed tidal evolution of Kleopatra's outer satellite
Authors: M. Brož, J. Ďurech, B. Carry, F. Vachier, F. Marchis, J. Hanuš, L. Jorda, P. Vernazza, D. Vokrouhlický, M. Walterová, R. Behrend
Abstract: The orbit of the outer satellite Alexhelios of (216) Kleopatra is already constrained by adaptive-optics astrometry, obtained with the VLT/SPHERE instrument. However, there is also a preceding occultation event in 1980 attributed to this satellite. Hereinafter, we try to link all observations, spanning 1980--2018. We find the nominal orbit exhibits an unexplained shift by $+60^\circ$ in the true l… ▽ More The orbit of the outer satellite Alexhelios of (216) Kleopatra is already constrained by adaptive-optics astrometry, obtained with the VLT/SPHERE instrument. However, there is also a preceding occultation event in 1980 attributed to this satellite. Hereinafter, we try to link all observations, spanning 1980--2018. We find the nominal orbit exhibits an unexplained shift by $+60^\circ$ in the true longitude. Using both periodogram analysis and an $\ell = 10$ multipole model suitable for the motion of mutually interacting moons about the irregular body, we confirmed that it is not possible to adjust the respective osculating period $P_2$. Instead, we were forced to use a model with tidal dissipation (and increasing orbital periods) to explain the shift. We also analyzed light curves, spanning 1977--2021, and searched for the expected spin deceleration of Kleopatra. According to our best-fit model, the observed period rate is $\dot P_2 = (1.8\pm 0.1)\cdot 10^{-8}\,{\rm d}\,{\rm d}^{-1}$ and the corresponding time lag $Δt_2 = 42\,{\rm s}$ of tides, for the assumed value of the Love number $k_2 = 0.3$. It is the first detection of tidal evolution for moons orbiting 100-km asteroids. The corresponding dissipation factor $Q$ is comparable with other terrestrial bodies, albeit at a higher loading frequency $2|ω-n|$. We also predict a secular evolution of the inner moon, $\dot P_1 = 5.0\cdot 10^{-8}$, as well as a spin deceleration of Kleopatra, $\dot P_0 = 1.9\cdot 10^{-12}$. In alternative models, with moons captured in the 3:2 mean-motion resonance or more massive moons, the respective values of $Δt_2$ are a factor of 2--3 lower. Future astrometric observations by direct imaging or occultations should allow to distinguish between these models, which is important for the internal structure and mechanical properties of (216) Kleopatra. △ Less
Submitted 25 October, 2021; originally announced October 2021.
-
Terrestrial planet formation by torque-driven convergent migration of planetary embryos
Authors: M. Brož, O. Chrenko, D. Nesvorný, N. Dauphas
Abstract: Massive cores of the giant planets are thought to have formed in a gas disk by accretion of pebble-size particles whose accretional cross-section is enhanced by aerodynamic gas drag [1][2]. A commonly held view is that the terrestrial planet system formed later (30-200 Myr after the dispersal of the gas disk) by giant collisions of tens of roughly Mars-size protoplanets [3]. Here we propose, inste… ▽ More Massive cores of the giant planets are thought to have formed in a gas disk by accretion of pebble-size particles whose accretional cross-section is enhanced by aerodynamic gas drag [1][2]. A commonly held view is that the terrestrial planet system formed later (30-200 Myr after the dispersal of the gas disk) by giant collisions of tens of roughly Mars-size protoplanets [3]. Here we propose, instead, that the terrestrial planets formed earlier by gas-driven convergent migration of protoplanets toward $\sim\!1\,{\rm au}$ (related ref. [4] invoked a different process to concentrate planetesimals). To investigate situations in which convergent migration occurs, we developed a radiation-hydrodynamic model with realistic opacities [5][6] to determine the thermal structure of the gas and pebble disks in the terrestrial planet zone. We find that protoplanets rapidly grow by mutual collisions and pebble accretion, and gain orbital eccentricities by gravitational scattering and the hot-trail effect [7][8]. The orbital structure of the terrestrial planet system is well reproduced in our simulations, including its tight mass concentration at 0.7-1 au and the small sizes of Mercury and Mars. The early-stage protosolar disk temperature exceeds 1500 K inside 0.4 au implying that Mercury grew in a highly reducing environment, next to the evaporation lines of iron and silicates, influencing Mercury's bulk composition [9]. In a late-stage cold gas disk, accretion of icy/hydrated pebbles would contribute to Earth's water budget. △ Less
Submitted 23 September, 2021; originally announced September 2021.
-
(216) Kleopatra, a low density critically rotating M-type asteroid
Authors: F. Marchis, L. Jorda, P. Vernazza, M. Brož, J. Hanuš, M. Ferrais, F. Vachier, N. Rambaux, M. Marsset, M. Viikinkoski, E. Jehin, S. Benseguane, E. Podlewska-Gaca, B. Carry, A. Drouard, S. Fauvaud, M. Birlan, J. Berthier, P. Bartczak, C. Dumas, G. Dudzinski, J. Durech, J. Castillo-Rogez, F. Cipriani, F. Colas , et al. (14 additional authors not shown)
Abstract: Context. The recent estimates of the 3D shape of the M/Xe-type triple asteroid system (216) Kleopatra indicated a density of 5 g.cm$^{-3}$. Such a high density implies a high metal content and a low porosity which is not easy to reconcile with its peculiar dumbbell shape. Aims. Given the unprecedented angular resolution of the VLT/SPHERE/ZIMPOL camera, we aim to constrain the mass and the shape of… ▽ More Context. The recent estimates of the 3D shape of the M/Xe-type triple asteroid system (216) Kleopatra indicated a density of 5 g.cm$^{-3}$. Such a high density implies a high metal content and a low porosity which is not easy to reconcile with its peculiar dumbbell shape. Aims. Given the unprecedented angular resolution of the VLT/SPHERE/ZIMPOL camera, we aim to constrain the mass and the shape of Kleopatra with high accuracy, hence its density. Methods. We combined our new VLT/SPHERE observations of Kleopatra recorded in 2017 and 2018 with archival data, as well as lightcurve, occultation, and delay-Doppler images, to derive its 3D shape model using two different algorithms (ADAM, MPCD). Furthermore, an N-body dynamical model allowed us to retrieve the orbital elements of the two moons as explained in the accompanying paper. Results. The shape of Kleopatra is very close to an equilibrium dumbbell figure with two lobes and a thick neck. Its volume equivalent diameter (118.75$\pm$1.40) km and mass (2.97$\pm$0.32) 10$^{18}$ kg imply a bulk density of (3.38$\pm$0.50) g cm$^{-3}$. Such a low density for a supposedly metal-rich body indicates a substantial porosity within the primary. This porous structure along with its near-equilibrium shape is compatible with a formation scenario including a giant impact followed by reaccumulation. Kleopatra's current rotation period and dumbbell shape imply that it is in a critically rotating state. The low effective gravity along the equator of the body, together with the equatorial orbits of the moons and possibly rubble-pile structure, opens the possibility that the moons formed via mass shedding. Conclusions. Kleopatra is a puzzling multiple system due to the unique characteristics of the primary. It deserves particular attention in the future, with the Extremely Large Telescopes and possibly a dedicated space mission. △ Less
Submitted 16 August, 2021; originally announced August 2021.
-
The young Hobson family: Possible binary parent body and low-velocity dispersal
Authors: David Vokrouhlický, Miroslav Brož, Bojan Novaković, David Nesvorný
Abstract: Asteroid families with ages younger than $1$ Myr offer an interesting possibility of studying the outcomes of asteroid disruptions that are little modified by subsequent evolutionary processes. We analyze a very young asteroid family associated with (18777) Hobson in the central part of the main belt. We aim at (i) understanding its peculiar size distribution, and (ii) setting an upper limit on th… ▽ More Asteroid families with ages younger than $1$ Myr offer an interesting possibility of studying the outcomes of asteroid disruptions that are little modified by subsequent evolutionary processes. We analyze a very young asteroid family associated with (18777) Hobson in the central part of the main belt. We aim at (i) understanding its peculiar size distribution, and (ii) setting an upper limit on the characteristic dispersal velocity at subkilometer sizes corresponding to the smallest visible Hobson members. We identified the Hobson family using an up-to-date asteroid catalog. A significant increase in the number of its known members allowed us to study their size distribution and compare it with computer simulations of catastrophic disruptions. Backward orbital integrations of the heliocentric orbits allowed us to confirm the previously suggested age of Hobson and helped to estimate limits of the ejection speed. The Hobson family has an unusual size distribution: two nearly equal-size bodies, followed by a population of smaller asteroids, whose distribution takes a characteristic power law. There are two possibilities to explain these data. Either a canonical impact onto a single parent body, requiring fine-tuned impact conditions that have not been studied so far, or an unconventional model for the parent object of the Hobson family, namely a binary with $\simeq 7-9$ km primary and a $\simeq 2.5$ km secondary. In the latter case, the primary was disrupted, leaving behind the largest remnant (18777) Hobson and a suite of subkilometer asteroids. The second largest asteroid, (57738) 2001 UZ160, is the nearly intact satellite of the parent binary. The excellent convergence of nominal orbits of Hobson members sets an upper limit of $\simeq (10-20)$ m s$^{-1}$ for the initial dispersal velocity of the known members, which is consistent with both formation models. △ Less
Submitted 11 August, 2021; originally announced August 2021.
-
An advanced multipole model for (216) Kleopatra triple system
Authors: M. Brož, F. Marchis, L. Jorda, J. Hanuš, P. Vernazza, M. Ferrais, F. Vachier, N. Rambaux, M. Marsset, M. Viikinkoski, E. Jehin, S. Benseguane, E. Podlewska-Gaca, B. Carry, A. Drouard, S. Fauvaud, M. Birlan, J. Berthier, P. Bartczak, C. Dumas, G. Dudziński, J. Ďurech, J. Castillo-Rogez, F. Cipriani, F. Colas , et al. (15 additional authors not shown)
Abstract: To interpret adaptive-optics observations of (216) Kleopatra, we need to describe an evolution of multiple moons, orbiting an extremely irregular body and including their mutual interactions. Such orbits are generally non-Keplerian and orbital elements are not constants. Consequently, we use a modified $N$-body integrator, which was significantly extended to include the multipole expansion of the… ▽ More To interpret adaptive-optics observations of (216) Kleopatra, we need to describe an evolution of multiple moons, orbiting an extremely irregular body and including their mutual interactions. Such orbits are generally non-Keplerian and orbital elements are not constants. Consequently, we use a modified $N$-body integrator, which was significantly extended to include the multipole expansion of the gravitational field up to the order $\ell = 10$. Its convergence was verified against the `brute-force' algorithm. We computed the coefficients $C_{\ell m},S_{\!\ell m}$ for Kleopatra's shape, assuming a~constant bulk density. For solar-system applications, it was also necessary to implement a variable distance and geometry of observations. Our $χ^2$ metric then accounts for the absolute astrometry, the relative astrometry (2nd moon with respect to 1st), angular velocities, and also silhouettes, constraining the pole orientation. This allowed us to derive the orbital elements of Kleopatra's two moons. Using both archival astrometric data and new VLT/SPHERE observations (ESO LP 199.C-0074), we were able to identify the true periods of the moons, $P_1 = (1.822359\pm0.004156)\,{\rm d}$, $P_2 = (2.745820\pm0.004820)\,{\rm d}$. They orbit very close to the 3:2 mean-motion resonance, but their osculating eccentricities are too small compared to other perturbations (multipole, mutual), so that regular librations of the critical argument are not present. The resulting mass of Kleopatra, $m_1 = (1.49\pm0.16)\cdot10^{-12}\,M_\odot$ or $2.97\cdot10^{18}\,{\rm kg}$, is significantly lower than previously thought. An implication explained in the accompanying paper (Marchis et al.) is that (216) Kleopatra is a critically rotating body. △ Less
Submitted 19 May, 2021; originally announced May 2021.
-
Optically thin circumstellar medium in $β$ Lyr A system
Authors: M. Brož, D. Mourard, J. Budaj, P. Harmanec, H. Schmitt, I. Tallon-Bosc, D. Bonneau, H. Božić, D. Gies, M. Šlechta
Abstract: beta Lyr A is a complex binary system with an extensive observational dataset: light curves (from FUV to FIR), interferometric squared visibility, closure phase, triple product measurements, spectral-energy distribution (SED), high-resolution spectroscopy, differential visibility amplitude, and also differential phase. In particular, we use spectra from Ondrejov 2m telescope from 2013 to 2015, to… ▽ More beta Lyr A is a complex binary system with an extensive observational dataset: light curves (from FUV to FIR), interferometric squared visibility, closure phase, triple product measurements, spectral-energy distribution (SED), high-resolution spectroscopy, differential visibility amplitude, and also differential phase. In particular, we use spectra from Ondrejov 2m telescope from 2013 to 2015, to measure the emission in Halpha, HeI, SiII, NeI, or CII lines, and differential interferometry by CHARA/VEGA from the 2013 campaign to measure wavelength-dependent sizes across Halpha and HeI 6678. This allows us to constrain not only optically thick objects (primary, secondary, accretion disk), but also optically thin objects (disk atmosphere, jets, shell). We extended our modelling tool Pyshellspec (based on Shellspec; a 1D LTE radiative transfer code) to include all new observables, to compute differential visibilities/phases, to perform a Doppler tomography, and to determine a joint chi^2 metric. After an optimisation of 38 free parameters, we derive a robust model of the beta Lyr A system. According to the model, the emission is formed in an extended atmosphere of the disk, two perpendicular jets expanding at ~700 km s^-1, and a symmetric shell with the radius ~70 R_S. The spectroscopy indicates a low abundance of carbon, 10^-2 of the solar value. We also quantify systematic differences between datasets and discuss alternative models, with higher resolution, additional asymmetries, or He-rich abundance. △ Less
Submitted 12 October, 2020; originally announced October 2020.
-
A pair of Jovian Trojans at the L4 Lagrange point
Authors: Timothy R. Holt, David Vokrouhlický, David Nesvorný, Miroslav Brož, Jonathan Horner
Abstract: Asteroid pairs, two objects that are not gravitationally bound to one another, but share a common origin, have been discovered in the Main belt and Hungaria populations. Such pairs are of major interest, as the study of their evolution under a variety of dynamical influences can indicate the time since the pair was created. To date, no asteroid pairs have been found in the Jovian Trojans, despite… ▽ More Asteroid pairs, two objects that are not gravitationally bound to one another, but share a common origin, have been discovered in the Main belt and Hungaria populations. Such pairs are of major interest, as the study of their evolution under a variety of dynamical influences can indicate the time since the pair was created. To date, no asteroid pairs have been found in the Jovian Trojans, despite the presence of several binaries and collisional families in the population. The search for pairs in the Jovian Trojan population is of particular interest, given the importance of the Trojans as tracers of planetary migration during the Solar system's youth. Here we report a discovery of the first pair, (258656) 2002~ES$_{76}$ and 2013~CC$_{41}$, in the Jovian Trojans. The two objects are approximately the same size and are located very close to the L4 Lagrange point. Using numerical integrations, we find that the pair is at least $360$~Myr old, though its age could be as high as several Gyrs. The existence of the (258656) 2002~ES$_{76}$--2013~CC$_{41}$ pair implies there could be many such pairs scattered through the Trojan population. Our preferred formation mechanism for the newly discovered pair is through the dissociation of an ancient binary system, triggered by a sub-catastrophic impact, but we can not rule out rotation fission of a single object driven by YORP torques. A by-product of our work is an up-to-date catalog of Jovian Trojan proper elements, which we have made available for further studies. △ Less
Submitted 14 September, 2020; originally announced September 2020.
-
Physical and dynamical characterization of the Euphrosyne asteroid Family
Authors: B. Yang, J. Hanus, M. Broz, O. Chrenko, M. Willman, P. Sevecek, J. Masiero, H. Kaluna
Abstract: The Euphrosyne asteroid family occupies a unique zone in orbital element space around 3.15 au and may be an important source of the low-albedo near-Earth objects. The parent body of this family may have been one of the planetesimals that delivered water and organic materials onto the growing terrestrial planets. We aim to characterize the compositional properties as well as the dynamical propertie… ▽ More The Euphrosyne asteroid family occupies a unique zone in orbital element space around 3.15 au and may be an important source of the low-albedo near-Earth objects. The parent body of this family may have been one of the planetesimals that delivered water and organic materials onto the growing terrestrial planets. We aim to characterize the compositional properties as well as the dynamical properties of the family. We performed a systematic study to characterize the physical properties of the Euphrosyne family members via low-resolution spectroscopy using the IRTF telescope. In addition, we performed smoothed-particle hydrodynamics (SPH) simulations and N-body simulations to investigate the collisional origin, determine a realistic velocity field, study the orbital evolution, and constrain the age of the Euphrosyne family. Our spectroscopy survey shows that the family members exhibit a tight taxonomic distribution, suggesting a homogeneous composition of the parent body. Our SPH simulations are consistent with the Euphrosyne family having formed via a reaccumulation process instead of a cratering event. Finally, our N-body simulations indicate that the age of the family is 280 Myr +180/-80 Myr, which is younger than a previous estimate. △ Less
Submitted 9 September, 2020; originally announced September 2020.
-
Binary asteroid (31) Euphrosyne: Ice-rich and nearly spherical
Authors: B. Yang, J. Hanus, B. Carry, P. Vernazza, M. Broz, F. Vachier, N. Rambaux, M. Marsset, O. Chrenko, P. Sevecek, M. Viikinkoski, E. Jehin, M. Ferrais, E. Podlewska Gaca, A. Drouard, F. Marchis, M. Birlan, Z. Benkhaldoun, J. Berthier, P. Bartczak, C. Dumas, G. Dudzinski, J. Durech, J. Castillo-Rogez, F. Cipriani , et al. (16 additional authors not shown)
Abstract: Asteroid (31) Euphrosyne is one of the biggest objects in the asteroid main belt and the Euphrosyne family occupies a highly inclined region in the outer main belt and contains a remarkably large number of members, which is interpreted as an outcome of a disruptive cratering event. The goals of this adaptive-optics imaging study were threefold: to characterize the shape of Euphrosyne, to constrain… ▽ More Asteroid (31) Euphrosyne is one of the biggest objects in the asteroid main belt and the Euphrosyne family occupies a highly inclined region in the outer main belt and contains a remarkably large number of members, which is interpreted as an outcome of a disruptive cratering event. The goals of this adaptive-optics imaging study were threefold: to characterize the shape of Euphrosyne, to constrain its density, and to search for the large craters that may be associated with the family formation event. We obtained disk-resolved images of Euphrosyne using SPHERE/ZIMPOL at ESO's 8.2-m VLT as part of our large program (ID: 199.C-0074, PI: Vernazza). We reconstructed its 3D-shape using the adam shape modeling algorithm based on the SPHERE images and the available lightcurves of this asteroid. We analyzed the dynamics of the satellite with the genoid meta-heuristic algorithm. Finally, we studied the shape of Euphrosyne using hydrostatic equilibrium models. Our SPHERE observations show that Euphrosyne has a nearly spherical shape with the sphericity index of 0.9888 and its surface lacks large impact craters. Euphrosyne's diameter is 268+/-6 km, making it one of the top 10 largest main belt asteroids. We detected a satellite of Euphrosyne -- S/2019 (31) 1-- that is about 4 km across, on an circular orbit. The mass determined from the orbit of the satellite together with the volume computed from the shape model imply a density of 1665+/-242 kg/m^3, suggesting that Euphrosyne probably contain a large fraction of water ice in its interior. We find that the spherical shape of Euphrosyne is a result of the reaccumulation process following the impact, as in the case of (10) Hygiea. However, our shape analysis reveals that, contrary to Hygiea, the axis ratios of Euphrosyne significantly differ from the ones suggested by fluid hydrostatic equilibrium following reaccumulation. △ Less
Submitted 15 July, 2020; originally announced July 2020.
-
Impacts into rotating targets: angular momentum draining and efficient formation of synthetic families
Authors: P. Ševeček, M. Brož, M. Jutzi
Abstract: About 10% of the observed asteroids have rotational periods lower than P = 3 h and they seem to be relatively close to the spin barrier. Yet, the rotation has often been neglected in simulations of asteroid collisions. To determine the effect of rotation, we perform a large number of SPH/N-body impact simulations with rotating targets. We developed a new unified SPH/N-body code with self-gravity,… ▽ More About 10% of the observed asteroids have rotational periods lower than P = 3 h and they seem to be relatively close to the spin barrier. Yet, the rotation has often been neglected in simulations of asteroid collisions. To determine the effect of rotation, we perform a large number of SPH/N-body impact simulations with rotating targets. We developed a new unified SPH/N-body code with self-gravity, suitable for simulations of both fragmentation phase and gravitational reaccumulation. The code has been verified against previous ones (Benz and Asphaug 1994), but we also tested new features, e.g. rotational stability, tensile stability, etc. Using the new code, we ran simulations with D_pb = 10 km and 100 km monolithic targets and compared synthetic asteroid families created by these impacts with families corresponding to non-rotating targets. The rotation affects mostly cratering events at oblique impact angles. The total mass ejected by these collision can be up to five times larger for rotating targets. We further compute the transfer of the angular momentum and determine conditions under which impacts accelerate or decelerate the target. While individual cratering collisions can cause both acceleration and deceleration, the deceleration prevails on average, collisions thus cause a systematic spin-down of asteroid population. △ Less
Submitted 8 August, 2019; originally announced August 2019.
-
Properties and nature of Be stars 31. The binary nature, light variability, physical elements, and emission-line changes of HD~81357
Authors: P. Koubský, P. Harmanec, M. Brož, L. Kotková, S. Yang, H. Božić, D. Sudar, Y. Frémat, D. Korčákov'a, V. Votruba, P. Škoda, M. Šlechta, D. Ruždjak
Abstract: Reliable determination of the basic physical properties of hot emission-line binaries with Roche-lobe filling secondaries is important for developing the theory of mass exchange in binaries. It is not easy, however, due to the presence of circumstellar matter. Here, we report the first detailed investigation of a new representative of this class of binaries, HD~81357, based on the analysis of spec… ▽ More Reliable determination of the basic physical properties of hot emission-line binaries with Roche-lobe filling secondaries is important for developing the theory of mass exchange in binaries. It is not easy, however, due to the presence of circumstellar matter. Here, we report the first detailed investigation of a new representative of this class of binaries, HD~81357, based on the analysis of spectra and photometry from several observatories. HD~81357 was found to be a double-lined spectroscopic binary and an ellipsoidal variable seen under an intermediate orbital inclination of $\sim(63\pm5)^\circ$, having an orbital period of 33\fd77445(41) and a~circular orbit. From an automated comparison of the observed and synthetic spectra, we estimate the component's effective temperatures to be 12930(540)~K and 4260(24)~K. The combined light-curve and orbital solutions, also constrained by a very accurate Gaia Data Release 2 parallax, give the following values of the basic physical properties: masses $3.36\pm0.15$ and $0.34\pm0.04$~\Mnom, radii $3.9\pm0.2$ and 13.97\pm0.05$~\Rnom, and a~mass ratio $10.0\pm0.5$. Evolutionary modelling of the system including the phase of mass transfer between the components indicated that HD~81357 is a~system observed in the final slow phase of the mass exchange after the mass-ratio reversal. Contrary to what has been seen for similar binaries like AU~Mon, no cyclic light variations were found on a~time scale an~order of magnitude longer than the orbital period. 243,1 15% △ Less
Submitted 7 August, 2019; originally announced August 2019.
-
The shape of (7) Iris as evidence of an ancient large impact?
Authors: J. Hanuš, M. Marsset, P. Vernazza, M. Viikinkoski, A. Drouard, M. Brož, B. Carry, R. Fetick, F. Marchis, L. Jorda, T. Fusco, M. Birlan, T. Santana-Ros, E. Podlewska-Gaca, E. Jehin, M. Ferrais, J. Grice, P. Bartczak, J. Berthier, J. Castillo-Rogez, F. Cipriani, F. Colas, G. Dudzinski, C. Dumas, J. Ďurech , et al. (13 additional authors not shown)
Abstract: Asteroid (7) Iris is an ideal target for disk-resolved imaging owing to its brightness (V$\sim$7-8) and large angular size of 0.33 arcsec during its apparitions. Iris is believed to belong to the category of large unfragmented asteroids that avoided internal differentiation, implying that its current shape and topography may record the first few 100 Myr of the solar system's collisional evolution.… ▽ More Asteroid (7) Iris is an ideal target for disk-resolved imaging owing to its brightness (V$\sim$7-8) and large angular size of 0.33 arcsec during its apparitions. Iris is believed to belong to the category of large unfragmented asteroids that avoided internal differentiation, implying that its current shape and topography may record the first few 100 Myr of the solar system's collisional evolution. We recovered information about the shape and surface topography of Iris from disk-resolved VLT/SPHERE/ZIMPOL images acquired in the frame of our ESO large program. We used the All-Data Asteroid Modeling (ADAM) shape reconstruction algorithm to model the 3D shape of Iris, using optical disk-integrated data and disk-resolved images from SPHERE as inputs. We analyzed the SPHERE images to infer the asteroid's global shape and the morphology of its main craters. We present the volume-equivalent diameter D$_{\rm eq}$=214$\pm$5 km, and bulk density $ρ$=2.7$\pm$0.3 g cm$^{-3}$ of Iris. Its shape appears to be consistent with that of an oblate spheroid with a large equatorial excavation. We identified eight putative surface features 20--40 km in diameter detected at several epochs, which we interpret as impact craters. Craters on Iris have depth-to-diameter ratios that are similar to those of analogous 10 km craters on Vesta. The bulk density of Iris is consistent with that of its meteoritic analog, namely LL ordinary chondrites. Considering the absence of a collisional family related to Iris and the number of large craters on its surface, we suggest that its equatorial depression may be the remnant of an ancient (at least 3 Gyr) impact. Iris's shape further opens the possibility that large planetesimals formed as almost perfect oblate spheroids. Finally, we attribute the difference in crater morphology between Iris and Vesta to their different surface gravities. △ Less
Submitted 25 February, 2019; originally announced February 2019.
-
Binary planet formation by gas-assisted encounters of planetary embryos
Authors: Ondřej Chrenko, Miroslav Brož, David Nesvorný
Abstract: We present radiation hydrodynamic simulations in which binary planets form by close encounters in a system of several super-Earth embryos. The embryos are embedded in a protoplanetary disk consisting of gas and pebbles and evolve in a region where the disk structure supports convergent migration due to Type I torques. As the embryos accrete pebbles, they become heated and thus affected by the ther… ▽ More We present radiation hydrodynamic simulations in which binary planets form by close encounters in a system of several super-Earth embryos. The embryos are embedded in a protoplanetary disk consisting of gas and pebbles and evolve in a region where the disk structure supports convergent migration due to Type I torques. As the embryos accrete pebbles, they become heated and thus affected by the thermal torque (Benítez-Llambay et al. 2015) and the hot-trail effect (Chrenko et al. 2017) which excites orbital eccentricities. Motivated by findings of Eklund & Masset (2017), we assume the hot-trail effect operates also vertically and reduces the efficiency of inclination damping. Non-zero inclinations allow the embryos to become closely packed and also vertically stirred within the convergence zone. Subsequently, close encounters of two embryos assisted by the disk gravity can form transient binary planets which quickly dissolve. Binary planets with a longer lifetime $\sim$$10^{4}$ yr form in 3-body interactions of a transient pair with one of the remaining embryos. The separation of binary components generally decreases in subsequent encounters and due to pebble accretion until the binary merges, forming a giant planet core. We provide an order-of-magnitude estimate of the expected occurrence rate of binary planets, yielding one binary planet per $\simeq$$2$--$5\times10^{4}$ planetary systems. Therefore, although rare, the binary planets may exist in exoplanetary systems and they should be systematically searched for. △ Less
Submitted 30 October, 2018; originally announced October 2018.
-
A study of 3-dimensional shapes of asteroid families with an application to Eos
Authors: Miroslav Brož, Alessandro Morbidelli
Abstract: In order to fully understand the shapes of asteroids families in the 3-dimensional space of the proper elements $(a_{\rm p}, e_{\rm p}, \sin I_{\rm p})$ it is necessary to compare observed asteroids with N-body simulations. To this point, we describe a rigorous yet simple method which allows for a selection of the observed asteroids, assures the same size-frequency distribution of synthetic astero… ▽ More In order to fully understand the shapes of asteroids families in the 3-dimensional space of the proper elements $(a_{\rm p}, e_{\rm p}, \sin I_{\rm p})$ it is necessary to compare observed asteroids with N-body simulations. To this point, we describe a rigorous yet simple method which allows for a selection of the observed asteroids, assures the same size-frequency distribution of synthetic asteroids, accounts for a background population, and computes a $χ^2$ metric. We study the Eos family as an example, and we are able to fully explain its non-isotropic features, including the distribution of pole latitudes $β$. We confirm its age $t = (1.3\pm0.3)\,{\rm Gyr}$; while this value still scales with the bulk density, it is verified by a Monte-Carlo collisional model. The method can be applied to other populous families (Flora, Eunomia, Hygiea , Koronis, Themis, Vesta, etc.). △ Less
Submitted 9 October, 2018; originally announced October 2018.
-
Dynamics of multiple protoplanets embedded in gas/pebble disks and its dependence on $Σ$ and $ν$ parameters
Authors: Miroslav Brož, Ondřej Chrenko, David Nesvorný, Michiel Lambrechts
Abstract: Protoplanets of Super-Earth sizes may get trapped in convergence zones for planetary migration and form gas giants there. These growing planets undergo accretion heating, which triggers a hot-trail effect that can reverse migration directions, increase eccentricities and prevent resonant captures (Chrenko et al. 2017). We study populations of embryos accreting pebbles using Fargo-Thorin 2D hydroco… ▽ More Protoplanets of Super-Earth sizes may get trapped in convergence zones for planetary migration and form gas giants there. These growing planets undergo accretion heating, which triggers a hot-trail effect that can reverse migration directions, increase eccentricities and prevent resonant captures (Chrenko et al. 2017). We study populations of embryos accreting pebbles using Fargo-Thorin 2D hydrocode. We find that embryos in a disk with high surface density ($Σ_0 = 990\,{\rm g}\,{\rm cm}^{-2}$) undergo `unsuccessful' two-body encounters which do not lead to a merger. Only when a 3rd protoplanet arrives to the convergence zone, three-body encounters lead to mergers. For a low-viscosity disk ($ν= 5\times10^{13}\,{\rm cm}^2\,{\rm s}^{-1}$) a massive coorbital is a possible outcome, for which a pebble isolation develops and the coorbital is stabilised. For more massive protoplanets ($5\,M_\oplus$), the convergence radius is located further out, in the ice-giant zone. After a series of encounters, there is an evolution driven by a dynamical torque of a tadpole region, which is systematically repeated several times, until the coorbital configuration is disrupted and planets merge. This may be a pathway how to solve the problem that coorbitals often form in simulations but they are not observed in nature. In contrast, the joint evolution of 120 low-mass protoplanets ($0.1\,M_\oplus$) reveals completely different dynamics. The evolution is no longer smooth, but rather a random walk. This is because the spiral arms, developed in the gas disk due to Lindblad resonances, overlap with each other and affect not only a single protoplanet but several in the surroundings. Our hydrodynamical simulations may have important implications for N-body simulations of planetary migration that use simplified torque prescriptions and are thus unable to capture protoplanet dynamics in its full glory. △ Less
Submitted 5 December, 2018; v1 submitted 8 October, 2018; originally announced October 2018.
-
Physical properties of beta Lyr A and its opaque accretion disk
Authors: D. Mourard, M. Broz, J. Nemravova, P. Harmanec, J. Budaj, F. Baron, J. Monnier, G. Schaefer, H. Schmitt, I. Tallon-Bosc, J. Armstrong, E. Baines, D. Bonneau, H. Bozic, J. M. Clausse, C. Farrington, D. Gies, J. Jurysek, D. Korcakova, H. McAlister, A. Meilland, N. Nardetto, P. Svoboad, M. Slechta, M. Wolf , et al. (1 additional authors not shown)
Abstract: Mass exchange and mass loss in close binaries can significantly affect their evolution, but a complete self-consistent theory of these processes is still to be developed. Processes such as radiative shielding due to a~hot-spot region, or a~hydrodynamical interaction of different parts of the gas stream have been studied previously. In order to test the respective predictions, it is necessary to ca… ▽ More Mass exchange and mass loss in close binaries can significantly affect their evolution, but a complete self-consistent theory of these processes is still to be developed. Processes such as radiative shielding due to a~hot-spot region, or a~hydrodynamical interaction of different parts of the gas stream have been studied previously. In order to test the respective predictions, it is necessary to carry out detailed observations of binaries undergoing the largescale mass exchange, especially for those that are in the rapid transfer phase. \bla is an archetype of such a system, having a long and rich observational history. Our goal for this first study is to quantitatively estimate the geometry and physical properties of the optically thick components, namely the Roche-lobe filling mass-losing star, and the accretion disk surrounding the mass-gaining star of \blae. A series of continuum visible and NIR spectro-interferometric observations by the NPOI, CHARA/MIRC and VEGA instruments covering the whole orbit of \bla acquired during a~two-week campaign in 2013 were complemented with \ubvr\ photometric observations acquired during a three-year monitoring of the system. We included NUV and FUV observations from OAO~A-2, IUE, and Voyager satellites. △ Less
Submitted 12 July, 2018; originally announced July 2018.
-
SPH/N-body simulations of small (D = 10 km) asteroidal breakups and improved parametric relations for Monte-Carlo collisional models
Authors: P. Ševeček, M. Brož, D. Nesvorný, B. Enke, D. Durda, K. Walsh, D. C. Richardson
Abstract: We report on our study of asteroidal breakups, i.e. fragmentations of targets, subsequent gravitational reaccumulation and formation of small asteroid families. We focused on parent bodies with diameters $D_{\rm pb} = 10$ km. Simulations were performed with a smoothed-particle hydrodynamics (SPH) code combined with an efficient N-body integrator. We assumed various projectile sizes, impact velocit… ▽ More We report on our study of asteroidal breakups, i.e. fragmentations of targets, subsequent gravitational reaccumulation and formation of small asteroid families. We focused on parent bodies with diameters $D_{\rm pb} = 10$ km. Simulations were performed with a smoothed-particle hydrodynamics (SPH) code combined with an efficient N-body integrator. We assumed various projectile sizes, impact velocities and impact angles (125 runs in total). Resulting size-frequency distributions are significantly different from scaled-down simulations with $D_{\rm pb} = 100$ km targets (Durda et al., 2007). We derive new parametric relations describing fragment distributions, suitable for Monte-Carlo collisional models. We also characterize velocity fields and angular distributions of fragments, which can be used as initial conditions for N-body simulations of small asteroid families. Finally, we discuss a number of uncertainties related to SPH simulations. △ Less
Submitted 28 March, 2018; originally announced March 2018.
-
CzeV - The Czech Variable Star Catalogue
Authors: Marek Skarka, Martin Mašek, Luboš Brát, Pavel Cagaš, Jakub Juryšek, Kateřina Hoňková, Miloslav Zejda, Ladislav Šmelcer, Martin Jelínek, František Lomoz, Martin Tylšar, Jaroslav Trnka, Ondřej Pejcha, Pavel Pintr, Martin Lehký, Jan Janík, Ladislav Červinka, Václav Přibík, David Motl, Filip Walter, Petr Zasche, Karel Koss, Petr Hájek, František Bílek, Jiří Liška , et al. (40 additional authors not shown)
Abstract: We present the first release of the Czech Variable star catalogue that currently contains 1228 stars whose variability was discovered by 60 Czech observers. The catalogue contains confirmed variable stars of various types, but also candidates. We give precise coordinates, cross identification with other catalogues, information about constellation, variability type, brightness, light elements, name… ▽ More We present the first release of the Czech Variable star catalogue that currently contains 1228 stars whose variability was discovered by 60 Czech observers. The catalogue contains confirmed variable stars of various types, but also candidates. We give precise coordinates, cross identification with other catalogues, information about constellation, variability type, brightness, light elements, name of the discoverer and year of discovery. In eighty-eight percent of stars the variability type is estimated, for more than 60 % of the stars the light ephemerides are given. △ Less
Submitted 26 September, 2017; originally announced September 2017.
-
Improved model of the triple system V746 Cas that has a bipolar magnetic field associated with the tertiary
Authors: P. Harmanec, M. Brož, P. Mayer, P. Zasche, L. Kotková, J. A. Nemravová, R. J. Dukes, D. Korčáková, M. Šlechta, E. Kıran, R. Kříček, J. Juryšek
Abstract: V746 Cas is known to be a triple system composed of a close binary with an alternatively reported period of either 25.4d or 27.8d and a third component in a 62000d orbit. The object was also reported to exhibit multiperiodic light variations with periods from 0.83d to 2.50d, on the basis of which it was classified as a slowly pulsating B star. Interest in further investigation of this system was r… ▽ More V746 Cas is known to be a triple system composed of a close binary with an alternatively reported period of either 25.4d or 27.8d and a third component in a 62000d orbit. The object was also reported to exhibit multiperiodic light variations with periods from 0.83d to 2.50d, on the basis of which it was classified as a slowly pulsating B star. Interest in further investigation of this system was raised by the detection of a variable magnetic field. Analysing spectra from four instruments, earlier published radial velocities, and several sets of photometric observations, we arrived at the following conclusions: (1) The optical spectrum is dominated by the lines of the B-type primary (Teff1~16500(100) K), contributing 70% of the light in the optical region, and a slightly cooler B tertiary (Teff3~13620(150) K). The lines of the low-mass secondary are below our detection threshold; we estimate that it could be a normal A or F star. (2) We resolved the ambiguity in the value of the inner binary period and arrived at a linear ephemeris of T_super.conj.=HJD 2443838.78(81)+25.41569(42)xE. (3) The intensity of the magnetic field undergoes a~sinusoidal variation in phase with one of the known photometric periods, namely 2.503867(19)d, which we identify with the rotational period of the tertiary. (4) The second photometric 1.0649524(40)d period is identified with the rotational period of the B-type primary, but this interpretation is much less certain and needs further verification. (5) If our interpretation of photometric periods is confirmed, the classification of the object as a slowly pulsating B star should be revised. (6) Applying an N-body model to different types of available observational data, we constrain the orbital inclination of the inner orbit to ~60 deg to 85 deg even in the absence of eclipses, and estimate the probable properties of the triple system and its components. △ Less
Submitted 19 September, 2017; originally announced September 2017.
-
Shape models of asteroids based on lightcurve observations with BlueEye600 robotic observatory
Authors: Josef Durech, Josef Hanus, Miroslav Broz, Martin Lehky, Raoul Behrend, Pierre Antonini, Stephane Charbonnel, Roberto Crippa, Pierre Dubreuil, Gino Farroni, Gilles Kober, Alain Lopez, Federico Manzini, Julian Oey, Raymond Poncy, Claudine Rinner, Rene Roy
Abstract: We present physical models, i.e. convex shapes, directions of the rotation axis, and sidereal rotation periods, of 18 asteroids out of which 10 are new models and 8 are refined models based on much larger data sets than in previous work. The models were reconstructed by the lightcurve inversion method from archived publicly available lightcurves and our new observations with BlueEye600 robotic obs… ▽ More We present physical models, i.e. convex shapes, directions of the rotation axis, and sidereal rotation periods, of 18 asteroids out of which 10 are new models and 8 are refined models based on much larger data sets than in previous work. The models were reconstructed by the lightcurve inversion method from archived publicly available lightcurves and our new observations with BlueEye600 robotic observatory. One of the new results is the shape model of asteroid (1663)~van~den~Bos with the rotation period of 749\,hr, which makes it the slowest rotator with known shape. We describe our strategy for target selection that aims at fast production of new models using the enormous potential of already available photometry stored in public databases. We also briefly describe the control software and scheduler of the robotic observatory and we discuss the importance of building a database of asteroid models for studying asteroid physical properties in collisional families. △ Less
Submitted 12 July, 2017; originally announced July 2017.
-
Eccentricity excitation and merging of planetary embryos heated by pebble accretion
Authors: Ondřej Chrenko, Miroslav Brož, Michiel Lambrechts
Abstract: Context: Planetary embryos can continue to grow by pebble accretion until they become giant planet cores. Simultaneously, these embryos mutually interact and also migrate due to torques arising from the protoplanetary disk. Aims: Our aim is to investigate how pebble accretion alters the orbital evolution of embryos undergoing the Type-I migration. In particular, we study whether they establish r… ▽ More Context: Planetary embryos can continue to grow by pebble accretion until they become giant planet cores. Simultaneously, these embryos mutually interact and also migrate due to torques arising from the protoplanetary disk. Aims: Our aim is to investigate how pebble accretion alters the orbital evolution of embryos undergoing the Type-I migration. In particular, we study whether they establish resonant chains, whether these chains are prone to instabilities and if giant planet cores form through embryo merging, thus occurring more rapidly than by pebble accretion alone. Methods: For the first time, we perform self-consistent global-scale radiative hydrodynamic simulations of a two-fluid protoplanetary disk consisting of gas and pebbles, the latter being accreted by embedded embryos. Accretion heating, along with other radiative processes, is accounted for to correctly model the Type-I migration. Results: We track the evolution of four super-Earth-like embryos, initially located in a region where the disk structure allows for a convergent migration. Generally, embryo merging is facilitated by rapidly increasing embryo masses and breaks the otherwise oligarchic growth. Moreover, we find that the orbital eccentricity of each embryo is considerably excited ($\simeq0.03$) due to the presence of an asymmetric underdense lobe of gas, a so-called `hot trail', produced by accretion heating of the embryo's vicinity. Eccentric orbits lead the embryos to frequent close encounters and make resonant locking more difficult. Conclusions: Embryo merging typically produces one massive core ($\gtrsim 10\,\mathrm{M_{E}}$) in our simulations, orbiting near $10\,\mathrm{AU}$. Pebble accretion is naturally accompanied by occurrence of eccentric orbits which should be considered in future efforts to explain the structure of exoplanetary systems. △ Less
Submitted 20 June, 2017; originally announced June 2017.
-
An advanced N-body model for interacting multiple stellar systems
Authors: M. Brož
Abstract: We construct an advanced model for interacting multiple stellar systems in which we compute all trajectories with a numerical N-body integrator, namely the Bulirsch--Stoer from the SWIFT package. We can then derive various observables: astrometric positions, radial velocities, minima timings (TTVs), eclipse durations, interferometric visibilities, closure phases, synthetic spectra, spectral-energy… ▽ More We construct an advanced model for interacting multiple stellar systems in which we compute all trajectories with a numerical N-body integrator, namely the Bulirsch--Stoer from the SWIFT package. We can then derive various observables: astrometric positions, radial velocities, minima timings (TTVs), eclipse durations, interferometric visibilities, closure phases, synthetic spectra, spectral-energy distribution, and even complete light curves. We use a modified version of the Wilson--Devinney code for the latter, in which the instantaneous true phase and inclination of the eclipsing binary are governed by the N-body integration. If one has all kinds of observations at disposal, a joint $χ^2$ metric and an optimisation algorithm (a~simplex or simulated annealing) allows to search for a global minimum and construct very robust models of stellar systems. At the same time, our N-body model is free from artefacts which may arise if mutual gravitational interactions among all components are not self-consistently accounted for. Finally, we present a number of examples showing dynamical effects that can be studied with our code and we discuss how systematic errors may affect the results (and how to prevent this from happening). △ Less
Submitted 14 June, 2017; originally announced June 2017.
-
Hektor - an exceptional D-type family among Jovian Trojans
Authors: Jakub Rozehnal, Miroslav Brož, David Nesvorný, Daniel D. Durda, Kevin Walsh, Derek C. Richardson, Erik Asphaug
Abstract: In this work, we analyze Jovian Trojans in the space of suitable resonant elements and we identify clusters of possible collisional origin by two independent methods: the hierarchical clustering and a so-called "randombox". Compared to our previous work (Brož and Rozehnal 2011), we study a twice larger sample. Apart from Eurybates, Ennomos and $1996\,\rm RJ$ families, we have found three more clus… ▽ More In this work, we analyze Jovian Trojans in the space of suitable resonant elements and we identify clusters of possible collisional origin by two independent methods: the hierarchical clustering and a so-called "randombox". Compared to our previous work (Brož and Rozehnal 2011), we study a twice larger sample. Apart from Eurybates, Ennomos and $1996\,\rm RJ$ families, we have found three more clusters --- namely families around asteroids (20961)~Arkesilaos, (624)~Hektor in the $L_4$ libration zone and (247341)~$2001\,\rm UV_{209}$ in $L_5$. The families fulfill our stringent criteria, i.e. a high statistical significance, an albedo homogeneity and a steeper size-frequency distribution than that of background. In order to understand their nature, we simulate their long term collisional evolution with the Boulder code (Morbidelli et al. 2009) and dynamical evolution using a modified SWIFT integrator (Levison and Duncan, 1994). Within the framework of our evolutionary model, we were able to constrain the the age of the Hektor family to be either 1 to 4 Gyr or, less likely, 0.1 to 2.5 Gyr, depending on initial impact geometry. Since (624) Hektor itself seems to be a bilobed--shape body with a satellite (Marchis et al. 2014), i.e. an exceptional object, we address its association with the D--type family and we demonstrate that the moon and family could be created during a single impact event. We simulated the cratering event using a Smoothed Particle Hydrodynamics (SPH, Benz and Asphaug, 1994). This is also the first case of a family associated with a D--type parent body. △ Less
Submitted 15 July, 2016; originally announced July 2016.
-
The origin of long-lived asteroids in the 2:1 mean-motion resonance with Jupiter
Authors: O. Chrenko, M. Brož, D. Nesvorný, K. Tsiganis, D. K. Skoulidou
Abstract: The 2:1 mean-motion resonance with Jupiter harbours two distinct groups of asteroids. The short-lived population is known to be a transient group sustained in steady state by the Yarkovsky semimajor axis drift. The long-lived asteroids, however, can exhibit dynamical lifetimes comparable to $4\,\mathrm{Gyr}$. They reside near two isolated islands of the phase space denoted $\mathrm{A}$ and… ▽ More The 2:1 mean-motion resonance with Jupiter harbours two distinct groups of asteroids. The short-lived population is known to be a transient group sustained in steady state by the Yarkovsky semimajor axis drift. The long-lived asteroids, however, can exhibit dynamical lifetimes comparable to $4\,\mathrm{Gyr}$. They reside near two isolated islands of the phase space denoted $\mathrm{A}$ and $\mathrm{B}$, with an uneven population ratio $\mathrm{B}/\mathrm{A} \simeq 10$. The orbits of $\mathrm{A}$-island asteroids are predominantly highly inclined, compared to island $\mathrm{B}$. The size-frequency distribution is steep but the orbital distribution lacks any evidence of a collisional cluster. These observational constraints are somewhat puzzling and therefore the origin of the long-lived asteroids has not been explained so far. With the aim to provide a viable explanation, we first update the resonant population and revisit its physical properties. Using an $N$-body model with seven planets and the Yarkovsky effect included, we demonstrate that the dynamical depletion of island $\mathrm{A}$ is faster, in comparison with island $\mathrm{B}$. Then we investigate (i) the survivability of primordial resonant asteroids and (ii) capture of the population during planetary migration, following a recently described scenario with an escaping fifth giant planet and a jumping-Jupiter instability. We also model the collisional evolution of the resonant population over past $4\,\mathrm{Gyr}$. Our conclusion is that the long-lived group was created by resonant capture from a narrow part of hypothetical outer main-belt family during planetary migration. Primordial asteroids surviving the migration were probably not numerous enough to substantially contribute to the observed population. △ Less
Submitted 16 May, 2015; originally announced May 2015.
-
The thermal emission from boulders on (25143) Itokawa and general implications for the YORP effect
Authors: P. Ševeček, M. Brož, D. Čapek, J. Ďurech
Abstract: Infrared radiation emitted from an asteroid surface causes a torque that can significantly affect rotational state of the asteroid. The influence of small topographic features on this phenomenon, called the YORP effect, seems to be of utmost importance. In this work, we show that a lateral heat diffusion in boulders of suitable sizes leads to an emergence of a local YORP effect which magnitude is… ▽ More Infrared radiation emitted from an asteroid surface causes a torque that can significantly affect rotational state of the asteroid. The influence of small topographic features on this phenomenon, called the YORP effect, seems to be of utmost importance. In this work, we show that a lateral heat diffusion in boulders of suitable sizes leads to an emergence of a local YORP effect which magnitude is comparable to the YORP effect due to the global shape. We solve a three-dimensional heat diffusion equation in a boulder and its surroundings by the finite element method, using the FreeFem++ code. The contribution to the total torque is inferred from the computed temperature distribution. Our general approach allows us to compute the torque induced by a realistic irregular boulder. For an idealized boulder, our result is consistent with an existing one-dimensional model. We also estimated (and extrapolated) a size distribution of boulders on (25143) Itokawa from close-up images of its surface. We realized that topographic features on Itokawa can potentially induce a torque corresponding to a rotational acceleration of the order 10^-7 rad day^-2 and can therefore explain the observed phase shift in light curves. △ Less
Submitted 2 April, 2015; originally announced April 2015.
-
Identification and Dynamical Properties of Asteroid Families
Authors: D. Nesvorny, M. Broz, V. Carruba
Abstract: Asteroids formed in a dynamically quiescent disk but their orbits became gravitationally stirred enough by Jupiter to lead to high-speed collisions. As a result, many dozen large asteroids have been disrupted by impacts over the age of the Solar System, producing groups of fragments known as asteroid families. Here we explain how the asteroid families are identified, review their current inventory… ▽ More Asteroids formed in a dynamically quiescent disk but their orbits became gravitationally stirred enough by Jupiter to lead to high-speed collisions. As a result, many dozen large asteroids have been disrupted by impacts over the age of the Solar System, producing groups of fragments known as asteroid families. Here we explain how the asteroid families are identified, review their current inventory, and discuss how they can be used to get insights into long-term dynamics of main belt asteroids. Electronic tables of the membership for 122 notable families are reported on the Planetary Data System node. △ Less
Submitted 5 February, 2015; originally announced February 2015.
-
A six-part collisional model of the main asteroid belt
Authors: Helena Cibulková, Miroslav Brož, Paula G. Benavidez
Abstract: In this work, we construct a new model for the collisional evolution of the main asteroid belt. Our goals are to test the scaling law of Benz and Asphaug (1999) and ascertain if it can be used for the whole belt. We want to find initial size-frequency distributions (SFDs) for the considered six parts of the belt (inner, middle, 'pristine', outer, Cybele zone, high-inclination region) and to verify… ▽ More In this work, we construct a new model for the collisional evolution of the main asteroid belt. Our goals are to test the scaling law of Benz and Asphaug (1999) and ascertain if it can be used for the whole belt. We want to find initial size-frequency distributions (SFDs) for the considered six parts of the belt (inner, middle, 'pristine', outer, Cybele zone, high-inclination region) and to verify if the number of synthetic asteroid families created during the simulation matches the number of observed families as well. We used new observational data from the WISE satellite (Masiero et al., 2011) to construct the observed SFDs. We simulate mutual collisions of asteroids with a modified version of the Boulder code (Morbidelli et al., 2009), where the results of hydrodynamic (SPH) simulations of Durda et al. (2007) and Benavidez et al. (2012) are included. Because material characteristics can significantly affect breakups, we created two models - for monolithic asteroids and for rubble-piles. To explain the observed SFDs in the size range D = 1 to 10 km we have to also account for dynamical depletion due to the Yarkovsky effect. The assumption of (purely) rubble-pile asteroids leads to a significantly worse fit to the observed data, so that we can conclude that majority of main-belt asteroids are rather monolithic. Our work may also serve as a motivation for further SPH simulations of disruptions of smaller targets (with a parent body size of the order of 1 km). △ Less
Submitted 23 July, 2014; originally announced July 2014.
-
An anisotropic distribution of spin vectors in asteroid families
Authors: J. Hanuš, M. Brož, J. Ďurech, B. D. Warner, J. Brinsfield, R. Durkee, D. Higgins, R. A. Koff, J. Oey, F. Pilcher, R. Stephens, L. P. Strabla, Q. Ulisse, R. Girelli
Abstract: Current amount of ~500 asteroid models derived from the disk-integrated photometry by the lightcurve inversion method allows us to study not only the spin-vector properties of the whole population of MBAs, but also of several individual collisional families. We create a data set of 152 asteroids that were identified by the HCM method as members of ten collisional families, among them are 31 newly… ▽ More Current amount of ~500 asteroid models derived from the disk-integrated photometry by the lightcurve inversion method allows us to study not only the spin-vector properties of the whole population of MBAs, but also of several individual collisional families. We create a data set of 152 asteroids that were identified by the HCM method as members of ten collisional families, among them are 31 newly derived unique models and 24 new models with well-constrained pole-ecliptic latitudes of the spin axes. The remaining models are adopted from the DAMIT database or the literature. We revise the preliminary family membership identification by the HCM method according to several additional criteria - taxonomic type, color, albedo, maximum Yarkovsky semi-major axis drift and the consistency with the size-frequency distribution of each family, and consequently we remove interlopers. We then present the spin-vector distributions for eight asteroidal families. We use a combined orbital- and spin-evolution model to explain the observed spin-vector properties of objects among collisional families. In general, we observe for studied families similar trends in the (a_p, β) space: (i) larger asteroids are situated in the proximity of the center of the family; (ii) asteroids with β>0° are usually found to the right from the family center; (iii) on the other hand, asteroids with β<0° to the left from the center; (iv) majority of asteroids have large pole-ecliptic latitudes (|β|\gtrsim 30°); and finally (v) some families have a statistically significant excess of asteroids with β>0° or β<0°. Our numerical simulation of the long-term evolution of a collisional family is capable of reproducing well the observed spin-vector properties. Using this simulation, we also independently constrain the age of families Flora (1.0\pm0.5 Gyr) and Koronis (2.5-4 Gyr). △ Less
Submitted 17 September, 2013; originally announced September 2013.
-
The Eos family halo
Authors: Miroslav Brož, Alessandro Morbidelli
Abstract: We study K-type asteroids in the broad surroundings of the Eos family because they seem to be intimately related, according to their colours measured by the Sloan Digital Sky Survey. Such `halos' of asteroid families have been rarely used as constraints for dynamical studies to date. We explain its origin as bodies escaping from the family `core' due to the Yarkovsky semimajor-axis drift and inter… ▽ More We study K-type asteroids in the broad surroundings of the Eos family because they seem to be intimately related, according to their colours measured by the Sloan Digital Sky Survey. Such `halos' of asteroid families have been rarely used as constraints for dynamical studies to date. We explain its origin as bodies escaping from the family `core' due to the Yarkovsky semimajor-axis drift and interactions with gravitational resonances, mostly with the 9/4 mean-motion resonance with Jupiter at 3.03 AU. Our N-body dynamical model allows us to independently estimate the age of the family 1.5 to 1.9 Gyr. This is approximately in agreement with the previous age estimate by Vokrouhlický et al. (2006) based on a simplified model (which accounts only for changes of semimajor axis). We can also constrain the geometry of the disruption event which had to occur at the true anomaly f =~ 120 to 180 deg. △ Less
Submitted 6 February, 2013; originally announced February 2013.
-
Asteroids' physical models from combined dense and sparse photometry and scaling of the YORP effect by the observed obliquity distribution
Authors: J. Hanuš, J. Ďurech, M. Brož, A. Marciniak, B. D. Warner, F. Pilcher, R. Stephens, R. Behrend, B. Carry, D. Čapek, P. Antonini, M. Audejean, K. Augustesen, E. Barbotin, P. Baudouin, A. Bayol, L. Bernasconi, W. Borczyk, J. -G. Bosch, E. Brochard, L. Brunetto, S. Casulli, A. Cazenave, S. Charbonnel, B. Christophe , et al. (95 additional authors not shown)
Abstract: The larger number of models of asteroid shapes and their rotational states derived by the lightcurve inversion give us better insight into both the nature of individual objects and the whole asteroid population. With a larger statistical sample we can study the physical properties of asteroid populations, such as main-belt asteroids or individual asteroid families, in more detail. Shape models can… ▽ More The larger number of models of asteroid shapes and their rotational states derived by the lightcurve inversion give us better insight into both the nature of individual objects and the whole asteroid population. With a larger statistical sample we can study the physical properties of asteroid populations, such as main-belt asteroids or individual asteroid families, in more detail. Shape models can also be used in combination with other types of observational data (IR, adaptive optics images, stellar occultations), e.g., to determine sizes and thermal properties. We use all available photometric data of asteroids to derive their physical models by the lightcurve inversion method and compare the observed pole latitude distributions of all asteroids with known convex shape models with the simulated pole latitude distributions. We used classical dense photometric lightcurves from several sources and sparse-in-time photometry from the U.S. Naval Observatory in Flagstaff, Catalina Sky Survey, and La Palma surveys (IAU codes 689, 703, 950) in the lightcurve inversion method to determine asteroid convex models and their rotational states. We also extended a simple dynamical model for the spin evolution of asteroids used in our previous paper. We present 119 new asteroid models derived from combined dense and sparse-in-time photometry. We discuss the reliability of asteroid shape models derived only from Catalina Sky Survey data (IAU code 703) and present 20 such models. By using different values for a scaling parameter cYORP (corresponds to the magnitude of the YORP momentum) in the dynamical model for the spin evolution and by comparing synthetics and observed pole-latitude distributions, we were able to constrain the typical values of the cYORP parameter as between 0.05 and 0.6. △ Less
Submitted 29 January, 2013; originally announced January 2013.
-
Constraining the cometary flux through the asteroid belt during the late heavy bombardment
Authors: M. Brož, A. Morbidelli, W. F. Bottke, J. Rozehnal, D. Vokrouhlický, D. Nesvorný
Abstract: In the Nice model, the late heavy bombardment (LHB) is related to an orbital instability of giant planets which causes a fast dynamical dispersion of a transneptunian cometary disk. We study effects produced by these hypothetical cometary projectiles on main-belt asteroids. In particular, we want to check whether the observed collisional families provide a lower or an upper limit for the cometary… ▽ More In the Nice model, the late heavy bombardment (LHB) is related to an orbital instability of giant planets which causes a fast dynamical dispersion of a transneptunian cometary disk. We study effects produced by these hypothetical cometary projectiles on main-belt asteroids. In particular, we want to check whether the observed collisional families provide a lower or an upper limit for the cometary flux during the LHB. We present an updated list of observed asteroid families as identified in the space of synthetic proper elements by the hierarchical clustering method, colour data, albedo data and dynamical considerations and we estimate their physical parameters. We selected 12 families which may be related to the LHB according to their dynamical ages. We then used collisional models and N-body orbital simulations to gain insight into the long-term dynamical evolution of synthetic LHB families over 4 Gyr. We account for the mutual collisions, the physical disruptions of comets, the Yarkovsky/YORP drift, chaotic diffusion, or possible perturbations by the giant-planet migration. Assuming a "standard" size-frequency distribution of primordial comets, we predict the number of families with parent-body sizes D_PB >= 200 km which seems consistent with observations. However, more than 100 asteroid families with D_PB >= 100 km should be created at the same time which are not observed. This discrepancy can be nevertheless explained by the following processes: i) asteroid families are efficiently destroyed by comminution (via collisional cascade), ii) disruptions of comets below some critical perihelion distance (q <~ 1.5 AU) are common. Given the freedom in the cometary-disruption law, we cannot provide stringent limits on the cometary flux, but we can conclude that the observed distribution of asteroid families does not contradict with a cometary LHB. △ Less
Submitted 26 January, 2013; originally announced January 2013.
-
Did the Hilda collisional family form during the late heavy bombardment?
Authors: M. Brož, D. Vokrouhlický, A. Morbidelli, D. Nesvorný, W. F. Bottke
Abstract: We model the long-term evolution of the Hilda collisional family located in the 3/2 mean-motion resonance with Jupiter. Its eccentricity distribution evolves mostly due to the Yarkovsky/YORP effect and assuming that: (i) impact disruption was isotropic, and (ii) albedo distribution of small asteroids is the same as for large ones, we can estimate the age of the Hilda family to be… ▽ More We model the long-term evolution of the Hilda collisional family located in the 3/2 mean-motion resonance with Jupiter. Its eccentricity distribution evolves mostly due to the Yarkovsky/YORP effect and assuming that: (i) impact disruption was isotropic, and (ii) albedo distribution of small asteroids is the same as for large ones, we can estimate the age of the Hilda family to be $4_{-1}^{+0}\,{\rm Gyr}$. We also calculate collisional activity in the J3/2 region. Our results indicate that current collisional rates are very low for a 200\,km parent body such that the number of expected events over Gyrs is much smaller than one. The large age and the low probability of the collisional disruption lead us to the conclusion that the Hilda family might have been created during the Late Heavy Bombardment when the collisions were much more frequent. The Hilda family may thus serve as a test of orbital behavior of planets during the LHB. We tested the influence of the giant-planet migration on the distribution of the family members. The scenarios that are consistent with the observed Hilda family are those with fast migration time scales $\simeq 0.3\,{\rm Myr}$ to $3\,{\rm Myr}$, because longer time scales produce a family that is depleted and too much spread in eccentricity. Moreover, there is an indication that Jupiter and Saturn were no longer in a compact configuration (with period ratio $P_{\rm S}/P_{\rm J} > 2.09$) at the time when the Hilda family was created. △ Less
Submitted 6 September, 2011; originally announced September 2011.
-
Eurybates - the only asteroid family among Trojans?
Authors: Miroslav Brož, Jakub Rozehnal
Abstract: We study orbital and physical properties of Trojan asteroids of Jupiter. We try to discern all families previously discussed in literature, but we conclude there is only one significant family among Trojans, namely the cluster around asteroid (3548) Eurybates. It is the only cluster, which has all of the following characteristics: (i) it is clearly concentrated in the proper-element space; (ii) si… ▽ More We study orbital and physical properties of Trojan asteroids of Jupiter. We try to discern all families previously discussed in literature, but we conclude there is only one significant family among Trojans, namely the cluster around asteroid (3548) Eurybates. It is the only cluster, which has all of the following characteristics: (i) it is clearly concentrated in the proper-element space; (ii) size-frequency distribution is different from background asteroids; (iii) we have a reasonable collisional/dynamical model of the family. Henceforth, we can consider it as a real collisional family. We also report a discovery of a possible family around the asteroid (4709) Ennomos, composed mostly of small asteroids. The asteroid (4709) Ennomos is known to have a very high albedo $p_V \simeq 0.15$, which may be related to a hypothetical cratering event which exposed ice (Fernández et al. 2003). The relation between the collisional family and the exposed surface of the parent body is a unique opportunity to study the physics of cratering events. However, more data are needed to confirm the existence of this family and its relationship with Ennomos. △ Less
Submitted 6 September, 2011; originally announced September 2011.
-
V2368 Oph: An eclipsing and double-lined spectroscopic binary used as a photometric comparison star for U Oph
Authors: P. Harmanec, H. Božić, P. Mayer, P. Eenens, M. Brož, M. Wolf, S. Yang, M. Šlechta, D. Ruždjak, D. Sudar, H. Ak
Abstract: The A-type star HR 6412 = V2368 Oph was used by several investigators as a photometric comparison star for the known eclipsing binary U Oph but was found to be variable by three independent groups, including us. By analysing series of new spectral and photometric observations and a critical compilation of available radial velocities, we were able to find the correct period of light and radial-velo… ▽ More The A-type star HR 6412 = V2368 Oph was used by several investigators as a photometric comparison star for the known eclipsing binary U Oph but was found to be variable by three independent groups, including us. By analysing series of new spectral and photometric observations and a critical compilation of available radial velocities, we were able to find the correct period of light and radial-velocity variations and demonstrate that the object is an eclipsing and double-lined spectroscopic binary moving in a highly eccentric orbit. We derived a linear ephemeris T min.I = HJD (2454294.67 +/- 0.01) + (38.32712 +/- 0.00004)d x E and estimated preliminary basic physical properties of the binary. The dereddened UBV magnitudes and effective temperatures of the primary and secondary, based on our light- and velocity-curve solutions, led to distance estimates that agree with the Hipparcos distance within the errors. We find that the mass ratio must be close to one, but the limited number and wavelength range of our current spectra does not allow a truly precise determination of the binary masses. Nevertheless, our results show convincingly that both binary components are evolved away from the main sequence, which makes this system astrophysically very important. There are only a few similarly evolved A-type stars among known eclipsing binaries. Future systematic observations and careful analyses can provide very stringent tests for the stellar evolutionary theory. △ Less
Submitted 11 May, 2011; originally announced May 2011.
-
A study of asteroid pole-latitude distribution based on an extended set of shape models derived by the lightcurve inversion method
Authors: Josef Hanus, Josef Durech, Miroslav Broz, Brian D. Warner, Frederick Pilcher, Robert Stephens, Julian Oey, Laurent Bernasconi, Silvano Casulli, Raoul Behrend, David Polishook, Tomas Henych, Martin Lehky, Fumi Yoshida, Takashi Ito
Abstract: Tens of thousands of sparse-in-time lightcurves from astrometric projects are publicly available. We investigate these data and use them in the lightcurve inversion method to derive new asteroid models. By having a greater number of models with known physical properties, we can gain a better insight into the nature of individual objects and into the whole asteroid population. We use sparse photome… ▽ More Tens of thousands of sparse-in-time lightcurves from astrometric projects are publicly available. We investigate these data and use them in the lightcurve inversion method to derive new asteroid models. By having a greater number of models with known physical properties, we can gain a better insight into the nature of individual objects and into the whole asteroid population. We use sparse photometry from selected observatories from the AstDyS database, either alone or in combination with dense lightcurves, to determine new asteroid models by the lightcurve inversion method. We present 80 new asteroid models derived from combined data sets where sparse photometry is taken from the AstDyS database and dense lightcurves are from the Uppsala Asteroid Photometric Catalogue (UAPC) and from several individual observers. For 18 asteroids, we present updated shape solutions based on new photometric data. The addition of new models increases the total number of models derived by the lightcurve inversion method to ~200. We also present a simple statistical analysis of physical properties of asteroids where we look for possible correlations between various physical parameters with an emphasis on the spin vector. We present the observed and de-biased distributions of ecliptic latitudes with respect to different size ranges of asteroids as well as a simple theoretical model of the latitude distribution and then compare its predictions with the observed distributions. From this analysis we find that the latitude distribution of small asteroids (D < 30 km) is clustered towards ecliptic poles and can be explained by the YORP thermal effect while the latitude distribution of larger asteroids (D > 60 km) exhibits an evident excess of prograde rotators, probably of primordial origin. △ Less
Submitted 20 April, 2011; originally announced April 2011.
