How Sublimation Delays the Onset of Dusty Debris Disk Formation Around White Dwarf Stars
Although numerous white dwarf stars host dusty debris disks, the temperature distribution of these stars differs significantly from the white dwarf population as a whole. Dusty debris disks exist exclusively around white dwarfs cooler than 27,000 K. This is all the more enigmatic given that the formation processes of dusty debris disks should favor younger, hotter white dwarfs, which likely host more dynamically unstable planetary systems. Here we apply a sophisticated material sublimation model to white dwarf systems to show that these statistics are actually a natural result of the interplay of thermal and tidal forces, and show how they define the circumstellar regions where dusty debris disks can form. We demonstrate that these processes tend to prevent stability against both sublimative destruction and reaccretion into planetesimals for rocky materials until white dwarfs cool to below ~25,000-32,000 K, in agreement with the observed limit of ~27,000 K. For pure water ice, this critical temperature is less than 2,700 K (requiring a cooling age older the universe); this precludes pure water ice-rich debris disks forming through the accepted two-step mechanism. The critical temperature is size-dependent; more massive white dwarfs could potentially host dusty debris disks at warmer temperatures.. Our model suggests that the location of the disks within the PG 0010+280, GD 56, GD 362, and PG 1541+651 systems are consistent with a forsterite-dominated olivine composition. We also find that very cool white dwarfs may simultaneously host multiple, independently formed dusty debris disks, consistent with observations of the LSPM J0207+3331 system.
Thiols in the ISM: first detection of HC(O)SH and confirmation of CHSH
The chemical compounds carrying the thiol group (-SH) have been considered essential in recent prebiotic studies regarding the polymerization of amino acids. We have searched for this kind of compounds toward the Galactic Centre quiescent cloud G+0.693-0.027. We report the first detection in the interstellar space of the trans-isomer of monothioformic acid (t-HC(O)SH) with an abundance of ~ 1 × 10. Additionally, we provide a solid confirmation of the gauche isomer of ethyl mercaptan (g-CHSH) with an abundance of ~ 3 × 10, and we also detect methyl mercaptan (CHSH) with an abundance of ~ 5 × 10. Abundance ratios were calculated for the three SH-bearing species and their OH-analogues, revealing similar trends between alcohols and thiols with increasing complexity. Possible chemical routes for the interstellar synthesis of t-HC(O)SH, CHSH and CHSH are discussed, as well as the relevance of these compounds in the synthesis of prebiotic proteins in the primitive Earth.
Slowdown and Heating of Interstellar Neutral Helium by Elastic Collisions Beyond the Heliopause
Direct sampling of interstellar neutral (ISN) atoms close to the Sun enables studies of the very local interstellar medium (VLISM) around the heliosphere. The primary population of ISN helium atoms has, until now, been assumed to reflect the pristine VLISM conditions at the heliopause. Consequently, the atoms observed at 1 au by the Interstellar Boundary Explorer (IBEX) were used to determine the VLISM temperature and velocity relative to the Sun, without accounting for elastic collisions with other species outside the heliopause. Here, we evaluate the effect of these collisions on the primary ISN helium population. We follow trajectories of helium atoms and track their collisions with slowed plasma and interstellar hydrogen atoms ahead of the heliopause. Atoms typically collide a few times in the outer heliosheath, and only ~1.5% of the atoms are not scattered at all. We use calculated differential cross sections to randomly choose scattering angles in these collisions. We estimate that the resulting primary ISN helium atoms at the heliopause are slowed down by ~0.45 km s and heated by ~1100 K compared to the pristine VLISM. The resulting velocity distribution is asymmetric and shows an extended tail in the antisunward direction. Accounting for this change in the parameters derived from IBEX observations gives the Sun's relative speed of 25.85 km s and temperature of 6400 K in the pristine VLISM. Finally, this paper serves as a source of the differential cross sections for elastic collisions with helium atoms.
Cluster Analysis of Presolar Silicon Carbide Grains: Evaluation of Their Classification and Astrophysical Implications
Cluster analysis of presolar silicon carbide grains based on literature data for C/C, N/N, Si/Si, and Si/Si including or not inferred initial Al/Al data, reveals nine clusters agreeing with previously defined grain types but also highlighting new divisions. Mainstream grains reside in three clusters probably representing different parent star metallicities. One of these clusters has a compact core, with a narrow range of composition, pointing to an enhanced production of SiC grains in asymptotic giant branch (AGB) stars with a narrow range of masses and metallicities. The addition of Al/Al data highlights a cluster of mainstream grains, enriched in N and Al, which cannot be explained by current AGB models. We defined two AB grain clusters, one with N and Al excesses, and the other with N and smaller Al excesses, in agreement with recent studies. Their definition does not use the solar N isotopic ratio as a divider, and the contour of the Al-rich AB cluster identified in this study is in better agreement with core-collapse supernova models. We also found a cluster with a mixture of putative nova and AB grains, which may have formed in supernova or nova environments. X grains make up two clusters, having either strongly correlated Si isotopic ratios or deviating from the 2/3 slope line in the Si 3-isotope plot. Finally, most Y and Z grains are jointly clustered, suggesting that the previous use of C/C = 100 as a divider for Y grains was arbitrary. Our results show that cluster analysis is a powerful tool to interpret the data in light of stellar evolution and nucleosynthesis modeling and highlight the need of more multi-element isotopic data for better classification.
Explaining the Galilean Satellites' Density Gradient by Hydrodynamic Escape
The Galilean satellites exhibit a monotonic decrease in density (and increase in ice mass fraction) with distance from Jupiter (Pollack & Fanale 1982). Whether this is because of the background conditions when they formed (Lunine & Stevenson 1982; Canup & Ward 2002; Mosqueira & Estrada 2003a; Ronnet et al. 2017), the process of accretion itself (Dwyer et al. 2013), or later loss due to tidal heating (Canup & Ward 2009), has been in dispute for forty years. We find that a hitherto largely neglected process - vapor loss driven by accretional heating (Kuramoto & Matsui 1994) - can reproduce the observed density trend for accretion timescales ≳300 kyr, consistent with gas-starved satellite formation models (Canup & Ward 2002, 2006). In this model both Io and Europa develop an early surface liquid water ocean. Vapor escape from this ocean causes the water inventories of Io and Europa to be completely and mostly lost, respectively. Isotopic fractionation arising from vapor loss means that Europa will develop a higher D/H ratio compared with Ganymede and Callisto. We make predictions that can be tested with in situ measurements of D/H of potential Europa plumes (Roth et al. 2014) by the Europa Clipper spacecraft, or infrared spectroscopic determinations (Clark et al. 2019) of D/H at all three bodies.
NuSTAR Observation of a Minuscule Microflare in a Solar Active Region
We present X-ray imaging spectroscopy of one of the weakest active region (AR) microflares ever studied. The microflare occurred at ~11:04 UT on 2018 September 9 and we studied it using the Nuclear Spectroscopic Telescope ARray (NuSTAR) and the Solar Dynamic Observatory's Atmospheric Imaging Assembly (SDO/AIA). The microflare is observed clearly in 2.5-7 keV with NuSTAR and in Fe XVIII emission derived from the hotter component of the 94 Å SDO/AIA channel. We estimate the event to be three orders of magnitude lower than a GOES A class microflare with an energy of 1.1 × 10 erg. It reaches temperatures of 6.7 MK with an emission measure of 8.0 × 10 cm. Non-thermal emission is not detected but we instead determine upper limits to such emission. We present the lowest thermal energy estimate for an AR microflare in literature, which is at the lower limits of what is still considered an X-ray microflare.
Accelerated Electrons Observed Down to <7 keV in a Solar Microflare
We report the detection of emission from a nonthermal electron distribution in a small solar microflare (GOES class A5.7) observed by the , with supporting observation by the (). The flaring plasma is well accounted for by a thick-target model of accelerated electrons collisionally thermalizing within the loop, akin to the "coronal thick-target" behavior occasionally observed in larger flares. This is the first positive detection of nonthermal hard X-rays from the Sun using a direct imager (as opposed to indirectly imaging instruments). The accelerated electron distribution has a spectral index of 6.3 ± 0.7, extends down to at least 6.5 keV, and deposits energy at a rate of ~2 × 10 erg s, heating the flare loop to at least 10 MK. The existence of dominant nonthermal emission in X-rays down to <5 keV means that emission is almost entirely nonthermal, contrary to what is usually assumed in spectroscopy. The ratio of nonthermal to thermal energies is similar to that of large flares, in contrast to what has been found in previous studies of small flares. We suggest that a coronal thick target may be a common property of many small microflares based on the average electron energy and collisional mean free path. Future observations of this kind will enable understanding of how flare particle acceleration changes across energy scales, and will aid the push toward the observational regime of nanoflares, which are a possible source of significant coronal heating.
Laboratory Observations of Ultra-Low Frequency Analogue Waves Driven by the Right-Hand Resonant Ion Beam Instability
The Right-Hand Resonant Instability (RHI) is one of several electromagnetic ion/ion beam instabilities responsible for the formation of parallel magnetized collisionless shocks and the generation of ultra-low frequency (ULF) waves in their foreshocks. This instability has been observed for the first time under foreshock-relevant conditions in the laboratory through the repeatable interaction of a preformed magnetized background plasma and a super-Alfvénic laser-produced plasma. This platform has enabled unprecedented volumetric measurements of waves generated by the RHI, revealing filamentary current structures in the transverse plane. These measurements are made in the plasma rest frame with both high spatial and temporal resolution, providing a perspective that is complementary to spacecraft observations. Direct comparison of data from both the experiment and the Wind spacecraft to 2D hybrid simulations demonstrates that the waves produced are analogous to the ULF waves observed upstream of the terrestrial bow shock.
Comet C/2017 S3 (PanSTARRS): Outbursts and Disintegration
The Solar Wind ANisotropies (SWAN) all-sky hydrogen Lyman-alpha camera on the SOlar and Heliospheric Observer (SOHO) satellite observed the hydrogen coma of comet C/2017 S3 (PanSTARRS) for the last month of its activity from 2018 July 4 to August 4 and what appears to have been its final disintegration just 11 days before its perihelion on August 15. The hydrogen coma indicated water production had a small outburst on July 8 at a heliocentric distance of 1.1AU and then a much larger one on July 20 at 0.8 AU. Over the following two weeks the water production dropped by more than a factor of ten after which it was no longer detectable. The behavior is reminiscent of comet C/1999 S4 (LINEAR) in 2000, which had a few small outbursts on its inbound orbit and a major outburst at a heliocentric distance of about 0.8 AU, which was close to its perihelion, followed by its complete disintegration that was documented by several sets of observations including SWAN. C/2017 S3 (PanSTARRS) however had a much larger water production rate than C/1999 S4 (LINEAR). Here we estimate the size of the nucleus of C/2017 S3 just before its final outburst and apparent disintegration was estimated using the total amount of water produced during its last weeks for a range of values of the refractory/ice ratio in the nucleus. We also determine the size distribution of the disintegrating particles as the comet faded.
Gas Accretion within the Dust Cavity in AB Aur
AB Aur is a Herbig Ae star hosting a well-known transitional disk. Because of its proximity and low inclination angle, it is an excellent object to study planet formation. Our goal is to investigate the chemistry and dynamics of the molecular gas component in the AB Aur disk, and its relation with the prominent horseshoe shape observed in continuum mm emission. We used the NOEMA interferometer to map with high angular resolution the J = 3-2 lines of HCO and HCN. By combining both, we can gain insight into the AB Aur disk structure. Chemical segregation is observed in the AB Aur disk: HCO shows intense emission toward the star position, at least one bright molecular bridge within the dust cavity, and ring-like emission at larger radii, while HCN is only detected in an annular ring that is coincident with the dust ring and presents an intense peak close to the dust trap. We use HCO to investigate the gas dynamics inside the cavity. The observed bright HCO bridge connects the compact central source with the outer dusty ring. This bridge can be interpreted as an accretion flow from the outer ring to the inner disk/jet system proving gas accretion through the cavity.
Cloud Coalescence: A Dynamical Instability Affecting Multiphase Environments
Mass and size distributions are the key characteristics of any astrophysical object, including the densest clumps comprising the cold phase of multiphase environments. In our recent papers, we showed how individual clouds of various sizes form and evolve in active galactic nuclei. In particular, we showed that large clouds undergo damped oscillations as a response to their formation process. Here we follow up this investigation, addressing how different size clouds interact. We find that smaller clouds become trapped in the advective flows generated by larger clouds. The explanation for this behavior leads to a rather remarkable conclusion: even in the absence of gravity, complexes of clouds are dynamically unstable. In an idealized environment (e.g., one free of turbulence and magnetic fields) a perfectly symmetric arrangement of static clouds will remain static, but any small spatial perturbation will lead to all clouds coalescing into a single, large cloud, given enough time. Using numerical simulations, we investigate the main factors that determine the rate of coalescence. In addition to the cloud separation distance, we find that the transient response of clouds to a disturbance is the primary factor. Turbulent motions in the flow can easily suppress this tendency for spatially well-separated clouds to coalesce, so it is as yet unclear if this phenomenon can occur in nature. Nevertheless, this Letter casts strong doubts on a recent hypothesis that large clouds are prone to fragmentation.
Multi-TeV Emission From the Vela Pulsar
Pulsed emission from the Vela pulsar at energies above 3 TeV has recently been detected by the H.E.S.S. II air-Cherenkov telescope. We present a model for the broad-band spectrum of Vela from infra-red (IR) to beyond 10 TeV. Recent simulations of the global pulsar magnetosphere have shown that most of the particle acceleration occurs in the equatorial current sheet outside the light cylinder and that the magnetic field structure is nearly force-free for younger pulsars. We adopt this picture to compute the radiation from both electron-positron pairs produced in polar cap cascades and from primary particles accelerated in the separatrix and current sheet. The synchrotron spectrum from pairs resonantly absorbing radio photons at relatively low altitude can account for the observed IR-optical emission. We set the parallel electric field in the current sheet to produce the GeV emission through curvature radiation, producing particles with energies of 30-60 TeV. These particles then produce Very-High-Energy emission up to around 30 TeV through inverse-Compton scattering of the IR-Optical emission. We present model spectra and light curves that can match the IR-Optical, GeV and make predictions for the multi-TeV emission.
TESS DISCOVERY OF A TRANSITING SUPER-EARTH IN THE MENSAE SYSTEM
We report the detection of a transiting planet around Men (HD 39091), using data from the (). The solar-type host star is unusually bright ( = 5.7) and was already known to host a Jovian planet on a highly eccentric, 5.7-year orbit. The newly discovered planet has a size of 2.04 ± 0.05 and an orbital period of 6.27 days. Radial-velocity data from the HARPS and AAT/UCLES archives also displays a 6.27-day periodicity, confirming the existence of the planet and leading to a mass determination of 4.82±0.85 . The star's proximity and brightness will facilitate further investigations, such as atmospheric spectroscopy, asteroseismology, the Rossiter-McLaughlin effect, astrometry, and direct imaging.
Variable Nature of Magnetically-Driven Ultra-Fast Outflows
Among a number of active galactic nuclei (AGNs) that drive ionized outflows in X-rays, a low-redshift ( = 0.184) quasar, PDS 456, is long known to exhibit one of the exemplary ultra-fast outflows (UFOs). However, the physical process of acceleration mechanism is yet to be definitively constrained. In this work, we model the variations of the Fe K UFO properties in PDS 456 over many epochs in X-ray observations in the context of magnetohydrodynamic (MHD) accretion-disk winds employed in our earlier studies of similar X-ray absorbers. We applied the model to the 2013/2014 spectra to determine the UFO's condition; namely, velocity, ionization parameter, column density and equivalent width (EW). Under some provisions on the dependence of X-ray luminosity on the accretion rate applicable to near-Eddington state, our photoionization calculations, coupled to a 2.5-dimensional MHD-driven wind model, can further reproduce the observed correlations of the UFO velocity and the anti-correlation of its EW with X-ray strength of PDS 456. This work demonstrates that UFOs, even without radiative pressure, can be driven as an extreme case purely by magnetic interaction while also producing the observed spectrum and correlations.
Discovery of Interstellar Isocyanogen (CNCN): Further Evidence that Dicyanopolyynes are Abundant in Space
It is thought that dicyanopolyynes could be potentially abundant interstellar molecules, although their lack of dipole moment makes it impossible to detect them through radioastronomical techniques. Recently, the simplest member of this chemical family, cyanogen (NCCN), was indirectly probed for the first time in interstellar space through the detection of its protonated form toward the dense clouds L483 and TMC-1. Here we present a second firm evidence of the presence of NCCN in interstellar space, namely the detection of the metastable and polar isomer isocyanogen (CNCN). This species has been identified in L483 and tentatively in TMC-1 by observing various rotational transitions in the λ 3 mm band with the IRAM 30m telescope. We derive beam-averaged column densities for CNCN of 1.6 × 10 cm in L483 and 9 × 10 cm in TMC-1, which imply fractional abundances relative to H in the range (5 - 9) × 10. While the presence of NCCN in interstellar clouds seems out of doubt owing to the detection of NCCNH and CNCN, putting tight constraints on its abundance is still hampered by the poor knowledge of the chemistry that links NCCN with NCCNH and especially with CNCN. We estimate that NCCN could be fairly abundant, in the range 10-10 relative to H, as other abundant nitriles like HCN and HCN.
Extremely Cr- and Ti-rich presolar oxide grains in a primitive meteorite: Formation in rare types of supernovae and implications for the astrophysical context of solar system birth
We report the identification of 19 presolar oxide grains from the Orgueil CI meteorite with substantial enrichments in Cr, with Cr/Cr ratios ranging from 1.2 to 56 times the solar value. The most enriched grains also exhibit enrichments at mass 50, most likely due in part to Ti, but close-to-normal or depleted Cr/Cr ratios. There is a strong inverse relationship between Cr enrichment and grain size; the most extreme grains are all <80 nm in diameter. Comparison of the isotopic data with predictions of nucleosynthesis calculations indicate that these grains most likely originated in either rare, high-density Type Ia supernovae (SNIa), or in electron-capture supernovae (ECSN) which may occur as the end stage of evolution for stars of mass 8-10 . This is the first evidence for preserved presolar grains from either type of supernova. An ECSN origin is attractive since these likely occur much more frequently than high-density SNIa, and their evolutionary timescales (~20 Myr) are comparable to those of molecular clouds. Self-pollution of the Sun's parent cloud from an ECSN may explain the heterogeneous distribution of n-rich isotopic anomalies in planetary materials, including a recently reported dichotomy in Mo isotopes in the solar system. The stellar origins of three grains with solar Cr/Cr, but anomalies in Cr or Cr, as well as of a grain enriched in Fe, are unclear.
Thermal History of CB Chondrules and Cooling Rate Distributions of Ejecta Plumes
It has been proposed that some meteorites, CB and CH chondrites, contain material formed as a result of a protoplanetary collision during accretion. Their melt droplets (chondrules) and FeNi metal are proposed to have formed by evaporation and condensation in the resulting impact plume. We observe that the skeletal olivine (SO) chondrules in CB chondrites have a blebby texture and an enrichment in refractory elements not found in normal chondrules. Because the texture requires complete melting, their maximum liquidus temperature of 1928 K represents a minimum temperature for the putative plume. Dynamic crystallization experiments show that the SO texture can be created only by brief reheating episodes during crystallization, giving a partial dissolution of olivine. The ejecta plume formed in a smoothed particle hydrodynamics simulation served as the basis for 3D modeling with the adaptive mesh refinement code FLASH4.3. Tracer particles that move with the fluid cells are used to measure the in situ cooling rates. Their cooling rates are ~10,000 K hr briefly at peak temperature and, in the densest regions of the plume, ~100 K hr for 1400-1600 K. A small fraction of cells is seen to be heating at any one time, with heating spikes explained by the compression of parcels of gas in a heterogeneous patchy plume. These temperature fluctuations are comparable to those required in crystallization experiments. For the first time, we find an agreement between experiments and models that supports the plume model specifically for the formation of CB chondrules.
Discovery of the Ubiquitous Cation NS in Space Confirmed by Laboratory Spectroscopy
We report the detection in space of a new molecular species which has been characterized spectroscopically and fully identified from astrophysical data. The observations were carried out with the 30m IRAM telescopea. The molecule is ubiquitous as its =2→1 transition has been found in cold molecular clouds, prestellar cores, and shocks. However, it is not found in the hot cores of Orion-KL and in the carbon-rich evolved star IRC+10216. Three rotational transitions in perfect harmonic relation ' = 2/3/5 have been identified in the prestellar core B1b. The molecule has a Σ electronic ground state and its =2→1 transition presents the hyperfine structure characteristic of a molecule containing a nucleus with spin 1. A careful analysis of possible carriers shows that the best candidate is NS. The derived rotational constant agrees within 0.3-0.7% with ab initio calculations. NS was also produced in the laboratory to unambiguously validate the astrophysical assignment. The observed rotational frequencies and determined molecular constants confirm the discovery of the nitrogen sulfide cation in space. The chemistry of NS and related nitrogen-bearing species has been analyzed by means of a time-dependent gas phase model. The model reproduces well the observed NS/NS abundance ratio, in the range 30-50, and indicates that NS is formed by reactions of the neutral atoms N and S with the cations SH and NH, respectively.
First Detection of Interstellar SH
We present the first detection of gas phase SH in the Horsehead, a moderately UV-irradiated nebula. This confirms the presence of doubly sulfuretted species in the interstellar medium and opens a new challenge for sulfur chemistry. The observed SH abundance is ~5×10, only a factor 4-6 lower than that of the widespread HS molecule. HS and SH are efficiently formed on the UV-irradiated icy grain mantles. We performed ice irradiation experiments to determine the HS and SH photodesorption yields. The obtained values are ~1.2×10 and <1×10 molecules per incident photon for HS and SH, respectively. Our upper limit to the SH photodesorption yield suggests that photo-desorption is not a competitive mechanism to release the SH molecules to the gas phase. Other desorption mechanisms such as chemical desorption, cosmic-ray desorption and grain shattering can increase the gaseous SH abundance to some extent. Alternatively, SH can be formed via gas phase reactions involving gaseous HS and the abundant ions S and SH. The detection of SH in this nebula could be therefore the result of the coexistence of an active grain surface chemistry and gaseous photo-chemistry.
High-resolution SOFIA/EXES Spectroscopy of SO Gas in the Massive Young Stellar Object MonR2 IRS3: Implications for the Sulfur Budget
Sulfur has been observed to be severely depleted in dense clouds leading to uncertainty in the molecules that contain it and the chemistry behind their evolution. Here, we aim to shed light on the sulfur chemistry in young stellar objects (YSOs) by using high-resolution infrared spectroscopy of absorption by the rovibrational band of SO obtained with the Echelon-Cross-Echelle Spectrograph on the Stratospheric Observatory for Infrared Astronomy. Using local thermodynamic equilibrium models we derive physical parameters for the SO gas in the massive YSO MonR2 IRS3. This yields a SO/H abundance lower limit of 5.6 ± 0.5 × 10, or >4% of the cosmic sulfur budget, and an intrinsic line width (Doppler parameter) of < 3.20 km s. The small line widths and high temperature ( = 234 ± 15 K) locate the gas in a relatively quiescent region near the YSO, presumably in the hot core where ices have evaporated. This sublimation unlocks a volatile sulfur reservoir (e.g., sulfur allotropes as detected abundantly in comet 67P/Churyumov-Gerasimenko), which is followed by SO formation by warm, dense gas-phase chemistry. The narrowness of the lines makes formation of SO from sulfur sputtered off grains in shocks less likely toward MonR2 IRS3.
Particle-in-cell Simulations of the Whistler Heat-flux Instability in Solar Wind Conditions
In collision-poor plasmas from space, e.g., solar wind or stellar outflows, the heat flux carried by the strahl or beaming electrons is expected to be regulated by the self-generated instabilities. Recently, simultaneous field and particle observations have indeed revealed enhanced whistler-like fluctuations in the presence of counter-beaming populations of electrons, connecting these fluctuations to the whistler heat-flux instability (WHFI). This instability is predicted only for limited conditions of electron beam-plasmas, and has not yet been captured in numerical simulations. In this Letter we report the first simulations of WHFI in particle-in-cell setups, realistic for the solar wind conditions, and without temperature gradients or anisotropies to trigger the instability in the initiation phase. The velocity distributions have a complex reaction to the enhanced whistler fluctuations conditioning the instability saturation by a decrease of the relative drifts combined with induced (effective) temperature anisotropies (heating the core electrons and pitch-angle and energy scattering the strahl). These results are in good agreement with a recent quasilinear approach, and support therefore a largely accepted belief that WHFI saturates at moderate amplitudes. In the anti-sunward direction the strahl becomes skewed with a pitch-angle distribution decreasing in width as electron energy increases, which seems to be characteristic of self-generated whistlers and not to small-scale turbulence.