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

Estimating distances from parallaxes IV: Distances to 1.33 billion stars in Gaia Data Release 2

Abstract: For the vast majority of stars in the second Gaia data release, reliable distances cannot be obtained by inverting the parallax. A correct inference procedure must instead be used to account for the nonlinearity of the transformation and the asymmetry of the resulting probability distribution. Here we infer distances to essentially all 1.33 billion stars with parallaxes published in the second \gaia\ data release. This is done using a weak distance prior that varies smoothly as a function of Galactic longitude and latitude according to a Galaxy model. The irreducible uncertainty in the distance estimate is characterized by the lower and upper bounds of an asymmetric confidence interval. Although more precise distances can be estimated for a subset of the stars using additional data (such as photometry), our goal is to provide purely geometric distance estimates, independent of assumptions about the physical properties of, or interstellar extinction towards, individual stars. We analyse the characteristics of the catalogue and validate it using clusters. The catalogue can be queried on the Gaia archive using ADQL at this http URL and downloaded from this http URL .

Binary companions of evolved stars in APOGEE DR14: Search method and catalog of ~5,000 companions

Abstract: Multi-epoch radial velocity measurements of stars can be used to identify stellar, sub-stellar, and planetary-mass companions. Even a small number of observation epochs can be informative about companions, though there can be multiple qualitatively different orbital solutions that fit the data. We have custom-built a Monte Carlo sampler (The Joker) that delivers reliable (and often highly multi-modal) posterior samplings for companion orbital parameters given sparse radial-velocity data. Here we use The Joker to perform a search for companions to 96,231 red-giant stars observed in the APOGEE survey (DR14) with $\geq 3$ spectroscopic epochs. We select stars with probable companions by making a cut on our posterior belief about the amplitude of the stellar radial-velocity variation induced by the orbit. We provide (1) a catalog of 320 companions for which the stellar companion properties can be confidently determined, (2) a catalog of 4,898 stars that likely have companions, but would require more observations to uniquely determine the orbital properties, and (3) posterior samplings for the full orbital parameters for all stars in the parent sample. We show the characteristics of systems with confidently determined companion properties and highlight interesting systems with candidate compact object companions.

Gaia Data Release 2: Observational Hertzsprung-Russell diagrams

Abstract: We highlight the power of the Gaia DR2 in studying many fine structures of the Hertzsprung-Russell diagram (HRD). Gaia allows us to present many different HRDs, depending in particular on stellar population selections. We do not aim here for completeness in terms of types of stars or stellar evolutionary aspects. Instead, we have chosen several illustrative examples. We describe some of the selections that can be made in Gaia DR2 to highlight the main structures of the Gaia HRDs. We select both field and cluster (open and globular) stars, compare the observations with previous classifications and with stellar evolutionary tracks, and we present variations of the Gaia HRD with age, metallicity, and kinematics. Late stages of stellar evolution such as hot subdwarfs, post-AGB stars, planetary nebulae, and white dwarfs are also analysed, as well as low-mass brown dwarf objects. The Gaia HRDs are unprecedented in both precision and coverage of the various Milky Way stellar populations and stellar evolutionary phases. Many fine structures of the HRDs are presented. The clear split of the white dwarf sequence into hydrogen and helium white dwarfs is presented for the first time in an HRD. The relation between kinematics and the HRD is nicely illustrated. Two different populations in a classical kinematic selection of the halo are unambiguously identified in the HRD. Membership and mean parameters for a selected list of open clusters are provided. They allow drawing very detailed cluster sequences, highlighting fine structures, and providing extremely precise empirical isochrones that will lead to more insight in stellar physics. Gaia DR2 demonstrates the potential of combining precise astrometry and photometry for large samples for studies in stellar evolution and stellar population and opens an entire new area for HRD-based studies.

Evidence for a systematic offset of -80 micro-arcseconds in the Gaia DR2 parallaxes

Abstract: We reprise the analysis of Stassun & Torres (2016), comparing the parallaxes of the eclipsing binaries reported in that paper to the parallaxes newly reported in the Gaia second data release (DR2). We find evidence for a systematic offset of $-82 \pm 33$ micro-arcseconds, in the sense of the Gaia parallaxes being too small, for brightnesses $(G \lesssim 12)$ and for distances (0.03--3 kpc) in the ranges spanned by the eclipsing binary sample. The offset does not appear to depend strongly on distance within this range, though there is marginal evidence that the offset increases (becomes slightly more negative) for distances $\gtrsim 1$ kpc, up to the 3 kpc distances probed by the test sample. The offset reported here is consistent with the expectation that global systematics in the Gaia DR2 parallaxes are below 100 micro-arcseconds.

New Parallaxes of Galactic Cepheids from Spatially Scanning the Hubble Space Telescope: Implications for the Hubble Constant

Abstract: We present new parallax measurements of 7 long-period (> 10 days) Milky Way Cepheids (SS CMa, XY Car, VY Car, VX Per, WZ Sgr, X Pup and S Vul) using astrometry from spatial scanning of WFC3 on HST. Observations were obtained at 6 month intervals over 4 years. The distances are 1.7--3.6 kpc with a mean precision of 45 microarcseconds and a best of 29 microarcseconds (SNR = 14). The accuracy of the parallaxes is demonstrated through independent analyses of >100 reference stars. This raises to 10 the number of long-period Cepheids with significant parallax measurements, 8 obtained from this program. We also present high-precision F555W, F814W, and F160W magnitudes of these Cepheids, allowing a direct, zeropoint-independent comparison to >1800 extragalactic Cepheids in the hosts of 19 SNeIa. This sample addresses two outstanding systematic uncertainties affecting prior comparisons of Milky Way and extragalactic Cepheids used to calibrate H_0: their dissimilarity of periods and photometric systems. Comparing the new parallaxes to their predicted values derived from reversing the distance ladder gives a ratio (or independent scale for H_0) of 1.037+/-0.036, consistent with no change and inconsistent at the 3.5 sigma level with a ratio of 0.91 needed to match the value predicted by Planck+LCDM. Using these data instead to augment the Riess et al. (2016) measurement of H_0 improves the precision to 2.3%, yielding 73.48+/-1.66 km/s/Mpc, and tension with Planck+LCDM increases to 3.7 sigma. The future combination of Gaia parallaxes and HST spatial scanning photometry of 50 Milky Way Cepheids can support a < 1% calibration of H_0.

BANYAN. XI. The BANYAN $\Sigma$ multivariate Bayesian algorithm to identify members of young associations within 150 pc.

Abstract: BANYAN {\Sigma} is a new Bayesian algorithm to identify members of young stellar associations within 150 pc of the Sun It includes 27 young associations with ages in the range ~1-800 Myr, modelled with multivariate Gaussians in 6-dimensional XYZUVW space It is the first such multi-associations classification tool to include the nearest sub-groups of the Sco-Cen OB star-forming region, the IC 2602, IC 2391, Pleiades and Platais 8 clusters, and the {\rho} Ophiuci, Corona Australis, and Taurus star-formation regions A model of field stars is built from a mixture of multivariate Gaussians based on the Besan\c{c}on Galactic model The algorithm can derive membership probabilities for objects with only sky coordinates and proper motion, but can also include parallax and radial velocity measurements, as well as spectrophotometric distance constraints from sequences in color-magnitude or spectral type-magnitude diagrams BANYAN {\Sigma} benefits from an analytical solution to the Bayesian marginalization integrals that makes it more accurate and significantly faster than its predecessor BANYAN II A contamination versus hit rate analysis is presented and demonstrates that BANYAN {\Sigma} achieves a better classification performance than other moving group classification tools, especially in terms of cross-contamination between young associations An updated list of bona fide members in the 27 young associations, augmented by the Gaia-DR1 release, are presented This new tool will make it possible to analyze large data sets such as the upcoming Gaia-DR2 to identify new young stars IDL and Python versions of BANYAN {\Sigma} are made available with this publication (shortened)

How to Constrain Your M dwarf II: the mass-luminosity-metallicity relation from 0.075 to 0.70$M_\odot$

Abstract: The mass-luminosity relation for late-type stars has long been a critical tool for estimating stellar masses. However, there is growing need for both a higher-precision relation and a better understanding of systematic effects (e.g., metallicity). Here we present an empirical relationship between Mks and mass spanning $0.075M_\odot

Disks ARound TTauri Stars with Sphere (DARTTS-S) I: Sphere / IRDIS Polarimetric Imaging of 8 prominent TTauri Disks

Abstract: We present the first part of our DARTTS-S (Disks ARound TTauri Stars with SPHERE) survey: Observations of 8 TTauri stars which were selected based on their strong (sub-)mm excesses using SPHERE / IRDIS polarimetric differential imaging (PDI) in the J and H bands. All observations successfully detect the disks, which appear vastly different in size, from $\approx$80 au in scattered light to $>$400 au, and display total polarized disk fluxes between 0.06% and 0.89% of the stellar flux. For five of these disks, we are able to determine the three-dimensional structure and the flaring of the disk surface, which appears to be relatively consistent across the different disks, with flaring exponents $\alpha$ between $\approx$1.1 and $\approx$1.6. We also confirm literature results w.r.t. the inclination and position angle of several of our disk, and are able to determine which side is the near side of the disk in most cases. While there is a clear trend of disk mass with stellar ages ($\approx$1 Myr to $>$10 Myr), no correlations of disk structures with age were found. There are also no correlations with either stellar mass or sub-mm flux. We do not detect significant differences between the J and H bands. However, we note that while a high fraction (7/8) of the disks in our sample show ring-shaped sub-structures, none of them display spirals, in contrast to the disks around more massive Herbig Ae/Be stars, where spiral features are common.

Confirmation of the ${\rm \it Gaia}$ DR2 parallax zero-point offset using asteroseismology and spectroscopy in the ${\rm \it Kepler}$ field

Abstract: We present an independent confirmation of the zero-point offset of ${\rm \it Gaia}$ Data Release 2 (DR2) parallaxes using asteroseismic data of evolved stars in the ${\rm \it Kepler}$ field. Using well-characterized red giant branch (RGB) stars from the APOKASC-2 catalogue we identify a ${\rm \it Gaia}$ astrometric pseudo-color ($ u_{\rm eff}$)- and ${\rm \it Gaia}$ $G$-band magnitude-dependent zero-point offset of $\varpi_{\rm seis} - \varpi_{Gaia} = 52.8 \pm 2.4 {\rm\ (rand.)} \pm 8.6 {\rm\ (syst.)} - (150.7 \pm 22.7)( u_{\rm eff} - 1.5) - (4.21 \pm 0.77)(G - 12.2) \mu{\rm as}$, in the sense that ${\rm \it Gaia}$ parallaxes are too small. The offset is found in high and low-extinction samples, as well as among both shell H-burning red giant stars and core He-burning red clump stars. We show that errors in the asteroseismic radius and temperature scales may be distinguished from errors in the ${\rm \it Gaia}$ parallax scale. We estimate systematic effects on the inferred global ${\rm \it Gaia}$ parallax offset, $c$, due to radius and temperature systematics, as well as choices in bolometric correction and the adopted form for ${\rm \it Gaia}$ parallax spatial correlations. Because of possible spatially-correlated parallax errors, as discussed by the ${\rm \it Gaia}$ team, our ${\rm \it Gaia}$ parallax offset model is specific to the ${\rm \it Kepler}$ field, but broadly compatible with the magnitude- and color-dependent offset inferred by the ${\rm \it Gaia}$ team and several subsequent investigations using independent methods.

Kinematic evidence for an embedded protoplanet in a circumstellar disc

Abstract: Discs of gas and dust surrounding young stars are the birthplace of planets. However, direct detection of protoplanets forming within discs has proved elusive to date. We present the detection of a large, localized deviation from Keplerian velocity in the protoplanetary disc surrounding the young star HD163296. The observed velocity pattern is consistent with the dynamical effect of a two Jupiter-mass planet orbiting at a radius $\approx$ 260au from the star.

Gaia Data Release 2: Summary of the variability processing & analysis results

Abstract: The Gaia Data Release 2 (DR2): we summarise the processing and results of the identification of variable source candidates of RR Lyrae stars, Cepheids, long period variables (LPVs), rotation modulation (BY Dra-type) stars, delta Scuti & SX Phoenicis stars, and short-timescale variables. In this release we aim to provide useful but not necessarily complete samples of candidates. The processed Gaia data consist of the G, BP, and RP photometry during the first 22 months of operations as well as positions and parallaxes. Various methods from classical statistics, data mining and time series analysis were applied and tailored to the specific properties of Gaia data, as well as various visualisation tools. The DR2 variability release contains: 228'904 RR Lyrae stars, 11'438 Cepheids, 151'761 LPVs, 147'535 stars with rotation modulation, 8'882 delta Scuti & SX Phoenicis stars, and 3'018 short-timescale variables. These results are distributed over a classification and various Specific Object Studies (SOS) tables in the Gaia archive, along with the three-band time series and associated statistics for the underlying 550'737 unique sources. We estimate that about half of them are newly identified variables. The variability type completeness varies strongly as function of sky position due to the non-uniform sky coverage and intermediate calibration level of this data. The probabilistic and automated nature of this work implies certain completeness and contamination rates which are quantified so that users can anticipate their effects. This means that even well-known variable sources can be missed or misidentified in the published data. The DR2 variability release only represents a small subset of the processed data. Future releases will include more variable sources and data products; however, DR2 shows the (already) very high quality of the data and great promise for variability studies.

Gaps and Rings in an ALMA Survey of Disks in the Taurus Star-forming Region

Abstract: Rings are the most frequently revealed substructure in ALMA dust observations of protoplanetary disks, but their origin is still hotly debated. In this paper, we identify dust substructures in 12 disks and measure their properties to investigate how they form. This subsample of disks is selected from a high-resolution ($\sim0.12''$) ALMA 1.33 mm survey of 32 disks in the Taurus star-forming region, which was designed to cover a wide range of sub-mm brightness and to be unbiased to previously known substructures. While axisymmetric rings and gaps are common within our sample, spiral patterns and high contrast azimuthal asymmetries are not detected. Fits of disk models to the visibilities lead to estimates of the location and shape of gaps and rings, the flux in each disk component, and the size of the disk. The dust substructures occur across a wide range of stellar mass and disk brightness. Disks with multiple rings tend to be more massive and more extended. The correlation between gap locations and widths, the intensity contrast between rings and gaps, and the separations of rings and gaps could all be explained if most gaps are opened by low-mass planets (super-Earths and Neptunes) in the condition of low disk turbulence ($\alpha=10^{-4}$). The gap locations are not well correlated with the expected locations of CO and N$_2$ ice lines, so condensation fronts are unlikely to be a universal mechanism to create gaps and rings, though they may play a role in some cases.

Gaia Data Release 2: Specific characterisation and validation of all-sky Cepheids and RR Lyrae stars

Abstract: The Gaia second Data Release (DR2) presents a first mapping of full-sky RR Lyrae stars and Cepheids observed by the spacecraft during the initial 22 months of science operations. The Specific Object Study (SOS) pipeline, developed to validate and fully characterise Cepheids and RR Lyrae stars (SOS Cep&RRL) observed by Gaia, has been presented in the documentation and papers accompanying the Gaia first Data Release. Here we describe how the SOS pipeline was modified to allow for processing the Gaia multiband (G, G_BP and G_RP) time series photometry of all-sky candidate variables and produce specific results for confirmed RR Lyrae stars and Cepheids that are published in the DR2 catalogue. The SOS Cep&RRL processing uses tools such as the period-amplitude and the period-luminosity relations in the G band. For the analysis of the Gaia DR2 candidates we also used tools based on the G_BP and G_RP photometry, such as the period-Wesenheit relation in (G,G_RP). Multiband time series photometry and characterisation by the SOS Cep&RRL pipeline are published in Gaia DR2 for 150,359 such variables (9,575 classified as Cepheids and 140,784 as RR Lyrae stars) distributed all over the sky. The sample includes variables in 87 globular clusters and 14 dwarf galaxies. To the best of our knowledge, as of 25 April 2018, variability of 50,570 of these sources (350 Cepheids and 50,220 RR Lyrae stars) is not known in the literature, hence likely they are new discoveries by Gaia. An estimate of the interstellar absorption is published for 54,272 fundamental-mode RR Lyrae stars from a relation based on the G-band amplitude and the pulsation period. Metallicities derived from the Fourier parameters of the light curves are also released for 64,932 RR Lyrae stars and 3,738 fundamental-mode classical Cepheids with period shorter than 6.3 days.

The APOGEE-2 Survey of the Orion Star Forming Complex II: Six-dimensional structure

Abstract: We present an analysis of spectrosopic and astrometric data from APOGEE-2 and Gaia DR2 to identify structures towards the Orion Complex. By applying a hierarchical clustering algorithm to the 6-dimensional stellar data, we identify spatially and/or kinematically distinct groups of young stellar objects with ages ranging from 1 to 12 Myr. We also investigate the star forming history within the Orion Complex, and identify peculiar sub-clusters. With this method we reconstruct the older populations in the regions that are presently largely devoid of molecular gas, such as Orion C (which includes the $\sigma$ Ori cluster), and Orion D (the population that traces Ori OB1a, OB1b, and Orion X). We report on the distances, kinematics, and ages of the groups within the Complex. The Orion D groups is in the process of expanding. On the other hand, Orion B is still in the process of contraction. In $\lambda$ Ori the proper motions are consistent with a radial expansion due to an explosion from a supernova; the traceback age from the expansion exceeds the age of the youngest stars formed near the outer edges of the region, and their formation would have been triggered when they were half-way from the cluster center to their current positions. We also present a comparison between the parallax and proper motion solutions obtained by Gaia DR2, and those obtained towards star-forming regions by Very Long Baseline Array.

Massive runaways and walkaway stars

Abstract: Anticipating the kinematic constraints from the Gaia mission, we perform an extensive numerical study of the evolution of massive binary systems to predict the peculiar velocities that stars obtain when their companion collapses and disrupts the system. Our aim is to (1) identify which predictions are robust against model uncertainties and assess their implications, (2) investigate which physical processes leave a clear imprint and may therefore be constrained observationally and (3) provide a suite of publicly available model predictions. We find that $22_{-8}^{+26}$% of all massive binary systems merge prior to the first core collapse in the system. Of the remainder, $86_{-9}^{+11}$% become unbound because of the core-collapse. Remarkably, this rarely produce runaway stars (i.e., stars with velocities above 30 km/s). These are outnumbered by more than an order of magnitude by slower unbound companions, or "walkaway stars". This is a robust outcome of our simulations and is due to the reversal of the mass ratio prior to the explosion and widening of the orbit, as we show analytically and numerically. We estimate a $10^{+5}_{-8}$% of massive stars to be walkaways and only $0.5^{+1.0}_{-0.4}$% to be runaways, nearly all of which have accreted mass from their companion. Our findings are consistent with earlier studies, however the low runaway fraction we find is in tension with observed fractions 10%. If Gaia confirms these high fractions of massive runaway stars resulting from binaries, it would imply that we are currently missing physics in the binary models. Finally, we show that high end of the mass distributions of runaway stars is very sensitive to the assumed black hole natal kicks and propose this as a potentially stringent test for the explosion mechanism. We discuss companions remaining bound which can evolve into X-ray and gravitational wave sources.

The White Dwarf Initial-Final Mass Relation for Progenitor Stars From 0.85 to 7.5 M$_\odot$

Abstract: We present the initial-final mass relation (IFMR) based on the self-consistent analysis of Sirius B and 79 white dwarfs from 13 star clusters. We have also acquired additional signal on eight white dwarfs previously analyzed in the NGC 2099 cluster field, four of which are consistent with membership. These reobserved white dwarfs have masses ranging from 0.72 to 0.97 M$_\odot$, with initial masses from 3.0 to 3.65 M$_\odot$, where the IFMR has an important change in slope that these new data help to observationally confirm. In total, this directly measured IFMR has small scatter ($\sigma$ = 0.06 M$_\odot$) and spans from progenitors of 0.85 to 7.5 M$_\odot$. Applying two different stellar evolutionary models to infer two different sets of white dwarf progenitor masses shows that when the same model is also used to derive the cluster ages, the resulting IFMR has weak sensitivity to the adopted model at all but the highest initial masses ($>$5.5 M$_\odot$). The non-linearity of the IFMR is also clearly observed with moderate slopes at lower masses (0.08 M$_{\rm final}$/M$_{\rm initial}$) and higher masses (0.11 M$_{\rm final}$/M$_{\rm initial}$) that are broken up by a steep slope (0.19 M$_{\rm final}$/M$_{\rm initial}$) between progenitors from 2.85 to 3.6 M$_\odot$. This IFMR shows total stellar mass loss ranges from 33\% of M$_{\rm initial}$ at 0.83 M$_\odot$ to 83\% of M$_{\rm initial}$ at 7.5 M$_\odot$. Testing this total mass loss for dependence on progenitor metallicity, however, finds no detectable sensitivity across the moderate range of -0.15 $<$ [Fe/H] $<$ +0.15.

A First Catalog of Variable Stars Measured by the Asteroid Terrestrial-impact Last Alert System (ATLAS)

Abstract: The Asteroid Terrestrial-impact Last Alert System (ATLAS) carries out its primary planetary defense mission by surveying about 13000 deg^2 at least four times per night. The resulting data set is useful for the discovery of variable stars to a magnitude limit fainter than r~18, with amplitudes down to 0.01 mag for bright objects. Here we present a Data Release One catalog of variable stars based on analyzing 142 million stars measured at least 100 times in the first two years of ATLAS operations. Using a Lomb-Scargle periodogram and other variability metrics, we identify 4.7 million candidate variables which we analyze in detail. Through Space Telescope Science Institute, we publicly release lightcurves for all of them, together with a vector of 169 classification features for each star. We do this at the level of unconfirmed candidate variables in order to provide the community with a large set of homogeneously analyzed photometry and avoid pre-judging which types of objects others may find most interesting. We use machine learning to classify the candidates into fifteen different broad categories based on lightcurve morphology. About 10% (430,000 stars) pass extensive tests designed to screen out spurious variability detections: we label these as `probable' variables. Of these, 230,000 receive specific classifications as eclipsing binaries, pulsating, Mira-type, or sinusoidal variables: these are the `classified' variables. New discoveries among the probable variables number more than 300,000, while 150,000 of the classified variables are new, including about 10,000 pulsating variables, 2,000 Mira stars, and 70,000 eclipsing binaries.

The Second APOKASC Catalog: The Empirical Approach.

Abstract: We present a catalog of stellar properties for a large sample of 6676 evolved stars with APOGEE spectroscopic parameters and \textit{Kepler} asteroseismic data analyzed using five independent techniques. Our data includes evolutionary state, surface gravity, mean density, mass, radius, age, and the spectroscopic and asteroseismic measurements used to derive them. We employ a new empirical approach for combining asteroseismic measurements from different methods, calibrating the inferred stellar parameters, and estimating uncertainties. With high statistical significance, we find that asteroseismic parameters inferred from the different pipelines have systematic offsets that are not removed by accounting for differences in their solar reference values. We include theoretically motivated corrections to the large frequency spacing ($\Delta u$) scaling relation, and we calibrate the zero point of the frequency of maximum power ($ u_{\rm max}$) relation to be consistent with masses and radii for members of star clusters. For most targets, the parameters returned by different pipelines are in much better agreement than would be expected from the pipeline-predicted random errors, but 22\% of them had at least one method not return a result and a much larger measurement dispersion. This supports the usage of multiple analysis techniques for asteroseismic stellar population studies. The measured dispersion in mass estimates for fundamental calibrators is consistent with our error model, which yields median random and systematic mass uncertainties for RGB stars of order 4\%. Median random and systematic mass uncertainties are at the 9\% and 8\% level respectively for RC stars.

CARMENES input catalogue of M dwarfs: III. Rotation and activity from high-resolution spectroscopic observations

Abstract: CARMENES is a spectrograph for radial velocity surveys of M dwarfs with the aim of detecting Earth-mass planets orbiting in the habitable zones of their host stars. To ensure an optimal use of the CARMENES Guaranteed Time Observations, in this paper we investigate the correlation of activity and rotation for approximately 2200 M dwarfs, ranging in spectral type from M0.0 V to M9.0 V. We present new high-resolution spectroscopic observations with FEROS, CAFE, and HRS of approximately 500 M dwarfs. For each new observation, we determined its radial velocity and measured its Halpha activity index and its rotation velocity. Additionally, we have multiple observations of many stars to investigate if there are any radial velocity variations due to multiplicity. The results of our survey confirm that early-M dwarfs are Halpha inactive with low rotational velocities and that late-M dwarfs are Halpha active with very high rotational velocities. The results of this high-resolution analysis comprise the most extensive catalogue of rotation and activity in M dwarfs currently available.

Stellar and substellar companions of nearby stars from Gaia DR2 - Binarity from proper motion anomaly

Abstract: The census of stellar and substellar companions of nearby stars is largely incomplete. It is, however, fundamentally important in the understanding of the stellar and planetary formation and evolution mechanisms. We aim to characterize the presence of physical companions in orbit around nearby stars. Orbiting secondary bodies influence the proper motion (PM) of their parent star through their gravitational reflex motion. Using the Hipparcos catalog and Gaia's second data release (GDR2), we determined the long-term PM of the stars common to these two catalogs. We then searched for a PM anomaly (PMa) between the long-term PM vector and the GDR2 (or Hipparcos) measurements, indicative of the presence of a perturbing secondary object. We focussed our analysis on the 6741 nearby stars located within 50 pc, and we also present a catalog of the PMa for >99% of the Hipparcos catalog (>117000 stars). 30% of the stars studied present a PMa greater than 3 sigmas. The PMa allows us to detect orbiting companions, or set stringent limits on their presence. We present a few applications to interesting targets. We set upper limits of 0.1-0.3 M_Jup to potential planets orbiting Proxima between 1 and 10 au (P_orb=3 to 100 years). We confirm that Proxima is gravitationally bound to alpha Cen. We recover the masses of the known companions of epsilon Eri, epsilon Ind, Ross 614 and beta Pic. We also detect a possible planet of a few Jovian masses orbiting tau Ceti. The combination of the GDR2 with Hipparcos results in very high accuracy PMa vectors. The detection of tangential velocity anomalies at a median accuracy of sigma(Delta v_tan)=1.0 m/s per parsec of distance opens the possibility to identify long period orbital companions otherwise inaccessible. Gaia's complementarity to radial velocity and transit techniques already appears to be remarkably powerful.

The Disk Substructures at High Angular Resolution Project (DSHARP): V. Interpreting ALMA maps of protoplanetary disks in terms of a dust model

Abstract: The Disk Substructures at High Angular Resolution Project (DSHARP) is the largest homogeneous high-resolution ($\sim 0.035$ arcsec, or $\sim$ 5 au) disk continuum imaging survey with ALMA so far. In the coming years, many more disks will be mapped with ALMA at similar resolution. Interpreting the results in terms of the properties and quantities of the emitting dusty material is, however, a very non-trivial task. This is in part due to the uncertainty in the dust opacities, an uncertainty which is not likely to be resolved any time soon. It is also partly due to the fact that, as the DSHARP survey has shown, these disk often contain regions of intermediate to high optical depth, even at millimeter wavelengths and at relatively large radius in the disk. This makes the interpretation challenging, in particular if the grains are large and have a large albedo. On the other hand, the highly structured features seen in the DSHARP survey, of which strong indications were already seen in earlier observations, provide a unique opportunity to study the dust growth and dynamics. To provide continuity within the DSHARP project, its follow-up projects, and projects by other teams interested in these data, we present here the methods and opacity choices used within the DSHARP collaboration to link the measured intensity $I_ u$ to dust surface density $\Sigma_d$.

Asteroseismology of 16000 Kepler Red Giants: Global Oscillation Parameters, Masses, and Radii.

Abstract: The Kepler mission has provided exquisite data to perform an ensemble asteroseismic analysis on evolved stars In this work we systematically characterize solar-like oscillations and granulation for 16,094 oscillating red giants, using end-of-mission long-cadence data We produced a homogeneous catalog of the frequency of maximum power (typical uncertainty $\sigma_{ u_{\rm max}}$=16\%), the mean large frequency separation ($\sigma_{\Delta u}$=06\%), oscillation amplitude ($\sigma_{\rm A}$=47\%), granulation power ($\sigma_{\rm gran}$=86\%), power excess width ($\sigma_{\rm width}$=88\%), seismically-derived stellar mass ($\sigma_{\rm M}$=78\%), radius ($\sigma_{\rm R}$=29\%), and thus surface gravity ($\sigma_{\log g}$=001 dex) Thanks to the large red giant sample, we confirm that red-giant-branch (RGB) and helium-core-burning (HeB) stars collectively differ in the distribution of oscillation amplitude, granulation power, and width of power excess, which is mainly due to the mass difference The distribution of oscillation amplitudes shows an extremely sharp upper edge at fixed $ u_{\rm max}$, which might hold clues to understand the excitation and damping mechanisms of the oscillation modes We find both oscillation amplitude and granulation power depend on metallicity, causing a spread of 15\% in oscillation amplitudes and a spread of 25\% in granulation power from [Fe/H]=-07 to 05 dex Our asteroseismic stellar properties can be used as reliable distance indicators and age proxies for mapping and dating galactic stellar populations observed by Kepler They will also provide an excellent opportunity to test asteroseismology using Gaia parallaxes, and lift degeneracies in deriving atmospheric parameters in large spectroscopic surveys such as APOGEE and LAMOST

The origin, early evolution and predictability of solar eruptions

Abstract: Coronal mass ejections (CMEs) were discovered in the early 1970s when space-borne coronagraphs revealed that eruptions of plasma are ejected from the Sun. Today, it is known that the Sun produces eruptive flares, filament eruptions, coronal mass ejections and failed eruptions; all thought to be due to a release of energy stored in the coronal magnetic field during its drastic reconfiguration. This review discusses the observations and physical mechanisms behind this eruptive activity, with a view to making an assessment of the current capability of forecasting these events for space weather risk and impact mitigation. Whilst a wealth of observations exist, and detailed models have been developed, there still exists a need to draw these approaches together. In particular more realistic models are encouraged in order to asses the full range of complexity of the solar atmosphere and the criteria for which an eruption is formed. From the observational side, a more detailed understanding of the role of photospheric flows and reconnection is needed in order to identify the evolutionary path that ultimately means a magnetic structure will erupt.

The spin of the second-born black hole in coalescing binary black holes

Abstract: Various binary black hole formation channels have been proposed since the first gravitational event GW150914 was discovered by the Advanced Laser Interferometer Gravitational-Wave Observatory (AdLIGO). For all evolutionary channels based on the evolution of isolated binaries, the immediate progenitor of the binary black hole is a close binary system composed of a black hole and a helium star. We perform detailed binary evolution and systematically explore the parameter space of initial binary properties, including initial black hole and helium star masses, initial rotation of the helium star as well as metallicity. We argue that the spin of the first-born black hole at its birth is negligible ($\lesssim 0.1$), hence the second-born black hole's spin dominates the measured effective spin, $\chi_{\rm eff}$, from gravitational wave events of double black hole mergers. We find that tides can be important only when orbital periods are shorter than 2 days. Upon core collapse, the helium star produces a black hole (the second-born black hole in the system) with a spin that can span the entire range from zero to maximally spinning. We show that the bimodal distribution of the spin of the second-born black hole obtained in recent papers is mainly due to oversimplifying assumptions. We find an anti-correlation between the merging timescale of the two black holes, T$_{\rm merger}$, and the effective spin $\chi_{\rm eff}$. Finally, we provide new prescriptions for the tidal coefficient E$_2$ for both H-rich and the helium-rich stars. We predict that, with future improvements to AdLIGO's sensitivity, the sample of merging binary black hole systems will show an overdensity of sources with positive but small $\chi_{\rm eff}$. originating from lower mass black hole mergers born at low redshift.

CO destruction in protoplanetary disk midplanes: inside versus outside the CO snow surface

Abstract: CO has long been thought to be the best tracer to measure gas masses as it is readily detected at (sub)mm wavelengths in many disks. Inferred gas masses from CO in recent ALMA observations of large samples of disks seem inconsistent with their inferred dust masses. The derived gas-to-dust mass ratios from CO are 1-2 orders of magnitude lower than $\sim$100 even if photodissociation and freeze-out are included. Herschel measurements of HD line emission imply gas masses in line with gas-to-dust mass ratios of 100. This suggests that at least one additional mechanism is removing gaseous CO. Here we test the suggestion that the bulk of the CO is chemically processed and that the carbon is sequestered into less volatile species such as CO2, CH3OH and CH4 in the dense, shielded midplane regions of the disk. Using our gas-grain chemical code we perform a parameter exploration and follow the CO abundance evolution over a range of conditions representative of shielded disk midplanes. We find that no chemical processing of CO takes place on 1-3 Myr timescales for low cosmic-ray ionisation rates, $< 5\times 10^{-18}$ s$^{-1}$. Assuming an ionisation rate of $10^{-17}$ s$^{-1}$, more than 90% of the CO is destroyed, but only in the parts of the disk below 30 K. This order of magnitude destruction of CO is robust against the choice of grain-surface reaction rate parameters except for temperatures below 20 K, where there is a strong dependence on the H tunnelling barrier. The low temperatures needed for CO chemical processing indicate that the exact disk temperature structure is important, with warm disks around luminous Herbig stars expected to have little to no CO conversion. For cold disks around T Tauri stars, a large fraction of the emitting CO layer is affected unless the disks are young (< 1 Myr). This can lead to inferred gas masses that are up to two orders of magnitude too low.

The Hyades tidal tails revealed by Gaia DR2

Abstract: Within a 200 pc sphere around the Sun we search for the Hyades tidal tails in the Gaia DR2 dataset. We use a modified convergent point method to search for stars with space velocities close to the space velocity of the Hyades cluster. We find a clear indication for the existence of the Hyades tidal tails, a preceding tail extending up to 170 pc from the centre of the Hyades with 292 stars (36 contaminants), and a following tail up to 70 pc with 237 stars (32 contaminants). A comparison with an N-body model of the cluster and its tails shows remarkably good coincidence. Five white dwarfs are found in the tails.

STiC -- A multi-atom non-LTE PRD inversion code for full-Stokes solar observations.

Abstract: The inference of the underlying state of the plasma in the solar chromosphere remains extremely challenging because of the nonlocal character of the observed radiation and plasma conditions in this layer. Inversion methods allow us to derive a model atmosphere that can reproduce the observed spectra by undertaking several physical assumptions. The most advanced approaches involve a depth-stratified model atmosphere described by temperature, line-of-sight velocity, turbulent velocity, the three components of the magnetic field vector, and gas and electron pressure. The parameters of the radiative transfer equation are computed from a solid ground of physical principles. To apply these techniques to spectral lines that sample the chromosphere, NLTE effects must be included in the calculations. We developed a new inversion code STiC to study spectral lines that sample the upper chromosphere. The code is based the RH synthetis code, which we modified to make the inversions faster and more stable. For the first time, STiC facilitates the processing of lines from multiple atoms in non-LTE, also including partial redistribution effects. Furthermore, we include a regularization strategy that allows for model atmospheres with a complex stratification, without introducing artifacts in the reconstructed physical parameters, which are usually manifested in the form of oscillatory behavior. This approach takes steps toward a node-less inversion, in which the value of the physical parameters at each grid point can be considered a free parameter. In this paper we discuss the implementation of the aforementioned techniques, the description of the model atmosphere, and the optimizations that we applied to the code. We carry out some numerical experiments to show the performance of the code and the regularization techniques that we implemented. We made STiC publicly available to the community.

Gaia Data Release 2: Variable stars in the colour-absolute magnitude diagram

Abstract: The ESA Gaia mission provides a unique time-domain survey for more than 1.6 billion sources with G ~ 21 mag. We showcase stellar variability across the Galactic colour-absolute magnitude diagram (CaMD), focusing on pulsating, eruptive, and cataclysmic variables, as well as on stars exhibiting variability due to rotation and eclipses. We illustrate the locations of variable star classes, variable object fractions, and typical variability amplitudes throughout the CaMD and illustrate how variability-related changes in colour and brightness induce `motions' using 22 months worth of calibrated photometric, spectro-photometric, and astrometric Gaia data of stars with significant parallax. To ensure a large variety of variable star classes to populate the CaMD, we crossmatch Gaia sources with known variable stars. We also used the statistics and variability detection modules of the Gaia variability pipeline. Corrections for interstellar extinction are not implemented in this article. Gaia enables the first investigation of Galactic variable star populations across the CaMD on a similar, if not larger, scale than previously done in the Magellanic Clouds. Despite observed colours not being reddening corrected, we clearly see distinct regions where variable stars occur and determine variable star fractions to within Gaia's current detection thresholds. Finally, we show the most complete description of variability-induced motion within the CaMD to date. Gaia enables novel insights into variability phenomena for an unprecedented number of stars, which will benefit the understanding of stellar astrophysics. The CaMD of Galactic variable stars provides crucial information on physical origins of variability in a way previously accessible only for Galactic star clusters or external galaxies.

Strong dipole magnetic fields in fast rotating fully convective stars

Abstract: M dwarfs are the most numerous stars in our Galaxy with masses between approximately 0.5 and 0.1 solar mass. Many of them show surface activity qualitatively similar to our Sun and generate flares, high X-ray fluxes, and large-scale magnetic fields. Such activity is driven by a dynamo powered by the convective motions in their interiors. Understanding properties of stellar magnetic fields in these stars finds a broad application in astrophysics, including, e.g., theory of stellar dynamos and environment conditions around planets that may be orbiting these stars. Most stars with convective envelopes follow a rotation-activity relationship where various activity indicators saturate in stars with rotation periods shorter than a few days. The activity gradually declines with rotation rate in stars rotating more slowly. It is thought that due to a tight empirical correlation between X-ray and magnetic flux, the stellar magnetic fields will also saturate, to values around ~4kG. Here we report the detection of magnetic fields above the presumed saturation limit in four fully convective M-dwarfs. By combining results from spectroscopic and polarimetric studies we explain our findings in terms of bistable dynamo models: stars with the strongest magnetic fields are those in a dipole dynamo state, while stars in a multipole state cannot generate fields stronger than about four kilogauss. Our study provides observational evidence that dynamo in fully convective M dwarfs generates magnetic fields that can differ not only in the geometry of their large scale component, but also in the total magnetic energy.

The Close Binary Fraction of Solar-type Stars is Strongly Anti-correlated with Metallicity.

Abstract: There is now strong evidence that the close binary fraction (P < 10$^4$ days; a < 10 AU) of solar-type stars ($M_1$ = 0.6-1.5M$_{\odot}$) decreases significantly with metallicity. Although early surveys showed that the observed spectroscopic binary (SB) fractions in the galactic disk and halo are similar (e.g., Carney-Latham sample), these studies did not correct for incompleteness. In this study, we examine five different surveys and thoroughly account for their underlying selection biases to measure the intrinsic occurrence rate of close solar-type binaries. We re-analyze: (1) a volume-limited sample of solar-type stars, (2) an SB survey of high-proper-motion stars, (3) various SB samples of metal-poor giants, (4) the APOGEE survey of radial velocity (RV) variables, and (5) Kepler eclipsing binaries (EBs). The observed APOGEE RV variability fraction and Kepler EB fraction both decrease by a factor of $\approx$4 across $-$1.0 < [Fe/H] < 0.5 at the 22$\sigma$ and 9$\sigma$ confidence levels, respectively. After correcting for incompleteness, all five samples exhibit a quantitatively consistent anti-correlation between the intrinsic close binary fraction (a < 10 AU) and metallicity: $F_{\rm close}$ = 53%$\pm$12%, 40%$\pm$6%, 24%$\pm$4%, and 10%$\pm$3% at [Fe/H] = $-$3.0, $-$1.0, $-$0.2 (mean field metallicity), and +0.5, respectively. We present fragmentation models that explain why the close binary fraction of solar-type stars strongly decreases with metallicity while the wide binary fraction, close binary fraction of OB stars, and initial mass function are all constant across $-$1.5 < [Fe/H] < 0.5. The majority of solar-type stars with [Fe/H] < $-$1.0 will interact with a stellar companion, which has profound implications for binary evolution in old and metal-poor environments such as the galactic halo, bulge, thick disk, globular clusters, dwarf galaxies, and high-redshift universe.