In this work, we investigate the rotational dynamics of the ginger-shaped near-Earth asteroid 4179 Toutatis, which was closely observed by Change-2 at a distance of $770pm120~$ meters from the asteroids surface during the outbound flyby citep{Huang20
13} on 13 December 2012. A sequence of high-resolution images was acquired during the flyby mission. In combination with ground-based radar observations collected over the last two decades, we analyze these flyby images and determine the orientation of the asteroid at the flyby epoch. The 3-1-3 Euler angles of the conversion matrix from the J2000 ecliptic coordinate system to the body-fixed frame are evaluated to be $-20.1^circpm1^circ$, $27.6^circpm1^circ$ and $42.2^circpm1^circ$, respectively. The least-squares method is utilized to determine the rotational parameters and spin state of Toutatis. The characteristics of the spin-state parameters and angular momentum variations are extensively studied using numerical simulations, which confirm those reported by citet{Takahashi2013}. The large amplitude of Toutatis precession is assumed to be responsible for its tumbling attitude as observed from Earth. Toutatis angular momentum orientation is determined to be described by $lambda_{H}=180.2^{+0.2^circ}_{-0.3^circ}$ and $beta_{H}=-54.75^{+0.15^circ}_{-0.10^circ}$, implying that it has remained nearly unchanged for two decades. Furthermore, using Fourier analysis to explore the change in the orientation of Toutatis axes, we reveal that the two rotational periods are 5.38 and 7.40 days, respectively, consistent with the results of the former investigation. Hence, our investigation provides a clear understanding of the state of the rotational dynamics of Toutatis.
On 13 December 2012, Change-2 conducted a successful flyby of the near-Earth asteroid 4179 Toutatis at a closest distance of 770 $pm$ 120 meters from the asteroids surface. The highest-resolution image, with a resolution of better than 3 meters, reve
als new discoveries on the asteroid, e.g., a giant basin at the big end, a sharply perpendicular silhouette near the neck region, and direct evidence of boulders and regolith, which suggests that Toutatis may bear a rubble-pile structure. Toutatis maximum physical length and width are (4.75 $times$ 1.95 km) $pm$10$%$, respectively, and the direction of the +$z$ axis is estimated to be (250$pm$5$^circ$, 63$pm$5$^circ$) with respect to the J2000 ecliptic coordinate system. The bifurcated configuration is indicative of a contact binary origin for Toutatis, which is composed of two lobes (head and body). Change-2 observations have significantly improved our understanding of the characteristics, formation, and evolution of asteroids in general.
We investigate the formation of terrestrial planets in the late stage of planetary formation using two-planet model. At that time, the protostar has formed for about 3 Myr and the gas disk has dissipated. In the model, the perturbations from Jupiter
and Saturn are considered. We also consider variations of the mass of outer planet, and the initial eccentricities and inclinations of embryos and planetesimals. Our results show that, terrestrial planets are formed in 50 Myr, and the accretion rate is about $60% - 80%$. In each simulation, 3 - 4 terrestrial planets are formed inside Jupiter with masses of $0.15 - 3.6 M_{oplus}$. In the $0.5 - 4$ AU, when the eccentricities of planetesimals are excited, planetesimals are able to accrete material from wide radial direction. The plenty of water material of the terrestrial planet in the Habitable Zone may be transferred from the farther places by this mechanism. Accretion could also happen a few times between two major planets only if the outer planet has a moderate mass and the small terrestrial planet could survive at some resonances over time scale of $10^8$ yr. In one of our simulations, com-mensurability of the orbital periods of planets is very common. Moreover, a librating-circulating 3:2 configuration of mean motion resonance is found.
We perform numerical simulations to study the Habitable zones (HZs) and dynamical structure for Earth-mass planets in multiple planetary systems. For example, in the HD 69830 system, we extensively explore the planetary configuration of three Neptune
-mass companions with one massive terrestrial planet residing in 0.07 AU $leq a leq$ 1.20 AU, to examine the asteroid structure in this system. We underline that there are stable zones of at least $10^5$ yr for low-mass terrestrial planets locating between 0.3 and 0.5 AU, and 0.8 and 1.2 AU with final eccentricities of $e < 0.20$. Moreover, we also find that the accumulation or depletion of the asteroid belt are also shaped by orbital resonances of the outer planets, for example, the asteroidal gaps at 2:1 and 3:2 mean motion resonances (MMRs) with Planet C, and 5:2 and 1:2 MMRs with Planet D. In a dynamical sense, the proper candidate regions for the existence of the potential terrestrial planets or HZs are 0.35 AU $< a < $ 0.50 AU, and 0.80 AU $< a < $ 1.00 AU for relatively low eccentricities, which makes sense to have the possible asteroidal structure in this system.
We perform numerical simulations to study the secular orbital evolution and dynamical structure in the quintuplet planetary system 55 Cancri with the self-consistent orbital solutions by Fischer and coworkers (2008). In the simulations, we show that
this system can be stable at least for $10^{8}$ yr. In addition, we extensively investigate the planetary configuration of four outer companions with one terrestrial planet in the wide region of 0.790 AU $leq a leq $ 5.900 AU to examine the existence of potential asteroid structure and Habitable Zones (HZs). We show that there are unstable regions for the orbits about 4:1, 3:1 and 5:2 mean motion resonances (MMRs) with the outermost planet in the system, and several stable orbits can remain at 3:2 and 1:1 MMRs, which is resemblance to the asteroidal belt in solar system. In a dynamical point, the proper candidate HZs for the existence of more potential terrestrial planets reside in the wide area between 1.0 AU and 2.3 AU for relatively low eccentricities.
We use a multi-dimensional hydrodynamics code to study the gravitational interaction between an embedded planet and a protoplanetary disk with emphasis on the generation of vortensity (Potential Vorticity or PV) through a Baroclinic Instability. We s
how that the generation of PV is very common and effective in non-barotropic disks through the Baroclinic Instability, especially within the coorbital region. Our results also complement previous work that non-axisymmetric Rossby-Wave Instabilities (RWIs, Lovelace et al. 1999) are likely to develop at local minima of PV distribution that are generated by the interaction between a planet and an inviscid barotropic disk. The development of RWIs results in non-axisymmetric density blobs, which exert stronger torques onto the planet when they move to the vicinity of the planet. Hence, large amplitude oscillations are introduced to the time behavior of the total torque acted on the planet by the disk. In current simulations, RWIs do not change the overall picture of inward orbital migration but cause a non-monotonic behavior to the migration speed. As a side effect, RWIs also introduce interesting structures into the disk. These structures may help the formation of Earth-like planets in the Habitable Zone or Hot Earths interior to a close-in giant planet.
