ترغب بنشر مسار تعليمي؟ اضغط هنا

A close-encounter method for simulating the dynamics of planetesimals

91   0   0.0 ( 0 )
 نشر من قبل Sebastian Lorek
 تاريخ النشر 2020
  مجال البحث فيزياء
والبحث باللغة English




اسأل ChatGPT حول البحث

The dynamics of planetesimals plays an important role in planet formation, because their velocity distribution sets the growth rate to larger bodies. When planetesimals form in protoplanetary discs, their orbits are nearly circular and planar due to the effect of gas drag. However, mutual close encounters of the planetesimals increase eccentricities and inclinations until an equilibrium between stirring and damping is reached. After disc dissipation, there is no more gas drag and mutual close encounters as well as encounters with planets stir the orbits again. The high number of planetesimals in protoplanetary discs renders it difficult to simulate their dynamics by means of direct N-body simulations of planet formation. Therefore, we developed a novel method for the dynamical evolution of planetesimals that is based on following close encounters between planetesimal-mass bodies and gravitational stirring by planet-mass bodies. To separate the orbital motion from the close encounters, we employ a Hamiltonian splitting scheme as used in symplectic N-body integrators. Close encounters are identified using a cell algorithm with linear scaling in the number of bodies. A grouping algorithm is used to create small groups of interacting bodies which are integrated separately. Our method allows simulating a high number of planetesimals interacting through gravity and collisions with low computational cost. The typical computational time is of the order of minutes or hours, up to a few days for more complex simulations, as compared to several hours or even weeks for the same setup with full N-body. The dynamical evolution of the bodies is sufficiently well reproduced. This will make it possible to study the growth of planetesimals through collisions and pebble accretion coupled to their dynamics for a much higher number of bodies than previously accessible with full N-body simulations.



قيم البحث

اقرأ أيضاً

We have used (a) HST ACS imaging and STIS spectroscopy, (b) ground-based PIONIER/VLT long-baseline interferometry, and (c) ground-based spectroscopy from different instruments to study the orbit of the extreme multiple system HD 93 129 Aa,Ab, which i s composed of (at least) two very massive stars in a long-period orbit with e>0.92 that will pass through periastron in 2017/2018. In several ways, the system is an eta Car precursor. Around the time of periastron passage the two very strong winds will collide and generate an outburst of non-thermal hard X-ray emission without precedent in an O+O binary since astronomers have been able to observe above Earths atmosphere. A coordinated multiwavelength monitoring in the next two years will enable a breakthrough understanding of the wind interactions in such extreme close encounters. Furthermore, we have found evidence that HD 93 129 Aa may be a binary system itself. In that case, we could witness a three-body interaction that may yield a runaway star or a stellar collision close to or shortly after the periastron passage. Either of those outcomes would be unprecedented, as they are predicted to be low-frequency events in the Milky Way.
For problems in astrophysics, planetary science and beyond, numerical simulations are often limited to simulating fewer particles than in the real system. To model collisions, the simulated particles (aka superparticles) need to be inflated to repres ent a collectively large collisional cross section of real particles. Here we develop a superparticle-based method that replicates the kinetic energy loss during real-world collisions, implement it in an $N$-body code and test it. The tests provide interesting insights into dynamics of self gravitating collisional systems. They show how particle systems evolve over several free fall timescales to form central concentrations and equilibrated outer shells. The superparticle method can be extended to account for the accretional growth of objects during inelastic mergers.
The discovery of 2012VP113 initiated the debate on the origin of the Sedna family of planetesimals in orbit around the Sun. Sednitos roam the outer regions of the Solar System between the Egeworth--Kuiper belt and the Oort cloud, in extraordinary wid e (a>150au) orbits with a large perihelion distance of q>30au compared to the Earths (a=1au and eccentricity e=(1-q/a) ~ 0.0167 or q=1au). This population is composed of a dozen objects, which we consider a family because they have similar perihelion distance and inclination with respect to the ecliptic i=10--30deg. They also have similar argument of perihelion omega=340+/-55deg. There is no ready explanation for their origin. Here we show that these orbital parameters are typical for a captured population from the planetesimal disk of another star.Assuming the orbital elements of Sednitos have not changed since they acquired their orbits, we reconstruct the encounter that led to their capture. We conclude that they might have been captured in a near miss with a 1.8MSun star that impacted the Sun at ~340au at an inclination with respect to the ecliptic of 17--34deg with a relative velocity at infinity of ~4.3km/s. We predict that the Sednitos-region is populated by 930 planetesimals and the inner Oort cloud acquired ~440 planetesimals through the same encounter.
Simulations of exoplanet albedo profiles are key to planning and interpreting future direct imaging observations. In this paper we demonstrate the use of the Planetary Spectrum Generator to produce simulations of reflected light exoplanet spectra. We use PSG to examine multiple issues relevant to all models of directly imaged exoplanet spectra and to produce sample spectra of the bright, nearby exoplanet $upsilon$ Andromedae d, a potential direct imaging target for next-generation facilities. We introduce a new, fast, and accurate subsampling technique that enables calculations of disk-integrated spectra one order of magnitude faster than Chebyshev-Gauss sampling for moderate- to high-resolution sampling. Using this method and a first-principles-derived atmosphere for $upsilon$ And d, we simulate phase-dependent spectra for a variety of different potential atmospheric configurations. The simulated spectra for $upsilon$ And d inclu
We present the Superparticle Model/Algorithm for Collisions in Kuiper belts and debris disks (SMACK), a new method for simultaneously modeling, in 3-D, the collisional and dynamical evolution of planetesimals in a debris disk with planets. SMACK can simulate azimuthal asymmetries and how these asymmetries evolve over time. We show that SMACK is stable to numerical viscosity and numerical heating over 10^7 yr, and that it can reproduce analytic models of disk evolution. We use SMACK to model the evolution of a debris ring containing a planet on an eccentric orbit. Differential precession creates a spiral structure as the ring evolves, but collisions subsequently break up the spiral, leaving a narrower eccentric ring.
التعليقات
جاري جلب التعليقات جاري جلب التعليقات
سجل دخول لتتمكن من متابعة معايير البحث التي قمت باختيارها
mircosoft-partner

هل ترغب بارسال اشعارات عن اخر التحديثات في شمرا-اكاديميا