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

MESA Isochrones and Stellar Tracks (MIST). I: Solar-Scaled Models

Abstract: This is the first of a series of papers presenting the Modules for Experiments in Stellar Astrophysics (MESA) Isochrones and Stellar Tracks (MIST) project, a new comprehensive set of stellar evolutionary tracks and isochrones computed using MESA, a state-of-the-art open-source 1D stellar evolution package. In this work, we present models with solar-scaled abundance ratios covering a wide range of ages ($5 \leq \rm \log(Age)\;[yr] \leq 10.3$), masses ($0.1 \leq M/M_{\odot} \leq 300$), and metallicities ($-2.0 \leq \rm [Z/H] \leq 0.5$). The models are self-consistently and continuously evolved from the pre-main sequence to the end of hydrogen burning, the white dwarf cooling sequence, or the end of carbon burning, depending on the initial mass. We also provide a grid of models evolved from the pre-main sequence to the end of core helium burning for $-4.0 \leq \rm [Z/H] < -2.0$. We showcase extensive comparisons with observational constraints as well as with some of the most widely used existing models in the literature. The evolutionary tracks and isochrones can be downloaded from the project website at this http URL

MESA Isochrones and Stellar Tracks (MIST) 0: Methods for the construction of stellar isochrones

Abstract: I describe a method to transform a set of stellar evolution tracks onto a uniform basis and then interpolate within that basis to construct stellar isochrones. The method accommodates a broad range of stellar types, from substellar objects to high-mass stars, and phases of evolution, from the pre-main sequence to the white dwarf cooling sequence. I discuss situations in which stellar physics leads to departures from the otherwise monotonic relation between initial stellar mass and lifetime and how these may be dealt with in isochrone construction. I close with convergence tests and recommendations for the number of points in the uniform basis and the mass between tracks in the original grid required in order to achieve a certain level of accuracy in the resulting isochrones. The programs that implement these methods are free and open-source; they may be obtained from the project webpage.

Physics of Core-Collapse Supernovae in Three Dimensions: a Sneak Preview

Abstract: Nonspherical mass motions are a generic feature of core-collapse supernovae, and hydrodynamic instabilities play a crucial role for the explosion mechanism. First successful neutrino-driven explosions could be obtained with self-consistent, first-principle simulations in three spatial dimensions (3D). But 3D models tend to be less prone to explosion than corresponding axisymmetric (2D) ones. This has been explained by 3D turbulence leading to energy cascading from large to small spatial scales, inversely to the 2D case, thus disfavoring the growth of buoyant plumes on the largest scales. Unless the inertia to explode simply reflects a lack of sufficient resolution in relevant regions, it suggests that some important aspect may still be missing for robust and sufficiently energetic neutrino-powered explosions. Such deficits could be associated with progenitor properties like rotation, magnetic fields or pre-collapse perturbations, or with microphysics that could lead to an enhancement of neutrino heating behind the shock. 3D simulations have also revealed new phenomena that are not present in 2D, for example spiral modes of the standing accretion shock instability (SASI) and a stunning dipolar lepton-emission self-sustained asymmetry (LESA). Both impose time- and direction-dependent variations on the detectable neutrino signal. The understanding of these effects and of their consequences is still in its infancy.

Pulkovo Compilation of Radial Velocities for 35495 Hipparcos Stars in a Common System

Abstract: The Pulkovo Compilation of Radial Velocities (PCRV) has been made to study the stellar kinematics in the local spiral arm. The PCRV contains weighted mean absolute radial velocities for 35 495 Hipparcos stars of various spectral types and luminosity classes over the entire celestial sphere mainly within 500 pc of the Sun. The median accuracy of the radial velocities obtained is 0.7 km s$^{-1}$. Results from 203 publications were used in the catalogue. Four of them were used to improve the radial velocities of standard stars from the IAU list. The radial velocities of 155 standard stars turned out to be constant within 0.3 km s$^{-1}$. These stars were used to analyze 47 768 mean radial velocities for 37 200 stars from 12 major publications ($\sim80\%$ of all the data used). Zero-point discrepancies and systematic dependences on radial velocity, $(B-V)$ color index, right ascension, and declination were found in radial velocity differences of the form "publication minus IAU list of standards". These discrepancies and dependences were approximated and taken into account when calculating the weighted mean radial velocities. 1128 stars whose independent radial-velocity determinations were available at least in three of these publications and agreed within 3 km s$^{-1}$ were chosen as the work list of secondary standards. Radial-velocity differences of the form "publication minus list of secondary standards" were used by analogy to correct the zero points and systematic dependences in the radial velocities from 33 more publications (~13\% of the data used). In addition, the radial velocities from 154 minor publications ($\sim7\%$ of the data used) pertaining to well-known instruments were used without any corrections.

A new route towards merging massive black holes

Abstract: Recent advances in gravitational-wave astronomy make the direct detection of gravitational waves from the merger of two stellar-mass compact objects a realistic prospect. Evolutionary scenarios towards mergers of double compact objects generally invoke common-envelope evolution which is poorly understood, leading to large uncertainties in merger rates. We explore the alternative scenario of massive overcontact binary (MOB) evolution, which involves two very massive stars in a very tight binary which remain fully mixed due to their tidally induced high spin. We use the public stellar-evolution code MESA to systematically study this channel by means of detailed simulations. We find that, at low metallicity, MOBs produce double-black-hole (BH+BH) systems that will merge within a Hubble time with mass ratios close to one, in two mass ranges, ~25...60msun and >~ 130msun, with pair instability supernovae (PISNe) being produced in-between. Our models are also able to reproduce counterparts of various stages in the MOB scenario in the local Universe, providing direct support for it. We map the initial parameter space that produces BH+BH mergers, determine the expected chirp mass distribution, merger times, Kerr parameters and predict event rates. We typically find that for Z~

Spectral Properties of Cool Stars: Extended Abundance Analysis of 1617 Planet Search Stars

Abstract: We present a catalog of uniformly determined stellar properties and abundances for 1626 F, G, and K stars using an automated spectral synthesis modeling procedure. All stars were observed using the HIRES spectrograph at Keck Observatory. Our procedure used a single line list to fit model spectra to observations of all stars to determine effective temperature, surface gravity, metallicity, projected rotational velocity, and the abundances of 15 elements (C, N, O, Na, Mg, Al, Si, Ca, Ti, V, Cr, Mn, Fe, Ni, & Y). Sixty percent of the sample had Hipparcos parallaxes and V-band photometry which we combined with the spectroscopic results to obtain mass, radius, and luminosity. Additionally, we used the luminosity, effective temperature, metallicity and alpha-element enhancement to interpolate in the Yonsei-Yale isochrones to derive mass, radius, gravity, and age ranges for those stars. Finally, we determined new relations between effective temperature and macroturbulence for dwarfs and subgiants. Our analysis achieved precisions of 25 K in Teff , 0.01 dex in [M/H], 0.028 dex for log g and 0.5 km/s in v sin ibased on multiple observations of the same stars. The abundance results were similarly precise, between 0.01 and - 0.04 dex, though trends with respect to Teff remained for which we derived empirical corrections. The trends, though small, were much larger than our uncertainties and are shared with published abundances. We show that changing our model atmosphere grid accounts for most of the trend in [M/H] between 5000 K and 5500 K indicating a possible problem with the atmosphere models or opacities.

The Gould's Belt Distances Survey (GOBELINS) II. Distances and structure towards the Orion Molecular Clouds

Abstract: We present the results of the Gould's Belt Distances Survey (GOBELINS) of young star forming regions towards the Orion Molecular Cloud Complex. We detected 36 YSOs with the Very Large Baseline Array (VLBA), 27 of which have been observed in at least 3 epochs over the course of 2 years. At least half of these YSOs belong to multiple systems. We obtained parallax and proper motions towards these stars to study the structure and kinematics of the Complex. We measured a distance of 388$\pm$5 pc towards the Orion Nebula Cluster, 428$\pm$10 pc towards the southern portion L1641, 388$\pm$10 pc towards NGC 2068, and roughly $\sim$420 pc towards NGC 2024. Finally, we observed a strong degree of plasma radio scattering towards $\lambda$ Ori.

The MUSCLES Treasury Survey I: Motivation and Overview

Abstract: Ground- and space-based planet searches employing radial velocity techniques and transit photometry have detected thousands of planet-hosting stars in the Milky Way. The chemistry of these atmospheres is controlled by the shape and absolute flux of the stellar spectral energy distribution, however, flux distributions of relatively inactive low-mass stars are poorly known at present. To better understand exoplanets orbiting low-mass stars, we have executed a panchromatic (X-ray to mid-IR) study of the spectral energy distributions of 11 nearby planet hosting stars, the {\it Measurements of the Ultraviolet Spectral Characteristics of Low-mass Exoplanetary Systems} (MUSCLES) Treasury Survey. The MUSCLES program consists of contemporaneous observations at X-ray, UV, and optical wavelengths. We show that energetic radiation (X-ray and ultraviolet) is present from magnetically active stellar atmospheres at all times for stars as late as M5. Emission line luminosities of \ion{C}{4} and \ion{Mg}{2} are strongly correlated with band-integrated luminosities. We find that while the slope of the spectral energy distribution, FUV/NUV, increases by approximately two orders of magnitude form early K to late M dwarfs ($\approx$~0.01~to~1), the absolute FUV and XUV flux levels at their corresponding habitable zone distances are constant to within factors of a few, spanning the range 10~--~70 erg cm$^{-2}$ s$^{-1}$ in the habitable zone. Despite the lack of strong stellar activity indicators in their optical spectra, several of the M dwarfs in our sample show spectacular flare emission in their UV light curves. Finally, we interpret enhanced $L(line)$/$L_{Bol}$ ratios for \ion{C}{4} and \ion{N}{5} as tentative observational evidence for the interaction of planets with large planetary mass-to-orbital distance ratios ($M_{plan}$/$a_{plan}$) with the transition regions of their host stars.

The H$\alpha$ emission of nearby M dwarfs and its relation to stellar rotation

Abstract: The high-energy emission from low-mass stars is mediated by the magnetic dynamo. Although the mechanisms by which fully convective stars generate large-scale magnetic fields are not well understood, it is clear that, as for solar-type stars, stellar rotation plays a pivotal role. We present 270 new optical spectra of low-mass stars in the Solar Neighborhood. Combining our observations with those from the literature, our sample comprises 2202 measurements or non-detections of H$\alpha$ emission in nearby M dwarfs. This includes 466 with photometric rotation periods. Stars with masses between 0.1 and 0.6 solar masses are well-represented in our sample, with fast and slow rotators of all masses. We observe a threshold in the mass-period plane that separates active and inactive M dwarfs. The threshold coincides with the fast-period edge of the slowly rotating population, at approximately the rotation period at which an era of rapid rotational evolution appears to cease. The well- defined active/inactive boundary indicates that H$\alpha$ activity is a useful diagnostic for stellar rotation period, e.g. for target selection for exoplanet surveys, and we present a mass-period relation for inactive M dwarfs. We also find a significant, moderate correlation between $L_{\mathrm{H}\alpha}/L_{\mathrm{bol}}$ and variability amplitude: more active stars display higher levels of photometric variability. Consistent with previous work, our data show that rapid rotators maintain a saturated value of $L_{\mathrm{H}\alpha}/L_{\mathrm{bol}}$. Our data also show a clear power-law decay in $L_{\mathrm{H}\alpha}/L_{\mathrm{bol}}$ with Rossby number for slow rotators, with an index of $-1.7 \pm 0.1$.

Gravitational Instabilities in Circumstellar Disks

Abstract: [Abridged] Star and planet formation are the complex outcomes of gravitational collapse and angular momentum transport mediated by protostellar and protoplanetary disks. In this review we focus on the role of gravitational instability in this process. We begin with a brief overview of the observational evidence for massive disks that might be subject to gravitational instability, and then highlight the diverse ways in which the instability manifests itself in protostellar and protoplanetary disks: the generation of spiral arms, small scale turbulence-like density fluctuations, and fragmentation of the disk itself. We present the analytic theory that describes the linear growth phase of the instability, supplemented with a survey of numerical simulations that aim to capture the non-linear evolution. We emphasize the role of thermodynamics and large scale infall in controlling the outcome of the instability. Despite apparent controversies in the literature, we show a remarkable level of agreement between analytic predictions and numerical results. We highlight open questions related to (1) the development of a turbulent cascade in thin disks, and (2) the role of mode-mode coupling in setting the maximum angular momentum transport rate in thick disks.

The Kepler Follow-Up Observation Program. I. A Catalog of Companions to Kepler Stars from High-Resolution Imaging

Abstract: We present results from high-resolution, optical to near-IR imaging of host stars of Kepler Objects of Interest (KOIs), identified in the original Kepler field. Part of the data were obtained under the Kepler imaging follow-up observation program over seven years (2009 - 2015). Almost 90% of stars that are hosts to planet candidates or confirmed planets were observed. We combine measurements of companions to KOI host stars from different bands to create a comprehensive catalog of projected separations, position angles, and magnitude differences for all detected companion stars (some of which may not be bound). Our compilation includes 2297 companions around 1903 primary stars. From high-resolution imaging, we find that ~10% (~30%) of the observed stars have at least one companion detected within 1" (4"). The true fraction of systems with close (< ~4") companions is larger than the observed one due to the limited sensitivities of the imaging data. We derive correction factors for planet radii caused by the dilution of the transit depth: assuming that planets orbit the primary stars or the brightest companion stars, the average correction factors are 1.06 and 3.09, respectively. The true effect of transit dilution lies in between these two cases and varies with each system. Applying these factors to planet radii decreases the number of KOI planets with radii smaller than 2 R_Earth by ~2-23% and thus affects planet occurrence rates. This effect will also be important for the yield of small planets from future transit missions such as TESS.

Measuring Turbulence in TW Hya with ALMA: Methods and Limitations

Abstract: We obtain high spatial and spectral resolution images of the CO J=2-1, CN N=2-1 and CS J=5-4 emission with ALMA in Cycle~2. The radial distribution of the turbulent broadening is derived with three approaches: two `direct' and one modelling. The first requires a single transition and derives \Tex{} directly from the line profile, yielding a \vturb{}. The second assumes two different molecules are co-spatial thus their relative linewidths allow for a calculation of \Tkin{} and \vturb{}. Finally we fit a parametric disk model where physical properties of the disk are described by power laws, to compare our `direct' methods with previous values. The two direct methods were limited to the outer $r > 40$~au disk due to beam smear. The direct method found \vturb{} ranging from $\approx$~\vel{130} at 40~au, dropping to $\approx$~\vel{50} in the outer disk, qualitatively recovered with the parametric model fitting. This corresponds to roughly $0.2 - 0.4~c_s$. CN was found to exhibit strong non-LTE effects outside $r \approx 140$~au, so \vturb{} was limited to within this radius. The assumption that CN and CS are co-spatial is consistent with observed linewidths only within $r \lesssim 100$~au, within which \vturb{} was found to drop from \vel{100} ($\approx~0.4~c_s$) to nothing at 100~au. The parametric model yielded a near constant \vel{50} for CS ($0.2 - 0.4~c_s$). We demonstrate that absolute flux calibration is and will be the limiting factor in all studies of turbulence using a single molecule. The magnitude of the dispersion is comparable with or below that predicted by the magneto-rotational instability theory. A more precise comparison would require to reach an absolute calibration precision of order 3\%, or to find a suitable combination of light and heavy molecules which are co-located in the disk.

Hot subluminous stars

Abstract: Hot subluminous stars of spectral type B and O are core helium-burning stars at the blue end of the horizontal branch or have evolved even beyond that stage. Strikingly, the distribution in the Hertzsprung-Russell diagram of He-rich vs. He-poor hot subdwarf stars of the globular clusters omega Cen and NGC~2808 differ from that of their field counterparts. The metal-abundance patterns of hot subdwarfs are typically characterized by strong deficiencies of some lighter elements as well as large enrichments of heavy elements. A large fraction of sdB stars are found in close binaries with white dwarf or very low-mass main sequence companions, which must have gone through a common-envelope phase of evolution.They provide a clean-cut laboratory to study this important but yet purely understood phase of stellar evolution. Substellar companions to sdB stars have also been found. For HW~Vir systems the companion mass distribution extends from the stellar into the brown dwarf regime. A giant planet to the pulsator V391 Peg was the first discovery of a planet that survived the red giant evolution of its host star. Several types of pulsating star have been discovered among hot subdwarf stars, the most common are the gravity-mode sdB pulsators (V1093 Her) and their hotter siblings, the p-mode pulsating V361 Hya stars. Another class of multi-periodic pulsating hot subdwarfs has been found in the globular cluster omega Cen that is unmatched by any field star. The masses of hot subdwarf stars are the key to understand the stars' evolution. A few pulsating sdB stars in eclipsing binaries have been found that allow mass determination. The results are in good agreement with predictions from binary population synthesis. New classes of binaries, hosting extremely low mass (ELM) white dwarfs (M<0.3 Msun), have recently been discovered, filling a gap in the mosaic of binary stellar evolution. (abbreviated)

Ambipolar diffusion in low-mass star formation. I. General comparison with the ideal MHD case

Abstract: We thank the anonymous referee for the suggestions and remarks that contributed to improve the quality of this manuscript. The research leading to these results has received funding from the European Research Council under the European Community’s Seventh Framework Programme (FP7/2007-2013 Grant Agreement no. 247060). B.C. gratefully acknowledges support from the French ANR Retour Postdoc program (ANR-11-PDOC-0031). We finally acknowledge financial support from the “Programme National de Physique Stellaire” (PNPS) of CNRS/INSU, France.

Nominal values for selected solar and planetary quantities: IAU 2015 Resolution B3

Abstract: In this brief communication we provide the rationale for, and the outcome of the International Astronomical Union (IAU) resolution vote at the XXIX-th General Assembly in Honolulu, Hawaii, in 2015, on recommended nominal conversion constants for selected solar and planetary properties. The problem addressed by the resolution is a lack of established conversion constants between solar and planetary values and SI units: a missing standard has caused a proliferation of solar values (e.g., solar radius, solar irradiance, solar luminosity, solar effective temperature and solar mass parameter) in the literature, with cited solar values typically based on best estimates at the time of paper writing. As precision of observations increases, a set of consistent values becomes increasingly important. To address this, an IAU Working Group on Nominal Units for Stellar and Planetary Astronomy formed in 2011, uniting experts from the solar, stellar, planetary, exoplanetary and fundamental astronomy, as well as from general standards fields to converge on optimal values for nominal conversion constants. The effort resulted in the IAU 2015 Resolution B3, passed at the IAU General Assembly by a large majority. The resolution recommends the use of nominal solar and planetary values, which are by definition exact and are expressed in SI units. These nominal values should be understood as conversion factors only, not as the true solar/planetary properties or current best estimates. Authors and journal editors are urged to join in using the standard values set forth by this resolution in future work and publications to help minimize further confusion.

The interior rotation of a sample of gamma Doradus stars from ensemble modelling of their gravity mode period spacings

Abstract: CONTEXT. Gamma Doradus stars (hereafter gamma Dor stars) are known to exhibit gravity- and/or gravito-intertial modes that probe the inner stellar region near the convective core boundary. The non-equidistant spacing of the pulsation periods is an observational signature of the stars' evolution and current internal structure and is heavily influenced by rotation. AIMS. We aim to constrain the near-core rotation rates for a sample of gamma Dor stars, for which we have detected period spacing patterns. METHODS. We combined the asymptotic period spacing with the traditional approximation of stellar pulsation to fit the observed period spacing patterns using chi-squared optimisation. The method was applied to the observed period spacing patterns of a sample of stars and used for ensemble modelling. RESULTS. For the majority of stars with an observed period spacing pattern we successfully determined the rotation rates and the asymptotic period spacing values, though the uncertainty margins on the latter were typically large. This also resulted directly in the identification of the modes corresponding with the detected pulsation frequencies, which for most stars were prograde dipole gravity and gravito-inertial modes. The majority of the observed retrograde modes were found to be Rossby modes. We further discuss the limitations of the method due to the neglect of the centrifugal force and the incomplete treatment of the Coriolis force. CONCLUSION. Despite its current limitations, the proposed methodology was successful to derive the rotation rates and to identify the modes from the observed period spacing patterns. It forms the first step towards detailed seismic modelling based on observed period spacing patterns of moderately to rapidly rotating gamma Dor stars.

Lupus disks with faint CO isotopologues: low gas/dust or large carbon depletion?

Abstract: An era has started in which gas and dust can be observed independently in protoplanetary disks, thanks to the recent surveys with ALMA. The first near-complete high-resolution disk survey in both dust and gas in a single star-forming region has been carried out in Lupus, finding surprisingly low gas/dust ratios. The goal of this work is to fully exploit CO isotopologues observations in Lupus, comparing them with physical-chemical model results, in order to obtain gas masses for a large number of disks. We have employed physical-chemical models to analyze continuum and CO isotopologues observations of Lupus disks, including isotope-selective processes and freeze-out. Employing also the ALMA 13CO-only detections, disk gas masses have been calculated for a total of 34 sources, expanding the sample of 10 disks studied by Ansdell et al. (2016), where also C18O was detected. We confirm that overall gas-masses are very low, often smaller than 1 $M_{\rm J}$, if volatile carbon is not depleted. Accordingly, global gas/dust ratios predominantly between 1 and 10. Low CO-based gas masses and gas/dust ratios may indicate rapid loss of gas, or alternatively chemical evolution, e.g. via sequestering of carbon from CO to more complex molecules, or carbon locked up in larger bodies. Current ALMA observations cannot distinguish between these two hypotheses. We have simulated both scenarios, but chemical model results do not allow us to rule out one of the two. Assuming that all Lupus disks have evolved mainly due to viscous processes over the past few Myr, the observed correlation between the current mass accretion rate and dust mass found by Manara et al. (2016) implies a constant gas-to-dust ratio, which is close to 100 based on the observed $M_{\rm disk}/\dot{M}_{\rm acc}$ ratio. This in turn points to a scenario in which carbon depletion is responsible for the low CO isotopologue line luminosities.

Physics Of Eclipsing Binaries. II. Towards the Increased Model Fidelity

Abstract: The precision of photometric and spectroscopic observations has been systematically improved in the last decade, mostly thanks to space-borne photometric missions and ground-based spectrographs dedicated to finding exoplanets. The field of eclipsing binary stars strongly benefited from this development. Eclipsing binaries serve as critical tools for determining fundamental stellar properties (masses, radii, temperatures and luminosities), yet the models are not capable of reproducing observed data well either because of the missing physics or because of insufficient precision. This led to a predicament where radiative and dynamical effects, insofar buried in noise, started showing up routinely in the data, but were not accounted for in the models. PHOEBE (PHysics Of Eclipsing BinariEs; this http URL) is an open source modeling code for computing theoretical light and radial velocity curves that addresses both problems by incorporating missing physics and by increasing the computational fidelity. In particular, we discuss triangulation as a superior surface discretization algorithm, meshing of rotating single stars, light time travel effect, advanced phase computation, volume conservation in eccentric orbits, and improved computation of local intensity across the stellar surfaces that includes photon-weighted mode, enhanced limb darkening treatment, better reflection treatment and Doppler boosting. Here we present the concepts on which PHOEBE is built on and proofs of concept that demonstrate the increased model fidelity.

Mass measurements in protoplanetary disks from hydrogen deuteride

Abstract: The total gas mass of a protoplanetary disk is a fundamental, but poorly determined, quantity. A new technique \citep{bergin+13} has been demonstrated to assess directly the bulk molecular gas reservoir of molecular hydrogen using the HD J=1-0 line at 112 $\mu$m. In this work we present a {\it Herschel} Space Observatory\footnote{Herschel is an ESA space observatory with science instruments provided by European-led Principal Investigator consortia and with important participation from NASA.} survey of six additional T Tauri disks in the HD line. Line emission is detected at $>$3$\sigma$ significance in two cases: DM Tau and GM Aur. For the other four disks, we establish upper limits to the line flux. Using detailed disk structure and ray tracing models, we calculate the temperature structure and dust mass from modeling the observed spectral energy distributions, and include the effect of UV gas heating to determine the amount of gas required to fit the HD line. The range of gas masses are 1.0-4.7$\times10^{-2}$ for DM Tau and 2.5-20.4$\times10^{-2}$ for GM Aur. These values are larger than those found using CO for GM Aur, while the CO-derived gas mass for DM Tau is consistent with the lower end of our mass range. This suggests a CO chemical depletion from the gas phase of up to a factor of five for DM Tau and up to two orders of magnitude for GM Aur. We discuss how future analysis can narrow the mass ranges further.

The impact of companions on stellar evolution

Abstract: Stellar astrophysicists are increasingly taking into account the effects of orbiting companions on stellar evolution. New discoveries, many thanks to systematic time-domain surveys, have underlined the role of binary star interactions in a range of astrophysical events, including some that were previously interpreted as due uniquely to single stellar evolution. Here, we review classical binary phenomena such as type Ia supernovae, and discuss new phenomena such as intermediate luminosity transients, gravitational wave-producing double black holes, or the interaction between stars and their planets. Finally, we examine the reassessment of well-known phenomena in light of interpretations that include both single and binary stars, for example supernovae of type Ib and Ic or luminous blue variables. At the same time we contextualise the new discoveries within the framework and nomenclature of the corpus of knowledge on binary stellar evolution. The last decade has heralded an era of revival in stellar astrophysics as the complexity of stellar observations is increasingly interpreted with an interplay of single and binary scenarios. The next decade, with the advent of massive projects such as the Large Synoptic Survey Telescope, the Square Kilometre Array, the James Webb Space Telescope and increasingly sophisticated computational methods, will see the birth of an expanded framework of stellar evolution that will have repercussions in many other areas of astrophysics such as galactic evolution and nucleosynthesis.

The MUSCLES Treasury Survey II: Intrinsic Lyman Alpha and Extreme Ultraviolet Spectra of K and M Dwarfs with Exoplanets

Abstract: The ultraviolet (UV) spectral energy distributions of low-mass (K- and M-type) stars play a critical role in the heating and chemistry of exoplanet atmospheres, but are not observationally well-constrained. Direct observations of the intrinsic flux of the Lyman alpha line (the dominant source of UV photons from low-mass stars) are challenging, as interstellar HI absorbs the entire line core for even the closest stars. To address the existing gap in empirical constraints on the UV flux of K and M dwarfs, the MUSCLES HST Treasury Survey has obtained UV observations of 11 nearby M and K dwarfs hosting exoplanets. This paper presents the Lyman alpha and extreme-UV spectral reconstructions for the MUSCLES targets. Most targets are optically inactive, but all exhibit significant UV activity. We use a Markov Chain Monte Carlo technique to correct the observed Lyman alpha profiles for interstellar absorption, and we employ empirical relations to compute the extreme-UV spectral energy distribution from the intrinsic Lyman alpha flux in ~100 A bins from 100-1170 A. The reconstructed Lyman alpha profiles have 300 km/s broad cores, while >1% of the total intrinsic Lyman alpha flux is measured in extended wings between 300 km/s to 1200 km/s. The Lyman alpha surface flux positively correlates with the MgII surface flux and negatively correlates with the stellar rotation period. Stars with larger Lyman alpha surface flux also tend to have larger surface flux in ions formed at higher temperatures, but these correlations remain statistically insignificant in our sample of 11 stars. We also present HI column density measurements for 10 new sightlines through the local interstellar medium.

Magnetic Inhibition of Convection and the Fundamental Properties of Low-Mass Stars. III. A Consistent 10 Myr Age for the Upper Scorpius OB Association

Abstract: When determining absolute ages of identifiably young stellar populations, results strongly depend on which stars are studied. Cooler (K, M) stars typically yield ages that are systematically younger than warmer (A, F, G) stars by a factor of two. I explore the possibility that these age discrepancies are the result of magnetic inhibition of convection in cool young stars by using magnetic stellar evolution isochrones to determine the age of the Upper Scorpius subgroup of the Scorpius-Centaurus OB Association. A median age of 10 Myr consistent across spectral types A through M is found, except for a subset of F-type stars that appear significantly older. Agreement is shown for ages derived from the Hertzsprung-Russell diagram and from the empirical mass-radius relationship defined by eclipsing multiple-star systems. Surface magnetic field strengths required to produce agreement are of order 2.5 kG and are predicted from a priori estimates of equipartition values. A region in the HR diagram is identified that plausibly connects stars whose structures are weakly influenced by the presence of magnetic fields with those whose structures are strongly influenced by magnetic fields. The models suggest this region is characterized by stars with rapidly thinning outer convective envelopes where the radiative core mass is greater than 75% of the total stellar mass. Furthermore, depletion of lithium predicted from magnetic models appears in better agreement with observed lithium equivalent widths than predictions from non-magnetic models. These results suggest that magnetic inhibition of convection plays an important role in the early evolution of low-mass stars and that it may be responsible for noted age discrepancies in young stellar populations.

Spectral Properties of Cool Stars: Extended Abundance Analysis of 1626 Planet Search Stars

Abstract: We present a catalog of uniformly determined stellar properties and abundances for 1626 F, G, and K stars using an automated spectral synthesis modeling procedure. All stars were observed using the HIRES spectrograph at Keck Observatory. Our procedure used a single line list to fit model spectra to observations of all stars to determine effective temperature, surface gravity, metallicity, projected rotational velocity, and the abundances of 15 elements (C, N, O, Na, Mg, Al, Si, Ca, Ti, V, Cr, Mn, Fe, Ni, & Y). Sixty percent of the sample had Hipparcos parallaxes and V-band photometry which we combined with the spectroscopic results to obtain mass, radius, and luminosity. Additionally, we used the luminosity, effective temperature, metallicity and alpha-element enhancement to interpolate in the Yonsei-Yale isochrones to derive mass, radius, gravity, and age ranges for those stars. Finally, we determined new relations between effective temperature and macroturbulence for dwarfs and subgiants. Our analysis achieved precisions of 25 K in Teff , 0.01 dex in [M/H], 0.028 dex for log g and 0.5 km/s in v sin ibased on multiple observations of the same stars. The abundance results were similarly precise, between 0.01 and - 0.04 dex, though trends with respect to Teff remained for which we derived empirical corrections. The trends, though small, were much larger than our uncertainties and are shared with published abundances. We show that changing our model atmosphere grid accounts for most of the trend in [M/H] between 5000 K and 5500 K indicating a possible problem with the atmosphere models or opacities.

Determining protoplanetary disk gas masses from CO isotopologues line observations

Abstract: Despite intensive studies of protoplanetary disks, there is still no reliable way to determine their total mass and their surface density distribution, quantities that are crucial for describing both the structure and the evolution of disks up to the formation of planets. The goal of this work is to use less abundant CO isotopologues, whose detection is routine for ALMA, to infer the gas mass of disks. Isotope-selective effects need to be taken into account in the analysis, because they can significantly modify CO isotopologues line intensities. CO isotope-selective photodissociation has been implemented in the physical-chemical code DALI and 800 disk models have been run for a range of disk and stellar parameters. Dust and gas temperature structures have been computed self-consistently, together with a chemical calculation of the main species. Both disk structure and stellar parameters have been investigated. Total fluxes have been ray-traced for different CO isotopologues and for various transitions for different inclinations. A combination of 13CO and C18O total intensities allows inference of the total disk mass, although with non-negligible uncertainties. These can be overcome by employing spatially resolved observations, i.e. the disk's radial extent and inclination. Comparison with parametric models shows differences at the factor of a few level, especially for extremely low and high disk masses. Finally, total line intensities for different CO isotopologue and for various low-J transitions are provided and are fitted to simple formulae. The effects of a lower gas-phase carbon abundance and different gas/dust ratios are investigated as well, and comparison with other tracers is made. Disk masses can be determined within a factor of a few by comparing CO isotopologue lines observations with the simulated line fluxes, modulo the uncertainties in the volatile elemental abundances.

Observables Processing for the Helioseismic and Magnetic Imager Instrument on the Solar Dynamics Observatory

Abstract: NASA's Solar Dynamics Observatory (SDO) was launched 11 February 2010 with 3 instruments on board, including the Helioseismic and Magnetic Imager (HMI). Since beginning normal operations on 1 May 2010, HMI has observed the Sun's entire visible disk almost continuously. HMI collects sequences of polarized filtergrams taken at a fixed cadence with two 4096 x 4096 cameras from which are computed arcsecond-resolution maps of photospheric observables: the line-of-sight (LoS) velocity and magnetic field, continuum intensity, line width, line depth, and the Stokes polarization parameters, I Q U V, at 6 wavelengths. Two processing pipelines implemented at the SDO Joint Science Operations Center (JSOC) at Stanford University compute observables from calibrated Level-1 filtergrams. One generates LoS quantities every 45s, and the other, primarily for the vector magnetic field, computes averages on a 720s cadence. Corrections are made for static and temporally changing CCD characteristics, bad pixels, image alignment and distortion, polarization irregularities, filter-element uncertainty and non-uniformity, as well as Sun-spacecraft velocity. This report explains issues affecting the resulting physical quantities, describes the impact of regular updates to the instrument calibration, and shows how the computations are optimized for actual HMI observations. During the 5 years of the SDO prime mission, regular calibration sequences have been used to regularly improve and update the instrument calibration and to monitor instrument changes. The instrument more than satisfies the original specifications for data quality and continuity. The procedures described here still have significant room for improvement. The most significant remaining systematic errors are associated with the spacecraft orbital velocity.

The MUSCLES Treasury Survey III: X-ray to Infrared Spectra of 11 M and K Stars Hosting Planets

Abstract: We present a catalog of panchromatic spectral energy distributions (SEDs) for 7 M and 4 K dwarf stars that span X-ray to infrared wavelengths (5 A - 5.5 {\mu}m). These SEDs are composites of Chandra or XMM-Newton data from 5 - ~50 A, a plasma emission model from ~50 - 100 A, broadband empirical estimates from 100 - 1170 A, HST data from 1170 - 5700 A, including a reconstruction of stellar Ly{\alpha} emission at 1215.67 A, and a PHOENIX model spectrum from 5700 - 55000 A. Using these SEDs, we computed the photodissociation rates of several molecules prevalent in planetary atmospheres when exposed to each star's unattenuated flux ("unshielded" photodissociation rates) and found that rates differ among stars by over an order of magnitude for most molecules. In general, the same spectral regions drive unshielded photodissociations both for the minimally and maximally FUV active stars. However, for ozone visible flux drives dissociation for the M stars whereas NUV flux drives dissociation for the K stars. We also searched for an FUV continuum in the assembled SEDs and detected it in 5/11 stars, where it contributes around 10% of the flux in the range spanned by the continuum bands. An ultraviolet continuum shape is resolved for the star {\epsilon} Eri that shows an edge likely attributable to Si II recombination. The 11 SEDs presented in this paper, available online through the Mikulski Archive for Space Telescopes, will be valuable for vetting stellar upper-atmosphere emission models and simulating photochemistry in exoplanet atmospheres.

The Radial Distribution of H2 and CO in TW Hya as Revealed by Resolved ALMA Observations of CO Isotopologues

Abstract: CO is widely used as a tracer of molecular gas. However, there is now mounting evidence that gas phase carbon is depleted in the disk around TW Hya. Previous efforts to quantify this depletion have been hampered by uncertainties regarding the radial thermal structure in the disk. Here we present resolved ALMA observations of 13CO 3-2, C18O 3-2, 13CO 6-5, and C18O 6-5 emission in TW Hya, which allow us to derive radial gas temperature and gas surface density profiles, as well as map the CO abundance as a function of radius. These observations provide a measurement of the surface CO snowline at ~30 AU and show evidence for an outer ring of CO emission centered at 53 AU, a feature previously seen only in less abundant species. Further, the derived CO gas temperature profile constrains the freeze-out temperature of CO in the warm molecular layer to < 21 K. Combined with the previous detection of HD 1-0, these data constrain the surface density of the warm H2 gas in the inner ~30 AU. We find that CO is depleted by two orders of magnitude from R=10-60 AU, with the small amount of CO returning to the gas phase inside the surface CO snowline insufficient to explain the overall depletion. Finally, this new data is used in conjunction with previous modeling of the TW Hya disk to constrain the midplane CO snowline to 17-23 AU.

Dust Evolution and the Formation of Planetesimals

Abstract: The solid content of circumstellar disks is inherited from the interstellar medium: dust particles of at most a micrometer in size. Protoplanetary disks are the environment where these dust grains need to grow at least 13 orders of magnitude in size. Our understanding of this growth process is far from complete, with different physics seemingly posing obstacles to this growth at various stages. Yet, the ubiquity of planets in our galaxy suggests that planet formation is a robust mechanism. This chapter focuses on the earliest stages of planet formation, the growth of small dust grains towards the gravitationally bound "planetesimals", the building blocks of planets. We will introduce some of the key physics involved in the growth processes and discuss how they are expected to shape the global behavior of the solid content of disks. We will consider possible pathways towards the formation of larger bodies and conclude by reviewing some of the recent observational advances in the field.

A Comparison of Flare Forecasting Methods, I: Results from the "All-Clear" Workshop

Abstract: Solar flares produce radiation which can have an almost immediate effect on the near-Earth environment, making it crucial to forecast flares in order to mitigate their negative effects. The number of published approaches to flare forecasting using photospheric magnetic field observations has proliferated, with varying claims about how well each works. Because of the different analysis techniques and data sets used, it is essentially impossible to compare the results from the literature. This problem is exacerbated by the low event rates of large solar flares. The challenges of forecasting rare events have long been recognized in the meteorology community, but have yet to be fully acknowledged by the space weather community. During the interagency workshop on "all clear" forecasts held in Boulder, CO in 2009, the performance of a number of existing algorithms was compared on common data sets, specifically line-of-sight magnetic field and continuum intensity images from MDI, with consistent definitions of what constitutes an event. We demonstrate the importance of making such systematic comparisons, and of using standard verification statistics to determine what constitutes a good prediction scheme. When a comparison was made in this fashion, no one method clearly outperformed all others, which may in part be due to the strong correlations among the parameters used by different methods to characterize an active region. For M-class flares and above, the set of methods tends towards a weakly positive skill score (as measured with several distinct metrics), with no participating method proving substantially better than climatological forecasts.

A cavity and further radial substructures in the disk around HD~97048

Abstract: Context: Gaps, cavities and rings in circumstellar disks are signposts of disk evolution and planet-disk interactions. We follow the recent suggestion that Herbig Ae/Be disks with a flared disk harbour a cavity, and investigate the disk around HD~97048. Aims: We aim to resolve the 34$\pm$ 4 au central cavity predicted by Maaskant et al. (2013) and to investigate the structure of the disk. Methods: We image the disk around HD~97048 using ALMA at 0.85~mm and 2.94~mm, and ATCA (multiple frequencies) observations. Our observations also include the 12CO J=1-0, 12CO J=3-2 and HCO+ J=4-3 emission lines. Results: A central cavity in the disk around HD~97048 is resolved with a 40-46 au radius. Additional radial structure present in the surface brightness profile can be accounted for either by an opacity gap at ~90 au or by an extra emitting ring at ~150 au. The continuum emission tracing the dust in the disk is detected out to 355 au. The 12CO J=3-2 disk is detected 2.4 times farther out. The 12CO emission can be traced down to $\approx$ 10 au scales. Non-Keplerian kinematics are detected inside the cavity via the HCO+ J=4-3 velocity map. The mm spectral index measured from ATCA observations suggests that grain growth has occurred in the HD~97048 disk. Finally, we resolve a highly inclined disk out to 150 au around the nearby 0.5~$M_{\odot}$ binary ISO-ChaI 126. Conclusions: The data presented here reveal a cavity in the disk of HD 97048, and prominent radial structure in the surface brightness. The cavity size varies for different continuum frequencies and gas tracers. The gas inside the cavity follows non-Keplerian kinematics seen in HCO+ emission. The variable cavity size along with the kinematical signature suggests the presence of a substellar companion or massive planet inside the cavity.