Earth and Planetary Astrophysics (astro-ph.EP)

Abstract: Disc accretion rate onto low mass protostar FU Ori suddenly increased hundreds of times 85 years ago and remains elevated to this day. We show that the sum of historic and recent observations challenges existing FU Ori models. We build a theory of a new process, Extreme Evaporation (EE) of young gas giant planets in discs with midplane temperatures exceeding 30, 000 K. Such temperatures are reached in the inner 0.1 AU during thermal instability bursts. In our 1D time-dependent code the disc and an embedded planet interact through gravity, heat, and mass exchange. We use disc viscosity constrained by simulations and observations of dwarf novae instabilities, and we constrain planet properties with a stellar evolution code. We show that dusty gas giants born in the outer self-gravitating disc reach the innermost disc in a $\sim$ 10,000 years with radius of $\sim 10 R_J$. We show that their EE rates are $\sim 10^{-5}$ Msun/yr; if this exceeds the background disc accretion activity then the system enters a planet-sourced mode. Like a stellar secondary in mass-transferring binaries, the planet becomes the dominant source of matter for the star, albeit for $\sim$ O(100) years. We find that a $\sim$ 6 Jupiter mass planet evaporating in a disc fed at a time-averaged rate of $\sim 10^{-6}$ Msun/yr appears to explain all that we currently know about FU Ori accretion outburst. More massive planets and/or planets in older less massive discs do not experience EE process. Future FUOR modelling may constrain planet internal structure and evolution of the earliest discs.

Abstract: Exomoons are expected to orbit gas giant exoplanets just as moons orbit solar system planets. Tidal heating is present in solar system satellites and it can heat up their interior depending on their orbital and interior properties. We aim to identify a Tidally Heated Exomoon's (THEM) orbital parameter space that would make it observable in infrared wavelengths with MIRI/JWST around $\epsilon$ Eridani b. We study the possible constraints on orbital eccentricity and interior properties that a successful THEM detection in infrared wavelengths can bring. We also investigate what exomoon properties need to be independently known in order to place these constraints. We use a coupled thermal-tidal model to find stable equilibrium points between the tidally produced heat and heat transported within a moon. For the latter, we consider a spherical and radially symmetric satellite with heat being transported via magma advection in a sub-layer of melt (asthenosphere) and convection in the lower mantle. We incorporate uncertainties in the interior and tidal model parameters to assess the fraction of simulated moons that would be observable with MIRI. We find that a $2 R_{Io}$ THEM orbiting $\epsilon$ Eridani b with an eccentricity of 0.02, would need to have a semi-major axis of 4 planetary Roche-radii for 100% of the simulations to produce an observable moon. These values are comparable with the orbital properties of gas giant solar system satellites. We place similar constraints for eccentricities up to 0.1. We conclude that if the semi-major axis and radius of the moon are known (eg. with exomoon transits), tidal dissipation can constrain the orbital eccentricity and interior properties of the satellite, such as the presence of melt and the thickness of the melt containing sub-layer.

Abstract: Cometary dust particles are best preserved remnants of the matter present at the onset of the formation of the Solar System. Space missions, telescopic observations and laboratory analyses advanced the knowledge on the properties of cometary dust. Cometary samples were returned from comet 81P/Wild2 by the Stardust mission. The chondritic (porous) anhydrous interplanetary dust particles and chondritic porous micrometeorites, and the ultracarbonaceous Antarctic micrometeorites (UCAMMs) also show strong evidence for a cometary origin. The composition of cometary dust is generally chondritic, but with high C and N compared with CI. The cometary organic matter is mixed with minor amounts of crystalline and amorphous minerals. The most abundant crystalline minerals are ferromagnesian silicates, refractory minerals and low Ni Fe sulfides are also present. The presence of carbonates in cometary dust is still debated, but a phyllosilicate-like phase was observed in a UCAMM. GEMS phases are usually abundant. Some of the organic matter present in cometary dust particle resembles the insoluble organic matter present in primitive meteorites, but amorphous carbon and exotic (e.g. N-rich) organic phases are also present. The H isotopic composition of the organic matter traces a formation at very low temperatures, in the protosolar cloud or in the outer regions of the protoplanetary disk. The presolar dust concentration in cometary dust can reach about 1%, which is the most elevated value observed in extraterrestrial samples. The differential size distribution of cometary dust in comet trails is well represented by a power-law distribution with a mean power index N typically ranging from -3 to -4. Polarimetric and light scattering studies suggest mixtures of porous agglomerates of sub-micrometer minerals with organic matter. Cometary dust particles have low tensile strength, and low density.

Abstract: The photospheric unsigned magnetic flux has been shown to be highly correlated with radial velocity (RV) variations caused by solar surface activity. This activity indicator is therefore a prime candidate to unlock the potential of RV surveys to discover Earth twins orbiting Sun-like stars. We show for the first time how a precise proxy of the unsigned magnetic flux ($\Delta\alpha B^2$) can be obtained from Sun-as-a-star intensity spectra by harnessing the magnetic information contained in over 4000 absorption lines in the wavelength range from 380 to 690 nm. This novel activity proxy can thus be obtained from the same spectra from which RVs are routinely extracted. We derived $\Delta\alpha B^2$ from 500 randomly selected spectra from the HARPS-N public solar data set, which spans from 2015 to 2018. We compared our estimates with the unsigned magnetic flux values from the Solar Dynamics Observatory (SDO) finding excellent agreement (median absolute deviation: 4.9 per cent). The extracted indicator $\Delta\alpha B^2$ correlates with SDO's unsigned magnetic flux estimates on the solar rotational timescale (Pearson correlation coefficient 0.67) and on the three-year timescale of our data set (correlation coefficient 0.91). We find correlations of $\Delta\alpha B^2$ with the HARPS-N solar RV variations of 0.49 on the rotational timescale and 0.78 on the three-year timescale. The Pearson correlation of $\Delta\alpha B^2$ with the RVs is found to be greater than the correlation of the classical activity indicators with the RVs. For solar-type stars, $\Delta\alpha B^2$ therefore represents the best simultaneous activity proxy known to date.

Abstract: Planetesimals are born fragile and are subject to destruction by wind erosion as they move through the gas of a protoplanetary disk. In microgravity experiments, we determined the shear stress necessary for erosion of a surface consisting of 1 mm dust pebbles down to 1 Pa ambient pressure. This is directly applicable to protoplanetary disks. Even pebble pile planetesimals with low eccentricities of 0.1 cannot survive inside of 1 au in a minimum-mass solar nebula, and safe zones for planetesimals with higher eccentricities are located even farther out.

Abstract: Young protostellar discs are likely to be both self-gravitating, and to support grain growth to sizes where the particles decoupled from the gas. This combination could lead to short-wavelength fragmentation of the solid component in otherwise non-fragmenting gas discs, forming Earth-mass solid cores during the Class 0/I stages of Young Stellar Object evolution. We use three-dimensional smoothed particle hydrodynamics simulations of two-fluid discs, in the regime where the Stokes number of the particles St>1, to study how the formation of solid clumps depends on the disc-to-star mass ratio, the strength of gravitational instability, and the Stokes number. Gravitational instability of the simulated discs is sustained by local cooling. We find that the ability of the spiral structures to concentrate solids increases with the cooling time, and decreases with the Stokes number, while the relative dynamical temperature between gas and dust of the particles decreases with the cooling time and the disc-to-star mass ratio, and increases with the Stokes number. Dust collapse occurs in a subset of high disc mass simulations, yielding clumps whose mass is close to linear theory estimates, namely 1-10 Earth masses. Our results suggest that if planet formation occurs via this mechanism, the best conditions correspond to near the end of the self-gravitating phase, when the cooling time is long and the Stokes number close to unity.