• ### SPH/N-body simulations of small (D = 10 km) asteroidal breakups and improved parametric relations for Monte-Carlo collisional models(1803.10666)

March 28, 2018 astro-ph.EP
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.
• ### Improved model of the triple system V746 Cas that has a bipolar magnetic field associated with the tertiary(1709.06340)

Sept. 19, 2017 astro-ph.SR
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.
• ### An advanced N-body model for interacting multiple stellar systems(1706.04469)

June 14, 2017 astro-ph.SR
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 $\chi^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).
• The goal of this report is to give a comprehensive overview of the rich field of forward physics, with a special attention to the topics that can be studied at the LHC. The report starts presenting a selection of the Monte Carlo simulation tools currently available, chapter 2, then enters the rich phenomenology of QCD at low, chapter 3, and high, chapter 4, momentum transfer, while the unique scattering conditions of central exclusive production are analyzed in chapter 5. The last two experimental topics, Cosmic Ray and Heavy Ion physics are presented in the chapter 6 and 7 respectively. Chapter 8 is dedicated to the BFKL dynamics, multiparton interactions, and saturation. The report ends with an overview of the forward detectors at LHC. Each chapter is correlated with a comprehensive bibliography, attempting to provide to the interested reader with a wide opportunity for further studies.
• ### The origin of long-lived asteroids in the 2:1 mean-motion resonance with Jupiter(1505.04329)

May 16, 2015 astro-ph.EP
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.
• ### The thermal emission from boulders on (25143) Itokawa and general implications for the YORP effect(1504.00566)

April 2, 2015 astro-ph.EP
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.
• ### Identification and Dynamical Properties of Asteroid Families(1502.01628)

Feb. 5, 2015 astro-ph.EP
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.
• ### An anisotropic distribution of spin vectors in asteroid families(1309.4296)

Sept. 17, 2013 astro-ph.EP
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, \beta) space: (i) larger asteroids are situated in the proximity of the center of the family; (ii) asteroids with \beta>0{\deg} are usually found to the right from the family center; (iii) on the other hand, asteroids with \beta<0{\deg} to the left from the center; (iv) majority of asteroids have large pole-ecliptic latitudes (|\beta|\gtrsim 30{\deg}); and finally (v) some families have a statistically significant excess of asteroids with \beta>0{\deg} or \beta<0{\deg}. 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).
• ### Asteroids' physical models from combined dense and sparse photometry and scaling of the YORP effect by the observed obliquity distribution(1301.6943)

Jan. 29, 2013 astro-ph.EP
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.
• ### Constraining the cometary flux through the asteroid belt during the late heavy bombardment(1301.6221)

Jan. 26, 2013 astro-ph.EP
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.
• The ALICE Collaboration has measured inclusive J/psi production in pp collisions at a center of mass energy sqrt(s)=2.76 TeV at the LHC. The results presented in this Letter refer to the rapidity ranges |y|<0.9 and 2.5<y<4 and have been obtained by measuring the electron and muon pair decay channels, respectively. The integrated luminosities for the two channels are L^e_int=1.1 nb^-1 and L^mu_int=19.9 nb^-1, and the corresponding signal statistics are N_J/psi^e+e-=59 +/- 14 and N_J/psi^mu+mu-=1364 +/- 53. We present dsigma_J/psi/dy for the two rapidity regions under study and, for the forward-y range, d^2sigma_J/psi/dydp_t in the transverse momentum domain 0<p_t<8 GeV/c. The results are compared with previously published results at sqrt(s)=7 TeV and with theoretical calculations.
• ### Did the Hilda collisional family form during the late heavy bombardment?(1109.1114)

Sept. 6, 2011 astro-ph.EP
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.
• ### V2368 Oph: An eclipsing and double-lined spectroscopic binary used as a photometric comparison star for U Oph(1105.2162)

May 11, 2011 astro-ph.SR
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.
• ### A unified solution for the orbit and light-time effect in the V505 Sgr system(1002.1145)

Feb. 5, 2010 astro-ph.SR
The multiple system V505 Sagittarii is composed of at least three stars: a compact eclipsing pair and a distant component, which orbit is measured directly using speckle interferometry. In order to explain the observed orbit of the third body in V505 Sagittarii and also other observable quantities, namely the minima timings of the eclipsing binary and two different radial velocities in the spectrum, we thoroughly test a fourth-body hypothesis - a perturbation by a dim, yet-unobserved object. We use an N-body numerical integrator to simulate future and past orbital evolution of 3 or 4 components in this system. We construct a suitable chi^2 metric from all available speckle-interferometry, minima-timings and radial-velocity data and we scan a part of a parameter space to get at least some of allowed solutions. In principle, we are able to explain all observable quantities by a presence of a fourth body, but the resulting likelihood of this hypothesis is very low. We also discuss other theoretical explanations of the minima timings variations. Further observations of the minima timings during the next decade or high-resolution spectroscopic data can significantly constrain the model.