On Using the Color-Magnitude Diagram Morphology of M67 to Test Solar Abundances

The open cluster M67 has solar metallicity and an age of about 4 Gyr. The turnoff (TO) mass is close to the minimum mass for which solar metallicity stars develop a convective core during main sequence evolution as a result of the development of hydrogen burning through the CNO cycle. The morphology of the color-magnitude diagram (CMD) of M67 around the TO shows a clear hook-like feature, a direct sign that stars close to the TO have convective cores. VandenBerg et al. investigated the possibility of using the morphology of the M67 TO to put constraints on the solar metallicity, particularly CNO elements, for which solar abundances have been revised downward by more than 30% over the last few years. Here, we extend their work, filling the gaps in their analysis. To this aim, we compute isochrones appropriate for M67 using new (low metallicity) and old (high metallicity) solar abundances and study whether the characteristic TO in the CMD of M67 can be reproduced or not. We also study the importance of other constitutive physics on determining the presence of such a hook, particularly element diffusion, overshooting and nuclear reaction rates. We find that using the new solar abundance determinations, with low CNO abundances, makes it more difficult to reproduce the characteristic CMD of M67. This result is in agreement with results by VandenBerg et al. However, changes in the constitutive physics of the models, particularly overshooting, can influence and alter this result to the extent that isochrones constructed with models using low CNO solar abundances can also reproduce the TO morphology in M67. We conclude that only if all factors affecting the TO morphology are completely under control (and this is not the case), M67 could be used to put constraints on solar abundances

A new view on exoplanet transits: Transit of Venus described using three-dimensional solar atmosphere Stagger-grid simulations

Stellar activity and, in particular, convection-related surface structures, potentially cause fluctuations that can affect the transit light curves. Surface convection simulations can help the interpretation of ToV. We used realistic three-dimensional radiative hydrodynamical simulation of the Sun from the Stagger-grid and synthetic images computed with the radiative transfer code Optim3D to provide predictions for the transit of Venus in 2004 observed by the satellite ACRIMSAT. We computed intensity maps from RHD simulation of the Sun and produced synthetic stellar disk image. We computed the light curve and compared it to the ACRIMSAT observations and also to the light curves obtained with solar surface representations carried out using radial profiles with different limb-darkening laws. We also applied the same spherical tile imaging method to the observations of center-to-limb Sun granulation with HINODE. We managed to explain ACRIMSAT observations of 2004 ToV and showed that the granulation pattern causes fluctuations in the transit light curve. We evaluated the contribution of the granulation to the ToV. We showed that the granulation pattern can partially explain the observed discrepancies between models and data. This confirms that the limb-darkening and the granulation pattern simulated in 3D RHD Sun represent well what is imaged by HINODE. In the end, we found that the Venus's aureole contribution during ToV is less intense than the solar photosphere, and thus negligible. Being able to explain consistently the data of 2004 ToV is a new step forward for 3D RHD simulations that are becoming essential for the detection and characterization of exoplanets. They show that the granulation have to be considered as an intrinsic incertitude, due to the stellar variability, on precise measurements of exoplanet transits of, most likely, planets with small diameters.Comment: Accepted for publication in Astronomy and Astrophysic

The Stagger-grid: A Grid of 3D Stellar Atmosphere Models - I. Methods and General Properties

We present the Stagger-grid, a comprehensive grid of time-dependent, 3D hydrodynamic model atmospheres for late-type stars with realistic treatment of radiative transfer, covering a wide range in stellar parameters. This grid of 3D models is intended for various applications like stellar spectroscopy, asteroseismology and the study of stellar convection. In this introductory paper, we describe the methods used for the computation of the grid and discuss the general properties of the 3D models as well as their temporal and spatial averages (). All our models were generated with the Stagger-code, using realistic input physics for the equation of state (EOS) and for continuous and line opacities. Our ~220 grid models range in Teff from 4000 to 7000K in steps of 500K, in log g from 1.5 to 5.0 in steps of 0.5 dex, and [Fe/H] from -4.0 to +0.5 in steps of 0.5 and 1.0 dex. We find a tight scaling relation between the vertical velocity and the surface entropy jump, which itself correlates with the constant entropy value of the adiabatic convection zone. The range in intensity contrast is enhanced at lower metallicity. The granule size correlates closely with the pressure scale height sampled at the depth of maximum velocity. We compare the models with widely applied 1D models, as well as with theoretical 1D hydrostatic models generated with the same EOS and opacity tables as the 3D models, in order to isolate the effects of using self-consistent and hydrodynamic modeling of convection, rather than the classical mixing length theory approach. For the first time, we are able to quantify systematically over a broad range of stellar parameters the uncertainties of 1D models arising from the simplified treatment of physics, in particular convective energy transport. In agreement with previous findings, we find that the differences can be significant, especially for metal-poor stars.Comment: Accepted for publication in A&A, 31 pages, 29 figure

Non-local Thermodynamic Equilibrium Stellar Spectroscopy with 1D and 〈3D〉 Models. II. Chemical Properties of the Galactic Metal-poor Disk and the Halo

From exploratory studies and theoretical expectations it is known that simplifying approximations in spectroscopic analysis (local thermodynamic equilibrium (LTE), 1D) lead to systematic biases of stellar parameters and abundances. These biases depend strongly on surface gravity, temperature and, in particular, for LTE versus non-LTE (NLTE), on metallicity of the stars. Here we analyze the [Mg/Fe] and [Fe/H] plane of a sample of 326 stars, comparing LTE and NLTE results obtained using 1D hydrostatic models and averaged $\langle 3{\rm{D}}\rangle$ models. We show that compared to the $\langle 3{\rm{D}}\rangle$ NLTE benchmark, the other three methods display increasing biases toward lower metallicities, resulting in false trends of [Mg/Fe] against [Fe/H], which have profound implications for interpretations by chemical evolution models. In our best $\langle 3{\rm{D}}\rangle$ NLTE model, the halo and disk stars show a clearer behavior in the [Mg/Fe]–[Fe/H] plane, from the knee in abundance space down to the lowest metallicities. Our sample has a large fraction of thick disk stars and this population extends down to at least [Fe/H] ~ −1.6 dex, further than previously proven. The thick disk stars display a constant [Mg/Fe] ≈ 0.3 dex, with a small intrinsic dispersion in [Mg/Fe] that suggests that a fast SN Ia channel is not relevant for the disk formation. The halo stars reach higher [Mg/Fe] ratios and display a net trend of [Mg/Fe] at low metallicities, paired with a large dispersion in [Mg/Fe]. These indicate the diverse origin of halo stars from accreted low-mass systems to stochastic/inhomogeneous chemical evolution in the Galactic halo

Asymmetries on red giant branch surfaces from CHARA/MIRC optical interferometry

Context. Red giant branch (RGB) stars are very bright objects in galaxies and are often used as standard candles. Interferometry is the ideal tool to characterize the dynamics and morphology of their atmospheres. Aims. We aim at precisely characterising the surface dynamics of a sample of RGB stars. Methods. We obtained interferometric observations for three RGB stars with the MIRC instrument mounted at the CHARA interfer- ometer. We looked for asymmetries on the stellar surfaces using limb-darkening models. Results. We measured the apparent diameters of HD197989 (Epsilon Cyg) = 4.61+-0.02 mas, HD189276 (HR7633) = 2.95+-0.01 mas, and HD161096 (Beta Oph) = 4.43+-0.01 mas. We detected departures from the centrosymmetric case for all three stars with the tendency of a greater effect for lower logg of the sample. We explored the causes of this signal and conclude that a possible explanation to the interferometric signal is the convection-related and/or the magnetic-related surface activity. However, it is necessary to monitor these stars with new observations, possibly coupled with spectroscopy, in order to firmly establish the cause.Comment: Accepted for publication as a Letter in Astronomy and Astrophysics, section 1. Letters to the Editor. The official date of acceptance is 06/03/201

A single low-energy, iron-poor supernova as the source of metals in the star SMSS J 031300.36-670839.3

The element abundance ratios of four low-mass stars with extremely low metallicities indicate that the gas out of which the stars formed was enriched in each case by at most a few, and potentially only one low-energy, supernova. Such supernovae yield large quantities of light elements such as carbon but very little iron. The dominance of low-energy supernovae is surprising, because it has been expected that the first stars were extremely massive, and that they disintegrated in pair-instability explosions that would rapidly enrich galaxies in iron. What has remained unclear is the yield of iron from the first supernovae, because hitherto no star is unambiguously interpreted as encapsulating the yield of a single supernova. Here we report the optical spectrum of SMSS J031300.36- 670839.3, which shows no evidence of iron (with an upper limit of 10^-7.1 times solar abundance). Based on a comparison of its abundance pattern with those of models, we conclude that the star was seeded with material from a single supernova with an original mass of ~60 Mo (and that the supernova left behind a black hole). Taken together with the previously mentioned low-metallicity stars, we conclude that low-energy supernovae were common in the early Universe, and that such supernovae yield light element enrichment with insignificant iron. Reduced stellar feedback both chemically and mechanically from low-energy supernovae would have enabled first-generation stars to form over an extended period. We speculate that such stars may perhaps have had an important role in the epoch of cosmic reionization and the chemical evolution of early galaxies.Comment: 28 pages, 6 figures, Natur

The STAGGER-grid: A grid of 3D stellar atmosphere models: V. Synthetic stellar spectra and broad-band photometry

Context. The surface structures and dynamics of cool stars are characterised by the presence of convective motions and turbulent flows which shape the emergent spectrum. Aims. We used realistic three-dimensional (3D) radiative hydrodynamical simulations from the STAGGER-grid to calculate synthetic spectra with the radiative transfer code OPTIM3D for stars with different stellar parameters to predict photometric colours and convective velocity shifts. Methods. We calculated spectra from 1000 to 200 000 Å with a constant resolving power of λ/Δλ = 20 000 and from 8470 and 8710 Å (Gaia Radial Velocity Spectrometer - RVS - spectral range) with a constant resolving power of λ/Δλ = 300 000. Results. We used synthetic spectra to compute theoretical colours in the Johnson-Cousins UBV (RI)C, SDSS, 2MASS, Gaia, SkyMapper, Strömgren systems, and HST-WFC3. Our synthetic magnitudes are compared with those obtained using 1D hydrostatic models. We showed that 1D versus 3D differences are limited to a small percent except for the narrow filters that span the optical and UV region of the spectrum. In addition, we derived the effect of the convective velocity fields on selected Fe I lines. We found the overall convective shift for 3D simulations with respect to the reference 1D hydrostatic models, revealing line shifts of between -0.235 and +0.361 km s-1. We showed a net correlation of the convective shifts with the effective temperature: lower effective temperatures denote redshifts and higher effective temperatures denote blueshifts. We conclude that the extraction of accurate radial velocities from RVS spectra need an appropriate wavelength correction from convection shifts. Conclusions. The use of realistic 3D hydrodynamical stellar atmosphere simulations has a small but significant impact on the predicted photometry compared with classical 1D hydrostatic models for late-type stars. We make all the spectra publicly available for the community through the POLLUX database.L.C. gratefully acknowledges support from the Australian Research Council (grants DP150100250, FT160100402)

New constraints on Planck-scale Lorentz Violation in QED from the Crab Nebula

We set constraints on O(E/M) Lorentz Violation in QED in an effective field theory framework. A major consequence of such assumptions is the modification of the dispersion relations for electrons/positrons and photons, which in turn can affect the electromagnetic output of astrophysical objects. We compare the information provided by multiwavelength observations with a full and self-consistent computation of the broad-band spectrum of the Crab Nebula. We cast constraints of order 10^{-5} at 95% confidence level on the lepton Lorentz Violation parameters.Comment: 23 pages, 9 figures. v2: added comments and references, matches version accepted by JCA

Measuring stellar granulation during planet transits

Context. Stellar activity and convection-related surface structures might cause bias in planet detection and characterization that use these transits. Surface convection simulations help to quantify the granulation signal. Aims. We used realistic three-dimensional (3D) radiative hydrodynamical (RHD) simulations from the Stagger grid and synthetic images computed with the radiative transfer code Optim3D to model the transits of three prototype planets: a hot Jupiter, a hot Neptune, and a terrestrial planet. Methods. We computed intensity maps from RHD simulations of the Sun and a K-dwarf star at different wavelength bands from optical to far-infrared that cover the range of several ground-and space-based telescopes which observe exoplanet transits. We modeled the transit using synthetic stellar-disk images obtained with a spherical-tile imaging method and emulated the temporal variation of the granulation intensity generating random images covering a granulation time-series of 13.3 h. We measured the contribution of the stellar granulation on the light curves during the planet transit. Results. We identified two types of granulation noise that act simultaneously during the planet transit: (i) the intrinsic change in the granulation pattern with timescale (e.g., 10 min for solar-type stars assumed in this work) is smaller than the usual planet transit (~hours as in our prototype cases); and (ii) the fact that the transiting planet occults isolated regions of the photosphere that differ in local surface brightness as a result of convective-related surface structures. First, we showed that our modeling approach returns granulation timescale fluctuations that are comparable with what has been observed for the Sun. Then, our statistical approach shows that the granulation pattern of solar and K-dwarf-type stars have a non-negligible effect of the light curve depth during the transit, and, consequentially on the determination of the planet transit parameters such as the planet radius (up to 0.90% and ~0.47% for terrestrial and gaseous planets, respectively). We also showed that larger (or smaller) orbital inclination angles with respect to values corresponding to transit at the stellar center display a shallower transit depth and longer ingress and egress times, but also granulation fluctuations that are correlated to the center-to-limb variation: they increase (or decrease) the value of the inclination, which amplifies the fluctuations. The granulation noise appears to be correlated among the different wavelength ranges either in the visible or in the infrared regions. Conclusions. The prospects for planet detection and characterization with transiting methods are excellent with access to large amounts of data for stars. The granulation has to be considered as an intrinsic uncertainty (as a result of stellar variability) on the precise measurements of exoplanet transits of planets. The full characterization of the granulation is essential for determining the degree of uncertainty on the planet parameters. In this context, the use of 3D RHD simulations is important to measure the convection-related fluctuations. This can be achieved by performing precise and continuous observations of stellar photometry and radial velocity, as we explained with RHD simulations, before, after, and during the transit period

Three-dimensional interferometric, spectrometric, and planetary views of Procyon

We used a new realistic 3D radiative-hydrodynamical model atmosphere of Procyon generated with the Stagger Code and synthetic spectra computed with the radiative transfer code Optim3D to re-analyze interferometric and spectroscopic data from the optical to the infrared of Procyon. We compute intensity maps in two optical filters centered at 500 and 800 nm (MARK III) and one infrared filter centered at 2200 nm (VINCI). We constructed stellar disk images accounting for the center-to-limb variations and used them to derive visibility amplitudes and closure phases. We provide 3D limb-darkening coefficients in the optical as well as in the infrared. We show that visibility curves and closure phases show clear deviations from circular symmetry from the 3rd lobe on. These deviations are detectable with current interferometers using closure phases. We derive new angular diameters at different wavelengths with two independent methods based on 3D simulations. We find a diameter_Vinci = 5.390 \pm 0.03 mas that this is confirmed by an independent asteroseismic estimation. The resulting Teff is 6591 K, which is consistent with the infrared flux method determinations. We find also a value of the surface gravity log g = 4.01 \pm 0.03 that is larger by 0.05 dex from literature values. Spectrophotometric comparisons with observations provide very good agreement with the spectral energy distribution and photometric colors, allowing us to conclude that the thermal gradient of the simulation matches fairly well Procyon. Finally, we show that the granulation pattern of a planet hosting Procyon-like star has a non-negligible impact on the detection of hot Jupiters in the infrared using interferometry closure phases. It is then crucial to have a comprehensive knowledge of the host star to directly detect and characterize hot Jupiters. In this respect, RHD simulations are very important to reach this aim.Comment: Accepted for publication on Astronomy and Astrophysics, 14 pages, 12 figure