scholarly journals Dynamics of multiple protoplanets embedded in gas and pebble discs and its dependence on Σ and ν parameters

2018 ◽  
Vol 620 ◽  
pp. A157 ◽  
Author(s):  
M. Brož ◽  
O. Chrenko ◽  
D. Nesvorný ◽  
M. Lambrechts

Protoplanets of super-Earth size 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 planetary eccentricities, and prevent resonant captures of migrating planets. In this work, we study populations of embryos that are accreting pebbles under different conditions, by changing the surface density, viscosity, pebble flux, mass, and the number of protoplanets. For modelling, we used the FARGO-THORIN two-dimensional (2D) hydrocode, which incorporates a pebble disc as a second pressure-less fluid, the coupling between the gas and pebbles, and the flux-limited diffusion approximation for radiative transfer. We find that massive embryos embedded in a disc with high surface density (Σ = 990 g cm−2 at 5.2 au) undergo numerous “unsuccessful” two-body encounters that do not lead to a merger. Only when a third protoplanet arrives in the convergence zone do three-body encounters lead to mergers. For a low-viscosity disc (ν = 5 × 1013 cm2 s−1), a massive co-orbital is a possible outcome, for which a pebble isolation develops and the co-orbital is further stabilised. For more massive protoplanets (5 M⊕), 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 co-orbital configuration is disrupted and planets merge. This may be a way to solve the problem that co-orbitals often form in simulations but they are not observed in nature. In contrast, the joint evolution of 120 low-mass protoplanets (0.1 M⊕) 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 disc due to Lindblad resonances, overlap with each other and affect not only a single protoplanet but several in the surrounding area. 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.

Author(s):  
Frédéric S. Masset

Planet migration is the variation over time of a planet’s semimajor axis, leading to either a contraction or an expansion of the orbit. It results from the exchange of energy and angular momentum between the planet and the disk in which it is embedded during its formation and can cause the semimajor axis to change by as much as two orders of magnitude over the disk’s lifetime. The migration of forming protoplanets is an unavoidable process, and it is thought to be a key ingredient for understanding the variety of extrasolar planetary systems. Although migration occurs for protoplanets of all masses, its properties for low-mass planets (those having up to a few Earth masses) differ significantly from those for high-mass planets. The torque that is exerted by the disk on the planet is composed of different contributions. While migration was first thought to be invariably inward, physical processes that are able to halt or even reverse migration were later uncovered, leading to the realization that the migration path of a forming planet has a very sensitive dependence on the underlying disk parameters. There are other processes that go beyond the case of a single planet experiencing smooth migration under the disk’s tide. This is the case of planetary migration in low-viscosity disks, a fashionable research avenue because protoplanetary disks are thought to have very low viscosity, if any, over most of their planet-forming regions. Such a process is generally significantly chaotic and has to be tackled through high-resolution numerical simulations. The migration of several low-mass planets is also is a very fashionable topic, owing to the discovery by the Kepler mission of many multiple extrasolar planetary systems. The orbital properties of these systems suggest that at least some of them have experienced substantial migration. Although there have been many studies to account for the orbital properties of these systems, there is as yet no clear picture of the different processes that shaped them. Finally, some recently unveiled processes could be important for the migration of low-mass planets. One process is aero-resonant migration, in which a swarm of planetesimals subjected to aerodynamic drag push a planet inward when they reach a mean-motion resonance with the planet, while another process is based on so-called thermal torques, which arise when thermal diffusion in the disk is taken into account, or when the planet, heated by accretion, releases heat into the ambient gas.


2020 ◽  
Vol 635 ◽  
pp. A204 ◽  
Author(s):  
S. Ataiee ◽  
W. Kley

Context. The most accurate method for modelling planetary migration and hence the formation of resonant systems is using hydrodynamical simulations. Usually, the force (torque) acting on a planet is calculated using the forces from the gas disc and the star, while the gas accelerations are computed using the pressure gradient, the star, and the planet’s gravity, ignoring its own gravity. For a non-migrating planet the neglect of the disc gravity results in a consistent torque calculation while for a migrating case it is inconsistent. Aims. We aim to study how much this inconsistent torque calculation can affect the final configuration of a two-planet system. We focus on low-mass planets because most of the multi-planetary systems, discovered by the Kepler survey, have masses around ten Earth masses. Methods. Performing hydrodynamical simulations of planet–disc interaction, we measured the torques on non-migrating and migrating planets for various disc masses as well as density and temperature slopes with and without considering the self-gravity of the disc. Using this data, we found a relation that quantifies the inconsistency, used this relation in an N-body code, and performed an extended parameter study modelling the migration of a planetary system with different planet mass ratios and disc surface densities, to investigate the impact of the torque inconsistency on the architecture of the planetary system. Results. Not considering disc self-gravity produces an artificially larger torque on the migrating planet that can result in tighter planetary systems. The deviation of this torque from the correct value is larger in discs with steeper surface density profiles. Conclusions. In hydrodynamical modelling of multi-planetary systems, it is crucial to account for the torque correction, otherwise the results favour more packed systems. We examine two simple correction methods existing in the literature and show that they properly correct this problem.


2020 ◽  
Vol 499 (1) ◽  
pp. 668-680
Author(s):  
Alejandro González-Samaniego ◽  
Enrique Vazquez-Semadeni

ABSTRACT We use two hydrodynamical simulations (with and without photoionizing feedback) of the self-consistent evolution of molecular clouds (MCs) undergoing global hierarchical collapse (GHC), to study the effect of the feedback on the structural and kinematic properties of the gas and the stellar clusters formed in the clouds. During this early stage, the evolution of the two simulations is very similar (implying that the feedback from low-mass stars does not affect the cloud-scale evolution significantly) and the star-forming region accretes faster than it can convert gas into stars, causing the instantaneous measured star formation efficiency (SFE) to remain low even in the absence of significant feedback. Afterwards, the ionizing feedback first destroys the filamentary supply to star-forming hubs and ultimately removes the gas from it, thus first reducing the star formation (SF) and finally halting it. The ionizing feedback also affects the initial kinematics and spatial distribution of the forming stars because the gas being dispersed continues to form stars, which inherit its motion. In the non-feedback simulation, the groups remain highly compact and do not mix, while in the run with feedback, the gas dispersal causes each group to expand, and the cluster expansion thus consists of the combined expansion of the groups. Most secondary star-forming sites around the main hub are also present in the non-feedback run, implying a primordial rather than triggered nature. We do find one example of a peripheral star-forming site that appears only in the feedback run, thus having a triggered origin. However, this appears to be the exception rather than the rule, although this may be an artefact of our simplified radiative transfer scheme.


2018 ◽  
Vol 2 (4) ◽  
pp. 74 ◽  
Author(s):  
Abinash Tripathy ◽  
Patryk Wąsik ◽  
Syama Sreedharan ◽  
Dipankar Nandi ◽  
Oier Bikondoa ◽  
...  

Functional ZnO nanostructured surfaces are important in a wide range of applications. Here we report the simple fabrication of ZnO surface structures at near room temperature with morphology resembling that of sea urchins, with densely packed, μm-long, tapered nanoneedles radiating from the urchin center. The ZnO urchin structures were successfully formed on several different substrates with high surface density and coverage, including silicon (Si), glass, polydimethylsiloxane (PDMS), and copper (Cu) sheets, as well as Si seeded with ZnO nanocrystals. Time-resolved SEM revealed growth kinetics of the ZnO nanostructures on Si, capturing the emergence of “infant” urchins at the early growth stage and subsequent progressive increases in the urchin nanoneedle length and density, whilst the spiky nanoneedle morphology was retained throughout the growth. ε-Zn(OH)2 orthorhombic crystals were also observed alongside the urchins. The crystal structures of the nanostructures at different growth times were confirmed by synchrotron X-ray diffraction measurements. On seeded Si substrates, a two-stage growth mechanism was identified, with a primary growth step of vertically aligned ZnO nanoneedle arrays preceding the secondary growth of the urchins atop the nanoneedle array. The antibacterial, anti-reflective, and wetting functionality of the ZnO urchins—with spiky nanoneedles and at high surface density—on Si substrates was demonstrated. First, bacteria colonization was found to be suppressed on the surface after 24 h incubation in gram-negative Escherichia coli (E. coli) culture, in contrast to control substrates (bare Si and Si sputtered with a 20 nm ZnO thin film). Secondly, the ZnO urchin surface, exhibiting superhydrophilic property with a water contact angle ~ 0°, could be rendered superhydrophobic with a simple silanization step, characterized by an apparent water contact angle θ of 159° ± 1.4° and contact angle hysteresis ∆θ < 7°. The dynamic superhydrophobicity of the surface was demonstrated by the bouncing-off of a falling 10 μL water droplet, with a contact time of 15.3 milliseconds (ms), captured using a high-speed camera. Thirdly, it was shown that the presence of dense spiky ZnO nanoneedles and urchins on the seeded Si substrate exhibited a reflectance R < 1% over the wavelength range λ = 200–800 nm. The ZnO urchins with a unique morphology fabricated via a simple route at room temperature, and readily implementable on different substrates, may be further exploited for multifunctional surfaces and product formulations.


2019 ◽  
Vol 488 (4) ◽  
pp. 5580-5593 ◽  
Author(s):  
Viraj Pandya ◽  
Joel Primack ◽  
Peter Behroozi ◽  
Avishai Dekel ◽  
Haowen Zhang ◽  
...  

ABSTRACT Hubble Space Telescope observations show that low-mass ($M_*=10^9\!-\!10^{10}\, \mathrm{M}_{\odot }$) galaxies at high redshift (z = 1.0–2.5) tend to be elongated (prolate) rather than disky (oblate) or spheroidal. This is explained in zoom-in cosmological hydrodynamical simulations by the fact that these galaxies are forming in cosmic web filaments where accretion happens preferentially along the direction of elongation. We ask whether the elongated morphology of these galaxies allows them to be used as effective tracers of cosmic web filaments at high redshift via their intrinsic alignments. Using mock light cones and spectroscopically confirmed galaxy pairs from the Cosmic Assembly Near-infared Deep Extragalactic Legacy Survey (CANDELS), we test two types of alignments: (1) between the galaxy major axis and the direction to nearby galaxies of any mass and (2) between the major axes of nearby pairs of low-mass, likely prolate, galaxies. The mock light cones predict strong signals in 3D real space, 3D redshift space, and 2D projected redshift space for both types of alignments (assuming prolate galaxy orientations are the same as those of their host prolate haloes), but we do not detect significant alignment signals in CANDELS observations. However, we show that spectroscopic redshifts have been obtained for only a small fraction of highly elongated galaxies, and accounting for spectroscopic incompleteness and redshift errors significantly degrades the 2D mock signal. This may partly explain the alignment discrepancy and highlights one of several avenues for future work.


2020 ◽  
Vol 637 ◽  
pp. A11 ◽  
Author(s):  
Thomas Baumann ◽  
Bertram Bitsch

Low-mass planets that are in the process of growing larger within protoplanetary disks exchange torques with the disk and change their semi-major axis accordingly. This process is called type I migration and is strongly dependent on the underlying disk structure. As a result, there are many uncertainties about planetary migration in general. In a number of simulations, the current type I migration rates lead to planets reaching the inner edge of the disk within the disk lifetime. A new kind of torque exchange between planet and disk, the thermal torque, aims to slow down inward migration via the heating torque. The heating torque may even cause planets to migrate outwards, if the planetary luminosity is large enough. Here, we study the influence on planetary migration of the thermal torque on top of previous type I models. We find that the formula of Paardekooper et al. (2011, MNRAS, 410, 293) allows for more outward migration than that of Jiménez & Masset (2017, MNRAS, 471, 4917) in most configurations, but we also find that planets evolve to very similar mass and final orbital radius using both formulae in a single planet-formation scenario, including pebble and gas accretion. Adding the thermal torque can introduce new, but small, regions of outwards migration if the accretion rates onto the planet correspond to typical solid accretion rates following the pebble accretion scenario. If the accretion rates onto the planets become very large, as could be the case in environments with large pebble fluxes (e.g., high-metallicity environments), the thermal torque can allow more efficient outward migration. However, even then, the changes for the final mass and orbital positions in our planet formation scenario are quite small. This implies that for single planet evolution scenarios, the influence of the heating torque is probably negligible.


2019 ◽  
Vol 629 ◽  
pp. A106 ◽  
Author(s):  
Sota Arakawa ◽  
Yuhito Shibaike

Jupiter has four large regular satellites called the Galilean satellites: Io, Europa, Ganymede, and Callisto. The inner three of the Galilean satellites orbit in a 4:2:1 mean motion resonance; therefore their orbital configuration may originate from the stopping of the migration of Io near the bump in the surface density distribution and following resonant trapping of Europa and Ganymede. The formation mechanism of the bump near the orbit of the innermost satellite, Io, is not yet understood, however. Here, we show that photophoresis in the circumjovian disk could be the cause of the bump using analytic calculations of steady-state accretion disks. We propose that photophoresis in the circumjovian disk could stop the inward migration of dust particles near the orbit of Io. The resulting dust-depleted inner region would have a higher ionization fraction, and thus admit increased magnetorotational-instability-driven accretion stress in comparison to the outer region. The increase of the accretion stress at the photophoretic dust barrier would form a bump in the surface density distribution, halting the migration of Io.


2020 ◽  
Vol 492 (3) ◽  
pp. 3440-3458 ◽  
Author(s):  
E Sanchis ◽  
G Picogna ◽  
B Ercolano ◽  
L Testi ◽  
G Rosotti

ABSTRACT We predict magnitudes for young planets embedded in transition discs, still affected by extinction due to material in the disc. We focus on Jupiter-sized planets at a late stage of their formation, when the planet has carved a deep gap in the gas and dust distributions and the disc starts to being transparent to the planet flux in the infrared (IR). Column densities are estimated by means of three-dimensional hydrodynamical models, performed for several planet masses. Expected magnitudes are obtained by using typical extinction properties of the disc material and evolutionary models of giant planets. For the simulated cases located at 5.2 au in a disc with a local unperturbed surface density of 127 $\mathrm{g} \, \mathrm{cm}^{-2}$, a 1MJ planet is highly extinct in the J, H, and Kbands, with predicted absolute magnitudes ≥ 50 mag. In the L and Mbands, extinction decreases, with planet magnitudes between 25 and 35 mag. In the Nband, due to the silicate feature on the dust opacities, the expected magnitude increases to ∼40 mag. For a 2MJ planet, the magnitudes in the J, H, and Kbands are above 22 mag, while for the L, M, and Nbands, the planet magnitudes are between 15 and 20 mag. For the 5MJ planet, extinction does not play a role in any IR band, due to its ability to open deep gaps. Contrast curves are derived for the transition discs in CQ Tau, PDS 70, HL Tau, TW Hya, and HD 163296. Planet mass upper limits are estimated for the known gaps in the last two systems.


Sign in / Sign up

Export Citation Format

Share Document