Origin of orbits of secondaries in the discovered trans-Neptunian binaries

The dependences of inclinations of orbits of secondaries in the discovered trans-Neptunian binaries on the distance between the primary and the secondary, on the eccentricity of orbits of the secondary around the primary, on the ratio of diameters of the secondary and the primary, and on the elements of heliocentric orbits of these binaries are studied. These dependences are interpreted using the model of formation of a satellite system in a collision of two rarefied condensations composed of dust and/or objects less than 1 m in diameter. It is assumed in this model that a satellite system forms in the process of compression of a condensation produced in such a collision. The model of formation of a satellite system in a collision of two condensations agrees with the results of observations: according to observational data, approximately 40% of trans-Neptunian binaries have a negative angular momentum relative to their centers of mass.     

Formation of Embryos of the Earth and the Moon from a Common Rarefied Condensation and Their Subsequent Growth

Embryos of the Moon and the Earth may have formed as a result of contraction of a common parental rarefied condensation. The required angular momentum of this condensation could largely be acquired in a collision of two rarefied condensations producing the parental condensation. With the subsequent growth of embryos of the Moon and the Earth taken into account, the total mass of as-formed embryos needed to reach the current angular momentum of the Earth–Moon system could be below 0.01 of the Earth mass. For the low lunar iron abundance to be reproduced with the growth of originally iron-depleted embryos of the Moon and the Earth just by the accretion of planetesimals, the mass of the lunar embryo should have increased by a factor of 1.3 at the most. The maximum increase in the mass of the Earth embryo due to the accumulation of planetesimals in a gas-free medium is then threefold, and the current terrestrial iron abundance is not attained. If the embryos are assumed to have grown just by accumulating solid planetesimals (without the ejection of matter from the embryos), it is hard to reproduce the current lunar and terrestrial iron abundances at any initial abundance in the embryos. For the current lunar iron abundance to be reproduced, the amount of matter ejected from the Earth embryo and infalling onto the Moon embryo should have been an order of magnitude larger than the sum of the overall mass of planetesimals infalling directly on the Moon embryo and the initial mass of the Moon embryo, which had formed from the parental condensation, if the original embryo had the same iron abundance as the planetesimals. The greater part of matter incorporated into the Moon embryo could be ejected from the Earth in its multiple collisions with planetesimals (and smaller bodies).     

Estimating the Parameters of Collisions between Fractal Dust Clusters in a Gas–Dust Protoplanetary Disk

Studying the origin and evolution of the Solar system is among the fundamental problems of modern natural science. This problem is interdisciplinary and requires the development of mathematical models for the physical structure and evolution of a gas–dust accretion disk from the initial stages of its formation to the formation of a planetary system. One of the key problems is the formation and growth of bodies in a protoplanetary disk, the basis for which is a study of the collisional processes of the solidbody component. We have performed a parametric analysis of the cluster–cluster collision processes occurring in a protoplanetary disk within the model of permeable particles being developed by us. The outcome of such collisions is shown to be affected significantly by the topological properties of colliding dust clusters with a fractal internal structure. The results of our parametric analysis show that for sufficiently “dense” fractal dust clusters, at low relative collision velocities, there exists a range in which the colliding clusters bounce. At the same time, for “porous” fractal clusters the bounce is impossible for any sets of collision parameters. As the relative collision velocities increase, the cluster coalescence processes begin to dominate due to a rearrangement of the fractal structure in the contact zone. However, as the kinetic energy of collisions increases further, a critical threshold is reached beyond which the collision energy exceeds the particle binding energy in clusters and the fractal dust cluster destruction processes are switched on during collisions. Thus, our parametric analysis imposes quite definite constraints on the dynamics and chronology of the evolution processes during the formation of primordial solid bodies and planetesimals. The proposed approach and the results obtained are fairly realistic and open prospects for more comprehensive model studies of the initial evolutionary phase of a protoplanetary disk.     

The satellites of Neptune and the origin of Pluto

Pluto and the chaotic satellite system of Neptune may have originated from a single encounter of Neptune with a massive solar system body. A series of numerical experiments has been carried out to try to set limits on the circumstances of such an encounter. These experiments show that orbits very much like those of Pluto, Triton, and Nereid can result from a single close encounter of such a body with Neptune. The implied mass range and encounter velocities limit the source of the encountering body to a former trans-Neptunian planet in the 2- to 5-Earth-mass range.     

Dynamics of rotating triple systems

We investigate the dynamical evolution of 100 000 rotating triple systems with equal-mass components. The system rotation is specified by the parameter ω=?c²E, where c and E are the angular momentum and total energy of the triple system, respectively. We consider ω=0.1,1, 2, 4, 6 and study 20 000 triple systems with randomly specified coordinates and velocities of the bodies for each ω. We consider two methods for specifying initial conditions: with and without a hierarchical structure at the beginning of the evolution. The evolution of each system is traced until the escape of one of the bodies or until the critical time equal to 1000 mean system crossing times. For each set of initial conditions, we computed parameters of the final motions: orbital parameters for the final binary and the escaping body. We analyze variations in the statistical characteristics of the distributions of these parameters with ω. The mean disruption time of triple systems and the fraction of the systems that have not been disrupted in 1000 mean crossing times increase with ω. The final binaries become, on average, wider at larger angular momenta. The distribution of their eccentricities does not depend on ω and generally agrees with the theoretical law f(e)=2e. The velocities of the escaping bodies, on average, decrease with increasing angular momentum of the triple system. The fraction of the angles between the escaping-body velocity vector and the triple-system angular momentum close to 90° increases with ω. Escapes in the directions opposite to rotation and prograde motions dominate at small and large angular momenta, respectively. For slowly rotating systems, the angular momentum during their disruption is, on average, evenly divided between the escaping body and the final binary, whereas in rapidly rotating systems, about 80% of the angular momentum is carried away by the escaping component. We compare our numerical simulations with the statistical theory of triple-system disruption.     

Problems of origin of asteroids and comets

Abstract— Various hypotheses of the origin of asteroids and comets are briefly discussed. Interaction of planetesimals in the asteroid zone (AZ) with the gas, their perturbations by proto-Jupiter, and sweeping them out by more massive Jupiter zone bodies when they penetrated the AZ are considered. If the gas was turbulent, it could prevent a settling of dust particles to the equatorial plane of the disk and formation of dust condensations due to gravitational instability. Then particles grew by sticking upon collision. Gas moved radially due to turbulent viscosity and its dissipation. Small particles moved more-or-less together with the gas. As a result of gas drag, larger particles and bodies moved relative to the gas in the direction of increasing gas pressure. Gas would remove much of the solid material from the AZ if most bodies larger than a few km disintegrated by collisions into fragments smaller than a few tens of meters. Most of these fragments would then move into the Martian zone, and the small mass of Mars would have no explanation. Resonant perturbations of asteroids by Jupiter are discussed. In the model of a small mass disk they could scan through the asteroid belt due to changes in Jupiter's distance from the Sun that occurred when this planet accreted the gas and ejected the bodies from the solar system. Such a scanning considerably accelerated the removal of asteroids from the AZ. Massive Jupiter zone bodies with large orbital eccentricities that crossed the AZ were probably efficient at sweeping out bodies. Larger bodies increased the random velocities of the remaining asteroids at close encounters to the present values ～ 5 km/s. Restrictions on the runaway growth of giant planets, on the relative velocities of bodies and the disk surface density that follow from the consideration of the origin of the asteroid belt and the cometary cloud are considered.     

The origin of the density distribution of disc galaxies: a new problem for the standard model of disc formation

We present new models for the formation of disc galaxies that improve upon previous models by following the detailed accretion and cooling of the baryonic mass, and by using realistic distributions of specific angular momentum. Under the assumption of detailed angular momentum conservation, the discs that form have density distributions that are more centrally concentrated than an exponential. We examine the influence of star formation, bulge formation, and feedback on the outcome of the surface brightness distributions of the stars. Low angular momentum haloes yield disc galaxies with a significant bulge component and with a stellar disc that is close to exponential, in good agreement with observations. High angular momentum haloes, on the other hand, produce stellar discs that are much more concentrated than an exponential, in clear conflict with observations. At large radii, the models reveal distinct truncation radii in both the stars and the cold gas. The stellar truncation radii result from our implementation of star formation threshold densities, and are in excellent agreement with observations. The truncation radii in the density distribution of the cold gas reflect the maximum specific angular momentum of the gas that has cooled. We find that these truncation radii occur at H  i surface densities of roughly 1 M_⊙ pc⁻², in conflict with observations. We examine various modifications to our models, including feedback, viscosity, and dark matter haloes with constant-density cores, but show that the models consistently fail to produce bulge less discs with exponential surface brightness profiles. This signals a new problem for the standard model of disc formation: if the baryonic component of the protogalaxies out of which disc galaxies form has the same angular momentum distribution as the dark matter, discs are too compact.     

Systems of colliding bodies in a gravitational field: Evolution and role of the internal rotations

To emphasize the rotational effects of a simple friction between colliding bodies in a keplerian field we investigate numerically the evolution of the rotational energies in a three dimensional system of spherical particles interacting through inelastic collisions in a deterministic model. All the particles are made of the same material but they possibly have different sizes. Each collision reduces the relative surface velocity and there are exchanges between orbital energy and rotational energy. Our results are compared with some previous papers and our aim is to supply other probabilists models with simple basic references about mean dynamical properties.The rotational energy of the colliding bodies tends to reach an equilibrium state that depends only on the rate of energy loss in the collision process. Internal rotations prevent the complete flattening of the system. With this model, light and small particles spin faster than the massive and big ones. We observe an excess of prograde rotations on counterclockwise orbits. The ratio of rotational and orbital energies is % MathType!MTEF!2!1!+-% feaafeart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr% 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq-Jc9% vqaqpepm0xbba9pwe9Q8fs0-yqaqpepae9pg0FirpepeKkFr0xfr-x% fr-xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamyramaaBa% aaleaacaWGYbaabeaakiaac+cacaWGfbWaaSbaaSqaaiaadUgaaeqa% aOGaeyisISRaaGymaiaaicdadaahaaWcbeqaaiabgkHiTiaaiodaaa% aaaa!3F83!\[E_r /E_k \approx 10^{ - 3} \] while the ratio of corresponding mean angular velocities is % MathType!MTEF!2!1!+-% feaafeart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr% 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq-Jc9% vqaqpepm0xbba9pwe9Q8fs0-yqaqpepae9pg0FirpepeKkFr0xfr-x% fr-xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaWaaaWaaeaacq% aHjpWDaiaawMYicaGLQmcacaGGVaWaaaWaaeaacqGHPoWvaiaawMYi% caGLQmcacqGHijYUcaaIYaaaaa!4008!\[\left\langle \omega \right\rangle /\left\langle \Omega \right\rangle \approx 2\] These values depends strongly on the dimensional scale of the model.     

Analytical formulation of impulsive collision avoidance dynamics

The paper deals with the problem of impulsive collision avoidance between two colliding objects in three dimensions and assuming elliptical Keplerian orbits. Closed-form analytical expressions are provided that accurately predict the relative dynamics of the two bodies in the encounter b-plane following an impulsive delta-V manoeuvre performed by one object at a given orbit location prior to the impact and with a generic three-dimensional orientation. After verifying the accuracy of the analytical expressions for different orbital eccentricities and encounter geometries the manoeuvre direction that maximises the miss distance is obtained numerically as a function of the arc length separation between the manoeuvre point and the predicted collision point. The provided formulas can be used for high-accuracy instantaneous estimation of the outcome of a generic impulsive collision avoidance manoeuvre and its optimisation.     

Masses of the Main Asteroid Belt and the Kuiper Belt from the Motions of Planets and Spacecraft

Dynamicalmass estimates for the main asteroid belt and the trans-Neptunian Kuiper belt have been found from their gravitational influence on the motion of planets. Discrete rotating models consisting ofmovingmaterial points have been used tomodel the total attraction fromsmall or as yet undetected bodies of the belts. The masses of the model belts have been included in the set of parameters being refined and determined and have been obtained by processing more than 800 thousand modern positional observations of planets and spacecraft. We have processed the observations and determined the parameters based on the new EPM2017 version of the IAA RAS planetary ephemerides. The large observed radial extent of the belts (more than 1.2 AU for the main belt and more than 8 AU for the Kuiper belt) and the concentration of bodies in the Kuiper belt at a distance of about 44 AU found from observations have been taken into account in the discrete models. We have also used individual mass estimates for large bodies of the belts as well as for objects that spacecraft have approached and for bodies with satellites. Our mass estimate for the main asteroid belt is (4.008 ± 0.029) × 10^?4/m_⊕ (3σ). The bulk of the Kuiper belt objects are in the ring zone from 39.4 to 47.8 AU. The estimate of its total mass together with the mass of the 31 largest trans-Neptunian Kuiper belt objects is (1.97 ± 0.30) × 10^?2m_⊕ (3σ), which exceeds the mass of the main asteroid belt almost by a factor of 50. The mass of the 31 largest trans-Neptunian objects (TNOs) is only about 40% of the total one.     

Satellite-sized planetesimals and lunar origin

Exploratory calculations using accretionary theory are made to demonstrate plausible sizes of second-largest, third-largest, etc., bodies at the close of planet formation in heliocentric orbits near the planets, assuming asteroid-like size distributions at the start of the calculation. Many satellite-sized bodies are found to be available for capture, cratering, or collisional fragmentation. In the case of Earth-sized planets, the models suggest second-largest bodies of 500 to 3000 km radius, and tens of bodies larger than 100 km radius. Many of these interact with the planet before suffering any fragmentation events with each other. Collision of a large body with Earth could eject iron-deficient crust and upper mantle material, forming a cloud of refractory, volatile-poor dust that could form the Moon. Other satellite systems may have been affected by major capture or collision events of chance character.     

Numerical simulations of collisions and gravitational encounters in systems of non-identical particles

Numerical simulations of 200 mutually colliding non-identical particles indicate that the equipartition of random kinetic energy is possible only in systems having a narrow distribution of particle masses. Otherwise the random energy is concentrated on heavy particles. The form of the velocity distribution versus particle mass depends also on the elastic properties of the particles, and on the relative importance of the particle size. If the coefficient of restitution is a weakly decreasing function of impact velocity, a large difference in the equilibrium velocities of largest and smallest particles is possible. On the other hand, if the elasticity drops to a low level even in the small velocity regime, the dispersion of velocities is maintained by finite size and differential rotation, and the velocities of smallest particles are, at most, slightly larger than those of the largest ones. The results of simulations are consistent with the predictions of the collisional theory of non-identical particles (Hämeen-Anttila, 1984). The application to Saturn's rings indicates that the geometric thickness of cm-sized particles is of the order of 50 m in the rarefied regions of the rings. Without the gravitational encounters a thickness of about 30 m is derived. These estimations are made by using the latest measurements (Bridges et al., 1984) for the restitution coefficient of icy particles.     

The cosmological dependence of galactic specific angular momenta

Hydrodynamical simulations of galaxy formation in spatially flat cold dark matter (CDM) cosmologies with and without a cosmological constant (Λ) are described. A simple star formation algorithm is employed and radiative cooling is allowed only after redshift z =1 so that enough hot gas is available to form large, rapidly rotating stellar discs if angular momentum is approximately conserved during collapse. The specific angular momenta of the final galaxies are found to be sensitive to the assumed background cosmology. This dependence arises from the different angular momenta contained in the haloes at the epoch when the gas begins to collapse and the inhomogeneity of the subsequent halo evolution. In the Λ-dominated cosmology, the ratio of stellar specific angular momentum to that of the dark matter halo (measured at the virial radius) has a median value of ∼0.24 at z =0. The corresponding quantity for the Λ=0 cosmology is over three times lower. It is concluded that the observed frequency and angular momenta of disc galaxies pose significant problems for spatially flat CDM models with Λ=0 but may be consistent with a Λ-dominated CDM universe.     

Can newly-forming stars act like Mini-QSOs?

The birth of stars takes place inside dense molecular clouds and is therefore difficult to observe with optical telescopes. Yet some stars form at the edges of clouds, and combined radio, infrared, optical and X-ray observations have recently revealed a great deal of structure: over a wide range of luminosities one observes high mass-outflow rates from pre-Main-Sequence and T Tauri-stars, with wind momenta exceeding the radiation momenta by large factors. The mass flows take the shape of two highly supersonic jets, perpendicular to a circumstellar disk. The jets have knotlike condensations, show strong linear polarization along the flow direction, and are often seen to reconverge. They resemble the twin-jets from the nuclei of active galaxies and may be driven by a similar type of engine. They may even hold a clue to the problem of how stars like our Sun got rid of the enormous angular momentum of their progenitor cloud. I propose that in both phenomena, the central engine is the maximally rotating core of a massive disk which produces a pair of thin, antipodal, magnetized relativistic jets.     

Planet formation: Mechanism of early growth

Experiments in vacuum (approx. 0.5 to 1 mbar) and in air quantify mechanics of collisions, rebound, and fragmentation at low velocities (1–50 m/sec), under the conditions usually postulated for the preplanetary environment in the primitive solar nebula. Such collisions have been little studied experimentally. Contrary to widespread assumptions, accretionary growth of the largest meteoroid- and asteroid-sized bodies in a given swarm results spontaneously from the simple mechanics of these collisions, without other ad hoc sticking mechanisms. The smaller bodies in the swarm are less likely to grow. Granular surfaces form, either by gravitational collapse of dust swarms or by rapid formation of regolith surfaces on solid planetesimals; these surfaces strongly promote further growth by retarding rebound. Growth of large bodies increases modal collision velocities, causing fragmentation of smaller bodies and eventual production of interstellar dust as a by-product planetesimal interactions.     

Matrix method for the Liapunov stability analysis of cyclic discrete mechanical systems

A matrix formalism is developed for the purpose of facilitating the Liapunov stability analysis of discrete, holonomic, mechanical systems with cyclic coordinates and with the Hamiltonian free of explicit time dependence. Matrix expressions are developmed for the kinetic energy, the Routhian, the Hamiltonian, and the quadratic approximation of the dynamic potential energy, with cyclic coordinates, cyclic-coordinate velocities, and cyclic-coordinate generalized momenta not explicitly involved in the last of these functions. The final result is an expression for the quadratic approximation of the dynamic potential energy that is calculated much more readily than by scalar analysis. From the condition for positive-definiteness of this function, Liapunov stability conditions are available. The method is applied to a dual-spin satellite to illustrate the procedure.     

Exploring the 7:4 mean motion resonance—II: Scattering evolutionary paths and resonance sticking

In our preliminary study, we have investigated basic properties and dynamical evolution of classical TNOs around the 7:4 mean motion resonance with Neptune (a∼43.7 AU), motivated by observational evidences that apparently present irregular features near this resonance (see [Lykawka and Mukai, 2005a. Exploring the 7:4 mean motion resonance—I. Dynamical evolution of classical trans-Neptunian objects (TNOs). Space Planet. Sci. 53, 1175-1187]; hereafter “Paper I”). In this paper, we aim to explore the dynamical long-term evolution in the scattered disk (but not its early formation) based on the computer simulations performed in Paper I together with extra computations. Specifically, we integrated the orbital motion of test particles (totalizing a bit more than 10,000) placed around the 7:4 mean motion resonance under the effect of the four giant planets for the age of the Solar System. In order to investigate chaotic diffusion, we also conducted a special simulation with on-line computation of proper elements following tracks in phase space over 4-5 Gyr. We found that: (1) A few percent (1-2%) of the test particles survived in the scattered disk with direct influence of other Neptunian mean motion resonances, indicating that resonance sticking is an extremely common phenomenon and that it helps to enhance scattered objects longevity. (2) In the same region, the so-called extended scattered TNOs are able to form via very long resonance trapping under certain conditions. Namely, if the body spends more than about 80% of its dynamical lifetime trapped in mean motion resonance(s) and there is the action of a k+1 or (k+2)/2 mean motion resonance (e.g., external mean motion resonances with Neptune described as (j+k)/j with j=1 and 2, respectively). According to this hypothetical mechanism, 5-15% of current scattered TNOs would possess thus probably constituting a significant part of the extended scattered disk. (3) Moreover, considering hot orbital initial conditions, it is likely that the trans-Neptunian belt (or Edgeworth-Kuiper belt) has been providing members to the scattered disk, so that scattered TNOs observed today would consist of primordial scattered bodies mixed with TNOs that came from unstable regions of the trans-Neptunian belt in the past.Considering the three points together, our results demonstrated that the scattered disk has been evolving continuously since early times until present.     

Creating an isotopically similar Earth–Moon system with correct angular momentum from a giant impact

The giant impact hypothesis is the dominant theory explaining the formation of our Moon. However, the inability to produce an isotopically similar Earth–Moon system with correct angular momentum has cast a shadow on its validity. Computer-generated impacts have been successful in producing virtual systems that possess many of the observed physical properties. However, addressing the isotopic similarities between the Earth and Moon coupled with correct angular momentum has proven to be challenging. Equilibration and evection resonance have been proposed as means of reconciling the models. In the summer of 2013, the Royal Society called a meeting solely to discuss the formation of the Moon. In this meeting, evection resonance and equilibration were both questioned as viable means of removing the deficiencies from giant impact models. The main concerns were that models were multi-staged and too complex. We present here initial impact conditions that produce an isotopically similar Earth–Moon system with correct angular momentum. This is done in a single-staged simulation. The initial parameters are straightforward and the results evolve solely from the impact. This was accomplished by colliding two roughly half-Earth-sized impactors, rotating in approximately the same plane in a high-energy, off-centered impact, where both impactors spin into the collision.     

Numerical investigation of collision orbits of lunar satellites

In the present study an investigation of the collision orbits of natural satellites of the Moon (considered to be of finite dimensions) is developed, and the tendency of natural satellites of the Moon to collide on the visible or the far side of the Moon is studied. The collision course of the satellite is studied up to its impact on the lunar surface for perturbations of its initial orbit arbitrarily induced, for example, by the explosion of a meteorite. Several initial conditions regarding the position of the satellite to collide with the Moon on its near (visible) or far (invisible) side is examined in connection to the initial conditions and the direction of the motion of the satellite. The distribution of the lunar craters-originating impact of lunar satellites or celestial bodies which followed a course around the Moon and lost their stability - is examined. First, we consider the planar motion of the natural satellite and its collision on the Moon's surface without the presence of the Earth and Sun. The initial velocities of the satellite are determined in such a way so its impact on the lunar surface takes place on the visible side of the Moon. Then, we continue imparting these velocities to the satellite, but now in the presence of the Earth and Sun; and study the forementioned impacts of the satellites but now in the Earth-Moon-Satellite system influenced also by the Sun. The initial distances of the satellite are taken as the distances which have been used to compute periodic orbits in the planar restricted three-body problem (cf. Gousidou-Koutita, 1980) and its direction takes different angles with the x-axis (Earth-Moon axis). Finally, we summarise the tendency of the satellite's impact on the visible or invisible side of the Moon.