Magnetic fields and stellar spots can alter the equivalent widths of absorption lines in stellar spectra, varying during the activity cycle. This also influences the information that we derive through spectroscopic analysis. In this study we analyse high-resolution spectra of 211 Sun-like stars observed at different phases of their activity cycles, in order to investigate how stellar activity affects the spectroscopic determination of stellar parameters and chemical abundances. We observe that equivalent widths of lines can increase as a function of the activity index log R$^\prime_{\rm HK}$ during the stellar cycle, which also produces an artificial growth of the stellar microturbulence and a decrease in effective temperature and metallicity. This effect is visible for stars with activity indexes log R$^\prime_{\rm HK}$$\geq$$-$5.0 (i.e., younger than 4-5 Gyr) and it is more significant at higher activity levels. These results have fundamental implications on several topics in astrophysics that are discussed in the paper, including stellar nucleosynthesis, chemical tagging, the study of Galactic chemical evolution, chemically anomalous stars, the structure of the Milky Way disk, stellar formation rates, photoevaporation of circumstellar disks, and planet hunting.

Low frequency ($\nu\lesssim 150\,$MHz) stellar radio emission is expected to originate in the outer corona at heights comparable to and larger than the stellar radius. Such emission from the Sun has been used to study coronal structure, mass ejections, space-weather conditions around the planets (Schwenn 2006). Searches for low-frequency emission from other stars have only detected a single active flare-star (Lynch et al. 2017) that is not representative of the wider stellar population. Here we report the detection of low-frequency radio emission from a quiescent star, GJ 1151--- a member of the most common stellar type (red dwarf or spectral class M) in the Galaxy. The characteristics of the emission are similar to those of planetary auroral emissions (Zarka 1998) (e.g. Jupiter's decametric emission), suggesting a coronal structure dominated by a global magnetosphere with low plasma density. Our results show that large-scale currents that power radio aurorae operate over a vast range of mass and atmospheric composition, ranging from terrestrial planets to main-sequence stars. The Poynting flux required to produce the observed radio emission cannot be generated by GJ 1151's slow rotation, but can originate in a sub-Alfv\'{e}nic interaction of its magnetospheric plasma with a short-period exoplanet. The emission properties are consistent with theoretical expectations (Zarka 2007; Lanza 2009; Saur et al. 2013; Turnpenney et al. 2018) for interaction with an Earth-size planet in a $\sim 1-5$ day-long orbit.

The recent detection of circularly polarized, long-duration (>8 hr) low-frequency (~150 MHz) radio emission from the M4.5 dwarf GJ 1151 has been interpreted as arising from a star-planet interaction via the electron cyclotron maser instability. The existence or parameters of the proposed planets have not been determined. Using 20 new HARPS-N observations, we put 99th-percentile upper limits on the mass of any close companion to GJ 1151 at Msini < 5.6 M earth. With no stellar, brown dwarf, or giant planet companion likely in a close orbit, our data are consistent with detected radio emission emerging from a magnetic interaction between a short-period terrestrial-mass planet and GJ 1151.

The abundance of short-period planetary systems with high orbital obliquities relative to the spin of their host stars is often taken as evidence that scattering processes play important roles in the formation and evolution of these systems. More recent studies have suggested that wide binary companions can tilt protoplanetary disks, inducing a high stellar obliquity that form through smooth processes like disk migration. DS Tuc Ab, a transiting planet with an 8.138 day period in the 40 Myr Tucana-Horologium association, likely orbits in the same plane as its now-dissipated protoplanetary disk, enabling us to test these theories of disk physics. Here, we report on Rossiter-McLaughlin observations of one transit of DS Tuc Ab with the Planet Finder Spectrograph on the Magellan Clay Telescope at Las Campanas Observatory. We confirm the previously detected planet by modeling the planet transit and stellar activity signals simultaneously. We test multiple models to describe the stellar activity-induced radial velocity variations over the night of the transit, finding the obliquity to be low: $\lambda = 12 \pm 13$ degrees, suggesting that this planet likely formed through smooth disk processes and its protoplanetary disk was not significantly torqued by DS Tuc B. The specific stellar activity model chosen affects the results at the $\approx 5$ degree level. This is the youngest planet to be observed using this technique; we provide a discussion on best practices to accurately measure the observed signal of similar young planets.

In recent years, dedicated extreme-precision radial velocity (EPRV) spectrographs have produced vast quantities of high-resolution, high-signal-to-noise time-series spectra for bright stars. These data contain valuable information for the dual purposes of planet detection via the measured RVs and stellar characterization via the co-added spectra. However, considerable data analysis challenges exist in extracting these data products from the observed spectra at the highest possible precision, including the issue of poorly-characterized telluric absorption features and the common use of an assumed stellar spectral template. In both of these examples, precision-limiting reliance on external information can be sidestepped using the data directly. Here we propose a data-driven method to simultaneously extract precise RVs and infer the underlying stellar and telluric spectra using a linear model (in the log of flux). The model employs a convex objective and convex regularization to keep the optimization of the spectral components fast. We implement this method in wobble, an open-source python package which uses TensorFlow in one of its first non-neural-network applications to astronomical data. In this work, we demonstrate the performance of wobble on archival HARPS spectra. We recover the canonical exoplanet 51 Pegasi b, detect the secular RV evolution of the M dwarf Barnard's Star, and retrieve the Rossiter-McLaughlin effect for the Hot Jupiter HD 189733b. The method additionally produces extremely high-S/N composite stellar spectra and detailed time-variable telluric spectra, which we also present here.

It is well established that the chemical structure of the Milky Way exhibits a bimodality with respect to the $\alpha$-enhancement of stars at a given [Fe/H]. This has been studied largely based on a bulk $\alpha$ abundance, computed as a summary of several individual $\alpha$-elements. Inspired by the expected subtle differences in their nucleosynthetic origins, here we probe the higher level of granularity encoded in the inter-family [Mg/Si] abundance ratio. Using a large sample of stars with APOGEE abundance measurements, we first demonstrate that there is additional information in this ratio beyond what is already apparent in [$\alpha$/Fe] and [Fe/H] alone. We then consider Gaia astrometry and stellar age estimates to empirically characterize the relationships between [Mg/Si] and various stellar properties. We find small but significant trends between this ratio and $\alpha$-enhancement, age, [Fe/H], location in the Galaxy, and orbital actions. To connect these observed [Mg/Si] variations to a physical origin, we attempt to predict the Mg and Si abundances of stars with the galactic chemical evolution model Chempy. We find that we are unable to reproduce abundances for the stars that we fit, which highlights tensions between the yield tables, the chemical evolution model, and the data. We conclude that a more data-driven approach to nucleosynthetic yield tables and chemical evolution modeling is necessary to maximize insights from large spectroscopic surveys.

Improved distance measurements to millisecond pulsars can enhance pulsar timing array (PTA) sensitivity to gravitational-waves, improve tests of general relativity with binary pulsars, improve constraints on fuzzy dark matter, and more. Here we report the parallax and distance measurements to seven Gaia DR2 objects associated with 7 International PTA pulsars: J0437-4715, J1012+5307, J1024-0719, J1732-5049, J1910+1256, J1843-1113, and J1949+3106. By multiplying the posteriors of the PTA-based parallax measurements with the Gaia parallax measurement to the pulsar's companion, we improve the distance measurements, and provide a tentative detection of a previously unknown binary companion to J1843-1113. We also find an order of magnitude improvement in the parallax measurement to J1949+3106. Finally, we recommend that future Gaia data releases use J0437-4715 as a calibration point, since its distance estimate in Gaia DR2 is relatively poor compared to pulsar timing measurements.

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

During its two year prime mission the Transiting Exoplanet Survey Satellite (TESS) will perform a time-series photometric survey covering over 80% of the sky. This survey comprises observations of 26 24 x 96 degree sectors that are each monitored continuously for approximately 27 days. The main goal of TESS is to find transiting planets around 200,000 pre-selected stars for which fixed aperture photometry is recorded every two minutes. However, TESS is also recording and delivering Full-Frame Images (FFIs) of each detector at a 30 minute cadence. We have created an open-source tool, eleanor, to produce light curves for objects in the TESS FFIs. Here, we describe the methods used in eleanor to produce light curves that are optimized for planet searches. The tool performs background subtraction, aperture and PSF photometry, decorrelation of instrument systematics, and cotrending using principal component analysis. We recover known transiting exoplanets in the FFIs to validate the pipeline and perform a limited search for new planet candidates in Sector 1. Our tests indicate that eleanor produces light curves with significantly less scatter than other tools that have been used in the literature. Cadence-stacked images, and raw and detrended eleanor light curves for each analyzed star will be hosted on MAST, with planet candidates on ExoFOP-TESS as Community TESS Objects of Interest (CTOIs). This work confirms the promise that the TESS FFIs will enable the detection of thousands of new exoplanets and a broad range of time domain astrophysics.

Observations of the Kepler-1625 system with the Kepler and Hubble Space Telescopes have suggested the presence of a candidate exomoon, Kepler-1625b I, a Neptune-radius satellite orbiting a long-period Jovian planet. Here we present a new analysis of the Hubble observations, using an independent data reduction pipeline. We find that the transit light curve is well fit with a planet-only model, with a best-fit $\chi^2_\nu$ equal to 1.01. The addition of a moon does not significantly improve the fit quality. We compare our results directly with the original light curve from Teachey & Kipping (2018), and find that we obtain a better fit to the data using a model with fewer free parameters (no moon). We discuss possible sources for the discrepancy in our results, and conclude that the lunar transit signal found by Teachey & Kipping (2018) was likely an artifact of the data reduction. This finding highlights the need to develop independent pipelines to confirm results that push the limits of measurement precision.

(Abridged) The Maunakea Spectroscopic Explorer (MSE) is an end-to-end science platform for the design, execution and scientific exploitation of spectroscopic surveys. It will unveil the composition and dynamics of the faint Universe and impact nearly every field of astrophysics across all spatial scales, from individual stars to the largest scale structures in the Universe. Major pillars in the science program for MSE include (i) the ultimate Gaia follow-up facility for understanding the chemistry and dynamics of the distant Milky Way, including the outer disk and faint stellar halo at high spectral resolution (ii) galaxy formation and evolution at cosmic noon, via the type of revolutionary surveys that have occurred in the nearby Universe, but now conducted at the peak of the star formation history of the Universe (iii) derivation of the mass of the neutrino and insights into inflationary physics through a cosmological redshift survey that probes a large volume of the Universe with a high galaxy density. MSE is positioned to become a critical hub in the emerging international network of front-line astronomical facilities, with scientific capabilities that naturally complement and extend the scientific power of Gaia, the Large Synoptic Survey Telescope, the Square Kilometer Array, Euclid, WFIRST, the 30m telescopes and many more.

To date, a handful of exoplanets have been photometrically mapped using phase-modulated reflection or emission from their surfaces, but the small amplitudes of such signals have limited previous maps almost exclusively to coarse dipolar features on hot giant planets. In this work, we uncover a signal using recently released data from the Transiting Exoplanet Survey Satellite (TESS), which we show corresponds to time-variable reflection from a terrestrial planet with a rotation period of 0.9972696 days. Using a spherical harmonic-based reflection model developed as an extension of the STARRY package, we are able to reconstruct the surface features of this rocky world. We recover a time-variable albedo map of the planet including persistent regions which we interpret as oceans and cloud banks indicative of continental features. We argue that this planet represents the most promising detection of a habitable world to date, although the potential intelligence of any life on it is yet to be determined.

The Maunakea Spectroscopic Explorer (MSE) is a planned 11.25-m aperture facility with a 1.5 square degree field of view that will be fully dedicated to multi-object spectroscopy. A rebirth of the 3.6m Canada-France-Hawaii Telescope on Maunakea, MSE will use 4332 fibers operating at three different resolving powers (R ~ 2500, 6000, 40000) across a wavelength range of 0.36-1.8mum, with dynamical fiber positioning that allows fibers to match the exposure times of individual objects. MSE will enable spectroscopic surveys with unprecedented scale and sensitivity by collecting millions of spectra per year down to limiting magnitudes of g ~ 20-24 mag, with a nominal velocity precision of ~100 m/s in high-resolution mode. This white paper describes science cases for stellar astrophysics and exoplanet science using MSE, including the discovery and atmospheric characterization of exoplanets and substellar objects, stellar physics with star clusters, asteroseismology of solar-like oscillators and opacity-driven pulsators, studies of stellar rotation, activity, and multiplicity, as well as the chemical characterization of AGB and extremely metal-poor stars.

The present work aims to examine in detail the depletion of lithium in solar twins to better constrain stellar evolution models and investigate its possible connection with exoplanets. We employ spectral synthesis in the region of the asymmetric 6707.75 \AA Li I line for a sample of 77 stars plus the Sun. As in previous works based on a smaller sample of solar twins, we find a strong correlation between Li depletion and stellar age. In addition, for the first time we show that the Sun has the lowest Li abundance in comparison with solar twins at similar age (4.6 $\pm$ 0.5 Gyr). We compare the lithium content with the condensation temperature slope for a sub-sample of the best solar twins and determine that the most lithium depleted stars also have fewer refractory elements. We speculate whether the low lithium content in the Sun might be related to the particular configuration of our Solar system.

The stellar Rotation $vs.$ Age relation is commonly considered as a useful tool to derive reliable ages for Sun-like stars. However, in the light of \kepler\ data, the presence of apparently old and fast rotators that do not obey the usual gyrochronology relations led to the hypothesis of weakened magnetic breaking in some stars. In this letter, we constrain the solar rotation evolutionary track using solar twins. Predicted rotational periods as a function of mass, age, [Fe/H] and given critical Rossby number ($Ro_{\rm crit}$) were estimated for the entire rotational sample. Our analysis favors the smooth rotational evolution scenario and suggests that, if the magnetic weakened breaking scenario takes place at all, it should arise after $Ro_{\rm crit}\gtrsim2.29$ or ages $\gtrsim$5.3 Gyr (at 95$\%$ confidence level).

The orbital properties of stars in the disk are signatures of their formation, but they are also expected to change over time due to the dynamical evolution of the Galaxy. Stellar orbits can be quantified by three dynamical actions, J_r, L_z, and J_z, which provide measures of the orbital eccentricity, guiding radius, and non-planarity, respectively. Changes in these dynamical actions over time reflect the strength and efficiency of the evolutionary processes that drive stellar redistributions. We examine how dynamical actions of stars are correlated with their age using two samples of stars with well-determined ages: 78 solar twin stars (with ages to ~5%) and 4376 stars from the APOKASC2 sample (~20%). We compute actions using spectroscopic radial velocities from previous surveys and parallax and proper motion measurements from Gaia DR2. We find weak gradients in all actions with stellar age, of (7.51 +/- 0.52, -29.0 +/- 1.83, 1.54 +/- 0.18) kpc km/s/Gyr for J_r, L_z, and J_z, respectively. There is, however, significant scatter in the action-age relation. We caution that our results will be affected by the restricted spatial extent of our sample, particularly in the case of J_z. Nevertheless, these action-age gradients and their associated variances provide strong constraints on the efficiency of the mechanisms that drive the redistribution of stellar orbits over time and demonstrate that actions are informative as to stellar age. The shallow action-age gradients combined with the large dispersion in each action at a given age, however, renders the prospect of age inference from orbits of individual stars bleak. Using the precision measurements of [Fe/H] and [$\alpha$/Fe] for our stars we investigate the abundance-action relationship and find weak correlations. Similar to our stellar age results, dynamical actions afford little discriminating power between low- and high-$\alpha$ stars.

We have investigated the thorium (Th) abundance in a sample of 53 thin disc solar twins covering a wide range of ages. These data provide constrains on the mantle energy budget of terrestrial planets that can be formed over the evolution of the Galaxy's thin disc. We have estimated Th abundances with an average precision of 0.025\,dex (in both [Th/H] and [Th/Fe]) through comprehensive spectral synthesis of a Th\,II line present at 4019.1290\,{\AA}, using very high resolution (R\,=\,115,000) high quality HARPS spectra obtained at the ESO La Silla Observatory. We have confirmed that there is a large energy budget from Th decay for maintaining mantle convection inside potential rocky planets around solar twins, from the Galactic thin disc formation until now, because the pristine [Th/H]$_{\rm ZAMS}$ is super-solar on average under a uniform dispersion of 0.056\,dex (varying from +0.037 up to +0.138\,dex based on linear fits against isochrone stellar age). Comparing to neodymium (Nd) and europium (Eu), two others neutron-capture elements, the stellar pristine abundance of Th follows Eu along the Galactic thin disc evolution, but it does not follow Nd, probably because neodymium has a significant contribution from the $s$-process (about 60\,per\,cent).

The compositions of stars are a critical diagnostic tool for many topics in astronomy such as the evolution of our Galaxy, the formation of planets, and the uniqueness of the Sun. Previous spectroscopic measurements indicate a large intrinsic variation in the elemental abundance patterns of stars with similar overall metal content. However, systematic errors arising from inaccuracies in stellar models are known to be a limiting factor in such studies, and thus it is uncertain to what extent the observed diversity of stellar abundance patterns is real. Here we report the abundances of 30 elements with precisions of 2% for 79 Sun-like stars within 100 parsecs. Systematic errors are minimized in this study by focusing on solar twin stars and performing a line-by-line differential analysis using high-resolution, high-signal-to-noise spectra. We resolve [X/Fe] abundance trends in galactic chemical evolution at precisions of $10^{-3}$ dex Gyr$^{-1}$ and reveal that stars with similar ages and metallicities have nearly identical abundance patterns. Contrary to previous results, we find that the ratios of carbon-to-oxygen and magnesium-to-silicon in solar metallicity stars are homogeneous to within 10% throughout the solar neighborhood, implying that exoplanets may exhibit much less compositional diversity than previously thought. Finally, we demonstrate that the Sun has a subtle deficiency in refractory material relative to >80% of solar twins (at 2$\sigma$ confidence), suggesting a possible signpost for planetary systems like our own.

It is well known that the magnetic activity of solar type stars decreases with age, but it is widely debated in the literature whether there is a smooth decline or if there is an early sharp drop until 1-2 Gyr followed by a relatively inactive constant phase. We revisited the activity-age relation using time-series observations of a large sample of solar twins whose precise isochronal ages and other important physical parameters have been determined. We measured the Ca II H and K activity indices using 9000 HARPS spectra of 82 solar twins. We measured an average solar activity of $S_{\rm MW}$ = 0.1712 $\pm$ 0.0017 during solar magnetic cycles 23$-$24 covered by HARPS observations and we also inferred an average of $S_{\rm MW}$ = 0.1694 $\pm$ 0.0025 for cycles 10$-$24, anchored on a S index vs. sunspot number correlation. Also, a simple relation between the average and dispersion of the activity levels of solar twins was found. This enabled us to predict the stellar variability effects on the age-activity diagram and, consequently, estimate the chromospheric age uncertainties due to the same phenomena. The age-activity relation is still statistically significant up to ages around 6$-$7 Gyr, in agreement with previous works using open clusters and field stars with precise ages. Our research confirms that Ca II H \& K lines remain a useful chromospheric evolution tracer until stars reach ages of at least 6$-$7 Gyr. We found an evidence that, for the most homogeneous set of old stars, the chromospheric activity indices seem to continue decreasing after the solar age towards the end of the main-sequence. Our results indicate that a significant part of the scatter observed in the age-activity relation of solar twins can be attributed to stellar cycle modulations effects. The Sun seems to have a normal activity level and variability for its age.

We produce light curves for all ~34,000 targets observed with K2 in Campaign 17 (C17), identifying 34 planet candidates, 184 eclipsing binaries, and 222 other periodic variables. The location of the C17 field means follow-up can begin immediately now that the campaign has concluded and interesting targets have been identified. The C17 field has a large overlap with C6, so this latest campaign also offers a rare opportunity to study a large number of targets already observed in a previous K2 campaign. The timing of the C17 data release, shortly before science operations begin with the Transiting Exoplanet Survey Satellite (TESS), also lets us exercise some of the tools and methods developed for identification and dissemination of planet candidates from TESS. We find excellent agreement between these results and those identified using only K2-based tools. Among our planet candidates are several planet candidates with sizes < 4 R_E and orbiting stars with KepMag < 10 (indicating good RV targets of the sort TESS hopes to find) and a Jupiter-sized single-transit event around a star already hosting a 6 d planet candidate.

Important insights into the formation and evolution of the Galactic disc(s) are contained in the chemical compositions of stars. We analysed high-resolution and high signal to noise HARPS spectra of 79 solar twin stars in order to obtain precise determinations of their atmospheric parameters, ages ($\sigma$$\sim$0.4 Gyr) and chemical abundances ($\sigma$$<$0.01~dex) of 12 neutron-capture elements (Sr, Y, Zr, Ba, La, Ce, Pr, Nd, Sm, Eu, Gd, and Dy). This valuable dataset allows us to study the [X/Fe]-age relations over a time interval of $\sim$10 Gyr and among stars belonging to the thin and thick discs. These relations show that i) the $s$-process has been the main channel of nucleosynthesis of $n$-capture elements during the evolution of the thin disc; ii) the thick disc is rich in $r$-process elements which suggests that its formation has been rapid and intensive. %; iii) a chemical continuity between the thin and thick discs is evident in the abundances of Ba. In addition, the heavy (Ba, La, Ce) and light (Sr, Y, Zr) $s$-process elements revealed details on the dependence between the yields of AGB stars and the stellar mass or metallicity. Finally, we confirmed that both [Y/Mg] and [Y/Al] ratios can be employed as stellar clocks, allowing ages of solar twin stars to be estimated with an average precision of $\sim$0.5~Gyr.

Previous studies on the rotation of Sun-like stars revealed that the rotational rates of young stars converge towards a well-defined evolution that follows a power-law decay. It seems, however, that some binary stars do not obey this relation, often by displaying enhanced rotational rates and activity. In the Solar Twin Planet Search program we observed several solar twin binaries, and found a multiplicity fraction of $42\% \pm 6\%$ in the whole sample; moreover, at least three of these binaries (HIP 19911, HIP 67620 and HIP 103983) clearly exhibit the aforementioned anomalies. We investigated the configuration of the binaries in the program, and discovered new companions for HIP 6407, HIP 54582, HIP 62039 and HIP 30037, of which the latter is orbited by a $0.06$ M$_\odot$ brown dwarf in a 1-month long orbit. We report the orbital parameters of the systems with well-sampled orbits and, in addition, the lower limits of parameters for the companions that only display a curvature in their radial velocities. For the linear trend binaries, we report an estimate of the masses of their companions when their observed separation is available, and a minimum mass otherwise. We conclude that solar twin binaries with low-mass stellar companions at moderate orbital periods do not display signs of a distinct rotational evolution when compared to single stars. We confirm that the three peculiar stars are double-lined binaries, and that their companions are polluting their spectra, which explains the observed anomalies.

HD3167 is a bright (V=8.9 mag) K0V star observed by the NASA's K2 space mission during its Campaign 8. It has been recently found to host two small transiting planets, namely, HD3167b, an ultra short period (0.96 d) super-Earth, and HD3167c, a mini-Neptune on a relatively long-period orbit (29.85 d). Here we present an intensive radial velocity follow-up of HD3167 performed with the FIES@NOT, HARPS@ESO-3.6m, and HARPS-N@TNG spectrographs. We revise the system parameters and determine radii, masses, and densities of the two transiting planets by combining the K2 photometry with our spectroscopic data. With a mass of 5.69+/-0.44 MEarth, radius of 1.574+/-0.054 REarth, and mean density of 8.00(+1.0)(-0.98) g/cm^3, HD3167b joins the small group of ultra-short period planets known to have a rocky terrestrial composition. HD3167c has a mass of 8.33 (+1.79)(-1.85) MEarth and a radius of 2.740(+0.106)(-0.100) REarth, yielding a mean density of 2.21(+0.56)(-0.53) g/cm^3, indicative of a planet with a composition comprising a solid core surrounded by a thick atmospheric envelope. The rather large pressure scale height (about 350 km) and the brightness of the host star make HD3167c an ideal target for atmospheric characterization via transmission spectroscopy across a broad range of wavelengths. We found evidence of additional signals in the radial velocity measurements but the currently available data set does not allow us to draw any firm conclusion on the origin of the observed variation.

The six planets of the Kepler-11 system are the archetypal example of a population of surprisingly low-density transiting planets revealed by the Kepler mission. We have determined the fundamental parameters and chemical composition of the Kepler-11 host star to unprecedented precision using an extremely high quality spectrum from Keck-HIRES (R$\simeq$67,000, S/N per pixel$\simeq$260 at 600 nm). Contrary to previously published results, our spectroscopic constraints indicate that Kepler-11 is a young main-sequence solar twin. The revised stellar parameters and new analysis raise the densities of the Kepler-11 planets by between 20-95% per planet, making them more typical of the emerging class of "puffy" close-in exoplanets. We obtain photospheric abundances of 22 elements and find that Kepler-11 has an abundance pattern similar to that of the Sun with a slightly higher overall metallicity. We additionally analyze the Kepler lightcurves using a photodynamical model and discuss the tension between spectroscopic and transit/TTV-based stellar density estimates.

We present the open-source radiative transfer code named HELIOS, which is constructed for studying exoplanetary atmospheres. In its initial version, the model atmospheres of HELIOS are one-dimensional and plane-parallel, and the equation of radiative transfer is solved in the two-stream approximation with non-isotropic scattering. A small set of the main infrared absorbers is employed, computed with the opacity calculator HELIOS-K and combined using a correlated-$k$ approximation. The molecular abundances originate from validated analytical formulae for equilibrium chemistry. We compare HELIOS with the work of Miller-Ricci & Fortney using a model of GJ 1214b, and perform several tests, where we find: model atmospheres with single-temperature layers struggle to converge to radiative equilibrium; $k$-distribution tables constructed with $\gtrsim 0.01$ cm$^{-1}$ resolution in the opacity function ($ \lesssim 10^3$ points per wavenumber bin) may result in errors $\gtrsim 1$-10 % in the synthetic spectra; and a diffusivity factor of 2 approximates well the exact radiative transfer solution in the limit of pure absorption. We construct "null-hypothesis" models (chemical equilibrium, radiative equilibrium and solar element abundances) for 6 hot Jupiters. We find that the dayside emission spectra of HD 189733b and WASP-43b are consistent with the null hypothesis, while it consistently under-predicts the observed fluxes of WASP-8b, WASP-12b, WASP-14b and WASP-33b. We demonstrate that our results are somewhat insensitive to the choice of stellar models (blackbody, Kurucz or PHOENIX) and metallicity, but are strongly affected by higher carbon-to-oxygen ratios. The code is publicly available as part of the Exoclimes Simulation Platform (ESP; exoclime.net).

[Methods]. We obtained high-precision radial velocities with HARPS on the ESO 3.6 m telescope and determined precise stellar elemental abundances (~0.01 dex) using MIKE spectra on the Magellan 6.5m telescope. [Results]. Our data indicate the presence of a planet with a minimum mass of 26 Earth masses around the solar twin HIP 68468. The planet is a super-Neptune, but unlike the distant Neptune in our solar system (30 AU), HIP 68468c is close-in, with a semi-major axis of 0.66 AU, similar to that of Venus. The data also suggest the presence of a super-Earth with a minimum mass of 2.9 Earth masses at 0.03 AU; if the planet is confirmed, it will be the fifth least massive radial velocity planet discovery to date and the first super-Earth around a solar twin. Both isochrones (5.9 Gyr) and the abundance ratio [Y/Mg] (6.4 Gyr) indicate an age of about 6 billion years. The star is enhanced in refractory elements when compared to the Sun, and the refractory enrichment is even stronger after corrections for Galactic chemical evolution. We determined a NLTE Li abundance of 1.52 dex, which is four times higher than what would be expected for the age of HIP 68468. The older age is also supported by the low log(R'HK) (-5.05) and low jitter. Engulfment of a rocky planet of 6 Earth masses can explain the enhancement in both lithium and the refractory elements. [Conclusions]. The super-Neptune planet candidate is too massive for in situ formation, and therefore its current location is most likely the result of planet migration that could also have driven other planets towards its host star, enhancing thus the abundance of lithium and refractory elements in HIP 68468. The intriguing evidence of planet accretion warrants further observations to verify the existence of the planets that are indicated by our data and to better constrain the nature of the planetary system around this unique star.

We report the discovery of K2-98 b (EPIC 211391664 b), a transiting Neptune-sized planet monitored by the K2 mission during its campaign 5. We combine the K2 time-series data with ground-based photometric and spectroscopic follow-up observations to confirm the planetary nature of the object and derive its mass, radius, and orbital parameters. K2-98 b is a warm Neptune-like planet in a 10-day orbit around a V=12.2~mag F-type star with $M_\star$=$ 1.074\pm0.042$, $R_\star$=$ 1.311 ^{+ 0.083} _{ - 0.048} $, and age of $5.2_{-1.0}^{+1.2}$~Gyr. We derive a planetary mass and radius of $M_\mathrm{p}$=$ 32.2 \pm 8.1 $ and $R_\mathrm{p}$=$4.3^{+0.3}_{-0.2}$. K2-98 b joins the relatively small group of Neptune-sized planets whose both mass and radius have been derived with a precision better than 25 %. We estimate that the planet will be engulfed by its host star in $\sim$3~Gyr, due to the evolution of the latter towards the red giant branch.

It is still unclear how common the Sun is when compared to other similar stars in regards to some of its physical properties, such as rotation. Considering that gyrochronology relations are widely used today to estimate ages of stars in the main sequence, and that the Sun is used to calibrate it, it is crucial to assess if these procedures are acceptable. We analyze the rotational velocities -- limited by the unknown rotation axis inclination angle -- of an unprecedented large sample of solar twins in order to study the rotational evolution of Sun-like stars, and assess if the Sun is a typical rotator. We use high-resolution ($R = 115000$) spectra obtained with the HARPS spectrograph and ESO's 3.6 m telescope at La Silla Observatory. The projected rotational velocities for 82 solar twins are estimated by line profile fitting with synthetic spectra. Macroturbulence velocities are inferred from a prescription that accurately reflects their dependence with effective temperature and luminosity of the stars. Our sample of solar twins include some spectroscopic binaries with enhanced rotational velocities, and we do not find any non-spectroscopic binaries with unusually high rotation velocities. We verified that the Sun does not have a peculiar rotation, but the solar twins exhibit rotational velocities that depart from the Skumanich relation. The Sun is a regular rotator when compared to solar twins with a similar age. Additionally, we obtain a rotational braking law that better describes the stars in our sample ($v \propto t^{-0.6}$) in contrast to previous, often-used scalings.

Solar twins are stars with similar stellar (surface) parameters to the Sun that can have a wide range of ages. This provide an opportunity to analyze the variation of their chemical abundances with age. Nissen (2015) recently suggested that the abundances of the s-process element Y and the $\alpha$-element Mg could be used to estimate stellar ages. This paper aims to determine with high precision the Y, Mg, and Fe abundances for a sample of 88 solar twins that span a broad age range ($0.3-10.0$\,Gyr) and investigate their use for estimating ages. We obtained high-quality Magellan Inamori Kyocera Echelle (MIKE) spectra and determined Y and Mg abundances using equivalent widths and a line-by-line differential method within a 1D LTE framework. Stellar parameters and iron abundances were measured in Paper I of this series for all stars, but a few (three) required a small revision. The [Y/Mg] ratio shows a strong correlation with age. It has a slope of -0.041$\pm$0.001 dex/Gyr and a significance of 41 $\sigma$. This is in excellent agreement with the relation first proposed by Nissen (2015). We found some outliers that turned out to be binaries where mass transfer may have enhanced the yttrium abundance. Given a precise measurement of [Y/Mg] with typical error of 0.02 dex in solar twins, our formula can be used to determine a stellar age with $\sim$0.8 Gyr precision in the 0 to 10 Gyr range.

Through our HARPS radial velocity survey for planets around solar twin stars, we have identified a promising Jupiter twin candidate around the star HIP11915. We characterize this Keplerian signal and investigate its potential origins in stellar activity. Our analysis indicates that HIP11915 hosts a Jupiter-mass planet with a 3800-day orbital period and low eccentricity. Although we cannot definitively rule out an activity cycle interpretation, we find that a planet interpretation is more likely based on a joint analysis of RV and activity index data. The challenges of long-period radial velocity signals addressed in this paper are critical for the ongoing discovery of Jupiter-like exoplanets. If planetary in nature, the signal investigated here represents a very close analog to the solar system in terms of both Sun-like host star and Jupiter-like planet.

The achievable level of precision on photospheric abundances of stars is a major limiting factor on investigations of exoplanet host star characteristics, the chemical histories of star clusters, and the evolution of the Milky Way and other galaxies. While model-induced errors can be minimized through the differential analysis of spectrally similar stars, the maximum achievable precision of this technique has been debated. As a test, we derive differential abundances of 19 elements from high-quality asteroid-reflected solar spectra taken using a variety of instruments and conditions. We treat the solar spectra as being from unknown stars and use the resulting differential abundances, which are expected to be zero, as a diagnostic of the error in our measurements. Our results indicate that the relative resolution of the target and reference spectra is a major consideration, with use of different instruments to obtain the two spectra leading to errors up to 0.04 dex. Use of the same instrument at different epochs for the two spectra has a much smaller effect (~0.007 dex). The asteroid used to obtain the solar standard also has a negligible effect (~0.006 dex). Assuming that systematic errors from the stellar model atmospheres have been minimized, as in the case of solar twins, we confirm that differential chemical abundances can be obtained at sub-0.01 dex precision with due care in the observations, data reduction and abundance analysis.

We are carrying out a search for planets around a sample of solar twin stars using the HARPS spectrograph. The goal of this project is to exploit the advantage offered by solar twins to obtain chemical abundances of unmatched precision. This survey will enable new studies of the stellar composition -- planet connection. Here we used the MIKE spectrograph on the Magellan Clay Telescope to acquire high resolution, high signal-to-noise ratio spectra of our sample stars. We measured the equivalent widths of iron lines and used strict differential excitation/ionization balance analysis to determine atmospheric parameters of unprecedented internal precision (DTeff=7K, Dlogg=0.019, D[Fe/H]=0.006dex, Dvt=0.016km/s). Reliable relative ages and highly precise masses were then estimated using theoretical isochrones. The spectroscopic parameters we derived are in good agreement with those measured using other independent techniques. The root-mean-square scatter of the differences seen is fully compatible with the observational errors, demonstrating, as assumed thus far, that systematic uncertainties in the stellar parameters are negligible in the study of solar twins. We find a tight activity-age relation for our sample stars, which validates the internal precision of our dating method. Furthermore, we find that the solar cycle is perfectly consistent both with this trend and its star-to-star scatter. We present the largest sample of solar twins analyzed homogeneously using high quality spectra. The fundamental parameters derived from this work will be employed in subsequent work that aims to explore the connections between planet formation and stellar chemical composition.

We study with unprecedented detail the chemical composition and stellar parameters of the solar twin 18 Sco in a strictly differential sense relative to the Sun. Our study is mainly based on high resolution (R ~ 110 000) high S/N (800-1000) VLT UVES spectra, which allow us to achieve a precision of about 0.005 dex in differential abundances. The effective temperature and surface gravity of 18 Sco are Teff = 5823+/-6 K and log g = 4.45+/-0.02 dex, i.e., 18 Sco is 46+/-6 K hotter than the Sun and log g is 0.01+/-0.02 dex higher. Its metallicity is [Fe/H] = 0.054+/-0.005 dex and its microturbulence velocity is +0.02+/-0.01 km/s higher than solar. Our precise stellar parameters and differential isochrone analysis show that 18 Sco has a mass of 1.04+/-0.02M_Sun and that it is ~1.6 Gyr younger than the Sun. We use precise HARPS radial velocities to search for planets, but none were detected. The chemical abundance pattern of 18 Sco displays a clear trend with condensation temperature, showing thus higher abundances of refractories in 18 Sco than in the Sun. Intriguingly, there are enhancements in the neutron-capture elements relative to the Sun. Despite the small element-to-element abundance differences among nearby n-capture elements (~0.02 dex), we successfully reproduce the r-process pattern in the solar system. This is independent evidence for the universality of the r-process. Our results have important implications for chemical tagging in our Galaxy and nucleosynthesis in general.

We present the first detailed chemical abundance analysis of the old 8.2 Gyr solar twin, HIP 102152. We derive differential abundances of 21 elements relative to the Sun with precisions as high as 0.004 dex ($\lesssim$1%), using ultra high-resolution (R = 110,000), high S/N UVES spectra obtained on the 8.2-m Very Large Telescope. Our determined metallicity of HIP 102152 is [Fe/H] = -0.013 $\pm$ 0.004. The atmospheric parameters of the star were determined to be 54 K cooler than the Sun, 0.09 dex lower in surface gravity, and a microturbulence identical to our derived solar value. Elemental abundance ratios examined vs. dust condensation temperature reveal a solar abundance pattern for this star, in contrast to most solar twins. The abundance pattern of HIP 02152 appears to be the most similar to solar of any known solar twin. Abundances of the younger, 2.9 Gyr solar twin, 18 Sco, were also determined from UVES spectra to serve as a comparison for HIP 102152. The solar chemical pattern of HIP 102152 makes it a potential candidate to host terrestrial planets, which is reinforced by the lack of giant planets in its terrestrial planet region. The following non-local thermodynamic equilibrium Li abundances were obtained for HIP 102152, 18 Sco, and the Sun: log $\epsilon$ (Li) = 0.48 $\pm$ 0.07, 1.62 $\pm$ 0.02, and 1.07 $\pm$ 0.02, respectively. The Li abundance of HIP 102152 is the lowest reported to date for a solar twin, and allows us to consider an emerging, tightly constrained Li-age trend for solar twin stars.