Showing papers in "arXiv: Solar and Stellar Astrophysics in 2019"

The Revised TESS Input Catalog and Candidate Target List.

Abstract: We describe the catalogs assembled and the algorithms used to populate the revised TESS Input Catalog (TIC), based on the incorporation of the Gaia second data release. We also describe a revised ranking system for prioritizing stars for 2-minute cadence observations, and assemble a revised Candidate Target List (CTL) using that ranking. The TIC is available on the Mikulski Archive for Space Telescopes (MAST) server, and an enhanced CTL is available through the Filtergraph data visualization portal system at the URL this http URL.

Modules for Experiments in Stellar Astrophysics (MESA): Pulsating Variable Stars, Rotation, Convective Boundaries, and Energy Conservation

Abstract: We update the capabilities of the open-knowledge software instrument Modules for Experiments in Stellar Astrophysics (MESA). RSP is a new functionality in MESAstar that models the non-linear radial stellar pulsations that characterize RR Lyrae, Cepheids, and other classes of variable stars. We significantly enhance numerical energy conservation capabilities, including during mass changes. For example, this enables calculations through the He flash that conserve energy to better than 0.001 %. To improve the modeling of rotating stars in MESA, we introduce a new approach to modifying the pressure and temperature equations of stellar structure, and a formulation of the projection effects of gravity darkening. A new scheme for tracking convective boundaries yields reliable values of the convective-core mass, and allows the natural emergence of adiabatic semiconvection regions during both core hydrogen- and helium-burning phases. We quantify the parallel performance of MESA on current generation multicore architectures and demonstrate improvements in the computational efficiency of radiative levitation. We report updates to the equation of state and nuclear reaction physics modules. We briefly discuss the current treatment of fallback in core-collapse supernova models and the thermodynamic evolution of supernova explosions. We close by discussing the new MESA Testhub software infrastructure to enhance source-code development.

CHIANTI - an atomic database for emission lines - Paper XV: Version 9, improvements for the X-ray satellite lines

Abstract: CHIANTI contains a large quantity of atomic data for the analysis of astrophysical spectra. Programs are available in IDL and Python to perform calculation of the expected emergent spectrum from these sources. The database includes atomic energy levels, wavelengths, radiative transition probabilities, rate coefficients for collisional excitation, ionization, and recombination, as well as data to calculate free-free, free-bound, and two-photon continuum emission. In Version 9, we improve the modelling of the satellite lines at X-ray wavelengths by explicitly including autoionization and dielectronic recombination processes in the calculation of level populations for select members of the lithium isoelectronic sequence and Fe XVIII-XXIII. In addition, existing datasets are updated, new ions added and new total recombination rates for several Fe ions are included. All data and IDL programs are freely available at this http URL or through SolarSoft and the Python code ChiantiPy is also freely available at this https URL.

Mind the gap: The location of the lower edge of the pair instability supernovae black hole mass gap

Abstract: Gravitational-wave detections are now starting to probe the mass distribution of stellar-mass black holes (BHs). Robust predictions from stellar models are needed to interpret these. Theory predicts the existence of a gap in the BH mass distribution because of pair-instability supernova. The maximum BH mass below the gap is the result of pulsational mass loss. We evolve massive helium stars through their late hydrodynamical phases of evolution using the open-source MESA stellar evolution code. We find that the location of the lower edge of the mass gap at 45$M_\odot$ is remarkably robust against variations in the metallicity ($\approx 3M_\odot$), the treatment of internal mixing ($\approx 1M_\odot$), stellar wind mass loss ($\approx 4M_\odot$), making it the most robust predictions for the final stages of massive star evolution. The reason is that the onset of the instability is dictated by the near-final core mass, which in turn sets the resulting BH mass. However, varying $^{12}C\left(\alpha,\gamma\right)^{16}O$ reaction rate within its $1\sigma$ uncertainties shifts the location of the gap between $40M_\odot$ and $56M_\odot$. We provide updated analytic fits for population synthesis simulations. Our results imply that the detection of merging BHs can provide constraints on nuclear astrophysics. Furthermore, the robustness against metallicity suggests that there is a universal maximum for the location of the lower edge of the gap, which is insensitive to the formation environment and redshift for first-generation BHs. This is promising for the possibility to use the location of the gap as a "standard siren" across the Universe.

Kinematic detection of a planet carving a gap in a protoplanetary disc

Abstract: We still do not understand how planets form, or why extra-solar planetary systems are so different from our own solar system. But the last few years have dramatically changed our view of the discs of gas and dust around young stars. Observations with the Atacama Large Millimeter/submillimeter Array (ALMA) and extreme adaptive-optics systems have revealed that most --- if not all --- discs contain substructure, including rings and gaps, spirals, azimuthal dust concentrations, and shadows cast by misaligned inner discs. These features have been interpreted as signatures of newborn protoplanets, but the exact origin is unknown. Here we report the kinematic detection of a few Jupiter-mass planet located in a gas and dust gap at 130 au in the disc surrounding the young star HD 97048. An embedded planet can explain both the disturbed Keplerian flow of the gas, detected in CO lines, and the gap detected in the dust disc at the same radius. While gaps appear to be a common feature in protoplanetary discs, we present a direct correspondence between a planet and a dust gap, indicating that at least some gaps are the result of planet-disc interactions.

Flare-productive active regions

Abstract: Strong solar flares and coronal mass ejections, here defined not only as the bursts of electromagnetic radiation but as the entire process in which magnetic energy is released through magnetic reconnection and plasma instability, emanate from active regions (ARs) in which high magnetic non-potentiality resides in a wide variety of forms. This review focuses on the formation and evolution of flare-productive ARs from both observational and theoretical points of view. Starting from a general introduction of the genesis of ARs and solar flares, we give an overview of the key observational features during the long-term evolution in the pre-flare state, the rapid changes in the magnetic field associated with the flare occurrence, and the physical mechanisms behind these phenomena. Our picture of flare-productive ARs is summarized as follows: subject to the turbulent convection, the rising magnetic flux in the interior deforms into a complex structure and gains high non-potentiality; as the flux appears on the surface, an AR with large free magnetic energy and helicity is built, which is represented by delta-sunspots, sheared polarity inversion lines, magnetic flux ropes, etc; the flare occurs when sufficient magnetic energy has accumulated, and the drastic coronal evolution affects magnetic fields even in the photosphere. We show that the improvement of observational instruments and modeling capabilities has significantly advanced our understanding in the last decades. Finally, we discuss the outstanding issues and future perspective and further broaden our scope to the possible applications of our knowledge to space-weather forecasting, extreme events in history, and corresponding stellar activities.

The CARMENES search for exoplanets around M dwarfs: Different roads to radii and masses of the target stars

Abstract: We determine the radii and masses of 293 nearby, bright M dwarfs of the CARMENES survey. This is the first time that such a large and homogeneous high-resolution (R>80 000) spectroscopic survey has been used to derive these fundamental stellar parameters. We derived the radii using Stefan-Boltzmann's law. We obtained the required effective temperatures $T_{\rm eff}$ from a spectral analysis and we obtained the required luminosities L from integrated broadband photometry together with the Gaia DR2 parallaxes. The mass was then determined using a mass-radius relation that we derived from eclipsing binaries known in the literature. We compared this method with three other methods: (1) We calculated the mass from the radius and the surface gravity log g, which was obtained from the same spectral analysis as $T_{\rm eff}$. (2) We used a widely used infrared mass-magnitude relation. (3) We used a Bayesian approach to infer stellar parameters from the comparison of the absolute magnitudes and colors of our targets with evolutionary models. Between spectral types M0V and M7V our radii cover the range $0.1\,R_{ ormalsize\odot}

The Evolution and Role of Solar Wind Turbulence in the Inner Heliosphere

Abstract: The first two orbits of the Parker Solar Probe (PSP) spacecraft have enabled the first in situ measurements of the solar wind down to a heliocentric distance of 0.17 au (or 36 Rs). Here, we present an analysis of this data to study solar wind turbulence at 0.17 au and its evolution out to 1 au. While many features remain similar, key differences at 0.17 au include: increased turbulence energy levels by more than an order of magnitude, a magnetic field spectral index of -3/2 matching that of the velocity and both Elsasser fields, a lower magnetic compressibility consistent with a smaller slow-mode kinetic energy fraction, and a much smaller outer scale that has had time for substantial nonlinear processing. There is also an overall increase in the dominance of outward-propagating Alfvenic fluctuations compared to inward-propagating ones, and the radial variation of the inward component is consistent with its generation by reflection from the large-scale gradient in Alfven speed. The energy flux in this turbulence at 0.17 au was found to be ~10% of that in the bulk solar wind kinetic energy, becoming ~40% when extrapolated to the Alfven point, and both the fraction and rate of increase of this flux towards the Sun is consistent with turbulence-driven models in which the solar wind is powered by this flux.

Magnetic connectivity of the ecliptic plane within 0.5 AU : PFSS modeling of the first PSP encounter.

Abstract: We compare magnetic field measurements taken by the FIELDS instrument on Parker Solar Probe (PSP) during its first solar encounter to predictions obtained by Potential Field Source Surface (PFSS) modeling. Ballistic propagation is used to connect the spacecraft to the source surface. Despite the simplicity of the model, our results show striking agreement with PSPs first observations of the heliospheric magnetic field from 0.5 AU (107.5 Rs) down to 0.16 AU (35.7 Rs). Further, we show the robustness of the agreement is improved both by allowing the photospheric input to the model to vary in time, and by advecting the field from PSP down to the PFSS model domain using in situ PSP/SWEAP measurements of the solar wind speed instead of assuming it to be constant with longitude and latitude. We also explore the source surface height parameter (RSS) to the PFSS model finding that an extraordinarily low source surface height (1.3-1.5Rs) predicts observed small scale polarity inversions which are otherwise washed out with regular modeling parameters. Finally, we extract field line traces from these models. By overlaying these on EUV images we observe magnetic connectivity to various equatorial and mid-latitude coronal holes indicating plausible magnetic footpoints and offering context for future discussions of sources of the solar wind measured by PSP.

First grids of low-mass stellar models and isochrones with self-consistent treatment of rotation : From 0.2 to 1.5 M_\odot at 7 metallicities from PMS to TAMS

Abstract: We present an extended grid of state-of-the art stellar models for low-mass stars including updated physics (nuclear reaction rates, surface boundary condition, mass-loss rate, angular momentum transport, torque and rotation-induced mixing prescriptions). We aim at evaluating the impact of wind braking, realistic atmospheric treatment, rotation and rotation-induced mixing on the structural and rotational evolution from the pre-main sequence to the turn-off. Using the STAREVOL code, we provide an updated PMS grid. We compute stellar models for 7 different metallicities, from [Fe/H] = -1 dex to [Fe/H] = +0.3 dex with a solar composition corresponding to $Z=0.0134$. The initial stellar mass ranges from 0.2 to 1.5\Ms with extra grid refinement around one solar mass. We also provide rotating models for three different initial rotation rates (slow, median and fast) with prescriptions for the wind braking and disc-coupling timescale calibrated on observed properties of young open clusters. The rotational mixing includes an up-to-date description of the turbulence anisotropy in stably stratified regions. The overall behaviour of our models at solar metallicity -- and its constitutive physics -- is validated through a detailed comparison with a variety of distributed evolutionary tracks. The main differences arise from the choice of surface boundary conditions and initial solar composition. The models including rotation with our prescription for angular momentum extraction and self-consistent formalism for angular momentum transport are able to reproduce the rotation period distribution observed in young open clusters over a broad mass-range. These models are publicly available and may be used to analyse data coming from present and forthcoming asteroseismic and spectroscopic surveys such as Gaia, TESS and PLATO.

The Evolution of Massive Helium Stars Including Mass Loss

Abstract: The evolution of helium stars with initial masses in the range 1.6 to 120 Msun is studied, including the effects of mass loss by winds. These stars are assumed to form in binary systems when their expanding hydrogenic envelopes are promptly lost just after helium ignition. Significant differences are found with single star evolution, chiefly because the helium core loses mass during helium burning rather than gaining it from hydrogen shell burning. Consequently presupernova stars for a given initial mass function have considerably smaller mass when they die and will be easier to explode. Even accounting for this difference, the helium stars with mass loss develop more centrally condensed cores that should explode more easily than their single-star counterparts. The production of low mass black holes may be diminished. Helium stars with initial masses below 3.2 Msun experience significant radius expansion after helium depletion, reaching blue supergiant proportions. This could trigger additional mass exchange or affect the light curve of the supernova. The most common black hole masses produced in binaries is estimated to be about 9 Msun. A new maximum mass for black holes derived from pulsational pair-instability supernovae is derived - 46 Msun, and a new potential gap at 10 - 12 Msun is noted. Models pertinent to SN 2014ft are presented and a library of presupernova models is generated.

Clusters and mirages: cataloguing stellar aggregates in the Milky Way

Abstract: Many open clusters listed in modern catalogues were initially reported by visual astronomers as apparent overdensities of bright stars. As observational techniques and analysis methods continue to improve, some of them have been shown to be chance alignments of stars and not true clusters. Recent publications have provided membership lists for over a thousand clusters, however, many nearby objects listed in the literature have so far evaded detection. We aim to update the Gaia DR2 cluster census by performing membership determinations for known clusters that had been missed by previous studies and for recently-discovered clusters. We investigate a sub-set of non-detected clusters that, according to their literature parameters, should be easily visible in Gaia. Confirming or disproving the existence of old, inner-disc, high-altitude clusters is especially important as their survival or disruption is linked to the dynamical processes that drive the evolution of the Milky Way. We employed the Gaia DR2 catalogue and a membership assignment procedure, as well as visual inspections of spatial, proper motion, and parallax distributions. We used membership lists provided by other authors when available. We derived membership lists for 150 objects, including ten that were already known prior to Gaia. We compiled a final list of members for 1481 clusters. Among the objects that we are still unable to identify with the Gaia data, we argue that many (mostly putatively old, relatively nearby, high-altitude objects) are not true clusters. At present, the only confirmed cluster located further than 500pc away from the Galactic plane within the Solar circle is NGC 6791. It is likely that the objects discussed in this study only represent a fraction of the non-physical groupings erroneously listed in the catalogues as genuine open clusters and that those lists need further cleaning.

Do Kepler superflare stars really include slowly-rotating Sun-like stars ? - Results using APO 3.5m telescope spectroscopic observations and Gaia-DR2 data -

Abstract: We report the latest view of Kepler solar-type (G-type main-sequence) superflare stars, including recent updates with Apache Point Observatory (APO) 3.5m telescope spectroscopic observations and Gaia-DR2 data. First, we newly conducted APO3.5m spectroscopic observations of 18 superflare stars found from Kepler 1-min time cadence data. More than half (43 stars) are confirmed to be "single" stars, among 64 superflare stars in total that have been spectroscopically investigated so far in this APO3.5m and our previous Subaru/HDS observations. The measurements of $v\sin i$ (projected rotational velocity) and chromospheric lines (Ca II H\&K and Ca II 8542A) support the brightness variation of superflare stars is caused by the rotation of a star with large starspots. We then investigated the statistical properties of Kepler solar-type superflare stars by incorporating Gaia-DR2 stellar radius estimates. As a result, the maximum superflare energy continuously decreases as the rotation period $P_{\mathrm{rot}}$ increases. Superflares with energies $\lesssim 5\times10^{34}$ erg occur on old, slowly-rotating Sun-like stars ($P_{\mathrm{rot}}\sim$25 days) approximately once every 2000--3000 years, while young rapidly-rotating stars with $P_{\mathrm{rot}}\sim$ a few days have superflares up to $10^{36}$ erg. The maximum starspot area does not depend on the rotation period when the star is young, but as the rotation slows down, it starts to steeply decrease at $P_{\mathrm{rot}}\gtrsim$12 days for Sun-like stars. These two decreasing trends are consistent since the magnetic energy stored around starspots explains the flare energy, but other factors like spot magnetic structure should also be considered.

A wide star-black-hole binary system from radial-velocity measurements

Abstract: All stellar mass black holes have hitherto been identified by X-rays emitted by gas that is accreting onto the black hole from a companion star. These systems are all binaries with black holes below 30 M$_{\odot}$$^{1-4}$. Theory predicts, however, that X-ray emitting systems form a minority of the total population of star-black hole binaries$^{5,6}$. When the black hole is not accreting gas, it can be found through radial velocity measurements of the motion of the companion star. Here we report radial velocity measurements of a Galactic star, LB-1, which is a B-type star, taken over two years. We find that the motion of the B-star and an accompanying H$\alpha$ emission line require the presence of a dark companion with a mass of $68^{+11}_{-13}$ M$_{\odot}$, which can only be a black hole. The long orbital period of 78.9 days shows that this is a wide binary system. The gravitational wave experiments have detected similarly massive black holes$^{7,8}$, but forming such massive ones in a high-metallicity environment would be extremely challenging to current stellar evolution theories$^{9-11}$.

YBC, a stellar bolometric corrections database with variable extinction coefficients: an application to PARSEC isochrones.

Abstract: We present the \texttt{YBC} database of stellar bolometric corrections (BCs), available at \url{this http URL}. We homogenize widely-used theoretical stellar spectral libraries and provide BCs for many popular photometric systems, including the Gaia filters. The database can be easily extended to additional photometric systems and stellar spectral libraries. The web interface allows users to transform their catalogue of theoretical stellar parameters into magnitudes and colours of selected filter sets. The BC tables can be downloaded or also be implemented into large simulation projects using the interpolation code provided with the database. We compute extinction coefficients on a star-by-star basis, hence taking into account the effects of spectral type and non-linearity dependency on the total extinction. We illustrate the use of these BCs in \texttt{PARSEC} isochrones. We show that using spectral-type dependent extinction coefficients is quite necessary for Gaia filters whenever $\Av\gtrsim 0.5$\,mag. BC tables for rotating stars and tables of limb-darkening coefficients are also provided.

Towards precise stellar ages: combining isochrone fitting with empirical gyrochronology

Abstract: We present a new age-dating technique that combines gyrochronology with isochrone fitting to infer ages for FGKM main-sequence and subgiant field stars. Gyrochronology and isochrone fitting are each capable of providing relatively precise ages for field stars in certain areas of the Hertzsprung-Russell diagram: gyrochronology works optimally for cool main-sequence stars, and isochrone fitting can provide precise ages for stars near the main-sequence turnoff. Combined, these two age-dating techniques can provide precise and accurate ages for a broader range of stellar masses and evolutionary stages than either method used in isolation. We demonstrate that the position of a star on the Hertzsprung- Russell or color-magnitude diagram can be combined with its rotation period to infer a precise age via both isochrone fitting and gyrochronology simultaneously. We show that incorporating rotation periods with 5% uncertainties into stellar evolution models improves age precision for FGK stars on the main sequence, and can, on average, provide age estimates up to three times more precise than isochrone fitting alone. In addition, we provide a new gyrochronology relation, calibrated to the Praesepe cluster and the Sun, that includes a variance model to capture the rotational behavior of stars whose rotation periods do not lengthen with the square-root of time, and parts of the Hertzsprung-Russell diagram where gyrochronology has not been calibrated. This publication is accompanied by an open source Python package, stardate, for inferring the ages of main-sequence and subgiant FGKM stars from rotation periods, spectroscopic parameters and/or apparent magnitudes and parallaxes.

The Asteroseismic Target List (ATL) for solar-like oscillators observed in 2-minute cadence with the Transiting Exoplanet Survey Satellite (TESS)

Abstract: We present the target list of solar-type stars to be observed in short-cadence (2-min) for asteroseismology by the NASA Transiting Exoplanet Survey Satellite (TESS) during its 2-year nominal survey mission. The solar-like Asteroseismic Target List (ATL) is comprised of bright, cool main-sequence and subgiant stars and forms part of the larger target list of the TESS Asteroseismic Science Consortium (TASC). The ATL uses Gaia DR2 and the Extended Hipparcos Compilation (XHIP) to derive fundamental stellar properties, calculate detection probabilities and produce a rank-ordered target list. We provide a detailed description of how the ATL was produced and calculate expected yields for solar-like oscillators based on the nominal photometric performance by TESS. We also provide publicly available source code which can be used to reproduce the ATL, thereby enabling comparisons of asteroseismic results from TESS with predictions from synthetic stellar populations.

How mergers magnetise massive stars

Abstract: Magnetic fields are ubiquitous in the Universe. The Sun's magnetic field drives the solar wind and causes solar flares and other energetic surface phenomena that profoundly affect space weather here on Earth. The first magnetic field in a star other than the Sun was detected in 1947 in the peculiar A-type star 78 Vir. It is now known that the magnetic fields of the Sun and other low-mass stars ( 1.5 solar masses; referred to as "massive" stars here) have relatively quiet, radiative envelopes where a solar-like dynamo cannot operate. However, about 10% of them, including 78 Vir, have strong, large-scale surface magnetic fields whose origin has remained a major mystery. The massive star $\tau$ Sco is a prominent member of this group and appears to be surprisingly young compared to other presumably coeval members of the Upper Scorpius association. Here, we present the first 3D magneto-hydrodynamical simulations of the coalescence of two massive main-sequence stars and 1D stellar evolution computations of the subsequent evolution of the merger product that can explain $\tau$ Sco's magnetic field, apparent youth and other observed characteristics. We argue that field amplification in stellar mergers is a general mechanism to form strongly-magnetised massive stars. These stars are promising progenitors of those neutron stars that host the strongest magnetic fields in the Universe, so-called magnetars, and that may give rise to some of the enigmatic fast radio bursts. Strong magnetic fields affect the explosions of core-collapse supernovae and, moreover, those magnetic stars that have rapidly-rotating cores at the end of their lives might provide the right conditions to power long-duration gamma-ray bursts and super-luminous supernovae.

A stellar census in globular clusters with MUSE: Binaries in NGC 3201

Abstract: We utilize multi-epoch MUSE spectroscopy to study binaries in the core of NGC 3201. Our sample consists of 3553 stars with 54883 spectra in total comprising 3200 main-sequence stars up to 4 magnitudes below the turn-off. Each star in our sample has between 3 and 63 (with a median of 14) reliable radial velocity (RV) measurements within five years of observations. We introduce a statistical method to determine the probability of a star showing RV variations based on the whole inhomogeneous RV sample. Using HST photometry and an advanced dynamical MOCCA simulation of this specific GC we overcome observational biases that previous spectroscopic studies had to deal with. This allows us to infer a binary frequency in the MUSE FoV and enables us to deduce the underlying true binary frequency of (6.75+-0.72) % in NGC 3201. The comparison of the MUSE observations with the MOCCA simulation suggests a significant fraction of primordial binaries. We can also confirm a radial increase of the binary fraction towards the GC centre due to mass segregation. We discovered that in our sample at least (57.5+-7.9) % of blue straggler stars (BSS) are in a binary system. For the first time in a study of GCs, we were able to fit Keplerian orbits to a significant sample of 95 binaries. We present the binary system properties of eleven BSS and show evidence that two BSS formation scenarios, the mass transfer in binary (or triple) star systems and the coalescence due to binary-binary interactions, are present in our data. We also describe the binary and spectroscopic properties of four sub-subgiant (or red straggler) stars. Furthermore, we discovered two new black hole (BH) candidates with minimum masses (Msini) of (7.68+-0.50) M_sun, (4.4+-2.8) M_sun, and refine the minimum mass estimate on the already published BH to (4.53+-0.21) M_sun. These BHs are consistent with an extensive BH subsystem hosted by NGC 3201.

Magnetic fields in M dwarfs from the CARMENES survey.

Abstract: M dwarfs are known to generate the strongest magnetic fields among main-sequence stars with convective envelopes, but the link between the magnetic fields and underlying dynamo mechanisms, rotation, and activity still lacks a consistent picture. In this work we measure magnetic fields from the high-resolution near-infrared spectra taken with the CARMENES radial-velocity planet survey in a sample of 29 active M dwarfs and compare our results against stellar parameters. We use the state-of-the-art radiative transfer code to measure total magnetic flux densities from the Zeeman broadening of spectral lines and filling factors. We detect strong kG magnetic fields in all our targets. In 16 stars the magnetic fields were measured for the first time. Our measurements are consistent with the magnetic field saturation in stars with rotation periods P<4d. The analysis of the magnetic filling factors reveal two different patterns of either very smooth distribution or a more patchy one, which can be connected to the dynamo state of the stars and/or stellar mass. Our measurements extend the list of M dwarfs with strong surface magnetic fields. They also allow us to better constrain the interplay between the magnetic energy, stellar rotation, and underlying dynamo action. The high spectral resolution and observations at near-infrared wavelengths are the beneficial capabilities of the CARMENES instrument that allow us to address important questions about the stellar magnetism.

A White Dwarf with Transiting Circumstellar Material Far Outside the Roche Limit

Abstract: We report the discovery of a white dwarf exhibiting deep, irregularly shaped transits, indicative of circumstellar planetary debris. Using Zwicky Transient Facility DR2 photometry of ZTF$\,$J013906.17+524536.89 and follow-up observations from the Las Cumbres Observatory, we identify multiple transit events that recur every ${\approx}\,107.2\,$d, much longer than the $4.5{-}4.9\,$h orbital periods observed in WD$\,$1145+017, the only other white dwarf known with transiting planetary debris. The transits vary in both depth and duration, lasting $15{-}25\,$d and reaching $20{-}45\,\%$ dips in flux. Optical spectra reveal strong Balmer lines, identifying the white dwarf as a DA with $T_{\mathrm{eff}}=10{,}530\pm140\,\mathrm{K}$ and $\log(g)=7.86\pm0.06$. A $\mathrm{Ca\,II\,K}$ absorption feature is present in all spectra both in and out of transit. Spectra obtained during one night at roughly $15\,\%$ transit depth show increased $\mathrm{Ca\,II\,K}$ absorption with a model atmospheric fit suggesting $[\mathrm{Ca/H}]=-4.6\pm0.3$, whereas spectra taken on three nights out of transit have $[\mathrm{Ca/H}]$ of -5.5, -5.3, and -4.9 with similar uncertainties. While the $\mathrm{Ca\,II\,K}$ line strength varies by only 2-sigma, we consider a predominantly interstellar origin for Ca absorption unlikely. We suggest a larger column density of circumstellar metallic gas along the line of site or increased accretion of material onto the white dwarf's surface are responsible for the Ca absorption, but further spectroscopic studies are required. In addition, high-speed time series photometry out of transit reveals variability with periods of 900 and 1030$\,$s, consistent with ZZ Ceti pulsations.

The CARMENES search for exoplanets around M dwarfs -- Photospheric parameters of target stars from high-resolution spectroscopy. II. Simultaneous multiwavelength range modeling of activity insensitive lines

Abstract: We present precise photospheric parameters of 282 M dwarfs determined from fitting the most recent version of PHOENIX models to high-resolution CARMENES spectra in the visible (0.52 - 0.96 $\mu$m) and near-infrared wavelength range (0.96 - 1.71 $\mu$m). With its aim to search for habitable planets around M dwarfs, several planets of different masses have been detected. The characterization of the target sample is important for the ability to derive and constrain the physical properties of any planetary systems that are detected. As a continuation of previous work in this context, we derived the fundamental stellar parameters effective temperature, surface gravity, and metallicity of the CARMENES M-dwarf targets from PHOENIX model fits using a $\chi^2$ method. We calculated updated PHOENIX stellar atmosphere models that include a new equation of state to especially account for spectral features of low-temperature stellar atmospheres as well as new atomic and molecular line lists. We show the importance of selecting magnetically insensitive lines for fitting to avoid effects of stellar activity in the line profiles. For the first time, we directly compare stellar parameters derived from multiwavelength range spectra, simultaneously observed for the same star. In comparison with literature values we show that fundamental parameters derived from visible spectra and visible and near-infrared spectra combined are in better agreement than those derived from the same spectra in the near-infrared alone.

Bright C2H emission in protoplanetary disks in Lupus: high volatile C/O>1 ratios

Abstract: ALMA surveys have shown that CO emission in protoplanetary disks is much fainter than expected. Accordingly, CO-based gas masses and gas/dust ratios are orders of magnitude lower than previously thought. This may be explained either as fast gas dispersal, or as chemical evolution and locking up of volatiles in larger bodies leading to the low observed CO fluxes. The latter processes lead to enhanced C/O ratios in the gas, which may be reflected in enhanced abundances of molecules like C2H. The goal of this work is to employ C2H observations to understand whether low CO fluxes are caused by volatile depletion, or by fast gas dissipation. We present ALMA Cycle 4 C2H observations of a subsample of nine sources in Lupus. The integrated C2H emission is determined and compared to previous CO observations and model predictions. Seven out of nine disks are detected in C2H, whose line emission is almost as bright as 13CO. All detections are significantly brighter than the typical sensitivity of the observations, hinting at a bimodal distribution of the C2H line intensities. When compared with DALI models, the observed C2H fluxes can be reproduced only if some level of volatile carbon and oxygen depletion is allowed and C/O>1 in the gas. Models with reduced gas/dust ratios fail instead to reproduce the observations. A steeper than linear correlation between C2H and CN emission line is found. This is linked to the fact that C2H emission lines are affected more strongly by C/O variations than CN lines. Ring-like structures are detected both in C2H and in continuum emission but, as for CN, they do not seem to be connected. Sz 71 shows ring shaped emission in both C2H and CN with the location of the peak intensity coinciding. Our new ALMA C2H observations favour volatile carbon and oxygen depletion rather than fast gas dispersal to explain the faint CO observations for most of the disks.

Prediction of the strength and timing of sunspot cycle 25 reveal decadal-scale space environmental conditions

Abstract: The Sun's activity cycle governs the radiation, particle and magnetic flux in the heliosphere creating hazardous space weather. Decadal-scale variations define space climate and force the Earth's atmosphere. However, predicting the solar cycle is challenging. Current understanding indicates a short window for prediction best achieved at previous cycle minima. Utilizing magnetic field evolution models for the Sun's surface and interior we perform the first century-scale, data-driven simulations of solar activity and present a scheme for extending the prediction window to a decade. Our ensemble forecast indicates cycle 25 would be similar or slightly stronger than the current cycle and peak around 2024. Sunspot cycle 25 may thus reverse the substantial weakening trend in solar activity which has led to speculation of an imminent Maunder-like grand minimum and cooling global climate. Our simulations demonstrate fluctuation in the tilt angle distribution of sunspots is the dominant mechanism responsible for solar cycle variability.

New study of the line profiles of sodium perturbed by H2

Abstract: The opacity of alkali atoms, most importantly of Na and K, plays a crucial role in the atmospheres of brown dwarfs and exoplanets. We present a comprehensive study of NaH2 collisional profiles at temperatures from 500 to 3000 K, the temperatures prevailing in the atmosphere of brown dwarfs and Jupiter-mass planets.The relevant H2 perturber densities reach several 10^19 cm^-3 in hot Teff > 1500 K Jupiter-mass planets and can exceed 10^20 cm^-3 for more massive or cooler objects. Accurate pressure-broadened profiles that are valid at high densities of H2 should be incorporated into spectral models. Unified profiles of sodium perturbed by molecular hydrogen were calculated in the semi-classical approach using up-to-date molecular data.New NaH2 collisional profiles and their effects on the synthetic spectra of brown dwarfs and hot Jupiters computed with petitCODE are presented.

Abundance Estimates for 16 Elements in 6 Million Stars from LAMOST DR5 Low-Resolution Spectra

Abstract: We present the determination of stellar parameters and individual elemental abundances for 6 million stars from $\sim$8 million low-resolution ($R\sim1800$) spectra from LAMOST DR5. This is based on a modeling approach that we dub $The$ $Data$--$Driven$ $Payne$ ($DD$--$Payne$), which inherits essential ingredients from both {\it The Payne} \citep{Ting2019} and $The$ $Cannon$ \citep{Ness2015}. It is a data-driven model that incorporates constraints from theoretical spectral models to ensure the derived abundance estimates are physically sensible. Stars in LAMOST DR5 that are in common with either GALAH DR2 or APOGEE DR14 are used to train a model that delivers stellar parameters ($T_{\rm eff}$, $\log g$, $V_{\rm mic}$) and abundances for 16 elements (C, N, O, Na, Mg, Al, Si, Ca, Ti, Cr, Mn, Fe, Co, Ni, Cu, and Ba) when applied to LAMOST spectra. Cross-validation and repeat observations suggest that, for ${\rm S/N}_{\rm pix}\ge 50$, the typical internal abundance precision is 0.03--0.1\,dex for the majority of these elements, with 0.2--0.3\,dex for Cu and Ba, and the internal precision of $T_{\rm eff}$ and $\log g$ is better than 30\,K and 0.07\,dex, respectively. Abundance systematics at the $\sim$0.1\,dex level are present in these estimates, but are inherited from the high-resolution surveys' training labels. For some elements, GALAH provides more robust training labels, for others, APOGEE. We provide flags to guide the quality of the label determination and to identify binary/multiple stars in LAMOST DR5. The abundance catalogs are publicly accessible via \href{url}{this http URL}.

Most Black Holes are Born Very Slowly Rotating

Abstract: The age of gravitational wave (GW) astronomy has begun, and black hole (BH) mergers detected by LIGO are providing novel constraints on massive star evolution. A major uncertainty in stellar theory is the angular momentum (AM) transport within the star that determines its core rotation rate and the resulting BH's spin. Internal rotation rates of low-mass stars measured from asteroseismology prove that AM transport is efficient, suggesting that massive stellar cores may rotate slower than prior expectations. We investigate AM transport via the magnetic Tayler instability, which can largely explain the rotation rates of low-mass stars and white dwarfs. Implementing an updated AM transport prescription into models of high-mass stars, we compute the spins of their BH remnants. We predict that BHs born from single stars rotate very slowly, with $a \sim 10^{-2}$, regardless of initial rotation rate, possibly explaining the low $\chi_{\rm eff}$ of most BH binaries detected by LIGO thus far. A limited set of binary models suggests slow rotation for many binary scenarios as well, although homogeneous evolution and tidal spin-up of post-common envelope helium stars can create moderate or high BH spins. We make predictions for the values of $\chi_{\rm eff}$ in future LIGO events, and we discuss implications for engine-powered transients.

On the competing effect of wind braking and interior coupling in the rotational evolution of solar-like stars

Abstract: Solar-like stars (M < 1.3 Msun) lose angular momentum through their magnetized winds. The resulting evolution of the surface rotation period, which can be directly measured photometrically, has the potential to provide an accurate indicator of stellar age, and is constrained by observations of rotation periods of coeval stars, such as members of Galactic open clusters. A prominent observational feature of the mass-rotation period diagrams of open clusters is a sequence of relatively slower rotators. The formation and persistence of this slow rotators sequence across several billion years imply an approximately coherent spin-down of the stars that belong to it. In particular, the sequence is observed to evolve coherently towards longer periods in progressively older clusters. Recent observations of the 700-Myr old Praesepe and the 1-Gyr old NGC 6811 clusters, however, seem to contradict this general pattern. While the 1 Msun stars on the slow rotators sequence of the older NGC 6811 have longer periods than their counterparts in the younger Praesepe, as expected, the two sequences essentially merge at lower masses (< 0.8 Msun). In other words, low-mass stars seem to have not been spinning down at all in the intervening 300 Myr. Here we show that this behavior is a manifestation of the variable rotational coupling in solar-like stars. The resurfacing of angular momentum from the interior can temporarily compensate for the loss due to wind braking at the surface. The internal redistribution of angular momentum has a steep mass dependence; as a result, the re-coupling occurs at different ages for stars of different masses. Our model explains the stalled surface spin-down of low-mass stars between Praesepe and NGC 6811, and predicts that the same behavior should be observable at other ages in other mass ranges.

The Solar Neighborhood XLV. The Stellar Multiplicity Rate of M Dwarfs Within 25 pc

Abstract: We present results of the largest, most comprehensive study ever done of the stellar multiplicity of the most common stars in the Galaxy, the red dwarfs. We have conducted an all-sky, volume-limited survey for stellar companions to 1120 M dwarf primaries known to lie within 25 pc of the Sun via trigonometric parallaxes. In addition to a comprehensive literature search, stars were explored in new surveys for companions at separations of 2" to 300". A reconnaissance of wide companions to separations of 300" was done via blinking archival images. I-band images were used to search our sample for companions at separations of 2" to 180". Various astrometric and photometric methods were used to probe the inner 2" to reveal close companions. We report the discovery of 20 new companions and identify 56 candidate multiple systems. We find a stellar multiplicity rate of 26.8 +/- 1.4% and a stellar companion rate of 32.4 +/- 1.4% for M dwarfs. There is a broad peak in the separation distribution of the companions at 4 -- 20 AU, with a weak trend of smaller projected linear separations for lower mass primaries. A hint that M dwarf multiplicity may be a function of tangential velocity is found, with faster moving, presumably older, stars found to be multiple somewhat less often. We calculate that stellar companions make up at least 17% of mass attributed to M dwarfs in the solar neighborhood, with roughly 11% of M dwarf mass hidden as unresolved companions. Finally, when considering all M dwarf primaries and companions, we find that the mass distribution for M dwarfs increases to the end of the stellar main sequence.

A Machine Learning Dataset Prepared From the NASA Solar Dynamics Observatory Mission

Abstract: In this paper we present a curated dataset from the NASA Solar Dynamics Observatory (SDO) mission in a format suitable for machine learning research. Beginning from level 1 scientific products we have processed various instrumental corrections, downsampled to manageable spatial and temporal resolutions, and synchronized observations spatially and temporally. We illustrate the use of this dataset with two example applications: forecasting future EVE irradiance from present EVE irradiance and translating HMI observations into AIA observations. For each application we provide metrics and baselines for future model comparison. We anticipate this curated dataset will facilitate machine learning research in heliophysics and the physical sciences generally, increasing the scientific return of the SDO mission. This work is a direct result of the 2018 NASA Frontier Development Laboratory Program. Please see the appendix for access to the dataset.