SOFIA upGREAT/FIFI-LS Emission-line Observations of Betelgeuse during the Great Dimming of 2019/2020

We report NASA-DLR SOFIA upGREAT circumstellar [O i] 63.2 μm and [C ii] 157.7 μm emission profiles and FIFI-LS [O i] 63.2 μm, [O i] 145.5 μm, and [C ii] 157.7 μm fluxes obtained shortly after Betelgeuse’s 2019/2020 Great Dimming event. Haas et al. noted a potential correlation between the [O i] 63.2 μm flux and V magnitude based on three Kuiper Airborne Observatory observations made with the CGS and FIFI instruments. The FIFI observation was obtained when V ≃ 0.88 and revealed a 3σ non-detection at a quarter of the previous CGS flux measurement made when V ≃ 0.35. A potential explanation could be a change in dust-gas drag heating by circumstellar silicates caused by variations in the photospheric radiation field. SOFIA observations provide a unique test of this correlation because the V-band brightness went to its lowest value on record, V ≃ 1.61, with the SOFIA observations being made when V FIFI−LS ≃ 1.51 and V upGREAT ≃ 1.36. The upGREAT spectra show a [O i] 63.2 μm flux larger than previous space observatory measurements obtained when V ≃ 0.58. The profile is consistent with formation in the slower, more turbulent inner S1 outflow, while the [C ii] 157.7 μm profile is consistent with formation farther out in the faster S2 outflow. Modeling of dust-gas drag heating, combined with 25 yr of Wing three-filter and V photometry, reveals that it is unlikely that the S1 circumstellar envelope and [O i] 63.2 μm fluxes are dominated by the dust-gas drag heating and that another heating source is also active. The [O i] 63.2 μm profile is hard to reconcile with existing outflow velocity models.

Modeling of Melt Flow and Heat Transfer in Stationary Gas Tungsten Arc Welding with Vertical and Tilted Torches

A 3D numerical simulation was conducted to study the transient development of temperature distribution in stationary gas tungsten arc welding with filler wire. Heat transfer to the filler wire and the workpiece was investigated with vertical (90°) and titled (70°) torches. Heat flux, current flux, and gas drag force were calculated from the steady-state simulation of the arc. The temperature in the filler wire was determined at three different time intervals: 0.12 s, 0.24 s, and 0.36 s. The filler wire was assumed not to deform during this short time, and was therefore simulated as solid. The temperature in the workpiece was calculated at the same intervals using heat flux, current flux, gas drag force, Marangoni convection, and buoyancy. It should be noted that heat transfer to the filler wire was faster with the titled torch compared to the vertical torch. Heat flux to the workpiece was asymmetrical with both the vertical and tilted torches when the filler wire was fully inserted into the arc. It was found that the overall trends of temperature contours for both the arc and the workpiece were in good agreement. It was also observed that more heat was transferred to the filler wire with the 70° torch compared with the 90° torch. The melted volume of the filler wire (volume above 1750 °K) was 12 mm3 with the 70° torch, compared to 9.2 mm3 with the 90° torch.

Circularization of tidal debris around white dwarfs: implications for gas production and dust variability

White dwarf (WD) pollution is thought to arise from the tidal disruption of planetary bodies. The initial fragment stream is extremely eccentric, while observational evidence suggest that discs are circular or nearly so. Here we propose a novel mechanism to bridge this gap and show that the fragments can rapidly circularise through dust or gas drag when they interact with a pre-existing compact disc. We assume that the tidal stream mainly consists of small cohesive fragments in the size range 10-1000 m, capable of resisting the WD tidal forces, whereas the compact discs span a wide mass range. We provide an analytical model, accompanied by N-body simulations, and find a large parameter space in fragment sizes and orbital separation that leads to full circularization. Partial circularization is possible for compact discs that are several orders of magnitudes less massive. We show that dust-induced circularization inherently produces gas as tidal fragments collisionally vaporize the pre-existing dust along their path. We show that ongoing gas production has a higher probability to occur during the early stages of tidal disruption events, resulting from the fact that smaller fragments are the first to circularize. Intermittent gas production however becomes more likely as the tidal stream matures. This could explain why only a small subset of systems with dusty compact discs also have an observed gaseous component. Additionally, the interaction yields fragment erosion by collisional shattering, sputtering, sublimation and possibly ram-pressure. Material scattered by the collisions might form a thin dusty halo that evolves through PR drag, in compatibility with observed infrared variability.

Massive black hole binary inspiral and spin evolution in a cosmological framework

Massive black hole (MBH) binary inspiral time scales are uncertain, and their spins are even more poorly constrained. Spin misalignment introduces asymmetry in the gravitational radiation, which imparts a recoil kick to the merged MBH. Understanding how MBH binary spins evolve is crucial for determining their recoil velocities, their gravitational wave (GW) waveforms detectable with LISA, as well as their retention rate in galaxies. Here we introduce a sub-resolution model for gas- and GW-driven MBH binary spin evolution using accreting MBHs from the Illustris cosmological hydrodynamics simulations. We also model binary inspiral via dynamical friction, stellar scattering, viscous gas drag, and GW emission. Our model assumes that the circumbinary disk always removes angular momentum from the binary. It also assumes differential accretion, which causes greater alignment of the secondary MBH spin in unequal-mass mergers. We find that 47% of the MBHs in our population merge by z = 0. Of these, 19% have misaligned primaries and 10% have misaligned secondaries at the time of merger in our fiducial model with initial eccentricity of 0.6 and accretion rates from Illustris. The MBH misalignment fraction depends strongly on the accretion disc parameters, however. Reducing accretion rates by a factor of 100, in a thicker disc, yields 79% and 42% misalignment for primaries and secondaries, respectively. Even in the more conservative fiducial model, more than 12% of binaries experience recoils of > 500km/s, which could displace them at least temporarily from galactic nuclei. We additionally find that a significant number of systems experience strong precession.

Inner dusty regions of protoplanetary discs – II. Dust dynamics driven by radiation pressure and disc winds

ABSTRACT We explore dust flow in the hottest parts of protoplanetary discs using the forces of gravity, gas drag, and radiation pressure. Our main focus is on the optically thin regions of dusty disc, where the dust is exposed to the most extreme heating conditions and dynamical perturbations: the surface of optically thick disc and the inner dust sublimation zone. We utilize results from two numerically strenuous fields of research. The first is the quasi-stationary solutions on gas velocity and density distributions from mangetohydrodynamical (MHD) simulations of accretion discs. This is critical for implementing a more realistic gas drag impact on dust movements. The second is the optical depth structure from a high-resolution dust radiation transfer. This step is critical for a better understanding of dust distribution within the disc. We describe a numerical method that incorporates these solutions into the dust dynamics equations. We use this to integrate dust trajectories under different disc wind models and show how grains end up trapped in flows that range from simple accretion on to the star to outflows into outer disc regions. We demonstrate how the radiation pressure force plays one of the key roles in this process and cannot be ignored. It erodes the dusty disc surface, reduces its height, resists dust accretion on to the star, and helps the disc wind in pushing grains outwards. The changes in grain size and porosity significantly affect the results, with smaller and porous grains being influenced more strongly by the disc wind and radiation pressure.

Physical models of streaming instabilities in protoplanetary discs

ABSTRACT We develop simple, physically motivated models for drag-induced dust–gas streaming instabilities, which are thought to be crucial for clumping grains to form planetesimals in protoplanetary discs. The models explain, based on the physics of gaseous epicyclic motion and dust–gas drag forces, the most important features of the streaming instability and its simple generalization, the disc settling instability. Some of the key properties explained by our models include the sudden change in the growth rate of the streaming instability when the dust-to-gas mass ratio surpasses one, the slow growth rate of the streaming instability compared to the settling instability for smaller grains, and the main physical processes underlying the growth of the most unstable modes in different regimes. As well as providing helpful simplified pictures for understanding the operation of an interesting and fundamental astrophysical fluid instability, our models may prove useful for analysing simulations and developing non-linear theories of planetesimal growth in discs.

Migrating low-mass planets in inviscid dusty protoplanetary discs

ABSTRACT Disc-driven planet migration is integral to the formation of planetary systems. In standard, gas-dominated protoplanetary discs, low-mass planets or planetary cores undergo rapid inwards migration and are lost to the central star. However, several recent studies indicate that the solid component in protoplanetary discs can have a significant dynamical effect on disc–planet interaction, especially when the solid-to-gas mass ratio approaches unity or larger and the dust-on-gas drag forces become significant. As there are several ways to raise the solid abundance in protoplanetary discs, for example through disc winds and dust trapping in pressure bumps, it is important to understand how planets migrate through a dusty environment. To this end, we study planet migration in dust-rich discs via a systematic set of high-resolution, two-dimensional numerical simulations. We show that the inwards migration of low-mass planets can be slowed down by dusty dynamical corotation torques. We also identify a new regime of stochastic migration applicable to discs with dust-to-gas mass ratios of ≳0.3 and particle Stokes numbers ≳0.03. In these cases, disc–planet interaction leads to the continuous development of small-scale, intense dust vortices that scatter the planet, which can potentially halt or even reverse the inwards planet migration. We briefly discuss the observational implications of our results and highlight directions for future work.

The aeolian-erosion barrier for the growth of metre-size objects in protoplanetary discs

ABSTRACT Aeolian erosion is a destructive process that can erode small-size planetary objects through their interaction with a gaseous environment. Aeolian erosion operates in a wide range of environments and under various conditions. Aeolian erosion has been extensively explored in the context of geophysics in terrestrial planets. Here we show that aeolian erosion of cobbles, boulders, and small planetesimals in protoplanetary discs can constitute a significant barrier for the early stages of planet formation. We use analytic calculations to show that under the conditions prevailing in protoplanetary discs small bodies ($10\!-\!10^4 \, \rm {m}$) are highly susceptible to gas-drag aeolian erosion. At this size-range aeolian erosion can efficiently erode the planetesimals down to tens-cm size and quench any further growth of such small bodies. It thereby raises potential difficulties for channels suggested to alleviate the metre-size barrier. Nevertheless, the population of ∼decimetre-size cobbles resulting from aeolian erosion might boost the growth of larger (>km size) planetesimals and planetary embryos through increasing the efficiency of pebble-accretion, once/if such large planetesimals and planetary embryos exist in the disc.

The solar nebula origin of (486958) Arrokoth, a primordial contact binary in the Kuiper Belt

The New Horizons spacecraft’s encounter with the cold classical Kuiper Belt object (486958) Arrokoth (provisional designation 2014 MU69) revealed a contact-binary planetesimal. We investigated how Arrokoth formed and found that it is the product of a gentle, low-speed merger in the early Solar System. Its two lenticular lobes suggest low-velocity accumulation of numerous smaller planetesimals within a gravitationally collapsing cloud of solid particles. The geometric alignment of the lobes indicates that they were a co-orbiting binary that experienced angular momentum loss and subsequent merger, possibly because of dynamical friction and collisions within the cloud or later gas drag. Arrokoth’s contact-binary shape was preserved by the benign dynamical and collisional environment of the cold classical Kuiper Belt and therefore informs the accretion processes that operated in the early Solar System.

The origin of the high metallicity of close-in giant exoplanets

Context. Recent studies suggest that in comparison to their host star, many giant exoplanets are highly enriched with heavy elements and can contain several tens of Earth masses of heavy elements or more. Such enrichment is considered to have been delivered by the accretion of planetesimals in late formation stages. Previous dynamical simulations, however, have shown that planets cannot accrete such high masses of heavy elements through “in situ” planetesimal accretion. Aims. We investigate whether a giant planet migrating inward can capture planetesimals efficiently enough to significantly increase its metallicity. Methods. We performed orbital integrations of a migrating giant planet and planetesimals in a protoplanetary gas disc to infer the planetesimal mass that is accreted by the planet. Results. We find that the two shepherding processes of mean motion resonance trapping and aerodynamic gas drag inhibit the planetesimal capture of a migrating planet. However, the amplified libration allows the highly-excited planetesimals in the resonances to escape from the resonance trap and to be accreted by the planet. Consequently, we show that a migrating giant planet captures planetesimals with total mass of several tens of Earth masses if the planet forms at a few tens of AU in a relatively massive disc. We also find that planetesimal capture occurs efficiently in a limited range of semi-major axis and that the total captured planetesimal mass increases with increasing migration distances. Our results have important implications for understanding the relation between giant planet metallicity and mass, as we suggest that it reflects the formation location of the planet – or more precisely, the location where runaway gas accretion occurred. We also suggest the observed metal-rich close-in Jupiters migrated to their present locations from afar, where they had initially formed.