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

Binary Population and Spectral Synthesis Version 2.1: construction, observational verification and new results

Abstract: The Binary Population and Spectral Synthesis (BPASS) suite of binary stellar evolution models and synthetic stellar populations provides a framework for the physically motivated analysis of both the integrated light from distant stellar populations and the detailed properties of those nearby. We present a new version 2.1 data release of these models, detailing the methodology by which BPASS incorporates binary mass transfer and its effect on stellar evolution pathways, as well as the construction of simple stellar populations. We demonstrate key tests of the latest BPASS model suite demonstrating its ability to reproduce the colours and derived properties of resolved stellar populations, including well- constrained eclipsing binaries. We consider observational constraints on the ratio of massive star types and the distribution of stellar remnant masses. We describe the identification of supernova progenitors in our models, and demonstrate a good agreement to the properties of observed progenitors. We also test our models against photometric and spectroscopic observations of unresolved stellar populations, both in the local and distant Universe, finding that binary models provide a self-consistent explanation for observed galaxy properties across a broad redshift range. Finally, we carefully describe the limitations of our models, and areas where we expect to see significant improvement in future versions.

Modules for Experiments in Stellar Astrophysics (MESA): Convective Boundaries, Element Diffusion, and Massive Star Explosions

Abstract: We update the capabilities of the software instrument Modules for Experiments in Stellar Astrophysics (MESA) and enhance its ease of use and availability. Our new approach to locating convective boundaries is consistent with the physics of convection, and yields reliable values of the convective core mass during both hydrogen and helium burning phases. Stars with $M<8\,{\rm M_\odot}$ become white dwarfs and cool to the point where the electrons are degenerate and the ions are strongly coupled, a realm now available to study with MESA due to improved treatments of element diffusion, latent heat release, and blending of equations of state. Studies of the final fates of massive stars are extended in MESA by our addition of an approximate Riemann solver that captures shocks and conserves energy to high accuracy during dynamic epochs. We also introduce a 1D capability for modeling the effects of Rayleigh-Taylor instabilities that, in combination with the coupling to a public version of the STELLA radiation transfer instrument, creates new avenues for exploring Type II supernovae properties. These capabilities are exhibited with exploratory models of pair-instability supernova, pulsational pair-instability supernova, and the formation of stellar mass black holes. The applicability of MESA is now widened by the capability of importing multi-dimensional hydrodynamic models into MESA. We close by introducing software modules for handling floating point exceptions and stellar model optimization, and four new software tools -- MESAWeb, MESA-Docker, pyMESA, and this http URL -- to enhance MESA's education and research impact.

ALMA unveils rings and gaps in the protoplanetary system HD 169142: signatures of two giant protoplanets

Abstract: The protoplanetary system HD 169142 is one of the few cases where a potential candidate protoplanet has been recently detected via direct imaging. To study the interaction between the protoplanet and the disk itself observations of the gas and dust surface density structure are needed. This paper reports new ALMA observations of the dust continuum at 1.3\,mm, $^{12}$CO, $^{13}$CO and C$^{18}$O $J=2-1$ emission from the system HD 169142 at angular resolution of $\sim 0".18 - 0".28$ ($\sim 20\,$au$ - 33\,$au). The dust continuum emission reveals a double-ring structure with an inner ring between $0".17-0".28$ ($\sim 20 - 35\,$au) and an outer ring between $0".48-0".64$ ($\sim 56 - 83\,$au). The size and position of the inner ring is in good agreement with previous polarimetric observations in the near-infrared and is consistent with dust trapping by a massive planet. No dust emission is detected inside the inner dust cavity ($R \lesssim 20\,$au) or within the dust gap ($\sim 35 - 56\,$au). In contrast, the channel maps of the $J=2-1$ line of the three CO isotopologues reveal the presence of gas inside the dust cavity and dust gap. The gaseous disk is also much larger than the compact dust emission extending to $\sim 1'.5$ ($\sim 180\,$au) in radius. This difference and the sharp drop of the continuum emission at large radii point to radial drift of large dust grains ($>$ \micron-size). Using the thermo-chemical disk code \textsc{dali}, the continuum and the CO isotopologues emission are modelled to quantitatively measure the gas and dust surface densities. The resulting gas surface density is reduced by a factor of $\sim 30-40$ inward of the dust gap. The gas and dust distribution hint at the presence of multiple planets shaping the disk structure via dynamical clearing (dust cavity and gap) and dust trapping (double ring dust distribution).

Photodissociation and photoionisation of atoms and molecules of astrophysical interest

Abstract: A new collection of photodissociation and photoionisation cross sections for 102 atoms and molecules of astrochemical interest has been assembled, along with a brief review of the basic processes involved. These have been used to calculate dissociation and ionisation rates, with uncertainties, in a standard ultraviolet interstellar radiation field (ISRF) and wavelength-dependent radiation fields. The new ISRF rates generally agree within 30% with our previous compilations, with a few notable exceptions. The reduction of rates in shielded regions was calculated as a function of dust, molecular and atomic hydrogen, atomic C, and self-shielding column densities. The relative importance of shielding types depends on the species in question and the dust optical properties. The new data are publicly available from the Leiden photodissociation and ionisation database. Sensitivity of rates to variation of temperature and isotope, and cross section uncertainties, are tested. Tests were conducted with an interstellar-cloud chemical model, and find general agreement (within a factor of two) with the previous iteration of the Leiden database for the ISRF, and order-of-magnitude variations assuming various kinds of stellar radiation. The newly parameterised dust-shielding factors makes a factor-of-two difference to many atomic and molecular abundances relative to parameters currently in the UDfA and KIDA astrochemical reaction databases. The newly-calculated cosmic-ray induced photodissociation and ionisation rates differ from current standard values up to a factor of 5. Under high temperature and cosmic-ray-flux conditions the new rates alter the equilibrium abundances of abundant dark cloud abundances by up to a factor of two. The partial cross sections for H2O and NH3 photodissociation forming OH, O, NH2 and NH are also evaluated and lead to radiation-field-dependent branching ratios.

The CARMENES search for exoplanets around M dwarfs: High-resolution optical and near-infrared spectroscopy of 324 survey stars

Abstract: The CARMENES radial velocity (RV) survey is observing 324 M dwarfs to search for any orbiting planets. In this paper, we present the survey sample by publishing one CARMENES spectrum for each M dwarf. These spectra cover the wavelength range 520--1710nm at a resolution of at least $R > 80,000$, and we measure its RV, H$\alpha$ emission, and projected rotation velocity. We present an atlas of high-resolution M-dwarf spectra and compare the spectra to atmospheric models. To quantify the RV precision that can be achieved in low-mass stars over the CARMENES wavelength range, we analyze our empirical information on the RV precision from more than 6500 observations. We compare our high-resolution M-dwarf spectra to atmospheric models where we determine the spectroscopic RV information content, $Q$, and signal-to-noise ratio. We find that for all M-type dwarfs, the highest RV precision can be reached in the wavelength range 700--900nm. Observations at longer wavelengths are equally precise only at the very latest spectral types (M8 and M9). We demonstrate that in this spectroscopic range, the large amount of absorption features compensates for the intrinsic faintness of an M7 star. To reach an RV precision of 1ms$^{-1}$ in very low mass M dwarfs at longer wavelengths likely requires the use of a 10m class telescope. For spectral types M6 and earlier, the combination of a red visual and a near-infrared spectrograph is ideal to search for low-mass planets and to distinguish between planets and stellar variability. At a 4m class telescope, an instrument like CARMENES has the potential to push the RV precision well below the typical jitter level of 3-4ms$^{-1}$.

A new generation of PARSEC-COLIBRI stellar isochrones including the TP-AGB phase

Abstract: We introduce a new generation of PARSEC-COLIBRI stellar isochrones that include a detailed treatment of the thermally-pulsing asymptotic giant branch (TP-AGB) phase, and covering a wide range of initial metallicities (0.0001

Different dust and gas radial extents in protoplanetary disks: consistent models of grain growth and CO emission

Abstract: ALMA observations of protoplanetary disks confirm earlier indications that there is a clear difference between the dust and gas radial extents The origin of this difference is still debated, with both radial drift of the dust and optical depth effects suggested in the literature In this work, the feedback of realistic dust particle distributions onto the gas chemistry and molecular emissivity is investigated, with a particular focus on CO isotopologues The radial dust grain size distribution is determined using dust evolution models that include growth, fragmentation and radial drift A new version of the code DALI is used to take into account how dust surface area and density influence the disk thermal structure, molecular abundances and excitation The difference of dust and gas radial sizes is largely due to differences in the optical depth of CO lines and millimeter continuum, without the need to invoke radial drift The effect of radial drift is primarily visible in the sharp outer edge of the continuum intensity profile The gas outer radius probed by $^{12}$CO emission can easily differ by a factor of $\sim 2$ between the models for a turbulent $\alpha$ ranging between typical values Grain growth and settling concur in thermally decoupling the gas and dust components, due to the low collision rate with large grains As a result, the gas can be much colder than the dust at intermediate heights, reducing the CO excitation and emission, especially for low turbulence values Also, due to disk mid-plane shadowing, a second CO thermal desorption (rather than photodesorption) front can occur in the warmer outer mid-plane disk The models are compared to ALMA observations of HD 163296 as a test case In order to reproduce the observed CO snowline of the system, a binding energy for CO typical of ice mixtures needs to be used rather than the lower pure CO value

Multiple Stellar Populations in Globular Clusters

Abstract: Globular Clusters (GCs) exhibit star-to-star variations in specific elements (e.g., He, C, N, O, Na, Al) that bare the hallmark of high temperature H burning. These abundance variations can be observed spectroscopically and also photometrically, with the appropriate choice of filters, due to the changing of spectral features within the band pass. This phenomenon is observed in nearly all of the ancient GCs, although, to date, has not been found in any massive cluster younger than 2~Gyr. Many scenarios have been suggested to explain this phenomenon, with most invoking multiple epochs of star-formation within the cluster, however all have failed to reproduce various key observations, in particular when a global view of the GC population is taken. We review the state of current observations, and outline the successes and failures of each of the main proposed models. The traditional idea of using the stellar ejecta from a 1st generation of stars to form a 2nd generation of stars, while conceptually straight forward, has failed to reproduce an increasing number of observational constraints. We conclude that the puzzle of multiple populations remains unsolved, hence alternative theories are needed.

On the generation of solar spicules and Alfv\'enic waves

Abstract: In the lower solar atmosphere, the chromosphere is permeated by jets known as spicules, in which plasma is propelled at speeds of 50 to 150 kilometers per second into the corona. The origin of the spicules is poorly understood, although they are expected to play a role in heating the million-degree corona and are associated with Alfv\'enic waves that help drive the solar wind. We compare magnetohydrodynamic simulations of spicules with observations from the Interface Region Imaging Spectrograph and the Swedish 1-m Solar Telescope. Spicules are shown to occur when magnetic tension is amplified and transported upward through interactions between ions and neutrals or ambipolar diffusion. The tension is impulsively released to drive flows, heat plasma (through ambipolar diffusion), and generate Alfv\'enic waves.

White Dwarf Rotation as a Function of Mass and a Dichotomy of Mode Linewidths: Kepler Observations of 27 Pulsating DA White Dwarfs Through K2 Campaign 8

Abstract: We present photometry and spectroscopy for 27 pulsating hydrogen-atmosphere white dwarfs (DAVs, a.k.a. ZZ Ceti stars) observed by the Kepler space telescope up to K2 Campaign 8, an extensive compilation of observations with unprecedented duration (>75 days) and duty cycle (>90%). The space-based photometry reveals pulsation properties previously inaccessible to ground-based observations. We observe a sharp dichotomy in oscillation mode linewidths at roughly 800 s, such that white dwarf pulsations with periods exceeding 800 s have substantially broader mode linewidths, more reminiscent of a damped harmonic oscillator than a heat-driven pulsator. Extended Kepler coverage also permits extensive mode identification: We identify the spherical degree of 61 out of 154 unique radial orders, providing direct constraints of the rotation period for 20 of these 27 DAVs, more than doubling the number of white dwarfs with rotation periods determined via asteroseismology. We also obtain spectroscopy from 4m-class telescopes for all DAVs with Kepler photometry. Using these homogeneously analyzed spectra we estimate the overall mass of all 27 DAVs, which allows us to measure white dwarf rotation as a function of mass, constraining the endpoints of angular momentum in low- and intermediate-mass stars. We find that 0.51-to-0.73-solar-mass white dwarfs, which evolved from 1.7-to-3.0-solar-mass ZAMS progenitors, have a mean rotation period of 35 hr with a standard deviation of 28 hr, with notable exceptions for higher-mass white dwarfs. Finally, we announce an online repository for our Kepler data and follow-up spectroscopy, which we collect at this http URL.

Precision Stellar Characterization of FGKM Stars using an Empirical Spectral Library

Abstract: Classification of stars, by comparing their optical spectra to a few dozen spectral standards, has been a workhorse of observational astronomy for more than a century. Here, we extend this technique by compiling a library of optical spectra of 404 touchstone stars observed with Keck/HIRES by the California Planet Search. The spectra are high-resolution (R~60000), high signal-to-noise (SNR~150/pixel), and registered onto a common wavelength scale. The library stars have properties derived from interferometry, asteroseismology, LTE spectral synthesis, and spectrophotometry. To address a lack of well-characterized late K-dwarfs in the literature, we measure stellar radii and temperatures for 23 nearby K-dwarfs, using SED modeling and Gaia parallaxes. This library represents a uniform dataset spanning the spectral types ~M5--F1 (Teff ~ 3000-7000K, Rstar ~ 0.1-1.6 Rsun). We also present "Empirical SpecMatch" (SpecMatch-Emp), a tool for parameterizing unknown spectra by comparing them against our spectral library. For FGKM stars, SpecMatch-Emp achieves accuracies of 100 K in effective temperature (Teff), 15% in stellar radius (Rstar), and 0.09 dex in metallicity [Fe/H]. Because the code relies on empirical spectra it performs particularly well for stars ~K4 and later which are challenging to model with existing spectral synthesizers, reaching accuracies of 70 K in Teff, 10% in Rstar, and 0.12 dex in [Fe/H]. We also validate the performance of SpecMatch-Emp, finding it to be robust at lower spectral resolution and SNR, enabling the characterization of faint late-type stars. Both the library and stellar characterization code are publicly available.

The binarity of the local white dwarf population

Abstract: White dwarfs (WDs) are powerful tools to study the evolutionary history of stars and binaries in the Galaxy. But do we understand their multiplicity from a theoretical point of view? This can be tested by a comparison with the sample of WDs within 20 pc, which is minimally affected by selection biases. From the literature, we compile the available information of the local WD sample with a particular emphasis on their multiplicity, and compare this to synthetic models of WD formation in single stars and binaries. As part of our population synthesis approach, we also study the effect of different assumptions concerning the star formation history, binary evolution, and the initial distributions of binary parameters. We find that the observed space densities of single and binary WDs are well reproduced by the models. The space densities of the most common WD systems (single WDs and unresolved WD-MS binaries) are consistent within a factor two with the observed value. We find a discrepancy only for the space density of resolved double WDs. We exclude that observational selection effects, fast stellar winds, or dynamical interactions with other objects in the Milky Way explain this discrepancy. We find that either the initial mass ratio distribution in the Solar neighbourhood is biased towards low mass-ratios, or more than ten resolved DWDs have been missed observationally in the 20 pc sample. Furthermore, we show that the low binary fraction of WD systems (~25%) compared to Solartype MS-MS binaries (~50%) is consistent with theory, and is mainly caused by mergers in binary systems, and to a lesser degree by WDs hiding in the glare of their companion stars. Lastly, Gaia will dramatically increase the size of the volume-limited WD sample, detecting the coolest and oldest WDs out to 50 pc. We provide a detailed estimate of the number of single and binary WDs in the Gaia sample.

An ALMA survey of DCN/H$^{13}$CN and DCO$^+$/H$^{13}$CO$^+$ in protoplanetary disks

Abstract: The deuterium enrichment of molecules is sensitive to their formation environment. Constraining patterns of deuterium chemistry in protoplanetary disks is therefore useful for probing how material is inherited or reprocessed throughout the stages of star and planet formation. We present ALMA observations at $\sim0.6"$ resolution of DCO$^+$, H$^{13}$CO$^+$, DCN, and H$^{13}$CN in the full disks around T Tauri stars AS 209 and IM Lup, the transition disks around T Tauri stars V4046 Sgr and LkCa 15, and the full disks around Herbig Ae stars MWC 480 and HD 163296. We also present ALMA observations of HCN in the IM Lup disk. DCN, DCO$^+$, and H$^{13}$CO$^+$ are detected in all disks, and H$^{13}$CN in all but the IM Lup disk. We find efficient deuterium fractionation for the sample, with estimates of disk-averaged DCO$^+$/HCO$^+$ and DCN/HCN abundance ratios ranging from $\sim0.02-0.06$ and $\sim0.005-0.08$, respectively, which is comparable to values reported for other ISM environments. The relative distributions of DCN and DCO$^+$ vary between disks, suggesting that multiple formation pathways may be needed to explain the diverse emission morphologies. In addition, gaps and rings observed in both H$^{13}$CO$^+$ and DCO$^+$ emission provide new evidence that DCO$^+$ bears a complex relationship with the location of the midplane CO snowline.

WIYN Open Cluster Study. LXXV. Testing the Metallicity Dependence of Stellar Lithium Depletion Using Hyades-Aged Clusters. 1. Hyades & Praesepe

Abstract: WIYN/Hydra spectroscopy (at R~15,000) of the moderately metal-rich Praesepe and Hyades open clusters was used to study their main sequence (MS) iron ([Fe/H]) and lithium (A(Li)) abundances. Self-consistent [Fe/H] and Li analyses of these clusters of consistent age, which we re-evaluate, confirms they have consistent [Fe/H] and provides a foundation to investigate the poorly understood G-dwarf and F-dwarf Li-depletions. Neither phenomenon agrees with standard stellar evolution theory, but possible explanations abound. We supplement our A(Li) with previously published results placed on a uniform abundance scale. This creates the largest self-consistently analyzed sample of A(Li) in both the Hyades (90) and Praesepe (110). For each star, high-precision UBVRI photometry was used to determine a ten color-based Teff and then to test for photometric peculiarities indicated by a large sigma-Teff (> 75 K). The stars with large sigma-Teff were predominantly found to be binaries or stars with peculiar (apparent) A(Li). When considering only proper-motion members that have low sigma-Teff and are also photometrically consistent with the cluster MS fiducial, each cluster has a more tightly defined Li morphology than previously observed and the two clusters' A(Li) are indistinguishable. This suggests that clusters of consistent age and metallicity may have consistent Li-depletion trends across a broad range of Teff; no additional major parameters are required, at least for these two clusters. We propose that the combined Hyades and Praesepe data offer more rigorous constraints than does either cluster alone, and we discuss newly-revealed features of the combined Li-Teff trend.

Gaia Data Release 1. Open cluster astrometry: performance, limitations, and future prospects

Abstract: Context. The first Gaia Data Release contains the Tycho-Gaia Astrometric Solution (TGAS). This is a subset of about 2 million stars for which, besides the position and photometry, the proper motion and parallax are calculated using Hipparcos and Tycho-2 positions in 1991.25 as prior information. Aims. We investigate the scientific potential and limitations of the TGAS component by means of the astrometric data for open clusters. Methods. Mean cluster parallax and proper motion values are derived taking into account the error correlations within the astrometric solutions for individual stars, an estimate of the internal velocity dispersion in the cluster, and, where relevant, the effects of the depth of the cluster along the line of sight. Internal consistency of the TGAS data is assessed. Results. Values given for standard uncertainties are still inaccurate and may lead to unrealistic unit-weight standard deviations of least squares solutions for cluster parameters. Reconstructed mean cluster parallax and proper motion values are generally in very good agreement with earlier Hipparcos-based determination, although the Gaia mean parallax for the Pleiades is a significant exception. We have no current explanation for that discrepancy. Most clusters are observed to extend to nearly 15 pc from the cluster centre, and it will be up to future Gaia releases to establish whether those potential cluster-member stars are still dynamically bound to the clusters. Conclusions. The Gaia DR1 provides the means to examine open clusters far beyond their more easily visible cores, and can provide membership assessments based on proper motions and parallaxes. A combined HR diagram shows the same features as observed before using the Hipparcos data, with clearly increased luminosities for older A and F dwarfs.

An extraordinary outburst in the massive protostellar system NGC6334I-MM1: quadrupling of the millimeter continuum

Abstract: Based on sub-arcsecond Atacama Large Millimeter/submillimeter Array (ALMA) and Submillimeter Array (SMA) 1.3 mm continuum images of the massive protocluster NGC 6334I obtained in 2015 and 2008, we find that the dust emission from MM1 has increased by a factor of 4.0$\pm$0.3 during the intervening years, and undergone a significant change in morphology. The continuum emission from the other cluster members (MM2, MM4 and the UCHII region MM3 = NGC 6334F) has remained constant. Long term single-dish maser monitoring at HartRAO finds that multiple maser species toward NGC 6334I flared beginning in early 2015, a few months before our ALMA observation, and some persist in that state. New ALMA images obtained in 2016 July-August at 1.1 and 0.87 mm confirm the changes with respect to SMA 0.87 mm images from 2008, and indicate that the (sub)millimeter flaring has continued for at least a year. The excess continuum emission, centered on the hypercompact HII region MM1B, is extended and elongated ($1.6" \times 1.0" \approx 2100 \times 1300$ au) with multiple peaks, suggestive of general heating of the surrounding subcomponents of MM1, some of which may trace clumps in a fragmented disk rather than separate protostars. In either case, these remarkable increases in maser and dust emission provide direct observational evidence of a sudden accretion event in the growth of a massive protostar yielding a sustained luminosity surge by a factor of $70\pm20$, analogous to the largest events in simulations by Meyer et al. (2017). This target provides an excellent opportunity to assess the impact of such a rare event on a protocluster over many years.

JHelioviewer - Time-dependent 3D visualisation of solar and heliospheric data

Abstract: Context. Solar observatories are providing the world-wide community with a wealth of data, covering large time ranges, multiple viewpoints, and returning large amounts of data. In particular, the large volume of SDO data presents challenges: it is available only from a few repositories, and full-disk, full-cadence data for reasonable durations of scientific interest are difficult to download practically due to their size and download data rates available to most users. From a scientist's perspective this poses three problems: accessing, browsing and finding interesting data as efficiently as possible. Aims. To address these challenges, we have developed JHelioviewer, a visualisation tool for solar data based on the JPEG2000 compression standard and part of the open source ESA/NASA Helioviewer Project. Since the first release of JHelioviewer, the scientific functionality of the software has been extended significantly, and the objective of this paper is to highlight these improvements. Methods. The JPEG2000 standard offers useful new features that facilitate the dissemination and analysis of high-resolution image data and offers a solution to the challenge of efficiently browsing petabyte-scale image archives. The JHelioviewer software is open source, platform independent and extendable via a plug-in architecture. Results. With JHelioviewer, users can visualise the Sun for any time period between September 1991 and today. They can perform basic image processing in real time, track features on the Sun and interactively overlay magnetic field extrapolations. The software integrates solar event data and a time line display. As a first step towards supporting science planning of the upcoming Solar Orbiter mission, JHelioviewer offers a virtual camera model that enables users to set the vantage point to the location of a spacecraft or celestial body at any given time.

Twelve years of spectroscopic monitoring in the Galactic Center: the closest look at S-stars near the black hole

Abstract: We study the young S-stars within a distance of 0.04 pc from the supermassive black hole in the center of our Galaxy. Given how inhospitable the region is for star formation, their presence is more puzzling the younger we estimate their ages. In this study, we analyse the result of 12 years of high resolution spectroscopy within the central arcsecond of the Galactic Center (GC). By co-adding between 55 and 105 hours of spectra we have obtained high signal to noise H- and K-band spectra of eight stars orbiting the central supermassive black hole. Using deep H-band spectra, we show that these stars must be high surface gravity (dwarf) stars. We compare these deep spectra to detailed model atmospheres and stellar evolution models to infer the stellar parameters. Our analysis reveals an effective temperature of 21000-28500 K, a rotational velocity of 60-170 km/s, and a surface gravity of 4.1-4.2. These parameters imply a spectral type of B0-B3V for these stars. The inferred masses lie within 8-14 Msun. We derive an age of 6.6^{+3.4}{-4.7} Myr for the star S2, which is compatible with the age of the clockwise rotating young stellar disk in the GC. We estimate the age of all other studied S-stars to be less than 15 Myr, which are compatible with the age of S2 within the uncertainties. The relatively low ages for these S-stars favor a scenario in which the stars formed in a local disk rather than the field-binary-disruption scenario throughout a longer period of time.

Nucleosynthesis in the Innermost Ejecta of Neutrino-Drive Supernova Explosions in Two Dimensions

Abstract: We examine the nucleosynthesis in the innermost, neutrino-processed ejecta (a few $10^{-3}\,M_\odot$) of self-consistent, two-dimensional explosion models of core-collapse supernovae for six progenitor stars with different initial masses. Three models have initial masses near the low-mass end of the supernova range, $8.8\,M_\odot$ (e8.8; electron-capture supernova), $9.6\,M_\odot$ (z9.6), and $8.1\,M_\odot$ (u8.1), with initial metallicities of 1, 0, and $10^{-4}$ times the solar metallicity, respectively. The other three are solar-metallicity models with initial masses of $11.2\,M_\odot$ (s11), $15\,M_\odot$ (s15), and $27\,M_\odot$ (s27). The low-mass models e8.8, z9.6, and u8.1 exhibit high production factors (nucleosynthetic abundances relative to the solar ones) of 100--200 for light trans-iron elements from Zn to Zr. This is associated with appreciable ejection of neutron-rich matter in these models. Remarkably, the nucleosynthetic outcomes for progenitors e8.8 and z9.6 are almost identical, including interesting productions of $^{48}$Ca and $^{60}$Fe, irrespective of their quite different (O-Ne-Mg and Fe) cores prior to collapse. In the more massive models s11, s15, and s27, several proton-rich isotopes of light trans-iron elements, including the $p$-isotope $^{92}$Mo (for s27) are made, up to production factors of $\sim$30. Both electron-capture and core-collapse supernovae near the low-mass end can therefore be dominant contributors to the Galactic inventory of light trans-iron elements from Zn to Zr and probably $^{48}$Ca and live $^{60}$Fe. The innermost ejecta of more massive supernovae may have only sub-dominant contributions to the chemical enrichment of the Galaxy except for $^{92}$Mo.

Partially Ionized Plasmas in Astrophysics

Abstract: Partially ionized plasmas are found across the Universe in many different astrophysical environments. They constitute an essential ingredient of the solar atmosphere, molecular clouds, planetary ionospheres and protoplanetary disks, among other environments, and display a richness of physical effects which are not present in fully ionized plasmas. This review provides an overview of the physics of partially ionized plasmas, including recent advances in different astrophysical areas in which partial ionization plays a fundamental role. We outline outstanding observational and theoretical questions and discuss possible directions for future progress.

A Systematic Survey of the Effects of Wind Mass Loss Algorithms on the Evolution of Single Massive Stars

Abstract: Mass loss is a key uncertainty in the evolution of massive stars. Stellar evolution calculations must employ parametric algorithms for mass loss, and usually only include stellar winds. We carry out a parameter study of the effects of wind mass loss on massive star evolution using the open-source stellar evolution code MESA. We provide a systematic comparison of wind mass loss algorithms for solar-metallicity, nonrotating, single stars in the initial mass range of $15-35\,M_\odot$. We consider combinations drawn from two hot phase algorithms, three cool phase algorithms, and two Wolf-Rayet algorithms. We consider linear wind efficiency scale factors of $1$, $0.33$, and $0.1$ to account for reductions in mass loss rates due to wind inhomogeneities. We find that the initial to final mass mapping for each zero-age main-sequence (ZAMS) mass has a $\sim 50\%$ uncertainty if all algorithm combinations and wind efficiencies are considered. The ad-hoc efficiency scale factor dominates this uncertainty. While the final total mass and internal structure of our models vary tremendously with mass loss treatment, final observable parameters are much less sensitive for ZAMS mass $\lesssim 30\,M_\odot$. This indicates that uncertainty in wind mass loss does not negatively affect estimates of the ZAMS mass of most single-star supernova progenitors from pre-explosion observations. Furthermore, we show that the internal structure of presupernova stars is sensitive to variations in both main sequence and post main-sequence mass loss. We find that the compactness parameter $\xi\propto M/R(M)$ varies by as much as $30\%$ for a given ZAMS mass evolved with different wind efficiencies and mass loss algorithm combinations. [abridged]

Common Envelope Wind Tunnel: Coefficients of Drag and Accretion in a Simplified Context for Studying Flows Around Objects Embedded Within Stellar Envelopes

Abstract: This paper examines the properties of flows around objects embedded within common envelopes in the simplified context of a "wind tunnel." We establish characteristic relationships between key common envelope flow parameters like the Mach number and density scale height. Our wind tunnel is a three-dimensional, cartesian geometry hydrodynamic simulation setup that includes the gravity of the primary and secondary stars and allows us to study the coefficients of drag and accretion experienced by the embedded object. Accretion and drag lead to a transformation of an embedded object and its orbit during a common envelope phase. We present two suites of simulations spanning a range of density gradients and Mach numbers -- relevant for flow near the limb of a stellar envelope to the deep interior. In one suite, we adopt an ideal gas adiabatic exponent of $\gamma=5/3$, in the other, $\gamma=4/3$. We find that coefficients of drag rise in flows with steeper density gradients and that coefficients of drag and accretion are consistently higher in the more compressible, $\gamma=4/3$ flow. We illustrate the impact of these newly derived coefficients by integrating the inspiral of a secondary object through the envelopes of $3M_\odot$ ($\gamma\approx5/3$) and $80M_\odot$ ($\gamma\approx4/3$) giants. In these examples, we find a relatively rapid initial inspiral because, near the stellar limb, dynamical friction drag is generated mainly from dense gas focussed from deeper within the primary-star's envelope. This rapid initial inspiral timescale carries potential implications for the timescale of transients from early common envelope interaction.

Protoplanetary Disks in $\rho$ Ophiuchus as Seen From ALMA

Abstract: We present a high angular resolution ($\sim 0.2^{\prime\prime}$), high sensitivity ($\sigma \sim 0.2$ mJy) survey of the 870 $\mu$m continuum emission from the circumstellar material around 49 pre-main sequence stars in the $\rho$ Ophiuchus molecular cloud. Because most millimeter instruments have resided in the northern hemisphere, this represents the largest high-resolution, millimeter-wave survey of the circumstellar disk content of this cloud. Our survey of 49 systems comprises 63 stars; we detect disks associated with 29 single sources, 11 binaries, 3 triple systems and 4 transition disks. We present flux and radius distributions for these systems; in particular, this is the first presentation of a reasonably complete probability distribution of disk radii at millimeter-wavelengths. We also compare the flux distribution of these protoplanetary disks with that of the disk population of the Taurus-Auriga molecular cloud. We find that disks in binaries are both significantly smaller and have much less flux than their counterparts around isolated stars. We compute truncation calculations on our binary sources and find that these disks are too small to have been affected by tidal truncation and posit some explanations for this. Lastly, our survey found 3 candidate gapped disks, one of which is a newly identified transition disk with no signature of a dip in infrared excess in extant observations.

Gaia Data Release 1. Testing the parallaxes with local Cepheids and RR Lyrae stars

Abstract: Parallaxes for 331 classical Cepheids, 31 Type II Cepheids and 364 RR Lyrae stars in common between Gaia and the Hipparcos and Tycho-2 catalogues are published in Gaia Data Release 1 (DR1) as part of the Tycho-Gaia Astrometric Solution (TGAS). In order to test these first parallax measurements of the primary standard candles of the cosmological distance ladder, that involve astrometry collected by Gaia during the initial 14 months of science operation, we compared them with literature estimates and derived new period-luminosity ($PL$), period-Wesenheit ($PW$) relations for classical and Type II Cepheids and infrared $PL$, $PL$-metallicity ($PLZ$) and optical luminosity-metallicity ($M_V$-[Fe/H]) relations for the RR Lyrae stars, with zero points based on TGAS. The new relations were computed using multi-band ($V,I,J,K_{\mathrm{s}},W_{1}$) photometry and spectroscopic metal abundances available in the literature, and applying three alternative approaches: (i) by linear least squares fitting the absolute magnitudes inferred from direct transformation of the TGAS parallaxes, (ii) by adopting astrometric-based luminosities, and (iii) using a Bayesian fitting approach. TGAS parallaxes bring a significant added value to the previous Hipparcos estimates. The relations presented in this paper represent first Gaia-calibrated relations and form a "work-in-progress" milestone report in the wait for Gaia-only parallaxes of which a first solution will become available with Gaia's Data Release 2 (DR2) in 2018.

Molecular gas in debris disks around young A-type stars

Abstract: According to the current paradigm of circumstellar disk evolution, gas-rich primordial disks evolve into gas-poor debris disks compose of second-generation dust. To explore the transition between these phases, we searched for $^{12}$CO, $^{13}$CO, and C$^{18}$O emission in seven dust-rich debris disks around young A-type stars, using ALMA in Band 6. We discovered molecular gas in three debris disks. In all these disks, the $^{12}$CO line was optically thick, highlighting the importance of less abundant molecules in reliable mass estimates. Supplementing our target list by literature data, we compiled a volume-limited sample of dust-rich debris disks around young A-type stars within 150 pc. We obtained a CO detection rate of 11/16 above a $^{12}$CO J=2$-$1 line luminosity threshold of $\sim 1.4 \times 10 ^4$ Jykms$^{-1}$pc$^2$ in the sample. This high incidence implies that the presence of CO gas in bright debris disks around young A-type stars is likely more the rule than the exception. Interestingly, dust-rich debris disks around young FG-type stars exhibit, with the same detectability threshold as for A-type stars, significantly lower gas incidence. While the transition from protoplanetary to debris phase is associated with a drop of dust content, our results exhibit a large spread in the CO mass in our debris sample, with peak values comparable to those in protoplanetary Herbig Ae disks. In the particularly CO-rich debris systems the gas may have primordial origin, characteristic of a hybrid disk.

Bombs and flares at the surface and lower atmosphere of the Sun

Abstract: A spectacular manifestation of solar activity is the appearance of transient brightenings in the far wings of the H$\alpha$ line, known as Ellerman bombs (EBs). Recent observations obtained by the Interface Region Imaging Spectrograph (IRIS) have revealed another type of plasma "bombs" (UV bursts) with high temperatures of perhaps up to $8\times 10^4$ K within the cooler lower solar atmosphere. Realistic numerical modeling showing such events is needed to explain their nature. Here, we report on 3D radiative magneto-hydrodynamic simulations of magnetic flux emergence in the solar atmosphere. We find that ubiquitous reconnection between emerging bipolar magnetic fields can trigger EBs in the photosphere, UV bursts in the mid/low chromosphere and small (nano-/micro-) flares ($10^6$ K) in the upper chromosphere. These results provide new insights on the emergence and build up of the coronal magnetic field and the dynamics and heating of the solar surface and lower atmosphere.

The Threatening Environment of the TRAPPIST-1 Planets

Abstract: Recently, four additional Earth-mass planets were discovered orbiting the nearby ultracool M8 dwarf TRAPPIST-1, making a remarkable total of seven planets with equilibrium temperatures compatible with the presence of liquid water on their surface. Temperate terrestrial planets around an M-dwarf orbit close to their parent star, rendering their atmospheres vulnerable to erosion by the stellar wind and energetic electromagnetic and particle radiation. Here, we use state-of-the-art 3D magnetohydrodynamic models to simulate the wind around TRAPPIST-1 and study the conditions at each planetary orbit. All planets experience a stellar wind pressure between $10^3$ and $10^5$ times the solar wind pressure on Earth. All orbits pass through wind pressure changes of an order of magnitude and most planets spend a large fraction of their orbital period in the sub-Alfv\'enic regime. For plausible planetary magnetic field strengths, all magnetospheres are greatly compressed and undergo much more dynamic change than that of the Earth. The planetary magnetic fields connect with the stellar radial field over much of the planetary surface, allowing direct flow of stellar wind particles onto the planetary atmosphere. These conditions could result in strong atmospheric stripping and evaporation and should be taken into account for any realistic assessment of the evolution and habitability of the TRAPPIST-1 planets.

The State-of-the-Art HST Astro-Photometric Analysis of the core of $\omega$ Centauri. I. The Catalog

Abstract: We have constructed the most-comprehensive catalog of photometry and proper motions ever assembled for a globular cluster (GC). The core of $\omega$Cen has been imaged over 650 times through WFC3's UVIS and IR channels for the purpose of detector calibration. There exist from 4 to over 60 exposures through each of 26 filters, stretching continuously from F225W in the UV to F160W in the infrared. Furthermore, the 11-year baseline between these data and a 2002 ACS survey has allowed us to more than double the proper-motion accuracy and triple the number of well-measured stars compared to our previous groundbreaking effort. This totally unprecedented complete spectral coverage for over 470,000 stars within the cluster's core, from the tip of the red-giant branch down to the white dwarfs, provides the best astro-photometric observational data base yet to understand the multiple-population phenomenon in any GC. In this first paper of the series we describe in detail the data-reduction processes and deliver the astro-photometric catalog to the astronomical community.

The Gaia-ESO Survey. Mg-Al anti-correlation in iDR4 globular clusters

Abstract: We use Gaia-ESO Survey iDR4 data to explore the Mg-Al anti-correlation in globular clusters, that were observed as calibrators, as a demonstration of the quality of Gaia-ESO Survey data and analysis. The results compare well with the available literature, within 0.1 dex or less, after a small (compared to the internal spreads) offset between the UVES and the GIRAFFE data of 0.10-0.15 dex was taken into account. In particular, we present for the first time data for NGC 5927, one of the most metal-rich globular clusters studied in the literature so far with [Fe/H]=-0.49 dex, that was included to connect with the open cluster regime in the Gaia-ESO Survey internal calibration. The extent and shape of the Mg-Al anti-correlation provide strong constraints on the multiple population phenomenon in globular clusters. In particular, we studied the dependency of the Mg-Al anti-correlation extension with metallicity, present-day mass, and age of the clusters, using GES data in combination with a large set of homogenized literature measurements. We find a dependency with both metallicity and mass, that is evident when fitting for the two parameters simultaneously, but no significant dependency with age. We confirm that the Mg-Al anti-correlation is not seen in all clusters, but disappears for the less massive or most metal-rich ones. We also use our dataset to see whether a normal anti-correlation would explain the low [Mg/$\alpha$] observed in some extragalactic globular clusters, but find that none of the clusters in our sample can reproduce it, and more extreme chemical compositions (like the one of NGC 2419) would be required. We conclude that GES iDR4 data already meet the requirements set by the main survey goals, and can be used to study in detail globular clusters even if the analysis procedures were not specifically designed for them.

Physical Properties and Catalogue of EW-type Eclipsing Binaries Observed by LAMOST

Abstract: Numerous EWs were discovered by several deep photometric survey and there are about 40785 EW-type binary systems listed in the international variable star index (VSX) by March 13, 2017. 7938 of them were observed by LAMOST by November 30, 2016 and their spectral types were given. Stellar atmospheric parameters of 5363 EW-type binary stars were determined based on good spectroscopic observations. In the paper, those EWs were catalogued and their properties are analyzed. The distributions of the orbital period (P), the effect temperature (T), the gravitational acceleration (Log(g)), the metallicity ([Fe/H]) and the radial velocity (RV) are presented for those observed EW-type systems. It is shown that about 80.6\% sample stars have metallicity below zero indicating that EW-type systems are old stellar population. This is in agreement with the the conclusion that the EW binaries are formed from moderately close binaries through angular momentum loss via magnetic braking that takes a few hundred million to a few billion years. The unusual high metallicities of a few percent of EWs may be caused by contaminating of material from the evolution of unseen neutron stars and black holes in the systems. The correlations between the orbital period and the effect temperature, the gravitational acceleration and the metallicity are presented and their scatters are mainly caused by (i) the presence of the third bodies and (ii) the wrong determined periods sometimes. It is shown that some EW contain evolved component stars and the physical properties of EWs are mainly depending on their orbital periods. It is found that the extremely short-period EWs may be older than their long-period cousins because they have lower metallicities. This reveals that they have a longer timescale of pre-contact evolution and their formation and evolution are mainly driven by angular momentum loss via magnetic braking.